- Chimera-users - RBVI Mailing Lists

Max 100 residues for local bug?
by Danilo Boškovic 16 Sep '19

16 Sep '19

Dear caretakers, I tried to do autodock vina but as soon as i did i got an error that i can only do max 100 residues. I use a local installation of AVina. Is this a bug or is there somewhere something im missing? 41965 build number 1.13.1 version Sincerely, D. B.

2 1

Where is the radius.nc map file in the Semliki Forest virus?
by sunyeping 16 Sep '19

16 Sep '19

Dear everyone, The Semliki Forest virus tutorial (https://www.cgl.ucsf.edu/chimera/data/tutorials/semliki/semliki.html) mentioned a radius.nc which is used for visualization of the layer. However, the tutorial doesn't describe where the reader can download it. Do you know where to get thid file? Thank you in advance. Yeping Sun

2 1

Re: [Chimera-users] Electrostatic potential molecular surfaces
by Elaine Meng 11 Sep '19

11 Sep '19

Hi Fernando, Our resources don’t allow spending much time on consulting with individual research projects, sorry, and I can only answer mostly as it pertains to Chimera as I don’t know the details of parameters developed by other people or other programs. Also Chimera questions should be sent to this mailing list instead of to me directly so that other people can also benefit, unless private data are discussed or attached. You found an old script that is really special-purpose, finding the surface area that is more red than blue after coloring by electrostatic potential. Note that this would include surface that looks white but has a tiny, tiny tint of pinkishness. Also, the surface area that is more blue than red is just the complement of the surface that is more red than blue. In other words, if the surface is 75% red>blue then it is 25% blue>red because very little if any will be pure white. They always add up to the total surface. So no, you cannot use these scripts to figure out how much of the surface is polar or nonpolar. It may not matter given the above, but which method, set of charges, and what parameter values to use for electrostatic potential (ESP) calculations are research decisions that I cannot make for you. There are always tradeoffs, and different people have different opinions, etc. I don’t know python, but your generalization to coloring by hydrophobicity by identifying cornflower blue vs.goldenrod might be wrong. Red vs. blue is simpler since they are RGB color components. However, even if the color comparison is done correctly, the logical problem is the same, they always add up to the whole surface, and nearly white areas will always be combined with one side or the other. However, you can get the total surface area given some range of hydrophobicity values independent of any coloring. Hydrophobicity in this case is just a constant lookup value for each type of amino acid, according to this table (kdHydrophobicity column): <http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/hydrophob.html> You would have to decide yourself which range of values is polar or nonpolar, but then you could just select all the residues that are less than or more than a certain value, e.g. command: select :/kdHydrophobicity>0.0 Alternatively you could just use command: select polar select ~polar …but you might not agree with what amino acids are considered polar, as shown in a table in this page: <http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/resprop/resprop.h…> Then if surface is shown you could add up the area for the selected residues, e.g. menu: Tools… Structure Analysis… Attribute Calculator, and in that tool calculate attribute named whatever you want for “molecules” with Formula: sum(areaSAS.residue) (… or SES instead of SAS, to be compared with the corresponding total area given in the Reply Log when you show the surface) with option to “Restrict formula domain to current selection” and “Show calculation results in Reply Log" turned on, other options turned off, click Apply. <http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/calculator/calcul…> Totally separate and not related to Chimera, but I happen to know of this web server you might be interested in, for comparing ESP of related proteins: <https://pipsa.h-its.org/pipsa/> Best regards, Elaine ----- Elaine C. Meng, Ph.D. UCSF Chimera(X) team Department of Pharmaceutical Chemistry University of California, San Francisco > On Sep 11, 2019, at 2:16 PM, Fernando Villa <fer.vdl1928(a)gmail.com> wrote: > > > Dear Elaine C. Meng, Ph.D > > I have some questions about electrostatic potential molecular surfaces and protein comparison. > > 1. The following scripts are added, edited from the redarea.py script of link http://plato.cgl.ucsf.edu/trac/chimera/wiki/Scripts, considering that the script allows reading the surface in Å2 colored more red than blue area. > > 1. Add scripts: > > Command: > > > open > > -blue_area_03.py > > -cornflowerblue_area_03.py > > -goldenrod_area_02.py > > -red_area_02.py > > (the scripts are attached) > > > 2. I used the following commands: > > > open 4MZ2 > > > delete solvent > > > addh > > > addcharge > > > Surface > > coulombic -10 red 0 white 10 blue > > > The electrostatic potential surface is presented by Coulombing coloring Surface: > > > > <image.png> > I execute the commands: > > > redarea #0 > > red area 3562, total area 5730, ratio 0.6217 > > > bluearea #0 > > blue area 2168, total area 5730, ratio 0.3783 > > > The dominant areas of both colors are obtained. > > The following command is input: > > > rangecolor kdHydrophobicity min cornflowerblue 0 white max goldenrod > > > > <image.png> > Next: > > > cornflowerbluearea # 0 > > cornflowerblue area 4109, total area 5730, ratio 0.7171 > > > > goldenrodarea # 0 > > Goldenrod area 1621, total area 5730, ratio 0.2829 > > > So, I believe (I'm not sure): > > Blue area + Red area = PSA ß Is correct? > > White area = NPSA ß Is correct? > > > So, the real blue and red areas would be (not sure): > > > (Redarea * PSA) / Total surf = Real Red area > > (Bluearea * PSA) / Total surf = Real Blue area > > So: > > Bluearea = 1554.6792 > > Redarea = 2554.3208 > > > Is this calculation correct? > > > Could I apply these calculations to compare them with other proteins of the same family? > > For example: > > Protein members of the immunoglobulin family. > > > > In the Chimera version 1.14, I Download the following programs pdb2pqr-windows-bin64-2.1.1 > > apbs1.5_win64 > > Downloaded from the website: http://www.poissonboltzmann.org > > Next, I assign charges and radii with the following parameters of PDB2PQR > > > <image.png> > And I saved the file. Pqr > > > <image.png> > Later I opened the file 001_B4.pqr with Chimera 1.14 > > > <image.png> > I use the following command: > > > surface > > And run the APBS interface with the following default parameters: > > <image.png> > > As a result, the color gradient of the molecular surface appeared as shown: > > > <image.png> > <image.png> > > > However, my question is: > > Are these parameters correct using the force field PARSE of PDB2PQR? Or, is it necessary to use a different force field? > > > Could I apply the scripts: > > -blue_area_03.py > > -cornflowerblue_area_03.py > > -goldenrod_area_02.py > > -red_area_02.py > > > to calculate the redarea and bluearea and compare it with other proteins? > > > Is there a script that gives me the correct areas of both red and blue colors? > > > Sorry for my long email, as I am new to this area. And I really need help. > > I apologize for so many questions. > > I would greatly appreciate your help. > > > Best regards, > > Fer. > > > > > <blue_area_03.py><red_area_02.py><goldenrod_area_02.py><cornflowerblue_area_03.py>

2 3

How to analyze cryo-EM reconstruction map of virus and prepare publishable figures?
by sunyeping 11 Sep '19

11 Sep '19

Dear all, I am new to cryo-EM structural biology. Now I have a cryo-EM reconstruction map for a virus particle. I wish to analyze the characteristics of the virus, such as the h, k and t number, the component of an cryo-EM asymmetric unit of the capid and prepare some publishable figures. I have no idea how to do these things. Could all these be done with Chimera? Is there any tutorials or manual for these calculation, visualization and drawing? If you could recommend some ways by which I can learn to do these, I would be greatly grateful. Best regards

2 1

APBS electrostatic potential map not visible
by Ahmad Khalifa 10 Sep '19

10 Sep '19

Hello, After PDB2PQR and APBS calculation, I get an .dx that I can't view (doesn't appear when I click show), let alone color. I tried to launch volume viewer, it does display a surface. What am I doing wrong? Regards.

