- Chimera-users - RBVI Mailing Lists

Re: [Chimera-users] Request for axes measurements
by Elaine Meng 14 Jun '15

14 Jun '15

Hi, You’re welcome, but I don’t have any further details than what is already given in the manual: <http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/structuremeas/str…> It’s recommended you send Chimera questions to chimera-users(a)cgl.ucsf.edu instead of me directly, in case others may be able to provide better answers! Best, Elaine ----- Elaine C. Meng, Ph.D. UCSF Computer Graphics Lab (Chimera team) and Babbitt Lab Department of Pharmaceutical Chemistry University of California, San Francisco On Mar 31, 2015, at 6:26 PM, Jayaraman T <jram(a)pusan.ac.kr> wrote: > > Hi Elaine Meng, > Yes i could figure now the axes as well as angle measurements between them and thanks for your timely help. > > I would be happy if you could tell me the principle on which these axes are calculated. I have noticed that it's from Eigen vector values that are calculated for the coordinates of the set of atoms. > > I'm really interested on the backend for this calculation particularly, might help a bit progressing more for my research. > > Thanks again.

2 6

on MD movie and morph conformation by chimera
by Smith Liu 12 Jun '15

12 Jun '15

Dear All, Will you please introduce to me on how to combine chimera MD movie function and morph conformation function to make conformation change movie files to insert into PPT files? Smith

2 1

mol2 for NH3
by MPI 11 Jun '15

11 Jun '15

Dear Users, I try to prepare a mol2 from pdb for ammonia (NH3) and ammonium (NH4+), respectively. I used two commands: addh and addcharge with am1 method. eg, addh addcharge all method am1 It turns out that Chimera gave two similar mol2 files for NH4+ but NOT for NH3. I wonder what command and parameters are needed to produce a correct mol2 for NH3. Here are the input pdb files for NH3 and NH4, respectively. # NH3 HETATM 2860 N NH3 A 402 53.512 -12.964 -36.364 0.88 45.66 N # NH4 HETATM 2858 N NH4 A 404 56.362 -12.467 -27.097 1.00 30.80 N Thanks, Dewey

2 5

Problems detecting a left-handed helix
by Felipe Vasquez 11 Jun '15

11 Jun '15

Hi, I am trying to visualize several of the secondary structure details on a protein of interest, and I recently had notice (by several research papers) about the existence of a left-handed helix on it. However, I could not detect the helix using UCSF Chimera (it is simply considered as a loop). Why could this happen and how could be resolved? I have tried to modify the ksdssp cut-off parameter, but it was unsuccessful, so I think now that maybe the specific nature of the helix (left-handed polyproline-like type II) could be a bit tricky for the program. What would you suggest? Thanks in advance for your help an patience. Best regards, *Andrés Felipe Vásquez J., BSc, MSc.* Grupo de Fisiología Molecular Subdirección de Investigación Científica y Tecnológica Dirección de Investigación en Salud Pública Instituto Nacional de Salud Avenida calle 26 No. 51-20 - Zona 6 CAN +57 (1) 2207700 ext. 1419 Bogotá, D.C., Colombia

2 3

Determination of angles formed by distances and ring axes
by Felipe Vasquez 10 Jun '15

10 Jun '15

Hi, I have problems trying to calculate the angle between two lines: one is really built as the distance between a residue atom from a protein and the centroid of a benzene ring (from a ligand), and the other is the line corresponding to the axis of the benzene ring. How could I accomplish this objective? I have read similar issues by chimera-users on the web, but I could not find enough information to complete the calculation (I could only estimate the first distance). In case of problems to clearly comprehend this particular task, please refer to the page 3, Fig 1, panel C in the attached research paper (phi angle). PD: Is it normal that the centroid of a benzene ring is above the ring itself? Because in the process of build the first distance, I don't know whether I have to calculate the distance between the centroid and the residue atom (which I found to be shorter), or instead calculate the distance between the center itself of the plane and the residue atom (which would be a bit longer). Thanks in advance for your time and assistance in this matter. Best regards, *Andrés Felipe Vásquez J., BSc, MSc.* Grupo de Fisiología Molecular Subdirección de Investigación Científica y Tecnológica Dirección de Investigación en Salud Pública Instituto Nacional de Salud Avenida calle 26 No. 51-20 - Zona 6 CAN +57 (1) 2207700 ext. 1419 Bogotá, D.C., Colombia

2 1

file size of a trajectory
by MPI 10 Jun '15

10 Jun '15

Dear Users, I use the lastest version of Chimera to read trajectory (.xtc) files output from Gromacs ver 4.6.5. When a .xtc file size larger than 2GB in most cases and 1GB in some cases, Chimera gets choked and crashed. When a .xtc file size < 1GB, Chimera works very well. I wonder if this is a limitation for reading Gromacs trajectory files. Is there a solution ? Thank you. Dewey

3 4

tetrahedral symmetry
by Morris, Kyle 10 Jun '15

10 Jun '15

Dear Chimera Developers, I have a volume reconstructed from Relion 1.3 imposing tetrahedral symmetry. Curiously, I am not able to find the symmetry group using: $ measure symmetry #0 Chimera returns C3, even with adjusting the parameters it uses to identify the symmetry group. The model is aligned as following the Xmipp convention for T symmetry (http://xmipp.cnb.csic.es/twiki/bin/view/Xmipp/Symmetry#T) which is the same as what Chimera calls T,z3? Any advice on getting the measure command to work in this convention would be greatly appreciated. Thanks! Best wishes, Kyle Morris University of Warwick, UK

2 1

Save Cluster Figure as Phython Object
by Beuerle, Matthias 09 Jun '15

09 Jun '15

Dear all, I have been using the Clustering tool in the MD Movie/Analysis/Cluster to examine my trajectories and I found it very useful. Now I would somehow like to save the box-chart matplotlib object in order to adjust the figure, so that I can use it for a publication, since the dialog buttons do not allow a lot of modifications. Is this possible somehow? Alternatively I would be very happy as well if one could use the output obtained from the 'Save' button to recreate a similar figure (in which I have failed so far). Help is very much appreciated. Kind regards, Matthias

2 1

joining two proteins
by López de Heredia Alonso, Miguel 08 Jun '15

08 Jun '15

Hi, I am new to Chimera, so perhaps someone could give me a hand. I have used Robetta for predicting the structure of a protein I am interested in and it gave back models for different parts of the protein. I would like to join them, so I am trying to join two models by a C-N peptide bond. I have selected the corresponding atoms in the last amino acid in one of the proteins and in the first amino acid of the other protein I want to join. With both atoms selected I took the "Join Models" option in "Build Structure" with the "C-N peptide bond" selected. And the "Apply" button does not activate. (See the screen shot attached) I have gone through the manual and I do not see what I am missing or what I am doing wrong. Any help is welcome. Thank you very much in advance. Kind regards, Miguel

2 1

APBS electrostatic potential
by Felipe Vasquez 06 Jun '15

06 Jun '15

Hi, I am interested on calculating the electrostatic potential of a molecular surface on a protein (specifically visualizing that potential on a previously cavity determined by CASTp). However, althought I did get reasonable results using Coulombic potential, I could not obtain similar results using APBS potential (I obtained a surface too electrostatic, charged negatively and the values have a huge range, not only from -10 to 10 as in the case of Coulombic potential). Should I modify some default values out of the calculation parameters in order to get a similar result to that obtain by Coulombic potential? And it is absolutely indispensable add hydrogens in an APBS estimation? Thanks in advance for your help. Best regards, *Andrés Felipe Vásquez J., BSc, MSc.* Grupo de Fisiología Molecular Subdirección de Investigación Científica y Tecnológica Dirección de Investigación en Salud Pública Instituto Nacional de Salud Avenida calle 26 No. 51-20 - Zona 6 CAN +57 (1) 2207700 ext. 1419 Bogotá, D.C., Colombia

2 2

surfaces and
by A G Szabo 03 Jun '15

03 Jun '15

I have been working with a very large structure of many subunits, the E.Coli Ribosome particle with the Signal Recognition Particle attached (PDB 2J28) I was able to create an image of the SRP part of the complex. I selected that subunit. And was able to hide all the other subunits. So with the subunit selected I tried to visualize only that image as a surface. Wouldn't let me, getting an error message that it was too large and I should somehow split it. Don't know how to do that. I know the whole dataset is very large, but the SRP portion is only 400 amino acids. A work around would be if I could somehow select the SRP subunit and have it as a dataset to reload as that subunit alone. Then not so big. Don't know how to do that either. Your suggestions. I couldn't find anything in the tutorials that might help. I was able to show the surface of another PDB file that was much smaller. Arthur G. Szabo Professor Emeritus Chemistry Wilfrid Laurier University Waterloo Ontario Canada.

2 1

Re: [Chimera-users] FW: Python code problem
by Mahendra B Thapa 03 Jun '15

03 Jun '15

Dear Eric I had not realized that not every amino acid consisted of CB atom. Thank you for the code which worked; just want to add a very minor addition ( """ ) to the left of the code, i.e. *sel = OSLSelection("""#0.1:%d.B@CB""" %i)* instead of *sel = OSLSelection(#0.1:%d.B@CB""" %i)* I tried another code as import string rc("open " + 'test.pdb') # 60 models, each model has chains A,B,C for i in [2196,2197]: for j in range(3429,3450): for *k* in ['A','B','C']: sel = OSLSelection("""#0.3:%d.%s@CA""" %(j,k)) sol = OSLSelection("""#0.2:%d.B@CA""" %i) if ((len(sel.residues()) != 1) or (len(sel.residues()) != 1)): continue if ((len(sol.residues()) != 1) and (len(sol.residues()) != 1)): continue rc("""distance #0.2:%d.B@CA #0.3:%d.%s@CA""" %(i,j,k)) I got the output as (without any error message) Distance between #0.2:2196.B@CA and #0.3:3450*.C*@CA: 104.315 The 'for' loop did not work for *k* (above) and only gave results for C only but not for A and B. Thank you for help, Mahendra Thapa University of Cincinnati,OH On Fri, May 29, 2015 at 1:55 PM, Thapa, Mahendra (thapamb) < thapamb(a)mail.uc.edu> wrote: > > ------------------------------ > *From:* Eric Pettersen > *Sent:* Friday, May 29, 2015 11:55:18 AM (UTC-06:00) Central America > *To:* Thapa, Mahendra (thapamb) > *Cc:* chimera-users(a)cgl.ucsf.edu > *Subject:* Re: [Chimera-users] Python code problem > > Hi Mahendra, > In the first case, as you saw, it tries exactly two distances; the first > involving :2011.B and the second involving :2420.B. In the second case it > will try every distance from :2011.B through :2419.B (note that the range > command does not include the ending number in the range). This means if > the structure is missing *any* of those residues then the distance > command will fail with the error you got, since the atom spec will only > select one atom. > So you either need to explicitly use just the residue numbers that exist, > or you need to test if the residue you are about to try to use exists. So > code to do that might be: > > from chimera.selection import OSLSelection > sel = OSLSelection(#0.1:%d.B@CB""" %i) > if len(sel.residues()) != 1: > continue > > --Eric > > Eric Pettersen > UCSF Computer Graphics Lab > http://www.cgl.ucsf.edu > > On May 29, 2015, at 6:37 AM, Mahendra B Thapa <thapamb(a)mail.uc.edu> > wrote: > > Dear Chimera users, > > The command I used was > chimera --nogui test.py > I tried two codes (following previous posts and tutorials in chimera) in > 'test.py' as follows > > (i) > import os > from chimera import runCommand as rc > from chimera import replyobj > rc("open " + 'test.pdb') > for i in [2011,2420]: > rc("""distance #0.2:2010.B@CB #0.1:%d.B@CB""" %i) > > *This code gave the result as expected*: > Executing test.py... > Opening test.pdb... > 60 models opened > Distance between #0.2:2010.B@CB and #0.1:2011.B@CB: 72.346 > Distance between #0.2:2010.B@CB and #0.1:2420.B@CB: 92.979 > Executed test.py > > (ii) > import os > from chimera import runCommand as rc > from chimera import replyobj > rc("open " + 'test.pdb') > for i in range(2011,2420): > rc("""distance #0.2:2010.B@CB #0.1:%d.B@CB""" %i) > > * This code gave the result as* > ' You selected %d.' % numSelected) > MidasError: Exactly two atoms/axes/planes must be selected. You selected > 1. > Error while processing test.py: > MidasError: Exactly two atoms/axes/planes must be selected. You selected > 1. > > File "/opt/UCSF/Chimera64-1.5.3/share/Midas/__init__.py", line 1266, in > distance > ' You selected %d.' % numSelected) > > What could be the error this code? Only difference between two codes is: > "for > * i in [2011,2420]:" & "for i in range(2011,2420):" * > That is, I am interested to use range operator in 2nd case instead of two > values as in the 1st case. > > Thank you for help, > Mahendra Thapa > University of Cincinnati,OH > > > _______________________________________________ > Chimera-users mailing list > Chimera-users(a)cgl.ucsf.edu > http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users > > > > >

2 1

savepos
by Dougherty, Matthew T 03 Jun '15

03 Jun '15

Hi, I have been automating the creation of positions using python scripts. Something relatively simple like p001, p002, etc. When I examine the positions using the python commands from Midas import positions pnames = positions.keys() at first things appear as I would expect, but later the sequence changes relative to the index, why is that? for example, using this to print the positions pnames = positions.keys() i=0 for pname in pnames: print i,pnames[i] i=i+1 I get these results: savepos p000 savepos p001 0 p000 1 p001 (as expected) savepos p002 0 p000 1 p001 2 p002 (as expected) savepos p003 0 p000 1 p001 2 p002 3 p003 (as expected) savepos p004 0 p004 1 p000 2 p001 3 p002 4 p003 (not as expected) savepos p005 0 p004 1 p005 2 p000 3 p001 4 p002 5 p003 (not as expected) I would expect pnames[0] to correspond to the first position created, pnames[1] to the second posiiton created, etc., assuming no ~savepos used. Want to create a situation where the default names would be pxxx, but could later be altered to more meaningful names but maintain the order of creation, or getting a session file from someone I would use their position names and tack on my default names to the end of the position list. Matthew Dougherty National Center for Macromolecular Imaging Baylor College of Medicine ===========================================================================

2 1

to color EM map
by Smith Liu 03 Jun '15

03 Jun '15

Dear All, Suppose I have an EM map for a protein complex of 4 subunits. I have fitted them in chimera. Will you please tell me how to give each part of the map a different color based on the subunits? For example, for subunit A, B,C, D part of the EM map, how to color the map accordingly in blue, red, yellow and purple? Smith

2 1

volume viewer/Data menu
by Dougherty, Matthew T 02 Jun '15

02 Jun '15

Hi, What is the pattern for the sequence of the Data list? Is there a way to sequence it the same as the Model Panel? Matthew Dougherty National Center for Macromolecular Imaging Baylor College of Medicine ===========================================================================

2 1

Unable to Change Lighting of Side Chains / Atoms
by David Nelson Bunck 02 Jun '15

02 Jun '15

Hello, Recently, for some reason, the settings changed on Chimera such that the atoms and side chains do not change color when I try to change the lighting settings. Thus when I get the protein (in ribbon format) to look appealing and shiny, the metal atoms in the protein remain black. Changing colors or lighting do not make any difference. The same is true for the side chains as well. This also occurs on several different sequences, so it is not just the one file. Sorry for the seemingly trivial question, but I have poked around for a while at the different settings and can't seem to find how to undo this change. Thanks for your help! david -- david bunck bunck(a)caltech.edu (626) 395-8919 Postdoctoral Scholar Heath Research Group Division of Chemistry and Chemical Engineering California Institute of Technology

2 2

Cage maker
by Anindito Sen 02 Jun '15

02 Jun '15

Dear All, Is it possible to generate a synthetic 3d density map of hexagons and/ pentagons using the Cage building option with file extension like mrc or spi or pif ( the usuals). If so, how to save it in such a format ? Thanks Andy Sent from my iPad

2 1

H-bond angle criteria (cut-off)
by Felipe Vasquez 02 Jun '15

02 Jun '15

Hi, I am using the H-bond function in UCSF Chimera and I noted that the relaxed parameter for angle criteria is 20 by default. However, since that I have no hydrogen atoms in the pdb structure, I assume that the angle theta (donor-hydrogen-acceptor) cannot be used by Chimera as an angle criteria, and so the psi or phi angles might be used instead. Which angle is actually used in the relaxed parameters? And most important, what is the default value of that angle (angle minus relaxing value)? Thanks in advance for your help. Best regards, *Andrés Felipe Vásquez J., BSc, MSc.* Grupo de Fisiología Molecular Subdirección de Investigación Científica y Tecnológica Dirección de Investigación en Salud Pública Instituto Nacional de Salud Avenida calle 26 No. 51-20 - Zona 6 CAN +57 (1) 2207700 ext. 1419 Bogotá, D.C., Colombia

2 1

Re: [Chimera-users] Mouse control
by Elaine Meng 02 Jun '15

02 Jun '15

Dear Prof. Szabo, I can’t see your image, but the following might help: menu: Actions… Focus - OR - command: focus … if you don’t have anything selected, either of the above should center and resize to frame all of the displayed stuff in the main window. By “selected” I mean when you Ctrl-click on something in the main window, it is highlighted with green outline. If you have something selected, you can use that same menu action to focus only on the selected part (the corresponding command would then be “focus sel”). <http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/menu.html#menuactions> <http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/focus.html> Then, if needed, you can enlarge by zooming with the mouse or in the Side View, or with the “scale” command. I hope this helps! For future questions, we recommend you send to chimera-users(a)cgl.ucsf.edu instead of only to me, so that others can see the answers, or give better answers than mine. Regards, Elaine —— Elaine C. Meng, Ph.D. UCSF Computer Graphics Lab (Chimera team) and Babbitt Lab Department of Pharmaceutical Chemistry University of California, San Francisco On Jun 2, 2015, at 1:27 PM, A G Szabo <agszabo(a)bell.net> wrote: > Elaine > You were helpful in early April I thought I would ask another question. > > I have been working with a PDB file, 2J28, that is the structure of a > full ribosome with a Signal Recognition Protein and piece of RNA bound to > the ribosome. I figured out how to get rid of all the various parts of the > Ribosome, and visualize only the SRP and piece of RNA attached to it. > > The problem is that I can't figure out how to centre that part of the > overall structure. If I move the SRP/RNA with the wheel on the mouse it > does move to the center of the screen but in the process it becomes very > small. > > I need to be able to transpose the SRP/RNA and then enlarge it. > > Can't figure it out and while there is some stuff in the tutorials and > user info, I couldn't find anything that pertained to what I needed to > do. > > This is the SRP - RNA structure from Protein Workshop after all the > ribosomal particles were hidden. > > > > The image in Chimera is even smaller. > > Thank you again for your help. > > Arthur G. Szabo > Professor Emeritus > Chemistry > Wilfrid Laurier University > Waterloo Canada. > > <winmail.dat>

1 1

Aligning on crystal axis
by Mohamed Noor 02 Jun '15

02 Jun '15

Hello. Is there any easy way to align the view along the crystal axes, such as the ab, ac and bc faces separately? Thanks.

2 1

Python code problem
by Mahendra B Thapa 29 May '15

29 May '15

Dear Chimera users, The command I used was chimera --nogui test.py I tried two codes (following previous posts and tutorials in chimera) in 'test.py' as follows (i) import os from chimera import runCommand as rc from chimera import replyobj rc("open " + 'test.pdb') for i in [2011,2420]: rc("""distance #0.2:2010.B@CB #0.1:%d.B@CB""" %i) *This code gave the result as expected*: Executing test.py... Opening test.pdb... 60 models opened Distance between #0.2:2010.B@CB and #0.1:2011.B@CB: 72.346 Distance between #0.2:2010.B@CB and #0.1:2420.B@CB: 92.979 Executed test.py (ii) import os from chimera import runCommand as rc from chimera import replyobj rc("open " + 'test.pdb') for i in range(2011,2420): rc("""distance #0.2:2010.B@CB #0.1:%d.B@CB""" %i) *This code gave the result as* ' You selected %d.' % numSelected) MidasError: Exactly two atoms/axes/planes must be selected. You selected 1. Error while processing test.py: MidasError: Exactly two atoms/axes/planes must be selected. You selected 1. File "/opt/UCSF/Chimera64-1.5.3/share/Midas/__init__.py", line 1266, in distance ' You selected %d.' % numSelected) What could be the error this code? Only difference between two codes is: "for * i in [2011,2420]:" & "for i in range(2011,2420):"* That is, I am interested to use range operator in 2nd case instead of two values as in the 1st case. Thank you for help, Mahendra Thapa University of Cincinnati,OH

2 1

Ramachandran plot for a selection
by Elena Lilkova 28 May '15

28 May '15

Dear Chimera community, First I would like to just say that I am new to Chimera. I made a MD simulation of a protein using NAMD. Now I am interested in generating a Ramachandran plot for the whole trajectory, but only for a specific selection of my protein, namely three specific residues. I fail to succeed. Could anyone probably suggest a solution? Elena Lilkova elilkova(a)phys.uni-sofia.bg PhD student, Atomic Physics Department, Faculty of Physics, University of Sofia "St. Kliment Ohridski" 5 James Bourchier Blvd. 1164 Sofia, Bulgaria Elena Lilkova elilkova(a)phys.uni-sofia.bg PhD student, Atomic Physics Department, Faculty of Physics, University of Sofia "St. Kliment Ohridski" 5 James Bourchier Blvd. 1164 Sofia, Bulgaria

3 3

Distances in biological assembly
by Mahendra B Thapa 27 May '15

27 May '15

Dear Chimera users The biological assembly of a virus in pdb format consists of 60 models; each model with 3 chains :A, B and C. I am interested to find atoms within 4 angstroms of the atom, say CA atom of residue 2059 of chain B of model 2. The chimera crashed in graphical interface (it could be due to size of the pdb file.), so I used the following command [ http://plato.cgl.ucsf.edu/pipermail/chimera-users/2007-November/002069.html ]: chimera --nogui file_test.cmd > data_rest.log where 'file_test.cmd' consists of open test.pdb findclash #2:2059.B@CA overlap -4 hbond 0 The output file 'data_rest.log' consisted of Opening test.pdb... 60 models opened Opened test.pdb containing 60 models, 677040 atoms, and 89700 residues No contacts No contacts *I want to know whether I have used right chimera command /code. *Any suggestion might help me. Thank you, Mahendra Thapa University of Cincinnati,OH

2 1

scale bar
by Dmitry 27 May '15

27 May '15

Dear developers, How are you? I have a question about scale bar. I am not completely sure that the scale bar I drew is correct. What are the scale settings for drawing the scale bar? By the default the scale bar is in Angstroms isn’t it? In ImageJ there is a possibility to set the scale 1 pixel = X value (nm/Å/etc) Is there the same possibility in Chimera? Thank you Sincerely, Dmitry

3 3

Starting multiple copies of Chimera on Mac
by Tom Goddard 26 May '15

26 May '15

It is possible to run multiple copies of the same Chimera on a Mac. Here’s how. From a Terminal (/Applications/Utilities/Terminal) type the command open -n -a Chimera The -n option tells the Mac to start another copy of the application. If you prefer to start the second copy by clicking an icon then you can make that icon using /Applications/Automator, choose to make an Application, select Run Shell Script, enter the command open -n -a Chimera as the script to run, then File / Save… to save this application as “Chimera+” or whatever you want to call it. This icon always starts a new copy of Chimera so you can have as many running as you like. Tom

1 0

Re: [Chimera-users] On the seeking of the orthologs
by James Starlight 25 May '15

25 May '15

Elaine, George, thank you very much for the suggestions and ideas!! One question- do you know some DB where is possible to seek for the orthologs (e,g based on HMM) agains *several* target sequences with flexible options like to detect minimal percent of the homology for orthologs which to be found or alternatively to make such search based on some *shared* features between those targets (e.g common structural motifs responsible for ligand binding etc) which might be specified by user initially? George, the idea to focus on particular specie seems interesting- initially I've thought in the same way but to focus on the human OR (because here I'm more comfortable with dynamical range patterns of those receptors based on the in vivo Matsunami data for instance). BTW briefly have a look on those DB I didn't find there engine to seek for the orthologs agains given target. Doest this option really exist within it? All the best, James 2015-05-22 21:10 GMT+02:00 George Tzotzos <gtzotzos(a)me.com>: > I would start from BioMart. I’m also working on OR (albeit insect ones). Try also MGI (Mouse Genome Informatics). Starting from mouse Ors you can scoop orthologs in other mammalian species. > > Cheers > > George > >> On 22 May 2015, at 18:44, Elaine Meng <meng(a)cgl.ucsf.edu> wrote: >> >> Hi James, >> Your questions are more about bioinformatics project design rather than Chimera use. >> >> Chimera provides basic BLAST searching of NCBI’s “nr” and “pdb” sequence databases, but these are the same databases generally used by BLAST whether you use it from Chimera or from another program, or stand-alone. Also, Chimera can fetch specific individual sequences from UniProt based on UniProt identifier, but that isn’t searching, it is just getting data you already knew you wanted. >> >> I believe NCBI’s “nr” and UniProtKB are both meant to be comprehensive (cover “all” sequences) but in their own ways, and you would have to go to those resources or their publications to learn the details. There are probably also specialized GPCR and/or olfactory receptor sequence databases. Then there are multiple ways of searching, with individual sequences or with derived patterns such as hidden Markov Models (HMMs). This is a really broad topic and in my opinion, beyond the scope of this forum. You’ve probably done some of this already, but I would recommend reading similar studies in quality journals by authors you respect, to see how they have designed their projects. >> >> Chimera doesn’t have HMM-related tools or evolutionary trace, it just allows conservation analysis on your whole alignment, or subsets of that alignment opened in separate windows. >> >> I hope this helps, >> Elaine >> ----- >> Elaine C. Meng, Ph.D. >> UCSF Computer Graphics Lab (Chimera team) and Babbitt Lab >> Department of Pharmaceutical Chemistry >> University of California, San Francisco >> >> On May 22, 2015, at 9:22 AM, James Starlight <jmsstarlight(a)gmail.com> wrote: >> >>> Dear Chimera users! >>> >>> I'd like to ask some suggestion regarding my bioinformatician task and >>> possible ways of its realisation using some bioinformatics methods >>> implemented in Chimera. >>> >>> 1)I need to seek for some sequence database (here I don't know >>> absolutely how to do it) of the olfactory receptors and to select 10 >>> sequences of the Ors corresponded to 10 different mammalian species -> >>> according to my idea it should be very different species associated >>> with different ecological niches so its olfactory spectrum should vary >>> dramatically. >>> Here the question: in what DB such search should be better with the >>> inclusion of flexible search criteria's like what I've listened above? >>> >>> 2)After allocating of those 10 sequences I should to find *all* >>> possible orthologs against them-> to find all OR in another species >>> which have conservation degree agains target above chosen threshold . >>> Eventually I'll need to make some analysis regarding overall >>> conservation of the OR within different species and to point out what >>> OR are most important (conservative) in different organisms in >>> accordance to some pre-defined criteria's (i.e ecological niche of the >>> species because it's associated with the odor spectrum). Here the >>> question: would the methods like blast-search be enough for such task >>> or I should to use also methods like evolutional tracing (to find all >>> important residues in different OR which tend to co-evaluate) based on >>> its alignment data? >>> >>> I would be very thankful for any suggestion! >>> >>> James >> >> >> _______________________________________________ >> Chimera-users mailing list >> Chimera-users(a)cgl.ucsf.edu >> http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users >

