- Chimera-users - RBVI Mailing Lists

Automate calculation and saving of residues' structural attributes
by André 21 Jan '22

21 Jan '22

Hello! I'm working on an automated pipeline (in Python) for the structural analysis of proteins. In particular, I'm using Chimera to calculate and save to .txt files some structural properties of residues (namely: areaSAS, areaSES, relSESA, bfactor, kdhydro, phi, psi). These files are the inputs of the pipeline. However, currently these files are exported manually, and, for the purpose of building a fully automated pipeline, it would be great to also automate this step of generating the attributes and exporting their respective .txt files. I searched a lot, and tried some ways to achieve such an automation, but nothing worked (if you'd like to know about these (unsuccessful) ideas that I had, I can share them in next messages). Given that, I'd like to know if anyone has done something like this, and/or have any guidance on how I could proceed, if that's something possible at all to be done with Chimera (or even ChimeraX). I really appreciate any help you can provide! Thank you so much, André.

2 2

Question: How ro show the H-bond
by Kenji MATSUI 16 Jan '22

16 Jan '22

Dear Chimera company Thank you for always answering my questions so carefully. Thank you for allowing me to use your software. I am Japanese, so I hope you will understand that there are some inappropriate expressions. I would like to ask you one question. I entered the data as per the manual. However, the hydrogen bond did not show up. What should I do? The figure of Parameters [image: image.png] The result of showing [image: osh4,ID 223-230 .png] Details on how to do it *Display interactions between ligand and protein.* 1.Favorites→Model Panel→group/ungroup( Go down at close.) [image: image.png] 2.Favorite→Command line: combine #0.1,0.2 →select combination at model Panel [image: image.png] 3.rename 4.Select →Residue→ LIG or command: select :LIG 5.Tool→Surface/Binding Analysis→ Find H Bond [image: image.png] Check three points 1."Only find H-bonds..." e.g. "with at least one end selected". 2. "If endpoint atom hidden, shown endpoint residue" or else it will not show all of the H-bonds it finds 3.Write information to reply log Apply→OK --------------------------------------------------------------------------- 松井健治東京農工大学大学院工学府産業技術専攻修士〒184-8588 東京都小金井市中町2-24-16 Mail: s214903z(a)st.go.tuat.ac.jp Kenji Matsui Graduate School of Tokyo University of Agriculture and Technology 2-24-16, Nakamachi, Koganei-shi, Tokyo 184-8588, Japan Mail: s214903z(a)st.go.tuat.ac.jp

2 4

zoom on the selection after orientation of the selected object
by Enrico Martinez 16 Jan '22

16 Jan '22

Dear Chimera users! I am using ChimeraX for the visualization of protein-ligand interactions comparing two different pdb structures A - ensemble for doclong; B- X-ray pose: open A.pdb open B.pdb' match #2 to #1.1 # close all docking poses with the exemption of the first pose: close #1.2-end view wait # hide the protein from X-ray pose keeping only ligand for visualisation: hide #2 & protein cartoon hide #2 & protein surface hide #2 & protein #orient view on the model 1.1 and zoom closer view #1.1 orient; turn z 90 zoom 1.5 wait The problem is that when I use zoom 1.5 after view #1.1 it hides completely the representation for the molecule on the current GUI. Otherwise, it works fine if I center view on the both A and B: view orient; turn z 90 zoom 1.5 How could I zoom the screen properly after orientation of the view on the first molecule? Many thanks in advance! Enrico

1 0

Re: [Chimera-users] ChimeraX in no-gui mode: visualisation of protein-ligand interactions
by Enrico Martinez 15 Jan '22

15 Jan '22

Thank you Elaine! I played with it today and it looks great! Some additional questions (sorry for it at the first time the thing seems abit tricky): 1)may I automatically place the 2D label on the top left corner of the screen instead of manually defining its X-Y using xpos ypos? 2dlab create title text \"name_of_the_complex" color red font gothic xpos .04 ypos .99 size 40 2) assuming that I open 2 molecules in the same session open A.pdb open B.pdb then I remove all the models with the exemption of the first close #1.2-end How could I superimpose A to B (align the structures) and then hide cartoon from B.pdb (keeping only ligand) ? 3) Is there any possibility to put your signature on the screen-shot? like "Observed by Enrico Martinez" ? :-) Many thanks in advance!! Enrico сб, 15 янв. 2022 г. в 00:11, Elaine Meng <meng(a)cgl.ucsf.edu>: > > See the "color" command for coloring the ligand, and "set" or "graphics" for coloring the background, and "label" (or "2dlabels" but that usually requires you to position them interactively) for labeling. All of these are included in the "making images" page I mentioned near the end of my previous message. > Best, > Elaine > > > > > On Jan 14, 2022, at 3:03 PM, Enrico Martinez <jmsstarlight(a)gmail.com> wrote: > > > > Thank you very much, Elaine! > > I really like ChimeraX and it looks so great for me with all of those > > presents :-) I will look for the details regarding the options. Just a > > few addtional questions regarding customizing representation of > > protein-ligand interactions using ghost preset: > > In particular: I would like to color ligand (in orange??) > > To add labels (like Try-193) on the displayed side-chains around the ligands > > To change background from black to indigo! > > :-) > > Many thanks in advance! > > Cheers > > Enrico > > > > пт, 14 янв. 2022 г. в 23:48, Elaine Meng <meng(a)cgl.ucsf.edu>: > >> > >> Hello, > >> It isn't my area of expertise, but I don't think you can save images in nogui mode except with Linux OSMesa. > >> > >> As for your other questions, which commands do certain things... > >>> may you suggest some commands that I may introduce into my batch file to > >>> 1) delete all models from the multi-model pdb file with the exemption > >>> of the first model. > >>> 2) switch from cartoon to surface representation of the protein > >>> 3) visualize protein-ligand non-covalent interactions in the current model > >> > >> > >> (1) close submodels except for first one. Maybe something like command: > >> > >> close #1.2-end > >> > >> <https://rbvi.ucsf.edu/chimerax/docs/user/commands/close.html> > >> > >> (2) switch from ribbon to molecular surface of protein. Well, you probably want to do some coloring and other stuff besides just hiding ribbon and showing surface, which could be something like commands: > >> > >> ~ribbon > >> surface protein > >> > >> <https://rbvi.ucsf.edu/chimerax/docs/user/commands/cartoon.html> > >> <https://rbvi.ucsf.edu/chimerax/docs/user/commands/surface.html> > >> > >> Or you could try a surface-showing preset, e.g. > >> > >> preset ghost > >> > >> <https://rbvi.ucsf.edu/chimerax/docs/user/commands/preset.html> > >> > >> (3) depends what you mean by showing interactions (lines? selection? coloring?). You should take a look at commands "contacts" and/or "hbonds" and choose the options that give what you want. > >> <https://rbvi.ucsf.edu/chimerax/docs/user/commands/clashes.html> > >> <https://rbvi.ucsf.edu/chimerax/docs/user/commands/hbonds.html> > >> > >> You will need to spend some time trying the commands since there are an infinite number of possibilities for combinations of coloring and display styles and lighting and graphics effects. See also commands "lighting" and "graphics" and "set" > >> > >> <https://rbvi.ucsf.edu/chimerax/docs/user/commands/lighting.html> > >> <https://rbvi.ucsf.edu/chimerax/docs/user/commands/graphics.html> > >> <https://rbvi.ucsf.edu/chimerax/docs/user/commands/set.html> > >> > >> ... and the "making images" page with several relevant commands > >> <https://rbvi.ucsf.edu/chimerax/docs/user/images.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 Jan 14, 2022, at 2:17 PM, Enrico Martinez via Chimera-users <chimera-users(a)cgl.ucsf.edu> wrote: > >>> > >>> Dear Chimera's users! > >>> First I wish Happy New Year to everyone! > >>> I am just switching from Chimera to ChimeraX that I am going to use in > >>> no-gui (batch) mode to prepare screen shots of protein-ligand > >>> interactions based on the protein-ligand docking poses. Basically I am > >>> loading a pdb file of the complex (consisted of protein and ligand in > >>> multi-model format) from the terminal of my macOSX. I have a template > >>> for a batch file created in the bash script on the fly, that I used 2 > >>> years ago with an early version of ChimeraX. > >>> > >>> # a batch file for chimeraX > >>> printf "open ${pdb} > >>> view; zoom 1.0; turn x 90; turn z -80; > >>> wait > >>> # options for appearance > >>> preset 'overall look' publication > >>> # set bg powder blue > >>> set silhouettes t; set silhouetteColor indigo; set silhouetteDepthJump 0.02 > >>> > >>> # make info of the screen > >>> 2dlab create title text \"${pdb_name}\" color red font gothic xpos .04 > >>> ypos .92 size 35 > >>> > >>> # options for light > >>> light depthcuestart .4 depthcueend .8; lighting full intensity 1.4 > >>> shadows true qualityOfShadows finer > >>> > >>> save image ${structure}/${pdb_name}.png format png width 800 height > >>> 600 supersample 4 transparentBackground false > >>> wait > >>> exit" > chimeraX.${pdb_name}.cxc" > >>> > >>> then I use it directly in terminal as > >>> chimerax-daily chimeraX.${pdb_name}.cxc" > >>> > >>> may you suggest some commands that I may introduce into my batch file to > >>> 1) delete all models from the multi-model pdb file with the exemption > >>> of the first model. > >>> 2) switch from cartoon to surface representation of the protein > >>> 3) visualize protein-ligand non-covalent interactions in the current model > >>> > >>> Yours with thanks > >>> Enrico > >> > > >

