<html><head><meta http-equiv="Content-Type" content="text/html; charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">Dear Elaine,<div class="">I'm trying to go a little further with the isosurfaces case… I would like to make a movie showing the movements of some lipids when they interact with the protein. The way you suggested me before works perfectly but now I would like to incorporate the perframe script commands into the first script that I designed, the .com script. And here is where I find find the problem because I can’t move the protein and the lipids as I did with the MD movie window. I am trying to combine both commands, perframe and coordset (I think this is the way to do it with commands) but I only obtain errors. Is it possible to do the protein and lipids movement part of the movie with commands instead of the MD movie windows?</div><div class=""><br class=""></div><div class="">Here is the .com script:</div><div class=""><br class=""></div><div class=""><div class=""><font color="#0432ff" class="">#protein in isosurface rep</font></div><div class="">movie record</div><div class="">molmap protein 4 grid 0.6 model 1</div><div class="">volume #1 color #dae4bceb8d98 </div><div class="">wait</div><div class="">~ribbon :.a </div><div class="">movie crossfade 60</div><div class=""><br class=""></div><div class=""><font color="#0432ff" class="">#PC and PA phosphates in sphere rep</font></div><div class="">repr sphere :DVPC & P</div><div class="">color #aaaa84bdc71c :DVPC & P</div><div class="">movie crossfade 30</div><div class="">2dlabel create l1 text 'Phosphatidylcholine (PC)' color #84bd4bdaaaaa ypos 0.9 xpos .7 size 68 typeface serif </div><div class="">wait 30</div><div class=""><br class=""></div><div class="">repr sphere :DVPA & P</div><div class="">movie crossfade 30</div><div class="">color #aaaab425684b :DVPA & P</div><div class="">movie crossfade 30</div><div class="">2dlabel create l2 text 'Phosphatidic Acid (PA)' color #4bda71c725ed ypos 0.9 xpos .7 size 68 typeface serif </div><div class="">movie crossfade 30</div><div class=""><br class=""></div><div class=""><font color="#0432ff" class="">#hide tails</font></div><div class="">~disp :DVPC & C </div><div class="">~disp :DVPC & H </div><div class="">~disp :DVPC & O </div><div class="">~disp :DVPC & N </div><div class=""><br class=""></div><div class="">~disp :DVPA & C </div><div class="">~disp :DVPA & H </div><div class="">~disp :DVPA & O</div><div class=""><br class=""></div><div class=""><font color="#0432ff" class="">#select DVPC and DVPA at z > 5 from the protein to make them transparent</font></div><div class="">select :.A z > 5 & :DVPC</div><div class="">transp 50,a :DVPC</div><div class="">select :.A z > 5 & :DVPA</div><div class="">transp 50,a :DVPA</div><div class="">movie crossfade 30</div><div class=""><br class=""></div><div class="">2dlabel delete *</div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><font color="#0432ff" class="">#select DVPC and DVPA lipids at z < 5 from the protein to create isosurfaces</font></div><div class="">molmap ":.A z<5 & :DVPC" 4 grid 0.8 model 2</div><div class="">~ disp element.H </div><div class="">volume #2 color #aaaa84bdc71c</div><div class="">movie crossfade 100</div><div class=""><br class=""></div><div class="">molmap ":.A z<5 & :DVPA" 4 grid 0.8 model 3</div><div class="">~ disp element.H </div><div class="">volume #3 color #aaaab425684b</div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><br class=""></div><div class="">At this point this is the representation that I obtain: </div><div class=""><br class=""></div><div class=""><img apple-inline="yes" id="A149130A-C3C5-4F7F-92BA-E28C5900C9EF" width="320" height="153" src="cid:769AECB8-B0A9-4FCA-92CE-AA776F2272CB" class=""></div></div><div class=""><br class=""></div><div class=""><br class=""></div><div class="">And now it is the time for the trajectories movements...</div><div class="">What I have been trying is:</div><div class="">perframe coordset #1-3 20,50,1 </div><div class="">or perframe "coordset #1-3 20,50,1” and other different variants but I only obtain errors.</div><div class=""><br class=""></div><div class="">Could you help me with this issue, please?,</div><div class="">clara</div><div class=""><br class=""></div><div class=""> </div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><div><br class=""><blockquote type="cite" class=""><div class="">El 22 abr 2019, a las 23:15, Elaine Meng <<a href="mailto:meng@cgl.ucsf.edu" class="">meng@cgl.ucsf.edu</a>> escribió:</div><br class="Apple-interchange-newline"><div class=""><div class="">Hi Clara<br class="">(oops sorry for using the wrong name before),<br class=""><br class="">You can include coloring in the script, e.g.<br class=""><br class="">volume #1 color blue<br class="">volume #2 color red<br class=""><br class="">Also to avoid selecting at every frame, you may be able to specify atoms directly in the molmap command, e.g.<br class=""><br class="">molmap ":.A z<5 & :DVPC” 4 grid 0.8 model 2<br class=""><br class="">The quotation marks should be plain-text quotation marks, not the fancy curved ones that may be shown in this message because the Mail app changes them automatically.<br class=""><br class="">I hope this helps,<br class="">Elaine<br class="">-----<br class="">Elaine C. Meng, Ph.D.<br class="">UCSF Chimera(X) team<br class="">Department of Pharmaceutical Chemistry<br class="">University of California, San Francisco<br class=""><br class=""><blockquote type="cite" class="">On Apr 22, 2019, at 10:43 AM, MARIA CLARA BLANES MIRA <<a href="mailto:c.blanes@goumh.umh.es" class="">c.blanes@goumh.umh.es</a>> wrote:<br class=""><br class="">Hello Elaine,<br class="">Just what I needed. Thank you very much. Furthermore, the graphics are great!.<br class=""><br class="">The only thing I can’t control is the color because they are changing in a rainbow way when I click on the play in MD Movie. I want the protein always for example blue and the the lipids near the protein always in red. I am trying different options in volume color zones but I can’t find the way.<br class="">This is what I introduced in the per-frame script:<br class=""><br class="">molmap protein 4 grid 0.6 model 1<br class="">select :.A z<5 & :DVPC <br class="">molmap sel 4 grid 0.8 model 2<br class=""><br class="">And it is great for me because they move but it remains the problem with the colours. So, is there any way to fix the colours?<br class="">Thank you very much for your attention!!,<br class="">clara<br class=""><br class=""><blockquote type="cite" class="">El 22 abr 2019, a las 18:22, Elaine Meng <<a href="mailto:meng@cgl.ucsf.edu" class="">meng@cgl.ucsf.edu</a>> escribió:<br class=""><br class="">Hi Maria,<br class="">I see that VMD Quicksurf makes a density map from the atoms and shows an isosurface. You can do the same thing in Chimera with the “molmap” command, but the density map is not automatically recalculated at every step of the trajectory. You would have to use a per-frame script and re-execute molmap at each frame to replace the old map with a new map every time the coordinates change.<br class=""><br class="">The molmap command has a resolution parameter and a grid spacing parameter just like Quicksurf.<br class=""><<a href="http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/molmap.html" class="">http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/molmap.html</a>><br class=""><br class="">You can test by showing one frame of your trajectory and trying various molmap resolution values, e.g. 2.5 angstroms with separate surfaces for protein and non-protein atoms:<br class=""><br class="">molmap protein 2.5 grid .5 model 1<br class="">molmap ~protein 2.5 grid .5 model 2<br class=""><br class="">Another issue is that the isosurface level isn’t controlled in the molmap command but could be controlled in a separate volume command, e.g.:<br class=""><br class="">volume #1 level 0.1<br class="">vol #2 level 0.13<br class=""><br class="">If you found commands that worked well enough for you, and you are using MD Movie to view the trajectory, in MD Movie menu: Per-Frame… Define Script and enter these commands as Chimera command script to be executed at each frame.<br class=""><<a href="http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/movie/framemovie.html" class="">http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/movie/framemovie.html</a>><br class=""><<a href="http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/movie/movie.html#per-frame" class="">http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/movie/movie.html#per-frame</a>><br class=""><br class="">I hope this helps,<br class="">Elaine<br class="">-----<br class="">Elaine C. Meng, Ph.D.<br class="">UCSF Chimera(X) team<br class="">Department of Pharmaceutical Chemistry<br class="">University of California, San Francisco<br class=""><br class=""><blockquote type="cite" class="">On Apr 21, 2019, at 11:57 AM, MARIA CLARA BLANES MIRA <<a href="mailto:c.blanes@goumh.umh.es" class="">c.blanes@goumh.umh.es</a>> wrote:<br class=""><br class="">Hello,<br class="">I would like to know if Chimera could show the surfaces like VMD does with the command QuickSurf (I attach the image with the protein in yellow and the lipids in green). <br class="">I have been trying with a lot of different variations of the surface command but I can not find the way to make the surface ’softer'….<br class="">I have also been trying with spheres representation but the image shows too much information. I need a simplest view, similar to the image because this way it is easier see them move in the .dcd file. It not so important for the protein, just for the lipids.<br class=""><br class="">Could you help me with this, please?<br class="">Thank you very much,<br class="">clara<br class=""></blockquote></blockquote><br class=""></blockquote><br class=""></div></div></blockquote></div><br class=""></div></body></html>