title: Bayesian mixture modeling using a hybrid sampler with application to protein subfamily identification.
abstract: Predicting protein function is essential to advancing our knowledge of biological processes. This article is focused on discovering the functional diversification within a protein family. A Bayesian mixture approach is proposed to model a protein family as a mixture of profile hidden Markov models. For a given mixture size, a hybrid Markov chain Monte Carlo sampler comprising both Gibbs sampling steps and hierarchical clustering-based split/merge proposals is used to obtain posterior inference. Inference for mixture size concentrates on comparing the integrated likelihoods. The choice of priors is critical with respect to the performance of the procedure. Through simulation studies, we show that 2 priors that are based on independent data sets allow correct identification of the mixture size, both when the data are homogeneous and when the data are generated from a mixture. We illustrate our method using 2 sets of real protein sequences.
authors.affiliation: Department of Biostatistics, University of Washington, Seattle, WA 98105, USA.
authors.forename: Youyi
authors.initials: Y
authors.lastname: Fong
authors.forename: Jon
authors.initials: J
authors.lastname: Wakefield
authors.forename: Kenneth
authors.initials: K
authors.lastname: Rice
mesh_terms: D000465:Algorithms; D000595:Amino Acid Sequence / Q000235:genetics; D000818:Animals; D001499:Bayes Theorem; D016000:Cluster Analysis; D003198:Computer Simulation; D006801:Humans; D016013:Likelihood Functions; D008390:Markov Chains; D015233:Models, Statistical; D009010:Monte Carlo Method; D011496:Protein Methyltransferases / Q000145:classification / Q000235:genetics; D011506:Proteins / Q000145:classification* / Q000235:genetics; D011968:Receptors, Histamine / Q000145:classification / Q000235:genetics; D016415:Sequence Alignment; D017386:Sequence Homology, Amino Acid
chemical_list: D011506:Proteins; D011968:Receptors, Histamine; D011496:Protein Methyltransferases