<div dir="ltr">Hi Patrick,<div><br></div><div>I dug up an old PDB parser I wrote for another purpose and made a new function to find the breaks using the method I suggested, you can maybe tweak it. </div><div><br></div><div>Lenna does bring up some legit concerns. This old script was made to compare distances, so perhaps you could hack it to try Lenna's Ca-Ca distance test.</div><div><br></div><div>Let me know, and I'll send it to you off list.</div><div><br></div><div>Dave</div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Feb 10, 2015 at 6:31 PM, PC <span dir="ltr"><<a href="mailto:patrick.cossins@inbox.com" target="_blank">patrick.cossins@inbox.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><u></u>
<div>
<div>Hi David,</div><div><br></div><div>Thank you, yes I am in the process of trying to do this.</div><div><br></div><div>Thank you for the suggestion of the spacing too, something I didn't think of.</div><div><br></div><div>Patrick </div><div><div class="h5"><br><br><blockquote style="PADDING-LEFT:5px;MARGIN-LEFT:5px;BORDER-LEFT:#0000ff 2px solid;MARGIN-RIGHT:0px"><div>-----Original Message-----<br><b>From:</b> <a href="mailto:davidsshin@lbl.gov" target="_blank">davidsshin@lbl.gov</a><br><b>Sent:</b> Tue, 10 Feb 2015 17:46:35 -0800<br><b>To:</b> <br><b>Subject:</b> Re: [Biopython] Extracting a PDB list<br><br></div><div><div><div dir="ltr">Hi Patrick,<div><br></div><div>You should be able to write a script to do this (shell script with some python or awk).</div><div><br></div><div>Off the top of my head, for each file you would:</div><div><br></div><div>for each file:</div><div> extract the lines with ^ATOM into a new file to make things easier</div><div> read each line into some list</div><div> subtract the residue number from each line from the next line in the list</div><div> if that value is > 1 </div><div> print something ( the file name, or some flag)</div><div> else there are no breaks... can do something else if you want</div><div>end</div><div><br></div><div>The only tough parts are using spaces to separate items. If say a protein had 1000 residues, then the 1000 will run into the chain ID. So that's something to consider. Using specific column numbers would be the better way. </div><div><br></div><div>That and I'm not sure about the uniformity of PDB files that are really old.</div><div><br></div><div>Let me know if that helps, if not, I can maybe help out further.</div><div><br></div><div>Dave</div><div><br></div><div><br></div><div><br></div><div><br><div>On Tue, Feb 10, 2015 at 2:24 PM, João Rodrigues <span dir="ltr"><<a href="mailto:j.p.g.l.m.rodrigues@gmail.com" target="_blank">j.p.g.l.m.rodrigues@gmail.com</a>></span> wrote:<br><blockquote><div dir="ltr">Hi,<div><br></div><div>Without manually checking every single one, there is no such list, at least that I know of. Your best bet could be to reduce your resolution as low as possible, usually those structures are of very good quality.</div><div><br></div><div>Cheers,</div><div><br></div><div>João</div></div><div><div><div><br><div>2015-02-10 22:35 GMT+01:00 PC <span dir="ltr"><<a href="mailto:patrick.cossins@inbox.com" target="_blank">patrick.cossins@inbox.com</a>></span>:<br><blockquote>Hi,<br>
<br>
I do know about PISCES lists but I want a list of PDB's without any chain breaks.<br>
<br>
Is there such a list or a way to obtain such a list?<br>
<br>
Thank you.<br>
<br>
____________________________________________________________<br>
FREE 3D EARTH SCREENSAVER - Watch the Earth right on your desktop!<br>
Check it out at <a href="http://www.inbox.com/earth" target="_blank">http://www.inbox.com/earth</a><br>
<br>
<br>
<br>
_______________________________________________<br>
Biopython mailing list - <a href="mailto:Biopython@mailman.open-bio.org" target="_blank">Biopython@mailman.open-bio.org</a><br>
<a href="http://mailman.open-bio.org/mailman/listinfo/biopython" target="_blank">http://mailman.open-bio.org/mailman/listinfo/biopython</a><br>
</blockquote></div><br></div>
</div></div><br>_______________________________________________<br>
Biopython mailing list - <a href="mailto:Biopython@mailman.open-bio.org" target="_blank">Biopython@mailman.open-bio.org</a><br>
<a href="http://mailman.open-bio.org/mailman/listinfo/biopython" target="_blank">http://mailman.open-bio.org/mailman/listinfo/biopython</a><br></blockquote></div><br><br clear="all"><div><br></div>-- <br><div><div>David Shin, Ph.D</div><div>Lawrence Berkeley National Labs</div><div>1 Cyclotron Road</div><div>MS 83-R0101</div><div>Berkeley, CA 94720</div><div>USA</div></div>
</div></div>
</div></div></blockquote>
</div></div><hr size="1px" noshade style="clear:both;margin-top:10px;min-height:1px">
<div style="font:12px Verdana,sans-serif;color:Black;background:white;padding:3px;line-height:1.3em">
<a href="http://www.inbox.com/earth" target="_blank"><img src="http://my.inbox.com/img/ftrs/earth.jpg" width="100" height="90" alt="3D Earth Screensaver Preview" align="left" border="0" style="margin-right:15px"></a>
<br><strong><font color="2086c3">Free 3D Earth Screensaver</font></strong><br>
<u></u>Watch the Earth right on your desktop!<u></u> <u></u>Check it out at <a href="http://www.inbox.com/earth" target="_blank">www.inbox.com/earth</a><u></u></div>
</div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature"><div>David Shin, Ph.D</div><div>Lawrence Berkeley National Labs</div><div>1 Cyclotron Road</div><div>MS 83-R0101</div><div>Berkeley, CA 94720</div><div>USA</div></div>
</div>