title: Accuracy estimation and parameter advising for protein multiple sequence alignment.
abstract: Abstract We develop a novel and general approach to estimating the accuracy of multiple sequence alignments without knowledge of a reference alignment, and use our approach to address a new task that we call parameter advising: the problem of choosing values for alignment scoring function parameters from a given set of choices to maximize the accuracy of a computed alignment. For protein alignments, we consider twelve independent features that contribute to a quality alignment. An accuracy estimator is learned that is a polynomial function of these features; its coefficients are determined by minimizing its error with respect to true accuracy using mathematical optimization. Compared to prior approaches for estimating accuracy, our new approach (a) introduces novel feature functions that measure nonlocal properties of an alignment yet are fast to evaluate, (b) considers more general classes of estimators beyond linear combinations of features, and (c) develops new regression formulations for learning an estimator from examples; in addition, for parameter advising, we (d) determine the optimal parameter set of a given cardinality, which specifies the best parameter values from which to choose. Our estimator, which we call Facet (for "feature-based accuracy estimator"), yields a parameter advisor that on the hardest benchmarks provides more than a 27% improvement in accuracy over the best default parameter choice, and for parameter advising significantly outperforms the best prior approaches to assessing alignment quality.
authors.affiliation: Department of Computer Science, University of Arizona, Tucson, AZ 85721, USA. kece@cs.arizona.edu
authors.forename: John
authors.initials: J
authors.lastname: Kececioglu
authors.forename: Dan
authors.initials: D
authors.lastname: DeBlasio
mesh_terms: D000595:Amino Acid Sequence; D030562:Databases, Protein; D017433:Protein Structure, Secondary; D011506:Proteins / Q000737:chemistry*; D016415:Sequence Alignment / Q000379:methods*
chemical_list: D011506:Proteins