Email updates

Keep up to date with the latest news and content from BMC Bioinformatics and BioMed Central.

Open Access Research article

Detailed protein sequence alignment based on Spectral Similarity Score (SSS)

Kshitiz Gupta123*, Dina Thomas1, SV Vidya4, KV Venkatesh23* and S Ramakumar45

Author Affiliations

1 Department of Computer Science & Engineering, Indian Institute of Technology, Bombay, Mumbai, India

2 Department of Chemical Engineering, Indian Institute of Technology, Bombay, Mumbai, India

3 School of Biosciences & Bioengineering, Indian Institute of Technology, Bombay, Mumbai, India

4 Department of Physics, Indian Institute of Science, Bangalore, India

5 Bioinformatics Center, Indian Institute of Science, Bangalore, India

For all author emails, please log on.

BMC Bioinformatics 2005, 6:105  doi:10.1186/1471-2105-6-105

Published: 23 April 2005

Abstract

Background

The chemical property and biological function of a protein is a direct consequence of its primary structure. Several algorithms have been developed which determine alignment and similarity of primary protein sequences. However, character based similarity cannot provide insight into the structural aspects of a protein. We present a method based on spectral similarity to compare subsequences of amino acids that behave similarly but are not aligned well by considering amino acids as mere characters. This approach finds a similarity score between sequences based on any given attribute, like hydrophobicity of amino acids, on the basis of spectral information after partial conversion to the frequency domain.

Results

Distance matrices of various branches of the human kinome, that is the full complement of human kinases, were developed that matched the phylogenetic tree of the human kinome establishing the efficacy of the global alignment of the algorithm. PKCd and PKCe kinases share close biological properties and structural similarities but do not give high scores with character based alignments. Detailed comparison established close similarities between subsequences that do not have any significant character identity. We compared their known 3D structures to establish that the algorithm is able to pick subsequences that are not considered similar by character based matching algorithms but share structural similarities. Similarly many subsequences with low character identity were picked between xyna-theau and xyna-clotm F/10 xylanases. Comparison of 3D structures of the subsequences confirmed the claim of similarity in structure.

Conclusion

An algorithm is developed which is inspired by successful application of spectral similarity applied to music sequences. The method captures subsequences that do not align by traditional character based alignment tools but give rise to similar secondary and tertiary structures. The Spectral Similarity Score (SSS) is an extension to the conventional similarity methods and results indicate that it holds a strong potential for analysis of various biological sequences and structural variations in proteins.