A maximum likelihood framework for protein design Export

@article{citeulike:721485, abstract = {BACKGROUND:The aim of protein design is to predict amino-acid sequences compatible with a given target structure. Traditionally envisioned as a purely thermodynamic question, this problem can also be understood in a wider context, where additional constraints are captured by learning the sequence patterns displayed by natural proteins of known conformation. In this latter perspective, however, we still need a theoretical formalization of the question, leading to general and efficient learning methods, and allowing for the selection of fast and accurate objective functions quantifying sequence/structure compatibility.RESULTS:We propose a formulation of the protein design problem in terms of model-based statistical inference. Our framework uses the maximum likelihood principle to optimize the unknown parameters of a statistical potential, which we call an inverse potential to contrast with classical potentials used for structure prediction. We propose an implementation based on Markov chain Monte Carlo, in which the likelihood is maximized by gradient descent and is numerically estimated by thermodynamic integration. The fit of the models is evaluated by cross-validation. We apply this to a simple pairwise contact potential, supplemented with a solvent-accessibility term, and show that the resulting models have a better predictive power than currently available pairwise potentials. Furthermore, the model comparison method presented here allows one to measure the relative contribution of each component of the potential, and to choose the optimal number of accessibility classes, which turns out to be much higher than classically considered.CONCLUSION:Altogether, this reformulation makes it possible to test a wide diversity of models, using different forms of potentials, or accounting for other factors than just the constraint of thermodynamic stability. Ultimately, such model-based statistical analyses may help to understand the forces shaping protein sequences, and driving their evolution.}, author = {Kleinman, Claudia and Rodrigue, Nicolas and Bonnard, Cecile and Philippe, Herve and Lartillot, Nicolas}, citeulike-article-id = {721485}, citeulike-linkout-0 = {http://dx.doi.org/10.1186/1471-2105-7-326}, citeulike-linkout-1 = {http://view.ncbi.nlm.nih.gov/pubmed/16808841}, citeulike-linkout-2 = {http://www.hubmed.org/display.cgi?uids=16808841}, day = {29}, doi = {10.1186/1471-2105-7-326}, issn = {1471-2105}, journal = {BMC Bioinformatics}, keywords = {bioinformatics, clique, phylogeny, robert-cedergren-center, stuff, udem}, month = {June}, number = {1}, pages = {326+}, posted-at = {2009-11-07 23:15:13}, priority = {2}, title = {A maximum likelihood framework for protein design}, url = {http://dx.doi.org/10.1186/1471-2105-7-326}, volume = {7}, year = {2006} }

TY - JOUR ID - citeulike:721485 L3 - citeulike-article-id:721485 N2 - BACKGROUND:The aim of protein design is to predict amino-acid sequences compatible with a given target structure. Traditionally envisioned as a purely thermodynamic question, this problem can also be understood in a wider context, where additional constraints are captured by learning the sequence patterns displayed by natural proteins of known conformation. In this latter perspective, however, we still need a theoretical formalization of the question, leading to general and efficient learning methods, and allowing for the selection of fast and accurate objective functions quantifying sequence/structure compatibility.RESULTS:We propose a formulation of the protein design problem in terms of model-based statistical inference. Our framework uses the maximum likelihood principle to optimize the unknown parameters of a statistical potential, which we call an inverse potential to contrast with classical potentials used for structure prediction. We propose an implementation based on Markov chain Monte Carlo, in which the likelihood is maximized by gradient descent and is numerically estimated by thermodynamic integration. The fit of the models is evaluated by cross-validation. We apply this to a simple pairwise contact potential, supplemented with a solvent-accessibility term, and show that the resulting models have a better predictive power than currently available pairwise potentials. Furthermore, the model comparison method presented here allows one to measure the relative contribution of each component of the potential, and to choose the optimal number of accessibility classes, which turns out to be much higher than classically considered.CONCLUSION:Altogether, this reformulation makes it possible to test a wide diversity of models, using different forms of potentials, or accounting for other factors than just the constraint of thermodynamic stability. Ultimately, such model-based statistical analyses may help to understand the forces shaping protein sequences, and driving their evolution. IS - 1 JF - BMC Bioinformatics SN - 1471-2105 TI - A maximum likelihood framework for protein design VL - 7 SP - 326 KW - bioinformatics KW - clique KW - phylogeny KW - robert-cedergren-center KW - stuff KW - udem AU - Kleinman, Claudia AU - Rodrigue, Nicolas AU - Bonnard, Cecile AU - Philippe, Herve AU - Lartillot, Nicolas PY - 2006/06/29/ UR - http://dx.doi.org/10.1186/1471-2105-7-326 ER -