The compositional adjustment of amino acid substitution matrices.
Amino acid substitution matrices are central to protein-comparison methods. In most commonly used matrices, the substitution scores take a log-odds form, involving the ratio of "target" to "background" frequencies derived from large, carefully curated sets of protein alignments. However, such matrices often are used to compare protein sequences with amino acid compositions that differ markedly from the background frequencies used for the construction of the matrices. Of course, the target frequencies should be adjusted in such cases, but the lack of an appropriate way to do this has been a long-standing problem. This article shows that if one demands consistency between target and background frequencies, then a log-odds substitution matrix implies a unique set of target and background frequencies as well as a unique scale. Standard substitution matrices therefore are truly appropriate only for the comparison of proteins with standard amino acid composition. Accordingly, we present and evaluate a rationale for transforming the target frequencies implicit in a standard matrix to frequencies appropriate for a nonstandard context. This rationale yields asymmetric matrices for the comparison of proteins with divergent compositions. Earlier approaches are unable to deal with this case in a fully consistent manner. Composition-specific substitution matrix adjustment is shown to be of utility for comparing compositionally biased proteins, including those of organisms with nucleotide-biased, and therefore codon-biased, genomes or isochores.
Yu
Yi-Kuo
YK
National Center for Biotechnology Information, National Library of Medicine, National Institutes of Health, Bethesda, MD 20894, USA.
Wootton
John C
JC
Altschul
Stephen F
SF
D000595:Amino Acid Sequence; D019943:Amino Acid Substitution / Q000235:genetics*; D000818:Animals; D001217:Aspartate-Ammonia Ligase / Q000737:chemistry / Q000235:genetics; D019143:Evolution, Molecular; D005787:Gene Frequency; D008969:Molecular Sequence Data; D020125:Mutation, Missense; D009169:Mycobacterium tuberculosis / Q000235:genetics; D010963:Plasmodium falciparum / Q000235:genetics; D011506:Proteins / Q000235:genetics; D015203:Reproducibility of Results; D016415:Sequence Alignment; D017386:Sequence Homology, Amino Acid
D011506:Proteins; D001217:Aspartate-Ammonia Ligase