1 0

On the seeking of the orthologs
by James Starlight 22 May '15

22 May '15

Dear Chimera users! I'd like to ask some suggestion regarding my bioinformatician task and possible ways of its realisation using some bioinformatics methods implemented in Chimera. 1)I need to seek for some sequence database (here I don't know absolutely how to do it) of the olfactory receptors and to select 10 sequences of the Ors corresponded to 10 different mammalian species -> according to my idea it should be very different species associated with different ecological niches so its olfactory spectrum should vary dramatically. Here the question: in what DB such search should be better with the inclusion of flexible search criteria's like what I've listened above? 2)After allocating of those 10 sequences I should to find *all* possible orthologs against them-> to find all OR in another species which have conservation degree agains target above chosen threshold . Eventually I'll need to make some analysis regarding overall conservation of the OR within different species and to point out what OR are most important (conservative) in different organisms in accordance to some pre-defined criteria's (i.e ecological niche of the species because it's associated with the odor spectrum). Here the question: would the methods like blast-search be enough for such task or I should to use also methods like evolutional tracing (to find all important residues in different OR which tend to co-evaluate) based on its alignment data? I would be very thankful for any suggestion! James

2 1

clustering
by Rafael Rocha 22 May '15

22 May '15

hello, I have a set of ligands that share Some atoms in common. How I can clustering only those atoms? Also, how I can change the variation of RMSD for this cluster analysis in the "ensemble cluster" tool of chimera? Thanks, Rafael.

2 2

chimera improvment suggestion
by Michał Kadlof 21 May '15

21 May '15

I have small improvement suggestion: Is it possible to display directories in Open File dialog with different font than files, for example bold? I'm loosing a lot of time on scrolling my directories tree to find particular file to open. Highlighting of directories would be very handy. -- best Michał

2 1

Re: [Chimera-users] How can I calculate RMSD with chimera script
by Eric Pettersen 19 May '15

19 May '15

Since runCommand() doesn't return any value, you need to call the underlying Python rmsd() function, which isn't hard to do. Somewhere at the top of your script put: from Midas import rmsd and then replace: rc("rmsd #0:{0}-{1} #1:{0}-{1}".format(s,k)) with: val = rmsd("#0:{0}-{1}".format(s,k), "#1:{0}-{1}".format(s,k)) --Eric On May 19, 2015, at 12:45 PM, Michał Kadlof <ramidas(a)gmail.com> wrote: > Yes it work, but then I have to write and run another script to parse output from chimera which cost HDD access time. I was wonder, if it possible to do it in single script without dumping output to separate file. > > -- > Michał > > 2015-05-19 18:53 GMT+02:00 Eric Pettersen <pett(a)cgl.ucsf.edu>: > Hi Michał, > Redirecting output from a program into a file is pretty easy in any Unix/Linux shell. This would redirect Chimera's non-error output into a file named "chimera-output.txt": > > chimera chimera-arguments-you-normally-use > chimera-output.txt > > You can also redirect the error output, but the syntax for that depends on the shell you are using. You should be able to easily Google it if you need to do that. > > --Eric > > Eric Pettersen > UCSF Computer Graphics Lab > http://www.cgl.ucsf.edu > > On May 19, 2015, at 1:16 AM, Michał Kadlof <ramidas(a)gmail.com> wrote: > >> Hello, >> >> I have a need of calculating about 1000 of RMSDs between two particular regions of reference model and testing model. I would like to do this with chimera script. I'm using chimera-1.10.1-Headless Linux 64-bit and have no access to any GUI. >> >> Currently my script looks like this: >> >> import sys, os, re >> from chimera import runCommand as rc >> >> << not important stuff >> >> >> rc("open %s" % (native)) >> for i in models: >> rc("open " + i) >> rc("rmsd #0:{0}-{1} #1:{0}-{1}".format(s,k)) >> rc("close #1") >> rc("close all") >> rc("stop now") >> >> This allow me to see output on the screen, but I don't know how to capture it for further processing. I found somewhere something like this: >> >> from chimera.tkgui import saveReplyLog >> << do something >> >> saveReplyLog('replyLog.txt') >> >> But it seem to work only in versions with GUI. >> >> Please help. >> >> -- >> best whishes >> Michał Kadlof >> >> University of Warsaw >> _______________________________________________ >> Chimera-users mailing list >> Chimera-users(a)cgl.ucsf.edu >> http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users > > > >

2 1

Load AMBER parm/traj from command line.
by ros 19 May '15

19 May '15

Hello! Is it possible to start up chimera and have it load an AMBER topology and trajectory file? For example, in Linux, you can start VMD with: vmd -parm7 name.prmtop -netcdf traj.nc And will load the topology with the corresponding trajectory frames. In chimera we do that using the 'Get Ensemble info' window in the MD movie tools. Is it possible to do it in a single line? Thank you! Rodrigo.

2 2

How can I calculate RMSD with chimera script
by Michał Kadlof 19 May '15

19 May '15

Hello, I have a need of calculating about 1000 of RMSDs between two particular regions of reference model and testing model. I would like to do this with chimera script. I'm using chimera-1.10.1-Headless Linux 64-bit and have no access to any GUI. Currently my script looks like this: import sys, os, re from chimera import runCommand as rc << not important stuff >> rc("open %s" % (native)) for i in models: rc("open " + i) rc("rmsd #0:{0}-{1} #1:{0}-{1}".format(s,k)) rc("close #1") rc("close all") rc("stop now") This allow me to see output on the screen, but I don't know how to capture it for further processing. I found somewhere something like this: from chimera.tkgui import saveReplyLog << do something >> saveReplyLog('replyLog.txt') But it seem to work only in versions with GUI. Please help. -- best whishes Michał Kadlof University of Warsaw

2 1

depth cue by color for tomogram
by Michael Elbaum 18 May '15

18 May '15

Hello all, I noticed a new sort of 3D color glasses (http://chromatek.com/) that work via blue-red chromatic aberration rather than shifted pairs. I think this would be a useful addition if it isn't available already somehow. The major advantage for a lecture situation is that 2D images are visible also without glasses. I tried to encode depth by color in Chimera but was not successful. I'm using solid style with transparency to see "into" the tomogram, loaded in Volume Viewer. I'd like the deep features to be in one color and the higher parts another, with sparse regions relatively transparent. In the main Viewing panel there's an option for depth cue on the Effects tab, but the cueing is relative to the Side View orientation rather than the real depth. Also there seems to be no option for a color gradient. thanks for any help, Michael Prof Michael Elbaum Dept of Materials and Interfaces Weizmann Institute of Science Rehovot 76100, ISRAEL tel +972 (0)8-9343537 fax +972 (0)8-9344138

3 5

Select multiscale model surface pieces with command lines
by TUBIANA Thibault CNRS 18 May '15

18 May '15

Hi everyone! I Want to script a movie, but to do so, I need to select some surface pieces of a virus multiscale model with command line. I didn't find any way to do that. I try to select them with the mouse (ctrl+click) and then rename the selection with "namesel", it works but when I save the session and reload it, theses selections are cleared.... I can't do that each times I load my session because I have plenty of them and it's to generate the movie on SSH on another computer more powerfull... I can't split them to create a new model (for exemple with coordinates) because I will have to update all my "savepos" position... Did someone know how to do that ? Thank you a lot for your help, Regards, -- Tubiana Thibault PhD Student Institut de Biologie Intégrative de la Cellule CNRS UMR 9198 CEA-Saclay, Bâtiment 532 - Pce 12A 91191 Gif-sur-Yvette 06.79.90.95.85 http://www.tubiana.me

2 1

Problem running Chimera remotely on CentOS7 system
by "Kestilä, Veli-Pekka" 12 May '15

12 May '15

Hi, I have following error when trying to run Chimera remotely on CentOS7 system from OSX this seems to affect also Windows with Xming libGL error: dlopen /usr/lib64/dri/swrast_dri.so failed (/opt/sci/UCSF-Chimera64-1.10.1/lib/libgcc_s.so.1: version `GCC_4.7.0' not found (required by /usr/lib64/dri/swrast_dri.so)) libGL error: unable to load driver: swrast_dri.so libGL error: failed to load driver: swrast X Error of failed request: GLXBadContextTag Major opcode of failed request: 150 (GLX) Minor opcode of failed request: 5 (X_GLXMakeCurrent) Serial number of failed request: 1359 Current serial number in output stream: 1359 On windows the error is: ibGL error: dlopen /usr/lib64/dri/swrast_dri.so failed (/opt/sci/UCSF-Chimera64-1.10.1/lib/libgcc_s.so.1: version `GCC_4.7.0' not found (required by /usr/lib64/dri/swrast_dri.so)) libGL error: unable to load driver: swrast_dri.so libGL error: failed to load driver: swrast X Error of failed request: GLXBadCurrentWindow Major opcode of failed request: 145 (GLX) Minor opcode of failed request: 5 (X_GLXMakeCurrent) Serial number of failed request: 1350 Current serial number in output stream: 1350 How I can get this to work? Intended use would be in classroom where we don't have permissions to run or install Chimera. Greetings, Veli-Pekka Kestilä -- Veli-Pekka Kestilä / Huone 4318 / Biokeskus 3 (Viikinkaari 1) / 050-3000014

2 1

Acylation of peptide
by Remmer Salas 08 May '15

08 May '15

Hi! I am new to Chimera and I am trying to modify the PDB file of a peptide by acylation. I have also a PDB file for the fatty acid which I got by removing the bound protein. I used the join models to form a peptide bond but prompted with this error: The following problems occurred while reading PDB file for anomt1.pdbStart residue of secondary structure not found: HELIX 1 1 LEU 2 LEU 10 1 8Model 1 (anomt1.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 3Use command 'help ksdssp' for more information.TypeError Exception in Tk callback Function: <bound method BuildStructureDialog._addParamBond of <BuildStructure.gui.BuildStructureDialog object at 0x04F9FBF0>> (type: <type 'instancemethod'>) Args: ()Traceback (innermost last): File "C:\Program Files\Chimera 1.10.1\bin\lib\site-packages\Pmw\Pmw_1_3_3\lib\PmwBase.py", line 1747, in __call__ return apply(self.func, args) File "C:\Program Files\Chimera 1.10.1\share\BuildStructure\gui.py", line 313, in _addParamBond self._appbDihedral.get(), phi=self._appbPhi.get()) File "C:\Program Files\Chimera 1.10.1\share\BuildStructure\__init__.py", line 869, in cnPeptideBond cn[1].residue.phi = phi File "C:\Program Files\Chimera 1.10.1\share\chimera\phipsi.py", line 91, in setPhi _setAngle(bond, phi, getPhi(res), "phi") File "C:\Program Files\Chimera 1.10.1\share\chimera\phipsi.py", line 134, in _setAngle br.angle = (newAngle - curAngle, br.biggerSide())TypeError: unsupported operand type(s) for -: 'float' and 'NoneType' TypeError: unsupported operand type(s) for -: 'float' and 'NoneType' File "C:\Program Files\Chimera 1.10.1\share\chimera\phipsi.py", line 134, in _setAngle br.angle = (newAngle - curAngle, br.biggerSide()) My first question is how to solve this error. Second, how can I optimize the generated PDB file using Chimera? I want also to convert this optimize PDB file to PQR. Any advice is greatly appreciated.

2 2

clip planes
by Bharat Lakhani 08 May '15

08 May '15

Dear All I want to visualize the core of my protein, Is it possible to slice the protein right into two half in order to visualize the core of the protein. I tried using chimera pre-capping model but it has not been work out well. Any suggestion on this would be great. Thanks

4 3

Coloring selected residues based on B-factor
by Mahendra B Thapa 07 May '15

07 May '15

Dear Chimera users, While using Chimera to study a protein with 511 residues, I colored residues 1- 274 'dark grey', residues 274-417 'tan' and residues 418-511 'dark khaki'. I want to color residues, say, 313,350,412,415,485,503, from blue to red with an increasing order of B-factors, say respective B-factors are 40,13,2,21,33,6 without altering the color of other residues. I manually put these B-factors values for the corresponding residues; rest residues have no b-factor values as the structure is obtained from a simulation. I tried some unsuccessful attempts visiting some previous posts (sorry if I missed such posts close to my query). Any comments and suggestions might help me. Thank you, Mahendra Thapa University of Cincinnati,OH

2 1

aborting a volume command
by Dougherty, Matthew T 07 May '15

07 May '15

I issued this command, which works correctly. volume #2 planes z,1,100,1,10 I try to abort it around slice 50 using freeze, which does not terminate. Is there a way to abort it? Matthew Dougherty National Center for Macromolecular Imaging Baylor College of Medicine ===========================================================================

2 1

Displaying and calculating protein surface enclosed pockets
by Matthias Fellner 07 May '15

07 May '15

Good day I am a PhD student at the University of Otago New Zealand. First I would like to thank you for this very useful software, I use it for all my protein crystallography pictures and a lot of analysis as well, I always cite your citation for sure when I publish using your software. Normally I get very far with your tutorials but this time I failed to achieve my goal and I wanted to write a direct question. I have a protein with an active site pocket. I can display the surface of the protein and then I see the pocket inside but I am really interested in the enclosed surface = the active site pocket. How can I display the "negative" (enclosed) surface area and possibly calculate its volume. Can this be done with the pdb file or is such an analysis only possible based on the electron density. I learned how to manipulate the density map in chimera, so if this is the only way I should also be able to do it if you explain the details. If I get this done I wanted to ask if there is a way to ignore water molecules from the displaying of surfaces and later calculations or if I should simply work with a pdb which has the waters removed. Thank you Matthias Fellner

4 5

listen selection command
by Mungo Carstairs (Staff) 07 May '15

07 May '15

Hi, I am trying out a REST listener for selections in Chimera and need help. I can talk to Chimera over REST. I can start an Http listener in Jalview - and verify it in my browser. I have sent this command to Chimera: listen start selection rest http://134.36.162.131:64745/ Then silence. I haven't managed to provoke Chimera into sending me a SelectionChanged notification. I am making selections in the Chimera window - as I can see by the 'Inspect Selection' panel. I have also started Chimera standalone (no REST service) and tried the listen command in the command line, but yet to see any output on a selection. What have I missed or how can I diagnose what I am doing wrong here? Thanks, Mungo Carstairs Jalview Computational Scientist The Barton Group Division of Computational Biology College of Life Sciences University of Dundee, Dundee, Scotland, UK. www.jalview.org<http://www.jalview.org/> www.compbio.dundee.ac.uk<http://www.compbio.dundee.ac.uk/> The University of Dundee is a registered Scottish Charity, No: SC015096

3 5

Re: [Chimera-users] [chimera-dev] VolumeEraser
by Tom Goddard 06 May '15

06 May '15

Hi Alex, Surfnet in Chimera does not produce a volume, only surfaces. It calculates the volume but then computes surfaces and deletes the volume. If you would to measure the enclosed volume of the disconnected pieces produced by surfnet, first make those pieces independently selectable with command sop split #0 where #0 is the id for the surfnet model (use Favorites / Model Panel to find the correct id number). Then select the piece of interest (ctrl click on it), then get its enclosed volume with command measure volume sel To type commands to Chimera use menu Favorites / Command Line. This question belongs on the chimera-users mailing list, not the chimera-dev list which is for programming issues. Tom > On May 6, 2015, at 8:18 AM, Alex wrote: > > Dear Sirs, > I recently started using chimera, please apologise my lack of knowledge. > In order to erase a surfnet volume I am trying to use the Erase volume module, but it doesn't work. Are the surfnet objects available to be used in VolumeEraser functions? Is there another way to study volume intersections? > If you know, please point out the article which describes the procedure for this kind of volume calculation. > > Best regards, > Alex > _______________________________________________ > Chimera-dev mailing list > Chimera-dev(a)cgl.ucsf.edu <mailto:Chimera-dev@cgl.ucsf.edu> > http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-dev <http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-dev>

1 0

PC shut down automatically while using Chimera
by kzhang＠mrc-lmb.cam.ac.uk 06 May '15

06 May '15

Dear Chimera Developers, Recently, I found something really annoying when I was using Chimera. The problem is that when I use chimera to display CryoEM maps(400X400X400, 32 bits), my computer shut down automatically. I don't think there's something wrong with the hardware and coot or pymol is fine. This is not just an accident. It happened quite a few times, especially when I suddenly zoom in, computer is suddenly shut down. Later there would be some error massage related to the instability of the power, CPU self protection "£$%^&# ... when I turned on the computer again. But I don't think there is anything wrong with the power as well, I've run other softwares without any problem. Also, it happens to both Window 7 and CentOS 6.5 in the same PC. So it's obviously hardware level issue related to Chimera. My guess is that there might be something in Chimera that activated the self-protection of CPU or similar things. Any ideas why this occurs and how can I solve it? Best wishes Kai

3 2

imod
by Dougherty, Matthew T 04 May '15

04 May '15

I have recently been getting some imod segmentation files and have been noticing three problems 1) when I read the MRC and imod files I am having a scaling problem, and the datasets do not integrate. For example, if the MRC voxel size is 31.22 31.22 31.22, and the origin index is 0 -114 0, I have to change the the voxel size to 3.122 3.122 3.122 , and the origin index to 0 0 0 so that graphics match correctly. Since there are number of people in the data chain before it gets to me, I thought it might check to determine if this is known problem before going upstream to unknown people to figure out if they are getting the imod output settings correct. 2) The graphics performance with imod segmentation files is sluggish in chimera, not sure what it is in imod. I assume that the format defines the vertices at the highest resolution corresponding with the MRC file. Is there a way to improve the performance on the chimera side? I am not using contouring, just reading the file in default mode. 3) The segmentation has a lot of components, such as vesicles and microtubules, that have different colors. Is there a way to break apart the components so they can be manipulated seperately? E.g., display/hide, change colors, translate and rotate. Matthew Dougherty National Center for Macromolecular Imaging Baylor College of Medicine ===========================================================================

2 2

running script + recording movie
by Dmitry 04 May '15

04 May '15

Dear colleagues, Does anybody know how to run the script in chimera - *.cmd (from windows) and in the same time record everything that script does as a movie? Sincerely, Dmitry

2 1

chimera meeting tomorrow
by Greg Couch 29 Apr '15

29 Apr '15

The agenda is available in advance! Thank you Conrad! Conrad says "feel free to add items" :-) https://plato.cgl.ucsf.edu/trac/chimera2/wiki/2015-4-30

1 1

Auto dock Vina Increasing Exhaustiveness
by Tunc Morova 29 Apr '15

29 Apr '15

Hello! I am using UCSF Chimera, Autodock Vina extension for obtaining conformations of my protein and ligand. In this interface, As you know, I can change the exhaustiveness value between 1-8. However, since my system is more than 30x30x30 A, I need to get my exhaustiveness value more than 8 for getting more accurate results. But I can’t make it more than 8 through Auto dock Vina interface in Chimera. How can I increase this exhaustiveness value ? Btw, Instead of using Opal web service, I am using vina locally. Thank you many times for your help Best Regards Tunc Morova

2 1

Converting sdf file into mol2 format
by Mahendra B Thapa 29 Apr '15

29 Apr '15

