<html><head><meta http-equiv="Content-Type" content="text/html charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">Hi Oliver,<div class=""><span class="Apple-tab-span" style="white-space:pre"> </span>And you can resort to a fairly simple Python script to re-establish bonds between the CAs after deleting all the other atoms:</div><div class=""><br class=""></div><div class="">from chimera import openModels, Molecule</div><div class="">for m in openModels.list(modelTypes=[Molecule]):</div><div class=""><span class="Apple-tab-span" style="white-space:pre"> </span>prev_res = None</div><div class=""><span class="Apple-tab-span" style="white-space:pre"> </span>for r in m.residues:</div><div class=""><span class="Apple-tab-span" style="white-space:pre"> </span>if prev_res == None:</div><div class=""><span class="Apple-tab-span" style="white-space:pre"> </span>prev_res = r</div><div class=""><span class="Apple-tab-span" style="white-space:pre"> </span>continue</div><div class=""><span class="Apple-tab-span" style="white-space:pre"> </span>if prev_res.id.chainId == r.chainId:</div><div class=""><span class="Apple-tab-span" style="white-space:pre"> </span>r.molecule.addBond(prev_res.atoms[0], r.atoms[0])</div><div class=""><span class="Apple-tab-span" style="white-space:pre"> </span>prev_res = r</div><div class=""><br class=""></div><div class="">—Eric</div><div class=""><br class=""><div><blockquote type="cite" class=""><div class="">On Feb 1, 2016, at 12:01 PM, Tom Goddard <<a href="mailto:goddard@sonic.net" class="">goddard@sonic.net</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><meta http-equiv="Content-Type" content="text/html charset=utf-8" class=""><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">Hi Oliver,<div class=""><br class=""></div><div class=""> Of course you can keep the entire backbone (command “delete ~@CA,C,N,O”) and then you keep your ribbon, but it is 4 times more atoms than CA only.</div><div class=""><br class=""></div><div class=""><span class="Apple-tab-span" style="white-space:pre"> </span>Tom</div><div class=""><br class=""><div class=""><blockquote type="cite" class=""><div class="">On Feb 1, 2016, at 11:57 AM, Tom Goddard wrote:</div><br class="Apple-interchange-newline"><div class=""><meta http-equiv="Content-Type" content="text/html charset=utf-8" class=""><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div class="">Hi Oliver,</div><div class=""><br class=""></div><div class=""> If you delete everything except CA (command “delete ~@CA”) you lose the ribbons. This is because Chimera needs a fake bond between consecutive CA atoms to make the ribbon. If Chimera reads a CA only PDB file it makes those fake bonds. But I don’t know another way to make those bonds. I’ve wanted this capability, requested it 4 years ago</div><div class=""><br class=""></div><span class="Apple-tab-span" style="white-space:pre"> </span><a href="http://plato.cgl.ucsf.edu/trac/chimera/ticket/9875" class="">http://plato.cgl.ucsf.edu/trac/chimera/ticket/9875</a><div class=""><br class=""></div><div class="">but think it will have to wait for ChimeraX. You could make a Python script that added a Chimera command that did all the steps you are doing now.</div><div class=""><br class=""></div><div class=""><span class="Apple-tab-span" style="white-space:pre"> </span>Tom</div><div class=""><br class=""></div><div class=""><br class=""><div class=""><blockquote type="cite" class=""><div class="">On Jan 30, 2016, at 1:44 PM, Oliver Clarke <<a href="mailto:olibclarke@gmail.com" class="">olibclarke@gmail.com</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div dir="ltr" class=""><div class="">Thanks for the helpful suggestions Elaine! I will look into them and see where it gets me - I have a feeling I have tried the perframe/scolor idea in the past without success but I may be misremembering.<br class=""><br class=""></div><div class="">One other thing - I'd like a one command way to convert a full atom model to a Calpha model, keeping the helix/sheet records at the start intact. This would be very helpful when running morphs of a large protein (where often I do not need to display sidechains, just ribbon), because it makes the calculations much faster (e.g. rather than 120,000 atoms I have ~17000 using CAs only). Currently what I do is I delete all the non-CA atoms in chimera, save the pdb, remove all the TER flags chimera has added between each CA atom using grep, and then replace the HELIX/SHEET records with those from the original file (because they get corrupted for some reason in the process. This works fine but it would be nice to have a one liner in chimera to handle the same at some point (although I realize it is a pretty niche functionality and maybe not worth it). (alternatively a "coarse grained" morph where parameters for each residue only take into account the CA would accomplish the same thing I guess.)<br class=""><br class=""></div><div class="">Cheers,<br class=""></div><div class="">Oli<br class=""></div></div><div class="gmail_extra"><br class=""><div class="gmail_quote">On Sat, Jan 30, 2016 at 9:23 PM, Elaine Meng <span dir="ltr" class=""><<a href="mailto:meng@cgl.ucsf.edu" target="_blank" class="">meng@cgl.ucsf.edu</a>></span> wrote:<br class=""><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Oliver,<br class="">
I haven’t tried these things in a situation similar to yours, but some initial ideas are:<br class="">
<br class="">
(1) use two “vop morph” commands separated by only a semicolon so that (if it is possible to run two map-morphs concurrently), they will start at the same frame<br class="">
<br class="">
(2) use “perframe” (possibly with “alias”) to specify running “scolor" zone command(s) at each frame<br class="">
<br class="">
I would try those things one at a time, for example, perframe scolor with a single vop morph. Even in that case, I’m not sure it will work. We may be able to come up with better answer or at least a more definite answer of whether it is possible next week. In the meanwhile, these are just some things to look into.<br class="">
Best,<br class="">
Elaine<br class="">
-----<br class="">
Elaine C. Meng, Ph.D.<br class="">
UCSF Computer Graphics Lab (Chimera team) and Babbitt Lab<br class="">
Department of Pharmaceutical Chemistry<br class="">
University of California, San Francisco<br class="">
<div class="HOEnZb"><div class="h5"><br class="">
On Jan 30, 2016, at 10:16 AM, Oliver Clarke <<a href="mailto:olibclarke@gmail.com" class="">olibclarke@gmail.com</a>> wrote:<br class="">
<br class="">
> Hello,<br class="">
><br class="">
> In a situation where I have maps for two conformations, and I've extracted out domain specific maps in both cases, is it possible to run morph maps on both groups of maps simultaneously? Or is there any way to apply color zone to the trajectory generated by morph maps? I would like to generate a morph map movie where the map is colored by domain throughout the morph, and can't figure out how to do so.<br class="">
><br class="">
> Cheers,<br class="">
> Oli<br class="">
<br class="">
</div></div></blockquote></div><br class=""></div>
_______________________________________________<br class="">Chimera-users mailing list<br class=""><a href="mailto:Chimera-users@cgl.ucsf.edu" class="">Chimera-users@cgl.ucsf.edu</a><br class=""><a href="http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users" class="">http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users</a><br class=""></div></blockquote></div><br class=""></div></div>_______________________________________________<br class="">Chimera-users mailing list<br class=""><a href="mailto:Chimera-users@cgl.ucsf.edu" class="">Chimera-users@cgl.ucsf.edu</a><br class=""><a href="http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users" class="">http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users</a><br class=""></div></blockquote></div><br class=""></div></div>_______________________________________________<br class="">Chimera-users mailing list<br class=""><a href="mailto:Chimera-users@cgl.ucsf.edu" class="">Chimera-users@cgl.ucsf.edu</a><br class="">http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users<br class=""></div></blockquote></div><br class=""></div></body></html>