<html><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div><div>On Aug 3, 2009, at 5:01 PM, Elaine Meng wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><span class="Apple-style-span" style="font-family: monospace; ">If you did want the exact charges from PDB2PQR, you could reformat <br>those charges (next-to-last column in the original PQR file) into an <br>attribute assignment file and reassign them with Define Attribute (or <br>command defattr), after reading in the edited PDB files (PQR files <br>with last 2 columns removed).</span></span></blockquote><br></div><div>In this vein, I wrote a Python script (pqr2chimera.py) that will read a PQR file and output a PDB file that Chimera can read (chimera.pdb) and a defattr file (chimera.defattr) that could be used to assign the PQR charges to the model, if for some reason you prefer them to the charges that Chimera would assign. The script is attached. You run it from a Unix shell like this:</div><div><br></div><div><span class="Apple-tab-span" style="white-space:pre"> </span>python pqr2chimera.py my_pqr_file.pqr</div><div><br></div><div>The two output files will appear in the directory you run it in. One caveat is that since the defattr output file uses "@/serialNumber=" as it's atom selector, defining the charge attribute in Chimera using that file is kind of slow -- it looks through the list of atoms every time to locate the one with the matching serial number. It probably would be faster if the selector were written in the ":res@atom" form, but since I'm not even sure you care about the defattr part I took the lazy road.</div><div><br></div><div>--Eric</div><br></body><div><span class="Apple-style-span" style="border-collapse: separate; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 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; -webkit-text-decorations-in-effect: none; text-indent: 0px; -webkit-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font face="Helvetica" size="5" style="font: normal normal normal 16px/normal Helvetica; "> Eric Pettersen</font></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font face="Helvetica" size="5" style="font: normal normal normal 16px/normal Helvetica; "> UCSF Computer Graphics Lab</font></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font face="Helvetica" size="5" style="font: normal normal normal 16px/normal Helvetica; "> http://www.cgl.ucsf.edu</font></div><br class="Apple-interchange-newline"></div></span></div><div><br></div><div></div></html>