[BiO BB] mapping swissprot sequence positions to pdb sequence positions
Iddo
idoerg at burnham.org
Thu Jul 15 20:02:01 EDT 2004
Sounds pretty good. Going through the hssp seqno would be my choice too.
Be aware that the gaps are marked with a "!", and I _think_ those "!"s
are numbered too. (To lazy to look at an HSSP file now).
Good luck,
Iddo
rich wrote:
>Hi, yes, thanks. I was thinking about using Align. I am
>using hssp entries which have a consecutive seqno
>position. I want to map these values onto the position
>within the swissprot entry for the protein sequence. I
>realise that there will be gaps etc in the pdb num but
>there won't be in the hssp seqno values.
>
>effectively i have snps and their positions and i want
>to map onto the seqno value from the hssp file. The
>seqno value is simply a consecutive seq position as
>opposed to the pdb seq position.
>
>I think the solution is to run Align on all of my
>sequences and then parse with Bio::Align, unless there
>is a quicker way that I haven't thought of
>
>cheers
>Rich
>
>thanks
>Rich
>
>Quoting Iddo Friedberg <idoerg at burnham.org>:
>
>
>
>>A pairwise alignment program, such as blast2seq
>>or EMBOSS's "needle" is
>>probably good for you.
>>
>>I'm not sure I understand your question, is the
>>sequence position coming
>>from the SwissProt file or from the PDB? If the
>>latter, be aware that:
>>
>>1) PDB sequences contain gaps.
>>2) Sequence numbers in PDB are non-consecutive
>>even if there are no gaps.
>>3) Sequence numbering can be alphanumeric and
>>contain insertion
>>codes(e.g. 183, 183A, 183B, 184)
>>4) Sequecnes can start with negative numbers
>>(e.g. for his-tagged proteins)
>>5) Mutations abound.
>>
>>Iddo
>>
>>rich wrote:
>>
>>
>>
>>>Hi,
>>>
>>>could anyone recommend a quick way of mapping
>>>
>>>
>>the
>>
>>
>>>positions in say a swissprot peptide sequence
>>>
>>>
>>file onto
>>
>>
>>>the pdb sequence file? I want to know the
>>>
>>>
>>actual
>>
>>
>>>sequence position values.
>>>
>>>thanks
>>>rich
>>>_______________________________________________
>>>BiO_Bulletin_Board maillist -
>>>
>>>
>>BiO_Bulletin_Board at bioinformatics.org
>>
>>https://bioinformatics.org/mailman/listinfo/bio_bulletin_board
>>
>>
>>>
>>>
>>>
>>>
>>_______________________________________________
>>BiO_Bulletin_Board maillist -
>>BiO_Bulletin_Board at bioinformatics.org
>>
>>
>>
>https://bioinformatics.org/mailman/listinfo/bio_bulletin_board
>
>
>>yzman at gwatson.eclipse.co.uk
>>
>>
>Received: from
>
>
>>mxa02.ex.eclipse.net.uk (mxa02.ex.eclipse.net.uk
>>[212.104.129.120])
>>
>>
> by
>
>
>>mda02.mail.eclipse.net.uk (Postfix) with ESMTP
>>id 4D5973CC179
>>
>>
> for
>
>
>><yzman at gwatson.eclipse.co.uk>; Mon, 12 Jul 2004
>>00:31:12 +0100 (BST)
>>
>>
>Received: from localhost
>
>
>>(localhost.localdomain [127.0.0.1])
>>
>>
> by
>
>
>>mxa02.ex.eclipse.net.uk (Postfix) with ESMTP id
>>BC0AE34EAB5
>>
>>
> for <yzman at gwatson.eclipse.co.uk>;
>
>
>>Mon, 12 Jul 2004 00:32:48 +0100 (BST)
>>
>>
>Received:
>
>
>>from mxa02.ex.eclipse.net.uk ([127.0.0.1])
>>
>>
> by
>
>
>>localhost (mxa02.ex.eclipse.net.uk [127.0.0.1])
>>(amavisd-new, port 10024)
>>
>>
> with LMTP id
>
>
>>15113-01-13 for <yzman at gwatson.eclipse.co.uk>;
>>
>>
>
>
>
>>Mon, 12 Jul 2004 00:32:48 +0100 (BST)
>>
>>
>Received:
>
>
>>from mx39.smf.ebay.com (mxsmfpool16.ebay.com
>>[66.135.209.213])
>>
>>
> by mxa02.ex.eclipse.net.uk
>
>
>>(Postfix) with ESMTP id 9158134E6E3
>>
>>
> for
>
>
>><yzman at gwatson.eclipse.co.uk>; Mon, 12 Jul 2004
>>00:32:47 +0100 (BST)
>>
>>
>Received: from
>
>
>>sjcbat01.sjc.ebay.com (sjcbat01.sjc.ebay.com
>>[10.6.37.40])
>>
>>
> by mx39.smf.ebay.com
>
>
>>(8.12.3/8.12.3) with ESMTP id
>>i6BNbVTE015225
>>
>>
> for
>
>
>><yzman at gwatson.eclipse.co.uk>; Sun, 11 Jul 2004
>>16:37:35 -0700
>>
>>
>Message-Id:
>
>
>><200407112337.i6BNbVTE015225 at mx39.smf.ebay.com>
>>
>>
>X-eBay-MailTracker:
>
>
>>10007.361.3.x
>>
>>
>To:
>
>
>>yzman at gwatson.eclipse.co.uk
>>
>>
>From:
>
>
>>dailystatus at ebay.co.uk
>>
>>
>Mime-Version:
>
>
>>1.0
>>
>>
>Content-Type: multipart/alternative;
>
>
>>boundary="==_MIME-Boundary-1_=="
>>
>>
>Subject: eBay
>
>
>>Watch List: 11-Jul-04 18:19:45 BST
>>
>>
>Date: Sun,
>
>
>>11 Jul 2004 16:36:33 PDT
>>
>>
>X-Virus-Scanned: by
>
>
>>Eclipse VIRUSshield at eclipse.net.uk
>>
>>
>
>
>
>
>_______________________________________________
>BiO_Bulletin_Board maillist - BiO_Bulletin_Board at bioinformatics.org
>https://bioinformatics.org/mailman/listinfo/bio_bulletin_board
>
>
>
>
--
Iddo Friedberg, Ph.D.
The Burnham Institute
10901 N. Torrey Pines Rd.
La Jolla, CA 92037 USA
Tel: +1 (858) 646 3100 x3516
Fax: +1 (858) 713 9930
http://ffas.ljcrf.edu/~iddo
More information about the BBB
mailing list