2 1

Modeller from a Python script
by Samo Lešnik 10 Sep '19

10 Sep '19

Dear Madam/Sir, I would like to run the Modeller missing loop function implemented in Chimera by using a Python script, as I want to go through a number of protein structures. Is this possible? I do not find the appropriate commands in the comand list. I would use a local instalation of Modeller. Thank you for your help. Best regards, Samo

2 1

DNA drawing inquiry
by Woorim Choi 09 Sep '19

09 Sep '19

2 1

Way to compute residue wise RMSD using python script in chimera:
by Kamalesh Damodaran 09 Sep '19

09 Sep '19

Dear sir/madam, I would like to superimpose the two structures and compute the residue wise RMSD values. With Gui support in chimera, I opened the two structures and used Matchmaker to superimpose two structures. In Matchmaker, I checked the box 'Show pairwise alignment(s)' and click Apply, subsequently it displays MultAlignViewer dialogue box, where --> Headers --> RMSD : backbone --> Save --> filenamegiven. This process would give me residue wise RMSD score only for backbone atoms in .hdr format. I wants to write python script to do the same, so that I will compare large number of structures, and my script is able to perform till the procedure of superimposition using matchmaker and display MultAlignViewer dialogue box. I don't know the way to store the RMSD : backbone and RMSD : calpha in .hdr format using python script through chimera. Can some one please help me in this. The reason why I specifically want to do in this approach is, because among the two structures, some residues of one structure is missing. Since this approach will compute RMSD only for the matching residues, I am interested in this. Hope my question is not confusing. Thanks in advance. Best regards, Kamalesh D.

2 3

cage builder
by Artur Biela 06 Sep '19

06 Sep '19

Hi all, I would like to build a cage model using cage builder, but the problem is that I need to use different polygons than available in the pop-up window of cage builder. My question is: is there any option to use for example octagons, decagons or even more-gons? Thank in advance Artur =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-= Artur BIELA, PhD Scientific Assistant Bionanoscience and Biochemistry Laboratory Malopolska Centre of Biotechnology Jagiellonian University 7A, Gronostajowa Street 30-387 Krakow Poland Department of Cell Biology and Imaging Institute of Zoology and Biomedical Research Jagiellonian University Gronostajowa 9 30-387 Krakow artur.biela(a)uj.edu.pl http://www.mcb.uj.edu.pl/en_GB/heddle-laboratory publons<https://publons.com/mashlets?el=badgeCont570&rid=M-3954-2019> =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=

3 3

How could I know the segment label
by Huiya Zhou 03 Sep '19

03 Sep '19

Hi, I use the “Segment Map” for the volume data to divide the density map into parts. How could I know the information of each region? For example, the label, the average value of each region. I also would like to know the regions which are connected with one specific region. How could I obtain information in python script? Thanks, Huiya

2 1

Support chimera
by Jose Alirio Mendoza Mesa 03 Sep '19

03 Sep '19

Hi there, my name is Jose Alirio Mendoza. Currently, I am carrying out my PHD thesis in which I am using the chimera program to show the tomography results. Nevertheless, I am being presented some problems when trying to measure my particle in Chimera. The file of my product is “mrc” and I have been investigating about chimera usage where I have realized that the file must be “pdb” for the program to recognize it. I do not really know how to convert my molecule to the required file, and I guess that is the problem by which I cannot measure the distance of my particle. When I click “tools – structure analysis – distance” it does not appear the code of my molecule as in picture below [image: imagen.png] This is my particle and what I get when using the tool for measuring. I have other three tomography for one publication. [image: imagen.png] There is not any code in the tool square and that is why I guess to no be able to measure the particle. Other tool I tried to use is the scale bar, however the value of it is not according to my calibration. [image: imagen.png] How may I change that calibration? In our research group we want the publish the articles citing chimera as visualization tool so we would really appreciate your help. Here is attached the file we have worked in chimera. pollo 13.mrc <https://drive.google.com/a/unal.edu.co/file/d/1cAs2g1kP88o8im65xuTbOz11o9pk…> -- *Qco. JOSE ALIRIO MENDOZA MESA* Candidato a Doctor en ciencias Universidad Nacional de Colombia. 3165000 Ext. 14475 *LABORATORIO DE CATALISIS HETEREOGENEA*

2 1

Re: [Chimera-users] code 256
by Greg Couch 30 Aug '19

30 Aug '19

