- ChimeraX-users - RBVI Mailing Lists

Suppress errors when using runscript?
by Vorländer,Matthias Kopano 03 Feb '25

03 Feb '25

Hi, I wondered if it is possible to suppress/ignore an error message when executing scripts with runscript. Specifically, I have a script to analyze alphafold (3) predictions and select residue interfaces residues using the interface select command. The selection for each of the five model is then saved in a named selection: # Select interface residues using the interfaces command. interfaces select last-opened&/A contacting last-opened&~/A bothSides true # Create a named selection for the interface residues. name frozen interface_residues$2_model_4 sel This works fine when all models show contacts , however, if one of the models/predictions does not have any residues involved in an interface, the script stops with the error “Nothing is selected by the specifier”. It would be great if the script does not pause at this message. Below the full .cxc script I am using: Best, Matthias # Generate AF contacts for all predictions of an AF3 run # # Usage: # runscript analyse_af3_prediction.cxc <AF_screen_directory> <custom_prediction_name> # The script assumes that five models were generated by AF. It expects the CIF files to be in a sub‐directory # matching the pattern "*model_X*.cif" and the corresponding JSON files in a sub‐directory matching "*data_X*.json" # for each model number X (4, 3, 2, 1, 0). # For the alignment and mseries slider command to work properly, the script assumes no other models are open. echo "Usage: runscript analyse_af3_prediction.cxc <AF_screen_directory> <custom_prediction_name>" ##### Model 4 ##### # Open the structure and its JSON scores for model 4. open $1/*model_4*.cif open $1/*data_4*.json structure last-opened hide show c # Generate AF contacts (using JSON scores) for chain A (or the first chain) alphafold contacts last-opened &/A name model_4_af_contacts_$2 # Create pseudobonds for the contacts and label them with residue names and numbers. sel last-opened&pbonds label sel pseudobonds text "{0.atoms[0].residue.name} {0.atoms[0].residue.number} to {0.atoms[1].residue.name} {0.atoms[1].residue.number}" # Save the pseudobond representation to a file. save $1/af_contacts_model_4_$2.txt format pseudobonds sel true echo "AlphaFold contacts (derived from JSON file) saved to $1/af_contacts_model_4_$2.txt" # Select interface residues using the interfaces command. interfaces select last-opened&/A contacting last-opened&~/A bothSides true # Create a named selection for the interface residues. name frozen interface_residues$2_model_4 sel show interface_residues$2_model_4 style interface_residues$2_model_4 ball # Save interface residue info to a file in the input directory. info residues interface_residues$2_model_4 saveFile $1/interface_residues$2_model_4_buriedArea_300.txt # Color the model for clarity. rainbow last-opened chains palette bupu color byhetero echo "Interface residues saved in named selection interface_residues$2_model_4 and info saved to $1/interface_residues$2_model_4_buriedArea_300.txt" ##### Model 3 ##### # Open the structure and its JSON scores for model 3. open $1/*model_3*.cif open $1/*data_3*.json structure last-opened hide show c # Generate AF contacts. alphafold contacts last-opened &/A # Create pseudobonds and label them. sel last-opened&pbonds label sel pseudobonds text "{0.atoms[0].residue.name} {0.atoms[0].residue.number} to {0.atoms[1].residue.name} {0.atoms[1].residue.number}" # Save the pseudobond representation. save $1/af_contacts_model_3_$2.txt format pseudobonds sel true echo "AlphaFold contacts (derived from JSON file) saved to $1/af_contacts_model_3_$2.txt" # Select interface residues using the interfaces command. interfaces select last-opened&/A contacting last-opened&~/A bothSides true # Create a named selection. name frozen interface_residues$2_model_3 sel show interface_residues$2_model_3 style interface_residues$2_model_3 ball # Save interface residue info to a file in the input directory. info residues interface_residues$2_model_3 saveFile $1/interface_residues$2_model_3_buriedArea_300.txt # Color the model. rainbow last-opened chains palette bupu color byhetero echo "Interface residues saved in named selection interface_residues$2_model_3 and info saved to $1/interface_residues$2_model_3_buriedArea_300.txt" ##### Model 2 ##### # Open the structure and its JSON scores for model 2. open $1/*model_2*.cif open $1/*data_2*.json structure last-opened hide show c # Generate AF contacts. alphafold contacts last-opened &/A # Create pseudobonds and label them. sel last-opened&pbonds label sel pseudobonds text "{0.atoms[0].residue.name} {0.atoms[0].residue.number} to {0.atoms[1].residue.name} {0.atoms[1].residue.number}" # Save the pseudobond representation. save $1/af_contacts_model_2_$2.txt format pseudobonds sel true echo "AlphaFold contacts (derived from JSON file) saved to $1/af_contacts_model_2_$2.txt" # Select interface residues using the interfaces command. interfaces select last-opened&/A contacting last-opened&~/A bothSides true # Create a named selection. name frozen interface_residues$2_model_2 sel show interface_residues$2_model_2 style interface_residues$2_model_2 ball # Save interface residue info to a file in the input directory. info residues interface_residues$2_model_2 saveFile $1/interface_residues$2_model_2_buriedArea_300.txt # Color the model. rainbow last-opened chains palette bupu color byhetero echo "Interface residues saved in named selection interface_residues$2_model_2 and info saved to $1/interface_residues$2_model_2_buriedArea_300.txt" ##### Model 1 ##### # Open the structure and its JSON scores for model 1. open $1/*model_1*.cif open $1/*data_1*.json structure last-opened hide show c # Generate AF contacts. alphafold contacts last-opened &/A # Create pseudobonds and label them. sel last-opened&pbonds label sel pseudobonds text "{0.atoms[0].residue.name} {0.atoms[0].residue.number} to {0.atoms[1].residue.name} {0.atoms[1].residue.number}" # Select interface residues using the interfaces command. interfaces select last-opened&/A contacting last-opened&~/A bothSides true # Create a named selection. name frozen interface_residues$2_model_1 sel show interface_residues$2_model_1 style interface_residues$2_model_1 ball # Save interface residue info to a file in the input directory. info residues interface_residues$2_model_1 saveFile $1/interface_residues$2_model_1_buriedArea_300.txt # Color the model. rainbow last-opened chains palette bupu color byhetero echo "Interface residues saved in named selection interface_residues$2_model_1 and info saved to $1/interface_residues$2_model_1_buriedArea_300.txt" ##### Model 0 ##### # Open the structure and its JSON scores for model 0. open $1/*model_0*.cif open $1/*data_0*.json structure last-opened hide show c # Generate AF contacts. alphafold contacts last-opened &/A # Create pseudobonds and label them. sel last-opened&pbonds label sel pseudobonds text "{0.atoms[0].residue.name} {0.atoms[0].residue.number} to {0.atoms[1].residue.name} {0.atoms[1].residue.number}" # Save the pseudobond representation. save $1/af_contacts_model_0_$2.txt format pseudobonds sel true echo "AlphaFold contacts (derived from JSON file) saved to $1/af_contacts_model_0_$2.txt" # Select interface residues using the interfaces command. interfaces select last-opened&/A contacting last-opened&~/A bothSides true # Create a named selection. name frozen interface_residues$2_model_0 sel show interface_residues$2_model_0 style interface_residues$2_model_0 ball # Save interface residue info to a file in the input directory. info residues interface_residues$2_model_0 saveFile $1/interface_residues$2_model_0_buriedArea_300.txt # Color the model. rainbow last-opened chains palette bupu color byhetero echo "Interface residues saved in named selection interface_residues$2_model_0 and info saved to $1/interface_residues$2_model_0_buriedArea_300.txt" ##### Final Steps ##### # Move all models (except the last-opened) to a common group (optional) mm #*&~last-opened&/A to last-opened&/A # Create an mseries slider for the five models. mseries slider #1-5 title $2 name $2_slider mseries #1-5 echo "Use the slider to go through the models!" echo "Pairwise contacts (based on the JSON fle) and all residues involved in an interface (based on a buried surface are of 300 Angstrom or more) have been saved in the input folder for each prediction" echo "You might want to increase stringency by calcualting the overlap between all predictions"

