Mathematical Problems in Engineering

Mathematical Problems in Engineering / 2012 / Article

Research Article | Open Access

Volume 2012 |Article ID 490649 |

Jian-Jun Shu, Kian Yan Yong, Weng Kong Chan, "An Improved Scoring Matrix for Multiple Sequence Alignment", Mathematical Problems in Engineering, vol. 2012, Article ID 490649, 9 pages, 2012.

An Improved Scoring Matrix for Multiple Sequence Alignment

Academic Editor: Gerhard-Wilhelm Weber
Received16 May 2011
Revised05 Oct 2011
Accepted05 Oct 2011
Published11 Jan 2012


The way for performing multiple sequence alignment is based on the criterion of the maximum-scored information content computed from a weight matrix, but it is possible to have two or more alignments to have the same highest score leading to ambiguities in selecting the best alignment. This paper addresses this issue by introducing the concept of joint weight matrix to eliminate the randomness in selecting the best multiple sequence alignment. Alignments with equal scores are iteratively rescored with the joint weight matrix of increasing level (nucleotide pairs, triplets, and so on) until one single best alignment is eventually found. This method for resolving ambiguity in multiple sequence alignment can be easily implemented by use of the improved scoring matrix.

1. Introduction

In the search for DNA regulatory elements such as binding sites, promoter, donor sites, TATA box, and genes, the multiple sequences containing these elements have to be aligned against each other. These elements are highly but not absolutely conserved and a weight matrix is used to represent and score the multiple sequences [1]. However, the current motif discovery algorithms based on the weight matrix technique for scoring multiple sequence alignment in terms of information content are not without their limitations [2]. From the analysis of these algorithms, the highest performance coefficient on the binding site level of search is only 30.2% using Motif Sampler [3], which is an algorithm modified from the widely adopted Gibbs sampling method [4]. This may be a result of randomness in selecting the best alignment from cases whereby there are multiple peaks. Hence, there are rooms for improvement, which is evident from many different approaches that have been developed [5โ€“9].

In this paper, a method of removing the randomness in selection is proposed. Randomness in selection occurs when there is more than one choice of alignments with the highest information content [10]. If one peak is randomly selected, the accuracy of multiple sequence alignment is compromised. This may be the reason that methods based on applied information theory cannot achieve much higher sensitivity, specificity, and performance. For example, by randomly selecting two peaks of similar information content, there is a 50% chance of selecting the wrong peak.

In order to overcome this problem, a simple method is proposed to eliminate the randomness of peak selection and to provide the best alignment, through the use of joint weight matrix (JWM) in this paper. Its flexibility means that a higher-level JWM can be used to work with cases with multiple peaks. The higher the level of JWM used, the lesser will the number of peaks be, until eventually a single peak is obtained. In this paper, JWM has been shown to reduce successfully the number of peaks in multiple sequence alignment.

2. Systems and Methods

The concept of JWM is presented here to demonstrate how two or more ambiguous selections can be reduced. Two sequences are used in this example. The longer one represents the DNA sequence and the shorter one represents a motif sequence, which is aligned to the former. The motif is assumed to be a perfect weight matrix with 100% base weightage at each position. The score is then either 1 for match or 0 for mismatch at each position for simplicity of demonstration.

Since the sequence is 7โ€‰bp (base pair) long and the motif is 4โ€‰bp, the total number of possible shift positions without introducing gaps is 7โˆ’4+1=4 in Figure 1. Table 1 shows the sequence alignment.



The score for the four possible alignments presents an ambiguous choice between positions 1 and 2, which are possible alignments with the highest score of 2. Since there is more than one peak or alignment, the second-level JWM is used to score the alignment. Table 2 shows the result using the second-level comparison.



The result clearly shows that between positions 1 and 2, the better match for the motif with the DNA is the position 2 with the matching score of 1 as compared with the position 1 with 0.

3. Algorithm

Here it is shown how JWM can be integrated with sequence alignment tool to remove the randomness of selection during the alignment process. The following are the additional steps added using JWM.

