[BioPython] HETATM records retrieval
O.Doehring at cs.ucl.ac.uk
O.Doehring at cs.ucl.ac.uk
Sat Jul 14 18:38:22 UTC 2007
Dear community,
how should HETATM records be retrieved via Biopython? I assume it should be
somewhere on the chain or residue level:
-
http://www.biopython.org/DIST/docs/api/public/Bio.PDB.Residue.Residue-class.html
- http://biopython.org/DIST/docs/api/private/Bio.PDB.Chain.Chain-class.html
Using the following basic sample code :
for model in self.structure.get_iterator():
for chain in model.get_iterator():
print chain.__repr__()
for residue in chain.get_iterator():
print residue.__repr__()
applied to protein 1DHR (http://www.pdb.org/pdb/files/1dhr.pdb) we get:
<Chain id= >
<Residue GLU het= resseq=5 icode= >
<Residue ALA het= resseq=6 icode= >
<Residue ARG het= resseq=7 icode= >
<Residue ARG het= resseq=8 icode= >
<Residue VAL het= resseq=9 icode= >
<Residue LEU het= resseq=10 icode= >
...
<Residue ALA het= resseq=238 icode= >
<Residue TYR het= resseq=239 icode= >
<Residue PHE het= resseq=240 icode= >
As one can see that are all ATOM records. Thanks.
Yours,
Orlando
More information about the Biopython
mailing list