Dear Chimera users, I am trying to convert sdf file ( downloaded from* diversity set* data from 'http://scs.illinois.edu/htsf/compound_collection/nci.php') into mol2 format. When I ran the command chimera --nogui processSDF.py as given the previous posting ( http://www.cgl.ucsf.edu/pipermail/chimera-users/2010-March/004974.html) I got the following error messages: Cannot find consistent set of bond orders for ring system containing atom #0.73:1@N1 --- Unknown redo atom type: N3 Traceback (most recent call last): Any suggestions for correcting these errors will be a good help for me. Thank you, Mahendra Thapa University of Cincinnati,OH

3 2

Create an artificial bond between ribbons for 3D printing
by David Bhella 28 Apr '15

28 Apr '15

Dear Chimera users, I am trying to prepare a ribbon diagram for 3D printing. As the molecule is a trimer, I would like to create a join between the ribbons such that the model is sent to the printer as a single part. So, I thought that I might be able to create a bond, but I suspect that the current tools may not be able to achieve this. Ideally I don’t want to have to show side chains, moreover the bond needs to be preserved when exported as STL format and finally it will be quite long. The only tool that would seem to come close to my requirement is the distance measure, however I can’t see how to convert that pseudo bond to a cylinder that will be exported. Any suggestions would be gratefully received. Thanks, D. Dr David Bhella MRC-University of Glasgow Centre for Virus Research Garscube Campus 464 Bearsden Road Glasgow G61 1QH Scotland (UK) Telephone: 0141-330-3685 Skype: d.bhella Virus structure group on Facebook: https://www.facebook.com/CVRstructure Molecular Machines - Images from Virus Research: http://www.molecularmachines.org.uk CVR website: http://www.cvr.ac.uk CVR on Facebook: https://www.facebook.com/centreforvirusresearch

3 3

Delete Bug?
by Korbin West 27 Apr '15

27 Apr '15

Hi all, I'm running some code on Chimera to select only certain chains and nearby aspects, but when I try to delete the extraneous information, Chimera freezes up. I'll detail the commands I use below, and for reference I'm interested in the H chain of 1CT1, but it doesn't seem to work with any of its chains. open 1CT1 zonesel :.h 3.5 extend true select invert del sel It then freezes on del sel. I've tried this with other pdb files and it works fine, but for some reason it stops on this one. For example, it works fine on 1PZ3 for its A and B chains. Is this some kind of bug or is there a work-around to this issue? We looked at the residues and couldn't find anything out of place, so we are not sure what is happening. Much thanks in advance. Best, Korbin West? Wabash College

3 3

Making a composite residue
by Repic Matej 23 Apr '15

23 Apr '15

Hello, I wonder if it is possible to construct a composite residue from two existing residues in Chimera? What I would like to do is concatenate residues ALA100 (mid chain) and GLY200 (end of chain) to make a composite residue AGL100. Chimera won't let me set the same residue number via the "resrenumber" command. Preferably, the composite residue would have all the atoms listed where ALA100 was in the original pdb. While I realize there are trivial text editor-based workarounds I would like to know if there is a native way. Thank you all, ---------------------------------------- Matej Repic Ecole Polytechnique Fédérale de Lausanne SB - ISIC - LCBC BCH 4108 CH - 1015 Lausanne ----------------------------------------

2 1

about volume viewer
by 张雷 21 Apr '15

21 Apr '15

Hi, the experts of Chimera, I am a new user of Chimera and have a question to consult you. When I display a 3D model in the volume viewer, I found that there was an original value of the level, as the attached figure shown. If I increased and decreased the values of level, the densities of 3D models were different. Therefore, how to deternime the right value of level for a protein. It is known that the average density of proteins is 1.35 g/cm3; If a protein has a molecular weight of 1350 kDa, and 1 Da = 1.66054E-27 kg; therefore, the protein volume is 1661 nm3.If I want to display the volume of 1661 nm3 in chimera, how can I decide the value of level? In another words, how can I calculate the value of level from the actual volume of a protein? Thank you! Best wishes! Lei Zhang Peking University, Beijing, China

3 3

Re: [Chimera-users] coloring individual segments of a protein
by Greg Couch 21 Apr '15

21 Apr '15

To get a different resolution than 72 dpi, you need to turn on the "Use print units" option in the Save Image dialog, and set the "Print resolution" and pick the image width and height. Or you can use a program like Adobe PhotoShop or GIMP to edit the DPI if you already have an image that you like. HTH, Greg On 04/21/2015 03:01 PM, A G Szabo wrote: > Greg > > Thanks for your reply > > > I had just sent Elaine another question that perhaps you might be able to > answer. > > I am currently writing a book and the editor would like the protein > structure images saved with at 300 dpi. When I saved an image of > Ribonuclease, PDB 5RSA, both as a png file and a Jpeg file the editor > indicated that the resolution was not sufficient and was only 72 dpi. > > Suggestions? > > Thanks > > ags

1 0

Residue Numbers vs. Residue Codes
by Korbin West 21 Apr '15

21 Apr '15

Hi all, I'm trying to determine uncommon ligands in pdb files using a python code, but have run into an odd dilemma. In my code, I open a pdb in Chimera, select the ligand, and write out the residue list. However, the formatting of this list changes depending on how I run the program. If I run the program through xterm, it will print out only residue numbers. If I run the program though Chimera's Idle, however, I will get the residue codes and the residue numbers. For example, I ran this with 3VKG, you can see the results in the attached files, where I show a sample code of what I'm using, the xterm results and the Chimera results. I want to do this for a large number of pdbs, so I'd like to do this --nogui through xterm, but I need the residue codes. Has anyone else ran into this problem or know of a work-around? Thanks in advance. Best, Korbin West Wabash College

2 1

coloring individual segments of a protein
by A G Szabo 21 Apr '15

21 Apr '15

How do I use Chimera to color a protein in the following way: for example: Helix 1 - red Helix 2 - yellow Helix 3 - magenta Strand 1 - dark blue Strand 2 - cyan Etc. A G Szabo

3 7

keyboard accelerators
by Dougherty, Matthew T 20 Apr '15

20 Apr '15

Hi, had a few questions about the documentation. A click in the graphics window may be required. Sometimes it is necessary to click in the graphics window to take back focus from some other dialog so that accelerator or command keystrokes will register. >> is there a way to send accelerators through the Viewing-Side View? I have been dedicating one monitor for the 3D graphics and the other for control, menus, etc. using various USB devices sending accelerators. Would prefer not put to the focus/mouse cursor on the 3D graphics. Accelerators could be shown on menu entries. The keyboard shortcuts could be shown on the corresponding menu entries to make them easier to learn.? >> How? Matthew Dougherty National Center for Macromolecular Imaging Baylor College of Medicine ===========================================================================

3 3

Help in MD movie
by Ankit Agrawal 18 Apr '15

18 Apr '15

Hi I have trajectory file from gromacs and opened in chimera. Now I am doing inertia calculation. I wrote a python script to show the ellipsoid. But this script is showing only the ellipsoid. I want ellipsoid as well as the molecular structure of that frame, i.e. both ellipsoid and the molecule. so that I can show the how molecule is changing its position. This is like a structure overlay but I want it with ellipsoid. thanks regards Ankit Agrawal

2 1

Reducing the map to a cubic box
by Sushree Tripathy 17 Apr '15

17 Apr '15

Hi, I went through the chimera user's guide and as per the syntax used "vop boxes #0 #0 useMarkerSize true" or "vop boxes #0 #0 size 36" (the simplest case) to create a cubic box from my density map (54, 36, 51). It is neither giving me any error nor giving me the output map. Please help. Thanks Sushree. Sushree Tripathy Graduate Student Department of Physics University at Buffalo

2 1

Generating a map from atomic molecule
by Tom Goddard 17 Apr '15

17 Apr '15

Hi Zhihai, The reason your map surface shows the triangulation is because the grid spacing is large. You made the map with molmap with 30A resolution. By default molmap sets the grid spacing to 1/3 of the requested resolution, so 10 Angstroms in this case. Apparently this is too coarse for your needs. So spedify a finer grid spacing in the molmap command. For instance molmap #0 30 gridSpacing 3 to produce 30 A resolution with 3 A grid spacing. Here is the molmap documentation http://www.cgl.ucsf.edu/chimera/current/docs/UsersGuide/midas/molmap.html <http://www.cgl.ucsf.edu/chimera/current/docs/UsersGuide/midas/molmap.html> Tom > On Apr 16, 2015, at 8:45 PM, 李智海 wrote: > > Hi Tom, > It seems that the low-resolution (about 30A) map generated from atomic model by “molmap" has many artifacts on its surface (see first attached figure), and I can see many rectangle signatures of it in Fourier space (see 2nd attached figure) so that I can not use it as the start map for further reconstruction. > Please tell me how can I get rid of this. Thank you! > > Best wishes! > > > Zhihai Li > > <molmap_30A.png><FFT.png> > > Zhihai Li > > Ph.D student in State Key Laboratory of Molecular Vaccinology and Molecular Diagnostics. > School of Life Sciences, Xiamen University. > Xiang'an Road(South), Xiang'an District, Xiamen City, > Fujian Province, China. > >> On Apr 15, 2015, at 1:28 AM, Tom Goddard wrote: >> >> Hi Zhihai, >> >> If you want the map computed from an atomic model but on a grid that matches a reference map then you should compute the map with molmap and then resample it on the reference map grid. If #0 is the atomic model and #1 the reference map. >> >> molmap #0 8 edgePadding 50 >> vop resample #0 onGrid #1 >> >> Here the 8 is the resolution for the computed map and 50 Angstroms controls how far beyond the molecule bounding to compute the map. >> >> Tom >> >> >> >>> On Apr 13, 2015, at 9:00 PM, 李智海 wrote: >>> >>> hi, >>> here I want to generate a density map from a atomic model in Chimera, then convert the generating map into .pif file to do as the initial map performing in Auto3dem. >>> I have met some questions to successfully finish these work. >>> >>> 1.How can I generate a mol-map from a fitting model with identical dimensions, center coordinates and pixel size of its reference map. I have tried to use “Fit in map”, but the generating map differed from its fitted map in the map size. >>> >>> 2. The mrc map created in Chimera has dimension of even number of columns. If I set the column number to odd number, is it ok (for example, if the origin is off by half pixel?). The reason I am asking such a question is that I will convert this mol-map to pif format, and pif format map has dimension of odd number columns. >>> >>> >>> Thanks in advanced! >>> >>> Best wishes! >>> >>> Zhihai Li >>> >>> >>> >>> _______________________________________________ >>> Chimera-users mailing list >>> Chimera-users(a)cgl.ucsf.edu <mailto:Chimera-users@cgl.ucsf.edu> >>> http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users <http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users> >> >

1 0

on script which can be read and executed by chimera
by Smith Liu 17 Apr '15

17 Apr '15

Dear All, Will you please introduce to me a web link on how to compile and edit scripts which can be read and executed by Chimera? I am looking forward to getting your reply. Smith

2 1

the "select" funtion of mouse is disable
by sunyeping 16 Apr '15

16 Apr '15

Dear all,I find the "select" funtion of mouse in Chimera is disable. I can not use "ctrl+left mouse button" to select any item in the structure. Have you met this problem before? My systems is windows 7 and the Chimera version is 1.10.1. How to solve this problem? Do I need to change some setting?Best regards,Yeping Sun Institute of Microbiology, Chinese Academy of Sciences

2 2

Re: [Chimera-users] Reducing the map to a cubic box
by Tom Goddard 15 Apr '15

15 Apr '15

In Chimera “atoms” means an atomic model, for example from a PDB file, not a density map. If you want to reduce your map to a cube by cropping and it is size 54,36,51 use command volume #0 region 0,0,0,35,35,35 That gets the lower left corner. To get the center part volume #0 region 9,0,7,44,35,42 The 6 numbers are imin,jmin,kmin,imax,jmax,kmax where i,j,k are indices along the 3 map axes. You can also do this with the Volume Viewer dialog menu entry Features / Region Bounds, or interactively with the mouse using Features / Subregion Selection. For more info about map manipulations http://www.cgl.ucsf.edu/chimera/data/tutorials/volumetour/volumetour.html Tom > On Apr 15, 2015, at 1:17 PM, Sushree Tripathy wrote: > > Hi, > > Thanks for your prompt reply. What I am trying to do is I have a map which I have loaded in chimera as #0. Now I want to reduce it to cubic size. Since the syntax is "vop boxes <> volume-spec atom-spec <https://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/frameatom_spec.html> size d", so I am choosing #0 for my volume and all the atoms in that volume. so #0 again. Please correct me if I am wrong. > > > Thanks > Sushree. > > Sushree Tripathy > Graduate Student > Department of Physics > University at Buffalo > > On Wed, Apr 15, 2015 at 4:13 PM, Tom Goddard wrote: > Hi Sushree, > > Why does your "vop boxes” command use model #0 for both the map and the markers? That can’t be right. If your marker(s) are in model #1 you would need “vop boxes #0 #1 size 36”. What are you trying to do? > > http://www.cgl.ucsf.edu/chimera/current/docs/UsersGuide/midas/vop.html#boxes <http://www.cgl.ucsf.edu/chimera/current/docs/UsersGuide/midas/vop.html#boxes> > > Tom > > > >> On Apr 15, 2015, at 1:03 PM, Sushree Tripathy wrote: >> >> Hi, >> >> I went through the chimera user's guide and as per the syntax used "vop boxes #0 #0 useMarkerSize true" or "vop boxes #0 #0 size 36" (the simplest case) to create a cubic box from my density map (54, 36, 51). It is neither giving me any error nor giving me the output map. >> >> Please help. >> Thanks >> Sushree. >> >> >> Sushree Tripathy >> Graduate Student >> Department of Physics >> University at Buffalo >> _______________________________________________ >> Chimera-users mailing list >> Chimera-users(a)cgl.ucsf.edu <mailto:Chimera-users@cgl.ucsf.edu> >> http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users <http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users> > >

1 0

Reducing the map to a cubic box
by Sushree Tripathy 15 Apr '15

15 Apr '15

Hi, I went through the chimera user's guide and as per the syntax used "vop boxes #0 #0 useMarkerSize true" or "vop boxes #0 #0 size 36" (the simplest case) to create a cubic box from my density map (54, 36, 51). It is neither giving me any error nor giving me the output map. Please help. Thanks Sushree. Sushree Tripathy Graduate Student Department of Physics University at Buffalo

2 1

generating a map from a atomic model with the same size of its reference map
by 李智海 14 Apr '15

14 Apr '15

hi, here I want to generate a density map from a atomic model in Chimera, then convert the generating map into .pif file to do as the initial map performing in Auto3dem. I have met some questions to successfully finish these work. 1.How can I generate a mol-map from a fitting model with identical dimensions, center coordinates and pixel size of its reference map. I have tried to use “Fit in map”, but the generating map differed from its fitted map in the map size. 2. The mrc map created in Chimera has dimension of even number of columns. If I set the column number to odd number, is it ok (for example, if the origin is off by half pixel?). The reason I am asking such a question is that I will convert this mol-map to pif format, and pif format map has dimension of odd number columns. Thanks in advanced! Best wishes! Zhihai Li

2 1

generating a map from a atomic model with the same size of its reference map
by 李智海 14 Apr '15

14 Apr '15

hi, here I want to generate a density map from a atomic model in Chimera, then convert the generating map into .pif file to do as the initial map performing in Auto3dem. I have met some questions to successfully finish these work. 1.How can I generate a mol-map from a fitting model with identical dimensions, center coordinates and pixel size of its reference map. I have tried to use “Fit in map”, but the generating map differed from its fitted map in the map size. 2. The mrc map created in Chimera has dimension of even number of columns. If I set the column number to odd number, is it ok (for example, if the origin is off by half pixel?). The reason I am asking such a question is that I will convert this mol-map to pif format, and pif format map has dimension of odd number columns. Thanks in advanced! Best wishes! Zhihai Li

1 0

Fitting pdb model into difference map
by Chinaza Egbuta 13 Apr '15

13 Apr '15

I generated an Fo-Fc difference map using CCP4i and was able to superimpose the PDB model structure onto the map in Coot but I am unable to do so with Chimera. The model and map are offset/rotated away from each other and I have not been able to fit the model into the map with Volume Viewer --> Tools --> Fit in Map. Please let me know how I can get around this issue. Thank you. CE

3 4

Extracting masked map region using selected sop piece?
by Oliver Clarke 13 Apr '15

13 Apr '15

Hi all, I find the sop command (particularly 'sop split' and 'sop hideDust metric volume rank') to be very useful for viewing a molecular envelope at a low threshold in an otherwise noisy map. I would like to save the volume enclosed by a single sop piece as a separate volume, but it does not seem like this is possible - when I select a sop piece and use 'volume nameregion’ to assign it a name, and then attempt to write just that region to an mrc, it saves the whole map. Is there some way to do this that I am missing? I realize that a similar outcome can be achieved with segger, but for larger maps this can be very slow (and potentially involve multiple cycles of selection and regrouping), whereas sop operations are almost instantaneous, so I feel it would be a valuable adjunct capability to that offered by segger. Also, the behaviour of sop split seems a little counterintuitive to me. When I use 'sop split' on a volume, it generates a new model with a different ID for the split map - fine so far. However, the original map seems to be completely invisible - it is listed as displayed in the model panel, but if I undisplay the split map in the model panel and toggle display of the original map in the model panel, nothing happens. If I go to the volume viewer, the original map is listed as displayed, and if I toggle it off and then on again there, it reappears. It seems like it would be maybe simpler for the user to just automatically undisplay the original map upon creation of the sop split map, but not to prevent the user from redisplaying it in the model panel. Here is a link to a screen recording demonstrating what I mean: https://www.dropbox.com/s/g40u9dhj0xslcbf/sop_split.mov?dl=0 <https://www.dropbox.com/s/g40u9dhj0xslcbf/sop_split.mov?dl=0> Best, Oliver.

3 3

Fitmap on overlap/correlation only within bounds of smallest map
by Oliver Clarke 13 Apr '15

13 Apr '15

Hi, I am trying to fit a set of large (whole protein) maps to a map corresponding to a single domain of the protein (extracted using ‘vop zone minimalbounds true’ with an atomic model of the domain in question. It works fine, except that it takes rather a long time (~2000 steps), and the final correlation is 0.33, which is strange because the overall correlation between the two complete maps is ~0.95, and I can see that the final overlap between the aligned maps is good. It seems like ‘fitmap’ is calculating the correlation/overlap based on all grid points in the map that I am aligning, which has a much larger box size than the map to which I am aligning it. Is there any way to restrict the metric analysis to just consider points within the smaller box of the reference map? I cannot find such an option in the documentation for fitmap. The use case for this is when one has multiple conformations of a protein (for example, EM volumes after 3D-classsification) and would like to align them on a single reference domain. Best, Oliver.

2 2

Ways to minimize data loss during resampling?
by Oliver Clarke 13 Apr '15

13 Apr '15

When one has multiple maps corresponding to the same protein in different conformations, it is convenient to save aligned copies of each map, for easy comparison later or for transfer to other programs. This requires interpolating the map to be aligned on the grid of the map to which one is aligning, e.g. using “vop resample x ongrid y” in Chimera. This step involves a certain degree of data loss, which is mostly insignificant but occasionally affects the quality of the resultant map in a way that is visually noticeable. Does anyone have any tips for minimizing this data loss, or is that impossible? E.g., if I oversample the map to which I am aligning by lets say 4 times, resample my map on that grid, and then save at gridstep 3 will that help, or is the data loss an unalterable consequence of the interpolation procedure? Best, Oliver.

2 1

Re: [Chimera-users] [chimera-dev] Using Chimera to display family shuffling information.
by Elaine Meng 10 Apr '15

10 Apr '15

On Apr 8, 2015, at 5:11 PM, Matthew Tiffany <matt.tiffany(a)icloud.com> wrote: > To Whom It May Concern: > I am interested in displaying results from DNA shuffling experiments on a crystal structure. Here is a graphical tool that I have found to do this: http://qpmf.rx.umaryland.edu/XoverDemo.html > > What I would like to do is input an alignment of parental sequences and separately a chimeric protein that is derived from the parental sequences by domain shuffling. Then color code a related structure based upon which domain of a parent the chimeric protein came from. The Xover program does this in a graph, but I would like to be able to do this in 3D space. It is similar to color coding by conservation in MultiAlign, but not exactly. Is there any existing plugin that you may be aware of that does this already or a potential backbone plugin that I can use as a starting point to try and code it myself. > > Regards, > Matthew Tiffany > Postdoctoral Fellow > Mark A. Kay Laboratory > Stanford University Hi Matthew, I’ve CC’d the chimera-users list since my answer doesn’t involve programming... Chimera itself includes a lot of sequence-structure interactivity and related features. Certainly you could do it semi-manually: (1) open your 3D structure and sequence alignment in Chimera. You can just the “File… Open” dailog. Chimera understands several sequence alignment formats. <http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/filetypes.html#alignment> <http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/multalignviewer/f…> (2) associate the structure with a sequence in the alignment. I don’t know exactly what your data look like, and so this requires some judgment calls and deciding amongst different strategies. My thought is at least initially, don’t worry about which sequence is associated with the structure, just try to get it associated with any sequence so that the 3D mapping along the alignment as a whole is correct. — approach A: just use sequence window menu Structure… Associations to force association with any one of the sequences, if it didn’t already associate automatically — approach B: add the sequence of the chimeric protein to your sequence alignment, sequence window menu Edit… Add Sequence, use From Structure (if it adds it poorly you may have to go through a few cycles of Edit… Delete Sequences/Gaps to remove it and then adding it back with different alignment parameters) … I’m leaning toward approach B as the better approach, but you can see whether one or the other works better for your situation… <http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/multalignviewer/m…> (3) then you can select the different parts of the structure for coloring by using the mouse in the sequence window to select the appropriate ranges of residues in the structure-associated sequence <http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/multalignviewer/m…> Another approach to #1-2 above would be to start with a sequence of just your structure (for example, open the structure, use Favorites… Sequence to show its sequence) and Blast against the PDB, then load in the Blast alignment. However, you already have your own alignment so I don’t think it would get you much further. Further, there is a newish “mda” command that performs multiple steps including blasting against PDB, automatically loading in the structures and Blast alignment, and conveniently laying out the structure hits N->C (based on alignment to query) in the 3D display. This also wouldn’t get you much further toward the color mapping of your chimeric structure, but the layout stuff is cool. It might not work if your shuffled cassettes are pretty short though. <http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/mda.html> I hope this helps, Elaine ----- Elaine C. Meng, Ph.D. UCSF Computer Graphics Lab (Chimera team) and Babbitt Lab Department of Pharmaceutical Chemistry University of California, San Francisco

1 0

Atomspec balloon don't show up
by LINA MARCELA PEDRAZA CASTILLO 09 Apr '15

09 Apr '15

Hi. (Sorry about my english, it's not my language) Dear users: I'm working on my thesis, and I use Chimera ucsf for everything, and suddenly, one day, it stops showing the atomspec balloons, (when the cursor is paused over an object (atom, bond, surface, *etc.*), whether to show the corresponding label information in a balloon in the graphics window, and PDB chain information, if available, in the status line <http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/chimerawindow.html#statusli…>.) OK, That's what's no happening anymore, and I need it, so... Can you help me please with this issue? I already read the docs and I did the modifications at preferences menu, and modify: Labels *Label font* (*Sans Serif* - *Normal -* *16* *Labels on top* (*true* *Show atomspec balloon* (*true* *Atomspec display style* (*simple* Please, help me with this. Thanks. I will appreciate it.

2 1

Mapping predicted secondary structure and solvent accessibility on structure or homology model?
by Oliver Clarke 09 Apr '15

09 Apr '15

Hi, I noticed that in the multalign sequence viewer, Chimera has the option to show predicted secondary structure, but it doesn’t seem to work, at least for my sequence (it seems to work initially, but when I show predicted helices/strands there is nothing to see). I think that incorporating modern secondary structure, solvent accessibility and disorder prediction (e.g. using JPRED, which I believe has an open API) would be a very valuable addition to Chimera for structural analysis - it would be particularly helpful for evaluating homology models and models that are based on low resolution crystallographic or CryoEM data, to pick up places where there might be a register error. Maybe this is something to consider for Chimera 2 (or maybe there is a way to load such data in now, if I could create it externally and load it in as a set of attributes…). Best, Oliver.

2 1

openGL error.
by Mosavverul Hassan 08 Apr '15

08 Apr '15

Hi, I am trying to get chimera run remotely on a linux server while logging into it from a mac terminal. I get the following error libGL error: failed to load driver: swrast X Error of failed request: GLXBadContextTag Major opcode of failed request: 150 (GLX) Minor opcode of failed request: 5 (X_GLXMakeCurrent) Serial number of failed request: 1386 Current serial number in output stream: 1386 glxinfo gives me the following name of display: localhost:13.0 libGL: OpenDriver: trying /usr/lib/x86_64-linux-gnu/dri/tls/swrast_dri.so libGL: OpenDriver: trying /usr/lib/x86_64-linux-gnu/dri/swrast_dri.so libGL: driver does not expose __driDriverGetExtensions_swrast(): /usr/lib/x86_64-linux-gnu/dri/swrast_dri.so: undefined symbol: __driDriverGetExtensions_swrast libGL: Can't open configuration file /mnt/raidc2/mhassan/.drirc: No such file or directory. libGL: Can't open configuration file /mnt/raidc2/mhassan/.drirc: No such file or directory. libGL error: failed to load driver: swrast display: localhost:13 screen: 0 direct rendering: No (If you want to find out why, try setting LIBGL_DEBUG=verbose) server glx vendor string: SGI server glx version string: 1.4 It seems to run fine from a linux computer but has problems with a mac. I am using X11 forwarding when I log in. I have X11 on my mac. Could someone please help me out? Regards, Mosavverul Hassan.

2 1

problems to start chimera
by torquato.bioq＠epm.br 08 Apr '15

08 Apr '15

Dear collegues !! I'm installed chimera-1.10.1-linux_x86_64.bin on my machine witch controlled by openSUSE 13.2. The installation procedure was fine, but when I try start chimera I got this messege: ricardo@ricardolinux:~$ chimera libGL error: unable to load driver: swrast_dri.so libGL error: failed to load driver: swrast X Error of failed request: GLXBadContextTag Major opcode of failed request: 155 (GLX) Minor opcode of failed request: 5 (X_GLXMakeCurrent) Serial number of failed request: 1457 Current serial number in output stream: 1457 Please can some one help T All the best Ricardo

2 1

plane normal axis
by Elaine Meng 07 Apr '15

07 Apr '15

Hi Oliver, If you have a plane defined in Axes/Planes/Centroids, you can write out its information including the normal (Save button on that tool), then place atoms at either end allowing you to define the axis. There are a few steps... In my little test, the plane output information included: Planes plane name, ID, center, normal, radius plane: p1 ( 19.023, 17.733, 7.080) (-0.380, 0.757, 0.531) 24.199 Then I used Tools… Build Structure, Start Structure, add atom at specific coordinates to add two atoms, one at the center of that plane and one to show the normal. The latter's coordinates are (center + normal) coordinates, and you can multiply the normal coordinates by some constant to make the axis longer. For example, I multiplied the normal coordinates by 10 before adding them to the center coordinates to give the position of the second atom as: x 15.223 y 25.303 z 12.390 Then I selected the two new atoms and used them to define an axis in Axes/Planes/Centroids, see attached image. I hope this helps, Elaine ---------- Elaine C. Meng, Ph.D. UCSF Computer Graphics Lab (Chimera team) and Babbitt Lab Department of Pharmaceutical Chemistry University of California, San Francisco On Apr 6, 2015, at 5:01 PM, Oliver Clarke <olibclarke(a)gmail.com> wrote: > Alternatively, would it be possible (or is there some way already) to define a separate axis that is the normal of a plane? In this case I’m working with an ion channel, so I’m looking for a consistent way to represent the change in orientation of particular extramembrane domains with respect to the channel axis. If I could define an axis for the membrane plane normal, I could at least measure the absolute orientation of the principal axes of particular domains or individual helices with respect to the channel axis. It seems like the plane normal is computed already (because there is an option in axes/planes/centroids) to align screen x/y/z to it, but I can’t see a way to generate it as a separate axis that I can use for measuring angles etc. > > Thanks for all your help as always! > > Oliver.

3 4

where is "Save each fit when doing multiple alignments" button
by sunyeping 07 Apr '15

07 Apr '15

Dear all,I am studying the chimera tutorial "fitting to segments" of Segger segmentation tutorials from Greg Pintilie at NCMI (http://ncmi.bcm.edu/ncmi/software/segger/docs_fitting) In the last paragraph, the tutorial states that to align the structure to multiple regions, one needs to check the "Save each fit when doing multiple alignments" button. But where is this button? I cannot find it in the "Fit to Segments", "Segment Map" or other dialogue boxes. Could you tell me where I can find it?Best regard,Yeping Sun Institute of Microbiology, Chinese Academy of Sciences

4 6

how to measure density diameter by chimera
by Smith Liu 07 Apr '15

07 Apr '15

Dear All, Will you please tell me the command or process to measure the diameter of the molecule density in the map by chimera? Smith

2 1

DTI stereo3D
by Dougherty, Matthew T 06 Apr '15

06 Apr '15

?When selecting DTI stereo3D, L&R images are placed side by side, and horizontally rescaled by half, known as SBS-half. I would like to generate SBS-full, ie no horizontal scaling. Is there a way to do that in chimera? Would prefer not to generate stereo pairs and rescale them back to 100% using a video editor. Matthew Dougherty National Center for Macromolecular Imaging Baylor College of Medicine ===========================================================================

2 2

segname
by Francesco Pietra 04 Apr '15

04 Apr '15

Hello: May I ask whether there is till a problem in having CHIMERA understanding segname? For those working with XPLOR (NAMD, for example) and multichain proteins, rendering results with CHIMERA becomes very tedious (having to use atom indexes) Thanks francesco pietra

2 6

Autodock vine problem
by Tunc Morova 04 Apr '15

04 Apr '15

Hello, I was trying to use built-in auto dock vina tool of chimera. I came a cross with an error saying service local`c:\”program filex(x86)”\”the scripps research Institute\Vina\Vina.exe` is unavailable. I have been checking forums for a longtime. Unfortunately, i couldn’t run the docking. couple days ago, I think i accidentally ran the docking locally but i am not sure. Is my pathing correct? how can my local auto dock vina be unavaible ? Please help me for finishing my thesis :D Thank you for your support -Tunc

3 3

Opening Gromacs .trr: error with "bham_lr"
by Alan Manning 03 Apr '15

03 Apr '15

Hi I want to import a Gromacs trajectory (a .trr file) in order to create an animation. I'm using the "MD Movie" tool, and I specify the input .tpr and .trr files. However, I get the following error when I try and make the movie: Don't know correct parameters for 'bham_lr' function Does anyone know what this is from? I'm using Gromacs 5.0.4. The complete Reply Log is as follows: using memory mapping VERSION 5.0.4 using floats version 100, generation 26 2613 atoms box o water Don't know correct parameters for 'bham_lr' function Thanks! Regards, Alan Manning

2 3

drawing ribbon
by Hartzell Jr., Criss 03 Apr '15

03 Apr '15

I have a protein with 10 transmembrane domains and would like to make a movie in which the ribbon is drawn in time from the N- to the C- terminus. I wrote a script that does that. But, the 10 transmembrane domains become hard see because they overlap one another. To solve this problem, I wanted to make the first 600 amino acids become 70% transparent after they had been drawn. Then the drawing would continue with the second half of the protein with no transparency. The problem is that the script that I have written to do this draws the first 600 amino acids 70% transparent immediately as they are drawn. Then the second half of the protein is drawn with 0 transparency. I set the wait time to 600 because I want to execute the transparecny command after the first 600 amino acids were drawn. But, instead, the drawing of the first 600 amino acids is tremendously slowed down by this command as if a wait of some duration occurs after each amino acid rather than after all 600 amino acids. I would appreciate your help. Thanks. movie record supersample 4; perframe "ribbon #2:$1" range 1,600; wait 600; transparency 70,r #2:1-600; perframe "ribbon #2:$1" range 601,959; movie stop; movie encode “d:\Dropbox\Chimera_files\movie1.mp4”; Criss Hartzell Professor Department of Cell Biology Emory University School of Medicine Atlanta, GA 30322 404-242-5719 (cell) For Fedex, UPS, other courier: Criss Hartzell 615 Michael St. 535 Whitehead Building Emory University Atlanta, GA 30322 ________________________________ This e-mail message (including any attachments) is for the sole use of the intended recipient(s) and may contain confidential and privileged information. If the reader of this message is not the intended recipient, you are hereby notified that any dissemination, distribution or copying of this message (including any attachments) is strictly prohibited. If you have received this message in error, please contact the sender by reply e-mail message and destroy all copies of the original message (including attachments).

2 2

New molecule from selection?
by Oliver Clarke 03 Apr '15

03 Apr '15

Hi, This may already be implemented, but I couldn’t find it - the nearest equivalent feature I know of is copy/combine, but that only works on whole models, not taking into account selections. I was wondering if it would be possible to implement a feature creating a new model comprised of all atoms in the current selection? This would be particularly helpful when combining many parts of multiple models, or when creating a separate model consisting of a subselection of the current model. This can be done at present using a combination of copy/combine and selecting and deleting those parts of the copied molecules that are not in the selection, but this gets messy for complicated selections. Best, Oliver.

3 3

select by color (or other display attributes)
by Oliver Clarke 02 Apr '15

02 Apr '15

Hi, one feature I would love to see in a future version of chimera is the capacity to select residues/atoms by color - this would be magic for modifying figures. Example use case: I’ve made a figure in which some complicated selection or set of selections is colored medium blue, and another selection is colored red. I’ve saved the session for the figure, but I (stupidly)didn’t save the corresponding selections, and now I have to remember exactly what selections I used. If there was a “select by color” function, I could simply select all the atoms that are colored medium blue (or some other color ID), and change them to another color, without having to know which residues were in the original selection. Analagously, it would also be handy, although less frequently useful, to be able to select on the basis of whether a residue is displayed, or represented in a particular manner - in particular, being able to restrict selection to displayed atoms/residues would be useful to avoid inadvertent modification of undisplayed models or parts thereof. Cheers, Oliver.

3 4

Smoothly change level of density map in animation?
by Oliver Clarke 02 Apr '15

02 Apr '15

Is there any way to alter the density level in a gradual fashion during an animation? I would like to first show the overall envelope of a cryoEM map, and then smoothly change to a higher contour level where I can see helices etc. Many other parameters, such as the presence of a MSMS surface or the position of a clipping plane, seem to automatically gradually transition between two key frames of an animation, but density level doesn’t seem to do so. Oliver.

2 2

mclip suggestion - match to arbitrary plane
by Oliver Clarke 01 Apr '15

01 Apr '15

Hi, Currently the mclip command can be used to match a model specific clipping plane to the global near/far clipping planes, or the per model clipping plane of another model. I think it would be a handy addition to a future chimera version to be able to match the per model clipping plane to an arbitrary plane, pre-defined in axes/planes/centroids - e.g. “mclip #1 match p2” would match the per model clipping plane of model #1 to the plane named “p2”. This would be useful because it allows one to easily generate a per-model clipping plane that intersects a certain set of atoms, which also facilitates precise description of a figure. Cheers, Oliver.

2 2

Request for axes measurements
by Jayaraman T 31 Mar '15

31 Mar '15

Hi, I'm Jayaraman, pursuing PhD in Protein Engineering at Pusan National University. I just want to calculate the Axes for specific set of helical residue CA atoms and then to calculate the angle between them. Is it possible with Chimera and if so, please do help me out that how to do it. It would be grateful and looking forward to hear at your earliest convenience. Thank you. - - With warm regards Jram | 201493155 PhD Scholar jram(a)pusan.ac.kr Bio molecular Engineering Lab Dept of Chemical Engineering Pusan National University Busan, South Korea +8210 9906 3680

2 1

Picking not working in RDS
by Hurt, Darrell (NIH/NIAID) [E] 31 Mar '15

31 Mar '15

Hi Chimera friends, I'm using the Microsoft Remote Desktop Mac client to connect to a Microsoft Remote Desktop Server (an "RDS" running Windows Server 2012 R2). That RDS is serving out a Windows 8 environment and I'm running the latest Chimera release in it very successfully. Except that selection by picking (Ctrl-mouse_button_1) and centering (Ctrl-mouse_button_3) does NOT work. When I try to select by picking, the cursor changes from the arrow to hand momentarily. I also get this behavior when connecting to the RDS via any other client. Any thoughts on what might be happening and how to fix it? Thanks, Darrell -- Darrell Hurt, Ph.D. Acting Chief Bioinformatics and Computational Biosciences Branch (BCBB) OCICB/OSMO/OD/NIAID/NIH 5601 Fishers Lane, 4A31 North Bethesda, MD 20852 Office: 240-669-2741 Mobile: 301-758-3559 Web: BCBB Home Page<http://www.niaid.nih.gov/about/organization/odoffices/omo/ocicb/Pages/bcbb.…> Twitter: @niaidbioit<https://twitter.com/niaidbioit> , @NIH3Dprint<https://twitter.com/nih3dprint> Disclaimer: The information in this e-mail and any of its attachments is confidential and may contain sensitive information. It should not be used by anyone who is not the original intended recipient. If you have received this e-mail in error please inform the sender and delete it from your mailbox or any other storage devices. National Institute of Allergy and Infectious Diseases shall not accept liability for any statements made that are sender's own and not expressly made on behalf of the NIAID by one of its representatives.

2 4

Toggle active vs displayed for checkboxes at bottom of screen?
by Oliver Clarke 30 Mar '15

30 Mar '15

Hi, this is more of a feature suggestion than anything, and I guess it is unlikely to be implemented before Chimera 2, but thought I’d put it out there anyway. Currently the checkboxes at the bottom of the main display control whether or not a model is active - this is very handy when reorienting models with respect to one another. I would love a similar means to toggle the display of individual models without going to the model panel. One way would be to alter the “Active models” text label (circled in screenshot) to “Displayed models”, and change the behaviour/meaning of the checkboxes accordingly, similar to the A/I/R/S selection mode changer at the bottom right of the main display. So you would click “Active models”, and the text would switch to “Displayed models”, and the configuration of the checkboxes would change so that the displayed models (rather than the active models) are checked. Is this something that would be possible in a future version? Best, Oliver.

3 3

Re: [Chimera-users] Running graphical Chimera from remote server
by Hurt, Darrell (NIH/NIAID) [E] 30 Mar '15

30 Mar '15

Hello again, I still have to try VirtualGL (as suggested on a previous thread). However, I want to share the following information: I learned that the 60 fps I was getting on glxgears is likely due to Ubuntu having "sync" enabled by default, which means that glxgears always runs at about 60 fps. So let's ignore that bit of information. http://askubuntu.com/questions/31488/how-to-test-if-my-video-card-has-3d-su pport My glxinfo (see below at "OpenGL renderer string") suggests that I have an NVIDIA graphics card installed. But "lspci -vnn | grep -i VGA" suggests that I have a simple Matrox graphics card: <output of lspci> 01:00.1 VGA compatible controller [0300]: Matrox Electronics Systems Ltd. MGA G200EH [102b:0533] (rev 01) (prog-if 00 [VGA controller]) <output of glxinfo> name of display: localhost:11.0 display: localhost:11 screen: 0 direct rendering: No (If you want to find out why, try setting LIBGL_DEBUG=verbose) server glx vendor string: SGI server glx version string: 1.4 server glx extensions: GLX_ARB_create_context, GLX_ARB_create_context_profile, GLX_ARB_multisample, GLX_EXT_import_context, GLX_EXT_visual_info, GLX_EXT_visual_rating, GLX_OML_swap_method, GLX_SGIS_multisample, GLX_SGIX_fbconfig client glx vendor string: Mesa Project and SGI client glx version string: 1.4 client glx extensions: GLX_ARB_create_context, GLX_ARB_create_context_profile, GLX_ARB_get_proc_address, GLX_ARB_multisample, GLX_EXT_import_context, GLX_EXT_visual_info, GLX_EXT_visual_rating, GLX_EXT_framebuffer_sRGB, GLX_EXT_create_context_es2_profile, GLX_MESA_copy_sub_buffer, GLX_MESA_multithread_makecurrent, GLX_MESA_swap_control, GLX_OML_swap_method, GLX_OML_sync_control, GLX_SGI_make_current_read, GLX_SGI_swap_control, GLX_SGI_video_sync, GLX_SGIS_multisample, GLX_SGIX_fbconfig, GLX_SGIX_pbuffer, GLX_SGIX_visual_select_group, GLX_EXT_texture_from_pixmap, GLX_INTEL_swap_event GLX version: 1.4 GLX extensions: GLX_ARB_create_context, GLX_ARB_create_context_profile, GLX_ARB_get_proc_address, GLX_ARB_multisample, GLX_EXT_import_context, GLX_EXT_visual_info, GLX_EXT_visual_rating, GLX_MESA_multithread_makecurrent, GLX_OML_swap_method, GLX_SGI_make_current_read, GLX_SGIS_multisample, GLX_SGIX_fbconfig, GLX_SGIX_pbuffer OpenGL vendor string: NVIDIA Corporation OpenGL renderer string: NVIDIA GeForce 9400M OpenGL Engine OpenGL version string: 1.4 (2.1 NVIDIA-8.16.81 310.40.00.20f04) OpenGL extensions: GL_ARB_depth_texture, GL_ARB_draw_buffers, GL_ARB_fragment_program, GL_ARB_fragment_program_shadow, GL_ARB_imaging, GL_ARB_multisample, GL_ARB_multitexture, GL_ARB_occlusion_query, GL_ARB_point_parameters, GL_ARB_point_sprite, GL_ARB_shadow, GL_ARB_texture_border_clamp, GL_ARB_texture_compression, GL_ARB_texture_cube_map, GL_ARB_texture_env_add, GL_ARB_texture_env_combine, GL_ARB_texture_env_crossbar, GL_ARB_texture_env_dot3, GL_ARB_texture_mirrored_repeat, GL_ARB_texture_non_power_of_two, GL_ARB_texture_rectangle, GL_ARB_transpose_matrix, GL_ARB_vertex_program, GL_ARB_window_pos, GL_EXT_abgr, GL_EXT_bgra, GL_EXT_blend_color, GL_EXT_blend_equation_separate, GL_EXT_blend_func_separate, GL_EXT_blend_minmax, GL_EXT_blend_subtract, GL_EXT_clip_volume_hint, GL_EXT_draw_range_elements, GL_EXT_fog_coord, GL_EXT_framebuffer_object, GL_EXT_multi_draw_arrays, GL_EXT_point_parameters, GL_EXT_rescale_normal, GL_EXT_secondary_color, GL_EXT_separate_specular_color, GL_EXT_shadow_funcs, GL_EXT_stencil_two_side, GL_EXT_stencil_wrap, GL_EXT_texture_compression_dxt1, GL_EXT_texture_compression_s3tc, GL_EXT_texture_edge_clamp, GL_EXT_texture_env_add, GL_EXT_texture_filter_anisotropic, GL_EXT_texture_lod_bias, GL_EXT_texture_mirror_clamp, GL_EXT_texture_rectangle, GL_APPLE_packed_pixels, GL_ATI_draw_buffers, GL_ATI_texture_env_combine3, GL_ATI_texture_mirror_once, GL_ATIX_texture_env_combine3, GL_IBM_texture_mirrored_repeat, GL_INGR_blend_func_separate, GL_NV_blend_square, GL_NV_depth_clamp, GL_NV_fog_distance, GL_NV_fragment_program_option, GL_NV_fragment_program2, GL_NV_light_max_exponent, GL_NV_multisample_filter_hint, GL_NV_point_sprite, GL_NV_texgen_reflection, GL_NV_texture_rectangle, GL_NV_vertex_program2_option, GL_NV_vertex_program3, GL_SGIS_generate_mipmap, GL_SGIS_texture_border_clamp, GL_SGIS_texture_edge_clamp, GL_SGIS_texture_lod, GL_SUN_multi_draw_arrays -- Darrell Hurt, Ph.D. Acting Chief Bioinformatics and Computational Biosciences Branch (BCBB) OCICB/OSMO/OD/NIAID/NIH 5601 Fishers Lane, 4A31 North Bethesda, MD 20852 Office: 240-669-2741 Mobile: 301-758-3559Web: BCBB Home Page <http://www.niaid.nih.gov/about/organization/odoffices/omo/ocicb/Pages/bcbb .aspx#niaid_inlineNav_Anchor> Twitter: @niaidbioit <https://twitter.com/niaidbioit> , @NIH3Dprint <https://twitter.com/nih3dprint> Disclaimer: The information in this e-mail and any of its attachments is confidential and may contain sensitive information. It should not be used by anyone who is not the original intended recipient. If you have received this e-mail in error please inform the sender and delete it from your mailbox or any other storage devices. National Institute of Allergy and Infectious Diseases shall not accept liability for any statements made that are sender's own and not expressly made on behalf of the NIAID by one of its representatives. On 3/29/15 7:10 PM, "Hurt, Darrell (NIH/NIAID) [E]" <darrellh(a)niaid.nih.gov> wrote: >Hi Chimera friends, > >I want to run a graphics-enabled Chimera from a remote Ubuntu server. The >Ubuntu server is quite beefy (and not virtualized), but it has a stock >video card. When I run "glxgears", I can see the gears and get about 60 >fps (reported), but the gears don't move and aren't interactive. > >I found this message from the listserv a few years ago that recommends >the use of VirtualGL: >http://plato.cgl.ucsf.edu/pipermail/chimera-users/2012-September/007953.ht >ml > >My X windowing system is Xquartz on a Mac. Is my problem related to >software GL or to the poor graphics card? > >Related question: what X windowing system do you recommend on a Windows >system? CYGWIN? > >Thanks, >Darrell > >-- >Darrell Hurt, Ph.D. >Acting Chief >Bioinformatics and Computational Biosciences Branch (BCBB) >OCICB/OSMO/OD/NIAID/NIH > >5601 Fishers Lane, 4A31 >North Bethesda, MD 20852 >Office: 240-669-2741 >Mobile: 301-758-3559 >Web: BCBB Home >Page<http://www.niaid.nih.gov/about/organization/odoffices/omo/ocicb/Pages >/bcbb.aspx#niaid_inlineNav_Anchor> >Twitter: @niaidbioit<https://twitter.com/niaidbioit> , >@NIH3Dprint<https://twitter.com/nih3dprint> > >Disclaimer: The information in this e-mail and any of its attachments is >confidential and may contain sensitive information. It should not be used >by anyone who is not the original intended recipient. If you have >received this e-mail in error please inform the sender and delete it from >your mailbox or any other storage devices. National Institute of Allergy >and Infectious Diseases shall not accept liability for any statements >made that are sender's own and not expressly made on behalf of the NIAID >by one of its representatives. > >_______________________________________________ >Chimera-users mailing list >Chimera-users(a)cgl.ucsf.edu >http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users

2 2

Rename models from command line
by Oliver Clarke 30 Mar '15

30 Mar '15

Hi all, Is there any way to rename models from the command line? I can do so using the model panel, but I have a script that generates a bunch of different models, and I would like to give them descriptive names apart from the model id in the script. Oliver.

3 4

MatchMaker-vs-MatchAlign
by Egor Tchesnokov 30 Mar '15

30 Mar '15

Hello, 1. I used MatchMaker to align two structures, worked well and I got an alignment file-1. Then I clicked on MatchAlign and got alignment file-2. The alignment-1 is different from alignment-2. Could you please let me know what is the differnce between MatchMaker and MatchAlign? 2. In the same function I can't find RR distance map. Is it applicable to MatchMaker/MatchAlign? Thank you Egor

2 2

on the section command
by Smith Liu 30 Mar '15

30 Mar '15

Dear All, Wil you please tell me how can I stop or pause the movie playing as a specific point which I want to watch carefully when the section command is in the execution? Smith

2 2

Running graphical Chimera from remote server
by Hurt, Darrell (NIH/NIAID) [E] 29 Mar '15

29 Mar '15

Hi Chimera friends, I want to run a graphics-enabled Chimera from a remote Ubuntu server. The Ubuntu server is quite beefy (and not virtualized), but it has a stock video card. When I run "glxgears", I can see the gears and get about 60 fps (reported), but the gears don't move and aren't interactive. I found this message from the listserv a few years ago that recommends the use of VirtualGL: http://plato.cgl.ucsf.edu/pipermail/chimera-users/2012-September/007953.html My X windowing system is Xquartz on a Mac. Is my problem related to software GL or to the poor graphics card? Related question: what X windowing system do you recommend on a Windows system? CYGWIN? Thanks, Darrell -- Darrell Hurt, Ph.D. Acting Chief Bioinformatics and Computational Biosciences Branch (BCBB) OCICB/OSMO/OD/NIAID/NIH 5601 Fishers Lane, 4A31 North Bethesda, MD 20852 Office: 240-669-2741 Mobile: 301-758-3559 Web: BCBB Home Page<http://www.niaid.nih.gov/about/organization/odoffices/omo/ocicb/Pages/bcbb.…> Twitter: @niaidbioit<https://twitter.com/niaidbioit> , @NIH3Dprint<https://twitter.com/nih3dprint> Disclaimer: The information in this e-mail and any of its attachments is confidential and may contain sensitive information. It should not be used by anyone who is not the original intended recipient. If you have received this e-mail in error please inform the sender and delete it from your mailbox or any other storage devices. National Institute of Allergy and Infectious Diseases shall not accept liability for any statements made that are sender's own and not expressly made on behalf of the NIAID by one of its representatives.

1 0

Tutorials
by Bobby Barnett 28 Mar '15

28 Mar '15

I was wondering if anyone had interactive tutorial session(s) that can be used to educate a small class of undergraduate students. Bobby Barnett Department of Chemistry University of Chemistry

2 5

Problems with reproducing 'write relative' output
by Dan McNamara 25 Mar '15

25 Mar '15

Hi all, I am encountering a strange inconsistency with the 'relative' option for the "write" command. Searching for similar problems on the mailing list came up short, so I thought I'd ask about it. To put it briefly, I run chimera with --nogui and perform ~780,000 total matches in which I write out new coordinates using 'relative.' My issue lies in automating this process with .com scripts, which creates different results from typing out the exact same commands into the command line of the GUI. This has been plaguing me for a year or so, and most recently I can cite using version 1.10.1 (build 40427). Below is a small example to illustrate my operations. The "B-new_#.pdb" output files will contain different coordinates if they originate from automated runs or manual execution. This problem is sporadic and I cannot narrow down the source of the discrepancy. From best I can tell, there are occasionally output files from automation (e.g. B-new_#.pdb) that are in the correct coordinate reference frame. However, many (most?) are not and it is a big obstacle. I hope that someone may have some insight on this problem. /////////// example /////////// open A.pdb open B.pdb open C.pdb ##write out selection from #0 [this file is always consistent] select #0:... write selected relative 0 #0 A-new.pdb ##move #2 onto #0 match #2:... #0:... ##move #1 onto #2 match #1:... #2:... ##write out selection from #1 and save it in its new position relative to #0 [inconsistent] select #1:... write selected relative 0 #1 B-new_1.pdb ~select all ##move #1 slightly match #1:... #2:...+1 ##write out selection from #1 and save it in its new position relative to #0 [inconsistent] select #1:... write selected relative 0 #1 B-new_2.pdb ~select all ##move #1 slightly more match #1:... #2:...+2 ##write out selection from #1 and save it in its new position relative to #0 [inconsistent] select #1:... write selected relative 0 #1 B-new_3.pdb /////////// example /////////// Best, Dan McNamara

2 1

black bar in saved image
by sunyeping 25 Mar '15

25 Mar '15

Dear all,I find in the saved volume image (png file), there are some black bar whick mask the image? Please see the following link:https://www.dropbox.com/s/nvy3douvafymqxd/image.png?dl=0Do you know how to get rid of them?Thank you very much! Yeping Sun Institute of Microbiology, Chinese Academy of Sciences

2 1

Creating MultiScale surface from a script
by Cruz, Phil (NIH/NIAID) [C] 24 Mar '15

24 Mar '15

Hello Chimera users, I need to create a MultiScale surface from a Chimera python script when the model is of large structures such as virus particles. On a Mac, the following lines work great: import MultiScale d = MultiScale.show_multiscale_model_dialog() d.make_multimers_cb(molecules=None) and produce the desired surface output. Unfortunately, where I really need this to work is on a Linux server using the special version of Chimera that doesn't require a graphics head. When I run the script there I get errors, presumably when the script tries to create the dialog: Traceback (most recent call last): File "/opt/UCSF/Chimera64-1.10.1/share/chimeraInit.py", line 698, in init midas_text.doRunScript("runscript", script) File "/opt/UCSF/Chimera64-1.10.1/share/Midas/midas_text.py", line 2206, in doRunScript execfile(scriptPath, scriptGlobals) File "Chimera_Molecular.py", line 301, in <module> d = MultiScale.show_multiscale_model_dialog() File "/opt/UCSF/Chimera64-1.10.1/share/MultiScale/__init__.py", line 2296, in show_multiscale_model_dialog return dialogs.display(MultiScale_Model_Dialog.name) File "/opt/UCSF/Chimera64-1.10.1/share/chimera/dialogs.py", line 76, in display dialog = find(name, create=1) File "/opt/UCSF/Chimera64-1.10.1/share/chimera/dialogs.py", line 61, in find return d() File "/opt/UCSF/Chimera64-1.10.1/share/chimera/baseDialog.py", line 665, in __init__ from tkgui import windowSystem File "/opt/UCSF/Chimera64-1.10.1/share/chimera/tkgui.py", line 73, in <module> import Togl ImportError: No module named Togl Is there a way to have the script create the multiscale surface without invoking the dialog (if that's what the issue is), or some other approach to get this to work in the required environment? Thanks for any suggestions. -- Phillip Cruz, Ph.D. Contractor, Medical Science & Computing Computational Structural Biologist Computational Biology Section Bioinformatics and Computational Biosciences Branch (BCBB) OCICB/OSMO/OD/NIAID/NIH 31 Center Dr., Room 3B62 Bethesda, MD 20892-0485 Office: 301-451-1089 http://bioinformatics.niaid.nih.gov<http://bioinformatics.niaid.nih.gov/> (Within NIH) http://exon.niaid.nih.gov<http://exon.niaid.nih.gov/> (Public) Disclaimer: The information in this e-mail and any of its attachments is confidential and may contain sensitive information. It should not be used by anyone who is not the original intended recipient. If you have received this e-mail in error please inform the sender and delete it from your mailbox or any other storage devices. National Institute of Allergy and Infectious Diseases shall not accept liability for any statements made that are sender's own and not expressly made on behalf of the NIAID by one of its representatives.

2 1

command related to the alignment of the ribbon form and the ball & stick form
by Smith Liu 24 Mar '15

24 Mar '15

Dear All, If by Chimera we open a PDB and a corresponding mrc file, and we display the molecule by both ribbon format and ball & stick format, we can find there is a minor disalignment between the ribbon representation and the ball & stick representation. I remember I have read somewhere by a ribbon related command the Chimera can align the ribbon representation and ball & stick presentation nicely, but I have forgotten where I read the tutorial or introduction. Will you please tell me the command to have chimera to get good alignment of the ribbon representation and ball & stick presentation, or tell me the weblink of the corresponding tutorial or reference? I am looking forward to getting your reply. Smith

2 1

Sampling the density of a map in the vicinity of each atom in a model
by Mario J. Borgnia 23 Mar '15

23 Mar '15

Hi, I have a density map in mrc format and a model in pdb format. I would like to calculate the average density within a 3D window centered at the coordinates of each atom in the model and include the result as an attribute for the atom. Despite some effort searching for the answer, I cannot find out how to do this easily. I could try and go deep into Python programming but, before I reinvent the wheel, I'd like to know whether there already is a simpler way to do this in Chimera. Thanks! Mario -- +---------------------------------+ | Mario J. Borgnia, Ph. D. | | | | Lab of Cell Biology | | NIH NCI/CCR | | 50 South Drive Rm 4306 MSC 8008 | | Bethesda MD 20892-8008 | | | | Tel: (301) 594 0563 | | mborgnia(a)nih.gov | +---------------------------------+

2 2

Re: [Chimera-users] Query on morph command in Chimera
by Elaine Meng 22 Mar '15

22 Mar '15

On Mar 22, 2015, at 9:38 AM, vinothkumar <vkumar(a)mrc-lmb.cam.ac.uk> wrote: > Dear Elaine, > I have three different volumes/maps of the same protein but differing in conformation. I am able to use the command line - ‘vop morph #0,1,2 constantVolume True model #3 frames 200’ which then interpolates and I could produce a movie. However, by default the threshold for the volumes and level goes to 4. Is there a way to insert in the command line to tell that all the volumes be set to level (1) and a certain threshold. > Thank you > Regards > Vinoth Dear Vinoth, You should adjust the threshold of the maps before the morphing. This could be done with the “volume” command if you are not using GUIs, e.g. command: volume all level 1 <http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/volume.html#general> Actually, my own experiments suggest it is only the threshold of the first input map that matters in setting the starting level of the morph. If you use the “constantVolume” morphing option, however, the threshold will be adjusted to keep the enclosed volume constant. If you really wanted to keep it exactly at 1 through the whole morph, you would not use that option. I hope this helps! For future questions, we recommend you send them to chimera-users(a)cgl.ucsf.edu (CC’d here) instead of to me or other individuals directly, unless you are including private data. Thanks, Elaine ----- Elaine C. Meng, Ph.D. UCSF Computer Graphics Lab (Chimera team) and Babbitt Lab Department of Pharmaceutical Chemistry University of California, San Francisco

4 4

problem related to open the chimera file by Coot
by Smith Liu 21 Mar '15

21 Mar '15

Dear All, I have fit a PDB file with a MRC by the Chimera. I saved both the fitted PDB file and the fitted MRC file by coot. When I open the fitted PDB file and the fitted MRC file by Chimera at the same time, they still fit perfectly. However when I open the Chimera fitted PDB file and the Chimera fitted MRC file by coot, the protein backbone (and the sidechains) does not fit the MRC file anymore. Will you please tell me what should I do in order to have the Chimera fitted PDB file and the Chimera fitted MRC file still fit in Coot? I am looking forward to getting your reply. Smith

2 2

Constraining map-to-map fitting based on symmetry?
by Oliver Clarke 20 Mar '15

20 Mar '15

Hi, If I have two maps of the same symmetry (say C4), and I align them using fit in map, is there any way to force chimera to keep the rotation axes of the two maps aligned, and to only search rotation around (and translation along) that axis, rather than arbitrary rotations? Cheers, Oliver.

2 3

Contour levels
by Nicholas Michelarakis 20 Mar '15

20 Mar '15

Hello all, Forgive me if this is a silly question but I was under the impression that when you opened two density maps in volume viewer and used the Fit in Map tool, with the Use only data above contour level option ticked, the correlation would change according to the contour level you selected in the volume viewer. In my case this does not happen. If this is indeed a bug, I would be happy to provide further information. Atm, I'm using Chimera 1.10. Thank you very much in advance. Best, Nick

2 1

Total number of atoms in all submodels
by Cruz, Phil (NIH/NIAID) [C] 19 Mar '15

19 Mar '15

Dear Chimera users, When I fetch a biological unit into Chimera from PDB, the symmetry related units are placed into separate submodels. I would like to find the total number of atoms in all submodels that make up the biological unit within a python script. Constructs (and variations that I've tried) such as the following only give the number of atoms in each sub model: numAtoms = len(openModels.list(modelTypes=[Molecule])[0].atoms) Is there some way to get the grand total of atoms? Alternatively, if the script could find the total number (or list) of the submodels, I could probably still handle it. My example is the virus capsid 1ej6. Thanks for any suggestions. -- Phillip Cruz, Ph.D. Contractor, Medical Science & Computing Computational Structural Biologist Computational Biology Section Bioinformatics and Computational Biosciences Branch (BCBB) OCICB/OSMO/OD/NIAID/NIH 31 Center Dr., Room 3B62 Bethesda, MD 20892-0485 Office: 301-451-1089 http://bioinformatics.niaid.nih.gov<http://bioinformatics.niaid.nih.gov/> (Within NIH) http://exon.niaid.nih.gov<http://exon.niaid.nih.gov/> (Public) Disclaimer: The information in this e-mail and any of its attachments is confidential and may contain sensitive information. It should not be used by anyone who is not the original intended recipient. If you have received this e-mail in error please inform the sender and delete it from your mailbox or any other storage devices. National Institute of Allergy and Infectious Diseases shall not accept liability for any statements made that are sender's own and not expressly made on behalf of the NIAID by one of its representatives.

2 2

Interaction interface between two residues and measure buriedArea
by ms 19 Mar '15

19 Mar '15

Dear all, I am trying to understand what is the area buried between two residues (or between one residue and two other ones) in several structures. It seemed to me the best way to do this was to use the command "measure buriedArea :res1 :res2" on the full structure (once all non-protein components have been cleaned out of it). The overall trend of these measures seems consistent with visual inspection. However I have a couple of questions about it: 1) I am unsure about what value is best to use/compare, buriedSAS and buriedSES. It seems to me that the difference between the two is essentially in the geometrical protocol used. Is there any general wisdom about what do they physically mean and how should I approach each value? 2) In some cases I receive *negative* values for buriedSAS and/or buriedSES. Looking on the ML it seems that it is a symptom of a buggy or failed area calculation (e.g. http://www.cgl.ucsf.edu/pipermail/chimera-users/2010-April/005119.html ). However it also says that I should see a "fallback" warning on the Reply Log -I see nothing of the sort. How do I troubleshoot this? Also in cases where the calculated value is positive, how can I be sure that the values are indeed meaningful? Any other feedback on the kind of measurement I am trying to do is more than welcome. Thanks a lot, Massimo

5 14

how to convert mtz file to a file which can be openne by Chimera
by Smith Liu 19 Mar '15

19 Mar '15

Dear All, Will you please tell me how to convert mtz file to a file which can be openne by Chimera? Smith

2 1

Chimera service unavailable
by SA 18 Mar '15

18 Mar '15

Dear developer, I am using Chimera PDB2PQR module. However, it is failed. The error message is below: self.backend = Backend(service, url) File "/soft/chimera/Chimera64-1.10.1/share/WebServices/opal_client.py", line 21, in __init__ raise NonChimeraError("Web service appears " NonChimeraError: Web service appears to be down. See Reply Log for more details. Service 'opal:pdb2pqr_1.8' is unavailable. See Reply Log for more details. For your indormation, I am using the latest linux version of Chimera. Thank you. Siti Postdocs UCSD

2 1

How to get prmtop file corresponding to a pdb file
by Mahendra B Thapa 18 Mar '15

18 Mar '15

Dear Chimera users, To get *Amber prmtop file* corresponding to a *pdb file using Chimera*, I did the following steps: 1] tools -> Amber -> write prmtop 2] In a pop-up window 'save file in chimera', I gave file name, say, aa.prmtop, then clicked save. 3] No charges message appeared, then I clicked 'Assign charges'. 4] In a pop-up window 'Add charge', I selected 'AM1-BCC' in other residues, then clicked ok. 5] In a pop up window 'specify net charge', I selected 'AMI-BCC' in charge method 6] I got 'chimera error' - failure running ANTECHEMBER for residue HIS. Check reply log for details. 7] In reply log, '*Failure running ANTECHEMBER for residue HIS*' seen. Any help in this issue will be a good help for me. Thank you. Mahendra Thapa University of Cincinnati,OH

