Animate a Reacton Mechanism?

I constructed a movie for the first part of a reaction mechanism (substrate docking). I would like the animated mechanism to continue as follows. The substrate donates a proton to a tyrosine residue of the enzyme, forming an intermediate. The intermediate rotates 180 (I have accomplished this part) and the proton rebounds to the intermediate forming a new product, which is released. I can animate the product release -- it's the reverse of substrate docking. The part of the animation I cannot envisage is how to smoothly show or morph the bond transitions (bond breaking and making) during the reaction progression. IAny advice on whether Chimera or other ancillary software can do this is appreciated.

Hi Kevin, I imagine you could save positions (savepos command) for the bond-broken and bond-formed structures of the complex and then use the reset command to go from one to the other over a specified number of frames. You could also save a position at one or more intermediate states, if you wanted to do something at that state (make the "bond" appear or disappear). The molecules that move relative to one another should be separate models, and since you can't make a bond between separate models, you would need to make a pseudobond. It could be shown as wire, or dashed line, or stick. Pseudobonds automatically "stretch" when the endpoints move relative to one another. You can create it ahead of time (not being a programmer, I would just make a distance monitor and then in Distances tool set its label to "None" and use "display options" button to set initial display style/color), save it in the session, and then with the setattr command control display style, color, and when it appears or disappears. With [a lot] more work you could gradually fade it in or out by changing the transparency of its color, but there is no single command to do the whole fade, it would take lots of individual color-changing commands. Example commands to change pseudobond color and and hide pseudobond: setattr p color 1,0,0,.5 setattr p display 0 ... where you could put the atomspec of the two endpoint atoms at the end of each command to control only that pseudobond, instead of all pseudobonds. With python there may be additional possibilities, such as adding a new atom to one model that is in the same exact place as the real bond partner atom in the other model, then adding an intramodel bond to that new atom, but I think that would be harder anyway, and more limiting, since the bond wouldn't automatically stretch like the pseudobond. To show it changing length would require changing the internal coordinates of the model. I hope this helps, Elaine ---------- Elaine C. Meng, Ph.D. UCSF Computer Graphics Lab (Chimera team) and Babbitt Lab Department of Pharmaceutical Chemistry University of California, San Francisco On Feb 19, 2012, at 4:09 PM, Kevin D. Walker wrote:
I constructed a movie for the first part of a reaction mechanism (substrate docking). I would like the animated mechanism to continue as follows. The substrate donates a proton to a tyrosine residue of the enzyme, forming an intermediate. The intermediate rotates 180 (I have accomplished this part) and the proton rebounds to the intermediate forming a new product, which is released.
I can animate the product release -- it's the reverse of substrate docking. The part of the animation I cannot envisage is how to smoothly show or morph the bond transitions (bond breaking and making) during the reaction progression. IAny advice on whether Chimera or other ancillary software can do this is appreciated.

On Feb 21, 2012, at 11:25 AM, Elaine Meng wrote:
setattr p color 1,0,0,.5 setattr p display 0
... where you could put the atomspec of the two endpoint atoms at the end of each command to control only that pseudobond, instead of all pseudobonds.
If you don't like the way the change in transparency looks, another alternative is to have the breaking bond become gradually thinner (attribute "radius") while the forming bond becomes thicker. It's just as awkward to script because you would need a setattr command each frame the motion is happening (luckily it's easy to have motions go on as other commands are executed, just don't use a "wait" command after the motion command). --Eric Eric Pettersen UCSF Computer Graphics Lab http://www.cgl.ucsf.edu

Hi Kevin, Another trick to animate bonds appearing and disappearing is to use the movie crossfade command. This allows you to fade between any two scenes. So you start with a bond showing, then use movie crossfade frames 25 and the following command undisplays the bond. This will cause the movie to fade between the bond shown and bond hidden state over 25 frames. In general it fades between the scene as it appeared just before the command to the scene just after the next command is executed. The fade does not appear on-screen -- only in the recorded movie. Tom
On Feb 21, 2012, at 11:25 AM, Elaine Meng wrote:
setattr p color 1,0,0,.5 setattr p display 0
... where you could put the atomspec of the two endpoint atoms at the end of each command to control only that pseudobond, instead of all pseudobonds.
If you don't like the way the change in transparency looks, another alternative is to have the breaking bond become gradually thinner (attribute "radius") while the forming bond becomes thicker.
It's just as awkward to script because you would need a setattr command each frame the motion is happening (luckily it's easy to have motions go on as other commands are executed, just don't use a "wait" command after the motion command).
--Eric
Eric Pettersen
UCSF Computer Graphics Lab
_______________________________________________ Chimera-users mailing list Chimera-users@cgl.ucsf.edu http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users

the 'crossfade' command omits bonds that are selected for later rotation in a command line script. See attached.

I want to change the geometry of a selected atom by executing a command line script. I would like to proceed as follows: 1) select [atom-id] (this is the easy part) ...then I want to apply a script for the following interpretation 2) modify structure; set Element to 'C'; bonds to '4'; geometry to 'tetrahedral' Is there a command script for this series? Kevin

On Feb 22, 2012, at 2:06 PM, Kevin D. Walker wrote:
I want to change the geometry of a selected atom by executing a command line script. I would like to proceed as follows:
1) select [atom-id] (this is the easy part)
...then I want to apply a script for the following interpretation
2) modify structure; set Element to 'C'; bonds to '4'; geometry to 'tetrahedral'
Is there a command script for this series?
Well, there is no command version of this part of Build Structure. Nonetheless everything is doable by calling the underlying Python. I've attached a Python script that does exactly what you requested above. It's got some explanatory comments. You would definitely want to look through it before running it because it offers some options for how you want the modified atom named. You would run the Python script from a regular command script with "open ~/c4.py" (assuming you named it c4.py and put it in your home directory). --Eric Eric Pettersen UCSF Computer Graphics Lab http://www.cgl.ucsf.edu
participants (4)
-
Elaine Meng
-
Eric Pettersen
-
Kevin D. Walker
-
Tom Goddard