Question: Finding Zinc Finger Domains Using Hmmer And Pfam
3
gravatar for justin
8.2 years ago by
justin30
justin30 wrote:

Reading this: http://selab.janelia.org/people/eddys/blog/?p=162.

I'm specifically researching zinc fingers.

Am I right in saying, I should either use HMMER 3.0 with Pfam 26 assuming Rob has rebuilt the seed alignments for the problematic families or use HMMER 2.3.2 with Pfam 23? (Before Pfam Cambridge used HMMER 3.0 to produce their profile HMMs)

What I'm really trying to do is run a hmmscan on protein sequences against pfam to identify all the domains. Then extract the proteins that have zinc finger domains in them.

Any other comments about this would be greatly appreciated.

Cheers

hmmer • 3.2k views
ADD COMMENTlink modified 6.9 years ago by Biostar ♦♦ 20 • written 8.2 years ago by justin30
3
gravatar for robert.d.finn
8.2 years ago by
robert.d.finn130
robert.d.finn130 wrote:

Hi,

I guess as I am specifically named in this post I should reply. So, although the curators at Pfam have added many additional Zn-finger families to Pfam, the current version of HMMER (3.0) still struggles with finding all of these short, relatively divergent domains. I can confirm that Sean Eddy, the author of HMMER, is actively working on fixing the problem. The problem arises due to the fact that HMMER 3.0 only produces local-local alignments (i.e. alignments can be generated from a partial match the model and a sub region of the sequence). In HMMER 2.3.2, Zn-fingers and such like, where detected via the global-local mode (i.e. the full length of the model had to be matched against the sequence or sub region of sequence). Forcing the match along the entire length of the model, provided just enough of a boost to the score, to enable clear distinction of the motif from the background or noise. In HMMER 3.1, glocal alignments will return and the problem will go away. From the initial results I have seen, it will even improve the detection rate of such motifs.

In the meantime, what to do? Well, if this were me and I was specifically interested in Zn-fingers, I would probably run both Pfam 23 and 26 with the respective matches and take the union of the hits. You may also consider running the sequences against InterPro as well. The different methods/signatures may detect additional matches.

I hope this helps,

Rob

ADD COMMENTlink written 8.2 years ago by robert.d.finn130

oh wow! didn't expect that, Awesome! Thanks for the clarification! Yeah I was planning to use InterPro as well :D

ADD REPLYlink modified 8.2 years ago • written 8.2 years ago by justin30

As of Hmmer4 and Nov 2017, this is still a problem for zinc fingers.

ADD REPLYlink written 3.0 years ago by Maximilian Haeussler1.3k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1138 users visited in the last hour