January 2023 - ChimeraX-users - RBVI Mailing Lists

Is RR Distance Maps implemented in ChimeraX?
by Michael Bruist 19 Jan '23

19 Jan '23

I use RR Distance Maps in Chimera for my teaching. I want to switch to ChimeraX, but I have not bee able to find this tool. Is it available in ChimeraX? Mike Bruist Michael F. Bruist Professor of Biochemistry Department of Chemistry and Biochemistry St. Joseph’s University 215-596-8530 mbruist(a)sju.edu

2 1

Reasign secondary structure for back-bone trace
by Enrico Martinez 19 Jan '23

19 Jan '23

Dear Chimera-X users! I am working on the visualization of a trajectory containing only back-bone atoms of the protein obtained from the PCA. When I load this multi-pdb ensemble to Chimera-X I can see only back-bone traces without any secondary structure information. I've tried to use dssp command to re-assign the secondary structure but it did not change the visualization state. Would it be possible to do it in some automatic fashion via some built-in command? Many thanks in advance Cheers Enrico

2 4

incomplete output from matchmaker embedded in python script
by Borchert, Milan 19 Jan '23

19 Jan '23

Hello ChimeraX community, I have been trying to use the Chimerax tool to sequentially check one query pdb against multiple (later up to 10000) pdb files. For this I have been working with the python scripting functionality. With the following code I was able to get the resulting data (number of matched atoms and resulting RMSD) saved as a html log file: import os from chimerax.core.commands import run as rc # change to folder with data files os.chdir("/home/ubuntu/Desktop/Chimerax_search/N_crassa_pdb") # gather the names of .pdb files in the folder file_names = [fn for fn in os.listdir(".") if fn.endswith(".pdb")] # give path to query file to be searched against pdb library query = "/home/ubuntu/Desktop/Chimerax_search/query/BDFB_007497.pdb" # open the query file as #1 in ChimeraX rc(session, "open /home/ubuntu/Desktop/Chimerax_search/query/BDFB_007497.pdb") # loop through the files for fn in file_names: # open first target pdb as #2 rc(session, "open " + fn) # match both structures rc(session, "match #1 to #2") # save log with ID matching input filename, close #2 and clear log rc(session, "log save " + fn.split(".")[0] + ".log") rc(session, "close #2") rc(session, "log clear") Instead of trying to gather the relevant data from thousands of rather large html output files I instead would prefer to only get then resulting values from matchmaker. I learned that the command actually has all the values I wanted as return values, so I tried to write a modified script, that saves these values in a variable and collects the values for every target pdb file in a dict: import os from chimerax.core.commands import run as rc # change to folder with data files os.chdir("/home/ubuntu/Desktop/Chimerax_search/N_crassa_pdb") # gather the names of .pdb files in the folder file_names = [fn for fn in os.listdir(".") if fn.endswith(".pdb")] # give path to query file to be searched against pdb library query = "/home/ubuntu/Desktop/Chimerax_search/query/BDFB_007497.pdb" # open the query file as #1 in ChimeraX rc(session, "open /home/ubuntu/Desktop/Chimerax_search/query/BDFB_007497.pdb") # create empty dict to collect matchmaker output output = {} # loop through the files, opening, processing, and closing each in turn for fn in file_names: id = fn.split("/")[-1].split(".")[0] rc(session, "open " + fn) match = rc(session, "match #1 to #2") output.update({id : match}) rc(session, "log save " + fn.split(".")[0] + ".log") rc(session, "close #2") rc(session, "log clear") # display dict print(output) For now I just used 4 target pdbs with my query and only tried printing the result to quickly check if it fits. However the result is a bit baffling to me. I do get the RMSD values, but not the other data: {'NCU01382': [{'full match atoms': <chimerax.atomic.molarray.Atoms object at 0x7f1944492130>, 'full ref atoms': <chimerax.atomic.molarray.Atoms object at 0x7f194448f610>, 'final match atoms': <chimerax.atomic.molarray.Atoms object at 0x7f1946694430>, 'final ref atoms': <chimerax.atomic.molarray.Atoms object at 0x7f1946694550>, 'full RMSD': 63.065617285514556, 'final RMSD': 1.365365439813019, 'transformation matrix': <chimerax.geometry.place.Place object at 0x7f194669c130>, 'aligned ref seq': <chimerax.atomic.molobject.Sequence object at 0x7f19f7b8f8b0>, 'aligned match seq': <chimerax.atomic.molobject.StructureSeq object at 0x7f19f7b8f610>} here is one of the resulting hits, for my first target pdb (NCU01382). As you can see, the RMSD scores are correctly given, but for other values the internal object names and no resulting value are displayed. Any help with figuring this out would be greatly appreciated. Best Wishes Milan Borchert

2 2

select displayed atoms via command line / transfer saved views
by Helge Paternoga 18 Jan '23

18 Jan '23

Dear ChimeraX users, I have to unrelated questions: 1) Is there a way to select displayed atoms via the command line? 2) Is it possible to transfer saved views from one session to another one? Many regards and thank you for reading, Helge

3 3

Clipper's strange behaviour
by André Graça 16 Jan '23

16 Jan '23

Hi everyone! Is there any known issue with Clipper plugin for ChimeraX? After installing clipper and restarting ChimeraX, the software is not responding after I execute the Open command. I expect that it takes some time to open an .mtz file as the map needs to be calculated, however I have waited more than 45 minutes and ChimeraX keeps unresponsive... I am using Chimera 1.5 and I installed the (hopefully compatible) version of Clipper 0.19.0 Did anyone faced the same problem? Let me know if you have a fix or trick ;) Thanks for any feedback in advance, André ____________________ [cid:image001.png@01D9272E.120B2740]

2 2

