- ChimeraX-users - RBVI Mailing Lists

Re: [chimerax-users] exporting crosslinks histogram
by Tom Goddard 17 May '22

17 May '22

Hi Pim, I added a Save Plot As... menu entry to the crosslinks histogram, can save as svg or pdf or png. I also added bins, maxLength, minLength and height options to the command to control how many bins, the x scale and height of y scale (in counts per bin). These changes will be in tomorrow's ChimeraX daily build which will become ChimeraX 1.5 eventually. Too late to make it into ChimeraX 1.4. Tom > Dear ChimeraX team, > > I am analysing multiple datasets generated using crosslinking-ms, loaded as separate pseudobond files on my model. The crosslinks histogram plots works just fine and I’d now like to put some plots neatly side-by-side. My questions: > > - can the histogram be saved or exported? Ideally as an svg or pdf file. But png or jpg would also already be useful. > I could export the .pb files and plot them elsewhere, but that seems like a workaround. > > - is it possible to specify the x and the y axis of the generated histograms, as well as the bin size? This would be helpful to show them side-by-side. > > Many thanks in advance, > > Pim > -running v1.2.5 on mac os catalina 10.15-

1 0

Re: [chimerax-users] [Chimera-users] Set scale factors from displayed contours?
by Tom Goddard 17 May '22

17 May '22

Hi Oliver, In tomorrow's daily ChimeraX builds you can use volume maximum #1,2 scale fromlevels and it will scale each map by 1/level where level is the displayed surface level before taking the maximum. This allows merging maps that have different intensity levels. Thanks for the suggestion. This will be in ChimeraX 1.5 but not in 1.4 and is not in Chimera. Tom > On May 6, 2022, at 2:38 PM, Elaine Meng via Chimera-users <chimera-users(a)cgl.ucsf.edu> wrote: > > Hi Oli, > Oh, I see, you already figured out the scale factors (essentially) by adjusting thresholds interactively. So the automation would simply be doing the math of taking the ratio of the thresholds and scaling the map(s) according to that, which sounds more feasible. > > I can make a feature request ticket (probably ChimeraX), but can't make promises as I would not be the one to triage it. > Best, > Elaine > ----- > Elaine C. Meng, Ph.D. > UCSF Chimera(X) team > Department of Pharmaceutical Chemistry > University of California, San Francisco > > >> On May 6, 2022, at 1:42 PM, Oliver Clarke via Chimera-users <chimera-users(a)cgl.ucsf.edu> wrote: >> >> Hi Elaine, >> >> I don't think I explained very well - I mean to say, calculate the scale factors from the displayed thresholds. >> >> E.g. if the user has set map1 to a threshold of 1.0, and map2 to a threshold of 2.0, in order to make them visually equivalent in a combined map, map1 needs to be scaled by a factor of 2. >> >> It would be great if this could be done automatically as a flag for the volume or vop command (as there is already a flag to provide per map scale factors). >> >> Cheers >> Oli >> >>> On May 6, 2022, at 4:39 PM, Elaine Meng <meng(a)cgl.ucsf.edu> wrote: >>> >>> Hi Oli, >>> Our volume expert Tom is on vacation until mid-month and might not be reviewing questions for a while... so I will give my not-so-expert thoughts: >>> >>> To me, the initial identification of "overlapping bits" sounds quite subjective and thus difficult to encode. If the user could designate the boundaries of this area (by drawing a box?), perhaps automatic scaling of one map to make its values as similar as possible to the values of another map inside that box could be done. However, drawing the box would probably be just as time-consuming as the process that you are already using. >>> >>> Also, new developments/agorithms would be going only into ChimeraX, not Chimera. ChimeraX has mostly the same "vop" options as Chimera, except included in the "volume" command instead of "vop": >>> <https://rbvi.ucsf.edu/chimerax/docs/user/commands/volume.html#vop> >>> >>> Best, >>> Elaine >>> ----- >>> Elaine C. Meng, Ph.D. >>> UCSF Chimera(X) team >>> Department of Pharmaceutical Chemistry >>> University of California, San Francisco >>> >>> >>>> On May 6, 2022, at 1:23 PM, Oliver Clarke via Chimera-users <chimera-users(a)cgl.ucsf.edu> wrote: >>>> >>>> Hi, >>>> >>>> When combining maps using vop maximum (or volume maximum in chimerax), we often want to adjust the scale factors of the maps such that the overlapping bits of two aligned maps are visually on the same scale. It is possible to calculate the appropriate scale factors from the displayed thresholds, but this becomes a little time consuming when combining many maps. >>>> >>>> Would it be possible to add a flag, something like "get scalefactor from display thresholds" that would take care of this automatically? >>>> >>>> Cheers >>>> Oli >>> >> >> >> _______________________________________________ >> Chimera-users mailing list: Chimera-users(a)cgl.ucsf.edu >> Manage subscription: https://www.rbvi.ucsf.edu/mailman/listinfo/chimera-users > > > _______________________________________________ > Chimera-users mailing list: Chimera-users(a)cgl.ucsf.edu > Manage subscription: https://www.rbvi.ucsf.edu/mailman/listinfo/chimera-users >