1 0

ChimeraX in no-gui mode: visualisation of protein-ligand interactions
by Enrico Martinez 15 Jan '22

15 Jan '22

Dear Chimera's users! First I wish Happy New Year to everyone! I am just switching from Chimera to ChimeraX that I am going to use in no-gui (batch) mode to prepare screen shots of protein-ligand interactions based on the protein-ligand docking poses. Basically I am loading a pdb file of the complex (consisted of protein and ligand in multi-model format) from the terminal of my macOSX. I have a template for a batch file created in the bash script on the fly, that I used 2 years ago with an early version of ChimeraX. # a batch file for chimeraX printf "open ${pdb} view; zoom 1.0; turn x 90; turn z -80; wait # options for appearance preset 'overall look' publication # set bg powder blue set silhouettes t; set silhouetteColor indigo; set silhouetteDepthJump 0.02 # make info of the screen 2dlab create title text \"${pdb_name}\" color red font gothic xpos .04 ypos .92 size 35 # options for light light depthcuestart .4 depthcueend .8; lighting full intensity 1.4 shadows true qualityOfShadows finer save image ${structure}/${pdb_name}.png format png width 800 height 600 supersample 4 transparentBackground false wait exit" > chimeraX.${pdb_name}.cxc" then I use it directly in terminal as chimerax-daily chimeraX.${pdb_name}.cxc" may you suggest some commands that I may introduce into my batch file to 1) delete all models from the multi-model pdb file with the exemption of the first model. 2) switch from cartoon to surface representation of the protein 3) visualize protein-ligand non-covalent interactions in the current model Yours with thanks Enrico

