September 2021 - ChimeraX-users - RBVI Mailing Lists

outline the interface region in chimeraX
by Zhao, Bingchun (NIH/VRC) [F] 09 Sep '21

09 Sep '21

Hi, I wonder if there is a way to outline the interface region in ChimeraX like below: [cid:image001.png@01D7A597.5344EDF0] Previously, I did this in photoshop, but can we label the interface regions with those lines in ChimeraX? Thank you very much! Bingchun Bingchun Zhao, PhD Postdoctoral, VF Biodefense Research Section VRC/NIAID/NIH Bldg 40, Room 2607B 40 Convent Drive Bethesda, MD 20892 Tel: 301-761-7649 Email: bingchun.zhao(a)nih.gov<mailto:bingchun.zhao@nih.gov>

3 6

Re: [chimerax-users] getting AlphaFold predictions
by Elaine Meng 09 Sep '21

09 Sep '21

Hi Yk, Forwarding to chimerax-users(a)cgl.ucsf.edu since I don't have a good answer for this currently. (We prefer you use that address for questions unless they include some data that needs to stay private.) Thanks! Elaine > On Sep 9, 2021, at 1:57 PM, Yunsik Kang <kangy(a)ohsu.edu> wrote: > > > Hi Elaine, > > Thank you again! > > I did open it and using ChimeraX I can see the pdb files and see the structure! It is amazing! > > One question. > > Since my protein is massive and this is just parts of my protein, is there a way ChimeraX can stitch it all together? > > Like Tristan says it is 1400 residues each, can I perhaps somehow stitch the whole thing together? > > Yk > > -----Original Message----- > From: Elaine Meng <meng(a)cgl.ucsf.edu> > Sent: Thursday, September 9, 2021 8:12 AM > To: Yunsik Kang <kangy(a)ohsu.edu> > Cc: meng Elaine <meng(a)cgl.ucsf.edu> > Subject: [EXTERNAL] Re: [chimerax-users] getting AlphaFold predictions > > Hi Yunsik, > I do not know the answer to your questions. We are not the developers of AlphaFold or the AlphaFold database. We only added tools in ChimeraX to connect to them. > > Maybe there is information at the AlphaFold database about what's in their download. I have never downloaded these myself. > <https://www.alphafold.ebi.ac.uk/download> > > All I know is that "tar" is a way of putting multiple files together and that you have to "untar" to get them apart. This is a standard computer file thing, not related to any specific kind of files. > > Usually if you doubleclick the tar file your computer will automatically untar it and make some new folder or folders containing the files. You would need to have enough space available on your computer for the large number of files. > > I don't know how they are named but maybe there is a README file that comes with it, or if not, you may need to look at the files with a text editor. Maybe they have the Uniprot IDs of the proteins somewhere in the names or header text of the file. > > Elaine > ----- > Elaine C. Meng, Ph.D. > UCSF Chimera(X) team > Department of Pharmaceutical Chemistry > University of California, San Francisco > > >> On Sep 8, 2021, at 7:36 PM, Yunsik Kang <kangy(a)ohsu.edu> wrote: >> >> Hi Elaine, >> >> Thank you for the fast response. My proteins are not found in the Alphafold search. >> >> I am new to this so bear with me. I downloaded the human TAR file as Tristan has said, but have no idea how to search or even start looking for my protein. >> >> Could you help? Or is this something I need professional help with? >> >> Would it help to tell you the protein I am interested in? >> >> Thank you again. >> Yunsik >> >> From: Elaine Meng <meng(a)cgl.ucsf.edu> >> Sent: Wednesday, September 8, 2021 3:44 PM >> To: Yunsik Kang <kangy(a)ohsu.edu> >> Subject: [EXTERNAL] Fwd: [chimerax-users] getting AlphaFold >> predictions >> >> Forwarding in case you did not get Tristan's reply below: >> >> >> Begin forwarded message: >> >> From: Tristan Croll via ChimeraX-users <chimerax-users(a)cgl.ucsf.edu> >> Subject: Re: [chimerax-users] getting AlphaFold predictions >> Date: September 8, 2021 at 3:17:28 PM PDT >> To: ChimeraX Users Help <ChimeraX-users(a)cgl.ucsf.edu> >> Reply-To: Tristan Croll <tic20(a)cam.ac.uk> >> >> The human protein will be in the EBI-Alphafold2 database, but not accessible to the search tool. For human proteins (and only human, so far) longer than ~2000 residues, the predictions are provided as multiple 1400-residue fragments starting at residue 1, 201, etc.. The only way to get those right now is to download the tarball of the complete set of human protein predictions, and find them in there. >> >> >> On 8 Sep 2021, at 23:08, Elaine Meng via ChimeraX-users <chimerax-users(a)cgl.ucsf.edu> wrote: >> >> Hi Yunsik, >> I don't know where/how you can run a prediction for such a long sequence, maybe others can comment on that. However: >> >> First make sure there is not already a predicted structure for that protein in the AlphaFold Database! E.g. try the Search button on the AlphaFold tool and see if the same UniProt entries as you are interested in are found. The database includes human and other well-studied species (maybe also the yeast you mentioned, depending on the exact species). If the prediction is already in the database, there is no need to run a new calculation. >> >> ChimeraX AlphaFold tool: >> <https://rbvi.ucsf.edu/chimerax/docs/user/tools/alphafold.html> >> >> I hope this helps, >> Elaine >> ----- >> Elaine C. Meng, Ph.D. >> UCSF Chimera(X) team >> Department of Pharmaceutical Chemistry University of California, San >> Francisco >> >> >> On Sep 8, 2021, at 2:45 PM, Yunsik Kang via ChimeraX-users <chimerax-users(a)cgl.ucsf.edu> wrote: >> >> Hello, >> >> My name is Yunsik Kang, and I am a postdoc in Marc Freeman’s lab at the Vollum Institute. >> >> I would love to use ChimeraX to predict the structure of my protein of interest. I watched all the YouTube videos and tried to run the program. >> >> Unfortunately, my protein is 5005 amino acids in humans and 2958 aa in yeast. I get a message “Please use the full AlphaFold system for long sequences.” >> >> My question is what is the best way to approach this problem? Should I cut the protein in half and run the program? In one of the videos, it mentions after 700 aa it will have problems. Will it work if I get Colab-Pro? Or would the server crash no matter what. >> >> I am not a structural biologist, but I hope the structure can help be predict me with my research. >> >> Thank you, >> Yunsik >> >> >> _______________________________________________ >> ChimeraX-users mailing list >> ChimeraX-users(a)cgl.ucsf.edu >> Manage subscription: >> https://www.rbvi.ucsf.edu/mailman/listinfo/chimerax-users >> >> _______________________________________________ >> ChimeraX-users mailing list >> ChimeraX-users(a)cgl.ucsf.edu >> Manage subscription: >> https://www.rbvi.ucsf.edu/mailman/listinfo/chimerax-users >> > >

