<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=Windows-1252">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Arial,Helvetica,sans-serif;">
<p>Hi again Chimera Team,</p>
<p><br>
</p>
<p>I have 2 questions, one specific and one quite general. If you see the emails below, I was speaking with Jaime about computing basic structure characteristics (pKa, net surface charge, proportions/predominance of 2ary structure etc).</p>
<p><br>
</p>
<p>The general question I have is what suggestions might you have for functions within chimera capable of computing any or all of these features? (ideally commandline based so that I can run this over many PDB structures).</p>
<p><br>
</p>
<p>The specific question that pertains to the below is that Jaime mentioned each residue has an int value that corresponds to a particular structure designation, ints from -1 to 6. Could you let me know what the correspondence is between these values and actual
structure?</p>
<p><br>
</p>
<div id="Signature">
<div id="divtagdefaultwrapper" style="font-size:12pt; color:#000000; background-color:#FFFFFF; font-family:Calibri,Arial,Helvetica,sans-serif">
Joe Healey<br style="font-size:10pt">
</div>
<div id="divtagdefaultwrapper" style="font-size:12pt; color:#000000; background-color:#FFFFFF; font-family:Calibri,Arial,Helvetica,sans-serif">
<br>
</div>
<div id="divtagdefaultwrapper" style="font-size:12pt; color:#000000; background-color:#FFFFFF; font-family:Calibri,Arial,Helvetica,sans-serif">
<span style="text-decoration:underline"> </span>
<br>
</div>
<div id="divtagdefaultwrapper" style="font-size:12pt; color:#000000; background-color:#FFFFFF; font-family:Calibri,Arial,Helvetica,sans-serif">
<span style="font-size:10pt">M.Sc. B.Sc. (Hons)</span><br>
</div>
<div id="divtagdefaultwrapper" style="font-size:12pt; color:#000000; background-color:#FFFFFF; font-family:Calibri,Arial,Helvetica,sans-serif">
<span style="font-size:10pt"></span>
<div><span style="font-size:10pt">PhD Student</span></div>
<div><span style="font-size:10pt">MOAC CDT, Senate House</span></div>
<div><span style="font-size:10pt">University of Warwick<br style="">
</span></div>
<div><span style="font-size:10pt">Coventry<br style="">
</span></div>
<div><span style="font-size:10pt">CV47AL<br style="">
Mob: +44 (0) 7536 042620 | </span><span style="font-size:10pt">Email: J.R.J.Healey@warwick.ac.uk</span></div>
<div><span style="font-size:10pt"><br style="">
</span></div>
<div><span style="font-size:10pt">Jointly working in:<br style="">
</span></div>
<div><span style="font-size:10pt"><a href="http://www2.warwick.ac.uk/fac/med/research/tsm/microinfect/staff/waterfieldlab/" id="LPNoLP">Waterfield Lab</a> (<span style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:13.3333330154419px; background-color:rgb(255,255,255)">WMS </span><span style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:13.3333330154419px; background-color:rgb(255,255,255)">Microbiology
and Infecti</span><span style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:13.3333330154419px; background-color:rgb(255,255,255)">on Un</span><span style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:13.3333330154419px; background-color:rgb(255,255,255)">it)</span></span></div>
<div><span style="font-size:10pt">and the <a href="http://www2.warwick.ac.uk/fac/sci/chemistry/research/gibson/gibsongroup/" id="LPNoLP">Gibson Lab</a> (Warwick Chemistry)<br style="">
</span></div>
<div><span style="font-size:10pt"><br style="">
</span></div>
<div><span style="font-size:10pt">Twitter: <a href="https://twitter.com/JRJHealey" id="LPNoLP">
@JRJHealey</a> | </span><span style="font-size:10pt">Website: </span><a href="http://www2.warwick.ac.uk/fac/sci/moac/people/students/2013/joseph_healey" style="font-size:10pt" id="LPNoLP">MOAC Page</a></div>
</div>
</div>
<br>
<br>
<div style="color: rgb(0, 0, 0);">
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>From:</b> Jaime Rodríguez-Guerra <jaime.rodriguezguerra@uab.cat><br>
<b>Sent:</b> 08 October 2016 09:55<br>
<b>To:</b> Healey, Joe<br>
<b>Subject:</b> Re: [Chimera-users] Recursive structure matching and acquisition of descriptive numbers</font>
<div> </div>
</div>
<div>
<div dir="ltr">Hi again Joe!
<div><br>
</div>
<div>No problem at all; I'll be pleased to help you.</div>
<div><br>
</div>
<div>For charges, I think Chimera can calculate them as part of the Surface/Binding analysis> Dock Prep wizard. You also have some information about 2nd structure inside chimera.Residue objects, specifically ssId attribute. However, these are int values, and
I don't really know which ss is each (ranging from -1 to 6, as I have observed).</div>
<div><br>
</div>
<div>For pKa, you can use PropKa (<a href="https://github.com/jensengroup/propka-3.1">https://github.com/jensengroup/propka-3.1</a>), and then I guess you'll be able to get the pI somewhere around the code. I've been designing a GUI to use this in Chimera,
but it's not yet ready... Sorry about that!</div>
<div><br>
</div>
<div>I don't know if you need anything else... but I think the Python community will have you covered one way or another! BioPython.PDB might be useful too (<a href="http://biopython.org/wiki/The_Biopython_Structural_Bioinformatics_FAQ">http://biopython.org/wiki/The_Biopython_Structural_Bioinformatics_FAQ</a>).</div>
<div><br>
</div>
<div>Cheers,</div>
<div>Jaime.</div>
<div><br>
</div>
</div>
<div class="gmail_extra"><br>
<div class="gmail_quote">2016-10-06 15:57 GMT+02:00 Healey, Joe <span dir="ltr"><<a href="mailto:J.R.J.Healey@warwick.ac.uk" target="_blank">J.R.J.Healey@warwick.ac.uk</a>></span>:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex; border-left:1px #ccc solid; padding-left:1ex">
<div>
<div id="m_1745486733791591556divtagdefaultwrapper" style="font-size:12pt; color:#000000; font-family:Calibri,Arial,Helvetica,sans-serif">
<p>Hi again Jaime,</p>
<p><br>
</p>
<p>Hope you don't mind but I had a couple of questions for you. You may have figured out from all the previous correspondance that I have a load of protein structure simulations. The RMSD script you helped me with was to sort of 'broadly' gauge whether the
simulations looked realistic or not.</p>
<p><br>
</p>
<p>For some of them, I'd like to build up a chart of their tertiary structure properties based on the simulations. Things like charge/hydrophobicity, isoelectric point, 2ndary structure dominance and so on.</p>
<p><br>
</p>
<p>Could I ask you for any advice you might about what tools are available to analyse PDBs? I'm sure chimera is probably capable of it too. I've mostly found webservers at the moment for specific tasks, but I'd like to do most of the analysis locally. My group
is full of bioinformaticians, but I'm the only one who really looks at proteins so I need the advice of a computational chemist I think!</p>
<p><br>
</p>
<p>Cheers,</p>
<p><br>
</p>
<p>Joe</p>
<span class="">
<p><br>
</p>
<p><br>
</p>
<div id="m_1745486733791591556Signature">
<div id="m_1745486733791591556divtagdefaultwrapper" style="font-size:12pt; color:#000000; background-color:#ffffff; font-family:Calibri,Arial,Helvetica,sans-serif">
Joe Healey<br style="font-size:10pt">
</div>
<div id="m_1745486733791591556divtagdefaultwrapper" style="font-size:12pt; color:#000000; background-color:#ffffff; font-family:Calibri,Arial,Helvetica,sans-serif">
<br>
</div>
<div id="m_1745486733791591556divtagdefaultwrapper" style="font-size:12pt; color:#000000; background-color:#ffffff; font-family:Calibri,Arial,Helvetica,sans-serif">
<span style="text-decoration:underline"> </span><br>
</div>
<div id="m_1745486733791591556divtagdefaultwrapper" style="font-size:12pt; color:#000000; background-color:#ffffff; font-family:Calibri,Arial,Helvetica,sans-serif">
<span style="font-size:10pt">M.Sc. B.Sc. (Hons)</span><br>
</div>
<div id="m_1745486733791591556divtagdefaultwrapper" style="font-size:12pt; color:#000000; background-color:#ffffff; font-family:Calibri,Arial,Helvetica,sans-serif">
<span style="font-size:10pt"></span>
<div><span style="font-size:10pt">PhD Student</span></div>
<div><span style="font-size:10pt">MOAC CDT, Senate House</span></div>
<div><span style="font-size:10pt">University of Warwick<br>
</span></div>
<div><span style="font-size:10pt">Coventry<br>
</span></div>
<div><span style="font-size:10pt">CV47AL<br>
Mob: <a href="tel:%2B44%20%280%29%207536%20042620" value="+447536042620" target="_blank">
+44 (0) 7536 042620</a> | </span><span style="font-size:10pt">Email: <a href="mailto:J.R.J.Healey@warwick.ac.uk" target="_blank">
J.R.J.Healey@warwick.ac.uk</a></span></div>
<div><span style="font-size:10pt"><br>
</span></div>
<div><span style="font-size:10pt">Jointly working in:<br>
</span></div>
<div><span style="font-size:10pt"><a href="http://www2.warwick.ac.uk/fac/med/research/tsm/microinfect/staff/waterfieldlab/" id="m_1745486733791591556LPNoLP" target="_blank">Waterfield Lab</a> (<span style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:13.3333330154419px; background-color:rgb(255,255,255)">WMS </span><span style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:13.3333330154419px; background-color:rgb(255,255,255)">Microbiology
and Infecti</span><span style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:13.3333330154419px; background-color:rgb(255,255,255)">on Un</span><span style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:13.3333330154419px; background-color:rgb(255,255,255)">it)</span></span></div>
<div><span style="font-size:10pt">and the <a href="http://www2.warwick.ac.uk/fac/sci/chemistry/research/gibson/gibsongroup/" id="m_1745486733791591556LPNoLP" target="_blank">Gibson Lab</a> (Warwick Chemistry)<br>
</span></div>
<div><span style="font-size:10pt"><br>
</span></div>
<div><span style="font-size:10pt">Twitter: <a href="https://twitter.com/JRJHealey" id="m_1745486733791591556LPNoLP" target="_blank">
@JRJHealey</a> | </span><span style="font-size:10pt">Website: </span><a href="http://www2.warwick.ac.uk/fac/sci/moac/people/students/2013/joseph_healey" id="m_1745486733791591556LPNoLP" target="_blank" style="font-size:10pt">MOAC Page</a></div>
</div>
</div>
</span></div>
<hr style="display:inline-block; width:98%">
<div id="m_1745486733791591556divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>From:</b> Healey, Joe<br>
<b>Sent:</b> 23 September 2016 11:51:40<br>
<b>To:</b> Jaime Rodríguez-Guerra<br>
<b>Cc:</b> Eric Pettersen; <a href="mailto:chimera-users@cgl.ucsf.edu" target="_blank">
chimera-users@cgl.ucsf.edu</a>
<div>
<div class="h5"><br>
<b>Subject:</b> Re: [Chimera-users] Recursive structure matching and acquisition of descriptive numbers</div>
</div>
</font>
<div> </div>
</div>
<div>
<div class="h5">
<div>
<div dir="ltr">
<div id="m_1745486733791591556x_divtagdefaultwrapper" style="font-size:12pt; color:#000000; background-color:#ffffff; font-family:Calibri,Arial,Helvetica,sans-serif">
<p>Not at all! Much appreciated, as I think I mentioned, this is my first go at writing anything remotely complex in python (and not a bad turn out for about a week I don't think...)</p>
<p><br>
</p>
<p>A lot of the mess comes from copying functional bits from StackOverflow in the interests of just getting it working first, and making it pretty second as I'm sure you can appreciate (you'd hate to see my LaTeX preamble)!</p>
<p><br>
</p>
<p>Not needed to close the file with 'with' is useful to know (takes some of the headache out of knowing where to place the close).</p>
<p><br>
</p>
<p>I'll certainly make the tidy ups you've suggested, for my own learning.</p>
<p><br>
</p>
<p><br>
</p>
<p>Thanks again,</p>
<p><br>
</p>
<p>Joe :)</p>
<p><br>
</p>
<div id="m_1745486733791591556x_Signature">
<div id="m_1745486733791591556x_divtagdefaultwrapper" style="font-size:12pt; color:#000000; background-color:#ffffff; font-family:Calibri,Arial,Helvetica,sans-serif">
Joe Healey<br style="font-size:10pt">
</div>
<div id="m_1745486733791591556x_divtagdefaultwrapper" style="font-size:12pt; color:#000000; background-color:#ffffff; font-family:Calibri,Arial,Helvetica,sans-serif">
<br>
</div>
<div id="m_1745486733791591556x_divtagdefaultwrapper" style="font-size:12pt; color:#000000; background-color:#ffffff; font-family:Calibri,Arial,Helvetica,sans-serif">
<span style="text-decoration:underline"> </span><br>
</div>
<div id="m_1745486733791591556x_divtagdefaultwrapper" style="font-size:12pt; color:#000000; background-color:#ffffff; font-family:Calibri,Arial,Helvetica,sans-serif">
<span style="font-size:10pt">M.Sc. B.Sc. (Hons)</span><br>
</div>
<div id="m_1745486733791591556x_divtagdefaultwrapper" style="font-size:12pt; color:#000000; background-color:#ffffff; font-family:Calibri,Arial,Helvetica,sans-serif">
<span style="font-size:10pt"></span>
<div><span style="font-size:10pt">PhD Student</span></div>
<div><span style="font-size:10pt">MOAC CDT, Senate House</span></div>
<div><span style="font-size:10pt">University of Warwick<br>
</span></div>
<div><span style="font-size:10pt">Coventry<br>
</span></div>
<div><span style="font-size:10pt">CV47AL<br>
Mob: <a href="tel:%2B44%20%280%29%207536%20042620" value="+447536042620" target="_blank">
+44 (0) 7536 042620</a> | </span><span style="font-size:10pt">Email: <a href="mailto:J.R.J.Healey@warwick.ac.uk" target="_blank">
J.R.J.Healey@warwick.ac.uk</a></span></div>
<div><span style="font-size:10pt"><br>
</span></div>
<div><span style="font-size:10pt">Jointly working in:<br>
</span></div>
<div><span style="font-size:10pt"><a href="http://www2.warwick.ac.uk/fac/med/research/tsm/microinfect/staff/waterfieldlab/" id="m_1745486733791591556LPNoLP" target="_blank">Waterfield Lab</a> (<span style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:13.3333330154419px; background-color:rgb(255,255,255)">WMS </span><span style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:13.3333330154419px; background-color:rgb(255,255,255)">Microbiology
and Infecti</span><span style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:13.3333330154419px; background-color:rgb(255,255,255)">on Un</span><span style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:13.3333330154419px; background-color:rgb(255,255,255)">it)</span></span></div>
<div><span style="font-size:10pt">and the <a href="http://www2.warwick.ac.uk/fac/sci/chemistry/research/gibson/gibsongroup/" id="m_1745486733791591556LPNoLP" target="_blank">Gibson Lab</a> (Warwick Chemistry)<br>
</span></div>
<div><span style="font-size:10pt"><br>
</span></div>
<div><span style="font-size:10pt">Twitter: <a href="https://twitter.com/JRJHealey" id="m_1745486733791591556LPNoLP" target="_blank">
@JRJHealey</a> | </span><span style="font-size:10pt">Website: </span><a href="http://www2.warwick.ac.uk/fac/sci/moac/people/students/2013/joseph_healey" id="m_1745486733791591556LPNoLP" target="_blank" style="font-size:10pt">MOAC Page</a></div>
</div>
</div>
</div>
<hr style="display:inline-block; width:98%">
<div id="m_1745486733791591556x_divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>From:</b> Jaime Rodríguez-Guerra <<a href="mailto:jaime.rodriguezguerra@uab.cat" target="_blank">jaime.rodriguezguerra@uab.cat</a><wbr>><br>
<b>Sent:</b> 23 September 2016 10:56:38<br>
<b>To:</b> Healey, Joe<br>
<b>Cc:</b> Eric Pettersen; Jaime Rodríguez-Guerra; <a href="mailto:chimera-users@cgl.ucsf.edu" target="_blank">
chimera-users@cgl.ucsf.edu</a><br>
<b>Subject:</b> Re: [Chimera-users] Recursive structure matching and acquisition of descriptive numbers</font>
<div> </div>
</div>
</div>
<font size="2"><span style="font-size:10pt">
<div class="m_1745486733791591556PlainText">Hi Joe!<br>
<br>
You're welcome! Glad to see you got it working. I've reviewed your<br>
code and would want to point out some details that might be helpful to<br>
you in the future, though. These are mostly style-related things, but<br>
those will contribute to a cleaner reading experience!<br>
<br>
1. You don't need backslashes if you're inside parenthesis. So, in all<br>
those parse_argument lines, you can safely delete them. Related to<br>
this, you can wrap imports with parenthesis, so that backslashes are<br>
not needed.<br>
2. Tuple unpacking does not need parenthesis in you're dealing with<br>
1D-tuples. Only needed if you are dealing with more than one dimension<br>
(ie, (animal, fruit), sport = [['tiger', 'pear'], 'soccer']<br>
3. Since you are opening files with the "with" context manager, you<br>
don't need to manually close the file later. It's closed automatically<br>
as soon as you leave that with block.<br>
4. Lines 166-170 could be replaced with a (cleaner, but not<br>
necessarily more performant) glob.glob() call. Check it out to see if<br>
it satisfies your requirements.<br>
5. While you are at it, take a look at the PEP8 docs and the Google<br>
Python style guide. You don't need to follow all the rules, but they<br>
are really helpful in bringing consistency to your own style!<br>
<br>
That's it! Hope you don't mind these pieces of advice :)<br>
<br>
Cheers,<br>
Jaime.<br>
<br>
2016-09-23 11:27 GMT+02:00 Healey, Joe <<a href="mailto:J.R.J.Healey@warwick.ac.uk" target="_blank">J.R.J.Healey@warwick.ac.uk</a>>:<br>
> That was exactly what I was after thanks! I figured the object orientation<br>
> should provide that somewhere but was barking up slightly the wrong tree!<br>
><br>
><br>
> Thank you all for all your help - definitely wouldn't have been able to do<br>
> it without you!<br>
><br>
><br>
> If it's of any interest, or it helps for any future questions where you<br>
> might want to refer back to code segments, I've put the<br>
> more-or-less-finished script in this paste:<br>
><br>
><br>
> <a href="http://pastebin.com/xJNSWJGq" target="_blank">http://pastebin.com/xJNSWJGq</a><br>
><br>
><br>
> Joe Healey<br>
><br>
><br>
> M.Sc. B.Sc. (Hons)<br>
> PhD Student<br>
> MOAC CDT, Senate House<br>
> University of Warwick<br>
> Coventry<br>
> CV47AL<br>
> Mob: <a href="tel:%2B44%20%280%29%207536%20042620" value="+447536042620" target="_blank">
+44 (0) 7536 042620</a> | Email: <a href="mailto:J.R.J.Healey@warwick.ac.uk" target="_blank">
J.R.J.Healey@warwick.ac.uk</a><br>
><br>
> Jointly working in:<br>
> Waterfield Lab (WMS Microbiology and Infection Unit)<br>
> and the Gibson Lab (Warwick Chemistry)<br>
><br>
> Twitter: @JRJHealey | Website: MOAC Page<br>
> ______________________________<wbr>__<br>
> From: Eric Pettersen <<a href="mailto:pett@cgl.ucsf.edu" target="_blank">pett@cgl.ucsf.edu</a>><br>
> Sent: 22 September 2016 20:54:35<br>
> To: Jaime Rodríguez-Guerra<br>
> Cc: <a href="mailto:chimera-users@cgl.ucsf.edu" target="_blank">chimera-users@cgl.ucsf.edu</a>; Healey, Joe<br>
> Subject: Re: [Chimera-users] Recursive structure matching and acquisition of<br>
> descriptive numbers<br>
><br>
> Thanks Jaime — exactly right. I realized that the “atoms1, atoms2” in my<br>
> original example didn’t really give any indication of which atoms were<br>
> which, which is why in my later examples I changed to "simAtoms, refAtoms”,<br>
> so it’s actually “atoms2” that are the reference atoms.<br>
><br>
> —Eric<br>
><br>
>> On Sep 22, 2016, at 11:45 AM, Jaime Rodríguez-Guerra<br>
>> <<a href="mailto:jaime.rodriguezguerra@uab.cat" target="_blank">jaime.rodriguezguerra@uab.cat</a><wbr>> wrote:<br>
>><br>
>> Hi!<br>
>><br>
>> Since you are getting lists of atoms back, you only need to do some<br>
>> attribute access (rather than method calling)!<br>
>><br>
>> Every chimera.Atom object has an attribute called molecule: a<br>
>> reference to its parent chimera.Molecule. chimera.Molecule objects<br>
>> have an attribute called 'name', but sometimes this is not very<br>
>> informative (depends on your input), so you need to resort to the pdb<br>
>> filename, stored in the first element of the openedAs tuple.<br>
>><br>
>> In code, this looks like this:<br>
>><br>
>> for atoms1, atoms2, rmsd, fullRmsd in match(CP_BEST, [ref, sims],<br>
>> defaults[MATRIX], "nw", defaults[GAP_OPEN], defaults[GAP_EXTEND]):<br>
>> molecule1 = atoms1[0].molecule # any atom from the list will do<br>
>> molecule2 = atoms2[0].molecule<br>
>> print <a href="http://molecule1.name" target="_blank">molecule1.name</a>, "\t",
<a href="http://molecule2.name" target="_blank">molecule2.name</a>, "\t", rmsd #<br>
>> molecule1.openedAs[0] will also work here<br>
>><br>
>> I don't know if atoms1 comes consistently from the reference molecule,<br>
>> but I'd guess it does.<br>
>><br>
>> Hope it helps!<br>
>> ______________________________<wbr>_________________<br>
>> Chimera-users mailing list: <a href="mailto:Chimera-users@cgl.ucsf.edu" target="_blank">
Chimera-users@cgl.ucsf.edu</a><br>
>> Manage subscription:<br>
>> <a href="http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users" target="_blank">
http://plato.cgl.ucsf.edu/<wbr>mailman/listinfo/chimera-users</a><br>
><br>
</div>
</span></font></div>
</div>
</div>
</div>
</blockquote>
</div>
<br>
</div>
</div>
</div>
</div>
</body>
</html>