4 5

Re: [Chimera-users] Question on pdb to mrc
by Elaine Meng 14 Jan '22

14 Jan '22

> On Jan 13, 2022, at 6:54 PM, Guan, Lan <Lan.Guan(a)ttuhsc.edu> wrote: > > Dear Elaine, > > I try to generate a reference model from crystal structure pdb file for cryo-EM particle alignment. > > I can get the map by “molmap #1 6”; #1 = my pdb file id; 6 = resolution I want to get in the output model. > > After imported into cryosparc, I can not use it. Error “ Cubic volumes only”. My question is how to change this Molmap output into a cubic volumes in Chimera? > I would appreciate for your input. > > Lan Hi Lan, (CC'ing chimera-users(a)cgl.ucsf.edu which is the address recommended for asking questions since others may have better answers) I don't know what cryosparc means by "cubic volumes only". My only guess is that maybe you did not save the molmap map correctly. The molmap command should make a new model which is the map, and then you can save the map with the Volume Viewer tool's File menu, or with the "volume" command. Saving data: <https://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/savemodel.html> Volume Viewer menu (File... Save map as...) <https://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/volumeviewer/vol…> volume command saving options <https://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/volume.html#output> Or, maybe cryosparc does not read any of the map formats that Chimera can write (listed in the link above). I hope this helps, Elaine ----- Elaine C. Meng, Ph.D. UCSF Chimera(X) team Department of Pharmaceutical Chemistry University of California, San Francisco