Step 1. Determine a weight matrix ๐‘ค(๐‘,๐‘–)=๐‘›(๐‘,๐‘–)โˆ‘๐‘โˆˆ{A,T,G,C}๐‘›(๐‘,๐‘–),(3.1) where ๐‘›(๐‘,๐‘–) is the number of each base ๐‘โˆˆ{A,T,G,C} at each position ๐‘–.

Step 2. Calculate the second-level JWM ๐‘ค2๎€ท๐‘1๐‘2๎€ธ๎€ท๐‘,๐‘–=๐‘ค1๎€ธ๐‘ค๎€ท๐‘,๐‘–2๎€ธ,๐‘–+1.(3.2) For the second-level JWM, the number of possible combinations of the four bases is 42=16. Hence, JWM is a matrix size of 16 by window length.

Step 3. From the weight matrix, the uncertainty of each combination of bases is ๐ป๐‘ (๐‘–)=โˆ’๐‘š๎“๐‘—=1๎“๐‘๐‘—โˆˆ{A,T,G,C}๐‘ค๐‘š๎€ท๐‘1๐‘2โ‹ฏ๐‘๐‘š๎€ธ,๐‘–log2๐‘ค๐‘š๎€ท๐‘1๐‘2โ‹ฏ๐‘๐‘š๎€ธ,๐‘–.(3.3) For the second-level JWM, ๐‘š is the value of 2.

Step 4. The information content for each base is then ๐‘…(๐‘–)=2๐‘š[]โˆ’๐ป๐‘ (๐‘–)โˆ’๐‘’๐‘›(๐‘–),(3.4) where ๐‘’[๐‘›(๐‘–)] is a small sample correction for ๐ป๐‘ (๐‘–) [11].

Step 5. The score for one shifting position is then ๐‘…shift(๎“sp)=๐‘–๐‘…(๐‘–).(3.5) The shift position (sp) ranges from negative to positive shifting parameter.

Step 6. Shift JWM as predetermined to get the alignment score plot of information content versus shifting position. From the alignment score plot, the highest peak is chosen among the ambiguous choice of the previous set of peaks to be generated.

Step 7. If there is still ambiguity after using the second-level JWM, a higher-level JWM (three or higher) should be calculated ๐‘ค๐‘š๎€ท๐‘1๐‘2โ‹ฏ๐‘๐‘š๎€ธ๎€ท๐‘,๐‘–=๐‘ค1๎€ธ๐‘ค๎€ท๐‘,๐‘–2๎€ธ๎€ท๐‘,๐‘–+1โ‹ฏ๐‘ค๐‘š๎€ธ,๐‘–+๐‘šโˆ’1.(3.6) Repeat the Steps 3 to 6 using the higher-level JWM in (3.6) when there is ambiguity in peak selection if using any lower-level JWM.

4. Implementation

An example of how JWM is used to eliminate or reduce ambiguity is shown using data from 16 randomly generated sequences of 15โ€‰bp (Tables 3 and 4) that bind to OxyR [12]. For illustration purpose, the centre 9th base is taken to be the start site of transcription, labeled as the position 0. The alignment score is obtained by using the window of 5 bases from โˆ’1 to +3 and the range of shifting position set from โˆ’8 to +6 with respect to the start site. The sequences are shifted one base at a time and the new alignment score is recalculated based on the simplified sequence logo [13] in Figure 2(a).







Window and shifting parameters are selected such that an ambiguous choice of more than one peak is resolved. By shifting one of the sequences from โˆ’8 to +6, the alignment score based on window from โˆ’1 to +3 show two peaks at shift positions โˆ’5 and 0 in Figure 2(b). From the simplified sequence logo, the information content prior to shifting of any sequence is ๐‘…shift(0)=0.0637+0.1950+0.5087+0.1504+0.2500=1.1678bits.(4.1)

One of the sequences is randomly selected and shifted about its position. The weight matrix is calculated for each new position and a set of ๐‘…shift(sp) is obtained by the end of the shift. The amount of shift required for the 16 sequences to produce ๐‘…shift is plotted in Figure 2(b), where two peaks are located at shift positions โˆ’5 and 0. The situation is ambiguous and a higher-level search is required by using JWM. The weight matrix is replaced by the second-level JWM in the new search. The new ๐‘…shift plot based on the higher-level JWM is shown in Figure 2(d).