2 1

Getting two color keys in the same ChimeraX session
by Elaine Meng 03 Feb '25

03 Feb '25

> Begin forwarded message: > > From: Eric Pettersen via Chimera-users <chimera-users(a)cgl.ucsf.edu> > Subject: [Chimera-users] Re: key command > Date: February 3, 2025 at 11:52:12 AM PST > To: Fabian Glaser <fglaser(a)technion.ac.il> > Cc: Chimera BB <chimera-users(a)cgl.ucsf.edu> > Reply-To: chimera-users(a)cgl.ucsf.edu, Eric Pettersen <pett(a)cgl.ucsf.edu> > > Hi Fabian, > There is no direct support for multiple keys, but it is pretty easy to “trick” ChimeraX into showing multiple keys. You would create your first key normally, getting it exactly the way you want it to look before creating the next key. You create the next key with the following Python code: > > from chimerax.color_key.model import ColorKeyModel > m = ColorKeyModel(session) > > If the above was in a file ending in ‘.py’, you could just use the ‘open’ command to run it and create the key. I’ve attached a file with that code in it for convenience. Once you have created the key, then all operations that affects keys (like the ‘key’ command) will affect the new key and leave the original key as is. If you had the Color Key tool up, it will continue to show the old values. Just close it and re-open it to show the new values. > > --Eric > > Eric Pettersen > UCSF Computer Graphics Lab > > >> On Feb 3, 2025, at 6:31 AM, Fabian Glaser via Chimera-users <chimera-users(a)cgl.ucsf.edu> wrote: >> >> Hi, >> >> We would be interested to add more than one key color bar in the same chimerax session. But we noticed the key command allows only one color key can be present at a time. >> >> Is there a way to add more than one color key in the same session? >> >> Thanks a lot, >> >> Fabian >> >> >> Fabian Glaser, PhD >> >> Technion Centre for Structural Biology (TCSB), Computational Section, Head >> Technion Human Health Initiative (THHI) >> Technion - Israel Institute of Technology, Haifa, Israel >> fglaser(a)technion.ac.il >> +972 (0)48293701 >> >> >> >> _______________________________________________ >> Chimera-users mailing list -- chimera-users(a)cgl.ucsf.edu >> To unsubscribe send an email to chimera-users-leave(a)cgl.ucsf.edu >> Archives: https://mail.cgl.ucsf.edu/mailman/archives/list/chimera-users@cgl.ucsf.edu/ > > _______________________________________________ > Chimera-users mailing list -- chimera-users(a)cgl.ucsf.edu > To unsubscribe send an email to chimera-users-leave(a)cgl.ucsf.edu > Archives: https://mail.cgl.ucsf.edu/mailman/archives/list/chimera-users@cgl.ucsf.edu/