4 3

Re: [Chimera-users] Protonating Histidine Residues using Chimera
by Elaine Meng 13 Jan '22

13 Jan '22

> On Jan 13, 2022, at 10:54 AM, Joshi, Keya <keyajoshi1996(a)ku.edu> wrote: > > Hello Prof. Meng, > Hope you are doing well and staying safe. My name is Keya Joshi. I am a second year PhD student in the department of Computational Biology at the University of Kansas. I am trying to use Chimera to protonate Histidine residue 816 in chain B. But I have no idea about how to do this using Chimera. > Can you please help me with this? > > With Regards > Keya Joshi Hi Keya, I hope you are well, too! (CC-ing the chimera-users(a)cgl.ucsf.edu address which is recommended for asking Chimera questions unless you are sending data that needs to remain private) In the structure that you sent, histidine 816 in chain B is already protonated (already has hydrogens). However, it is protonated to the neutral state. Maybe you meant you wanted one more hydrogen to make it the positively charged state. Two ideas: (1) You can just use the AddH tool (menu: Tools... Structure Editing... AddH) or "addh" command to add hydrogens to the whole protein. See help for information about the options: <https://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/addh/addh.html> <https://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/addh.html> On the PDB file that you sent, the default options of AddH (= simply using command "addh" without any other options) do add a hydrogen to 816 in chain B to put it in the positively charged state, but it also adds lots of other hydrogens including on carbons and on all of the other residues. I don't know if that is a problem for you, or whether it is an improvement on the original structure. Of course you can delete some hydrogens afterward, but it depends what you want. (2) To just add one hydrogen to that sidechain nitrogen, instead use Build Structure (menu: Tools... Structure Editing... Build Structure), and in that tool change from "Start Structure" to "Modify Structure". Select (Ctrl-click on) the nitrogen of His 816 that needs the hydrogen. Selection is shown with a green outline on that atom. Then in the tool, you want that selected atom to be N (stay the same element) but with 3 bonds and trigonal geometry as shown in the screenshot. Also as shown in the screenshot, choose the option to leave the residue name unchanged. Then click Apply to add that one hydrogen. ... after I click Apply it looks like this: about selection: <https://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/selection.html> Build Structure tool: <https://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/editing/editing.…> You can find Chimera features using Chimera menu: Help... Search Documentation, e.g. to search for "hydrogens" or "building" Our newer program ChimeraX also has an "addh" command. I hope this helps, Elaine ----- Elaine C. Meng, Ph.D. UCSF Chimera(X) team Department of Pharmaceutical Chemistry University of California, San Francisco

1 0

Autodock Vina plugin
by Enrico Martinez 12 Jan '22

12 Jan '22

