September 2021 - ChimeraX-users - RBVI Mailing Lists

Question about use of ChimeraX
by Jonathan Choong 07 Sep '21

07 Sep '21

Hi there, Hope you are doing well. I just had a question regarding ChimeraX that I was hoping to get some clarity on. I was wondering whether it was possible to use ChimeraX commands (such as open, run, save) in the system command line without actually opening the application interface? Thanks, hope to hear back from you. Kind regards, Jonathan

2 1

Re: [chimerax-users] [Chimera-users] matchmaker for more than one chain
by Elaine Meng 03 Sep '21

03 Sep '21

Hi Mani, There is a different e-mail address for ChimeraX questions: it is chimerax-users(a)cgl.ucsf.edu CC'd here. Responses below... > On Sep 2, 2021, at 2:22 AM, manikandank via Chimera-users <chimera-users(a)cgl.ucsf.edu> wrote: > > Hi ChimeraX-users, > I have used recently the mmaker cmd to align two pdb each have multiple chains as below. > > mmaker #12/L,M,N to #5/S,A,B matrix blosum-62 showAlignment true > > But the sequence view shows only the alignment of chain L and S. > > 1. I would like to have aligned sequences of all chains used in mmaker alignement. Your command above does not use all the chains for fitting, it only chooses the single pair of chains that gives the best sequence-alignment score to superimpose the two models. The other chains might look superimposed but it is only because they moved along with the ones that were actually used for the fit. To use all three pairs of chains, you need to also use the "pairing ss" option: mmaker #12/L,M,N to #5/S,A,B matrix blosum-62 showAlignment true pair ss Then you would get three pairwise sequence-alignment windows: L/S, M/A, and N/B. matchmaker or mmaker command <https://rbvi.ucsf.edu/chimerax/docs/user/commands/matchmaker.html> pairing options <https://rbvi.ucsf.edu/chimerax/docs/user/commands/matchmaker.html#options> > > 2. Is it possible to open the saved RMSD header values in a file later just to color the pdb? Yes, but with several steps so it may not be that convenient. Probably easier just to save and restore the ChimeraX session (which includes the sequence alignment windows). Also remember if you move the structures relative to each other, the RMSD values will change. However, if you really want to save the RMSD values to file after you showed sequence alignments and header, then you could save them as an attribute file (that could be reopened/applied later): <https://rbvi.ucsf.edu/chimerax/docs/user/commands/save.html#attributes> The attribute name is seq_rmsd as explained here: <https://rbvi.ucsf.edu/chimerax/docs/user/tools/sequenceviewer.html#headers> > > 3. Do mmaker align entire pdb based on selected chains and Can calculate RMSD for all residues in the pdb? For the first question, 3A: yes, depending on the "pairing" option as explained above. For the second question, 3B, I'm not sure what you want. Remember that the main purpose of matchmaker is to superimpose the structures well, and then it happens to also give the RMSD over the positions that it chose to use, which are shown in the pinkish boxes on the sequence alignment. By default it uses fit iteration and discards loops that don't superimpose well in space. Matchmaker only gives the RMSD over the positions it uses. You could tell it to use all positions, but that would change the superposition. If you like the default superposition (that does not use all of the positions), you can use a separate command "rmsd" to measure RMSD over any sets of atoms that you want, without moving the structures or changing the superposition. Or, if you really want to use all of the positions in the sequence alignment for the fitting, you would have to use matchmaker option "cutoff none" ... but that usually does not give the best superposition. cutoff option <https://rbvi.ucsf.edu/chimerax/docs/user/commands/matchmaker.html#cutoffDis…> rmsd command (calculate RMSD in current positions, no fitting) <https://rbvi.ucsf.edu/chimerax/docs/user/commands/rmsd.html> Matchmaker only gives alpha-carbon RMSD. The rmsd command can give RMSD for any sets of atoms that you want, alpha-carbons only, or backbone only, or all atoms, as long as you give equal numbers of atoms from the two structures. However, it can be tricky to specify the correct sets of atoms if the residue numbering or residue types are different between the two structures. So you might want to first use "mmaker" with whichever options give you the best superposition, and then to measure RMSD over all positions in the alignment (which might not be used in the matchmaker fit), use the "rmsd" command afterward. > > 4. I also would like to change the displayed text for the sequence name. Presently it puts the filename as sequence name. > > [similar to 'Edit sequence name' in the sequence viewer in chimera] Sorry, there is no way yet to change the displayed sequence name in ChimeraX. > > Please let me know how this can be done. > thank you very much > with kind regards > Mani. I hope this helps, Elaine ----- Elaine C. Meng, Ph.D. UCSF Chimera(X) team Department of Pharmaceutical Chemistry University of California, San Francisco

2 1

bundle dependency not installed
by Lagerwaard, I.M. (Ilse) 02 Sep '21

02 Sep '21

Hello, I'm currently developing a bundle with a dependency on pandas. Even though I included this dependency in my bundle_info.xml file, I get "ModuleNotFoundError: No module named 'pandas'". For now, I have installed pandas on the ChimeraX console via the command prompt. For future distribution, however, it would be nice if it were immediately installed together with the bundle. I would greatly appreciate your help with this issue. Thanks in advance, Ilse

2 1

Re: [chimerax-users] output mmCIF file
by Elaine Meng 02 Sep '21

02 Sep '21

> On Aug 31, 2021, at 1:50 PM, Liu, Xiangan <Xiangan.Liu(a)uth.tmc.edu> wrote: > > Hi, > I have issue to output the models I build in pdb formation to mmCIF formation. I did not provide the sequences in orginal pdb files. Could please let me know what is the problem and how to fix it? > Thanks, > Xiangan > ---------- > File "/usr/lib/ucsf-chimerax/lib/python3.8/site-packages/chimerax/mmcif/__init__.py", line 89, in save > mmcif_write.write_mmcif(session, path, **kw) > File "src/mmcif_write.pyx", line 145, in chimerax.mmcif.mmcif_write.write_mmcif > File "src/mmcif_write.pyx", line 149, in chimerax.mmcif.mmcif_write.write_mmcif > File "src/mmcif_write.pyx", line 542, in chimerax.mmcif.mmcif_write.save_structure > KeyError: ('a1', 'RPGMMDSQEFS') > > KeyError: ('a1', 'RPGMMDSQEFS') > > File "src/mmcif_write.pyx", line 542, in chimerax.mmcif.mmcif_write.save_structure Hi Xiangan, We recommend using the chimerax-users(a)cgl.ucsf.edu address for questions (CC'd on this message, so you don't need to do anything this time). The address you used is just for commercial licensing. Actually, if you get a traceback like this, you can usually use ChimeraX menu: Help... Report a Bug instead of sending e-mail. You may need to attach your session file so that we can reproduce the problem. I don't know the answer, so I'm using this address to share your question with the others. Best, Elaine ----- Elaine C. Meng, Ph.D. UCSF Chimera(X) team Department of Pharmaceutical Chemistry University of California, San Francisco

3 3