3 2

visualization of electrostatic potential
by Christos Deligkaris 18 Mar '15

18 Mar '15

Dear Chimera Users, I am interested in visualizing the electrostatic potential calculated by the APBS tool Chimera provides and I have several questions regarding this. 1) We have been having troubles running the PDB2PQR/APBS tools (I am not sure if we had trouble with both, it has been a few weeks) that Chimera provides, and I believe we fixed this by changing the version number at the very end of the web server address. Is there anything else we could have done to fix this? 2) I noticed that the numbers that are displayed when I visualize the APBS-calculated potential on the MSMS are very different from the numbers shown when I use the volume viewer (I mapped the values to atoms and then displayed those as labels). Why are those numbers not the same? Is it because the volume data will map the potential near the atoms and display that, whereas the potential on the surface (or actually 1.4 Angstroms from the surface) will display the potential value on the surface? Should I assume that the units are the same? 3) I am interested in finding the "patches" of the molecular surface the have fairly positive potential. I am familiar with the PatchFinderPlus server but I am interested in finding several patches, not just the largest one. I have not found any tool related to this in Chimera, please let me know if I am wrong. So,I found the redarea.py script avalaible from Chimera's website and I am thinking about modifying it in order to do this. The script iterates over all surface pieces. Is it possible to calculate the distance between vertices and atoms? Is it possible to obtain the potential value at a vortex? I think this would help me get started... Thank you very much, Christos Deligkaris, PhD Assistant Professor of Physics, Drury University 900 N Benton Ave, Springfield MO, 65802 Office Phone: (417) 873-7234 www2.drury.edu/christos <http://www2.drury.edu/christos/index.html> @DeligkarisGroup <https://twitter.com/DeligkarisGroup> +ChristosDeligkaris <http://google.com/+ChristosDeligkaris>

