June 2024 - ChimeraX-users - RBVI Mailing Lists

Re: CHimeraX and Chimera website down
by Elaine Meng 19 Jun '24

19 Jun '24

Hi Yaikhomba, Yes, we had a planned outage for server maintenance the last few days. Usually our email is out at the same time so it will not help to send email either, sorry! There was an announcement about it to this chimerax-users email list and also a news item on the ChimeraX home page, although of course you can't see the home page when the website is down. You can change the map origin with the ChimeraX "volume" command option "origin" or "originIndex", see the help for what these mean and how to use them: <https://rbvi.ucsf.edu/chimerax/docs/user/commands/volume.html#dimensions> I hope this helps, Elaine ----- Elaine C. Meng, Ph.D. UCSF Chimera(X) team Resource for Biocomputing, Visualization, and Informatics Department of Pharmaceutical Chemistry University of California, San Francisco > On Jun 19, 2024, at 8:49 AM, M Y <yai25ram(a)gmail.com> wrote: > > Hi Elaine > > I hope you are well. > > I had trimmed the box size of a cryoEM-map using the RELION software (to 375 pixels from 600 pixels). To fit in a PDB model to the trimmed map, I also rescaled the trimmed map to 1.048 A/pixel from 1.07 A/pixel. However, these 2 maps do not overlay with one another (image attached). > > I wanted to change the origin of the trimmed map such that 2 maps overlay with one another. How could I do this in ChimeraX? > > I tried looking at the Chimera or the ChimeraX website. But, the website is either down or I am not able to access it (Screenshot attached). > > Your suggestions will be helpful. > > Many thanks > Yaikhomba > <image.png> > <image.png>

1 0

Creating bond between two selected atoms
by Eleanor Streeper 17 Jun '24

17 Jun '24

Hi Chimera X, I hope you are well. I am a student at Swarthmore College and am currently trying to use Chimera X to covalently attach glycans to asparagine residues. I was able to get 3 out of 5 of the glycans attached by selecting both atoms > structure editing > build structure > adjust bonds > add reasonable bonds. I did the same process for the other 2 glycans and no bond was made. The log says bond sel, 0 bonds. The glycans are combined with the protein using the combine command. I believed it may be an issue with an unreasonable binding distance but when I moved the glycan closer, the bond still didn't work. What else should I try? Also, do you know why it is working for some glycans but not others? Best, Ellie Streeper

2 4

Updating pdb SEQRES record with a loaded .fasta file
by Samuel Haysom 14 Jun '24

14 Jun '24

Hello, I would like to use a .fasta file that I have loaded into ChimeraX to update the SEQRES records of a PDB file that is also loaded into ChimeraX. Currently the loaded PDB file only has SEQRES records for the residues that are modelled, while the .fasta file has the full, sequence with no gaps. I have tried loading in the .fasta file and then using the $ sequence associate command to associate it with the relevant pdb chain. However, when I then save out the pdb file (using the $ save command) and then examine the pdb file it still only records the modelled residues. Is there any way to force the written pdb file to have the full sequence as uploaded in the .fasta file? I am running chimeraX version 1.7 (2023-12-19). Thanks for your help, Sam Samuel Haysom Research Scientist • Biophysics Nxera Pharma UK Limited +44 (0)1223 949201 • Steinmetz Building, Granta Park, Great Abington Cambridge, CB21 6DG, United Kingdom The information in this e-mail is confidential and may be legally privileged. It is intended for the exclusive attention of the addressee stated above and should not be copied or disclosed to any other. If you have received this transmission in error, please contact the sender immediately, make no use of its contents and delete this message from your system. Contact and registered office address: Nxera Pharma UK Limited, Granta Park, Great Abington, Cambridge, CB21 6DG. Company No: 06267989.

2 3

RBVI downtime
by Eric Pettersen 13 Jun '24

13 Jun '24

