[Chimera-users] Coulombic surface coloring for a trajectory
goddard at sonic.net
Wed Nov 22 10:03:04 PST 2017
To work around this obscure Chimera bug you could compute the surface, coulombic color it, then copy it, and delete the original surface. The copy will not update automatically when coordinates change so the coloring won’t be erased. So I think the script would be
sop split #0
The first line closes the surface copy from the preceding frame, then we compute the new surface, color it, copy it (using the sop split command which will create model #1), then close the computed surface.
As Elaine told you in the original email, the molecular surface calculation fails often enough that you probably won’t get all the surfaces for a trajectory. This would work in a simpler and totally reliable way in our new ChimeraX, but we do not have the Coulombic electrostatics in ChimeraX yet so that is not an option.
> On Nov 22, 2017, at 6:54 AM, RYDER WHITTAKER HAWKINS wrote:
> Ryder Whittaker Hawkins
> MSc Candidate
> PI: Dr. Luc Vallières
> ryder.whittaker at crchudequebec.ulaval.ca
> Work: (418) 525-4444
> Cell: (613) 979-7729
> Dept. of Molecular Medicine
> Université Laval (CHUL)
> 2705 Laurier Blvd
> Québec, QC
> G1V 4G2
> De : Tom Goddard
> Date d'envoi : 21 novembre 2017 12:37
> À : chimera List
> Cc : RYDER WHITTAKER HAWKINS
> Objet : Re: [Chimera-users] Coulombic (Mailing list question)
> If you wanted to try molecular surfaces I think the perframe script would include the surface command
> coulombic ...
> As Elaine points out if you leave out the surface command the surface does automatically update but only after the coloring which erases the coloring. I'm not sure if the surface command above will do the trick, maybe it will think the surface is already there and do nothing so you may need to delete the surface and recreate it in the perframe script.
>> On Nov 21, 2017, at 9:28 AM, Elaine Meng wrote:
>> Hello Ryder,
>> Your Answer of the Day:
>> The error is probably “No molecular surfaces selected with atom spec” and you would get that error when you haven’t shown a surface, since the Coulombic coloring applies only to surfaces.
>> (Also, for “classical” coloring by electrostatic potential, you probably want white between the red and the blue.)
>> So, the short answer may be to show a surface first (e.g. command “surface”). However, this does not work well with trajectories, for two reasons. The molecular surface automatically recalculates when coordinates change, and when I put “coulombic” into a per-frame script it is poorly synchronized so that the surface update wipes out the Coulombic coloring so that it’s never visible, at least in my test system. Secondarily, molecular surface calculation is prone to failures and usually at least a few of the trajectory frames have a surface-calculation error or crazy-looking surface. However, you could try it on your system, e.g. show surface and then apply a per-frame script with something like the following command (if your molecular surface is #0, can check in the Model Panel):
>> coulombic -10 red 0 white 10 blue #0
>> Assuming you have the same prohibitive problems I observed, here is one possible workaround but it is more effort in deciding waht parameters to use. At each frame calculate a molmap surface (isosurface of density simulated from atomic positions), adjust isosurface level so it is somewhat similar to the molecular surface shape, apply coulombic coloring. E.g. per-frame script something like the following, but with variations depending on your system:
>> molmap protein 4 gridSpacing 0.25 mod #2
>> volume #2 level 0.2
>> coulombic atoms protein -10 red 0 white 10 blue #2
>> I decided on the “level” value of 0.2 by showing a molecular surface along with the molmap surface, changing the molmap surface to mesh representation (in the Volume Viewer dialog that automatically appears when you use molmap) and dragging the bar on the Volume Viewer histogram so that the mesh surface was in a similar location to the molecular surface. The value is shown in the “Level” field in the Volume Viewer dialog. This only needs to be done once, i.e. the isosurface level that's reasonable for one frame should be equally reasonable for the other frames. The two kinds of surfaces will never be exactly the same since they are calculated differently, but you can shoot for similar enough. Then you would close (delete) the molecular surface and change the molmap surface back to “surface” style in Volume Viewer before proceeding with your trajectory playback with the per-frame script.
>> See the help pages of the commands for details on the options, but in a nutshell the molmap command simulates a density map from the protein atoms (you could specify the atoms differently, I just mean my example uses “protein”) and makes that map a new model #2. The volume command sets the isosurface level of model #2 to 0.2. The coulombic command colors the molmap surface #2 by the potential calculated from the protein atoms.
>> I hope this helps,
>> Elaine C. Meng, Ph.D.
>> UCSF Chimera(X) team
>> Department of Pharmaceutical Chemistry
>> University of California, San Francisco
>>> On Nov 21, 2017, at 7:41 AM, RYDER WHITTAKER HAWKINS wrote:
>>> My Question of the Day:
>>> - I am running an MD movie in Chimera
>>> - I want to display the solvent excluded surface colored by coulombic potential.
>>> - I would like to implement this command as a script at each frame.
>>> What is the command? I am getting error after error...
>>>> coulombic -10 red 10 blue
>>> How could I put this in a Per-Frame analysis?
>>> Thanks for any suggestions. The documentation is not so clear for this coulombic function.
>>> Ryder Whittaker Hawkins
>>> MSc Candidate
>>> PI: Dr. Luc Vallières
>>> ryder.whittaker at crchudequebec.ulaval.ca
>>> Work: (418) 525-4444
>>> Cell: (613) 979-7729
>>> Dept. of Molecular Medicine
>>> Université Laval (CHUL)
>>> 2705 Laurier Blvd
>>> Québec, QC
>>> G1V 4G2
>> Chimera-users mailing list: Chimera-users at cgl.ucsf.edu
>> Manage subscription: http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users
> Chimera-users mailing list: Chimera-users at cgl.ucsf.edu
> Manage subscription: http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users
More information about the Chimera-users