<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div><div>On Apr 10, 2013, at 8:25 AM, Ibai Gazpio wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div>Dear,<br><br>We are developing a chimera extension and I would like to know if<br>there is a way to detect or select those models which have an internal<br>sequence gap, represented by the program as a pseudobond in order to<br>distinguish them from those who haven't. Is there a function to do<br>that? Thank you very much and congratulations for the program.</div></blockquote><br></div><div>Hi Ibai,</div><div><span class="Apple-tab-span" style="white-space:pre">        </span>There are a couple of ways I can think of.  One is go through the model's sequences (m.sequences()) and check the sequences' "residues" list.  There will be Nones in the list when residues in the SEQRES records don't exist in the structure.  You will need to ignore/strip leading/trailing series of Nones and then see if any Nones remain, which would indicate an internal gap.</div><div><span class="Apple-tab-span" style="white-space:pre">    </span>The other way is to get the pseudobond group that depicts the gap.  This will do that:</div><div><br></div><div>import chimera</div><div>pbg = chimera.misc.getPseudoBondGroup(chimera.LONGBOND_PBG_NAME, create=False)</div><div><br></div><div>Then if pbg is not None, go through its pseudoBonds list and check one atom of each to see if its "molecule" attr is your model.  Of course, if you know that you only have one model open then simply check if the pseudoBonds list is empty or not.</div><div><span class="Apple-tab-span" style="white-space:pre">   </span>Let me know if you need more info.  I'm happy to hear you're developing an extension!</div><div><br></div><div>--Eric</div><br><div>
<span class="Apple-style-span" style="border-collapse: separate; border-spacing: 0px 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 16px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: auto; -khtml-text-decorations-in-effect: none; text-indent: 0px; -apple-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><p style="margin: 0.0px 0.0px 0.0px 0.0px"><font face="Helvetica" size="5" style="font: 16.0px Helvetica"><span class="Apple-converted-space">                       <span class="Apple-converted-space"> </span></span>Eric Pettersen</font></p><p style="margin: 0.0px 0.0px 0.0px 0.0px"><font face="Helvetica" size="5" style="font: 16.0px Helvetica"><span class="Apple-converted-space">                       <span class="Apple-converted-space"> </span></span>UCSF Computer Graphics Lab</font></p><p style="margin: 0.0px 0.0px 0.0px 0.0px"><font face="Helvetica" size="5" style="font: 16.0px Helvetica"><span class="Apple-converted-space">                        </span><a href="http://www.cgl.ucsf.edu">http://www.cgl.ucsf.edu</a></font></p><br class="Apple-interchange-newline"></span>
</div>
<br></body></html>