1 0

Removing text labels generated with command distance
by Shubham Devesh Ramgoolam 16 May '22

16 May '22

Hi, I hope you are doing well. I have a python script, to be executed in nogui mode, which does the following: -open a pdb file -define centroids for specific residues -draw pseudobonds between specific residues -> I am using the command distance for this purpose The issue arises when I have to save the session (.cxs) when in nogui mode, I am getting this error: QGuiApplication::font(): no QGuiApplication instance and no application font set. I believe it is because of the text labels created by the command distance. So, I included the following in my code: run(session, ‘hide #2.*’) and I can now save the session. I noticed in the model panel that the pdb file is assigned to #1 and the pseudobond to #2. The text labels are created as a submodel of the pseudobond (so #2.1, 2.2,..), hence why the ‘hide #2.*’ Will this fix always work? Is the model assignment also dependent on the pdb file being opened? Regards, Shubham Sent from Mail<https://go.microsoft.com/fwlink/?LinkId=550986> for Windows

4 4

Where is the renumber sequence tool in ChimeraX?
by Phan, Isabelle 16 May '22

16 May '22

The sequence of my structure includes a tag and I want to renumber from -7 so that the first residue is the start Met. I tried the old chimera command resrenumber = FAIL. I googled and searched the docs, read the list of commands one by one, and still could not figure it out. The API has a renumber_sequence() function but I don't understand the description (what’s "the given residues”? The chain?) and was unable to find an example on how to call this from the interface. I opened the python shell and typed: from chimerax.core.commands import run run(renumber_residues('#1/A:*', -7)) # FAIL Sorry if this is trivial, I am still completely baffled by ChimeraX despite using it for years. Thanks, Isabelle CONFIDENTIALITY NOTICE: This e-mail, including any attachments, is for the sole use of the intended recipient(s) and may contain confidential and privileged information protected by law. Any unauthorized review, use, disclosure or distribution is prohibited. If you are not the intended recipient, please contact the sender by reply e-mail and destroy all copies of the original message.

2 5

align 400 copies of the same density map into a circle
by SHIQING LIAO 16 May '22

16 May '22

Hi there, I am a new user and would really appreciate any input on the task I want to do using Chimerax. Currently I have a mrc map, and it is a part of a cylindrical tube-shaped structure (like the area I drew in the picture). The map looks like a curved plane from the side of the cylinder, if it helps. If I piece 400 identical copies of map with some rotation and transition, I would ideally get a cynlider shaped strutcure, without top and bottom (like the picture below). The reason why I want to piece the mrc together is that I don't have the full look of the tube, but I know the rough radius (that's how 400 was derived), and I want to make a model for the full tube, using the model of the composing protein. I managed to manually fit the models of individual protein into one map. For the other mrc copies, I tried to do the same again but it took a lot of time and not very accurate, so I'm looking for any coding solutions that might help. I'm currently thinking if it is possible for me to place a marker at the center of one mrc map. Then make 400 copies of the map with the marker. Then write a code to tell each 1/400 mrc map to rotate 0.9 degrees around the marked center and move to a coordinate. Then I can repeat the same for the protein models I mapped to one mrc map. But I have no idea how to achieve it, and any suggestions would be greatly appreciated! Thank you so much for helping me!!! Best, Shiqing [image: cf23fde5-632f-4568-9700-9307f24e0f6c.png]

2 1

crossfade issue
by Omid Haji-Ghassemi 16 May '22

16 May '22

