- ChimeraX-users - RBVI Mailing Lists

Re: [chimerax-users] recent query
by Rangarajan, Amith 26 Jan '23

26 Jan '23

Dear Eric, thanks for taking time to answer my query. It was very detailed and i am beginning to understand how bundles work on chimerax. i am working on a methods paper that uses sidechain center of masses as a representation of mutation position and generates clusters to intrepret the impact of the mutations. i tried scripting this, but i realised that pandas library is not present in chimerax (?) and i did not proceed further. is there a way to include scikit learn , pandas libraries in chimerax script? thanks once again for all your support ! best amith

2 3

Macbook Pro, fans, and ChimeraX
by Lewis James martin 26 Jan '23

26 Jan '23

Hi all, Ever since I got my Macbook Pro (16-inch, 2019, 2.3GHz 8‑core Intel Core i9) I noticed the fans can go full bore for seemingly small tasks. ChimeraX is one of those - swinging a protein view around is enough to set them off, and Activity Monitor reports around 200% usage for the "UCSF ChimeraX" process. I can see that GPU usage on the inbuilt AMD Radeon Pro 5500M can get up to about 50%, too. The bottom can get quite hot - it's Summer here now, but this happens in Winter too. A couple of questions: - Is this expected (for ChimeraX, not generally)? Do other people experience this? - Does ChimeraX use the GPU for rendering? If so, what's driving the CPU spike when I move the protein on screen? Is it just from calculating the coordinate rotations? Thanks, Lewis

2 3

Re: [chimerax-users] [Chimera-users] Rotation along symmetry axys
by Tom Goddard 25 Jan '23

25 Jan '23

Hi Leonardo, As Elaine said, rotating the view in ChimeraX does not change what is saved in the map (*.mrc) file. The data needs to be rearranged in the 3D array. That can be done with the ChimeraX "volume resample" command (in Chimera it is called vop resample) but that interpolates the data on a new grid, and it would be better to just rotate the 3D data array so you don't get interpolation errors. I made a ChimeraX command to do that "volume rotate_90" available on the ChimeraX Recipes web site. https://rbvi.github.io/chimerax-recipes/rotate_90/rotate_90.html We no longer develop Chimera so I provide this for ChimeraX. Tom > On Jan 25, 2023, at 6:14 AM, Leonardo Talachia Rosa via Chimera-users <chimera-users(a)cgl.ucsf.edu> wrote: > > Dear Chimera team, > I hope you are well. > > I am processing a phage dataset which presents D5 pseudossymetry. The C5 axys is aligned to the z axys. I would like to rotate the map 90 degree in the z axys, so that Z axys now separate the molecule in a C2 pseudosymmetry, in order to play with some symmetry relaxation features to refine differences in the two map vortexes. > Is there a way in chimera to rotate the map in relation to the xyz axys? The way I see, all rotate commands I know rotate the axys themselves alongside the map. > > Thank you for your time, > Kind wishes, > Leonardo > > -- > ************************* > Leonardo TALACHIA ROSA, PhD > PostDoctoral Researcher > Microbiology and Structural Biology > Instituto de Química, Universidade de São Paulo - IQ-USP > Av. Prof. Lineu Prestes, 748 - Butantã, São Paulo - SP > ************************** > _______________________________________________ > Chimera-users mailing list: Chimera-users(a)cgl.ucsf.edu > Manage subscription: https://www.rbvi.ucsf.edu/mailman/listinfo/chimera-users

1 0

color palette in transprency ?
by vincent Chaptal 25 Jan '23

25 Jan '23

Hi, I'm rendering a series of models using this command line (which works very well): rainbow #1 structures palette cyan:blue I would like to do the same but in a series of transparent models, with a range from 0 to 80 (or similar). Is there a command for this, I couldn't find it? Thanks a lot. Vincent -- Vincent Chaptal, PhD Director of GdR APPICOM Drug Resistance and Membrane Proteins Lab MMSB -UMR5086 7 passage du Vercors 69007 LYON FRANCE +33 4 37 65 29 01 http://www.appicom.cnrs.fr http://mmsb.cnrs.fr/en/

2 4

vector
by Dr.S.Somasundaram 25 Jan '23

25 Jan '23

2 2

Pseudobonds change style
by Samo Lešnik 25 Jan '23

25 Jan '23

Dear ChimeraX developers, I am wondering if it is possible to change the style of pseudobonds, specifically "Distances", from "dotted disks" to "whole undotted lines"? It is very easy to do this in Chimera, with just changing the style in the graphical interface of Distances. Thank you very much, dr. Samo Lešnik

2 2

Re: [chimerax-users] downloading pdb files from another webserver
by Rangarajan, Amith 24 Jan '23

24 Jan '23

hi, i mostly resort to using the resource from dunbrack's lab to download structures renumbered as per uniprot numbering. i currently do it using a simple .py script: [cid:14eb2e7f-deaa-4341-92eb-3f9c56db264b] i want to make this as a chimerax script that i can call from command line with open {scriptname.py} { pdbid} wondering if some one has tried to do similar things and some tips regarding which modules to use might be helpful. thanks a lot , best amith

