- ChimeraX-users - RBVI Mailing Lists

Side by side view in ChimeraX (like pymol grid_mode)
by Thibault TUBIANA 06 May '22

06 May '22

Hello all :-) In ChimeraX, is there a way to display multiple models into the same windows like the "grid_mode" mode in Pymol ? Thank you! Best regards, Thibault Tubiana.

2 2

Chimerax scripting error and passing user supplied arguments
by Ryan Woltz 06 May '22

06 May '22

Dear community, I'm sorry if this is a repeated question, I've seen a few posts about taking in a file for the command line and doing some things then exiting. I got most of this down and I'm sure at this point it is just a matter or correcting a few lines of code. I'm trying to automate a protein modeling protocol and I have 1000 decoys at this point. Now I need to insert the 5 metals into 5 different pockets into all my decoys for a final relax before selecting my top model. I have a template pdb with all the metals in place. For 4 of these metals I have used rosetta align mover to align a decoy onto the template around a single metal binding site, save the pdb, then add on the metal coordinates then repeat for the next binding site. So this method works because the chain is continuous around these metal binding sites. However the final metal binding site is in the center and coordinated by 4 chains coming together and the positioning needs to be very precise. Chimera or chimerax solves this very well but I don't know how to put this into a script. I have a few solutions if someone can help me edit this to work. My script (chimera-data-processing.cxc) looks like this right now. cd results/ open S_0001_X-4-all-chain-A.pdb ###X = 1-1000 open hsk2-cam-homology-from-map-refined-6cno-alignment-template-1-all-chain-A.pdb info align #1:240-244,765-769,1290-1294,1815-1819@o to_atoms #2:240-244,765-769,1290-1294,1815-1819@o matchAtomNames true; wait 300 save S_0001_X-4-all-chain-A-sf-aligned.pdb #1 format pdb relModel #2; wait 30 cd ../ close quit This script is perfect if run for have a single model with the exact name written into the script after typing: chimerax-daily --nogui chimera-data-processing.cxc I have two main solutions A) run a while loop that counts down and passes a filename to chimerax: while ( X < 1000 ) chimerax-daily --nogui chimera-data-processing.cxc S_0001_X-4-all-chain-A.pdb hsk2-cam-homology-from-map-refined-6cno-alignment-template-1-all-chain-A.pdb Where in the chimerax script I take the arg1 (S_0001*.pdb) and arg2 (hsk2*.pdb)(assuming chimera-data-processing.cxc is arg0) and modify to look like this: open $arg1 open $arg2 ... save $arg1-sf-aligned.pdb #1 format pdb relModel #2; wait 30 Option B is to use the foreach command described in the open command documentation ( https://www.cgl.ucsf.edu/chimerax/docs/user/commands/open.html) However, a copy and past of that command: *open myscript.cxc foreach ~/data/*.cif * *or* *open chimera-data-processing.cxc foreach results/S_0001_*-4-* all-chain-A.pdb results in this error: Missing filename suffix foreach I also tried this: chimerax-daily --nogui chimera-data-processing.cxc foreach results/S_0001_*-4-all-chain-A.pdb But it just runs the command. this seems to ignore everything after the first filename. I also read that if I used foreach inside the script I could use $file to replace the filename it is currently working on, would there be a way to fave $file1 and $file2? to automate this I'd also like the change the name of the template file I'm aligning to. I also attempted the python method using sys.arg = arg[1] then running things but I couldn't get it to work because I have no experience in python. What is the suggested way to complete this task? Again this is in a much larger script with user inputs so ideally it would look like this in a script (very shorthanded just for the concepts). I can also run a python command script in my bash sctipt pdb transformations if (y = 5) ###alignment metal number) while ( x < 1000) ###arg0 = command script arg1 = pdb to align arg2 = template pdb with metals chimerax-daily --nogui chimera-command-script.cxc $filenamebase$x-$alignmentnumber-*.pdb $template.pdb ###strings for variables are provided when script is initialized. ###code to extract metal coords for S_0001_$x-4-*.pdb done fi chimera-command-script.cxc: open arg1 open arg2 align model #1 to model #2 on specific atoms save arg1-new-filename.pdb close exit if someone would like a follow-up email with my attempt at chimera python sctipt using runcommand I'm very happy to supply that as a foundation to help fix, I just thought this one is easier to fix since it is really close, the script works but only for 1 pdb at a time. Apologies for bad scripting terminology and confusion with python, I've had to teach coding to myself so I might have misused terms. Thanks for any help that can be given, Ryan