Movie clipping planes in chimera x
by Greg Martin 13 Jan '23

13 Jan '23

Hi, I’m trying to make a movie where I highlight different parts of a large complex. For one scene, I’m trying to interpolate between a zoomed out view of the whole complex (view p1), and then a zoomed in view of an interface (view p2), which shows a series of interactions. What I want is a smooth transition of both the view and clipping plane from p1 to p2. However, when I give the command: view p1 p2 100 it interpolates the two positions/zooms fine, in 100 frames, but then after the movement adjusts the clipping plane all at once. Is there a way to have the plane adjusted smoothly during transition between p1 and p2? I’ve also tried fly, but this disregards clipping planes altogether. Views are saved with view name p1 view name p2 Thanks in advance, Greg -- Gregory Martin, PhD. Staff Scientist Lab of Andrew Ward Department of Integrative Structural & Computational Biology The Scripps Research Institute La Jolla, CA 92037

3 3

Inquiry: Map sequence conservation onto a protein structure
by Gomez Barillas, Andrea 13 Jan '23

13 Jan '23

Dear ChimeraX developers, I am trying to map sequence conservation onto a protein structure. I am using ChimeraX version (version 1.4 (2022-06-03). For this purpose, I was trying to use a tutorial from 2014 for Chimera for Mapping Sequence Conservation onto Structures (Showing Sequence Conservation in Chimera (ucsf.edu)<https://www.cgl.ucsf.edu/chimera/data/tutorials/systems/outline.html>). It seems that some the tools mentioned in this tutorial are no longer available, or maybe work differently. Specifically, the tutorial's "Case 4" (calculating values outside of Chimera but showing them in Chimera) has you set an attribute, and then render/color the structure based on the attribute. The Tools > Structure Analysis > Define attribute mentioned in the tutorial is no longer available. I was as able to assign values directly to structure residues in a text file and read it in using the command defattr instead. The next step would be to use the tool Render by Attribute for coloring but I cannot find this option anymore. Is there a different command I can use for this step, or a new functionality built into ChimeraX I could use for my purpose of mapping sequence conservation onto the structure? One relevant point is that I already have my sequence identify values calculated from another program, and I want to import those values. Thank you in advance for your help, Andrea Showing Sequence Conservation in Chimera<https://www.cgl.ucsf.edu/chimera/data/tutorials/systems/outline.html> Example 2B Chimera session: ConSurf-1hd2-chimera19.py [show this last!] (result of submitting 1hd2 to the ConSurf server and choosing to show results in Chimera; alignment has 150 sequences including that of the query, 1hd2) . Given a protein structure, the ConSurf server estimates the evolutionary conservation of amino acid positions based on the phylogenetic relations between homologous ... www.cgl.ucsf.edu

2 1

Inquiry: Map sequence conservation onto a protein structure
by Gomez Barillas, Andrea 13 Jan '23

13 Jan '23

Dear ChimeraX developers, I am trying to map sequence conservation onto a protein structure. I am using ChimeraX version (version 1.4 (2022-06-03). For this purpose, I was trying to use a tutorial from 2014 for Chimera for Mapping Sequence Conservation onto Structures (Showing Sequence Conservation in Chimera (ucsf.edu)<https://www.cgl.ucsf.edu/chimera/data/tutorials/systems/outline.html>). It seems that some the tools mentioned in this tutorial are no longer available, or maybe work differently. Specifically, the tutorial's "Case 4" (calculating values outside of Chimera but showing them in Chimera) has you set an attribute, and then render/color the structure based on the attribute. The Tools > Structure Analysis > Define attribute mentioned in the tutorial is no longer available. I was as able to assign values directly to structure residues in a text file and read it in using the command defattr instead. The next step would be to use the tool Render by Attribute for coloring but I cannot find this option anymore. Is there a different command I can use for this step, or a new functionality built into ChimeraX I could use for my purpose of mapping sequence conservation onto the structure? One relevant point is that I already have my sequence identify values calculated from another program, and I want to import those values. Thank you in advance for your help, Andrea

1 0

ChimeraX Log window remains blank
by Scott White 13 Jan '23

13 Jan '23

Hi, Sorry if this is a dumb question, but recently when starting ChimeraX the Log window is not on Display. Opening the Log window (Tools> log) opens it as a separate window, and whether docked or free-floating the Log window remains blank no matter what I do. Even just "open 1nov" on the command line will download and open the PDB entry 1NOV, but no text appears in the Log window. I've tried a fresh install of ChimeraX 1.5 to no avail. I've tried opening an older version of ChimeraX, which despite having worked countless times before, now shows the same behaviour regarding the Log window as my abberant ChimeraX 1.5. I've tried rebooting the computer - no fix. I've tried "log clear" and "log show". It's a Windows 10, university owned PC. I don't have the issue on other computers. Is it a software clash? Thanks in advance for any help or advice. Scott _______________________________________ Dr Scott A. White SFHEA School of Biosciences, Room T716 University of Birmingham, B15 2TT E: s.a.white(a)bham.ac.uk<mailto:s.a.white@bham.ac.uk> T: +44-(0)121-414-7534 Office Hours Generally Tues 10am - 11am, Wed 3 pm - 4 pm Please contact me for an appointment. ______________________________________

2 2

Changing MD trajectory default frame rate
by Visvaldas K. 12 Jan '23

12 Jan '23

Dear all, Is there a way to change frame rate when I hit that red "record" button when playing MD trajectory? I could do "movie record", but then it becomes interactive (I hit triangle "play", and then "movie encode /Users/visvaldas/Desktop/movie.mp4 framerate 5"), so inadvertently can record some undesired frames. Of course I probably can reduce framerate using external programs, but how to that non-interactively in ChimeraX? Best regards, Vis

3 3