Hi all, The RBVI compute cluster will be down for scheduled maintenance Monday and Tuesday of next week (June 17/18) which means that the Chimera and ChimeraX web sites (including the ChimeraX toolshed) and mailing lists will be down during that time. This also means that web services from Chimera/ChimeraX that access our compute cluster will not work (Blast Protein, Modeller tools, sequence realignment). We are sorry for any inconvenience. --Eric Eric Pettersen UCSF Computer Graphics Lab

1 0

Resizing an ion model representation
by Ioannis Skalidis 13 Jun '24

13 Jun '24

Dear all, I want to show a metal ion surrounded by some side-chains for a figure panel, but the ion ball representation is quite small and gets lost in the figure panel. Is there a way to select only the specific ion and increase its size? I tried a bit with the "size" command, including the ballScale option but I had no luck. Any suggestions? Thanks a lot, Ioannis. Dr. Ioannis Skalidis | Structural Biochemist | Twitter: @imjohnskal <https://twitter.com/imjohnskal>

2 1

Feature Request: Visual feedback on selecting menu items
by Cam McMenamie 12 Jun '24

12 Jun '24

Hello! I'm just suggesting a UI feature that I think would be very helpful. Currently, there are several menu items that do not show any feedback if they are "selected" or not, for example, in the screenshot below (hard to tell which options have been selected / toggled). This is especially the case for mutually exclusive selections such as in the 'coloring' section in the "molecule display" tab. [cid:34830b00-bed7-4734-b191-dc4fd0fad018] Thanks in advance, Cam

2 1

Filtering PAE Contacts–Interfaces
by David S. Fay 10 Jun '24

10 Jun '24

Hi, Was attempting (unsuccessfully) to filter out lower scoring (red/yellow) PAE contacts when viewing interaction interfaces between two or more proteins. I normally use a command such as: alphafold contacts /A to /B distance 5 But then I’d like to limit what I see to regions with high pLDDT scores (blue lines). I tried playing with range, but keep getting errors. Thanks! David

3 3

Animate RELION jobs
by James Hooker 10 Jun '24

10 Jun '24

Hi all, Is there an easy way to generate animations of RELION classification jobs in ChimeraX? I came across https://github.com/Stefan-Zukin/relionReport, but haven’t been able to get it to work. Best, James James Hooker jhooker(a)mrc-lmb.cam.ac.uk PhD candidate Löwe group, Structural Studies division MRC Lab of Molecular Biology Cambridge, United Kingdom

2 1

Re: [Chimera-users] ChimeraX_InfoRequest
by Elaine Meng 10 Jun '24

10 Jun '24

Hello Dr. Altan Kara, Resolution is just the number of pixels width and height. To get higher resolution, specify larger numbers for the pixel width and height in the save image options (GUI or save command). Chimera has a separate setting for physical size (e.g. inches or cm) and then the pixels per physical size unit (dpi), but the result is still the total number of pixels. For example 2 inches wide with 600 pixels/inch is 1200 pixels wide. ChimeraX does not have separate settings for the physical size. For the same example, in ChimeraX you would just choose 1200 pixels wide. The result will not be marked with a physical size, however. You would have to use some other program like Photoshop if you want to mark it as 2 inches wide. If you know the physical size you want and the resolution at that size, just calculate the total pixel width and height and enter that when you save the image. For the future, if you want to ask a ChimeraX question it is better to use the ChimeraX email address: chimerax-users(a)cgl.ucsf.edu I hope this helps, Elaine ----- Elaine C. Meng, Ph.D. UCSF Chimera(X) team Resource for Biocomputing, Visualization, and Informatics Department of Pharmaceutical Chemistry University of California, San Francisco > On Jun 7, 2024, at 5:17 AM, Altan KARA (MAM-GMBE) via Chimera-users <chimera-users(a)cgl.ucsf.edu> wrote: > > To whom it may concern, > > How we can create high resolution images with ChimeraX? > No matter what I tried still the outputted images are not in high resolution.

2 2