3 2

Re: [Chimera-users] FW: How to get prmtop file corresponding to a pdb file
by Mahendra B Thapa 18 Mar '15

18 Mar '15

Dear Dr. Victor, Thank you for reminding me about the protonation state of HIS. As a trial, I changed all HIS to HID and again repeated my previous steps. This time, I got some python track back errors. I suspect the source of error could be the initial pdb file which I used in chimera ; the pdb file was obtained from other software (like online server), so some atom names were different from that of AMBER. My goal was to get AMBER prmtop corresponding to pdb files which were generated by programs other than AMBER so that I could use ambertools. So, may I need to convert all such pdb files compatible with amber pdb format ? Thank you for help, Mahendra Thapa On Wed, Mar 18, 2015 at 5:08 AM, Thapa, Mahendra (thapamb) < thapamb(a)mail.uc.edu> wrote: > > ------------------------------ > *From:* Victor Muñoz > *Sent:* Wednesday, March 18, 2015 3:08:18 AM (UTC-06:00) Central America > *To:* Thapa, Mahendra (thapamb) > *Cc:* chimera-users(a)cgl.ucsf.edu > *Subject:* Re: [Chimera-users] How to get prmtop file corresponding to a > pdb file > > Dear Mahendra, > > The error is due to HIS not being the standard residue name for Histidine > in the Amber force field. In Amber, histidine residues are named either > HIE, HID or HIP depending on their protonation state ( > http://ambermd.org/Questions/HIS.html) Change the name of the histidines > on the PDB file accordingly to their protonation state and it should work > fine. > > Kind regards, > > Victor > > 2015-03-18 6:39 GMT+01:00 Mahendra B Thapa <thapamb(a)mail.uc.edu>: > >> Dear Chimera users, >> >> To get *Amber prmtop file* corresponding to a *pdb file using Chimera*, >> I did the following steps: >> 1] tools -> Amber -> write prmtop >> 2] In a pop-up window 'save file in chimera', I gave file name, say, >> aa.prmtop, then clicked save. >> 3] No charges message appeared, then I clicked 'Assign charges'. >> 4] In a pop-up window 'Add charge', I selected 'AM1-BCC' in other >> residues, then clicked ok. >> 5] In a pop up window 'specify net charge', I selected 'AMI-BCC' in >> charge method >> 6] I got 'chimera error' - failure running ANTECHEMBER for residue HIS. >> Check reply log for details. >> 7] In reply log, '*Failure running ANTECHEMBER for residue HIS*' seen. >> >> Any help in this issue will be a good help for me. >> >> Thank you. >> Mahendra Thapa >> University of Cincinnati,OH >> >> _______________________________________________ >> Chimera-users mailing list >> Chimera-users(a)cgl.ucsf.edu >> http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users >> >> >

2 1

on color zone
by Smith Liu 17 Mar '15

17 Mar '15

Dear All, I have a mrc file for a protein complex with a 3-fragment sheets (from A to B, from B to C, from C to D, and they were connected by H-bonds). Will you please introduce to me on how to use color zone to split to get the mrc file for each fragment of the sheets? I have tried to do it by myself, however I find the splitting was not done based on fragment, i.e., I cannot get the mrc file for each fragment of the sheets. I am looking forward to getting your reply. Smith

2 1

Chimera tutorial at VIZBI conference in Boston on Tuesday, March 24
by Samuel Hertig 16 Mar '15

16 Mar '15

Dear Chimera users, In case you are interested in learning the basics of Chimera and if you can make it to Boston on March 24, the RBVI will be teaching a small Chimera tutorial as part of the Visualizing Biological Data (VIZBI) conference. Details and registration: http://www.vizbi.org/2015/Program/#Tuesday Best regards, Sam Hertig, PhD Postdoctoral researcher Resource for Biocomputing, Visualization, and Informatics University of California, San Francisco

1 0

UCSF Chimera error: X Error of failed request: BadDrawable (invalid Pixmap or Window parameter)
by Muhammad Sufian 16 Mar '15

16 Mar '15

Dear all Chimera users, I tried to install following version of Chimera; chimera-1.10.1-linux_x86_64.bin but after installation, following error appeared after its execution; sufian@system01:/usr/local/chimera/bin> ./chimera X Error of failed request: BadDrawable (invalid Pixmap or Window parameter) Major opcode of failed request: 72 (X_PutImage) Resource id in failed request: 0x2200064 Serial number of failed request: 1357 Current serial number in output stream: 1363 I also tried to install Chimera 1.9, but got same error. My system specifications are; Processor: x86_64 Operating-system: GNU/Linux (OpenSUSE 12.1) Kernel: Linux 3.4.63-2.44-desktop Kindly guide me how to resolve this issue. Regards. Muhammad SufianM.Phil. scholar (Molecular Medicine),Computational Biology Laboratory,Lab no. P-103, PCMD-Extension,Dr. Panjwani Center for Molecular Medicine and Drug ResearchInternational Center for Chemical and Biological SciencesUniversity of Karachi, Karachi, Pakistan

2 5

a question related to segment map
by Smith Liu 16 Mar '15

16 Mar '15

Dear All, After I segment a map into diffrent color, will you please tell me in detail how can I erase some specific parts of the map which I do not nee for some purpose? Smith

2 1

Selecting just 1 carbohydrate chain
by Tein Leong (PG Research) 12 Mar '15

12 Mar '15

Hello, Thank you for the amazing programme. I’m using it for the first time and I can do most of the things I needed to do. I am writing to ask if it is possible to “color” just one carbohydrate chain ? I’m coloring an antibody crystal structure which is made up of two heavy chains. Each chain has a carbohydrate attached. I wanted to color both carbohydrate in different colors but I’m unable to do so. When I go to “Select” > “Residue” , the programme can only select all carbohydrate residues. Please advice. Thanks ! Best wishes, Tein Foong Leong PhD Student Jacqui Wood Cancer Centre, Level 7 Ninewells Hospital University of Dundee DD1 9SY The University of Dundee is a registered Scottish Charity, No: SC015096

2 1

binned color ranges for scolor?
by Oliver Clarke 11 Mar '15

11 Mar '15

Hi all, I was wondering whether it is possible to use binned ranges for scolor? Specifically, I would like to color one volume by the density values of a second volume, but rather than using a continuous color palette, I would like to highlight several discrete bins with particular colors - e.g coloring all values between 2 and 4 red, and all values between 6 and 8 green, while leaving other voxels uncolored (or colored gray, for example). Is this currently possible? I couldn’t figure it out - both cmap and map range seem to give continuous coloring. Best, Oliver.

2 1

How to lower the map size
by 李智海 09 Mar '15

09 Mar '15

Hi, Now I have a density map with a dimension of 820*820*820, it is too large for me to perform image processes with Chimera. So I was wondering whether there is a way to decrease the map size rather than re-reconstruct the map with bined particle image. please help me out! Thank you in advance! Zhihai Li

2 1

Re: [Chimera-users] Request for support - UCSF Chimera
by Eric Pettersen 04 Mar '15

04 Mar '15

On Mar 4, 2015, at 10:14 AM, Nishant Ray Velagapudi <nishray(a)uw.edu> wrote: > Hello Mr. Petterson, > > I was looking for some help with using UCSF chimera and I was wondering if I can directly email you my question. If this is not the appropriate address, my apologies, and could you point me in the appropriate direction for support? Hi Nishant, You should probably send questions like these to the chimera-users mailing list. That way other people can benefit from the answer, and it also allows others to answer in the cases where I'm not really the best person for the question. > My issue is that I want to automate the process of creating and saving structure based sequence alignments. It has been simple to open the file and align them, but I am not sure how to save the .fasta alignment file out. Here is what I have so far: > > for fn in file_names: > replyobj.status("Processing: " + fn) > rc("open " + fn) > rc(mm #0 #1 show true) > rc("close " + fn) > > The bold line is the problematic one - is there any way to write out the resulting alignment file in a programmatic fashion? Yes. You can locate the MultAlign Viewer (MAV) instance using the code in this previous chimera-users mail: http://plato.cgl.ucsf.edu/pipermail/chimera-users/2014-June/010017.html That message also describes how to save the alignment to a file, in conjunction with the info in this message: http://plato.cgl.ucsf.edu/pipermail/chimera-users/2014-June/010015.html You can search the chimera-users list from this page to see if your question has been asked previously: http://www.cgl.ucsf.edu/chimera/docs/feedback.html I hope this helps. --Eric Eric Pettersen UCSF Computer Graphics Lab http://www.cgl.ucsf.edu

1 0

concatenate
by sunyeping 04 Mar '15

04 Mar '15

Dear all,I am dealing with a EM map which is a double helical filament and composed of two protein chains winding each other. Each of the two chains is composed of multiple same protein subunits. Now I am trying to color one chain of the map with one color, and the other chain with another color. How should I do that?Besides, it is known that multiple such double helixes can constitute a longer helix by a head-to-tail way. I can open several helixes in chimera, but how to concatenate they correctly into a long helix?Could you help me with these questions?With many thanks! Yeping Sun Institute of Microbiology, Chinese Academy of Sciences

2 1

Command for zooming in
by Anindito Sen 04 Mar '15

04 Mar '15

Dear All, What is the command for zooming into a portion of a density map ? This will be just using the command line for the zoom-in instead of using the middle mouse button. Thanks Andy

2 1

Select residue by array index, not residue id
by Thomas Hrabe 03 Mar '15

03 Mar '15

Hi everyone, is there a way to select a residue by it’s array index, starting at 0 -> sequence length rather than using the residue id? Thank you in advance for your help, Thomas

2 1

Re: [Chimera-users] rmsd calculation for pre aligned structures
by Elaine Meng 02 Mar '15

02 Mar '15

Hi Oliver, In that case, you can just show the sequence of one protein, then use the sequence window menu: Structure… Associations to associate both structures with that one sequence, then sequence window Headers menu to show the RMSD of interest, etc. This would work even if the sequences of the structures were slightly different, as long as they were both similar enough to the sequence in Multalign Viewer that Chimera could figure out the association accurately. Elaine On Mar 2, 2015, at 3:12 PM, Oliver Clarke <olibclarke(a)gmail.com> wrote: > Hi Elaine, thanks for the reply! This sounds like exactly the right solution. However, when I try to add a sequence in the sequence dialog, Chimera hangs on “running Needleman-Wunsch”. > > I suspect this is because my sequences are so long (~5000 residues). > > But I don’t actually need to align them - they are already aligned as they are actually the same sequence, because I am looking at different conformations of the same protein. > > I guess I can create a dummy sequence alignment, with two identical sequences, but it would be nice to have some way of telling chimera not to align the sequences, or to allow for manual alignment in cases such as these - perhaps a “Do not align” or “align manually” checkbox in the add sequence dialog would do the trick. > > Thanks again, > Oliver.

1 0

rmsd calculation for pre aligned structures
by Oliver Clarke 02 Mar '15

02 Mar '15

Hi, I have two structures that are pre-aligned in another program, and I would like to calculate a per-residue rmsd between the two structures and assign this to an attribute that I can use to color one of the structures. Is there any way to do this at present? The mavRMSD attribute only seems to be calculated after matchmaker, not after using the regular “rmsd” command. Perhaps there is a way to do a dummy run of matchmaker, somehow forcing both structures to remain fixed? Best, Oliver.

3 2

explicit parentheses in atom selections
by Oliver Clarke 02 Mar '15

02 Mar '15

Hi all, Would it be possible at some point to allow the use of explicit parentheses when grouping selections on the command line? It would make it much easier/clearer when using multiple boolean operations (e.g "ligand z<5 &~ ligand &~ solvent” could be “ligand z<5 &~(ligand | solvent)”). Is there a reason why this kind of syntax is not used at the moment? Also, just a random thought for a future version - it would be really handy to have something like the “Macro recorder” function of ImageJ or similar programs. This records whatever the user does using the GUI (while the recorder is active) and spits out the equivalent scripting commands to a dialog, where they can be saved and reused, or turned into a script. This can be a very powerful tool for learning the command line functionality of a complex application in an interactive manner. Cheers, Oliver.

2 2

Embed chimera into my gui code
by hy liao 26 Feb '15

26 Feb '15

Hi, I am creating a GUI using tkinter, and I would like to embed chimera into it. Specifically, I would like to display a 3D molecule density map using chimera and inside my GUI. Chimera would be activated while the GUI is on. Is there documentation to accomplish this? Thanks! Stu

2 1

Inverse volume-erasing
by kzhang＠mrc-lmb.cam.ac.uk 25 Feb '15

25 Feb '15

Dear Chimera users, Dose anyone know how to do inverse volume-erasing? Normally we use the eraser to erase density we don't want to keep. However, what shall I do if I want to keep the density within the eraser while getting rid of other parts? Many thanks Kai

3 2

Re: [Chimera-users] Inverse volume-erasing
by kzhang＠mrc-lmb.cam.ac.uk 25 Feb '15

25 Feb '15

Thanks! This is a good solution. But I wonder if there's an interactive way just like the eraser, but doing the opposite thing. It might help on difficult cases when you want a real time feedback what you are erasing and keeping just in Chimera without external command. Kai > difference ma p > <kzhang(a)mrc-lmb.cam.ac.uk>于2015年2月25日星期三写道： >> Dear Chimera users, >> Dose anyone know how to do inverse volume-erasing? Normally we use the eraser to erase density we don't want to keep. However, what shall I do if >> I want to keep the density within the eraser while getting rid of other parts? >> Many thanks >> Kai >> _______________________________________________ >> Chimera-users mailing list >> Chimera-users(a)cgl.ucsf.edu <javascript:;> >> http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users

1 0

about model coordinates
by Tat Cheung Cheng 20 Feb '15

20 Feb '15

Hi, I recently got a ~4.4A EM map and built a improved model using MDFF. The map and model is well aligned when loaded in Chimera. But it is off when I hit "center" in the volume viewer, as the origin is shifted. As I try to improve the model in Coot, the map and model isn't aligned at all in Coot, seemingly by the same offset when "center" button was hit. And also the values in the original index in the volume viewer is different from the values in the header in the density map. So I wonder what determines the origin index? and is it possible to find the shift and transform the coordinates in the model so as to align the model with the map in coot? Thanks a lot. Tc

2 1

problem installation on linux 64 bit
by Chiara Parravicini 20 Feb '15

20 Feb '15

<chimera-users(a)cgl.ucsf.edu>Dear Chimera Users, After installing Chimera from the precompiled version chimera-1.10.1-linux_x86_64.bin (the installation completed without any error), I get the following error when I run the binary: libGL error: unable to load driver: nouveau_dri.so libGL error: driver pointer missing libGL error: failed to load driver: nouveau libGL error: unable to load driver: swrast_dri.so libGL error: failed to load driver: swrast X Error of failed request: GLXBadContextTag Major opcode of failed request: 156 (GLX) Minor opcode of failed request: 5 (X_GLXMakeCurrent) Serial number of failed request: 1455 Current serial number in output stream: 1455 I would appreciate any help or suggestion! Many Thanks in advance Chiara

3 2

Restricting "select by conservation" to a subset of sequences in a multiple alignment
by Oliver Clarke 18 Feb '15

18 Feb '15

Hi, I very much like the sequence tools within Chimera - I find the capacity to load in multiple alignments and select by conservation particularly handy. I have two suggestions that I would love to see in a future version of chimera (or maybe they are already in there and I don’t know about them!): 1. Selection of sub-groups of sequences within a multiple alignment for selecting residues by conservation. This would enable the easy selection/identification of residues that are specifically conserved within individual families with different functions or substrates - for example, a position that is an arginine, and always an arginine, in one subgroup of proteins that we know have a particular function, and always a glutamate in another group of proteins that we know possess a different but related function. It is possible to do this currently by loading in multiple different sequence alignments and saving and intersecting selections, but having a popup menu in the select/render by attribute tab where you could select or deselect sequences that contribute to the analysis of conservation would make it much easier. 2. Optional addition of sequence similarity, as well as identity, to the analysis of sequence conservation would make analysis of alignments including quite divergent sequences much more meaningful. Even better would be if the definition of similar groups for this purpose was user-adjustable and had allowance for overlapping groups of similar residues. Cheers, Oliver Clarke.

3 4

MD movie export from per frame script
by Ankit Agrawal 18 Feb '15

18 Feb '15

Hi all I want to write a script which can create movie (.mp4) from per frame script. Because I want to show all models in movie. Please help Thanks Ankit Agrawal Department of Chemical Science IISER Mohali, India

2 2

showing shapes in only selected Chimera animation scenes
by Michael Oldham 17 Feb '15

17 Feb '15

Chimera users, I am making an animation that needs to include either a box generated with an opened .bild file using a .box command or a rectangle generated by "shape rectangle” I wish the box or rectangle to appear in some scenes but not in others. My problem is, if I choose to show the box or rectangle in the first scene then it is shown in all subsequent scenes even if they were created in which the box or rectangle were not shown. Any suggestions on how to show the box or rectangle in only selected scenes of the animation. Thanks Mike Oldham Rockefeller University

2 1

Fwd: Problem running chimera through a command script
by Christos Deligkaris 17 Feb '15

17 Feb '15

Dear all, I have a python script that creates a .com command script for chimera and then attempts to run the chimera script. Here is the python script: ************** #!/usr/bin/env python import re,os,sys,math, decimal, subprocess from copy import deepcopy from openeye.oechem import * from decimal import * script_file=open("chimera_script_for_" + "ligand" +".com", 'w') script_file.write("open "+ "127D" +"\n") script_file.write("open " + "127D" +"\n") script_file.write("open "+ "ligand.pdb" +"\n") script_file.write("match #0:"+ "6" + " " + "#1:"+ "18" +"\n") script_file.write("write format pdb 2 "+ "flipped_" + "ligand" + "\n") script_file.write("stop") script_file.close print script_file.name chimera=subprocess.call(["chimera", " --nogui ", " chimera_script_for_ligand.com"]) print chimera #delete the script, leaving just the flipped ligand #os.remove(script_file.name) ******************* When I run the script, it successfully generates the .com chimera script. However, I get the following error message: ******************* [christos@biophysics Desktop]$ ./script.py chimera_script_for_ligand.com Traceback (most recent call last): File "/share/apps/PROGRAMS/CHIMERA/CHIMERA-1.9/share/__main__.py", line 69, in <module> value = chimeraInit.init(sys.argv) File "/share/apps/PROGRAMS/CHIMERA/CHIMERA-1.9/share/chimeraInit.py", line 638, in init splash.create() File "/share/apps/PROGRAMS/CHIMERA/CHIMERA-1.9/share/chimera/splash.py", line 25, in create sync=chimera.debug) File "/share/apps/PROGRAMS/CHIMERA/CHIMERA-1.9/lib/python2.7/lib-tk/Tkinter.py", line 1747, in __init__ self.tk = _tkinter.create(screenName, baseName, className, interactive, wantobjects, useTk, sync, use) _tkinter.TclError: no display name and no $DISPLAY environment variable 1 ********************** The python script successfully creates the .com script: ****************************** [christos@biophysics Desktop]$ more chimera_script_for_ligand.com open 127D open 127D open ligand.pdb match #0:6 #1:18 write format pdb 2 flipped_ligand stop *************************** and if I manually execute this chimera script, everything works out well: ************************** [christos@biophysics Desktop]$ chimera --nogui chimera_script_for_ligand.com Fetching 127D from web site www.rcsb.org 0 Kbytes received Fetch 127D: finished Done fetching 127D; verifying... Opening 127D... #0, chain A: dna (5'-D(*cp*gp*cp*gp*ap*ap*tp*tp*cp*gp*cp*G)- 3') #0, chain B: dna (5'-D(*cp*gp*cp*gp*ap*ap*tp*tp*cp*gp*cp*G)- 3') Opening VRML model in 127D - Nucleotides... 127D opened Opened VRML model in 127D - Nucleotides Opened 127D containing 1 model, 639 atoms, and 146 residues Fetching 127D from web site www.rcsb.org 88 Kbytes received 80 Kbytes received Fetch 127D: finished Done fetching 127D; verifying... Opening 127D... #1, chain A: dna (5'-D(*cp*gp*cp*gp*ap*ap*tp*tp*cp*gp*cp*G)- 3') #1, chain B: dna (5'-D(*cp*gp*cp*gp*ap*ap*tp*tp*cp*gp*cp*G)- 3') Opening VRML model in 127D - Nucleotides... 127D opened Opened VRML model in 127D - Nucleotides Opened 127D containing 1 model, 639 atoms, and 146 residues Opening ligand.pdb... ligand.pdb opened Opened ligand.pdb containing 1 model, 32 atoms, and 1 residues Executing match ['#0:6', '#1:18'], no iteration RMSD between 21 atom pairs is 0.630 angstroms Wrote flipped_ligand into /home/christos/Desktop ***************************** Why do I get the error message when I attempt to run chimera through script.py? Is this a python issue or a chimera issue? I appreciate any thoughts on solving this problem....thank you, Christos Deligkaris, PhD Assistant Professor of Physics, Drury University 900 N Benton Ave, Springfield MO, 65802 Office Phone: (417) 873-7234 www2.drury.edu/christos <http://www2.drury.edu/christos/index.html> @DeligkarisGroup <https://twitter.com/DeligkarisGroup> +ChristosDeligkaris <http://google.com/+ChristosDeligkaris>

2 5

Announcement: Release of DOCK 6.7
by Scott Brozell 12 Feb '15

12 Feb '15

We are pleased to announce the release of DOCK 6.7. DOCK is a suite of programs for molecular docking. The source code for DOCK 6.7 is available for download and free for academic users at http://dock.compbio.ucsf.edu/. In version 6.7 the default values for several input parameters were updated based mainly on a performance assessment using large data sets and employing multiple metrics including pose reproduction, cross-docking, and database enrichment For full information on what is new in DOCK 6.7, please visit: http://dock.compbio.ucsf.edu/DOCK_6/new_in_6.7.txt Sincerely, The DOCK Team Please visit us at the DOCK Web site. http://dock.compbio.ucsf.edu

1 0

UnitCell.molecule_center(): Attribute does not exist
by Alejandro Virrueta 11 Feb '15

11 Feb '15

Hello all, I'm trying to execute the script to fill out the crystal unit cell found here: http://www.cgl.ucsf.edu/chimera/docs/ProgrammersGuide/faq.html#q13 However, I get the error that UnitCell does not have the 'molecule_center' attribute (code line ~12). I checked with dir(), and yeah, it's not there: ['ModelessDialog', 'Unit_Cell_Dialog', '__builtins__', '__doc__', '__file__', '__name__', '__package__', '__path__', 'chimera', 'dialogs', 'find_model_by_name', 'outline_box', 'place_molecule_copies', 'remove_extra_copies', 'show_unit_cell_dialog'] Where did the attribute go? Is there some other module I could use? Could someone explain (or direct me to some webpage that explains) the process of expanding the ASU and I can write my own script. Any help is greatly appreciated! Thanks! Alex

2 1

Re: [Chimera-users] Problem using structureViz plugin for Cytoscape - Chimera
by Scooter Morris 11 Feb '15

11 Feb '15

Hi Juan, I've uploaded a new version of structureViz (1.0.1) that should fix your problem. Don't hesitate to let us know if there are other issues/enhancements you would like to see. -- scooter On Feb 9, 2015, at 6:39 AM, Juan Eiros Zamora <j.eiros-zamora14(a)imperial.ac.uk <mailto:j.eiros-zamora14@imperial.ac.uk>> wrote: > Hi all, > > I am using Chimera 1.10.1 and Cytoscape 3.2.0, and have installed the > structureViz2 plugin. > > I use structureviz to open chimera through cytoscape and load an MD > trajectory of a protein to create a RIN through Analysis --> Residue > interaction network... > > What I am trying to do is colour the nodes in cytoscape based on which > chain they belong to. My protein has three chains, and in chimera I > color each differently. In cytoscape, I go to apps --> structureviz--> > synchronize residue colors. Here I am given two options. The first one > is “Apply colors from current network view to associated Chimera > models”: This one works fine, but it is not what I want. What I want > to do is the second option: “Apply colors from associated Chimera > models to current network view”, but when I clik on OK it does nothing. > > Anybody having the same problem or some suggestions on how I might > solve this? > > Much appreciated, > > Juan Eiros >

1 0

Manipulate individual hexagons/pentagons in cage builder
by vamsee 10 Feb '15

10 Feb '15

Hello Elaine, I am trying to build a cage but with incorrectly attached pentagons and hexagons to represent mis-assembly. Though I don't believe it is possible using the cage builder, I was wondering if the wonderfully brilliant chimera team would know of any other ways of doing this or even with the cage builder. Thank you yet again, Vamsee

3 2

pymol get_view and set_view in Chimera?
by Jan Kosinski 10 Feb '15

10 Feb '15

Is there an equivalent of pymol get_view and set_view in Chimera? I am aware of savepos command, but I am looking for a feature like Pymol get_view command (http://www.pymolwiki.org/index.php/Get_View) which allows to copy the current view to a script as text and restore it using set_view (http://www.pymolwiki.org/index.php/Set_View) Thanks in advance, Jan

2 4

adding morphs to scene animation in Chimera
by Michael Oldham 10 Feb '15

10 Feb '15

Chimera users, I want to use a combination of various scene selections and morphs in a Chimera animation. Can someone tell me how morphs and scenes be combined to make one movie. Thanks Mike Oldham Rockefeller University

2 2

Problem using structureViz plugin for Cytoscape - Chimera
by Juan Eiros Zamora 09 Feb '15

09 Feb '15

Hi all, I am using Chimera 1.10.1 and Cytoscape 3.2.0, and have installed the structureViz2 plugin. I use structureviz to open chimera through cytoscape and load an MD trajectory of a protein to create a RIN through Analysis --> Residue interaction network... What I am trying to do is colour the nodes in cytoscape based on which chain they belong to. My protein has three chains, and in chimera I color each differently. In cytoscape, I go to apps --> structureviz--> synchronize residue colors. Here I am given two options. The first one is “Apply colors from current network view to associated Chimera models”: This one works fine, but it is not what I want. What I want to do is the second option: “Apply colors from associated Chimera models to current network view”, but when I clik on OK it does nothing. Anybody having the same problem or some suggestions on how I might solve this? Much appreciated, Juan Eiros

2 1

radial density map
by Jian Guan 09 Feb '15

09 Feb '15

Hi all, I have a mrc map of virus, T=7 icosahedron. Does anyone know how to make radial density map (averaged radial density VS radius of virion)? I tried chimera method of shortcut rd. It took long time but return one result. For the rd result, is it make along x, y or z axis? Or along assigned axis, like 5 fold symmetric axis? I have asymmetric extra density around 5-fold axis. What will happen in my case if use rd? Are there any other ways to make the plot besides of rd? Thank you so much. Sincerely yours, Jian

2 1

CHIMERA movie with parm7/dcd
by Francesco Pietra 09 Feb '15

09 Feb '15

Hello: In the recent past, issues in running movie with parm7 files were solved. Now, with ambertools14 problems come out again, unless I dropped into a silly mistake. That is, parm7/dcd was refused because no top file was found. Renaming the prm7 extension prmtop led to File "/opt/UCSF/Chimera64-1.7/share/Trajectory/formats/Amber/Amber.py", line 354, in SetupFrame raise IOError, "Trajectory/Restart header line is greater than 80 characters -- invalid!!" Thanks for your kind attention francesco pietra

2 1

command line
by Dieter Blaas 09 Feb '15

09 Feb '15

Hi, sorry for asking again an old question but I lost the solution: how do I show single solid planes of several volumes? I tried: 'vol all style solid plane one axis z' and many combinations thereof but it always tells me: 'volume invalid planes argument "one": Value must be one of x,y,z, got "one" when I try ' vol all style solid plane z' it tells me: 'volume invalid planes argument "z"' and so on.... Thanks a lot for your patience, best, 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 2

Re: [Chimera-users] Ramachandran plot statistics
by Nadia Raboana 08 Feb '15

08 Feb '15

Hi, It works perfectly, thank you for helping me. Best regards, Nadia Sent from Windows Mail From: Elaine Meng Sent: ‎February‎ ‎8‎, ‎2015 ‎12‎:‎52‎ ‎AM To: Nadia Raboana CC: chimera-users(a)cgl.ucsf.edu Subject: Re: [Chimera-users] Ramachandran plot statistics Hi Nadia, You would have to decide for yourself what probability cutoff you want to use as the boundary between allowed and disallowed. First assign the probabilities to the residues: in the Ramachandran Plot dialog, click the button at the bottom to “Assign Residue Probabilities” (or use the ramachandran command option “assign”): <http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/ramachandran/rama…> <http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/ramachandran.html> Then in Select by Attribute (Chimera menu: Select… By Attribute Value) choose attributes of “residues” with Attribute: ramaProb … that will show a histogram of the probability values. Then you can move the green bars to enclose the range of probabilities you want to select. After you use that dialog to select everything above your cutoff probability, then you can get a count of the selected residues by clicking the green magnifying glass at the bottom right corner of the Chimera window to show the Selection Inspector. <http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/render/render.htm…> <http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/inspection.html> Then you would divide by the total number of amino acid residues to get the percentage. You could count the total number of amino acid residues, for example, by using command “select protein” and then clicking the green magnifying glass. I hope this helps, Elaine ----- Elaine C. Meng, Ph.D. UCSF Computer Graphics Lab (Chimera team) and Babbitt Lab Department of Pharmaceutical Chemistry University of California, San Francisco On Feb 7, 2015, at 4:49 AM, Nadia Raboana <nadia_raboana(a)hotmail.fr> wrote: > Hello, > I am using UCSF Chimera to generate Ramachandran plot related to my protein structure. I would like to know how to get the statistics about the analysis. I would precisely know the percentage values of residues in most favored regions, allowed regions and disallowed regions. > Thank you for answering me. > Best regards, > Nadia Raboanatahiry

1 0

Change selection mode programatically
by Mario J. Borgnia 07 Feb '15

07 Feb '15

Hi, I am trying to find a way to split a complex selection pattern among a number of lines in a script. I was expecting that this would be an option to the select command, but apparently it is not. I know that I can combine "select" commands by changing the selection mode to "append" in the GUI, but I am struggling to find the way of doing this from command line or from a script. Is there a way to change the selection mode to "append" without using the GUI? Thanks Mario -- +---------------------------------+ | Mario J. Borgnia, Ph. D. | | | | Lab of Cell Biology | | NIH NCI/CCR | | 50 South Drive Rm 4306 MSC 8008 | | Bethesda MD 20892-8008 | | | | Tel: (301) 594 0563 | | mborgnia(a)nih.gov | +---------------------------------+

2 2

Ramachandran plot statistics
by Nadia Raboana 07 Feb '15

07 Feb '15

Hello,I am using UCSF Chimera to generate Ramachandran plot related to my protein structure. I would like to know how to get the statistics about the analysis. I would precisely know the percentage values of residues in most favored regions, allowed regions and disallowed regions.Thank you for answering me.Best regards, Nadia RaboanatahirySchool of Life Science and TechnologyHuazhong University of Science and TechnologyWuhan, CHINA

2 1

Help in writing python code
by Ankit Agrawal 06 Feb '15

06 Feb '15

Hi I am doing trajectory analysis which has 40000 frames. So for all frames I wanted to calculate the all three inertia axes. That I did. But problem is that when I run the python code it starts making all the three axes for each frames that becomes messy. So I want a little change in my code so that during trajectory analysis it should create axes but should delete or hide the previous axes also. So during run I will see only the axes which are related to that frame. Here is the python code in per-frame analysis. from chimera import runCommand frame = mdInfo['frame'] runCommand("measure inertia #1:0-120 color yellow") runCommand("measure inertia #1:120-237 color green") Thanks. Ankit

2 2

Help in executing python script in per-frame analysis
by Ankit Agrawal 05 Feb '15

05 Feb '15

Hi I am analysing a protein molecule (MD simulations). I wrote a python script to calculate inertia for the protein and ran it over 40000 frames. I got a data in log file but it is creating a ellipsoid around the molecule. I want 3 axes of the ellipsoid to be shown in that protein molecule not the ellipsoid as it covers the whole protein molecule. So that I can easily visualize the movement of all 3 axes of ellipsoid over the 40000 frames (means how the axes are changing during trajectory run). Thanks Ankit Department of Chemical Science IISER-Mohali, India

2 1

multiple Chimera sessions on Mac
by Vladimir Yarov-Yarovoy 03 Feb '15

03 Feb '15

Hello, Is it possible to open multiple Chimera sessions on Mac OS X Yosemite? I tried to start multiple chimera sessions from different terminal windows that was recommended in this thread - http://www.cgl.ucsf.edu/pipermail/chimera-users/2008-October/003231.html - but it doesn't work. Thank you, Vladimir

2 3

Command line problem
by Luca Pellegrini 03 Feb '15

03 Feb '15

Hi, I am having the following problem when I try to use the command line. When I select ‘Command line’ from the ‘Favorites’ menu, the command line window appears, but I can’t type anything in it. Or rather, I can type in it, because if I hit return the command gets executed, but I can’t see what I type. Any idea of what’s going on? I suspect that it might be something to do with Yosemite 10.10.2, since this problem arose after upgrading the MacOSX operating system. I am using Chimera 1.10.1, build 40415. Thanks for any advice. Kind regards, Luca Luca Pellegrini Department of Biochemistry University of Cambridge 80 Tennis Court Road Cambridge CB2 1GA UK Email: lp212(a)cam.ac.uk Tel: 0044-1223-760469 Fax: 0044-1223-766002 Sanger building, room 3.59

2 2

Error: UnboundLocalError: local variable 'i' referenced before assignment
by Letizia Meregalli 03 Feb '15

03 Feb '15

Hi, I have a tRNA molecule and I would like to measure the angles between two nitrogenous bases. I created two planes, each for each base. When I select both planes, in order to measure the angles between them (Structure Measurements Axes/Planes/Centroids) I get following message: <type 'exceptions.UnboundLocalError'> Exception in Tk callback Function: <bound method SortableTable._browseCmd of <CGLtk.Table.SortableTable instance at 0x117990638>> (type: <type 'instancemethod'>) Args: ('1',) Traceback (innermost last): File "/Applications/Chimera.app/Contents/Resources/lib/python2.7/site-packages/Pmw/Pmw_1_3/lib/PmwBase.py", line 1747, in __call__ return apply(self.func, args) File "/Applications/Chimera.app/Contents/Resources/share/CGLtk/Table.py", line 1043, in _browseCmd self.userBrowseCmd(self.selected()) File "/Applications/Chimera.app/Contents/Resources/share/StructMeasure/Geometry.py", line 268, in _tableCB self.feedback(sels) File "/Applications/Chimera.app/Contents/Resources/share/StructMeasure/Geometry.py", line 113, in feedback diff = " ".join(words2[i:]) <type 'exceptions.UnboundLocalError'>: local variable 'i' referenced before assignment UnboundLocalError: local variable 'i' referenced before assignment File "/Applications/Chimera.app/Contents/Resources/share/StructMeasure/Geometry.py", line 113, in feedback diff = " ".join(words2[i:]) See reply log for Python traceback. I really don’t know how can I solve this problem.. I’d be really grateful if you could help me or give me some hints! thank you a lot Letizia

2 1

Help Needed
by angshu＠klyuniv.ac.in 03 Feb '15

03 Feb '15

Hi, I need a small help. I need to find out the amino acid residues from one chain of a protein which are within 5A of a specific amino acid residue of some another chain of the same protein. For example: Say I have a protein 1XYZ.pdb containing A, B, C chains. I want to find which amino acids from chain B and C are within 5A distance from ALA 78 belonging to chain A. Please let me know. Thanks. Angshuman

2 1

Using Chimera to measure voxel densities
by Noble, Alex 30 Jan '15

30 Jan '15

Hi all, I have a pdb fit into an EM density in Chimera and I would like to be able to see what the nearest voxel densities are at specific atomic coordinates. Is there any way to do this in Chimera? Best, -Alex Noble

3 4

Question
by Monete Rajão Gomes 28 Jan '15

28 Jan '15

Hi, I have some problem these days about superimpose templates with matchmaker. I followed this tutorial http://plato.cgl.ucsf.edu/pipermail/chimera-users/2011-August/006686.html. There are regions on aminoacids sequences which are identical. So, when I superimpose them, they stay separated. I am sending a screenshot attached. It would be great if you could help me. Best, Monete

2 2

Morph conformations - estimated time to completion in status bar?
by Oliver Clarke 26 Jan '15

26 Jan '15

Hi all - silly little feature suggestion, but I was wondering if it would be possible to add an live-updating “x % complete” or “estimated time to completion” entry in the status bar when running Morph Conformations? When morphing two structures with very large numbers of atoms, morph conformations can take a long time (hours) to complete, and it would be handy to know whether it is nearly done or not (currently it just shows “computing interpolation 1” until it is finished). On a related note, is there any way to speed up morph conformations for large structures? E.g calculating transformations using C-alpha atoms only? Best, Oliver.

2 1

modelling a dimer
by David Norman (Staff) 23 Jan '15

23 Jan '15

I am trying to model a symmetrical dimer using modeller within Chimera. I cannot work out how to present the alignment and pdb. files in oder to do this. Any hints please. Thanks David David Norman Ph.D Nucleic Acids Research Group School of Life Sciences MSI/WTB Complex Dow St. Dundee, DD1 5EH, Scotland, UK. phone +44(0)1382 384798 mobile +44(0)7808572788 E.mail d.g.norman(a)dundee.ac.uk Web site: http://www.lifesci.dundee.ac.uk/people/david-norman The University of Dundee is a Scottish Registered Charity, No. SC015096 The University of Dundee is a registered Scottish Charity, No: SC015096

2 1

how to find non-bonded interactions between a protein and a ligand?
by Sam Paul D. 22 Jan '15

22 Jan '15

Hi, I need to find non-bonded interactions between a protein and a ligand through command line. Could anyone please suggest which chimera command to use with the syntax? With regards, Sam

3 4

Find subsequence programatically
by Thomas Hrabe 22 Jan '15

22 Jan '15

Hi Chimera Team, How can I automate the Sequence View -> Edit -> Find Subsequence Tool programmatically with the Chimera shell or better within python? Thank you, Thomas

2 2

Coot capabilities in Chimera
by Mario J. Borgnia 22 Jan '15

22 Jan '15

Dear all, WIth the advent of atomic resolution maps in Cryo-EM, we find ourselves refining structures. The blunt question is, should I move to Coot and leave behind Chimera for that or are there overlapping functions and what will I be missing if I choose to stay with Chimera. I have some ideas about the answers but would like to read what the developers and power users have to say on the issue. Thanks!!!

2 1

How to predict/model the structure of two known structure proteins
by Phan Trong Hoang 22 Jan '15

22 Jan '15

Dear UCSF developers, I am Hoang, working in IPK in Germany. I am following your series online tutorials about how to use UCSF Chimera. This program is very useful for me. However, I did not see your video showing how to predict/model the structure of a fusion protein that is made by combination of two known structure proteins in PDB (for example gi|534286646|pdb|4KRN| and gi|390136394|pdb|4DDF|A). Is this function available in your software? If yes can you show me how to manipulate with it. I would like appreciate it very much Thank you for your attention and I am looking forward to hearing from you. Best regards Hoang Phan

2 3

SDF file
by Hernando J Sosa 22 Jan '15

22 Jan '15

I have an sdf file with a small ligand molecule in 2D (it looks flat when opened in chimera). What would be the best way to make it 3D (i.e. assign/minimize the proper bond angles & distances). Is it possible to do it in chimera? Thanks

3 3

How can one create a sculpted view of a cellular assembly?
by Kenward Vaughan 21 Jan '15

21 Jan '15

As I slowly approach the operational realization of my computational/visualization lab, I received your lovely holiday card. This seriously elevated my dreams of larger scale possibilities for chemistry and biology students through the cellPack resources. In playing with it, I found myself wanting a more robust version of a clipping plane which would have several characteristics: 3 dimensional - take a cell and hack out 1 or 2 adjoining octants relative to its center to give a sliced view similar to what is found in some texts (also often used in geology books to depict the earth's internal structure--see http://www.the-science-site.com/earths-interior.html for an example). Being able to leave the nucleus alone while the surface strips away the cytoplasm would be cool, much like the referenced image above showing the core of the earth. Likely a difficult one to actually implement? Have the clipped surface not be flat, but rather allow a component to stay visible as the surface is moved past it until some determined relationship at which point it vanishes (perhaps when it has passed entirely through the clipping surface?). The clipped surface would thus show structure of a sort. Have the choice of the clipping surface being static with the assembly (it could rotate with the assembly, keeping the same internal view), or separately static (how a clipping plane behaves now--as I know it--where rotating the structure causes things to move through the unmoving plane. One approach I can imagine would be the use of a spherical or ellipsoidal (?) surface of specific size at a specific distance from the center of the cellular assembly. Is this already implemented in a way that can be manipulated to achieve at least the major aspects of this effect (e.g. logical ANDing of clipping by 3 intersecting planes)? Sorry for the long question... ;-/ 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

3 3

Ringer/Chimera/IDLE
by Alejandro Virrueta 21 Jan '15

21 Jan '15

Hello, I am modifying some Ringer files for my own needs, but I am running into a weird issue. If I run Ringer via typical means: $RINGER/ringer/ringer -i ringer_in_${pdb_ID}.txt -o ringer_out_${pdb_ID}.txt it works fine. However, when I step through it in what I think is my version of IDLE ("$CHIMERA_HOME/bin/chimera --nogui --nostatus --script `which idle2.7` /usr/bin/idle2.7", with `which idle2.7` returning "/usr/bin/idle2.7"), I run into this error: Traceback (most recent call last): File "/Users/av376/ringer-2.0/ringer/main_copy.py", line 233, in <module> main() File "/Users/av376/ringer-2.0/ringer/main_copy.py", line 161, in main move_set, pdb = structureFile(parameters, path) File "/Users/av376/ringer-2.0/ringer/main_copy.py", line 65, in structureFile pdb = chimera.openModels.open(parameters.pdbfileName) File "/Applications/Chimera.app/Contents/Resources/share/chimera/__init__.py", line 1817, in open checkForChanges=False, noprefs=noprefs) File "/Applications/Chimera.app/Contents/Resources/share/chimera/__init__.py", line 1644, in add makePseudoBondsToMetals(realMolecules) File "/Applications/Chimera.app/Contents/Resources/share/chimera/__init__.py", line 1961, in makePseudoBondsToMetals cmPBG = mol.metalComplexGroup(issueHint=True) File "/Applications/Chimera.app/Contents/Resources/share/chimera/__init__.py", line 119, in _getMetalPbg am.color = preferences.get(MOLECULE_DEFAULT, MOL_COMPLEX_COLOR) File "/Applications/Chimera.app/Contents/Resources/share/chimera/preferences/base.py", line 608, in get forPrefSave=forPrefSave) File "/Applications/Chimera.app/Contents/Resources/share/chimera/preferences/base.py", line 179, in get return self._options[name].get(forPrefSave=forPrefSave) File "/Applications/Chimera.app/Contents/Resources/share/chimera/preferences/base.py", line 80, in get return self.prefToVal(self.value) File "/Applications/Chimera.app/Contents/Resources/share/chimera/tkoptions.py", line 1229, in _prefToColor return getColorByName(pref) File "/Applications/Chimera.app/Contents/Resources/share/chimera/colorTable.py", line 87, in getColorByName color = chimera.MaterialColor(r/255.0, g/255.0, b/255.0) ValueError: missing default material Any suggestions? Cheers, Alex

2 1

Alternate conformations in residues
by Alejandro Virrueta 21 Jan '15

21 Jan '15

Hello all, Is there a way to get the all the information for an alternate conformation from a residue object? I have res.altLocs resulting in set([a,b]), but how can I get chi1-4 for conformation b for example? Thanks! Alex

2 1

Alternate conformations in residues
by Alejandro Virrueta 21 Jan '15

21 Jan '15

Hello all, Is there a way to get the all the information for an alternate conformation from a residue object? I have res.altLocs resulting in set([a,b]), but how can I get chi1-4 for conformation b for example? Thanks! Alex

1 0

Mesh defined by distance to selection
by R. Michael Sivley 21 Jan '15

21 Jan '15

I have a set of residues in a protein structure that each mark the center of a sphere over which a value has been calculated. I can select a subset of residues by value, but the values actually correspond to the sphere around that residue. I would like to visualize this by creating a mesh defined by being within 10A of any selected residue. If a selected residue isn't within 10A of any other selected residue, I would like this to produce two meshes, but if they are within 10A, I would like a single continuous mesh whose boundaries are defined as being 10A from either residue. Zone-selecting residues within 10A of the selection and then showing their surface as a mesh is a close approximation of what I'm trying to convey, but doesn't give me exactly what I'm looking for. I found this in the Chimera help pages, and it seems to be exactly what I need: http://www.cgl.ucsf.edu/chimera/data/tutorials/volumetour/volumetour.html#z… I think I'm missing a first step, though, because I don't understand how to actually create the mesh. I select my residues, open the Volume Viewer, and define the zone radius, but then I'm stuck. How do I actually instantiate the mesh given those values? I've tried showing the surface for the selected residues, but the changes I make in the Volume Viewer don't appear to have any effect on the surface, so I'm assuming these are two separate features. There is no "apply" button in the Volume Viewer, and closing the window does not create a mesh. I've included a screenshot of where I'm stuck, should it be helpful. Could you help me identify the missing step in this process? Thanks! Mike Sivley PhD Student in Biomedical Informatics Center for Human Genetics Research Vanderbilt University

2 1

UHBD map, generated with APBS
by Feldman-Salit, Anna 20 Jan '15

20 Jan '15

My current "strange" problem: I use APBS stand-alone to generate electrostatic potential for a protein in PQR format. What I do is to ask for output in UHBD (GRD not DX) format. Unfortunately, the generated GRD file is not read by CHIMERA, which upsets me... Strange, but VMD seems to be less sensitive to the GRD generated with APBS to visualize the map. Another point: I am very glad that CHIMERA computes potential. Do you think one could save it not only in DX format? Maybe you could suggest me a way to convert DX to UHBD format? Thank you! Have a lovely working day! Greetings, Anna --- Dr. Anna Feldman-Salit Molecular and Cellular Modeling (MCM) Group HITS gGmbH, Schloß-Wolfsbrunnenweg 35 69118 Heidelberg, Germany phone: +49-6221-533-260 fax: +49 6221 533298

2 1

Simply, keep up the good work
by Kevin D. Walker 20 Jan '15

20 Jan '15

I really enjoy using the Chimera program, and I am trying to convince my lab members to make the switch. I hope you get these types of emails all of the time. Your program is duly acknowledged during my seminars. ---------------------------------- Kevin D Walker Associate Professor Michigan State University Department of Chemistry and Department of Biochemistry & Molecular Biology 578 S Shaw Lane; Room 208 East Lansing, MI 48824 Phone: 517 355 9715 x 257 Fax: 517 353 1793 websites: 1) http://bmb.natsci.msu.edu/about/directory/faculty/kevin-d-walker/current-re…

2 1

Creating a cubic box for a density map
by Juan Luis Loredo Varela 18 Jan '15

18 Jan '15

Dear Chimera users When I use the Segment map option to save a specific fragment as .mrc file, Chimera uses the minimal box containing the segment. Does anyone know how to modify the dimensions of a density map to have it in a cubic box that can be used in Flex-EM? Thanks Juan -- Juan Luis Loredo Varela PhD Student Antson Group York Structural Biology Laboratory L0 Department of Chemistry University of York

3 2

Scripting image writing in Mac OSX
by Mike Sivley 16 Jan '15

16 Jan '15

I have a Chimera script that takes a set of attribute files and generates a collection of scenes images. This script is usually run on Scientific Linux using the Headless Linux 64-bit distribution of Chimera. I'm now trying to run this script on Mac OSX Mavericks and hitting issues with the GUI. If I pass the --nogui flag to Chimera, then the script fails, stating that the GUI is necessary to save images. If I remove the flag, the GUI starts up correctly, but scripted control stops once the new attributes are defined. The user is prompted by the Render by Attribute window, and the script will not progress until the user confirms the changes and closes the GUI. How can I get the script to run start-to-finish without stalling for user-interaction, the way it does in linux? Thanks! Mike Sivley PhD Student in Biomedical Informatics Center for Human Genetics Research Vanderbilt University

2 1

Re: [Chimera-users] dotted lines
by Elaine Meng 16 Jan '15

16 Jan '15

On Jan 16, 2015, at 8:46 AM, Thomas Evangelidis <tevang3(a)gmail.com> wrote: > Hi Elaine, > I don't know why, but none of the 3 approached works, even with the latest Chimera version. > […] Hi Thomas, Your pseudobonds are in stick style. You would have to change the stick to wire before doing those other things; in Chimera there are no “dotted or dashed” sticks, only wires (lines). E.g. in Selection Inspector, inspect: Pseudobond, change bond style to wire. That only changes the specific pseudobond(s) you have selected; if you first select one distance pseudobond, you can press keyboard up arrow to expand the selection to all the distance pseudobonds. Can also be done with setattr command, e.g.: setattr p drawMode 0 Hard to see on white background. Can change color with setattr too: setattr p color black When you first create distances they are normally already in wire style, which is why I didn’t mention the stick/wire issue before. I hope this helps, Elaine ----- Elaine C. Meng, Ph.D. UCSF Computer Graphics Lab (Chimera team) and Babbitt Lab Department of Pharmaceutical Chemistry University of California, San Francisco

2 2

vop add
by Fco Javier Chichón 16 Jan '15

16 Jan '15

Dear all I would like to merge several volumes. Not all the voxels of the two volumes fit, so when I used vop add command, Chimera add the densities instead of average them. Anyone know how to do that? Thanks -- ======================================================= Francisco Javier Chichón García Centro.: Centro Nacional de Biotecnología Dept.: Macromolecular Structure Lab.: S1-S0.1 Tel.: (+34) 91585-5497:4690 =======================================================

2 3

dotted lines
by Thomas Evangelidis 15 Jan '15

15 Jan '15

Greetings, I have created a figure with manually added Hbonds (Tools->Structure Analysis->Distances). Is there any way to make these lines dotted/dashed? thanks, Thomas -- ====================================================================== Thomas Evangelidis PhD student University of Athens Faculty of Pharmacy Department of Pharmaceutical Chemistry Panepistimioupoli-Zografou 157 71 Athens GREECE email: tevang(a)pharm.uoa.gr tevang3(a)gmail.com website: https://sites.google.com/site/thomasevangelidishomepage/ =============================================================== *Physics is the only real science. The rest are just stamp collecting.* *- Ernest Rutherford*

2 1

Distance measurement in density map
by Maarten Schledorn 14 Jan '15

14 Jan '15

Dear reader, I'm looking at EM tomography data summarised in *.mrc density files. Since there are no atoms defined in these files, the usual distance measurement tool is inappropriate (I suppose). Is there another way to assess the dimensions of the oligomerised protein I'm looking at? With my sincere thanks in advance, Maarten

2 2

shown/hidden state of axes/planes not saved with scene
by Oliver Clarke 14 Jan '15

14 Jan '15

Hi all, Planes seem to be shown by default when I load a scene, regardless of whether they were displayed when I saved the scene. This is not normally an issue, but it is a little problematic when preparing animations using many different scenes - I will often have a plane defined at the start so I can use it to align the x/y/z axes to the normal of that plane (e.g. aligning the fourfold axis of an ion channel), but I do not want the plane displayed during the animation, and I don’t want to have to delete it from each scene and re-arrange the animation timeline if I can avoid doing so. Would it be possible to alter this behavior such that chimera knows the displayed/undisplayed states of axes/planes/centroids when saving a scene or session? Best, Oliver Clarke.

2 1

New to Chimera
by Mike 12 Jan '15

12 Jan '15

This is my first day with Chimera. I'm trying to find how to draw a small molecule, minimize the structure of the small molecule and then have it interact with my protein of interest. However, I don't see any way for me to draw a small molecule (potential drug). Can you provide any direction ? Thanks, Mike

2 3

need guidence
by Maria Saqib 08 Jan '15

08 Jan '15

Dear user i am using chimera for some images of protien .I have a problem while labelling i want to change font colors , How can i label it with different font colors? Waiting for your reply Maria

2 2

VR in Chimera?
by Andrey Malyutin 07 Jan '15

07 Jan '15

Good morning, With recent resurgence in virtual reality technology, such as Oculus Rift and Samsung Gear VR, is it possible that we'll see support for VR headsets in future releases of Chimera? Thank you, Andrey

2 1

Distance Grid SES calculation available?
by Hurt, Darrell (NIH/NIAID) [E] 06 Jan '15

06 Jan '15

Hi Chimera friends, This question is probably for Tom Goddard, but I'm happy to hear from anyone. I'm really interested in using Tom's code for "distance grid SES" calculations as described here: http://www.cgl.ucsf.edu/chimera/data/surface-oct2013/surface.html The above page suggests that this code will be part of Chimera 2, but I'm wondering if it is available now for surface calculations. Too many of my surface calculations are failing for our production server using Chimera. This would help immensely. Molmap is a decent alternative in some instances, but for my purposes, distance grid calculations are going to be more MSMS-like, which will be a better solution (cf. the description on the page comparing the internal pockets formed by molmap vs MSMS surfaces). Is there some way to call the routine, especially in a script? Python maybe? Thanks, Darrell -- Darrell Hurt, Ph.D. Section Head, Computational Biology Bioinformatics and Computational Biosciences Branch (BCBB) OCICB/OSMO/OD/NIAID/NIH 5601 Fishers Lane, 4A31 North Bethesda, MD 20852 Office: 240-669-2741 Mobile: 301-758-3559 Web: BCBB Home Page<http://www.niaid.nih.gov/about/organization/odoffices/omo/ocicb/Pages/bcbb.…> Twitter: @niaidbioit<https://twitter.com/niaidbioit> , @NIH3Dprint<https://twitter.com/nih3dprint> Disclaimer: The information in this e-mail and any of its attachments is confidential and may contain sensitive information. It should not be used by anyone who is not the original intended recipient. If you have received this e-mail in error please inform the sender and delete it from your mailbox or any other storage devices. National Institute of Allergy and Infectious Diseases shall not accept liability for any statements made that are sender's own and not expressly made on behalf of the NIAID by one of its representatives.

2 2

Parsing a new 3D file format
by ros 05 Jan '15

05 Jan '15

Hello! I would like to ask for advice and starting directions to be able to parse a new file format in Chimera so I can use Chimera's tools and generate images. The file format that I want to be able to open is generated from the Atoms in Molecules program AIMAll ( http://aim.tkgristmill.com/ ) The AIMAll program of course can also do visualization and is very good at it. I am attaching an example of how a molecule looks. The problem is that if you want to work with big biomolecules (like DNA), the program is not very good to cloak/hide/select atoms so you have to click atom by atom to hide it and make some sense of the visualization. So, I will try to make a parser (or something that you would recommend) to open the file in Chimera and use Chimera's masking rules to select/deselect, hide/show atoms, which is very powerful. The text output generated by AIMAll has an XYZ list for each atom like this: Nuclear Charges and Cartesian Coordinates: ------------------------------------------------------------------------------- Atom Charge X Y Z ------------------------------------------------------------------------------- C1 6.0 0.0000000000E+00 2.3265459100E+00 0.0000000000E+00 C2 6.0 0.0000000000E+00 0.0000000000E+00 1.7647056500E+00 C3 6.0 0.0000000000E+00 0.0000000000E+00 -1.7647056500E+00 H4 1.0 1.6866221200E+00 3.4814514700E+00 0.0000000000E+00 H5 1.0 -1.6866221200E+00 3.4814514700E+00 0.0000000000E+00 C6 6.0 -2.0148478600E+00 -1.1632729600E+00 0.0000000000E+00 H7 1.0 -2.1717143500E+00 -3.2013833400E+00 0.0000000000E+00 H8 1.0 -3.8583364700E+00 -2.8006813000E-01 0.0000000000E+00 C9 6.0 2.0148478600E+00 -1.1632729600E+00 0.0000000000E+00 H10 1.0 2.1717143500E+00 -3.2013833400E+00 0.0000000000E+00 H11 1.0 3.8583364700E+00 -2.8006813000E-01 0.0000000000E+00 H12 1.0 0.0000000000E+00 0.0000000000E+00 3.8040484300E+00 H13 1.0 0.0000000000E+00 0.0000000000E+00 -3.8040484300E+00 That is easy, I think. But as you can see from the attached image, I want to be able to visualize the green spheres (which are called critical points and represents special properties of the electron density, extracted from quantum mechanical calculations). Those points are represented in the output file as: CP# 16 Coords = 5.97894312844235E-20 1.21974356425640E+00 -9.60854463947958E-01 Type = (3,-1) BCP C1 C3 Rho = 2.4067790431E-01 GradRho = 3.5016504866E-18 -2.7948129921E-14 1.7716730860E-14 HessRho_EigVals = -4.4278076226E-01 -4.3836581782E-01 3.5682199801E-01 HessRho_EigVec1 = 1.0000000000E+00 0.0000000000E+00 0.0000000000E+00 HessRho_EigVec2 = 0.0000000000E+00 -6.1321358320E-01 7.8991714842E-01 HessRho_EigVec3 = 0.0000000000E+00 7.8991714842E-01 6.1321358320E-01 DelSqRho = -5.2432458207E-01 Bond Ellipticity = 1.0071370222E-02 V = -2.5219393995E-01 G = 6.0556397216E-02 CP# 18 Coords = 1.05289204670721E+00 3.04112062239339E+00 1.96757752864282E-18 Type = (3,-1) BCP H4 C1 Rho = 2.9524693796E-01 GradRho = 1.0061396161E-16 -1.6479873022E-16 1.6660693347E-18 HessRho_EigVals = -7.9444808957E-01 -7.8708168641E-01 4.1039750300E-01 HessRho_EigVec1 = -5.7017191910E-01 8.2152539989E-01 0.0000000000E+00 HessRho_EigVec2 = 0.0000000000E+00 0.0000000000E+00 1.0000000000E+00 HessRho_EigVec3 = 8.2152539989E-01 5.7017191910E-01 0.0000000000E+00 DelSqRho = -1.1711322730E+00 Bond Ellipticity = 9.3591342377E-03 V = -3.7782718105E-01 G = 4.2522056404E-02 etc etc etc First line indicates the critical point number followed by the XYZ coordinate of the point, and the second line shows which atoms are involved between that critical point. The rest is the electronic properties of the critical point (nothing to visualize there) So, I am guessing that the critical points can be represented as dummy atoms in Chimera following the XYZ coordinates? The idea is that first the molecule coordinates are parsed and the molecule is displayed and then display the critical points of each atom. If possible, show lines joining atoms that have a critical point between them... I know it is a difficult task, but maybe you can provide some starting directions? Anything will help because right now is double clicking 1000+ atoms! Thank you so much and have a good day, Rodrigo. ------------------------- Rodrigo Galindo-Murillo, PhD Department of Medicinal Chemistry University of Utah 30 South 2000 East, Room 307 Salt Lake City, UT 84112-5820 http://home.chpc.utah.edu/~cheatham/ (801) 587-9652 (801) 585-6208 (Fax)

2 1

WebGL viewer
by Thomas Hrabe 05 Jan '15

05 Jan '15

Hi everyone, I am trying to implement a EM density viewer in a HTML page based on the chimera webgl & html export feature. The original page chimera exports when you save a webgl / html scene has been modified to load multiple JSON objects with densities as they come from chimera. The bar at the bottom of the page essentially performs a coarse density sliding similar to chimera. Almost all works great but what I am missing is how to connect the rotation that was applied to one model and to propagate it to all other models. What is happening now is that when one rotates one model and slides the bar, the next model will not be rotated but will be in it’s original position. Here’s my page I am currently working on http://localize.pytom.org/php/displayWebGL.php?jobID=tutorial I tried catching the mouse event with a loop, but it seems it fails to propagate the rotation to the other models. Check the block at line 916 in the code try { var rotMat = vsphere(mouse_position, new_position); var camera = this.camera; var matrix = new Mat4; matrix.$translate(camera.target[0], camera.target[1], camera.target[2]); matrix.$mulMat4(rotMat); matrix.$translate(-camera.target[0], -camera.target[1], -camera.target[2]); for (var i = 0, models = this.scene.models, l = models.length; i < l; ++i) { var elem = models[i]; elem.matrix = matrix.mulMat4(elem.matrix); } } Thank you in advance for your help, Thomas

2 1

mmaker usage with md-movie perframe
by Bala subramanian 02 Jan '15

02 Jan '15

Friends, I have a crystal structure of a protein (370aa). I have a trajectory of a similar protein (not same) with (210aa). I opened the crystal structure and the trajectory via md-movie tool. In the perframe analysis, i gave the following command (#1 my traj, #0 crystal.pdb) *mmaker #1 #0 computeSS false* mmaker uses different number of atom pairs for rmsd calculation. If there are say 150aa that match between my target and crystal, then mmaker should make the rmsd firt for 150 atom pairs in every frame. Why does it use different no. of atom pairs. Do i need to give some extra options ? Thanks, Bala -- C. Balasubramanian

3 4

Add reasonable bonds from command line
by Enrico Ronca 29 Dec '14

29 Dec '14

Dear Sir. I'm Enrico Ronca a post-doc researcher at the National Research Council of Italy in Perugia. I'm contacting you because I often work with metal atoms and not always chimera show me bonds. I solved the problem removing all bonds and pseudo bonds and adding "reasonable bond between selected atoms" in the adjust bonds menu of chimera. Is there a way to run the same procedure from the command line so that I could automate it by a script? Thank you in advance for your reply. I'm looking forward to hearing from you Enrico

2 1

Re: [Chimera-users] Customizing Chimera
by Eric Pettersen 25 Dec '14

25 Dec '14

On Dec 22, 2014, at 1:08 AM, Michał Kadlof wrote: > Hello, > > I'm currently work on proteins with disulfide bonds. I'm very often do the same sequence of tasks in Chimera: > > 1. open structure > 2. Select -> chemistry -> functional group -> disulfide > 3. extend selection to have selected all atoms of Cys-Cys bridge > 4. Show selected atoms > > Is there any way to wrap it into one-click button on toolbar? Hi Michał, While it is literally possible to put what you want on the toolbar, it's a bit easier to make a menu item instead. To do that, you create a command alias for the actions, and put that in a file that Chimera reads at startup so that it is always available. Details: The alias would be alias ^showbridge sel disulfide; sel up; disp sel This creates an alias named "showbridge" which is only expanded if typed at the beginning of a line -- that's what the ^ indicates. Such beginning-of-line aliases are also placed in an Aliases menu that appears after the Favorites menu (if you have any such aliases). To get the alias read by Chimera at startup you would put it in any of the files listed in the "Files to read at startup" section of the Command Line preferences (Favorites->Preferences). This is also described in more detail here: Command Line. Finally, if you really want it on the toolbar instead, that will involve making your own (simple) plug-in in Python. I can provide more detail about that if you like. --Eric Eric Pettersen UCSF Computer Graphics Lab

2 1

Can't open the CASTp site
by Lei Jinmin 18 Dec '14

18 Dec '14

Hi, I have to use CASTp tool for research sake, but I can't open the site anymore( http://sts-fw.bioengr.uic.edu/castp/index.php ). I remember it worked a couple of weeks ago. Does anyboday know when it begins working again? Thank you in advance! Best, Jinmin Lei Researh assitant Aalto University, Finland

3 2

constrain measurement of rotation to a particular axis?
by Elaine Meng 18 Dec '14

18 Dec '14

Begin forwarded message: > From: Tom Goddard <goddard(a)sonic.net> > Subject: Re: [Chimera-users] constrain measurement of rotation to a particular axis? > Date: December 18, 2014 11:07:57 AM PST > > I think the sensible thing is to report the hinge angles without trying to use a common rotation axis. I can think of different ways to define the common rotation axis (you have to align all the structure in some way, then say how much the observed rotation about the axis for a particular structure translates into a rotation about a slightly different common axis), all of them arbitrary and hard to explain, and even harder to compute, so I wouldn’t go that way. > > Tom > >> From: Oliver Clarke <olibclarke(a)gmail.com> >> Subject: [Chimera-users] constrain measurement of rotation to a particular axis? >> Date: December 17, 2014 2:36:45 PM PST >> To: Mailing List <chimera-users(a)cgl.ucsf.edu> >> >> Hi all, >> >> When using the “measure rotation” command, is there any way to constrain the rotation to be measured about a particular axis? >> >> I am thinking specifically of a case where I have multiple structures of a protein with different inter domain hinge angles - I can use chimera to measure the hinge angle, but for each pair of structures it chooses a slightly different rotation axis, whereas I would like to keep the same axis so I can directly compare the rotation observed between pairs of structures. >> >> Or perhaps there is there another way to do this that would make more sense? Suggestions are welcome... >> >> Best, >> Oliver. >>

1 0

How to change colour of H-bond between atoms
by Mahendra B Thapa 18 Dec '14

18 Dec '14

Dear Chimera Users, There is a hydrogen bond between, say, H20 of a ligand (LIG,res no: 240; this H-atom is attached with N atom of LIG) and O atom of GLU 120 residue. There are other hydrogen bonds too, but I am interested to change the color of that H-bond only to, say, red. I tried some like bondcolor red :240@H20 But, it is not working. Any suggestion might help me. Also, is there any way to change width of the H-bond? Thank you, Mahendra Thapa University of Cincinnati,OH

2 1

constrain measurement of rotation to a particular axis?
by Oliver Clarke 17 Dec '14

17 Dec '14

Hi all, When using the “measure rotation” command, is there any way to constrain the rotation to be measured about a particular axis? I am thinking specifically of a case where I have multiple structures of a protein with different inter domain hinge angles - I can use chimera to measure the hinge angle, but for each pair of structures it chooses a slightly different rotation axis, whereas I would like to keep the same axis so I can directly compare the rotation observed between pairs of structures. Or perhaps there is there another way to do this that would make more sense? Suggestions are welcome... Best, Oliver.

1 0

Re: [Chimera-users] Request for Source code
by Elaine Meng 16 Dec '14

16 Dec '14

Dear Jayaraman T, You can turn off helical correction and mass weighting using command options, for example with this command: define axis perHelix true helicalCorrection false massWeighting false See: <http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/define.html> I’m just trying to convince you this is the easy way! However, you already also have the source code in your download, as Eric described in his reply: <http://plato.cgl.ucsf.edu/pipermail/chimera-users/2014-December/010574.html> Best, Elaine On Dec 15, 2014, at 6:57 PM, Jayaraman T <jaai.nasa(a)gmail.com> wrote: > Yes Elaine, > > We can use Chimera but there are helical correction and mass weighting options exist. And am trying to get it done for all these helices without these options. > > Anyway let me follow up and write you back if there is a need of further assistance. > > I think it would help us more and thanks again for your kindness. > > Good day. > > On Tue, Dec 16, 2014 at 6:51 AM, Elaine Meng <meng(a)cgl.ucsf.edu> wrote: > Dear Jayaraman T, > Chimera download already includes the Python code, but why not just use Chimera? It seems like that would be much, much easier than trying to find and understand all the parts of the code to do it. You can script it with Chimera commands “define” to create helix axes cylinders and “angle” to measure the angle between any pair of those. > > See: > <http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/define.html> > <http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/angle.html> > > …and information on looping through structures: > <http://www.rbvi.ucsf.edu/chimera/docs/ProgrammersGuide/basicPrimer.html> > > I hope this helps, > Elaine > ----- > Elaine C. Meng, Ph.D. > UCSF Computer Graphics Lab (Chimera team) and Babbitt Lab > Department of Pharmaceutical Chemistry > University of California, San Francisco > > On Dec 14, 2014, at 6:25 PM, Jayaraman T <jaai.nasa(a)gmail.com> wrote: > > > Dear, > > I'm Jayaraman T pursuing PhD here in Pusan National University. Since we are interested to calculate the interhelical angle and axes for a set of proteins, we found Chimera is really helpful for getting us the results. And i would like to do certain analysis on the basis of these interhelical angles, so i kindly request you to provide me the source code and other requirements for calculating the angle between helices. It would be grateful to you and please help me out at your earliest convenience. > > Thanks > > > > -- > Regards > > > Jram | 201493155 > PhD Student > Bio molecular Engineering Lab > School of Chemical Engineering > Pusan National University. > 30 Jangjeon-dong, Geumjeong-gu, > Busan, South Korea - 609-735. > +8210 9906 3680 > jram(a)pusan.ac.kr >

1 0

Request for Source code
by Jayaraman T 15 Dec '14

15 Dec '14

Dear, I'm Jayaraman T pursuing PhD here in Pusan National University. Since we are interested to calculate the interhelical angle and axes for a set of proteins, we found Chimera is really helpful for getting us the results. And i would like to do certain analysis on the basis of these interhelical angles, so i kindly request you to provide me the source code and other requirements for calculating the angle between helices. It would be grateful to you and please help me out at your earliest convenience. Thanks -- Regards Jram | 201493155 PhD Student Bio molecular Engineering Lab School of Chemical Engineering Pusan National University. 30 Jangjeon-dong, Geumjeong-gu, Busan, South Korea - 609-735. +8210 9906 3680 jram(a)pusan.ac.kr

3 2

different behavior of the Win and OS-X version seen in case of surface calculations
by Marek Maly 15 Dec '14

15 Dec '14

Hello, I found already several times the differences between OS-X and Win version of Chimera (1.10.). Problem is regarding the ability to calculate molecular surface. There are some cases when Win version of Chimera 1.10 is not able to calculate molecular surface of some more complicated molecular structure (e.g. protein) but on the other hand the OS-X version of Chimera (again 1.10) is able to do it (although writes warning that there was some partial problems). Moreover when I create PY session of such molecule with surface calculated/visualized under OS-X and then I want to open this session under Win version it is opened but without the surface and writes many errors into the log - see below. In both versions I am using default settings (no special changes in parameters important for the surface calc.). How is this different behavior of the same program under Win and OS-X possible ? Do you have any recommendation how to eliminate this problem ? Thank you in advance, Best wishes, Marek Here is the err list after the attempt to load into Win Chimera 1.10 session containing "problematic" molecule with surface calculated under OS-X version of the Chimera 1.10 ***** C:\Program Files\Chimera 1.10\bin\mscalc.exe 1.400000 2.000000 0 MSMSLIB 1.3 started on Local PC Copyright M.F. Sanner (March 2000) Compilation flags Surface calculation failed, mscalc returned code 5. Surface calculation frequently fails for large, multi-chain structures. The calculation may be successful if the chains are treated individually, by using the "split" command before generating a surface. If splitting is not desired or the structure is already a single chain, changing molecular surface parameters in the Selection Inspector or (before surface creation) the New Surfaces category of Preferences may allow the calculation to succeed. More details are given at http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/surfprobs.html Traceback (most recent call last): File "gB_2G_S16_A.py", line 2669, in <module> restoreCoreModels() File "gB_2G_S16_A.py", line 2654, in restoreCoreModels restoreSurfaces(surfInfo) File "C:\Program Files\Chimera 1.10\share\SimpleSession\versions\v40.py", line 55, in restoreSurfaces elif len(customColors) != len(s.triangleData()[0]): TypeError: 'NoneType' object has no attribute '__getitem__' Error restoring core models: TypeError: 'NoneType' object has no attribute '__getitem__' File "C:\Program Files\Chimera 1.10\share\SimpleSession\versions\v40.py", line 55, in restoreSurfaces elif len(customColors) != len(s.triangleData()[0]): See reply log for Python traceback. Traceback (most recent call last): File "gB_2G_S16_A.py", line 3025, in <module> restoreRemainder() File "gB_2G_S16_A.py", line 3020, in restoreRemainder restoreSilhouettes(silhouettes) File "C:\Program Files\Chimera 1.10\share\SimpleSession\versions\v62.py", line 265, in restoreSilhouettes model = idLookup(modelID) File "C:\Program Files\Chimera 1.10\share\SimpleSession\versions\v23.py", line 162, in idLookup return globals.sessionMap[sesID] KeyError: 27136 Error restoring post-model state: KeyError: 27136 File "C:\Program Files\Chimera 1.10\share\SimpleSession\versions\v23.py", line 162, in idLookup return globals.sessionMap[sesID] See reply log for Python traceback. Warning: Could not restore surface color on surface model with id 0.0 because that surface was not restored. C:\Program Files\Chimera 1.10\bin\lib\site-packages\matplotlib\axes.py:2369: UserWarning: Attempting to set identical left==right results in singular transformations; automatically expanding. left=1, right=1 + 'left=%s, right=%s') % (left, right)) -- Tato zpráva byla vytvořena převratným poštovním klientem Opery: http://www.opera.com/mail/

