<html>
<head>
<style>
body {
font-family: Verdana, sans-serif;
font-size: 0.8em;
color:#484848;
}
h1, h2, h3 { font-family: "Trebuchet MS", Verdana, sans-serif; margin: 0px; }
h1 { font-size: 1.2em; }
h2, h3 { font-size: 1.1em; }
a, a:link, a:visited { color: #2A5685;}
a:hover, a:active { color: #c61a1a; }
a.wiki-anchor { display: none; }
fieldset.attachments {border-width: 1px 0 0 0;}
hr {
width: 100%;
height: 1px;
background: #ccc;
border: 0;
}
span.footer {
font-size: 0.8em;
font-style: italic;
}
</style>
</head>
<body>
Issue #3146 has been updated by Peter Cock.
<ul>
<li><strong>Description</strong> updated (<a title="View differences" href="https://redmine.open-bio.org/journals/diff/15344?detail_id=1617">diff</a>)</li>
</ul>
<p>It would appear the PDB parser is more noisy about warnings, but were you able to run the DSSP part as well Vincent?</p>
<hr />
<h1><a href="https://redmine.open-bio.org/issues/3146#change-15344">Bug #3146: DSSP ungraceful failure</a></h1>
<ul><li>Author: Patrick empty</li>
<li>Status: New</li>
<li>Priority: Low</li>
<li>Assignee: Biopython Dev Mailing List</li>
<li>Category: Main Distribution</li>
<li>Target version: 1.53</li>
<li>URL: </li></ul>
<p>The DSSP annotator should probably fail gracefully when the PDBParser and DSSP disagree about the existence of a residue at a certain position. Here, DSSP reports values for residue 115 of chain A, while the PDBParser throws a key error.</p>
<p>from Bio.PDB import PDBParser<br />parser = PDBParser()<br />from Bio.PDB.DSSP import DSSP<br />structure=parser.get_structure("2p0i","pdb2p0i.ent")<br />model=structure<sup><a href="#fn0">0</a></sup><br />dssp=<abbr title="model, "pdb2p0i.ent"">DSSP</abbr><br />Traceback (most recent call last):<br /> File "<stdin>", line 1, in <module><br /> File "/usr/lib/pymodules/python2.6/Bio/PDB/DSSP.py", line 175, in <i>init</i><br /> res=chain[res_id]<br /> File "/usr/lib/pymodules/python2.6/Bio/PDB/Chain.py", line 71, in <i>getitem</i><br /> return Entity.__getitem__(self, id)<br /> File "/usr/lib/pymodules/python2.6/Bio/PDB/Entity.py", line 38, in <i>getitem</i><br /> return self.child_dict[id]<br />KeyError: (' ', 115, ' ')</p>
<blockquote><blockquote><blockquote>
<p>model['A'][114]</p>
</blockquote></blockquote></blockquote>
<p><Residue GLN het= resseq=114 icode= ></p>
<blockquote><blockquote><blockquote>
<p>model['A'][116]</p>
</blockquote></blockquote></blockquote>
<p><Residue TYR het= resseq=116 icode= ></p>
<blockquote><blockquote><blockquote>
<p>model['A'][115]</p>
</blockquote></blockquote></blockquote>
<p>Traceback (most recent call last):<br /> File "<stdin>", line 1, in <module><br /> File "/usr/lib/pymodules/python2.6/Bio/PDB/Chain.py", line 71, in <i>getitem</i><br /> return Entity.__getitem__(self, id)<br /> File "/usr/lib/pymodules/python2.6/Bio/PDB/Entity.py", line 38, in <i>getitem</i><br /> return self.child_dict[id]<br />KeyError: (' ', 115, ' ')</p>
<fieldset class="attachments"><legend>Files</legend>
<a href="https://redmine.open-bio.org/attachments/download/1791/pdb2p0i.ent">pdb2p0i.ent</a>
(2.24 MB)<br />
</fieldset>
<hr />
<span class="footer"><p>You have received this notification because you have either subscribed to it, or are involved in it.<br />To change your notification preferences, please click here and login: <a class="external" href="http://redmine.open-bio.org">http://redmine.open-bio.org</a></p></span>
</body>
</html>