[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