3 2

Chimera Published Paper
by Kurt Degabriele 08 Dec '14

08 Dec '14

Dear sir/madame, I am a pharmacy student currently undergoing my 4th year pharmacy project at the University of Malta as part of my 5 year study course; I am currently doing research within the Med. Chem. field and would like to ask if you can send me a copy of the following paper so I can include it in my study: UCSF Chimera—A visualization system for exploratory research and analysisEric F. Pettersen, Thomas D. Goddard,Conrad C. Huang, Gregory S. Couch,Daniel M. Greenblatt, Elaine C. Meng andThomas E. Ferrin* Article first published online: 1 JUL 2004DOI: 10.1002/jcc.20084 Regards Kurt

2 1

Re: [Chimera-users] How to use a colour-valued attribute?
by Elaine Meng 04 Dec '14

04 Dec '14

On Dec 3, 2014, at 1:30 AM, Mungo Carstairs (Staff) <g.m.carstairs(a)dundee.ac.uk> wrote: > Thanks all. > Tom, your suggestion of comma- instead or or-delimiter improves speed several times! > Elaine, your point about specifying the chain is noted. > > Conrad, I didn't explain too clearly: my 'real' command is the 'non-contiguous' example (think Zappo colouring with some custom features colours on top). The 'contiguous' example was just what I was using to explore behaviour of setattr / defattr / rangecolor. > (Also by the way the REST API is working fine for me.) > > Instead of 'color' I am now using 'setattr' followed by 'rangecolor'. > I tried both 'setattr' and 'defattr' (passing a temporary file) and didn't see any significant performance difference between them, at the scales I am testing so far (up to 3 models of 1600 residues). > > So I have what I need for now. But I'm still curious to know if / how you get Chimera to 'act on' (display) a colour-valued attribute. > > Best regards, > Mungo Hi Mungo, Setting the color attribute with setattr or defattr does change the color immediately. The problem in what you sent earlier: setattr r color green #0:1-84.B … is that there is no residue attribute named “color". There are atom and molecule-model attributes by that name (setattr a color …) but the only per-residue color setting is ribbon color, name “ribbonColor” … thus you would want something like: setattr r ribbonColor green #0:54.a,76.a … and that would only color the ribbon, not the atoms/bonds. <http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/atom_spec.html#respr…> Actually I am mystified as to why this would be any faster than using the color command with the same atomspec (#0:54.a,76.a in the example above), but maybe the others who know more about the “guts” of Chimera can comment. Best, Elaine ----- Elaine C. Meng, Ph.D. UCSF Computer Graphics Lab (Chimera team) and Babbitt Lab Department of Pharmaceutical Chemistry University of California, San Francisco

