Hidden Markov models that use predicted local structure for fold recognition: alphabets of backbone geometry.
An important problem in computational biology is predicting the structure of the large number of putative proteins discovered by genome sequencing projects. Fold-recognition methods attempt to solve the problem by relating the target proteins to known structures, searching for template proteins homologous to the target. Remote homologs that may have significant structural similarity are often not detectable by sequence similarities alone. To address this, we incorporated predicted local structure, a generalization of secondary structure, into two-track profile hidden Markov models (HMMs). We did not rely on a simple helix-strand-coil definition of secondary structure, but experimented with a variety of local structure descriptions, following a principled protocol to establish which descriptions are most useful for improving fold recognition and alignment quality. On a test set of 1298 nonhomologous proteins, HMMs incorporating a 3-letter STRIDE alphabet improved fold recognition accuracy by 15% over amino-acid-only HMMs and 23% over PSI-BLAST, measured by ROC-65 numbers. We compared two-track HMMs to amino-acid-only HMMs on a difficult alignment test set of 200 protein pairs (structurally similar with 3-24% sequence identity). HMMs with a 6-letter STRIDE secondary track improved alignment quality by 62%, relative to DALI structural alignments, while HMMs with an STR track (an expanded DSSP alphabet that subdivides strands into six states) improved by 40% relative to CE.
Karchin
Rachel
R
Center for Biomolecular Science and Engineering, Baskin School of Engineering, University of California, Santa Cruz 95064, USA. rachelk@soe.ucsc.edu
Cline
Melissa
M
Mandel-Gutfreund
Yael
Y
Karplus
Kevin
K
D000465:Algorithms; D019295:Computational Biology / Q000379:methods; D008390:Markov Chains; D017510:Protein Folding; D017433:Protein Structure, Secondary; D011506:Proteins / Q000737:chemistry*; D015203:Reproducibility of Results; D016415:Sequence Alignment / Q000379:methods
D011506:Proteins