This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Obtaining SNP locations on protein sequences.

Hello. I'm a researcher in the field in proteomics. Currently, I'm trying to obtain SNP locations on protein sequences, along with protein sequences.

My question is> Which database/service currently hosts protein SNP location, where I can batch download data for every protein? Is this even possible?

Thank you very much.

snp protein-sequence
• 2,870 views • link
updated by Ram Ram • written by Bioaln Bioaln

2 answers

This information is available in uniprot ftp://ftp.uniprot.org/pub/databases/uniprot/current_release/knowledgebase/complete/uniprot_sprot.xml.gz

(...)
  <feature description="In dbSNP:rs11542705." id="VAR_048095" type="sequence variant">
    <original>M</original>
    <variation>I</variation>
    <location>
      <position position="155"/>
    </location>
  </feature>
(...)
• 0 views • link
updated by Ram Ram • written by Pierre Lindenbaum Pierre Lindenbaum

Sorry for late reply, I tried to parse this with python, and I successfully get all of the variants. Is there any way I can obtain only nsSNPs? Or ones, connected with pathogenic effects?

• 0 views • link
written by Bioaln Bioaln

You can use a tool we have recently published at bioinformatics: I-PV. It will print your protein sequence along with SNVs, their polyphen and sift scores, Indels, aminoacid sequence, their chemical properties and corresponding codons. You will be able to see possible point mutations at each location and distribution of a set of amino acids to another set of amino acids. I have uploaded a set of introductory videos at the I-PV's website. Here is one of them: http://i-pv.org/intro_ipv_alt4.html

You will need the fasta files of your mRNA (NM_...) and protein sequence. You will also need a text file of conservation scores separated by newline character. (You can upload a dummy conservation file of random numbers if you like). Lastly, you will need the variant file for your SNVs where you can download from Biomart for your protein of interest. Or alternatively you can use a vcf file.

The resulting image will be interactive and you can still plot/hide data on it using the highlight tool or drop down menus. To have an idea what the output looks like and whether if it fits what you want take a look at some examples:

I hope this helps,

Good luck with your research,

• 0 views • link
updated by Ram Ram • written by Ibrahim Tanyalcin Ibrahim Tanyalcin

Thanks for the answer, were those made with Circos?

• 0 views • link
written by Bioaln Bioaln

Dear Bioaln,

The software is built on top of circos correct. It is a combination of circos and javascript. However, you do not have to generate datatracks yourself, they are automatically generated from the fasta files you provide. The output will open in a browser, and when you click on the SNPs, it will take you to the corresponding page for further information at dbSNP, like this example:

http://i-pv.org/gifs/snpToDbsnp.gif

I hope this helps,

• 0 views • link
written by Ibrahim Tanyalcin Ibrahim Tanyalcin

Log in to answer this question.

More posts like this