Question: How Can I Convert Hmmer Hmm Descriptions To Psi-Blast Pssms?
gravatar for Alex Richter
9.2 years ago by
Alex Richter210
San Diego
Alex Richter210 wrote:

I'm attempting to generate an rpsblast database from a set of protein cluster alignments. For most of them, I could use blastpgp -B <alignment> to generate the PSSM without problem, but several of the alignments (namely the large ones, or the ones with huge sequences) crash (seg-fault) blastpgp (legacy v2.2.23, if it matters).

However, these alignments can be hmmbuilt into HMMER profiles without complaint. I know there has to be a way to convert the HMMs to psiblast profiles, since does it for CDD with PFAM and TIGRFAM. I'm just hoping it something better than exporting the alignment, then blastpgp -- particularly since that's what's crashing.

Anyone got a script to convert them directly?

hmm • 2.9k views
ADD COMMENTlink written 9.2 years ago by Alex Richter210

Was the script working or were you able to create your own conversion script? If so, please share your findings because I'm looking into the exact same thing right now.

ADD REPLYlink written 9.1 years ago by Michael Schubert6.9k

Can someone explain the script of Yun He? I've never used Perl and I don't understand where I have to put what. If I change the

&read_hmmprof("", \@hmmprof);

to my own hmm profile it only prints the effective alphabet and nothing else.

ADD REPLYlink modified 11 weeks ago by RamRS24k • written 9.0 years ago by Niek De Klein2.5k
gravatar for iw9oel_ad
9.2 years ago by
iw9oel_ad6.0k wrote:

This post by Yun He from 2005 on the Sequence Search mailing list offers a method via hmmconvert and a script. I haven't tried it myself, nor have I tried to check to see if it's correct.

ADD COMMENTlink written 9.2 years ago by iw9oel_ad6.0k

Looks interesting, although it's missing a couple columns in the PSSM, and I'm not entirely sure where the $your_query_sequence is coming from. I'll have to see if I can get it to work. If nothing else, it'll give me a place to start.

ADD REPLYlink modified 11 weeks ago by RamRS24k • written 9.2 years ago by Alex Richter210
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: 1821 users visited in the last hour