- Chimera-users - RBVI Mailing Lists

3D to 2D
by Sameer K Singh 30 Sep '16

30 Sep '16

Hi, I was wondering if it's possible to extract information from 3D plots in Chimera and render them as 2d plots.Specifically I'm trying to extract information from a .bild file displayed in chimera and display it as a 2D plot of some kind in an external program (eg: excel) Thank you very much. Best, Sameer ________________ Sameer Singh, PhD Albert Einstein College of Medicine 1300 Morris Park Avenue New York, NY - 10461 Phone: 718-430-2048 Cell: 917-848-9319

2 1

differentiating lower and upper case chain ID's
by Melissa Jurica 29 Sep '16

29 Sep '16

Hi- Is it possible to differentiate between chain ID’s that are in upper vs. lower case (i.e. A vs. a) in the command line? Spliceosome complexes have too many chains (ex. 5GM6). Thanks, Melissa

2 4

Surface area list
by esserlo＠helix.nih.gov 28 Sep '16

28 Sep '16

Hi, with the help of "split #0" I created 10 segments of my protein and calculated the surface of each fragment. I can go to tools and measure volume and area but I'd prefer to write out the results as a file. Can Chimera do this ? I am asking because I intend to do more complex things than just 10 surfaces and it would be rather tedious if not impossible to do all this by hand. Thanks for any advice. Lothar

3 3

Problem with running Chimera
by Patty Tran 27 Sep '16

27 Sep '16