Dear Chimera Users, First Happy New Year to you and your families! I am trying to perform protein-ligand docking using UCSF Chimera 1.16v installed on my iMAC Big Sur. Using the pdb of the receptor and mol2 of the ligand (that I have already tested with VINA 1.2.3!) I got the following error: Model 1 (receptor.pdb) appears to be a protein without secondary structure assignments. Automatically computing assignments using 'ksdssp' and parameter values: energy cutoff -0.5 minimum helix length 3 minimum strand length 3 Use command 'help ksdssp' for more information. No SEQRES records for receptor.pdb (#1) chain A; guessing terminii instead Chain-initial residues that are actual N terminii: #1 SER 1.A Chain-initial residues that are not actual N terminii: Chain-final residues that are actual C terminii: Chain-final residues that are not actual C terminii: #1 THR 304.A 278 hydrogen bonds Removing spurious proton from 'C' of #1 THR 304.A Hydrogens added Wrote test.receptor.pdb into /Users/gleb/Desktop/docking_tuto adding gasteiger charges to peptide Traceback (most recent call last): File "/Applications/Chimera.app/Contents/Resources/share/vina/ws.py", line 395, in prepareReceptor execfile(scriptPath, d) File "/Applications/Chimera.app/Contents/Resources/lib/python2.7/site-packages/AutoDockTools/Utilities24/prepare_receptor4.py", line 172, in <module> delete_single_nonstd_residues=delete_single_nonstd_residues) File "/Applications/Chimera.app/Contents/Resources/lib/python2.7/site-packages/AutoDockTools/MoleculePreparation.py", line 533, in __init__ version=version, delete_single_nonstd_residues=delete_single_nonstd_residues) File "/Applications/Chimera.app/Contents/Resources/lib/python2.7/site-packages/AutoDockTools/MoleculePreparation.py", line 140, in __init__ self.addCharges(mol, charges_to_add) File "/Applications/Chimera.app/Contents/Resources/lib/python2.7/site-packages/AutoDockTools/MoleculePreparation.py", line 226, in addCharges chargeCalculator.addCharges(mol.allAtoms) File "/Applications/Chimera.app/Contents/Resources/lib/python2.7/site-packages/MolKit/chargeCalculator.py", line 80, in addCharges babel.assignHybridization(atoms) File "/Applications/Chimera.app/Contents/Resources/lib/python2.7/site-packages/PyBabel/atomTypes.py", line 127, in assignHybridization a.babel_atomic_number = self.get_atomic_number(a.babel_type) File "/Applications/Chimera.app/Contents/Resources/lib/python2.7/site-packages/PyBabel/atomTypes.py", line 103, in get_atomic_number (name,_name) ) ValueError: Could not find atomic number for Hn Hn Receptor preparation for AutoDock Vina failed; please look in Reply Log to see error messages.cannot prepare receptor for AutoDock Vina; please look in Reply Log and/or run Chimera with --debug flag to see errors. So here are my questions: 1) May this error be related to the version of OSX installed on my imac (I can not run MGL tools on the Big Sur for instance) ?? 2) Alternatively may I use the ADT scripts integrated in the Chimera like prepare_receptor4.py for my scripting purposes (outside of the UCSF Chimera) on the same iMAC machine ? Many thanks in advance! Cheers Enrico

6 13

Error on installation on a student's computer
by Kenward Vaughan 11 Jan '22

11 Jan '22

Hello all, As usual I have my students installing Chimera as part of their continuing adventure in organic chemistry, but have one of them reporting an error on attempted startup after installation. So far I have little information (we haven't started classes yet--next week--but she has sent the error message obtained. It is apparently on Windows. I suggested that she try starting from the command line to do the debug run, but this is from a distance at this time and I don't know how savvy she is with such stuff. Is there enough information in just the number noted by this message for someone to know what might be happening with her setup? Many thanks, and cheers! Kenward -- In a completely rational society, the best of us would aspire to be _teachers_ and the rest of us would have to settle for something less, because passing civilization along from one generation to the next ought to be the highest honor and the highest responsibility anyone could have. - Lee Iacocca Any fool can know. The point is to understand. - Albert Einstein

2 3

Problem with the app
by Nicola De Trino 11 Jan '22

11 Jan '22

I have download the chimera-1.16-win64.exe <https://www.cgl.ucsf.edu/chimera/cgi-bin/secure/chimera-get.py?file=win64/c…> for Windows 64 bit. I have installed but when i try to launch the program, it appears an error message : *chimera exited abnormally; the exit code was: -1073731819. * I try to install again the program but same result. I don't know how to fix the problem. <http://www.avg.com/email-signature?utm_medium=email&utm_source=link&utm_cam…> Mail priva di virus. www.avg.com <http://www.avg.com/email-signature?utm_medium=email&utm_source=link&utm_cam…> <#DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2> -- ------------------------ Indirizzo istituzionale di posta elettronica degli studenti e dei laureati dell'Università degli Studi di TorinoOfficial University of Turin email address for students and graduates

2 2