3 3

Memory problems with Orbital Viewer
by Michael Patzschke 24 Jan '23

24 Jan '23

Dear colleagues, I am using ChimeraX mostly for small molecular complexes, typically up to 100 atoms. I am running the code on a MacBook Pro with 16 GB of memory under macOS 11.8.8. When I try to render molecular orbitals from an orca calculation I get even for very small molecules the warning message: Estimated peak memory usage (0.0GB) is above or close to the available memory (0.0GB). Exceeding available memory might affect the stability of your computer. You may attempt to continue, but it is recommended that you lower your resolution, decrease padding, or use fewer threads. Press "Ok" to continue or "Abort" to cancel Is there any way to fix this problem? Kind regards, Michael Patzschke

2 4

electrostatic surface of docking model
by Citterico, Matteo 23 Jan '23

23 Jan '23

Hello everyone, I am trying to visualize the electrostatic potential surface of a docking model generated with ClusPro2.0 (job: https://cluspro.bu.edu/models.php?job=898454) When I click on the electrostatic button on ChimeraX 1.3 to visualize the charge, it gives the error: "Using Amber 20 recommended default charges and atom types for standard residues. Hydrogen copy of model.000.00.pdb #/A PHE 244 H1 bonded to atom that should not have hydrogens (copy of model.000.00.pdb #/A PHE 244 C)." How can I solve this problem? Thank you for the help and always a big up to the Chimera staff! Matteo

2 6

Coot view in ChimeraX
by Cryo EM 23 Jan '23

23 Jan '23

Hi everyone, Does anyone know of a command for getting a coot view of map in ChimeraX? I remember seeing it in chimera, so I guess there should be something of such sort in ChimeraX too? Regards and thanks

2 1

Error in opening Chimera
by Ashfaque Rahman 22 Jan '23

22 Jan '23

Why is this error showing? What is the solution? [image: image.png]

2 1

Non contiguous residues selection from command line
by LEONARDO SORCI 22 Jan '23

22 Jan '23

Hello there! I wondered if there is a way to quickly select noncontiguous residues from the command line (useful, for example, when selecting a set of residues defining an active site, a template for docking, etc.). In standard chimera, I would use the command SELECT with numbers, commas, and hyphens (e.g., SELECT: 172,192,193,208:214). Thanks, Leonardo Polytechnic University of Marche, Ancona, Italy

2 2

How to select best pose from docking model?
by Luigi Marongiu, PhD 22 Jan '23

22 Jan '23

Dear all, I am new to docking analysis. I have obtained a list of poses from VINA analysis, but what of them shall I choose? Is the first the best choice? The problem is that, according to the pose, the ligand binds the protein in a different position. Different ligands might or might not bind on the same position, thus choosing the right pose is essential to determine whether there is a position that allows binding of multiple molecules. What is the procedure to determine whether such a multiple binding site exists? What is the protocol to choose the right pose? Thank you

2 1

Re: [chimerax-users] overwrite session aliases from the startup options
by Helge Paternoga 20 Jan '23

20 Jan '23

Hi Elaine, Thanks for the tip with "alias delete all", I had not considered that one. My solution so far was to open the startup screen and copy all new aliases, then paste them in the command line and hit enter which overwrites everything with the newest ones. So a toggle option to "re-run startup script upon session loading" or something similar could alleviate the issue. Best, Helge Am 18.01.23 um 21:00 schrieb chimerax-users-request(a)cgl.ucsf.edu: > Dear Helge, > I don't have an elegant solution, but faced with this problem I would > > (1) open old session > > (2) use "alias delete" to get rid of the old versions of the aliases (or "alias delete all" if you are OK with getting rid of all of them) > <https://rbvi.ucsf.edu/chimerax/docs/user/commands/alias.html> > > (3) re-save session > > I understand that it is a pain to do this if you have a large repertoire of saved sessions, however. > I hope this helps, > Elaine

2 4

ChimeraX image saving automation
by Attana, Fedaa 19 Jan '23

19 Jan '23

Hi all, I have started to use Chimera recently and found it very cool. I liked the way I can edit all of my structure with one command when no model is selected. However, to save images, I have to manually show each structure (model) and run the same saving command with a new name. I wonder if there is an easier way to do this or if somebody did automate this process so one can save images of individual models with a script? Thanks in advance, Fedaa

3 2

scalemol extension not working
by André Graça 19 Jan '23

19 Jan '23

Hi everyone, I hope to at least get a reply from Tom, but this may need a bit of troubleshooting (?). I looked into a very old post of yours, Tom, about rescaling an atomic model (see it at the end of my email). I put everything in place but I get an error and I wonder if anything has changed with the PDB format since 2005, or some incompatibility with new chimera versions... If there is actually a way of doing the same with ChimeraX, that would be incredible =) Cheers! André Error and forwarded message comes now: __________________________________________________ AttributeError Exception in Tk callback Function: <bound method MidasUI.processCommand of <Midas.midas_ui.MidasUI instance at 0x00000000079092C8>> (type: <type 'instancemethod'>) Module: <module 'Midas.midas_ui' from 'C:\Program Files\Chimera 1.16\share\Midas\midas_ui.pyc'> (line: 618) Args: (<Tkinter.Event instance at 0x00000000198D8AC8>,) Event type: KeyPress (type num: 2) Traceback (innermost last): File "C:\Program Files\Chimera 1.16\bin\lib\site-packages\Pmw\Pmw_1_3_3\lib\PmwBase.py", line 1747, in __call__ return apply(self.func, args) File "C:\Program Files\Chimera 1.16\share\Midas\midas_ui.py", line 628, in processCommand midas_text.makeCommand(cmdText) File "C:\Program Files\Chimera 1.16\share\Midas\midas_text.py", line 69, in makeCommand f(c, args) File "C:\Program Files\Chimera 1.16\share\scalemol\ChimeraExtension.py", line 23, in scale_molecule x, y, z = a.coord().xyz.data() AttributeError: '_chimera.Point' object has no attribute 'xyz' ================================================ Event contents: char: delta: 0 height: ?? keycode: 13 keysym: Return keysym_num: 65293 num: ?? serial: 14638 state: 8 time: 154527203 type: 2 widget: .121252488L.126915336L.428746824L.428747592L.428747720L width: ?? x: 228 x_root: 288 y: -22 y_root: 941 ______________________________________________________________________ ------- Start of forwarded message ------- Date: Wed, 27 Apr 2005 11:05:55 -0700 (PDT) From: Thomas Goddard <goddard at cgl.ucsf.edu<http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users>> Hi Jinghua, By "virus cage pdb file" I guess you mean you placed atoms at icosahedral positions, and want to show them connected up to show some icosahedral pattern. I think the easiest approach is to open your old session, run a command to change the diameter of the cage, then save a new session. Below is some Python code that provides a "scalemol" Chimera command. Make a directory chimera/share/scalemol in your Chimera distribution and put the code below in file chimera/share/scalemol/ChimeraExtension.py. Then start Chimera, open the cage PDB model, use menu entry Favorites/Command line, and type a command: scalemol 0 1.5 to make model number 0 bigger by a factor of 1.5. If you have more than one model open you can see the model numbers using Favorites / Model Panel, the model number is shown in the left hand column. Tom - ---- scalemol/ChimeraExtension.py code follows: # ----------------------------------------------------------------------------- # def scale_molecule(cmdname, args): from Midas.midas_text import error fields = args.split() try: model_num = int(fields[0]) factor = float(fields[1]) except: error('Syntax error: scalemol <model-number> <factor>') return import chimera mlist = chimera.openModels.list(id = model_num) if len(mlist) == 0: error('scalemol: No model number %d' % model_num) return for model in mlist: for a in model.atoms: x, y, z = a.coord().xyz.data() c = chimera.Coord() c.x, c.y, c.z = (factor*x, factor*y, factor*z) a.setCoord(c) # ----------------------------------------------------------------------------- # import Midas.midas_text Midas.midas_text.addCommand('scalemol', scale_molecule) ------- End of forwarded message -------

