Get results of keyword search on Pfam via python script

Question

I'm interested in all proteins that are in any way associated with Danio rerio. I decided to look them up at Pfam data base and when I just make a keyword search, I get a a nice list which looks like this http://pfam.xfam.org/search/keyword?query=danio+rerio&submit=Submit. Important thing to note is that most of the associations of these protein families can be found in seq_info column.

My question is: how to extract that list from a python script? Moreover, I want to make a script that when run, will fetch the latest release, filter out the proteins that cannot be found in D. rerio and do some additional processing.

I've been searching for ways to directly download files from Pfam ftp server and mine for the proteins that are associated with D. rerio, however, I didn't find a single way to do it. I simply cannot see where do they store the information about species in which the protein can be found.

I would like to retrieve everything about protein families of interest, because the next step I will need to do is to find the DNA binding proteins and for finding a way to filter out the ones of no interest, I would liketo have as much data as possible.

Answer 1

If I understand your question correctly, then you can use the Pfam proteom search to get the data you are looking for.

Go to the

Pfam website -> Browse -> select the letter D from proteoms category -> D. rerio

There you can find a download link under the section Domain composition:

You can download a list of all regions for this proteome from our FTP site.

The download link seems to be bound to the species ID and if you want to automate your data retrieval you may be able to simply download the file from there in your python script.