Hello! I installed two versions of Chimera (Microsoft Windows 64-bit chimera-1.11.2-win64.exe <https://www.cgl.ucsf.edu/chimera/cgi-bin/secure/chimera-get.py?file=win64/c…> and Microsoft Windows 64-bit chimera-1.11-win64.exe) <https://www.cgl.ucsf.edu/chimera/cgi-bin/secure/chimera-get.py?file=win64/c…> I have only been able to run the program once and that was when I clicked "start program" after installing the 1.11 version. I have been unable to find the Chimera desktop icons and the program itself in the list of programs. I thought I could find the programs under my C:/Program Files/Chimera but there are no folders. Could someone provide instruction on how to uninstall the programs and install them so I could find and use the program easily? Thank you. Patricia

2 1

Measuring angles and distances between planes defined by two aromatic rings
by George Tzotzos 26 Sep '16

26 Sep '16

Having generated a 100ns trajectory, I’m interested in measuring the time evolution of the distance between to aromatic ring centroids as well as the angle between the planes of these two rings. The idea is to determine if I’m dealing with pi-pi stacking interactions. I’ve check past correspondence and tried various permutations of commands to no avail. Maybe my best option is to do this for selected snapshots and define for each snapshot two centroids (e.g. c1 and c2) and measure the distance c1 c2. I haven’t managed to find a solution for measuring the angles between the two planes defined by the ring atoms. I’d appreciate any suggestions on possible solutions to the problem described above. I’m attaching a snapshot of measuring distances between two atoms (not really a good solution). Thanks in advance George

2 8

Interpolate colors in Morph after Color zone
by Haselbach, David 26 Sep '16

26 Sep '16

Hi everyone, Want to make a movie of a movement of a single domain in my complex through morphing through several cryo EM densities. However as it is very complex I wanted to color the domain of interest. I did so with color zone and got a nice result for all the intermediate states. However as soon as I start morphing the color is gone even though interpolate colors is on. Is there are trick to achieve that? I tried segmenting the map and then interpolating through the segmented maps which works in principle fine but I don't really like the edge that is created by the segmentation. Kind regards David Dr. David Haselbach Max-Planck-Institute for biophysical Chemistry Department for structural Dynamics Am Fassberg 11 37077 Germany Tel. +495512011302

3 2

segmentation and resolution
by Kristin Parent 23 Sep '16

23 Sep '16

To whom it may concern, I am currently working on segmenting a high resolution cryo-EM map (~4Å by the ³gold standard FSC²) to phase crystallography data. After the segmentation, my colleague doing the phasing tells me that map appears more like 12Å resolution and I am worried I am doing something wrong during the process and somehow distorting my original data. I believe the 4Å of the original map based on structural features that I can clearly see. And the segmented density does not appear as defined. Therefore it appears that there has been some changes. In order to begin the process, I did some manipulations before I started. This is an icosahedral virus map created through Auto3dem. I am only using an small quadrant of the original map, and since it was in pif format originally, I used proc3d to convert to mrc and invert the contrast to be compatible with segger. Before segmentation the map visually looked the same as the original pif map when both were displayed in chimera. The procedure I followed after segmentation was: "Save selected regions to .mrc file... save density map masked by the selected <https://www.cgl.ucsf.edu/chimera/docs/UsersGuide/selection.html> segmentation regions to an MRC file (map dimensions set to the minimal box containing the regions)². I have made sure that the correct binning has been applied during segmentation. Any help would be much appreciated. Kristin -- Kristin N. Parent, Ph.D. Assistant Professor Dept. of Biochemistry and Molecular Biology Michigan State University 603 Wilson Road Room 519 Biochemistry Building East Lansing, MI 48824 Office phone: (517) 432-8434 https://kparentlab.natsci.msu.edu/

2 1

Re: [Chimera-users] Recursive structure matching and acquisition of descriptive numbers
by Jaime Rodríguez-Guerra 23 Sep '16

23 Sep '16

Hi! Since you are getting lists of atoms back, you only need to do some attribute access (rather than method calling)! Every chimera.Atom object has an attribute called molecule: a reference to its parent chimera.Molecule. chimera.Molecule objects have an attribute called 'name', but sometimes this is not very informative (depends on your input), so you need to resort to the pdb filename, stored in the first element of the openedAs tuple. In code, this looks like this: for atoms1, atoms2, rmsd, fullRmsd in match(CP_BEST, [ref, sims], defaults[MATRIX], "nw", defaults[GAP_OPEN], defaults[GAP_EXTEND]): molecule1 = atoms1[0].molecule # any atom from the list will do molecule2 = atoms2[0].molecule print molecule1.name, "\t", molecule2.name, "\t", rmsd # molecule1.openedAs[0] will also work here I don't know if atoms1 comes consistently from the reference molecule, but I'd guess it does. Hope it helps!

3 4

Saving scenes
by Tom Goddard 22 Sep '16

22 Sep '16

Hi Josh, Chimera session files (menu File / Save Session) save the data, views, colors, display styles, what is hidden and shown. The savepos and reset commands save and restore the positions of models and clip planes and view direction and zoom, just in memory (e.g. “savepos view1”, “reset view1”). Positions saved this way are also saved in sessions so after reloading a session you can use the reset command to get those positions back. Savepos does not save the colors, display styles, or what is hidden and shown. “Scenes” are created by the “Animation” tool (menu Tools / Utilities / Animation) and the “scene” command and include positions, camera view, colors, display styles, what is hidden and shown. Scenes try to capture exactly how the scene looks so they are like sessions but no file is used (although the are preserved in session files). So to save the current scene use command “scene s1 save” and to restore it “scene s1 reset”. Like savepos/reset you can save as many scenes as you want with different names. Scenes are a more recent feature and are not completely reliable — they don’t record everything and sometimes the cause errors saving sessions. Tom > On Sep 22, 2016, at 5:38 AM, Broyde, Joshua E. wrote: > > Dear Tom, > I hope you are well. I was wondering how to save a particular view (i.e. with specified zoom, atoms displayed and angle of view) in Chimera? Currently, I simply save this as a different session, but I was wondering if there is a better or easier way? > > Josh >

1 0

Re: [Chimera-users] Chimera-users Digest, Vol 161, Issue 19
by Healey, Joe 22 Sep '16

22 Sep '16

Ah yes now I see. mmaker is what I want. I've figured out the appropriate syntax for the command I'd want to execute in chimera syntax I think. In the code snippet below, it seems to me I have 2 options? I can synthesise a string and iterate the mmaker function in chimera syntax with something like: chimera.openModels.open(top_hit, type="PDB") # Open model structures for model_path in model_list: chimera.openModels.open(model_path,type="PDB") print("Opened: " + model_path) rmsd_file = '{0}.tsv'.format(rmsd) with open(rmsd_file, "a") as rmsd_output_file: for j in len(model_list): match_cmd = 'mmaker #0 #{0} alg nw iter false'.format(j) runCommand(match_cmd) # <somewhere in here retrieve the rmsd and append to a result file> or the other option would be to do the above but avoiding the string synthesis step with something like: chimera.module.mmaker() If the module can be accessed directly? If so which module is the mmaker function within? Though I still don't know how to access the RMSD value from the mmaker function in either of these manners. Would executing mmaker in either such way return a tuple or list where the RMSD is a particular index? Looking at your Match script the rmsd was the 3rd value of the tuple in the line: atoms1, atoms2, rmsd = match(CP_SPECIFIC_SPECIFIC, [(c1, c2)], defaults[MATRIX],"nw", defaults[GAP_OPEN], defaults[GAP_EXTEND],iterate=defaults[ITER_CUTOFF])[0] Would the same be true for mmaker? (Sorry to ask so many questions in one go. I don't know what your policy is on whether all related questions should be in the same issue, or whether I should be raising separate ones) Thanks again for all your help Joe Healey M.Sc. B.Sc. (Hons) PhD Student MOAC CDT, Senate House University of Warwick Coventry CV47AL Mob: +44 (0) 7536 042620 | Email: J.R.J.Healey(a)warwick.ac.uk Jointly working in: Waterfield Lab<http://www2.warwick.ac.uk/fac/med/research/tsm/microinfect/staff/waterfield…> (WMS Microbiology and Infection Unit) and the Gibson Lab<http://www2.warwick.ac.uk/fac/sci/chemistry/research/gibson/gibsongroup/> (Warwick Chemistry) Twitter: @JRJHealey<https://twitter.com/JRJHealey> | Website: MOAC Page<http://www2.warwick.ac.uk/fac/sci/moac/people/students/2013/joseph_healey> ________________________________ From: chimera-users-bounces(a)cgl.ucsf.edu <chimera-users-bounces(a)cgl.ucsf.edu> on behalf of chimera-users-request(a)cgl.ucsf.edu <chimera-users-request(a)cgl.ucsf.edu> Sent: 20 September 2016 19:30 To: chimera-users(a)cgl.ucsf.edu Subject: Chimera-users Digest, Vol 161, Issue 19 Send Chimera-users mailing list submissions to chimera-users(a)cgl.ucsf.edu To subscribe or unsubscribe via the World Wide Web, visit http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users or, via email, send a message with subject or body 'help' to chimera-users-request(a)cgl.ucsf.edu You can reach the person managing the list at chimera-users-owner(a)cgl.ucsf.edu When replying, please edit your Subject line so it is more specific than "Re: Contents of Chimera-users digest..." Today's Topics: 1. Re: Recursive structure matching and acquisition of descriptive numbers (Eric Pettersen) ---------------------------------------------------------------------- Message: 1 Date: Tue, 20 Sep 2016 11:30:38 -0700 From: Eric Pettersen <pett(a)cgl.ucsf.edu> To: "Healey, Joe" <J.R.J.Healey(a)warwick.ac.uk> Cc: "chimera-users(a)cgl.ucsf.edu BB" <chimera-users(a)cgl.ucsf.edu> Subject: Re: [Chimera-users] Recursive structure matching and acquisition of descriptive numbers Message-ID: <30002D29-1B25-4596-8868-EB3370B72E42(a)cgl.ucsf.edu> Content-Type: text/plain; charset="windows-1252" Hi Joe, I guess the first thing is that the command equivalent for MatchMaker is ?matchmaker? or ?mmaker?, not ?match?. The ?match? command is for when you know a priori the exact set of atoms you want matched in the two structures. ?matchmaker? is for when you let Chimera figure out the set of atoms to use based on the sequences and secondary structure. The documentation for the matchmaker command is here: https://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/matchmaker.html#opti… <https://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/matchmaker.html#opti…> I looked at your script and didn?t see anything obviously wrong. What exactly isn?t working? And unfortunately the answer is no, you can?t run a GUI operation and easily find the command equivalent. ChimeraX will have that capability though. ?Eric > On Sep 20, 2016, at 5:07 AM, Healey, Joe <J.R.J.Healey(a)warwick.ac.uk> wrote: > > Hi again, > > I'm making decent headway with the problem I mentioned before with Jaime's help and the use of pychimera. I'm down to the nitty-gritty of doing the actual RMSD/Match function alignment though and could use a bit more information. > > Is there are resource with more information on the commandline implementation of match anywhere other than the link below? > > https://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/match.html#options <https://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/match.html#options> > > I took a look at Eric's suggestion of the match.py script here: http://plato.cgl.ucsf.edu/pipermail/chimera-users/2009-March/003610.html <http://plato.cgl.ucsf.edu/pipermail/chimera-users/2009-March/003610.html> > > But I can't quite unwrap the options specified in the command for my uses. > > Perhaps another way to ask this is, can I run a command in GUI chimera, and find the EXACT command that was executed? > > I've put the semi-complete script in pastebin here http://pastebin.com/6k9cew55 <http://pastebin.com/6k9cew55> incase it gives you some indication of what I'm attempting to do. > > Many thanks > > > Joe Healey > > > M.Sc. B.Sc. (Hons) > PhD Student > MOAC CDT, Senate House > University of Warwick > Coventry > CV47AL > Mob: +44 (0) 7536 042620 | Email: J.R.J.Healey(a)warwick.ac.uk <mailto:J.R.J.Healey@warwick.ac.uk> > > Jointly working in: > Waterfield Lab <http://www2.warwick.ac.uk/fac/med/research/tsm/microinfect/staff/waterfield…> (WMS Microbiology and Infection Unit) > and the Gibson Lab <http://www2.warwick.ac.uk/fac/sci/chemistry/research/gibson/gibsongroup/> (Warwick Chemistry) > > Twitter: @JRJHealey <https://twitter.com/JRJHealey> | Website: MOAC Page <http://www2.warwick.ac.uk/fac/sci/moac/people/students/2013/joseph_healey> > From: Greg Couch <gregc(a)cgl.ucsf.edu <mailto:gregc@cgl.ucsf.edu>> > Sent: 08 September 2016 22:43:55 > To: Healey, Joe; chimera-users(a)cgl.ucsf.edu <mailto:chimera-users@cgl.ucsf.edu> BB > Subject: Re: [Chimera-users] Recursive structure matching and acquisition of descriptive numbers > > chimera --nogui can't access graphics. Headless chimera has software OpenGL rendering, so it can be used to create images. Make sure to use the windowsize command to set the size of the framebuffer. Don't know if pychimera cares. > > -- Greg > > On 09/08/2016 11:53 AM, Healey, Joe wrote: >> Brilliant, that sounds like the most likely option. I had stumbled across pychimera a couple of days ago while trying to find out if chimera could simply be imported as a module. I initially steered away from it as I don't have sudo rights for the server where I'm trying to set all this up (couldnt run pip), but I've since sorted it with our sysadmin. >> >> On a related point, is running chimera with --nogui equivalent to running headless chimera? I'm just wondering if pychimera is happy with either >> >> >> Joe Healey >> >> >> M.Sc. B.Sc. (Hons) >> PhD Student >> MOAC CDT, Senate House >> University of Warwick >> Coventry >> CV47AL >> Mob: +44 (0) 7536 042620 | Email: J.R.J.Healey(a)warwick.ac.uk <mailto:J.R.J.Healey@warwick.ac.uk> >> >> Jointly working in: >> Waterfield Lab <http://www2.warwick.ac.uk/fac/med/research/tsm/microinfect/staff/waterfield…> (WMS Microbiology and Infection Unit) >> and the Gibson Lab <http://www2.warwick.ac.uk/fac/sci/chemistry/research/gibson/gibsongroup/> (Warwick Chemistry) >> >> Twitter: @JRJHealey <https://twitter.com/JRJHealey> | Website: MOAC Page <http://www2.warwick.ac.uk/fac/sci/moac/people/students/2013/joseph_healey> >> From: Greg Couch <gregc(a)cgl.ucsf.edu> <mailto:gregc@cgl.ucsf.edu> >> Sent: 07 September 2016 19:49:19 >> To: Healey, Joe; chimera-users(a)cgl.ucsf.edu <mailto:chimera-users@cgl.ucsf.edu> BB >> Subject: Re: [Chimera-users] Recursive structure matching and acquisition of descriptive numbers >> >> Yes, Chimera's Python code is only "warrantied" to work with its own Python. See question 3b of the Chimera Programming FAQ,https://www.cgl.ucsf.edu/chimera/docs/ProgrammersGuide/faq.html#q3b <https://www.cgl.ucsf.edu/chimera/docs/ProgrammersGuide/faq.html#q3b>, for how to install other software into Chimera. Other parts of the FAQ may help you too. >> >> That said, Jaime Rodr?guez-Guerra took my previous comments on this subject, http://www.cgl.ucsf.edu/pipermail/chimera-users/2015-January/010647.html <http://www.cgl.ucsf.edu/pipermail/chimera-users/2015-January/010647.html>, and came up with pychimera, https://pypi.python.org/pypi/pychimera <https://pypi.python.org/pypi/pychimera>. I haven't checked to see how it handles all of my caveats, and it can't completely without using a modified Python, but I suspect that it would be good enough for you. >> >> HTH, >> >> Greg >> >> On 09/07/2016 07:36 AM, Healey, Joe wrote: >>> Hi Eric, >>> >>> That looks good thank you. I think I'm going about this slightly the wrong way after some investigation however. >>> >>> Am I right in thinking that chimera can't be imported in to python inside a normal bash shell, and rather it must be done the other way round? So file manipulation and such must be done by calls to the system from within the chimera python shell, rather than a python interpreter calling chimera modules in? >>> >>> It certainly makes the chimera steps I plan to do much easier, but I was going to include it in a pipeline with some other programs and it makes that a little more taxing (though not insurmountable!). >>> >>> Additionally, I plan to invoke this inside a script on a headless linux box (hence why running it through a python script was ideal), so whatever solution I use will need to be 'command line friendly'. >>> >>> In effect all I need is the actual algorithm for the matching process to be accessed somehow from a python script. >>> >>> Do you have any advice about how best to proceed? >>> >>> Many thanks, >>> >>> Joe Healey >>> >>> >>> M.Sc. B.Sc. (Hons) >>> PhD Student >>> MOAC CDT, Senate House >>> University of Warwick >>> Coventry >>> CV47AL >>> Mob: +44 (0) 7536 042620 | Email: J.R.J.Healey(a)warwick.ac.uk <mailto:J.R.J.Healey@warwick.ac.uk> >>> >>> Jointly working in: >>> Waterfield Lab <http://www2.warwick.ac.uk/fac/med/research/tsm/microinfect/staff/waterfield…> (WMS Microbiology and Infection Unit) >>> and the Gibson Lab <http://www2.warwick.ac.uk/fac/sci/chemistry/research/gibson/gibsongroup/> (Warwick Chemistry) >>> >>> Twitter: @JRJHealey <https://twitter.com/JRJHealey> | Website: MOAC Page <http://www2.warwick.ac.uk/fac/sci/moac/people/students/2013/joseph_healey> >>> >>> From: Eric Pettersen <pett(a)cgl.ucsf.edu> <mailto:pett@cgl.ucsf.edu> >>> Sent: 31 August 2016 18:23 >>> To: Healey, Joe >>> Cc: chimera-users(a)cgl.ucsf.edu <mailto:chimera-users@cgl.ucsf.edu> >>> Subject: Re: [Chimera-users] Recursive structure matching and acquisition of descriptive numbers >>> >>> Hi Joe, >>> This old chimera-users post has a Python script attachment that demonstrates how to call the underlying MatchMaker function to get the RMSD value programmatically. The mailing list gives the attachment a ?.bin? extension but it?s really a .py file. Just change the extension if you need to. Let me know if you have any questions about how to adapt the script to your needs? >>> >>> http://plato.cgl.ucsf.edu/pipermail/chimera-users/2009-March/003610.html <http://plato.cgl.ucsf.edu/pipermail/chimera-users/2009-March/003610.html> >>> >>> ?Eric >>> >>> Eric Pettersen >>> UCSF Computer Graphics Lab >>> >>>> On Aug 31, 2016, at 3:02 AM, Healey, Joe <J.R.J.Healey(a)warwick.ac.uk <mailto:J.R.J.Healey@warwick.ac.uk>> wrote: >>>> >>>> Hi, >>>> >>>> I have several hundred structure simulations that were performed on our HPC resource. Each simulation has approximately 5 models associated with it. I also have HMM model homologies for each of these proteins. >>>> >>>> I would like to use chimera's matchmaker function to compare the simulated structured to their nearest 'real' counterpart with a resolved structure. In order to do this I plan to script chimera to pull in the simulated models as well as the near homolog and then run the matchmaker algorithm. >>>> >>>> My question is: can the RMSD value that is returned in the reply log be 'accessed directly' such that for each match that is performed, I can obtain the RMSD and thus a ranking for the best fitting structures - or will it require parsing the reply log? >>>> >>>> Thanks >>>> >>>> Joe Healey >>>> >>>> >>>> M.Sc. B.Sc. (Hons) >>>> PhD Student >>>> MOAC CDT, Senate House >>>> University of Warwick >>>> Coventry >>>> CV47AL >>>> Mob: +44 (0) 7536 042620 | Email: J.R.J.Healey(a)warwick.ac.uk <mailto:J.R.J.Healey@warwick.ac.uk> >>>> >>>> Jointly working in: >>>> Waterfield Lab <http://www2.warwick.ac.uk/fac/med/research/tsm/microinfect/staff/waterfield…> (WMS Microbiology and Infection Unit) >>>> and the Gibson Lab <http://www2.warwick.ac.uk/fac/sci/chemistry/research/gibson/gibsongroup/> (Warwick Chemistry) >>>> >>>> Twitter: @JRJHealey <https://twitter.com/JRJHealey> | Website: MOAC Page <http://www2.warwick.ac.uk/fac/sci/moac/people/students/2013/joseph_healey>_______________________________________________ >>>> Chimera-users mailing list: Chimera-users(a)cgl.ucsf.edu <mailto:Chimera-users@cgl.ucsf.edu> >>>> Manage subscription: http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users <http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users> >>> >>> >>> _______________________________________________ >>> Chimera-users mailing list: Chimera-users(a)cgl.ucsf.edu <mailto:Chimera-users@cgl.ucsf.edu> >>> Manage subscription: http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users <http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users> >> > >

2 6