1 0

Grid mode and multiple viewports
by Roden Deng Luo 03 Feb '25

03 Feb '25

Hi, I find there is a `tile` command. After I enter the tile mode, the zoom is global. I guess `tile` is implemented by transforming the active models in the global scene explicitly. I wonder if there is a way to have the zoom independent for each model, similar to the grid mode in PyMol. A further question, more from a developer's perspective, is whether having multiple 3D rendering viewports in ChimeraX is possible. I imagine each viewport can have its own active models in the scene and its own camera. Different cameras and other operations might be synced upon users' request. This video (https://youtu.be/7KZ1P-nPav8?si=Z1hiBBuRiWop2WFC) demos the feature (in an old version of Houdini). Thanks, Roden -- This message and its contents, including attachments are intended solely for the original recipient. If you are not the intended recipient or have received this message in error, please notify me immediately and delete this message from your computer system. Any unauthorized use or distribution is prohibited. Please consider the environment before printing this email.

4 8

aligning maps generated by chimeraX
by msteves＠berkeley.edu 01 Feb '25

01 Feb '25

Hello, I'm new to Chimera, but would like to use it to generate surface maps of electrostatic potential and hydrophobicity, then output for further analysis. I've generated maps using coulombic and mlp commands, and successfully saved these to .MRC format. However, when I go and read these files, the dimensions of the maps are different, so I'm not sure how to properly overlay these maps in my later analysis. Is there a way to have Chimera generate maps only within a given spatial range, so I end up with maps of the same size? Or another way you would recommend going about this? In later analysis I'd also like to make use of the SAS calculated by Chimera, and I imagine I will also have difficulties properly aligning the surface to the maps, so any advice on this would also be greatly appreciated! Thanks, Megan Steves

2 1

Question about Surface Volume Overlap Calculation
by John Watters 31 Jan '25

31 Jan '25

Hello, I was wondering if anyone had any tips/tricks/useful functions for quantifying the volume of the overlap between two surfaces in ChimeraX? I am looking specifically to make a metric to compare steric clashes of molecules in different orientations. Thank you for any help you can provide! Best, John

3 4

Add help for reusable ChimeraX script that takes command line arguments?
by Vorländer,Matthias Kopano 30 Jan '25

30 Jan '25

Hi, It’s me again 😉 I wondered how I can add a help page to a re-usable ChimeraX script that is run via runscript. I saved the following script in my ChimeraX presets: #Required arguments: first is the uniprot ID, second is the model ID that will be used for mapping! alphafold fetch $1 #get the missense data open $1 fromDatabase alpha_missense format amiss #associate missense data to structure mutationscores structure $2 #label each residue with matric mutationscores label #1 amiss height 3 #define a new score called avg by avergaing over all substitution scores mutationscores define avg fromScore amiss setAttribute true combine mean #color by avg missense score color byattribute r:avg $2 target csab palette palette bluered range full #scale cartoon by acg missense score cartoon byattribute r:avg #!$2 and then defined an alias in my startup commands so I don’t need to type the full path alias map_missense runscript /Users/matthias.vorlaender/Documents/ChimeraX_presets/map_missense.cxc $1 $2 Is there a way to display the expected user inputs when the script is called without arguments? Cheers, Matthias

2 3

Feasibility of selecting chain ID via uniprot ID?
by Vorländer,Matthias Kopano 30 Jan '25

30 Jan '25

Hi again : ) When comparing structures from different sources (such as different PDB depositions, or an Alphafold structure), it is always tedious to figure out which chain belongs to which protein (especially for large structures like spliceosomes) for consistent coloring or alignment. I wondered if it would be technically feasible to use the uniport ID to match a specific chain directly, since the association between chain ID and uniport ID should be present at least in structures from curated databases such as the PDB and Alphafold. So the ambition would be to be able to use ‘color IF4A3_HUMAN red’ directly, without having to look chain IDs. An added bonus would be to be able to skip the species identifier, in case you have models containing IF4A3_MOUSE and IF4A3_HUMAN, and match both with IF4A3 …. Would be curious to hear your thoughts. Thanks a lot, and best wishes, Matthias

