Hi David, I’ve attached a Python script mean_distance.py that can compute the distance of each atom in an ensemble from the mean structure. Then you can color using that distance. Here’s an example set of Chimera commands. open 1mtx show open mean_distance.py rangecolor distance_from_mean 0 blue 3 white 6 red The script uses all open molecules as the ensemble. The coloring can also be done with menu Tools / Depiction / Render by Attribute. Tom
On May 3, 2018, at 10:10 AM, Morgan, David Gene wrote:
Tom,
Hope you are doing well. Jim Powers and I here at IU keep thinking we need to invite you and/or the team back for another chimera workshop, but we haven't done the ground work to see how possible that would really be. But I keep hoping!
In the meantime, I have a question about whether chimera can do something, and if not, I have a feature request... I was watching a grad student give a talk showing lots of superposed atomic structures. The structures included some molecular dynamics simulations and some crystal structures of the same protein in subtly different conformations (packing into a viral capsid, where the different subunits in the asymmetric unit have slightly different conformations). It occurred to me that _after_ showing the mess of spaghetti he had, the points he needed to make would be better served with a different visual: If he were to average the ensemble of MD atomic positions to a single chain (and color-code by the ensemble RMS difference from the average) and do the same for all his ensemble of X-ray structures, it would be much easier to talk about and show the consistencies and inconsistencies in the atomic co-ordinates from all the data. This same sort of averaging an ensemble of atomic co-ordinates is sometimes done with NMR structure ensembles (but I don't know about the color coding with the RMS).
I have looked at lots of packages that manipulate atomic co-ordinates, but haven't found anything that seems to actually do what I want to do. For example, gromacs is supposed to do the averaging but it really wants to work with MD data and at a minimum, I'm not smart enough to make it work with a series of standard PDB files. My NMR colleague suggested molmol for the averaging, and I need to explore that more. But there doesn't seem to be a recent linux version nor does it seem to be "scriptable" (i.e., if I wanted to do this for a dozen different ensembles, I would have to painfully mouse click my way through the process). On top of all that, I ran into numerous on-line discussions about how averaging was "the wrong thing to do." All those people had valid reasons for saying so, but right now, I'm really talking about a more simple visual that conveys some of the information contained in any sort of ensemble of atomic co-ordinates.
So my first question is: does chimera have a utility that would do this sort of thing? I looked at various commands but did not see anything that looked like it would do this. If not, can I put a feature request into the queue?
Looking forward to hearing from you and chatting some more.
-- politics is more difficult than physics. A. Einstein
David Gene Morgan Electron Microscopy Center 047D Simon Hall IU Bloomington 812 856 1457 (office) 812 856 3221 (3200) http://iubemcenter.indiana.edu
