- ChimeraX-users - RBVI Mailing Lists

Showing C-alphas
by Casey 30 Jul '21

30 Jul '21

I am trying to show a single sphere at c-alpha carbons for some point mutants that i designed. When I use the command: show /a:473,476,477 atoms I can see all the atoms in chain A for residues 473, 476, 477. Using: show /a:473,476,477@cb atoms I can show "two spheres" for the c-alpha and c-beta carbon. When I use: show /a:473,476,477@ca atoms nothing happens and i don't see any spheres. If I show the full atom and then do: hide /a:473,476,477@ca* atoms Just the c-alpha sphere gets hidden. Why is it that i can't just show the c-alpha spheres? Unfortunately for my figure i make silhouettes of my model and so i can't just hide the other atoms by making them fully transparent

2 1

Re: [chimerax-users] Automating Map segment
by Tom Goddard 30 Jul '21

30 Jul '21

To group by connectivity looks like you need to add d._group_con.value = True To figure out this stuff you would need to look at the Segger user interface code, it is in the ChimeraX distribution (under site-packages/chimerax/segger/segment_dialog.py) or on GitHub https://github.com/tomgoddard/segger/blob/master/Segger/src/segment_dialog.… <https://github.com/tomgoddard/segger/blob/master/Segger/src/segment_dialog.…> Tom > On Jul 30, 2021, at 10:42 AM, Sam W Marchant <smarchant(a)UCSD.EDU> wrote: > > Thanks so much im trying to group by connectivity and I am having trouble finding the right line of code to add that ability. Any thoughts? > >> On Jul 29, 2021, at 3:56 PM, Tom Goddard <goddard(a)sonic.net <mailto:goddard@sonic.net>> wrote: >> >> Hi Sam, >> >> I did an initial port of Segger from Chimera to ChimeraX (also called Segment Map) but that code was written by Greg Pintilie and he has not worked the ChimeraX version so there is no new command to run it in a script. It would probably not be too hard to write some Python to run it. The Chimera script mentioned in the mailing list is attached to my message about it. >> >> https://plato.cgl.ucsf.edu/pipermail/chimera-users/2017-October/013979.html <https://urldefense.proofpoint.com/v2/url?u=https-3A__plato.cgl.ucsf.edu_pip…> >> >> That won't work in ChimeraX, but here I attach the equivalent for ChimeraX. >> >> Segger was developed to segment molecules in EM maps, so I have doubts that it will do much good for your light microscopy. But you can decide that by running it by hand on some time steps and see if you like the result. >> >> Here is a different approach we used to quantify protrusions of crawling neutrophils in light sheet microscopy by placing markers that might be of interest. >> >> https://www.cgl.ucsf.edu/chimera/data/protrusion-sep2015/protrusions.html <https://urldefense.proofpoint.com/v2/url?u=https-3A__www.cgl.ucsf.edu_chime…> >> >> Tom >> >> >> >>> On Jul 29, 2021, at 3:24 PM, Sam W Marchant via ChimeraX-users <chimerax-users(a)cgl.ucsf.edu <mailto:chimerax-users@cgl.ucsf.edu>> wrote: >>> >>> Hi, >>> I am looking to apply the Map Segmenting to multiple time series of light sheet imaging and so I was wondering if there is either a command or a python function to automate the process. It looks here <https://urldefense.proofpoint.com/v2/url?u=http-3A__plato.cgl.ucsf.edu_pipe…> that you said it was going to be added to ChimeraX but I cannot find any info on how to do it. If not you mentioned a script in this answer, could you potentially send that to me? >>> Thanks, >>> Sam Marchant. >>> _______________________________________________ >>> ChimeraX-users mailing list >>> ChimeraX-users(a)cgl.ucsf.edu <mailto:ChimeraX-users@cgl.ucsf.edu> >>> Manage subscription: >>> https://www.rbvi.ucsf.edu/mailman/listinfo/chimerax-users <https://urldefense.proofpoint.com/v2/url?u=https-3A__plato.cgl.ucsf.edu_mai…> >> >> <segscript.py> >

1 0

How to Align Chromosome Models with Different Number of Atoms?
by Cardiff Jiang 30 Jul '21

30 Jul '21

Dear ChimeraX User List, How can I compare two chromosome 3D models (PDB files) that have different numbers of "atoms"? Matchmaker tool reports, "Reference and/or match model contains no nucleic or amino acid chains. Use the command-line 'align' command to superimpose small molecules/ligands." align #1 toAtoms #2 reports, "Unequal number of atoms to pair, 49 and 63." align #1 toAtoms #2 matchNumbering true reports, "Pairing dropped 49 atoms and 63 reference atoms. No atoms paired for alignment." align #1:1-40 to #2:1-40 reports, "No atoms paired for alignment." How should I compare two chromosome 3D models below? Thank you. [image: image.png] Best regards, Zichen "Cardiff" Jiang

4 9

Automating Map segment
by Sam W Marchant 29 Jul '21

29 Jul '21

Hi, I am looking to apply the Map Segmenting to multiple time series of light sheet imaging and so I was wondering if there is either a command or a python function to automate the process. It looks here<http://plato.cgl.ucsf.edu/pipermail/chimera-users/2017-October/013980.html> that you said it was going to be added to ChimeraX but I cannot find any info on how to do it. If not you mentioned a script in this answer, could you potentially send that to me? Thanks, Sam Marchant.