3 3

Encode: saving a movie on a Windows laptop, and colouring residues in the multiple sequence
by Anne-Frances Miller 30 Jan '25

30 Jan '25

Dear ChimeraX My class is having a blast with this! I have assigned the class to make a movie, based on the morph and encode commands. On Macs, your instructions are flawless. However I have students using Windows10 and Windows 11. For these, could you please provide an example movie encode command including path and good choice of format please? Also, when I open a multiple sequence alignment in ChimeraX, is there a way to apply Taylor or other colouring to the MSA? This used to be possible in Chimera's MAV. Thank you, Anne-Frances Novel Meritorious Research (and teaching) —•—•—•—•—•—•—•—•—•—•—•—•—•— A.-F. Miller, Ph. D. (she/elle/sie) Distinguished Prof. of Arts and Sciences Professor of Chemistry University of Kentucky https://uky.zoom.us/j/3250841836 https://millerlablive.wordpress.com/ Visiting Fellow, Einstein Foundation of Berlin

3 4

Possible to plot attribute value over sequence?
by Vorländer,Matthias Kopano 30 Jan '25

30 Jan '25

Dear ChimeraX team, I wondered if it is possible to generate a line or scatter plot of an attribute value (B factor, Alphafold pLDDT scores, or AlphafoldMissense score…) over the residue number? For example something like this (taken from https://journals.plos.org/plosone/article?id=10.1371/journal.pone.0089116) : [cid:image001.png@01DB7242.B1C3F8A0] Specifically, for the AlphafoldMissense data, such a plot would be more informative then the default histogram that is displayed when fetching the scores when looking for functional hotspots in your POI. Many thanks, MAtthias

3 5

Chai discovery data
by Reza Salavati, Prof. 29 Jan '25

29 Jan '25

Hello, I have downloaded the Chai discovery data in .npy and .json file formats. How can I open these files in Chimera X? Thank you Reza

4 4