Markov chain Monte Carlo computation of confidence intervals for substitution-rate variation in proteins.
We suggest a method implemented in a computer program, immodestly dubbed TSUNAMI, that allows us to compare two homologous protein subfamilies with respect to the distribution of substitution rates along sequences. This study furthers our earlier work on a wavelet model of rate variation (1). The current approach allows sensitive detection of subtle discordances in the selection patterns between two protein subfamilies. In addition to performing fast computation of the maximum posterior probability estimates of the relative substitution rates, the method can select the most appropriate number of wavelet parameters for a particular dataset. TSUNAMI is based on a Markov chain Monte Carlo technique, and appears to be more applicable to larger datasets than is the full likelihood-based approach.
Columbia Genome Center, Columbia University, 1050 St. Nicholas Avenue, Unit 109, New York, NY 10032, USA. ar345@columbia.edu
A
A
Rzhetsky
P
P
Morozov
D019943:Amino Acid Substitution; D003198:Computer Simulation; D016001:Confidence Intervals; D003627:Data Interpretation, Statistical; D019143:Evolution, Molecular; D016013:Likelihood Functions; D008390:Markov Chains; D008957:Models, Genetic; D009010:Monte Carlo Method; D011506:Proteins / Q000235:genetics*; D012984:Software
D011506:Proteins