4 6

volume threshold-level
by Dmitry Semchonok 05 May '22

05 May '22

Dear colleagues, Is there a tool in chimeraX that can set the correct volume threshold level when the Molecular weight of the molecule is known? Thank you! Kind regards, Dmitry

3 2

transparency vs. opacity
by Pranav Shah 04 May '22

04 May '22

Hi Elaine, I am trying to control the opacity/transparency of maps..however it seems to me that these two different concepts within the framework although they are used to achieve the same effect... My understanding is that transparency controls the surface being displayed whereas opacity controls the opacity of the color assigned to the surface. Is there room to unify these two concepts in ChimeraX? On an operational level, how do I control the opacity of a color in chimeraX commandline? Currently I am assigning colors based on hexcodes for specific shades... but there isnt a way to take into account the opacity of that shade with hexcodes only... Best, Pranav -- Pranav Shah Postdoctoral Research Fellow. Division of Structural Biology, Wellcome Trust Centre for Human Genetics, University of Oxford, Roosevelt Drive, Oxford OX3 7BN, UK

2 1

fitmap
by Dmitry A. Semchonok 04 May '22

04 May '22

2 1

Color by RMSD for nucleic acid
by Rundlet, Emily 03 May '22

03 May '22

Hi ChimeraX developers, Is there a way to color nucleic acids by seq_rmsd? When I use the following, no atoms are colored. mm #1 to #2 matrix Nucleic showAlignment true color byattribute seq_rmsd #2 Thank you in advance for your help! -Emily ________________________________ Email Disclaimer: www.stjude.org/emaildisclaimer Consultation Disclaimer: www.stjude.org/consultationdisclaimer

3 3

security bug in ChimeraX 1.1 and earlier on WIndows
by Greg Couch 03 May '22

03 May '22

FYI, ChimeraX uses the Qt user interface toolkit, and recently a security bug was found in it on Windows. It affects ChimeraX 1.1 and earlier. See https://kb.cert.org/vuls/id/411271 for details. Our recommended fix is to update to the current version of ChimeraX. -- Greg

1 0

Help pages tool hyperlink from Tool index
by William Hofsøy 03 May '22

03 May '22

Hi, I am currently writing some help pages for my tools and noticed that on the Tool index page in the chimerax help pages my tools get listed but not hyperlinked to the pages I write. Is there a way I can fix this? I also noticed that my commands are not listed here. --William

2 1

Dali server alignment
by HONG ZHAN 02 May '22

02 May '22

Dear all, I am trying to align molecules based on Dali server output. However, I can’t download the superposition results from the Dali server webpage (or anyone know?), but I can generate an alignment file with 2D structural alignment from the Dali server. Is there a way to use matchmaker to superpose molecules based on output alignment from Dali server? How to run structure conservation calculation to color the molecule? If I just give Dali’s alignment file, it still calculate conservation based on sequence not structure. Best, Hong

3 6

Re: [chimerax-users] ChimeraX question
by Tom Goddard 29 Apr '22

29 Apr '22

Hi Ryan, ChimeraX has a "morph" command. Here's an example using it open 2vaa open 2vab matchmaker #2 to #1 morph #1,2 save morph.pdb model #3 allCoordsets true Morphing of course only contains the atoms that are in both structures. So if you have an unbound and a bound state, and the bound tapasin is only in one structure then it will not be included in the morph since there is nothing to morph it to in the unbound state. You can delete chains that the two structures do not have in common before morphing, or you can use the same true option to the morph command "morph #1,2 same true" which will ignore chains which do not appear in both structures. Tom > On Apr 29, 2022, at 9:28 AM, Kissinger, Ryan (NIH/NIAID) [C] <ryan.kissinger(a)nih.gov> wrote: > > Hi Tom, > > I’m trying to create a morph animation in ChimeraX and I am having a lot of trouble finding decent tutorials online. I have pdb entries for unbound MHC1 and Tapasin and bound MHC1 and Tapasin. If possible I would like to morph between the bound and unbound forms. I would like to do this with ribbon and the atomic models. Any advice? > > Best Regards, > > Ryan

2 1