[BioPython] Extracting residue list from PDB
Gad Abraham
gabraham at cs.rmit.edu.au
Fri Oct 28 22:14:27 EDT 2005
On Fri, Oct 28, 2005 at 09:07:53AM -0700, Iddo Friedberg wrote:
> Hi Gad,
>
> The reason the chains seem to be of the wrong length, is that they are
> generated from the ATOM records, rather than the SEQRES records. Those
> disagree often enough in PDB files.
> This problem does not exists in mmCIF, I believe.
>
> Using your code, I got only six chains, so I cannot comment on the
> second problem.
>
I see. I'm consistently getting 10 chains (lengths 161, 804, 97, 176,
11, 71, 87, 662, 133, 286) for 1N62.
It seems that parsing the SEQRES is the way go.
Thanks,
Gad
--
+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
gabraham at cs.rmit.edu.au
http://yallara.cs.rmit.edu.au/~gabraham
+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
More information about the BioPython
mailing list