Like I said before in http://plato.cgl.ucsf.edu/pipermail/chimera-users/2019-August/016113.html, you need to use the daily build of Chimera. Go the Chimera download web page, http://www.cgl.ucsf.edu/chimera/download.html. Scroll down past the 1.13.1 downloads, and pick the 64-bit Linux daily build version to download. Then install it. When 1.14 comes out, it will not have this problem. -- Greg On 8/29/2019 11:22 PM, Aleksandra Jarmołowicz wrote: > > I am using Ubuntu 19.04. So what can I do to use Chimera with this > version? Before update of the system programme runs well. > > Aleksandra Jarmołowicz > Ph.D. Student > > Department of RNA Structural Genomics > Institute of Bioorganic Chemistry > Polish Academy of Sciences > Z. Noskowskiego str. 12/14 > 61-704 Poznań > Phone: +48 61 852 8503 ext 143 > ajarmolowicz(a)ibch.poznan.pl <mailto:ajarmolowicz@ibch.poznan.pl>

1 0

Buried_area python function
by Kevin Wang 29 Aug '19

29 Aug '19

Hi, I am trying to use the buried_area function below. What are the arguments that I need to pass in? def buried_area(operation, atoms1, atoms2, probeRadius = 1.4, vertexDensity = 2.0): My current code is as follows: #calculate buried area import os import chimera import xlsxwriter from Measure import measure from chimera import * workbook = xlsxwriter.Workbook('contactArea.xlsx') worksheet = workbook.add_worksheet() Path = "/Users/kevinwang/desktop/VEGFR_NO_LIGAND/results/clustering/pdbfit/split" i = 0 filelist = os.listdir(Path) for x in filelist: if x.startswith("trp.pdb."): #if (x == "trp.pdb.00001"): opened = chimera.openModels.open(x, type ="PDB") chimera.update.checkForChanges() try: runCommand('surface :286-596') runCommand('surface :899-1230') y = measure.buried_area() worksheet.write(i,0,x) worksheet.write(i,1,y) except: worksheet.write(i,0,x) worksheet.write(i,1,"ERROR") i = i + 1 print("Processed: " + x) print("Total #: " + str(i)) runCommand('close all') chimera.update.checkForChanges() workbook.close() Best, Kevin

2 1

copy transformation between volumes
by Oleksiy Kovtun 29 Aug '19

29 Aug '19

Dear community, I need your advice on the following problem. I align map #1 to the reference volume using fit command, and then need to confer this transformation to map #2. This could be done by using matrixcopy command. However in this particular case #1 has twice larger voxel size than #2. It is understod, that since transformation matrix describes shifts in non-scaled pixels, I can not use matrixcopy command. Is it correct to scale translations (column 4) in the fit matrix for #1 and apply it to #2 with matrixset command? Thank you and best regards, Oleksiy Dr Oleksiy Kovtun Research Fellow Briggs group MRC Laboratory of Molecular Biology Francis Crick Avenue Cambridge CB2 0QH UK Tel: + 44 1223 267551 Email: okovtun(a)mrc-lmb.cam.ac.uk

2 1

Python Script to Find Contact Area Between 2 Proteins
by Kevin Wang 29 Aug '19

29 Aug '19

Hi, I am writing a python script to find the area between 2 proteins using the contact_area python command below. def contact_area(p1, p2, d, color = None, offset = None, slab = None, smooth = False, optimize = True): Is there a specific object for p1 and p2 and if so how do I get it? My current code is as follows: import os import chimera from Measure import contactarea from chimera import * Path = "/Users/kevinwang/desktop/VEGFR_NO_LIGAND/results/clustering/pdbfit/split" filelist = os.listdir(Path) for x in filelist: #if x.startswith("trp.pdb."): if (x == "trp.pdb.00001"): print x opened = chimera.openModels.open(x, type ="PDB") runCommand('surface') y = contactarea.contact_area('0.1', '0.2', 2.5) text_file = open("buriedArea.txt", "w") text_file.write(y) text_file.close() The error message I am getting is as follows: AttributeError: 'str' object has no attribute 'geometry' File "/Applications/Chimera.app/Contents/Resources/share/Measure/contactarea.py", line 8, in contact_area v1, t1 = p1.geometry How should I be passing in my arguments? Best, Kevin

2 3

code 256
by Aleksandra Jarmołowicz 29 Aug '19

29 Aug '19

Dear Sir or Madame, During chimera-1.13.1-linux_x86_64 installation process I have found some difficulties. I have not found any solutions on your website. After running instalation an error appears : ERROR in chimera_final_install: unable to install deskop icon: result code from installer: 256 I would be very greatful for any help. Kind regards, Aleksandra Jarmołowicz -- Aleksandra Jarmołowicz Ph.D. Student Department of RNA Structural Genomics Institute of Bioorganic Chemistry Polish Academy of Sciences Z. Noskowskiego str. 12/14 61-704 Poznań Phone: +48 61 852 8503 ext 143 ajarmolowicz(a)ibch.poznan.pl

2 1

Can multiple marker sets be batch-saved into individual .cmm files?
by Jun Ha Song 28 Aug '19

28 Aug '19