2 1

Commands for multiscale models?
by Hurt, Darrell (NIH/NIAID) [E] 03 Dec '14

03 Dec '14

Hi Chimera friends, I was reading this post from 2009 that indicates that there are no Chimera commands for multiscale models: http://www.cgl.ucsf.edu/pipermail/chimera-users/2009-June/004009.html Is it still the case that there are no Chimera commands for manipulating the multiscale models? I saw the blue-tongue virus example script for doing something similar using Python itself: https://www.cgl.ucsf.edu/chimera/docs/ContributedSoftware/multiscale/blueto… I want to script the loading of an arbitrary VIPERdb model and then create multiscale model surfaces for every chain (not just the copies). Can it be done with Chimera commands? Or do I need to work up my Python-fu? Thanks, Darrell -- Darrell Hurt, Ph.D. Section Head, Computational Biology Bioinformatics and Computational Biosciences Branch (BCBB) OCICB/OSMO/OD/NIAID/NIH 31 Center Drive, Room 3B62B, MSC 2135 Bethesda, MD 20892-2135 Bethesda office: 301-402-0095 Rockville office: 240-669-2741 Mobile: 301-758-3559 Web: BCBB Home Page<http://www.niaid.nih.gov/about/organization/odoffices/omo/ocicb/Pages/bcbb.…> Twitter: @niaidbioit<https://twitter.com/niaidbioit> , @NIH3Dprint<https://twitter.com/nih3dprint> Disclaimer: The information in this e-mail and any of its attachments is confidential and may contain sensitive information. It should not be used by anyone who is not the original intended recipient. If you have received this e-mail in error please inform the sender and delete it from your mailbox or any other storage devices. National Institute of Allergy and Infectious Diseases shall not accept liability for any statements made that are sender's own and not expressly made on behalf of the NIAID by one of its representatives.