Dear ChimeraX community, I am having some issues with the crossfade command while making a movie. For some reason the command stretches out the image while running the transition. I am just doing a simple transition for a movie as part of a more complicated movie. To illustrate the problem I made a simple script that reproduces the issue. movie record supersample 3 size 640,480 wait 60 crossfade; cartoon hide #1/D-F; wait 60 movie encode framerate 60 output test.mov Be great if someone can let me know what I am doing wrong. Best, Omid --------------------------------------------------- --------------------------------------------------- Dr. Omid Haji-Ghassemi Assistant Professor, Biological Sciences Department University of Calgary https://hajighassemilab.org

3 4

Re: [chimerax-users] [Chimera-users] Making a movie from map series while keeping map's colors
by Elaine Meng 16 May '22

16 May '22

Hi Michal, Since this is a ChimeraX question I'm putting it on chimerax-users(a)cgl.ucsf.edu -- otherwise the Chimera users will be confused! In the future you may want to use that address for ChimeraX questions. If I remember correctly, "vseries" re-renders the surfaces at the time of playback, unless you specify storing more of them with the "cacheFrames" option -- so try using "vseries" with "cacheFrames 20": <https://rbvi.ucsf.edu/chimerax/docs/user/commands/vseries.html> See previous answer to a similar question: <https://www.rbvi.ucsf.edu/pipermail/chimerax-users/2022-February/003271.html> I hope this helps, Elaine ----- Elaine C. Meng, Ph.D. UCSF Chimera(X) team Department of Pharmaceutical Chemistry University of California, San Francisco > On May 16, 2022, at 12:45 AM, Michal Arie Shmuelov via Chimera-users <chimera-users(a)cgl.ucsf.edu> wrote: > > Hi everyone, > I have 20 maps (from 3D variability analysis) that I open as series in chimeraX. > I’m coloring the map series using the command “color zone #X near #X distance X” (coloring the map according to pdb model). > Next I’m trying to make a movie of this map series while retaining the map colors. > For some reason when the morph/ vseries starts to run the colors disappear. > Can someone help me to solve this problem? > What are the commands that I should write? > Thanks, > Michal >

2 1

Python code for obtaining chain information
by Shubham Devesh Ramgoolam 16 May '22

16 May '22

Hi, I hope you are doing well. I have a python script in which I need generate a list of all the chains present in a structure: e.g [A,B] However, I do not know which module to import in order to write a code for this. Regards, Shubham Sent from Mail<https://go.microsoft.com/fwlink/?LinkId=550986> for Windows

2 1

Showing Bond/PBond to Cartoon/Ribbons
by Anderson, Jacob 13 May '22

13 May '22

There has been a prior post on this problem, but I wanted to ask about some asymmetry I see in ChimeraX when displaying residues emanating from ribbons/cartoons. https://www.rbvi.ucsf.edu/pipermail/chimerax-users/2021-October/002850.html For example, sometimes when I select a residue on a ribbon/cartoon and write "show sel atoms,bonds", I get a floating residue. However, other times I get a nice connection between the beta carbon and the ribbon (the desired behavior). How can users get a consistent behavior where there is a connection from the beta carbon to the ribbon? Good Example - connected to ribbon [cid:6ddcc2ee-f9b2-4060-8e3f-fde8ccc74881] Two floating examples: [cid:46516984-036d-4a39-9713-9056a161548a]

3 6

no gui mode and redirecting output of some commands
by Shubham Devesh Ramgoolam 12 May '22

12 May '22

Hi, I hope you are doing well. I have a python script (myscript.py) as follows: from chimerax.core.commands import run import os os.chdir("C:\\Users\\Shubham\\Desktop") run(session, "open C:\\Users\\Shubham\\Desktop\\PDB_Script\\1ebt.pdb") # location of my pdb file run(session, "sel aromatic & protein") run(session, "log clear") run (session, "info residue sel") run(session, "log save aromatic_list.html") I am using the following prompt for calling chimerax from command prompt: Gui enabled: chimerax --script C:\\Users\\Shubham\\Desktop\\myscript.py Gui disabled: chimerax --nogui --script C:\\Users\\Shubham\\Desktop\\myscript.py With gui mode enabled, the html file, aromatic_list.html, is created. However, when disabling gui, nothing is created. I suspect this is because I am relying on the log which may not be functional when gui is disabled? Is there a way of redirecting the output of certain commands (e.g info, distance, angle) to a desired text file when disabling gui? I tried with sys.stdout but it did not work (empty text files were created.). Regards, Shubham Sent from Mail<https://go.microsoft.com/fwlink/?LinkId=550986> for Windows

2 6