1 / 36

Remote Homology Detection of Beta-Structural Motifs Using Random Fields

Remote Homology Detection of Beta-Structural Motifs Using Random Fields. Matt Menke, Tufts Bonnie Berger, MIT Lenore Cowen, Tufts ISMB 3Dsig 2010 July 10, 2010. Inferring structural similarity from homology is hard at the SCOP superfamily/fold level. Profile HMMs.

tory
Download Presentation

Remote Homology Detection of Beta-Structural Motifs Using Random Fields

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Remote Homology Detection of Beta-Structural Motifs Using Random Fields Matt Menke, Tufts Bonnie Berger, MIT Lenore Cowen, Tufts ISMB 3Dsig 2010 July 10, 2010

  2. Inferring structural similarity from homology is hard at the SCOP superfamily/fold level

  3. Profile HMMs

  4. HMM is trained from Sequence Alignment of Known Structures But: cannot capture pariwise long-range beta-sheet interactions!

  5. HMMs cannot capture statistical preferences from residues close in space but far, and a variable distance apart in seq. Pectate Lyase C (Yoder et al. 1993)

  6. Only look at Pairs: Generalize to Markov Random Fields Liu et al. 2009 Zhao et al. 2010 Menke et al. 2010 (This work) B3 T2 B2 B1 Look at Just Pairs or Generalize to Markov Random Fields [Bradley, Cowen, Menke, King, Berger, PNAS, 2001, 98:26, 14,819-14,824; Cowen, Bradley, Menke, King, Berger (2002), J Comp Biol, 9, 261-276]

  7. Let’s look at what this would mean for propeller folds

  8. Goal: capture HMM sequence information and pairwise information in beta-structural motifs at the same time! SCOP (http://scop.mrc-lmb.cam.ac.uk/scop

  9. Structural Motifs Using Random Fields SMURF

  10. Structural Motifs Using Random Fields Can we get the benefit of pairwise correlations without having to throw away all sequence info?

  11. The template is learned from solved structures in the PDB

  12. The template is learned from solved structures in the PDB:Aligned with Matt

  13. Digression: Matt structural alignment program Menke, Berger, Cowen, (PLOS Combio 2008) Specifically designed to align more distant homologs AFP chaining using dynamic programming with “translations and twists” (flexibility)

  14. The template is learned from solved structures in the PDB:Aligned with Matt

  15. Two beta tables are learned from amphapathic beta sheets that are not propellers from solved structures in the PDB. Two pairwise Exposed Residue Buried Residue http://bcb.cs.tufts.edu/propellers/si/

  16. Computing a Score • Sequences are scored by computing their best “threading” or “parse” against the template as a sum of HMM(score) + pairwise(score) • No longer polynomial time (multi-dimensional dynamic programming) • Tractable on propellers because paired beta-strands don’t interleave too much

  17. Let’s look at what this would mean for propeller folds

  18. Let’s look at what this would mean for propeller folds • Training set for HMM score: leave-superfamily-out cross validation • Training set for pairwise score: amphapathic beta-sheets from NON-propellers

  19. Results on Propellers

  20. Results on Propellers • Note that this is “6 (or 7)” bladed propeller versus non-propeller– distinguishing the number of blades in the propeller seems to be a much harder problem….

  21. Different propeller closures 1jof 2trc

  22. So: what new sequences fold into propellers? • We predict a double propeller motif in the N-terminal region of a hybrid 2-component sensor protein.

  23. What are these proteins? • First found in a benign bacteria in human gut. • May be involved in adapting to changes in diet/efficiently processing different sugars • Found in other bacterial species: help sense and adapt to environmental changes. • Big stretch (I am not a biologist): help to study human obesity epidemic??

  24. Popular Domains • HisKA histidine kinase domain • GGDEF adenylyl cyclase signalling domain • SpoIIE sporulation domain • Gaf domain • PAS domain • HATPase domain

  25. Species distribution

  26. Distinguishing Number of Blades • The automatic SMURF consensus 7-bladed template only learns 6 blades. • Sequence motifs are similar– the same Pfam motif occurs in propellers with different numbers of blades • The fix: throw out propellers with a “funky” 7th blade by hand and build a new template. Now 6-bladed propellers don’t like the 7-bladed template • Double propellers we found are probably 7-7 (but 7-6 is also plausible).

  27. Predict propellers with Smurf! • http://smurf.cs.tufts.edu • Accepts sequences in FASTA format • 6,7,8-bladed templates, as well as all 9 double-propeller template http://bcb.cs.tufts.edu/propellers/si pairwise tables long list of predicted propeller sequences

  28. What’s Next for SMURF? Long-range dependencies Deeply interleaved β-strand pairs

  29. Conclusions • Combining an HMM score with a pairwise score can help recognize beta-structures • Computing this score exactly with a random field is highly computationally intensive • We will begin to look at when it is feasible and when we should use heuristics. • Also: add side-chain packing, other model refinements.

  30. More Questions • When should we over-weight the HMM versus the pair portion of the score? -- the case of 8-bladed propellers • Are there other ways to incorporate pairwise dependencies into HMMs?

  31. An Hmm is only as good as its training data • An Hmm is only as good as its training data– or is it? • Idea: we augment the training set, using the simplest model of evolution! • See Kumar and Cowen’s ISMB proceedings paper!

  32. Acknowledgements • National Institutes of Health Thank you!

More Related