<span class="gmail_quote"></span>I strongly beleive this can be done using a 'grep' and 'substitution' commands in linux. <br>Since the ligand data starts with HETATM, they can  be easily grepable. <br><br>No need any 'parser' - a simple perl script can do this ! 
<br>Cheers, <br><span class="sg">Shameer Khadar<br>NCBS-TIFR <br>B'lore </span><div><span class="e" id="q_10e6452923636c0a_2"><br> <br><div><span class="gmail_quote">On 10/20/06, <b class="gmail_sendername">Tomek Jarzynka
</b> <<a href="mailto:tomee@genesilico.pl" target="_blank" onclick="return top.js.OpenExtLink(window,event,this)">tomee@genesilico.pl</a>> wrote:
</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Hi,<br><br>I would like to create a script that would take a PDB file as input,<br>try and identify the ligand structures and delete them from the PDB
<br>file. I figured this would be possible with Biopython by looking at<br>the atom id (whether it contains 'H_') and subclassing Select to<br>not allow those items to be written to a file.<br>Is there any more elegant way of doing this, perhaps another PDB parser
<br>framework, or maybe someone has already done similar work?<br><br>Thanks in advance,<br>--<br>Tomasz K. Jarzynka  /  +48 601 706 601  /  tomee(a-t)genesilico(d-o-t)pl<br>Laboratory of Bioinformatics and Protein Engineering | 
<a href="http://www.genesilico.pl" target="_blank" onclick="return top.js.OpenExtLink(window,event,this)">www.genesilico.pl</a><br>International Institute of Molecular and Cell Biology | <a href="http://www.iimcb.gov.pl" target="_blank" onclick="return top.js.OpenExtLink(window,event,this)">
www.iimcb.gov.pl</a><br><br>"You can have either freedom of speech or quality of communication.
<br>        -- Orson Scott Card"<br>_______________________________________________<br>BioPython mailing list  -  <a href="mailto:BioPython@lists.open-bio.org" target="_blank" onclick="return top.js.OpenExtLink(window,event,this)">
BioPython@lists.open-bio.org</a><br><a href="http://lists.open-bio.org/mailman/listinfo/biopython" target="_blank" onclick="return top.js.OpenExtLink(window,event,this)">
http://lists.open-bio.org/mailman/listinfo/biopython</a><br></blockquote></div><br>

</span></div>