Ottenere sequenze proteiche accedendo a Uniprot (con Python)
Sep 21 2020
Ho un elenco di ID di proteine che sto cercando di accedere alle sequenze di proteine da Uniprot con Python. Mi sono imbattuto in questo post: Sequenza proteica da uniprot protein id python ma fornisce un elenco di elementi e non la sequenza effettiva:
Codice
import requests as r
from Bio import SeqIO
from io import StringIO
cID='P04637'
baseUrl="http://www.uniprot.org/uniprot/"
currentUrl=baseUrl+cID+".fasta"
response = r.post(currentUrl)
cData=''.join(response.text)
Seq=StringIO(cData)
pSeq=list(SeqIO.parse(Seq,'fasta'))
che dà output:
produzione
[SeqRecord(seq=Seq('MQAALIGLNFPLQRRFLSGVLTTTSSAKRCYSGDTGKPYDCTSAEHKKELEECY...SSS', SingleLetterAlphabet()), id='sp|O45228|PROD_CAEEL', name='sp|O45228|PROD_CAEEL', description='sp|O45228|PROD_CAEEL Proline dehydrogenase 1, mitochondrial OS=Caenorhabditis elegans OX=6239 GN=prdh-1 PE=2 SV=2', dbxrefs=[])]
Ero solo curioso di sapere come posso effettivamente ottenere la sequenza stessa.
Risposte
3 Pallie Sep 21 2020 at 20:06
[record.seq for record in pSeq]
modifica: ti consigliamo str(pSeq[0].seq)