The new alignment score using JWM shows clearly that the shift position 0 has higher information content, as compared with the shift position โˆ’5. Hence, the best alignment is the original position 0. Instead of randomly selecting one of the peaks, it is rational to select the peak with higher information content.

5. Discussion and Conclusions

In the selection of the best multiple sequence alignment using the conventional weight matrix, it is assumed that the probability of each base is independent of its neighboring one. Output from a multiple sequence alignment program is not always the same. This can be attributed to several factors. One of the important contributing factors is the conventional scoring matrix. The best alignment at each stage is decided by the highest score with the conventional scoring matrix. However, there are cases whereby there is more than one of such score. This creates an ambiguity in selecting the best alignment. A random choice can be made, but it may result in a less than optimal alignment.

The following shows examples of ambiguities found using the conventional scoring matrix. The benchmark database (Table 5) consists of DNA sequences containing amelogenin protein in the study of its origin and evolutionary path [14]. Cases of ambiguity during multiple sequence alignment using the conventional scoring matrix are shown in Figure 3. For example, the ambiguity is found when the sequence 2 (DMSPARC), at position 0 and 3 of window, with window placed at the 18th base from the start (first base on the left).

No.NamesSequencesWindow positionAmbiguity positions

180, 3
180, 3
180, 3
11โˆ’5, โˆ’4
20โˆ’4, 0
7โˆ’1, 0
20โˆ’4, 0
20โˆ’4, 0
20โˆ’4, 0
20โˆ’4, 0
16โˆ’3, 0
19โˆ’3, 0
160, 3
7โˆ’4, 5
160, 3

The examples above show that ambiguities are frequent enough to be of concern during multiple sequence alignment, which may result in a suboptimal alignment. This problem can be overcome by using the proposed joint weight matrix for scoring. The proposed scoring matrix allows a closer look at each alignment by considering two or more bases for each scoring element. By comparing two bases at one time, the probability of the next base is affected by what appears before it. In fact, there are 16 probabilities of a pair of bases as compared with just 4 probabilities if only one base is considered. This increases the depth of search to reduce the number of peaks. Under Section 4, it is shown how the second-level JWM can identify the highest peak when a conventional weight matrix could not. This reduces the error that may occur when โ€œconflicts are resolvedโ€ by making a โ€œpseudorandom choiceโ€ [10].