2 1

Re: [chimerax-users] command line way for inverting selection
by Rangarajan, Amith 29 Jul '21

29 Jul '21

Hi, I am looking for a way to invert selection without reaching for the tool bar. I would like to know if there is a way to invert selection via the select command ( or some other command) thanks a lot, Amith

3 3

ATOM has no alt loc B
by James Loy 28 Jul '21

28 Jul '21

Hi folks- I'm trying to add hydrogens to an mmCIF file (5byl) using the addh command. However, every time I try, I encounter the following error: Traceback (most recent call last): File "/usr/lib/ucsf-chimerax/lib/python3.8/site-packages/ChimeraX_main.py", line 712, in init run(sess, cmd) File "/usr/lib/ucsf-chimerax/lib/python3.8/site-packages/chimerax/core/commands/run.py", line 36, in run results = command.run(text, log=log, return_json=return_json) File "/usr/lib/ucsf-chimerax/lib/python3.8/site-packages/chimerax/core/commands/cli.py", line 2852, in run result = ci.function(session, **kw_args) File "/usr/lib/ucsf-chimerax/lib/python3.8/site-packages/chimerax/addh/cmd.py", line 63, in cmd_addh add_h_func(session, structures, template=template, in_isolation=in_isolation, **prot_schemes) File "/usr/lib/ucsf-chimerax/lib/python3.8/site-packages/chimerax/addh/cmd.py", line 170, in hbond_add_hydrogens add_hydrogens(session, atoms, type_info_for_atom, naming_schemas, hydrogen_totals, File "/usr/lib/ucsf-chimerax/lib/python3.8/site-packages/chimerax/addh/hbond.py", line 214, in add_hydrogens _attach_hydrogens(atom, altloc_hpos_info, bonding_info) File "/usr/lib/ucsf-chimerax/lib/python3.8/site-packages/chimerax/addh/hbond.py", line 895, in _attach_hydrogens add_altloc_hyds(atom, altloc_hpos_info, invert, bonding_info, total_hydrogens, naming_schema) File "/usr/lib/ucsf-chimerax/lib/python3.8/site-packages/chimerax/addh/cmd.py", line 862, in add_altloc_hyds h = new_hydrogen(atom, i+1, total_hydrogens, naming_schema, File "/usr/lib/ucsf-chimerax/lib/python3.8/site-packages/chimerax/addh/cmd.py", line 882, in new_hydrogen metal_pos = metal.get_alt_loc_coord(alt_loc) File "atomic_cpp/cymol.pyx", line 587, in chimerax.atomic.cymol.CyAtom.get_alt_loc_coord ValueError: Atom /A MG 702 MG has no alt loc B I'm using ucsf-chimerax_1.2.5-1_amd64.deb installed on Ubuntu 20.04. I've done some digging, and I cannot figure out what is going on. Any help would be appreciated. Thanks in advance! -James -- James Loy, Ph. D. 614.915.9792 LinkedIn <https://www.linkedin.com/in/james-m-loy-ph-d-ba2256a/>

3 4

The PDB code for the demo picture on the website
by BAOCHENG LIU 28 Jul '21

28 Jul '21

Hi chimerax group, I am a postdoc in UCLA using chimerax. Thanks for providing all scientists this powerful software that makes life easier. I have a question here when I was learning the cartoon/ribbon style command line. So, I saw the 4 demo pictures on the website https://www.rbvi.ucsf.edu/chimerax/docs/user/commands/cartoon.html#style and just realized the structure just looks similar to one of my structures. Could you please tell me the PDB code for that structure so I can compare it with mine in detail? Hope it is released. Thank you. Best regards, Baocheng

2 1

resfit command for nuclei acid
by Wen-Jin Winston Wu 27 Jul '21

27 Jul '21

Hi Tom, I have just tried to use the "resfit" command ChimeraX Does the nice command of "resfit" in ChimeraX support inspection on nuclei acid? It seems that it only support protein residues as with the command " resfit #1/D:DT map #3", it returned with the message that "None of the 9 specified residues have backbone atoms N, CA and C". I have not checked out the daily build version yet, but thought that I can ask first. Thanks, Best wishes, Winston

2 2

Custom colour for atoms with bfactor==0
by Lorena Zuzic 27 Jul '21

27 Jul '21

Hi, I wanted to know if there is a way of colouring my structure by bfactor, but with an additional condition where all bfactor==0 atoms are coloured in grey. Currently, I am using a white to red gradient for my B-factors: color bfactor range -2.2, 2.2 In this case, bfactor==0 values are coloured white. However, I would like them to be grey and outside the gradient (as they represent an unmapped portion of the protein). Thanks! Lorena

2 2

create mask from "shapes"
by Roberto Marabini 27 Jul '21

27 Jul '21

Hi, It is possible to create a mask from a shape?. For example from the result of "*shape icosahedron".* If the command does not exist I guess I may create it as follows * the input would be the 3D map to be masked and the shape to be used as mask * request all the (varray, tarray) from the shape. * from varray,tarray build the triangles that define the surface of the shape * for each triangle create the normal vector * multiply each normal vector by the position of each voxel (dot product) * keep those voxels that have negative dot products for all the triangles. thanks for the help

2 2