Home > Community > DMPfold error while running protein structure prediction
Upvote

25

Downvote
+ Bioinformatics
+ Biochemistry
Posted by
Melissa Zimmerman

DMPfold error while running protein structure prediction

Nan.s  Follow

Thanks, I did this but got this erro:

blastpgp -a 12 -b 0 -j 3 -h 0.001 -d c600m2-tail-before-lysin-sequence.a3m -i c600m2-tail-before-lysin-sequence.temp.fasta[blastpgp] WARNING: Unable to open BLOSUM62[blastpgp] WARNING: BlastScoreBlkMatFill returned non-zero status[blastpgp] WARNING: SetUpBlastSearch failed.BLASTP 2.2.26 [Sep-21-2011]Reference: Altschul, Stephen F., Thomas L. Madden, Alejandro A. Schaffer, Jinghui Zhang, Zheng Zhang, Webb Miller, and David J. Lipman (1997), "Gapped BLAST and PSI-BLAST: a new generation of protein database searchprograms",  Nucleic Acids Res. 25:3389-3402.Reference for compositional score matrix adjustment: Altschul, Stephen F., John C. Wootton, E. Michael Gertz, Richa Agarwala, Aleksandr Morgulis,Alejandro A. Schaffer, and Yi-Kuo Yu (2005) "Protein database searchesusing compositionally adjusted substitution matrices", FEBS J. 272:5101-5109.Reference for composition-based statistics starting in round 2:Schaffer, Alejandro A., L. Aravind, Thomas L. Madden,Sergei Shavirin, John L. Spouge, Yuri I. Wolf,  Eugene V. Koonin, and Stephen F. Altschul (2001), "Improving the accuracy of PSI-BLAST protein database searches with composition-based statistics and other refinements",  Nucleic Acids Res. 29:2994-3005.Query= C600M2_00100 tail protein [Salmonella phage BPS17S6]         (377 letters)Database: c600m2-tail-before-lysin-sequence.a3m            46,978 sequences; 11,323,713 total lettersSearchingSegmentation fault (core dumped)

Could you please let me know how this could be solved? Thanks

More

Upvote

VOTE

Downvote
Jake Missy  Follow
Ok, I got it resolved by changing the correct blastpgp (legacy blast) executable. Thanks!More
Upvote

VOTE

Downvote
Andy Fletcher  Follow

The comment by Devon Ryan is correct - looks like your blastpgp command is not working properly. For a protein of that length it's more likely to be a problem with the BLAST executable running on your system than a memory issue.

Try running the blastpgp command (https://github.com/psipred/DMPfold/blob/master/bin/runpsipredandsolvwithdb#L26) on its own and see what happens:

~/path/to/blastpgp -a 12 -b 0 -j 3 -h 0.001 -d yourfile.a3m -i yourfile.fasta

More

Upvote

VOTE

Downvote