<html><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div><div>On Jun 4, 2012, at 2:42 PM, Anne Bowen wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"> <meta http-equiv="Content-Type" content="text/html; charset=us-ascii"> <div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "> <div apple-content-edited="true"><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: -webkit-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; font-size: medium; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: -webkit-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; font-size: medium; "> <div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "> Hi,</div> <div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "> <br> </div> <div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "> I currently have a "Per-Frame Command" set up (in the MD Movie dialog) that finds the H bonds (findhbond makePseudobonds true lineWidth 2.0 color green intramodel true intermodel false selRestrict any) </div> <div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "> every frame of a trajectory. Whenever a hbond is detected, I would like to turn the entire structure a different color just for that frame, and then have it go back to being colored by heteroatom. Is there an easy way to accomplish this? I tried to look through the documentation to determine if "findhbond" could return true or false if an h bond was found and it did not seem like it was possible. Is the only way to do this to change the findhbond script itself? Any help or suggestions is greatly appreciated!</div></span></span></div></div></blockquote><br></div><div>Hi Anne,</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>You will have to resort to Python to do what you want to do here. <i>Very simple</i> Python, but Python nonetheless. The reason it can be simple is that there is a Python command named runCommand that allows you to execute regular Chimera commands and therefore the bulk of the script is just Python calls of commands you already know, with a little extra logic/glue to perform the test for the H-bonds and execute the corresponding color command.</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>So first, in the "Per-Frame" dialog change the "Interpret script as" to "Python". Then put this in as your script:</div><div><br></div><div>from chimera import runCommand</div><div>runCommand("findhbond makePseudobonds true lineWidth 2.0 color green intramodel true intermodel false selRestrict any")</div><div>if len(chimera.misc.getPseudoBondGroup("hydrogen bonds").pseudoBonds) > 0:</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>runCommand("color red")</div><div>else:</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>runCommand("~color; color byhet")</div><div><br></div><div><br></div><div>I chose red as the "other color"; change it to whatever you like. The "~color" is needed before the "color byhet" so that the carbon atoms will go to their non-red color.</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>I'm guessing you'd like it if MD Movie's Analysis->Plot menu item allowed for plotting number of H-bonds?</div><div><br></div><div>--Eric</div><div><br></div><div apple-content-edited="true"> <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; "><div style="word-wrap: break-word; -khtml-nbsp-mode: space; -khtml-line-break: after-white-space; "><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"></div></span> </div><br></body></html>