2 1

missing or invalid sequences
by Liang, Junyi 19 Jan '23

19 Jan '23

I would like to predicate interaction between two proteins, and tried chimerax version 1.5 after noticing a youtube video https://www.youtube.com/watch?v=6lXeCPuTePs. [https://www.bing.com/th?id=OVP.95Ygyv4AXQTgCS0WB9iXtgHgFo&pid=Api]<https://www.youtube.com/watch?v=6lXeCPuTePs> Running AlphaFold to Predict Protein Complexes from ChimeraX<https://www.youtube.com/watch?v=6lXeCPuTePs> How to predict the structure of a photoreaction center complex of 3 proteins from ChimeraX using AlphaFold-Multimer. Requires a ChimeraX daily build newer than December 2, 2021 or version 1.4. It is not in ChimeraX version 1.3. Web page describing steps https://www.rbvi.ucsf.edu/chimerax/data/alphafold-dec2021/af_multimer.html www.youtube.com I input the two protein sequences in same way as shown in the video and hit predict. But your software says "missing or invalid sequences". How come? Please consider the environment before printing this e-mail Cleveland Clinic is currently ranked as the No. 2 hospital in the country by U.S. News & World Report (2017-2018). Visit us online at http://www.clevelandclinic.org for a complete listing of our services, staff and locations. Confidentiality Note: This message is intended for use only by the individual or entity to which it is addressed and may contain information that is privileged, confidential, and exempt from disclosure under applicable law. If the reader of this message is not the intended recipient or the employee or agent responsible for delivering the message to the intended recipient, you are hereby notified that any dissemination, distribution or copying of this communication is strictly prohibited. If you have received this communication in error, please contact the sender immediately and destroy the material in its entirety, whether electronic or hard copy. Thank you.

2 5

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