Reduced amino acid alphabets exhibit an improved sensitivity and selectivity in fold assignment.
BACKGROUND
Many proteins with vastly dissimilar sequences are found to share a common fold, as evidenced in the wealth of structures now available in the Protein Data Bank. One idea that has found success in various applications is the concept of a reduced amino acid alphabet, wherein similar amino acids are clustered together. Given the structural similarity exhibited by many apparently dissimilar sequences, we undertook this study looking for improvements in fold recognition by comparing protein sequences written in a reduced alphabet.


RESULTS
We tested over 150 of the amino acid clustering schemes proposed in the literature with all-versus-all pairwise sequence alignments of sequences in the Distance mAtrix aLIgnment database. We combined several metrics from information retrieval popular in the literature: mean precision, area under the Receiver Operating Characteristic curve and recall at a fixed error rate and found that, in contrast to previous work, reduced alphabets in many cases outperform full alphabets. We find that reduced alphabets can perform at a level comparable to full alphabets in correct pairwise alignment of sequences and can show increased sensitivity to pairs of sequences with structural similarity but low-sequence identity. Based on these results, we hypothesize that reduced alphabets may also show performance gains with more sophisticated methods such as profile and pattern searches.


BACKGROUND
A table of results as well as the substitution matrices and residue groupings from this study can be downloaded from (http://www.rpgroup.caltech.edu/publications/supplements/alphabets).
Peterson
Eric L
EL
Department of Physics, California Institute of Technology, Pasadena, CA 91125, USA.
Kondev
Jané
J
Theriot
Julie A
JA
Phillips
Rob
R
D000465:Algorithms; D000595:Amino Acid Sequence; D000596:Amino Acids / Q000737:chemistry*; D019295:Computational Biology / Q000379:methods*; D030562:Databases, Protein; D011487:Protein Conformation; D017510:Protein Folding; D016415:Sequence Alignment / Q000379:methods*; D020539:Sequence Analysis, Protein / Q000379:methods*
D000596:Amino Acids