ChimeraX segfaults when left on loop to make movies/images
by Artie Kushner 07 Jun '24

07 Jun '24

Hi, I hope the summer is starting off really well for you. tldr: I leave ChimeraX running movie/image generating commands from a python script in a tmux session and it segfaults on multiple machines. I was wondering if someone has encountered this before or can lend advice as to why this happens. I'd like to produce some movies (or at least pictures) in bulk (around 2000 individual molecules). I preprocess each one (load locally or open, add some representation, and either start recording a movie or "save XXXX.png"). After about 8-10 successful movies or images are produced I get a ChimeraX segfault and the process dies. This is quite perplexing to me because - I tried it on multiple machines: Ubuntu 22.04, installed from the .deb <https://www.cgl.ucsf.edu/chimerax/download.html> and Arch, installed from aur distribution. - both with `--offscreen` and `--nogui` flags and without them ( in fact, following up some threads in the mailing list it seems like it's not possible to produce movies with --nogui, just pngs) - there seems to be *no* excessive memory or page cache issues, although i haven't profiled it. Only one core is occupied when rendering both mp4/png's. Below I'm attaching a crash log from one of the machines (Ubuntu 22.04) and the description of my setup, which is a little awkward, but I don't think it's the cause since it produces some samples successfully before crashing. I hope i'm wrong. Best regards and any tips are really appreciated! I was really looking forward to having some crisp images. --------------------------------------------------------------------------------------- I start the process by calling a python script at the command line: ` chimerax --script "/home/student/dev/riboxyz/chimerax/ribrepr.py" --nogui --offscreen` when i tried generated images. This was after I gave up on getting the movies, which i did by just opening a ChimeraX gui on our workstation and calling a "loop.py" after a load the commands from "ribrepr.py" file from the gui command line. I do both in a tmux session. ```loop.py import os from chimerax.core.commands import run from concurrent.futures import ThreadPoolExecutor *DEST_DIR* = "/home/rtviii/dev/riboxyz/chimerax/movies" # This expands to ~1800 PDB ids: rcsb_id:= {4UG0, 7K00 and so on} for rcsb_id in os.listdir('/home/rtviii/dev/RIBETL_DATA'): movie_dir = os.path.join(DEST_DIR, "{}.mp4".format(rcsb_id)) print("Attempting to write ", movie_dir) if "{}.mp4".format(rcsb_id) in os.listdir(DEST_DIR): print("Skipping {}. Exists".format(rcsb_id)) continue else: run(session,"ribmovie {}".format(rcsb_id)) # <---- ribmovie command corresponds to the `*produce_and_save_movie*` in ribrepr.py below ``` The contents of the `ribrepr.py` are just to register a command for custom visualization of a given molecule + loop over molecules for image generation. ``` ribrepr.py import enum import json import os import sys from chimerax.core.commands import register, CmdDesc from chimerax.atomic import ResiduesArg, Chains, ChainArg, StructureArg from chimerax.atomic import Structure, AtomicStructure, Chain from chimerax.core.commands import run, runscript from chimerax.core.commands import CmdDesc, register, StringArg # https://mail.cgl.ucsf.edu/mailman/archives/list/chimera-users@cgl.ucsf.edu/… # just a data folder, it exists. *RIBETL_DATA* = os.environ.get("RIBETL_DATA", "/home/rtviii/dev/RIBETL_DATA") class *PolymerClass*(str, enum.Enum): tRNA = "tRNA" bS1m = "bS1m" # ... another ~200 enum members, redacted for brevity CHIMERAX_COLORS = [ ["tan", "#d2b48c"], ["sienna", "#a0522d"], #... more colors redacted for brevity ] # helper function to map colors def *get_polymer_color*(polymer_class:str): polyix = list(map(lambda x: x.value, (PolymerClass))) class_index = polyix.index(polymer_class) % len(polyix) if polymer_class == None: return "gray" return CHIMERAX_COLORS[(class_index % len(CHIMERAX_COLORS))][1] def *ribosome_representation*(session, structure: AtomicStructure): for _ in list(PolymerClass): print(_.name) from chimerax.core.commands import run from chimerax.core.colors import hex_color from chimerax.atomic import Residue, Atom, Chain # i get the PDB ID from the StructureArg rcsb_id = str(structure.name).upper().split('.')[0] # <-- the structure gets opened with the basename ex "(5AFI.cif)" run(session, "set bgColor white") run(session, "sym #1 assembly 1") # take only one assembly if multiple are available run(session, "hide #2") with open(os.path.join(RIBETL_DATA, rcsb_id, "{}.json".format(rcsb_id)), "r") as f: profile = json.load(f) # load my annotations, just a json file #--------------------------------------------[ This whole bit just colors by chain annotation basically] polymers = {} polymer_chains = [ *profile["proteins"], *profile["rnas"], *profile["other_polymers"], ] [ polymers.update(x) for x in [{chain["auth_asym_id"]: chain} for chain in polymer_chains] ] c: Chain for c in structure.chains: aaid = c.chain_id polyclass = None if len(polymers[aaid]["nomenclature"]) < 1: continue else: polyclass = polymers[aaid]["nomenclature"][0] if polymers[aaid]["entity_poly_polymer_type"] == "RNA": run(session, "surf /{}".format(aaid)) run(session, "transp /{} 100".format(aaid)) run(session, "color /{} gray".format(aaid)) else: run(session, "show /{} cartoon".format(aaid, get_polymer_color(polyclass))) run(session, "color /{} {}".format(aaid, get_polymer_color(polyclass))) #--------------------------------------------[ This whole bit just colors by chain annotation basically] run(session, "graphics silhouettes true width 1") run(session, "light soft") # # ! Didn't work out for now. This is how i tried to make movies. (Loop over every structure (loop.py above) and call this command). Chimerax segfaults when looped on this. # def *produce_and_save_movie*(session, target:str): # print("GOT TARGET", target) # RCSB_ID = target # run(session, "open /home/rtviii/dev/RIBETL_DATA/{}/{}.cif".format(RCSB_ID, RCSB_ID)) # run(session, "sym #1 assembly 1") # take only one assembly if multiple are available # run(session, "ribrep #2") # run(session, "movie record") # run(session, "turn y 2 180") # run(session, "wait 180") # run(session, "movie encode /home/rtviii/dev/riboxyz/chimerax/movies/{}.mp4".format(RCSB_ID)) # run(session, "close all") # # ! This is just to not download the mmcif files i have locally. def *register_ribetl_command*(logger): def ribetl(session, rcsb_id:str): rcsb_id = rcsb_id.upper() run(session, "open /home/rtviii/dev/RIBETL_DATA/{}/{}.cif".format(rcsb_id, rcsb_id)) desc = CmdDesc( required= [("rcsb_id", StringArg)], ) register("ribetl", desc, ribetl, logger=logger) # # ! Register the above repr sequence as a command. def *register_ribrepr_command*(logger): from chimerax.core.commands import CmdDesc, register from chimerax.atomic import AtomicStructureArg, Chain, Residue, Atom desc = CmdDesc( required = [("structure", AtomicStructureArg)], required_arguments = ["structure"], synopsis = "representation ", ) register("ribrep", desc, ribosome_representation, logger=logger) *register_ribrepr_command*(session.logger) *register_ribetl_command*(session.logger) # This loop expands into a list of ~1800 PDB : 5AFI, 4UG0 and so on. for rcsb_id in os.listdir(RIBETL_DATA): run(session,"open {}".format(rcsb_id)) run(session,"ribrep #1") run(session,"save /home/student/dev/riboxyz/chimerax/{}.png width 800 height 800 transparentBackgroun d true".format(rcsb_id)) run(session,"close all") ``` Let me know if i can elucidate this situation further with hardware specs or whatnot. I'm quite at loss as to how to proceed with this.

2 2