The higher-level of JWM can be used depending on the level of accuracy required. For example, the second-level JWM may be able to reduce the number of peaks from 5 to 3. The randomness is reduced when one is choosing the best peak from 3 instead of 5 possible sites. However, if the application requires a level of match to be of greater accuracy, a higher-level of JWM may be needed to proceed. The higher-level of JWM can further filter out more peaks till only one obvious choice is left. Although the higher-level of JWM may require more computation time and additional scan, this may be compensated by the faster convergence of results as a better alignment is selected early in the iterations. This is true especially for cases whereby a large number of iterations are required before a satisfactory convergence can be found [15]. JWM can be used to improve applications using conventional weight matrix system in bioinformatics. Besides aligning DNA sequences, JWM can also be implemented in protein sequence alignment.


  1. G. D. Stormo and G. W. Hartzell, โ€œIdentifying protein-binding sites from unaligned DNA fragments,โ€ Proceedings of the National Academy of Sciences of the United States of America, vol. 86, no. 4, pp. 1183โ€“1187, 1989. View at: Google Scholar
  2. J. Hu, B. Li, and D. Kihara, โ€œLimitations and potentials of current motif discovery algorithms,โ€ Nucleic Acids Research, vol. 33, no. 15, pp. 4899โ€“4913, 2005. View at: Publisher Site | Google Scholar
  3. G. Thijs, K. Marchal, M. Lescot et al., โ€œA Gibbs sampling method to detect overrepresented motifs in the upstream regions of coexpressed genes,โ€ Journal of Computational Biology, vol. 9, no. 2, pp. 447โ€“464, 2002. View at: Publisher Site | Google Scholar
  4. C. E. Lawrence, S. F. Altschul, M. S. Boguski, J. S. Liu, A. F. Neuwald, and J. C. Wootton, โ€œDetecting subtle sequence signals: a gibbs sampling strategy for multiple alignment,โ€ Science, vol. 262, no. 5131, pp. 208โ€“214, 1993. View at: Google Scholar
  5. Y. Liu, X. S. Liu, L. Wei, R. B. Altman, and S. Batzoglou, โ€œEukaryotic regulatory element conservation analysis and identification using comparative genomics,โ€ Genome Research, vol. 14, no. 3, pp. 451โ€“458, 2004. View at: Publisher Site | Google Scholar
  6. J. J. Shu and L. S. Ouw, โ€œPairwise alignment of the DNA sequence using hypercomplex number representation,โ€ Bulletin of Mathematical Biology, vol. 66, no. 5, pp. 1423โ€“1438, 2004. View at: Publisher Site | Google Scholar
  7. A. V. Favorov, M. S. Gelfand, A. V. Gerasimova, D. A. Ravcheev, A. A. Mironov, and V. J. Makeev, โ€œA Gibbs sampler for identification of symmetrically structured, spaced DNA motifs with improved estimation of the signal length,โ€ Bioinformatics, vol. 21, no. 10, pp. 2240โ€“2245, 2005. View at: Publisher Site | Google Scholar
  8. L. Kuo and T. Y. Yang, โ€œAn improved collapsed Gibbs sampler for Dirichlet process mixing models,โ€ Computational Statistics & Data Analysis, vol. 50, no. 3, pp. 659โ€“674, 2006. View at: Publisher Site | Google Scholar
  9. J. J. Shu and Y. Li, โ€œHypercomplex cross-correlation of DNA sequences,โ€ Journal of Biological Systems, vol. 18, no. 4, pp. 711โ€“725, 2010. View at: Publisher Site | Google Scholar
  10. T. D. Schneider and D. N. Mastronarde, โ€œFast multiple alignment of ungapped DNA sequences using information theory and a relaxation method,โ€ Discrete Applied Mathematics, vol. 71, no. 1โ€“3, pp. 259โ€“268, 1996. View at: Publisher Site | Google Scholar
  11. T. D. Schneider, G. D. Stormo, L. Gold, and A. Ehrenfeucht, โ€œInformation content of binding sites on nucleotide sequences,โ€ Journal of Molecular Biology, vol. 188, no. 3, pp. 415โ€“431, 1986. View at: Google Scholar
  12. M. B. Toledano, I. Kullik, F. Trinh, P. T. Baird, T. D. Schneider, and G. Storz, โ€œRedox-dependent shift of OxyR-DNA contacts along an extended DNA-binding site: a mechanism for differential promoter selection,โ€ Cell, vol. 78, no. 5, pp. 897โ€“909, 1994. View at: Publisher Site | Google Scholar
  13. T. D. Schneider and R. M. Stephens, โ€œSequence logos: a new way to display consensus sequences,โ€ Nucleic Acids Research, vol. 18, no. 20, pp. 6097โ€“6100, 1990. View at: Google Scholar
  14. S. Delgado, D. Casane, L. Bonnaud, M. Laurin, J. Y. Sire, and M. Girondot, โ€œMolecular evidence for precambrian origin of amelogenin, the major protein of vertebrate enamel,โ€ Molecular Biology and Evolution, vol. 18, no. 12, pp. 2146โ€“2153, 2001. View at: Google Scholar
  15. J.-J. Shu, Q.-W. Wang, and K.-Y. Yong, โ€œDNA-based computing of strategic assignment problems,โ€ Physical Review Letters, vol. 106, no. 18, Article ID 188702, 2011. View at: Publisher Site | Google Scholar

Copyright ยฉ 2012 Jian-Jun Shu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Related articles

No related content is available yet for this article.
 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

No related content is available yet for this article.

Article of the Year Award: Outstanding research contributions of 2021, as selected by our Chief Editors. Read the winning articles.