Dear RBVI, I need to save around 340 market sets into individual .cmm files, 1 .cmm for each marker set I have in my Chimera session. I have been looking at UCSF Chimera documentation for some time ( https://www.cgl.ucsf.edu/chimera/docs/ContributedSoftware/volumepathtracer/…) and I wanted to ask if multiple marker sets (each with its own unique model ID) can be saved in single batch command line. Is there a possibility of specifying commands like *save* or *write* or *export* with atom specification of all the model that I want to save to make Chimera save all my 340 marker sets as their individual .cmm files? I was hoping that I could simply use something like *save #1-340 individual UseModelIDasFilename format .cmm directory E:\LBNL\Individual_Models* to make it happen, but all the commands that I thought of as potential candidate for batch-save command doesn't mention whether it can work with volume tracer's "Save current marker set as...". or, if I can use Volume Tracer's "Save current marker set as..." using a script, I could handle this by making myself very simple macro, but I don't know whether volume tracer's codes would allow some quick script to substitute its "Save current marker set as..." function. So, in summary, I wanted to know if... Q1 : Is there a command in UCSF Chimera that allows batch saving multiple market sets as individual .cmm files? Q2 : Can Volume Tracer's "Save current marker set as..." function be used by using some kind of script? Sincerely, Junha Song.

2 1

Get map contour level with python script or command
by Kyle Morris 28 Aug '19

28 Aug '19

Dear chimera users and dev, Is there a chimera or python command that can report the current contour level of a volume that is in model number #n back to the reply log or equivalent? I want to do this after executing 'volume #n sdLevel 3’ such that I can programmatically retrieve the contour level associated with a particular sigma level of a currently loaded volume. Apologies if in my searching I have missed this on the boards, I think this is a novel question though. Thanks as always for your help. Best wishes, Kyle

3 2

How to read a volume in python script
by Huiya Zhou 28 Aug '19

28 Aug '19

Hello, I am trying to use UCSF Chimera to analysis cry-EM density maps. I would like to read a density map (for example: a .mrc file) in python script and get the array (or a three dimensional matrix) of all voxels of this density map. How could I realize this? By the way, I tried to learn some basic chimera functions which can be used in python on this website: https://www.cgl.ucsf.edu/chimera/docs/ProgrammersGuide/Examples/index.html <https://www.cgl.ucsf.edu/chimera/docs/ProgrammersGuide/Examples/index.html> However, it seems this website can’t be opened and it shows me a blank page. Best, Huiya

3 2

Chimera Error
by Ryan Harlich 27 Aug '19

27 Aug '19

Hi Chimera, I am having trouble opening a pdb file. I get the following error: ValueError: can’t normalize zero length array File \Chimera1.13.1\share\chimera\tkgui.py, line 3030, In showLabelBalloon Objs = self.viewer.pick(event.x, event.y) All other PDBs are opening so I am not sure why this one won’t open. Best, Ryan Harlich Sent from Mail<https://go.microsoft.com/fwlink/?LinkId=550986> for Windows 10

2 1

Subtracting the difference in EM volumes
by Y. Mutum 27 Aug '19

27 Aug '19

Hi Is there a way to subtract the non-overlapping bits and retain only the overlapping volumes between 2 EM density maps in Chimera? I have tried the vop subtract command, but it only seems to subtract the overlapping bits and not the inverse of it (non-overlapping bits)? Sorry, if I have missed it. Thanks Yilog

2 1

Divided the Volume Data into two parts
by Huiya Zhou 23 Aug '19

23 Aug '19

Hello, I would like to divided a density map (volume data) into two parts. For example, the size of the whole density map is 320*320*320. I want to get one part which covers the region from (i,j,k) to (x,y,z), e.g. from (0,0,0) to (100,100,100). The rest is the other part. The transparency for the first part is 0 and for the second part is 0.5. And I still want to get the histogram of the first part. How could I achieve this in a python script? Actually, I want to get the histogram of different regions. I plan to run a loop for different regions (Change the (i,i,k) and (x,y,z) for the first part). Could I save all of the histograms in one file which is not a window of the Chimera but a picture. Best, Huiya

2 1

Selecting bonds/pseudobonds with a cutt-off distance
by Yasser Almeida Hernández 23 Aug '19

23 Aug '19

Hi all, how can I select bonds/pseudobonds with a cutt-off distance, say bonds/pseudobonds with more than/less than X? Thanks in advance, Yasser

2 2

回复： scipy problem
by William 22 Aug '19

22 Aug '19

Dear Eric, I have tried the instructions following the FAQ, but there is another problem. The scipy installed successfully. However, it installed a higher version of the numpy when installing scipy. And the older version of numpy is still in. The chimera can not work. Anyway, thanks a lot. ------------------ 原始邮件 ------------------ 发件人: "Eric Pettersen"<pett(a)cgl.ucsf.edu>; 发送时间: 2019年8月20日(星期二) 凌晨1:46 收件人: "William"<3121437640(a)qq.com>; 抄送: "chimera-users"<chimera-users(a)cgl.ucsf.edu>; 主题: Re: [Chimera-users] scipy problem Hi William, I don’t know why Jaime’s instructions don’t work. Maybe he can follow up. Regardless, we have our own instructions for installing pip (and pip-able packages like scipy). It’s in the Chimera Programming FAQ. Question 3b. If those instructions don’t work, please let us know. —Eric Eric Pettersen UCSF Computer Graphics Lab

2 1

How to read a volume in python script
by Huiya Zhou 22 Aug '19

22 Aug '19

Hello, I am trying to use UCSF Chimera to analysis cry-EM density maps. I would like to read a density map (for example: a .mrc file) in python script and get the array (or a three dimensional matrix) of all voxels of this density map. How could I realize this? By the way, I tried to learn some basic chimera functions which can be used in python on this website: https://www.cgl.ucsf.edu/chimera/docs/ProgrammersGuide/Examples/index.html <https://www.cgl.ucsf.edu/chimera/docs/ProgrammersGuide/Examples/index.html> However, it seems this website can’t be opened and it shows me a blank page. Best, Huiya

2 1

copy/combine molecular models
by Rui Chen 22 Aug '19

22 Aug '19

Hello, I am working on a homodimer. I modelled loops of each chain one by one. I chose the best loops for each chain. Now I want to save them as a pub file. Please see the following screenshot. How to decide the new model's ID? I know the coordinate system is the reference coordinate, how do I decide which chain to use as the reference chain? Should I choose rename them uniquely? Looking forward to hearing from you. Thank you. [image: Screen Shot 2019-08-22 at 12.01.28 AM.png] Best regards, Rui

2 2

how to guess ligand's net charge?
by Thomas Evangelidis 21 Aug '19

21 Aug '19

Greetings, How does Chimera guess the net charge of a ligand when prompted to select the charge type (AM1-bcc or Gasteiger)? I am doing high throughput protein-ligand MD simulations and would like to know which open source library Chimera uses for this task. Thanks in advance. Thomas -- ====================================================================== Dr. Thomas Evangelidis Research Scientist IOCB - Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences <https://www.uochb.cz/web/structure/31.html?lang=en>, Prague, Czech Republic & CEITEC - Central European Institute of Technology <https://www.ceitec.eu/>, Brno, Czech Republic email: tevang3(a)gmail.com, Twitter: tevangelidis <https://twitter.com/tevangelidis>, LinkedIn: Thomas Evangelidis <https://www.linkedin.com/in/thomas-evangelidis-495b45125/> website: https://sites.google.com/site/thomasevangelidishomepage/

2 1

color by zone in chimeraX
by Pranav Shah 20 Aug '19

20 Aug '19

Hi Elaine and Co., I am trying to color a virus map with the relevant pdb file, the series of commands I am running are as follows, open virus.mrc open virus.pdb color /A plum target c color /B lightgreen target c color /C dodger blue target c color zone #1 near #2 dist 3 However, the resultant surface is rendered in splotches of red and blue, which leads me to think the it is coloring the surface using atom type. Could you tell me how I can get Chimera X to color the surface using the colors of the chains only? FYI, color by zone works flawlessly in chimera. Cheers, Pranav -- Pranav Shah Postdoctoral Research Fellow. Hogle Lab Department of Biological Chemistry and Molecular Pharmacology, Building C2 - First Floor 240 Longwood Avenue Boston, MA 02115 (617) 432-4360 (fax) (617) 432-3839 (lab)

2 2

scipy problem
by William 19 Aug '19

19 Aug '19

Hello, I tried to install scipy in chimera python following the steps below: When I tried the third step, it raised midaserror. The error is as follows： Is there anything wrong with my path? Thanks.

3 2

Match showmatrix true motion in original coordinates
by Ahmad Khalifa 13 Aug '19

13 Aug '19

Hi, Running match showmatrix true, gives two matrices. What's the difference between the two? I already saved the file after fitting it so it should have only one set of coordinates, its original coordinates?! Regards Motion of #76 in #2 coordinate system Matrix rotation and translation 0.99998889 -0.00169318 0.00440013 -1.57825672 0.00165962 0.99996960 0.00761849 -4.96792547 -0.00441290 -0.00761111 0.99996130 169.74075587 Axis -0.85023004 0.49200930 0.18717833 Axis point 0.00000000 21533.64775327 2718.66675667 Rotation angle (degrees) 0.51315730 Shift along axis 30.66940614 Motion of #76 original coordinates in #2 coordinate system Matrix rotation and translation 0.78916125 0.61416829 0.00467203 -143.02960322 -0.61417583 0.78916918 0.00023085 242.10381668 -0.00354524 -0.00305162 0.99998906 98.04867640 Axis -0.00267220 0.00668954 -0.99997405 Axis point 281.92931592 330.08522286 0.00000000 Rotation angle (degrees) 37.89294229 Shift along axis -96.04436567

2 1

Chimera bug report submission
by Tanara Arenhart 12 Aug '19

12 Aug '19

The following bug report has been submitted: Platform: Windows64 (WGL) Microsoft Windows Unknown 10.0 (Build 17134) x64 Chimera Version: 1.13.1_b41965 Description Hi! I'm trying to dock my protein with a ligand, but everytime appears this message: "receptor and ligand cannot be the same model". I don't know why, because i'm doing the same things I always do. What could be happening? Thank you! Open models #0 3dbs.pdb #1 smiles:CC1CCC2C3(CCC(... #2 AutoDock Vina Search Volume Reply log: Failed to translate SMILES to 3D structure via NCI web service (SMILES: receptor and ligand cannot be the same model) Failed to translate SMILES to 3D structure via Indiana University web service (SMILES: receptor and ligand cannot be the same model) Web services failed to translate SMILES string to 3D structure. Translated SMILES to 3D structure via NCI web service (SMILES: CC1CCC2C3(CCC(C(C3CCC2(C14CC5=C(O4)C(=CC(=O)C5=O)C)C)(C)C)O)C) No incomplete side chains Chain-initial residues that are actual N terminii: 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: 0 hydrogen bonds Hydrogens added Charge model: AMBER ff14SB Assigning partial charges to residue STP (net charge +0) with gasteiger method Running ANTECHAMBER command: C:/Program Files/Chimera 1.13.1/bin/amber16/bin\antechamber -i c:\users\nara_\appdata\local\temp\tmpi9opye\ante.in.mol2 -fi mol2 -o c:\users\nara_\appdata\local\temp\tmpi9opye\ante.out.mol2 -fo mol2 -c gas -nc 0 -j 5 -s 2 -dr n (STP) (STP) Welcome to antechamber 17.3: molecular input file processor. (STP) (STP) Info: Finished reading file (c:\users\nara_\appdata\local\temp\tmpi9opye\ante.in.mol2). (STP) Running: "C:/Program Files/Chimera 1.13.1/bin/amber16/bin/bondtype" -j part -i ANTECHAMBER_BOND_TYPE.AC0 -o ANTECHAMBER_BOND_TYPE.AC -f ac (STP) (STP) (STP) Running: "C:/Program Files/Chimera 1.13.1/bin/amber16/bin/atomtype" -i ANTECHAMBER_AC.AC0 -o ANTECHAMBER_AC.AC -p gaff (STP) Running: "C:/Program Files/Chimera 1.13.1/bin/amber16/bin/atomtype" -i ANTECHAMBER_GAS.AC -o ANTECHAMBER_GAS_AT.AC -d "C:/Program Files/Chimera 1.13.1/bin/amber16/dat/antechamber/ATOMTYPE_GAS.DEF" (STP) Charges for residue STP determined Total charge for #0: 0.000 #0, chain A: phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit γ isoform Traceback (most recent call last): File "C:\Program Files\Chimera 1.13.1\share\chimera\preferences\base.py", line 447, in save WindowsError: [Error 32] O arquivo j� est� sendo usado por outro processo Cannot replace previous preferences file: C:\Users\nara_\AppData\Roaming\.chimera\preferences: WindowsError: [Error 32] O arquivo j� est� sendo usado por outro processo File "C:\Program Files\Chimera 1.13.1\share\chimera\preferences\base.py", line 447, in save See reply log for Python traceback. Residues with incomplete side chains: MET 143.A GLU 322.A Replacing each by 'swapaa same (residue atom spec) lib Dunbrack preserve True' MET 143.A: phi none, psi -122.2 trans Chi angles for MET 143.A: GLU 322.A: phi -130.9, psi 172.4 trans Chi angles for GLU 322.A: Applying GLU rotamer (chi angles: -62.4 -179.4 -3.8) to GLU 322.A Applying MET rotamer (chi angles: -65.8 -61.8 -69.4) to MET 143.A Terminii for 3dbs.pdb (#0) chain A determined from SEQRES records Chain-initial residues that are actual N terminii: MET 143.A Chain-initial residues that are not actual N terminii: GLU 267.A, CYS 357.A, LEU 460.A, PHE 497.A, HIS 525.A, ARG 544.A, GLU 981.A Chain-final residues that are actual C terminii: Chain-final residues that are not actual C terminii: ALA 253.A, GLU 322.A, CYS 435.A, GLY 489.A, ASN 522.A, ALA 528.A, HIS 967.A, VAL 1091.A 718 hydrogen bonds Adding 'H' to GLU 267.A Adding 'H' to CYS 357.A Adding 'H' to LEU 460.A Adding 'H' to PHE 497.A Adding 'H' to HIS 525.A Adding 'H' to ARG 544.A Adding 'H' to GLU 981.A Removing spurious proton from 'C' of VAL 1091.A Hydrogens added Charge model: AMBER ff14SB Total charge for #0: -8.000 receptor and ligand cannot be the same model receptor and ligand cannot be the same model OpenGL Vendor: Intel OpenGL Renderer: Intel(R) HD Graphics OpenGL Version: 4.4.0 - Build 20.19.15.4835 Manufacturer: Dell Inc. Model: Inspiron 5452 TotalPhysicalMemory: 4203823104 ―――――――――――――――――――― Multisampling: False Shadows: False Shadow texture size: 2048 Silhouettes: False Depth cue: True Subdivision quality: 1.50 Single-layer transparency: True Transparent background: False Shaders supported: True Using shader: True Window size: 1366 663 Camera mode: mono Orthographic projection: False Center of rotation: front center Near/far clipping: False Key light: True Fill light: True Back light: False Ambient light: 0.20 Specular sharpness: 30.00 Specular reflectivity: 1.00 ―――――――――――――――――――― AntialiasLines: native AntialiasPoints: native BlendEquation: native BlendFuncSeparate: native BrokenAttribLocation: native, disabled ChoosePixelFormat: extension ColorTable: not supported CompileAndExecute: native CompiledVertexArray: extension CubeMap: native CullVertex: not supported DrawElementsInstanced: native, disabled DrawRangeElements: native FBConfig: not supported FBOShadows: native FastMultisampling: native, disabled FramebufferMultisample: native FramebufferObject: native LimitVertexAttribDivisor: native, disabled Multisample: native Multitexture: native PackedDepthStencil: native PalettedTexture: not supported PointParameters: native SeamlessCubeMap: native SeparateSpecularColor: native Shading: native Shadows: native StereoMultisample: native StereoRubberBanding: native Texture3D: native TextureColorTable: not supported TextureEdgeClamp: native TrustColorLogicBlend: native TrustNormals: native VertexArray: native VertexAttrib: native VertexAttribDivisor: native VertexBufferObject: native, disabled WindowPos: not supported

2 1

A question about the Volume Viewer
by Grosman, Claudio F 09 Aug '19

09 Aug '19

Hello everyone, I have (what I think is) a very simple question: what is the particular value of "level" (that is, the threshold) that appears automatically when one opens a 3D density map? After playing around with the threshold, I always feel like the initial value is the best one, but I wonder what exactly this value is. In the online documentation of Chimera, this issue is addressed. The text reads that the initial value of this threshold is such that 1% of the voxels have density values above it. Is this really the case? Only voxels with intensities above the threshold are displayed. Am I seeing only 1% of the voxels, then? I am afraid, I am confused. Thank you very much, in advance. Claudio

2 2

Selecting one ID 0 when two ID 0s are present
by Ryan Harlich 09 Aug '19

09 Aug '19

Hi Chimera, I have a problem where there are two ID 0. One is the normal density map the other is a solid (or what seems to be a movie of the protein). I want to select the density map. In the GUI interface I would open volume viewer and select radio button surface for category Style. Then a little below select All instead of One. I was wondering how to do this with command line? What I am trying to do is get a chain of the pdb structure and zone the density map. However the files are really big so I am trying to do it via ssh on another computer so I will be doing it with the –nogui flag. So if I can do everything via command line that would be great. Best, Ryan Harlich Sent from Mail<https://go.microsoft.com/fwlink/?LinkId=550986> for Windows 10

2 5

Double chain IDs
by Albert WEIXLBAUMER 08 Aug '19

08 Aug '19

Dear Chimera developers, My apologies if this question has been asked before but I was not able to find an answer so far. I am dealing with PDB files that contain too many individual chains and as such need two-character chain IDs. Ideally I do not want to split the PDB-files (to be able and work with single-character chain IDs). In addition, the PDB-file also contains more than 100,000 atoms and thus uses alphanumeric atom numbers. It seems Chimera is not aware of the two-character chain IDs because selecting a certain chain in the GUI shows only single character chain IDs and it leads to the selection of the two distinct chains (e.g. chain AA and chain BA are both treated as chain A). Is there a workaround? I’ve also noticed that saving such a PDB-file in Chimera changes all ATOM records to HETATM - I don’t know if this is related to the double-character chain IDs or alphanumeric atom numbers but this causes problems downstream because it’s very difficult to change this back to ATOM records without also affecting true HETATM entries. Thank you for your help! Best regards, Albert -- Albert Weixlbaumer Centre for Integrative Biology Department of Integrated Structural Biology Institute of Genetics and Molecular and Cellular Biology IGBMC - UMR 7104 - U 1258 1, rue Laurent Fries BP 10142 67404 ILLKIRCH CEDEX FRANCE Phone: +33 3 6948 5103 Email: Albert.Weixlbaumer(a)igbmc.fr<mailto:Albert.Weixlbaumer@igbmc.fr> OpenPGP KeyID: 0x7F644168

2 1

Making PDBs with over 10,000 atoms
by Ryan Harlich 07 Aug '19

07 Aug '19

Hi Chimera, I am trying to make a PDB with over 10,000 atoms. Currently we have this marked as a limitation to our project. I am researching to see if there is a way around this. Our current idea is to make the next 10,000 atoms another chain. Although when I apply this to a sample I get a break point between chain A and chain B. I do want them connected and to be different chains. So: 9,998.A 9,999.A 0.B 1.B connect. In addition, I would need helices and sheets to span over these two chains. So helix or sheet would start at 9,998.A and end at 1.B. Best, Ryan Harlich Sent from Mail<https://go.microsoft.com/fwlink/?LinkId=550986> for Windows 10

2 1

Superimpose protein chains
by Peter Snyder 07 Aug '19

07 Aug '19

How can I superimpose two (or three) chains from a single structure (ie two homologous subunits that form a multimeric protein)? Thank you Pete

2 2

Is "volume #0 voxelSize 1" a isosurface display changes
by Ryan Harlich 06 Aug '19

06 Aug '19

Hi Chimera, Is "volume #0 voxelSize 1" a isosurface display changes? So if I volume save this would it save this change? Best, Ryan Harlich Sent from Mail<https://go.microsoft.com/fwlink/?LinkId=550986> for Windows 10

2 3

Morphing, ligands, and secondary structure
by Yan Kung 06 Aug '19

06 Aug '19

Hi Chimera, I’m a new-ish user, and I’ve found that morphing structures is usually very straightforward and easy to use. But I’m now running into an issue with morphing and ligands (and secondary structure). I’d like to make a single morph composed of states A-B-C-B-A (so 4 movements between 5 states). These movements are: (1) A-B: domain opens, (2) B-C: a ligand flies out, (3) C-B: the same ligand flies back in, and (4) B-A: domain closes. First problem: the ligand does not appear in the resulting morph, a problem suggested by the documentation, which I have tried (and failed) to solve using the tips I found in archived posts from this mailing list. Also, because the ligand has to fly in and out, I can’t simply superimpose the ligand from any of the individual states over the morph. So I tried a workaround: make each movement a separate morph. This works for movements 2 and 3 (where the ligand flies out and in). But the ligand still does not appear in movements 1 and 4 (where the ligand remains in place and the domain moves). I’m not sure why the ligand’s appearance would depend on whether it flies in/out, or whether a domain moves… But it seems to. Moreover, a second problem emerges that has to do with secondary structure and is difficult to explain. (Essentially: a few protein regions are on the border between helices and loops such that they appear as helices in some states but as loops in others, and vice versa. If possible, I’d like the secondary structure to remain constant across all states, or else it’s distracting. Even still, I’ve tried to embrace this and use ksdssp to show the helices/loop switching back and forth, but this is distracting as they flicker between helices and loops, each region at its own separate frame. Specifying ksdssp only in certain frames to minimize the flickering doesn’t work either, as different regions switch back and forth at different times.) The only solution I can think of that would maintain the same secondary structure throughout is to just make a single morph of all movements (and not use ksdssp). But this puts me back where I started! Any ideas? (Separately, I realize that I could make a simpler morph of just A-B-C, then record it as “roundtrip”, i.e. forward and backwards, to achieve the same thing as A-B-C-B-A. But this doesn’t solve the problems above.) Thanks very much, Yan -- Yan Kung Assistant Professor Department of Chemistry Bryn Mawr College http://www.brynmawr.edu/people/yan-kung Follow: @KungLab<https://twitter.com/KungLab>

2 1

Nvidia 3D Vision on a large display
by Tom Goddard 02 Aug '19

02 Aug '19

Krish Raman asked what kind of large displays might work with Nvidia 3D Vision stereo glasses with Chimera. Hi Krish, I believe any display capable of a high refresh rate, at least 120 Hz (even 100Hz would do) will be able to show stereo with Nvidia 3D vision. Conventional 60 Hz displays could only show stereo at 30 Hz in each eye (since the display first shows the left eye, then shows the right eye in sequence), and 30 Hz is too slow causing flicker. Video gaming displays are usually advertised as 144 Hz and I believe any of those would work. It is hard to find big displays at 144 Hz, but here is an example, 49 inch ultra wide Samsung CHG90 https://www.amazon.com/Samsung-49-Inch-Curved-Monitor-LC49HG90DMNXZA/dp/B07… <https://www.amazon.com/Samsung-49-Inch-Curved-Monitor-LC49HG90DMNXZA/dp/B07…> If one person is using this sitting at a desk it is very wide and would probably give very nice immersive stereo because of the width. It is not ideal for your use with an audience standing further back because of its extreme 32:0 aspect ratio. You could look for a video gaming TV, such as these https://www.tomsguide.com/us/best-4k-gaming-tv,review-4837.html <https://www.tomsguide.com/us/best-4k-gaming-tv,review-4837.html> Keep in mind your Nvidia 3DVision system may not work with any of these. I would strongly recommend testing at a store with your 3DVision before buying. Stereo is an an extremely obscure niche product and there is no telling if it will work without testing. Another option is to use a high refresh rate projector. That is what we use at my lab -- (ours is a $50,000 Christie, so I don't recommend it). Also I am sure you are aware that Nvidia discontinued support for 3D Vision and it will not work with any Nvidia drivers released after April 2019. So the 3D Vision system is going to be dead soon. Tom

1 0

crash on using Preset Publication 2
by Dieter Blaas 02 Aug '19

02 Aug '19

Hi, I have been using the Publication Preset 2 for making figures of a number of pdb structures without problems. However, I now have a pdb that shows up well in Chimera and does not have any obvious problems. However, when trying rendering it with the Preset Publication 2 (and only with this one) Chimera crashes with the error message below: *** stack smashing detected +++: /home/blaas/.local/UCSF-Chimera64-1.10.2/bin/python2.7 terminated Since I have prepared the other figures with this setting I'd like using it as well for the next figure but how? Thanks a lot, best regards, Dieter ------------------------------------------------------------------------ Dieter Blaas, Max F. Perutz Laboratories Medical University of Vienna, Inst. Med. Biochem., Vienna Biocenter (VBC), Dr. Bohr Gasse 9/3, A-1030 Vienna, Austria, Tel: 0043 1 4277 61630, Fax: 0043 1 4277 9616, e-mail: dieter.blaas(a)meduniwien.ac.at ------------------------------------------------------------------------

2 1

Chimera Viewdock Bug
by Sacco, Michael 02 Aug '19

02 Aug '19

Hello all, I would like to report a bug. When using Chimera in dual monitor mode on windows OS and loading the viewdock module, Chimera becomes unresponsive. However, when the second monitor is unplugged, it functions normally. I have confirmed this on multiple computers. After force closing Chimera a pop-up appears with the message "Chimera exited abnormally; the exit code was: 1.". Thanks everyone! Michael

2 1

OBJ files are exported too small and won't open on other programs.
by Josh Moon 01 Aug '19

01 Aug '19

Hello! I am a relatively new and inexperienced user of Chimera. I am having a difficult time working with obj files. The Scenes I'm exporting as obj's end up having a very small file size (67 bytes in one case) and are being processed very quickly. The resulting file won't open successfully on any program that should understand an obj. Stl's are fine. Conversely, when I try to open an obj file into Chimera, it asks me to manually "Designate the File Type" as though it doesn't know what to do with it. Can someone please suggest what features I'm missing so that I can successfully export an OBJ? I'm also interested in how best to handle color information. Thank you! Josh ----------------------------------------- Josh Moon, Ph.D. (he/him/his) Educational Technology Specialist Kalamazoo College (269) 337-7415

2 1

inquire about molecular docking
by Mohamed Hegazy 31 Jul '19

31 Jul '19

Hey , My name is Mohamed Hegazy, beginner in molecular docking , recently I performed a molecular docking procedure using chimera and I found that the initial structure of the ligand differs from the minimized one , and the binding score for a ligand would change during my different runs. The binding score you obtain from my docking cannot be a constant number. Sent from Mail for Windows 10 Mohamed Said Hegazy --- This email has been checked for viruses by Avast antivirus software. https://www.avast.com/antivirus

2 4

AreaSAS Issue
by Rucco,Dominic J 31 Jul '19

31 Jul '19

Hello, I'll start by saying your software has helped me visualize protein structure for quite some time and overall has been easy to use! I do have a question however: I was recently using the "render by attribute" feature to display areaSAS measurements but have run into an issue. I can no longer select this attribute from the dropdown menu. I noticed this when I attempted to render it on the protein mCherry (PDB 2H5Q). Even when I have gone back to other PDB entries I know have worked for areaSAS rendering, it doesn't appear as an option anymore. The only time I can get it to appear is when I load a session that already has the rendering performed on the protein. Any help would be appreciated, thank you, Dominic Dominic Rucco | Ph.D. Candidate George and Josephine Butler Polymer Research Laboratory Center for Macromolecular Science and Engineering Department of Chemistry University of Florida Leigh Hall 314 | PO Box 117200 | Gainesville, FL 32611-7200 drucco(a)ufl.edu<mailto:drucco@ufl.edu> | www.savin.chem.ufl.edu<http://savin.chem.ufl.edu/>

2 1

--nogui gets permission denied when fetching pdb id
by Ryan Harlich 30 Jul '19

30 Jul '19

Hi Chimera, I am trying to run a script with the –nogui option and run command “open 1a1x” and I get permission denied. It works with the GUI so I am not sure why it is not working with –nogui. Best, Ryan Harlich Sent from Mail<https://go.microsoft.com/fwlink/?LinkId=550986> for Windows 10

3 5

Minimization with sulfur atom
by sara20＠mail.ustc.edu.cn 30 Jul '19

30 Jul '19

Dear users I have complex structure with ligand having sulfur atom in it. I am trying to minimize it but it fails giving the following error: (LIG) /home/usr/.local/UCSF-Chimera64-1.12/bin/amber16/bin/antechamber: Fatal Error! (LIG) Cannot properly run "/home/usr/.local/UCSF-Chimera64-1.12/bin/amber16/bin/sqm -O -i sqm.in -o sqm.out". I assume it's because of sulfur atom in ligand because I tried running it with out sulfur and it worked fine. I checked MMKT Atom database for sulfur file and it's already there: name = 'sulphur' symbol = 'S' mass = [(31.97207070, 95.02), (32.97145843, 0.75), (33.96786665, 4.21), (35.9670806, 0.02)] color = 'yellow' vdW_radius = 0.18 b_coherent = 2.847*fm b_incoherent = 0.188*fm b_coherent_deut = 2.847*fm b_incoherent_deut = 0.188*fm Any help would be highly appreciated. Thanks Sara

2 1

Applying session parameters to another pdb file.
by David Sáez 29 Jul '19

29 Jul '19

Dear Chimera experts, I have saved a session of Chimera showing the perfect angle and size for the molecules of interest. Now I would like to apply the same parameters to visualize another (but similar, only differnet coordinates) pdb file. Is there a way to do that? Thanks in advance. David

2 2

Editing sequences within Chimera
by Joshua Plank 26 Jul '19

26 Jul '19

Hi, I am wondering if it is possible to edit the amino acid sequence somehow? I have a .pdb file of the protein I'm interested in, but would like to change two of the amino acids to something else. Is this possible? thanks, Josh -- *Joshua Plank* *PhD Student* *Department of Electrical and Computer Engineering* *University of California, Riverside* *email: jplan002(a)ucr.edu <jplan002(a)ucr.edu>* *phone: 505-977-9919*

3 4

Measuring the distance between residues and a surface
by Stanley, Julian A 25 Jul '19

25 Jul '19

Hello! I'm fairly new to Chimera, so if you're able to point me towards any resources that can help me answer questions like this one independently in the future, please do! Here's the question: I would like to measure the distance between the ends arginine and lysine residues and a surface, using Chimera's python modules. Typically, I would load my .pdb file, generate a surface, and then use the "measure distance" command and parse the output from the Reply Log. This approach works fine, except since "measure distance" prints to the reply log, I need to run my script in GUI mode. If possible, I'd prefer to run everything from the command line only. I thought that I might be able to do this with the distance() command, but I can't figure out how to pass a surface to that command. Does anyone have any tips/advice? In summary, what I usually do: ``` from chimera import runCommand as rc # load the file rc("open " + pdbFile) # Generate a surface rc("select protein") rc("split") rc("surface allComponents false") rc("~select") # Measure distance between arg and lys and the surface rc("select #0") rc("~select #0:@") rc("measure distance :arg@cd|:lys@ce selection multiple true show true") ``` And then I parse out the distances from the Reply Log. What I would prefer to do: It would be better to have the output to go to STDOUT when I run the script. I imagine that I'd use something like the `distance()` command, of the premade scripts (like here: http://plato.cgl.ucsf.edu/trac/chimera/attachment/wiki/Scripts/atomdist.py) but they are designed to measure the distance between atoms or residues, not between residues and a surface. Thanks for your help! Sincerely, Julian

3 2

Modeller Restraints
by Danilo Boškovic 24 Jul '19

24 Jul '19

Dear all, how do i add restraints specific to my templates that i use for comparative modelling in chimera? Specifically what are the names assigned to those templates? Do i need the names to specify the template residues that i want to restrain? Are the names same as in the tab when i choose the templates for the run? Thanks, D. B.

3 5