[Biojava-l] PDBFileParser question using PDBID 470D
Steve Darnell
darnells at dnastar.com
Mon Nov 29 20:51:06 UTC 2010
Greetings,
After parsing PDBID 470D with biojava-3.0-alpha5, Chain A returns an
empty SEQRES sequence (Chain.getSeqResSequence) and empty SEQRES group
list (Chain.getSeqResGroups) but the one-letter ATOM sequence is
properly translated and the ATOM group list contains the appropriate
number of groups (LoadChemCompInfo set to true).
This is an exotic sequence, but my expectation is that the SEQRES group
list would have members in it (and one-letter sequence translated if
LoadChemCompInfo is true). Am I mistaken and the current behavior is
the intended result?
Best regards,
Steve Darnell
--
SEQRES records exist in 470D:
SEQRES 1 A 12 C43 G48 C43 G48 A44 A44 U36 U36 C43 G48 C43 G48
SEQRES 1 B 12 C43 G48 C43 G48 A44 A44 U36 U36 C43 G48 C43 G48
Sample println output (ln 1 record type, ln 2 get${TYPE}Sequence, ln 3
get${TYPE}Groups):
SEQRES
''
[]
ATOM
'CGCGAAUUCGCG'
[PDB: C43 1 trueatoms: 21, PDB: G48 2 trueatoms: 27, PDB: C43 3
trueatoms: 24, ...]
More information about the Biojava-l
mailing list