Hidden Markov models that use predicted local structure for fold recognition: Alphabets of backbone geometry

Citations Over TimeTop 10% of 2003 papers

Abstract

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.

Related Papers

• → Homology-extended sequence alignment(2005)119 cited
• → Fast structure alignment for protein databank searching(1992)138 cited
• → Optimization of multiple‐sequence alignment based on multiple‐structure alignment(2005)52 cited
• → ALIGN_MTX—An optimal pairwise textual sequence alignment program, adapted for using in sequence-structure alignment(2009)10 cited
• → Using Structure to Explore the Sequence Alignment Space of Remote Homologs(2011)11 cited