Question: Creating An Rpsblast Database From An Hmm Profile Library
gravatar for Michael Schubert
10.2 years ago by
Cambridge, UK
Michael Schubert7.0k wrote:

I have downloaded the HMM library file from CATH and would like to convert it to a rpsblast database.

Therefor, I think I need to do the following steps:

  1. Extract the HMM profiles from the library (this seems to be a simple split at //)
  2. Convert each profile into a PSSM (with something like this)
  3. Use formatrpsdb to create a database from the PSSMs

Is there a smarter way to do is? I do not want to use hmmsearch directly because I would like to integrate the CATHdb into an already exisiting rpsblast framework.

blast conversion hmmer • 2.5k views
ADD COMMENTlink modified 10.2 years ago by Neilfws49k • written 10.2 years ago by Michael Schubert7.0k
gravatar for Neilfws
10.2 years ago by
Sydney, Australia
Neilfws49k wrote:

Your approach sounds reasonable. Another solution is outlined in this HMMER tutorial - scroll down to "Using an HMM model to create an alignment". It is for HMMER 2.3.2; I don't know if CATH has moved to HMMER3 to build profiles?

Basically, they suggest using hmmalign on a set of sequences retrieved from hmmsearch. The resulting alignment can then be used to build a PSSM. However, your solution might be quicker, if you can go straight from HMMER profile to PSSM.

ADD COMMENTlink written 10.2 years ago by Neilfws49k
Please log in to add an answer.


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