[Chimera-users] triggering events when findhbond finds a bond

Eric Pettersen pett at cgl.ucsf.edu
Mon Jun 4 15:55:51 PDT 2012


On Jun 4, 2012, at 2:42 PM, Anne Bowen wrote:

> Hi,
>
> 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)
>  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!

Hi Anne,
	You will have to resort to Python to do what you want to do here.   
Very simple 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.
	So first, in the "Per-Frame" dialog change the "Interpret script as"  
to "Python".  Then put this in as your script:

from chimera import runCommand
runCommand("findhbond makePseudobonds true lineWidth 2.0 color green  
intramodel true intermodel false selRestrict any")
if len(chimera.misc.getPseudoBondGroup("hydrogen bonds").pseudoBonds)  
 > 0:
	runCommand("color red")
else:
	runCommand("~color; color byhet")


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.
	I'm guessing you'd like it if MD Movie's Analysis->Plot menu item  
allowed for plotting number of H-bonds?

--Eric

                         Eric Pettersen
                         UCSF Computer Graphics Lab
                         http://www.cgl.ucsf.edu


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://plato.cgl.ucsf.edu/pipermail/chimera-users/attachments/20120604/288e45fe/attachment.html>


More information about the Chimera-users mailing list