Applying A Pssm To Psi-Blast
1
1
Entering edit mode
11.9 years ago

I'm able to generate a PSSM and a sequence from combinatorial experimental data. I want to BLAST the results to find proteins displaying these motifs. However, I have information about each position's sensitivity to substitutions, in addition to a consensus sequence. I know how to create a PSSM from this information, but I don't know how to incorporate this into a BLAST search. Does anyone know a cmdline argument to blast or a way of hacking a PSSM into a PSI-BLAST search to add position information to a sequence?

pssm blast • 3.7k views
ADD COMMENT
0
Entering edit mode
11.9 years ago

You might find some information in this post : http://www.biostars.org/post/show/1699/how-to-scale-a-pssm-for-ncbi-blast/

As well as here : http://etutorials.org/Misc/blast/Part+V+BLAST+Reference/Chapter+13.+NCBI-BLAST+Reference/13.8+blastpgp+Parameters+PSI-BLAST+and+PHI-BLAST/ and here http://gogarten.uconn.edu/MCB372/Laboratories/assign3.html

What you probably need to do:

  • generate a .chk checkpoint file similar to the one generated by 'blastpgp' containing your PSSM
  • feed it to blastpgp with something similar to:

    blastpgp -i test1.fa -d genome.faa -R test1.chk -o genome.test1.br -I T -a 2

ADD COMMENT

Login before adding your answer.

Traffic: 1868 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6