1 0

Save pdb file
by Zhao, Bingchun (NIH/VRC) [F] 09 Sep '21

09 Sep '21

Hi, First, I wanna thank you for all your work to develop ChimeraX, it’s so powerful. I do have a question and would like to have your input. I opened a EM density map (.mrc) and two pdb files in chimerax (version 1.2), then fit all the pdb models into the EM map. Then, I hide the EM map, only show the pdb models after fitting, and then I would like to save these all pdb files in a new pdb file with all position information after fitting. Would you please give me any direction how to do this? I tried saving them, but only partial pdb files saved, also tried combine commend, but it said no such commend. Thank you very much! Best, bingchun Bingchun Zhao, PhD Postdoctoral, VF Biodefense Research Section VRC/NIAID/NIH Bldg 40, Room 2607B 40 Convent Drive Bethesda, MD 20892 Tel: 301-761-7649 Email: bingchun.zhao(a)nih.gov<mailto:bingchun.zhao@nih.gov>

2 2

error saving a surface
by Diego Amaya 08 Sep '21

08 Sep '21

Hi, I'm running a python script in --nogui and --script mode and I'm getting an error when doing rc(session, 'save surface.stl format stl'). The error message is : *QGuiApplication::font(): no QGuiApplication instance and no application font set.Segmentation fault* I'm working with debian 10 and I used the production release 1.2.5-1 Thanks in advance for the help. -- Diego A. Amaya Ramírez -- *Aviso legal:* El contenido de este mensaje y los archivos adjuntos son confidenciales y de uso exclusivo de la Universidad Nacional de Colombia. Se encuentran dirigidos sólo para el uso del destinatario al cual van enviados. La reproducción, lectura y/o copia se encuentran prohibidas a cualquier persona diferente a este y puede ser ilegal. Si usted lo ha recibido por error, infórmenos y elimínelo de su correo. Los Datos Personales serán tratados conforme a la Ley 1581 de 2012 y a nuestra Política de Datos Personales que podrá consultar en la página web www.unal.edu.co <http://www.unal.edu.co/>.* *Las opiniones, informaciones, conclusiones y cualquier otro tipo de dato contenido en este correo electrónico, no relacionados con la actividad de la Universidad Nacional de Colombia, se entenderá como personales y de ninguna manera son avaladas por la Universidad.

2 2

capture reply from run(session, command)
by Sean Connell 08 Sep '21

08 Sep '21

Hello I was wondering if it is possible to capture the output that goes to the log when you run the python command run(session, command). For example if I run: run(session, "mmaker #2/A to #1/A”) Can I store in a variable the reply in the log file that has information about the RSMD? Sean

2 2

distance between centers of mass using python API
by Diego Amaya 08 Sep '21

08 Sep '21

Hello, I'm trying to compute the distance between the centers of mass of two residues using the python API. I tried two different ways : 1. *rc(session, 'define centroid :43 massWeighting true name AA1')rc(session, 'define centroid :35 massWeighting true name AA2')rc(session, 'distance AA1 AA2') * 2. *rc(session, ' measure center :43 mark true name AA1')rc(session, 'measure center :35 mark true name AA2')rc(session, 'distance AA1 AA2')* But I get this error in the third line in both cases : *chimerax.core.errors.UserError: Missing or invalid "objects" argument: invalid objects specifier * What is wrong ? Thanks in advance. -- Diego A. Amaya Ramírez -- *Aviso legal:* El contenido de este mensaje y los archivos adjuntos son confidenciales y de uso exclusivo de la Universidad Nacional de Colombia. Se encuentran dirigidos sólo para el uso del destinatario al cual van enviados. La reproducción, lectura y/o copia se encuentran prohibidas a cualquier persona diferente a este y puede ser ilegal. Si usted lo ha recibido por error, infórmenos y elimínelo de su correo. Los Datos Personales serán tratados conforme a la Ley 1581 de 2012 y a nuestra Política de Datos Personales que podrá consultar en la página web www.unal.edu.co <http://www.unal.edu.co/>.* *Las opiniones, informaciones, conclusiones y cualquier otro tipo de dato contenido en este correo electrónico, no relacionados con la actividad de la Universidad Nacional de Colombia, se entenderá como personales y de ninguna manera son avaladas por la Universidad.

2 5

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

ChimeraX-users September 2021