title: Maximum-Likelihood Phylogenetic Inference with Selection on Protein Folding Stability.
abstract: Despite intense work, incorporating constraints on protein native structures into the mathematical models of molecular evolution remains difficult, because most models and programs assume that protein sites evolve independently, whereas protein stability is maintained by interactions between sites. Here, we address this problem by developing a new mean-field substitution model that generates independent site-specific amino acid distributions with constraints on the stability of the native state against both unfolding and misfolding. The model depends on a background distribution of amino acids and one selection parameter that we fix maximizing the likelihood of the observed protein sequence. The analytic solution of the model shows that the main determinant of the site-specific distributions is the number of native contacts of the site and that the most variable sites are those with an intermediate number of native contacts. The mean-field models obtained, taking into account misfolded conformations, yield larger likelihood than models that only consider the native state, because their average hydrophobicity is more realistic, and they produce on the average stable sequences for most proteins. We evaluated the mean-field model with respect to empirical substitution models on 12 test data sets of different protein families. In all cases, the observed site-specific sequence profiles presented smaller Kullback-Leibler divergence from the mean-field distributions than from the empirical substitution model. Next, we obtained substitution rates combining the mean-field frequencies with an empirical substitution model. The resulting mean-field substitution model assigns larger likelihood than the empirical model to all studied families when we consider sequences with identity larger than 0.35, plausibly a condition that enforces conservation of the native structure across the family. We found that the mean-field model performs better than other structurally constrained models with similar or higher complexity. With respect to the much more complex model recently developed by Bordner and Mittelmann, which takes into account pairwise terms in the amino acid distributions and also optimizes the exchangeability matrix, our model performed worse for data with small sequence divergence but better for data with larger sequence divergence. The mean-field model has been implemented into the computer program Prot_Evol that is freely available at http://ub.cbm.uam.es/software/Prot_Evol.php.
authors.affiliation: Department of Cell Biology and Immunology, Centro de Biología Molecular Severo Ochoa (CSIC-UAM), Universidad Autónoma de Madrid, Madrid, Spain.
authors.forename: Miguel
authors.initials: M
authors.lastname: Arenas
authors.affiliation: Department of Cell Biology and Immunology, Centro de Biología Molecular Severo Ochoa (CSIC-UAM), Universidad Autónoma de Madrid, Madrid, Spain.
authors.forename: Agustin
authors.initials: A
authors.lastname: Sánchez-Cobos
authors.affiliation: Department of Cell Biology and Immunology, Centro de Biología Molecular Severo Ochoa (CSIC-UAM), Universidad Autónoma de Madrid, Madrid, Spain ubastolla@cbm.csic.es.
authors.forename: Ugo
authors.initials: U
authors.lastname: Bastolla
mesh_terms: D008956:Models, Chemical; D008957:Models, Genetic; D010802:Phylogeny; D017510:Protein Folding; D055550:Protein Stability; D011506:Proteins / Q000737:chemistry* / Q000235:genetics*
chemical_list: D011506:Proteins
keywords: folding stability; maximum-likelihood estimate; misfolded state; structurally constrained substitution models
references.citation: Pac Symp Biocomput. 2001;:191-202
references.pmid: 11262940
references.citation: Mol Biol Evol. 2001 May;18(5):691-9
references.pmid: 11319253
references.citation: Mol Biol Evol. 2001 May;18(5):750-6
references.pmid: 11319259
references.citation: Trends Genet. 2001 May;17(5):262-72
references.pmid: 11335036
references.citation: Proteins. 2001 Aug 1;44(2):79-96
references.pmid: 11391771
references.citation: J Mol Biol. 2001 Sep 7;312(1):289-307
references.pmid: 11545603
references.citation: Proteins. 2002 Jan 1;46(1):105-9
references.pmid: 11746707
references.citation: Mol Biol Evol. 2002 Mar;19(3):352-6
references.pmid: 11861895
references.citation: Q Rev Biophys. 2002 Aug;35(3):205-86
references.pmid: 12599750
references.citation: Mol Biol Evol. 2003 Oct;20(10):1692-704
references.pmid: 12885968
references.citation: J Mol Evol. 2003;57 Suppl 1:S103-19
references.pmid: 15008407
references.citation: Mol Biol Evol. 2004 Jun;21(6):1095-109
references.pmid: 15014145
references.citation: J Mol Biol. 2004 Nov 5;343(5):1451-66
references.pmid: 15491623
references.citation: Biophys J. 1987 Dec;52(6):1083-5
references.pmid: 3427197
references.citation: Mol Biol Evol. 1987 Jul;4(4):406-25
references.pmid: 3447015
references.citation: Biophys Chem. 1989 Nov;34(3):187-99
references.pmid: 2611345
references.citation: Comput Appl Biosci. 1992 Jun;8(3):275-82
references.pmid: 1633570
references.citation: Mol Biol Evol. 1993 Nov;10(6):1396-401
references.pmid: 8277861
references.citation: Proc Natl Acad Sci U S A. 1995 Feb 28;92(5):1282-6
references.pmid: 7877968
references.citation: Proteins. 1995 Mar;21(3):167-95
references.pmid: 7784423
references.citation: Fold Des. 1997;2(5):261-9
references.pmid: 9261065
references.citation: Proteins. 1997 Dec;29(4):461-6
references.pmid: 9408943
references.citation: Proc Natl Acad Sci U S A. 1998 Apr 28;95(9):4976-81
references.pmid: 9560213
references.citation: Mol Biol Evol. 1998 Jul;15(7):910-7
references.pmid: 9656490
references.citation: Proteins. 1998 Aug 15;32(3):289-95
references.pmid: 9715905
references.citation: Mol Biol Evol. 1998 Dec;15(12):1600-11
references.pmid: 9866196
references.citation: Mol Biol Evol. 1999 Feb;16(2):173-9
references.pmid: 10028285
references.citation: J Theor Biol. 1999 Sep 7;200(1):49-64
references.pmid: 10479539
references.citation: Proc Natl Acad Sci U S A. 1999 Sep 14;96(19):10689-94
references.pmid: 10485887
references.citation: BMC Evol Biol. 2004 Oct 28;4:42
references.pmid: 15511291
references.citation: Mol Biol Evol. 2009 Oct;26(10):2387-95
references.pmid: 19597162
references.citation: Biochim Biophys Acta. 2010 Jun;1804(6):1231-64
references.pmid: 20117254
references.citation: PLoS Comput Biol. 2010 May;6(5):e1000767
references.pmid: 20463869
references.citation: Annu Rev Phys Chem. 2011;62:301-26
references.pmid: 21453060
references.citation: Proteins. 2011 May;79(5):1396-407
references.pmid: 21337623
references.citation: Nucleic Acids Res. 2012 Jan;40(Database issue):D290-301
references.pmid: 22127870
references.citation: BMC Evol Biol. 2011;11:361
references.pmid: 22171550
references.citation: Protein Sci. 2012 Jun;21(6):769-85
references.pmid: 22528593
references.citation: PLoS Comput Biol. 2012;8(6):e1002572
references.pmid: 22761562
references.citation: Mol Biol Evol. 2013 Apr;30(4):772-80
references.pmid: 23329690
references.citation: Proteins. 2013 Jul;81(7):1102-12
references.pmid: 23280507
references.citation: Bioinformatics. 2013 Dec 1;29(23):3020-8
references.pmid: 24037213
references.citation: Mol Biol Evol. 2014 Mar;31(3):736-49
references.pmid: 24307688
references.citation: BMC Evol Biol. 2014;14:78
references.pmid: 24716445
references.citation: Proteins. 2005 Jan 1;58(1):22-30
references.pmid: 15523667
references.citation: Proc Natl Acad Sci U S A. 2005 Jan 18;102(3):606-11
references.pmid: 15644440
references.citation: Mol Biol Evol. 2005 Mar;22(3):630-8
references.pmid: 15537801
references.citation: Gene. 2005 Mar 14;347(2):207-17
references.pmid: 15733531
references.citation: Proc Natl Acad Sci U S A. 2005 Jul 5;102(27):9541-6
references.pmid: 15980155
references.citation: Nat Rev Genet. 2005 Sep;6(9):678-87
references.pmid: 16074985
references.citation: Chem Rev. 2006 May;106(5):1559-88
references.pmid: 16683745
references.citation: BMC Evol Biol. 2006;6:43
references.pmid: 16737532
references.citation: Mol Biol Evol. 2007 Aug;24(8):1586-91
references.pmid: 17483113
references.citation: Proteins. 2008 Dec;73(4):872-88
references.pmid: 18536008