2 3

volume viewer step value
by Ludovic Renault 03 Dec '14

03 Dec '14

Hi, When I load different EM volumes into Chimera, depending on the size of the box the display step value changes from 1 to 4. Is it possible to force chimera to always load the 3d map with a step value of 1? Or at least change some parameters in preferences to have step 1 loaded with my routine volumes (up to 300 pixel boxes for example). Thanks, Ludo

2 1

show surface fails
by Hartzell Jr., Criss 02 Dec '14

02 Dec '14

I am struggling to show an electrostatic surface on the solvent accessible surface for several related proteins. It works fine in Pymol with APBS but with Chimera I continue to get "surface calculation failed, mscalc returned code 5". I have read posts going back at least 5 years about this problem and I have tried all of the various work-arounds (vdwdef, split, probe radius and vertex density). Has this problem been solved? I am using version 1.8.1 on a PC. Is there a way to identify which part of the protein is causing the problem? Criss Hartzell Professor Department of Cell Biology Emory University School of Medicine Atlanta, GA 30322 ________________________________ This e-mail message (including any attachments) is for the sole use of the intended recipient(s) and may contain confidential and privileged information. If the reader of this message is not the intended recipient, you are hereby notified that any dissemination, distribution or copying of this message (including any attachments) is strictly prohibited. If you have received this message in error, please contact the sender by reply e-mail message and destroy all copies of the original message (including attachments).

3 3

regarding distance measurement
by abhishek jamwal 02 Dec '14

02 Dec '14

Dear All, I have created a stereoview using chimera where I am showing distance between sidechain atoms of a protein. when I save an image of this stereoview (tiff, 600dpi), the connector lines showing the distances are fading, they are not of same quality in the image as in the original view ! Is there a way to save image without disturbing quality of the connector lines showing distances ? thnx abhishek

2 1

How to use a colour-valued attribute?
by Mungo Carstairs 02 Dec '14

02 Dec '14

Hi, This feels like a basic question but it has stumped me for a day now so here we go... I'm exploring ways of programmatically setting colours to find the most efficient method. This is for the 'custom coloured' case (residue-based colour schemes are not a problem). The simple version is to send as many 'color' commands as there are colours in my (Jalview) sequence: color #ffafaf #0:1.B|#0:5.B|#0:7-8.B|....;color #ffff00 #0:39.B|#0:44.B|#0:47.B|#0:77.B; This works, but doesn't seem to scale up well for large/many structures - I start to see Chimera using 100% CPU for minutes at a time. I'm trying instead to see whether setting a colour-value attribute is faster. I can set the attributes ok with commands like setattr r color green #0:1-84.B etc and confirm they are set with list residues spec #0:1-84.B attribute color but I can't figure out how to get Chimera to use the attribute value as the display colour. 'Render by Attribute' only seems to consider numeric (range-valued) attributes, not colour-valued. I thought perhaps 'my' colour was just hidden in the colour hierarchy, but removing other colours doesn't seem to help. What am I missing? Thanks, Mungo Carstairs Division of Computational Biology University of Dundee http://www.jalview.org The University of Dundee is a registered Scottish Charity, No: SC015096

4 3

Display atoms after multiscale model
by vamsee 02 Dec '14

02 Dec '14

Hello, I am working on a virus capsid and trying to figure out the buried surface area between subunits (similar to what ViperDB has. I have to do this for a few capsids and I am trying to automate this to the maximum possible extent. I have been able to follow the steps in this link Buried surface area calculation <http://www.cgl.ucsf.edu/pipermail/chimera-users/2006-December/001131.html> and manually calculate the information. Are there commands available to do the same? Also, what is the difference between the surfaces created using a multiscale model and surfaces computed from the atom-spec? Would this difference in the surface affect the calculation of the buried surface area? And if it does affect it which one be a more accurate representation? Thank you Vamsee

2 3