- Chimera-users - RBVI Mailing Lists

Python script for reading volume values at atom positions
by Kyle Morris 19 Jan '17

19 Jan '17

Dear chimera dev and users, Can anyone confirm with the latest chimera version (i’m running v1.11.2 build 41376) if this python script for measuring map values at atoms works for them? The script was reported previously on the users list here: http://plato.cgl.ucsf.edu/pipermail/chimera-users/2007-May/001587.html For me I get a python import error: from VolumeViewer import Data_Region, full_region, Rendering_Options Thanks for your insight! Best wishes, Kyle

3 6

beta sheet to loop
by chemocev marker 19 Jan '17

19 Jan '17

Hi I am looking at some region of the protein in chimera and it looks like a loop but the from the procheck analysis it shows a beta-hair pin like structure. In the electron density, its shows a hydrogen bonding between the main chain atom. Why is this difference and how to overcome on it. best Jiri

2 1

forget last question
by A G Szabo 18 Jan '17

18 Jan '17

Forget the last question regarding visualizing nucleotides as balls and sticks. Figured it out. ags

1 0

quick question
by A G Szabo 18 Jan '17

18 Jan '17

I am looking at a protein polynucleotide structure. How do I turn the rectangles that represent the individual nucleotide into balls and sticks? Thanks Arthur G. Szabo

1 0

Autodock vina + modeller + ATOM syntax error
by Zirpel, Bastian 18 Jan '17

18 Jan '17

Hi everyone, I'm using Chimera 1.11.2 on windows. I am experience following error after my workflow. I fetched a known pdb model from the servers and generated several structural models of similar sequences using the integrated MODELLER script. Since my pdb model contains a covalently attached cofactor in the active site, I let MODELLER include non-water HETATM into my structural models as well. Afterwards I tried to do substrate/product binding studies with the known pdb model as well as with the generated structural models using the AutoDockVina script. However, when I run the script with the known pdb model, everything works fine. But when I'm running the script with the structural models, I receive following error msg: "Parse error on line 8553 in file "receptor.pdbqt": ATOM syntax incorrect: "Ho" is not a valid AutoDock type. Note that AutoDock atom types are case-sensitive." It should be noted that this error does not occur if I won't include non-water HETATM by MODELLER. Makes sense, since having a look at the generated pdbqt file, the FAD lines are generating the error. " HETATM 8554 HO3' FAD 525 62.339 60.023 -29.788 1.00 0.00 0.238 Ho HETATM 8555 HO2' FAD 525 66.741 58.106 -28.406 1.00 0.00 0.397 Ho HETATM 8556 HO3B FAD 525 60.634 66.984 -32.274 1.00 0.00 0.221 Ho HETATM 8557 HO4' FAD 525 65.182 62.428 -29.602 1.00 0.00 0.217 Ho HETATM 8558 1H6A FAD 525 68.751 67.392 -38.565 1.00 0.00 0.157 HD HETATM 8559 2H6A FAD 525 68.156 65.783 -38.183 1.00 0.00 0.157 HD HETATM 8560 HO2B FAD 525 62.183 68.217 -34.408 1.00 0.00 0.253 Ho " I found that AutoDockVina needs a "HD" instead of a "Ho" entry at these positions? Is this correct? So my question is, why is AutoDock vina generating a pdbqt file in the first place which it cannot read properly in the next step? How can I prevent this? I cannot just change the pdbqt entries as AutoDockVina will overwrite it again. Thanks in advance for any help! Bastian

2 1

model panel: long 'clicking time'
by Julian von der Ecken 17 Jan '17

17 Jan '17

Dear Chimera Team, there is a long delay between on and off of the ’S’-checkbox (show button) in the model panel. I mean the model it self vanished instantly but the ckeckbox takes much longer to be set to on/off. Mac OS X 64-bit Version tested (1.11 2016-06-30, 1.11.2 2016-12-02; ): very slow. In older versions (1.10.2 2015-07-23) it was much faster. It is kind unfortunate when you click very fast through all models manually or you just want to alternatively show and hide one model very fast to see differences. Is it possible to set the speed somewhere/somehow? Thanks a lot for your advise Julian ————————————————————————————————— Julian von der Ecken, Ph.D. Max-Planck-Institute of Molecular Physiology Department of Structural Biochemistry Otto-Hahn-Str. 11 44227 Dortmund, Germany room: B2.27 tel: (+49) (0)231 133-2306 e-mail: julian.von-der-ecken(a)mpi-dortmund.mpg.de<mailto:julian.von-der-ecken@mpi-dortmund.mpg.de> institute: http://www.mpi-dortmund.mpg.de group: http://www.mpi-dortmund.mpg.de/forschungsgruppen/raunser

2 2

Transparency of surface colored with electrostatic potential file
by Sarah Piper 17 Jan '17

17 Jan '17

Hi, I was wondering if you can make a surface transparent, which contains an electrostatic potential file (.dx), calculated with APBS? It would be nice to be able to see the ribbon structure of my .pdb file within the coloured surface. I’m using mol maps or SES to create the surface (my surface calculations (MSMS) usually fails, probably because of too large .pdb files) and then color this surface (Surface Color) with the submitted potential file from APBS. If I then want to make the surface with the potential colouring transparent, it doesn’t work. It just works to make the ’normal’ surface color itself transparent (e.g. transp 70,s #1). Does anyone have an idea if I can make the potential file or surface with electrostatic potential transparent as well? Many thanks for your help. Best regards, Sarah

2 2

Visualization Martini coarse grained model in Chimera
by MOHD HOMAIDUR RAHMAN 17 Jan '17

17 Jan '17

Dear All It is possible to visualize martini coarse grained bead of proteins and lipid in Chimera? if yes please guide me for same. Thanks & regards Rahman ''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' *M*d* H*omaidur* R**ahman *( Research Scholar ) Lab No. 510, Computational Biophysics Lab. Department of Biotechnology Indian Institute of Technology-Madras Chennai-600 036, India *Mobile No = +91- 7845991785* '''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''

3 6

on Unit Cell function
by Smith Liu 17 Jan '17

17 Jan '17

Dear All, Once we load a PDB in chimera, in Tools-Higher order structure-Unit Cell, if we click "Make copies", the copies of the PDB would be created based on space group, with each copy has a specific # (0,1,2,3..., for example). Will you please let me know is any rule available for how the number of # for each copy of PDB was defined? It seems the PDB copies with consecutive specific # were not most closely related in the Chimera display. I am looking forward to getting a reply from you. Smith

2 1

Fixing the Frame size
by chemocev marker 17 Jan '17

17 Jan '17

Hi I am making some images and I would like to know how to fix the zoom size so that in all the images size of the secondary structure should be similar. best Jiri

2 1

Preference/default for H-bond representation
by Markus Heller 16 Jan '17

16 Jan '17

Hello, Is there a way to set the preference for H-bonds to be represented as dashed lines, not as solid lines? I don't seem to be able to find anything in the Preferences or the documentation. Thanks and Cheers Markus

3 5

virtualization project - UCSF
by Scholz, Sebastian 16 Jan '17

16 Jan '17

Dear Sir or Madam, we, the COMPAREX AG, are a German Reseller. For one of our clients we conduct an evaluation license and have a list of questions provided, in the answer we need your help. This involves essentially the question whether the software can be virtualized and, if so, on what terms. These are crucial questions for the future orientation of its IT strategy for our customer. The Software: (for every Software is one questionnaire) Chimera 1.x Please note, that our questions refer to the latest software version. I therefore kindly ask for a timely response to the questionnaire, which you can found in the appendix (in German and English), but not later than January 23rd 2017. Kindly forward our request to whom it may concern, if you are not the responsible contact person. Thank you for your support. I look forward hearing from you soon. Kind regards / Mit freundlichen Grüßen Sebastian Scholz International Sales Team ______________________________________ COMPAREX AG Software Inside Sales Blochstraße 1 04329 Leipzig Direct contact: Sebastian Scholz mobile: +49 341 2568 6643 email: Sebastian.Scholz(a)comparex.de<mailto:Sebastian.Scholz@comparex.de> COMPAREX AG | Blochstraße 1 | 04329 Leipzig | Amtsgericht Leipzig | HRB: 15064 | USt-ID: DE 141626867 | Vorstand: Mag. Wilfried Pruschak (Vorstandsvorsitzender, CEO) | Marc Betgem (CSO) | Dr. Thomas Reich (CFO) | Aufsichtsratsvorsitzender: Mag. Wilhelm Doupnik | phone: +49 341 2568 000 | fax: +49 341 2568 999 | email: info(a)comparex.de<mailto:info@comparex.de> | web: www.comparex.de<http://www.comparex.de/>

1 0

Docking license
by Farah Yousef 13 Jan '17

13 Jan '17

Hello, I need an assisdtance in what I faced in submitting after I filled the required info for DOCK 6. Pages cannot be opened though I checked my net connection and I know DOCK program is free of charge and available for academic projects. Hope I get help in obtaining the license for DOCK 6, please. Best wishes,Farah Yousef

3 2

Web services (Modeller, Blast, etc.) running
by Elaine Meng 13 Jan '17

13 Jan '17

Hello users, I believe the RBVI web services called from Chimera are now running, thanks to the efforts of my coworkers. I just now verified running Modeller, Blast (requires Chimera 1.11.2 or later), and sequence realignment with Clustal Omega or Muscle. 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

Deleting a surface on command line, not just hiding it
by Rebecca Swett 12 Jan '17

12 Jan '17

Hi all, I'm working on an animation, and I'd like to know if there's a way I can delete or close a calculated surface from the command line. For my purposes, I'd like to show the surface representation of the protein, remove it, then show a dynamics trajectory, and show the surface again at a separate point (while the trajectory is stopped). Currently, without being able to close the surface, it tries to recalculate the surface at each trajectory frame, which doesn't work for my purposes. Let me know if this is possible. Cheers, Rebecca This email message and any attachments are confidential and intended for use by the addressee(s) only. If you are not the intended recipient, please notify me immediately by replying to this message, and destroy all copies of this message and any attachments. Thank you.

2 2

Re: [Chimera-users] dipoles
by Sagar, Vatsala (NIH/NEI) [E] 12 Jan '17

12 Jan '17

Hi Tom, So if you align the pdb files of the rotated proteins, they have an exact match. But when you make the PQR files, it adds hydrogens differently to the rotated proteins, resulting in a charge difference in the said histidine. Should the program be adding hydrogens differently to the rotated proteins, and is there a way to remedy this> Thanks a lot, Vatsala Sagar From: "Sagar, Vatsala (NIH/NEI) [F]" <vatsala.sagar(a)nih.gov <mailto:vatsala.sagar@nih.gov>> Date: Friday, January 6, 2017 at 10:27 AM To: "chimera-bugs(a)cgl.ucsf.edu <mailto:chimera-bugs@cgl.ucsf.edu>" <chimera-bugs(a)cgl.ucsf.edu <mailto:chimera-bugs@cgl.ucsf.edu>> Subject: FW: [Chimera-users] dipoles Hi Tom, I am attaching the PQR files and raw PDB files of the crystal structure of human gamma D with two monomers rotated 90 degrees to one another. Please take a look at it and tell me what you think, Thanks, Vatsala From: Tom Goddard <goddard(a)sonic.net <mailto:goddard@sonic.net>> Date: Thursday, January 5, 2017 at 7:07 PM To: "chimera-users(a)cgl.ucsf.edu <mailto:chimera-users@cgl.ucsf.edu> BB" <chimera-users(a)cgl.ucsf.edu <mailto:chimera-users@cgl.ucsf.edu>>, "Sagar, Vatsala (NIH/NEI) [F]" <vatsala.sagar(a)nih.gov <mailto:vatsala.sagar@nih.gov>> Subject: Re: [Chimera-users] dipoles Hi Vatsala, I guess you are refering to the dipole.py script on the Chimera Python scripts web page: http://plato.cgl.ucsf.edu/trac/chimera/wiki/Scripts <http://plato.cgl.ucsf.edu/trac/chimera/wiki/Scripts> Yes it should give the same answer for a molecule and a copy of that molecule rotated by 90 degrees. Without have the structure and rotated PQR files I don’t know what is going wrong. Tom > On Jan 4, 2017, at 11:17 AM, Sagar, Vatsala (NIH/NEI) [F] wrote: > > Hi, > Do you know if there is a method to calculate dipoles in Chimera? I have generated a PQR structure in predicted protonation state and uploaded it onto Chimera. I then used a downloaded app called dimer.py to calculate the dipole moments of two of the same monomers rotated 90degrees. Technically they should both give the same dipole as the structures align 100% but according to the PDB coordinates are only rotated relative to one another. The dipole.py scipt is calculating them as two different dipoles. Can you help me with this? > Thanks, > Vatsala Sagar > _______________________________________________ > Chimera-users mailing list: Chimera-users(a)cgl.ucsf.edu <mailto:Chimera-users@cgl.ucsf.edu> > Manage subscription: http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users <http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users>

2 1

Re: [Chimera-users] dipoles
by Tom Goddard 10 Jan '17

10 Jan '17

Hi Vatsala, If you run the Chimera pdb2pqr command or graphical interface on a molecule and a rotated copy of the same molecule and it gives different protonation states for histidine that seems like a bug. As described in the Chimera documentation this uses the PDB2PQR web service provided by the National Biomedical Computation Resource (NBCR). https://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/pdb2pqr.html <https://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/pdb2pqr.html> You might want to contact the developers of PDB2PQR to pursue this matter further http://www.poissonboltzmann.org/support/home/ <http://www.poissonboltzmann.org/support/home/> Chimera can also add hydrogens and compute charges using the Chimera DockPrep tool so maybe that is an option for you to avoid this problem in pdb2pqr. Tom > On Jan 10, 2017, at 7:35 AM, Sagar, Vatsala (NIH/NEI) [E] <vatsala.sagar(a)nih.gov> wrote: > > Hi Tom, > So if you align the pdb files of the rotated proteins, they have an exact match. But when you make the PQR files, it adds hydrogens differently to the rotated proteins, resulting in a charge difference in the said histidine. Should the program be adding hydrogens differently to the rotated proteins, and is there a way to remedy this> > Thanks a lot, > Vatsala Sagar > > From: "Sagar, Vatsala (NIH/NEI) [F]" <vatsala.sagar(a)nih.gov <mailto:vatsala.sagar@nih.gov>> > Date: Friday, January 6, 2017 at 10:09 AM > To: Tom Goddard <goddard(a)sonic.net <mailto:goddard@sonic.net>>, "chimera-users(a)cgl.ucsf.edu <mailto:chimera-users@cgl.ucsf.edu> BB" <chimera-users(a)cgl.ucsf.edu <mailto:chimera-users@cgl.ucsf.edu>>, "Sagar, Vatsala (NIH/NEI) [F]" <vatsala.sagar(a)nih.gov <mailto:vatsala.sagar@nih.gov>> > Subject: Re: [Chimera-users] dipoles > > Hi Tom, > I am attaching the PQR files and raw PDB files of the crystal structure of human gamma D with two monomers rotated 90 degrees to one another. Please take a look at it and tell me what you think, > Thanks, > Vatsala > From: Tom Goddard <goddard(a)sonic.net <mailto:goddard@sonic.net>> > Date: Thursday, January 5, 2017 at 7:07 PM > To: "chimera-users(a)cgl.ucsf.edu <mailto:chimera-users@cgl.ucsf.edu> BB" <chimera-users(a)cgl.ucsf.edu <mailto:chimera-users@cgl.ucsf.edu>>, "Sagar, Vatsala (NIH/NEI) [F]" <vatsala.sagar(a)nih.gov <mailto:vatsala.sagar@nih.gov>> > Subject: Re: [Chimera-users] dipoles > > Hi Vatsala, > > I guess you are refering to the dipole.py script on the Chimera Python scripts web page: > > http://plato.cgl.ucsf.edu/trac/chimera/wiki/Scripts <http://plato.cgl.ucsf.edu/trac/chimera/wiki/Scripts> > > Yes it should give the same answer for a molecule and a copy of that molecule rotated by 90 degrees. Without have the structure and rotated PQR files I don’t know what is going wrong. > > Tom > > >> On Jan 4, 2017, at 11:17 AM, Sagar, Vatsala (NIH/NEI) [F] wrote: >> >> Hi, >> Do you know if there is a method to calculate dipoles in Chimera? I have generated a PQR structure in predicted protonation state and uploaded it onto Chimera. I then used a downloaded app called dimer.py to calculate the dipole moments of two of the same monomers rotated 90degrees. Technically they should both give the same dipole as the structures align 100% but according to the PDB coordinates are only rotated relative to one another. The dipole.py scipt is calculating them as two different dipoles. Can you help me with this? >> Thanks, >> Vatsala Sagar >> _______________________________________________ >> Chimera-users mailing list: Chimera-users(a)cgl.ucsf.edu <mailto:Chimera-users@cgl.ucsf.edu> >> Manage subscription: http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users <http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users> >

1 0

Map conservation of my sequence compared to sequence alignment
by mfellner＠msu.edu 09 Jan '17

09 Jan '17

Good day I like the mapping sequence conservation tool in chimera, however I ran into a problem that IÂ could not overcome with anything IÂ found online. I have the consensus mapped onto my protein but that is not what IÂ want. IÂ want to show how the residues of my protein are conserved. For example at a particular residue the consensus is a C but in my structure its a W. at the moment the W is colored to be highly conserved based on the consenus of the C. The W is maybe 15% conserved, how can IÂ map the conservation of my sequence compared to all the other sequences. IÂ hope IÂ explained what my issue is Thank you for this great tool Matthias Fellner

3 2

dipoles
by Sagar, Vatsala (NIH/NEI) [F] 06 Jan '17

06 Jan '17

Hi, Do you know if there is a method to calculate dipoles in Chimera? I have generated a PQR structure in predicted protonation state and uploaded it onto Chimera. I then used a downloaded app called dimer.py to calculate the dipole moments of two of the same monomers rotated 90degrees. Technically they should both give the same dipole as the structures align 100% but according to the PDB coordinates are only rotated relative to one another. The dipole.py scipt is calculating them as two different dipoles. Can you help me with this? Thanks, Vatsala Sagar

2 3

Mutation in Chimera
by Islam, Tanzila 06 Jan '17

06 Jan '17

Hi I am learning Chimera and I have gone through soe online tutorials. The tutorials showed how to introduce point mutation by Rotamers command. But it only allows to replace 1 ammino acid residue at a time. I want to insert 4-5 ammino acid residues in a gap of a sequence. Is there any standard procedure, command in Chimera that will let me do this type of residue insertion? This is really crucial for my research work. Thanks Tanzila Islam Graduate Student Washington State University

2 6

Problem opening DOI from Protein Data Bank
by Tones,Jason R 05 Jan '17

05 Jan '17

Hello, I am a student and recently downloaded chimera to use for a class, when I clicked to open the DOI of a search on the protein database I picked the wrong file under chimera to use to open it and it wouldn't open. I tried uninstalling and reinstalling and it keeps going to the file that won't open what I'm trying to do. Thank you for any help or suggestions you may have for me. Jason

2 1

Measuring hexamer asymmetry
by Tom Goddard 05 Jan '17

05 Jan '17

Hi Elena, Ok. The “measure rotation” command shows the rotation axis defined by the rigid motion of one subunit onto another — it has nothing to do with the hexamer axis. I think you are asking for is that if you take subunit 1, rotate it 60 degrees about the hexamer axis, then compare the orientation and translation of subunit 2 to that moved copy of subunit 1. Then do the same for subunits 3,4,5,6. To do this take subunit 1 and make 5 rotated copies by 60,120,180,240,300 degrees and save that as a new hexamer PDB file. Now open the asymmetric hexamer and this symmetric hexamer. If the subunits are chain ids A,B,C,D,E,F then use “measure rotation #0:.B #1:.B”, “measure rotation #0:.C #1:.C", …. This will still show a rotation axis that is not the hexamer axis which indicates the rotation needed to make the exact symmetry copy of subunit 1 rotate onto the asymmetric subunit 2, and it will give the translation, rotation axis and angle in the reply log. Maybe you are not interested in how the symmetric copies of subunit 1 rotate when being aligned to asymmetric subunits 2,3,4…. Then you can just ignore this rotation information and get the translation from the reply log. Tom > On Jan 5, 2017, at 7:38 AM, Elena Zehr wrote: > > Hi Tom, > > Thank you very much for you response. Yes within the asymmetric hexamer subunits have different conformations. To test your advice on RMSD, I first generated hexamer by using six copies of the same subunit to align to the original hexametric structure, then I closed the original hexametric structure. The asymmetric hexamer now consisted of subunits with identical conformation, please note the hexamer is still asymmetric. Then I opened a monomer with the same conformation as the subunits within the newly-generated hexamer, and calculated rotation and translation matrices between the first protomer within the hexamer and the rest of the subunits (#2-6). The central axis for each of the calculations deviated from each other. > > My issue is the following: to make the calculations reliable I wish to have the same central axis running through the center of the hexamer and not the axes between #1 and #2 or #1 and #4, etc. which deviate from each other due to the asymmetry of the hexametric structure. Is it possible? For example I wish to calculate the central axis of the whole hexamer and then use it as a reference to calculate rotation and translation for each subunit (#2-5) around the ring with respect to the subunit #1 and the central axis of the whole hexamer. Is it possible to define parameters for the central axis in the command line together with “measure rotation #1 #2 or #1 and #4” ? > > > Many thanks, > Elena > > > On Jan 4, 2017, at 4:47 PM, Tom Goddard wrote: > >> Hi Elena, >> >> Elaine asked me to look at your problem measuring orientation of subunits of a hexamer. Does your hexamer have cyclic symmetry? If so I think the reason the displayed axes don’t appear to align with the hexamer symmetry axis is because your alignment of one subunit to another is not good. You can get the root mean square deviation of C-alpha positions for the alignment by looking in the Reply Log (menu Favorites) after running the match maker command, for example, >> >> RMSD between 308 pruned atom pairs is 0.361 angstroms; (across all 308 pairs: 0.361) >> >> From your picture with blue and tan complexes it does not appear they are aligned at all. If the sequences of the subunits are very different then match maker may not be able to align them. >> >> Tom >> >> >> >> Begin forwarded message: >> >>> From: Elena Zehr >>> Subject: Re: user-defined rotation axis >>> Date: January 4, 2017 at 1:09:10 PM PST >>> To: Elaine Meng >>> >>> Dea Elaine, >>> >>> Thank you for getting back to me so soon. Sorry for sending you my question again. >>> >>> I am getting the following (see attached image) Where the central axis between protomer A and B is different than the one between protomer A and F. I wish to have the same central axis that would be the central axis of the hexamer, otherwise I assume my measurements are not accurate. I get the same if I use the hexamers or separate monomers >>> >>> >> <axis.tiff> >>> On Jan 4, 2017, at 1:15 PM, Elaine Meng wrote: >>> >>> Begin forwarded message: >>> >>>> From: Elena Zehr >>>> Subject: Re: user-defined rotation axis >>>> Date: January 4, 2017 at 1:09:10 PM PST >>>> To: Elaine Meng >>>> >>>> Dea Elaine, >>>> >>>> Thank you for getting back to me so soon. Sorry for sending you my question again. >>>> >>>> I am getting the following (see attached image) Where the central axis between protomer A and B is different than the one between protomer A and F. I wish to have the same central axis that would be the central axis of the hexamer, otherwise I assume my measurements are not accurate. I get the same if I use the hexamers or separate monomers >>>> >>>> >>> <axis.tif> >> >

1 0

Showing color key programatically colouring from a volume...
by Pablo Conesa 31 Dec '16

31 Dec '16

Dear all We have this simple script to open a volume (#0) coloured using values from a second volume (#1): open betaGal.mrc open extra/MG_Chimera_resolution.vol volume #1 voxelSize 3.54 vol #1 hide scolor #0 volume #1 cmap rainbow reverseColors True For scolor command we are using a volume. https://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/scolor.html#volume It would be nice to achieve something like this at start up (note the color key is based on the volume): I've did it dragging the mouse over the canvas while "Color key" utility is open. And, is there a posibility to round up the values. Some kind of: Psaudocode: Give the resulting color key for each round values display the color key I guess that if there is a way to get the color key, then rest can be done through commands. Cheers, Pablo.

2 3

color zone on command line
by Kyle Morris 24 Dec '16

24 Dec '16

Hi Chimera dev, Following on from this post back from 2009: http://plato.cgl.ucsf.edu/pipermail/chimera-users/2009-January/003514.html <http://plato.cgl.ucsf.edu/pipermail/chimera-users/2009-January/003514.html> How would one also split map in this python script? Thanks! Kyle

2 1

Re: [Chimera-users] Chimera-users Digest, Vol 164, Issue 26
by Healey, Joe 23 Dec '16

23 Dec '16

I might be able to offer a bit of help here, the Chimera team and Jaime helped me with basically this exact problem. I have a script here https://github.com/jrjhealey/bioinfo-tools/blob/master/strucfit.py that you could take from. You'll probably want to look specifically at lines 216-240. (Note - you'll also need to look at the relevant module imports for match and pychimera etc. The script will need Jaime's implementation of pychimera too (he's super helpful!) https://github.com/insilichem/pychimera The script gives output as below, as the script automatically finds the best matching PDB file (that's what the first 150 lines or so is looking for), there is additional information there. Best Model I'm Hit Comparing RMSD Assorted homology metrics 2k4r PAK_01796_model5.pdb 2.81843275676 37.5 13 0.00035 21.5 2k4r PAK_01796_model4.pdb 12.6869844869 37.5 13 0.00035 21.5 2k4r PAK_01796_model3.pdb 14.4133021316 37.5 13 0.00035 21.5 2k4r PAK_01796_model1.pdb 0.459177533295 37.5 13 0.00035 21.5 2k4r PAK_01796_model2.pdb 15.1861411912 37.5 13 0.00035 21.5 Hope that's useful! Happy to offer more help (though I'm sure the chimera team have it covered!) Joe M.Sc. B.Sc. (Hons) MSRB PhD Student MOAC CDT, Senate House University of Warwick Coventry CV47AL Mob: +44 (0) 7536 042620 | Email: J.R.J.Healey(a)warwick.ac.uk Jointly working in: Waterfield Lab<http://www2.warwick.ac.uk/fac/med/research/tsm/microinfect/staff/waterfield…> (WMS Microbiology and Infection Unit) and the Gibson Lab<http://www2.warwick.ac.uk/fac/sci/chemistry/research/gibson/gibsongroup/> (Warwick Chemistry) Twitter: @JRJHealey<https://twitter.com/JRJHealey> | Website: MOAC Page<http://www2.warwick.ac.uk/fac/sci/moac/people/students/2013/joseph_healey> | ORCID: orcid.org/0000-0002-9569-6738 ________________________________ From: chimera-users-bounces(a)cgl.ucsf.edu <chimera-users-bounces(a)cgl.ucsf.edu> on behalf of chimera-users-request(a)cgl.ucsf.edu <chimera-users-request(a)cgl.ucsf.edu> Sent: 22 December 2016 20:00 To: chimera-users(a)cgl.ucsf.edu Subject: Chimera-users Digest, Vol 164, Issue 26 Send Chimera-users mailing list submissions to chimera-users(a)cgl.ucsf.edu To subscribe or unsubscribe via the World Wide Web, visit http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users or, via email, send a message with subject or body 'help' to chimera-users-request(a)cgl.ucsf.edu You can reach the person managing the list at chimera-users-owner(a)cgl.ucsf.edu When replying, please edit your Subject line so it is more specific than "Re: Contents of Chimera-users digest..." Today's Topics: 1. Saving Output From Match Align (Muk, Sanychen) ---------------------------------------------------------------------- Message: 1 Date: Thu, 22 Dec 2016 16:11:00 +0000 From: "Muk, Sanychen" <smuk(a)coh.org> To: "chimera-users(a)cgl.ucsf.edu" <chimera-users(a)cgl.ucsf.edu> Subject: [Chimera-users] Saving Output From Match Align Message-ID: <A69D81CB04F0734E96302AA9E36F9B233E521BEE(a)ppwexch2kx03.coh.org> Content-Type: text/plain; charset="iso-8859-1" Hello, I'm working on a script to superimpose pdb structures onto one another and then run the match -> align tool to align them to get the RMSD between residues of the pdb structures. I know that from the multialign interviewer window, I can then go to structure acess match and then save the attributes as a text file that will contain the residue number and RMSD. Is there a command that I can use to output this attributes file automatically? Thank you, and happy holidays. Best, Sanychen Muk --------------------------------------------------------------------- *SECURITY/CONFIDENTIALITY WARNING: This message and any attachments are intended solely for the individual or entity to which they are addressed. This communication may contain information that is privileged, confidential, or exempt from disclosure under applicable law (e.g., personal health information, research data, financial information). Because this e-mail has been sent without encryption, individuals other than the intended recipient may be able to view the information, forward it to others or tamper with the information without the knowledge or consent of the sender. If you are not the intended recipient, or the employee or person responsible for delivering the message to the intended recipient, any dissemination, distribution or copying of the communication is strictly prohibited. If you received the communication in error, please notify the sender immediately by replying to this message and deleting the message and any accompanying files from your system. If, due to the security risks, you do not wi! sh to receive further communications via e-mail, please reply to this message and inform the sender that you do not wish to receive further e-mail from the sender. (fpc5p) ---------------------------------------------------------------------

1 0

Saving Output From Match Align
by Muk, Sanychen 22 Dec '16

22 Dec '16

Hello, I'm working on a script to superimpose pdb structures onto one another and then run the match -> align tool to align them to get the RMSD between residues of the pdb structures. I know that from the multialign interviewer window, I can then go to structure acess match and then save the attributes as a text file that will contain the residue number and RMSD. Is there a command that I can use to output this attributes file automatically? Thank you, and happy holidays. Best, Sanychen Muk --------------------------------------------------------------------- *SECURITY/CONFIDENTIALITY WARNING: This message and any attachments are intended solely for the individual or entity to which they are addressed. This communication may contain information that is privileged, confidential, or exempt from disclosure under applicable law (e.g., personal health information, research data, financial information). Because this e-mail has been sent without encryption, individuals other than the intended recipient may be able to view the information, forward it to others or tamper with the information without the knowledge or consent of the sender. If you are not the intended recipient, or the employee or person responsible for delivering the message to the intended recipient, any dissemination, distribution or copying of the communication is strictly prohibited. If you received the communication in error, please notify the sender immediately by replying to this message and deleting the message and any accompanying files from your system. If, due to the security risks, you do not wish to receive further communications via e-mail, please reply to this message and inform the sender that you do not wish to receive further e-mail from the sender. (fpc5p) ---------------------------------------------------------------------

2 1

ChimeraX daily builds available!
by Elaine Meng 22 Dec '16

22 Dec '16

Prerelease daily builds of UCSF ChimeraX are now available to the public! UCSF ChimeraX home: http://www.rbvi.ucsf.edu/chimerax/ Download: http://www.rbvi.ucsf.edu/chimerax/download.html Please keep in mind that ChimeraX features are currently quite limited, that ChimeraX may not work on older machines and/or graphics cards, and that due to the upcoming holidays, responses may be slow if you contact us. On behalf of the Chimera/ChimeraX team, 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

change default color maps
by Michał Kadlof 21 Dec '16

21 Dec '16

I found Palette Editor Tool, but I don't know how to use it. How can I apply chosen pallete in Surface Color Tool or load colors from palette into rainbow tool? -- pozdrawiam serdecznie Michał Kadlof <m.kadlof(a)cent.uw.edu.pl>

2 1

Crossfade and disp
by Daniel Asarnow 21 Dec '16

21 Dec '16

Hello Chimera experts, Here is another movie making question. I noticed that disp and ~disp trigger frames to be recorded, without a call to wait. For example: disp #2 movie crossfade 25 ~disp #2 disp #3 This would produce a movie that crossfades from model #2 to blank, with #3 suddenly appearing. To get crossfades between two ligands, I had to use: disp #2 movie crossfade 25 perframe "~disp #2" frames 1 disp #3 movie crossfade 25 Based on the documentation I wasn't sure if this was intended behavior or a bug. Did I miss a more elegant way to crossfade between two models? Thanks, -da

3 3

Multiscale models command
by Ian Holmes 20 Dec '16

20 Dec '16

Hello Chimera gurus, I am using multiscale models to view a virus capsid (PDB 2QA0) in Chimera. I would like to create a Python script to render the capsid as a multiscale model and then highlight (i.e. select) certain residues. However, I can't figure out the command-line syntax to create a multiscale model. The docs also suggest the Midas Python module as an alternative way of scripting Chimera in Python, but I haven't yet found that module. Here is the rough order of what I am doing from the graphical user interface: Open model 2QA0 Tools menu -> Higher-order structure -> Multiscale Models panel (The following actions are all in the Multiscale Models panel) Click "Make models" button in "Models from molecules and matrices" section Click "All" button in "Select chains:" section Click "Hide" button in "Act on selected chains" section, "Selected chains" subsection Click "All" button in "Select chains:" section (again) Select "Ribbon" from "Style" menu in "Act on selected chains" section (then back to the main view) Zoom out (with mouse wheel) Favorites menu -> Command Line Type "select :10,20,30" (amino acids I want to highlight) Here is what I've got so far in terms of equivalent Python code: import os from chimera import runCommand as rc rc("open 2QA0") # make multiscale model here? rc("select :10,20,30") # select some residues Your help would be appreciated - thanks! Ian

3 3

(no subject)
by Kartheek 19 Dec '16

19 Dec '16

Hi UCSF CHIMERA users, I am trying to make a movie in chimera using MD movie and pdb option. However I have been stumbled with the following error "Residue 256.C not in first model on line 20792 of /tmp/tmpvr4a2t.pdb". I do not have any end records, how can I rectify the error ?? -- Sincerely, P.Kartheek, Research Scholar, IIIT-H, Hyderabad, INDIA.

2 2

Surface clip coloring
by Michał Kadlof 16 Dec '16

16 Dec '16

How can color volume by value with binning several values into single color? I achived the effect partialy by opening several times the same map. set different thresholds and colors to each one and playing with Per-Model Clipping tool, but it have some drawbacks. Planes can't be at exact same level because I observe artifacts. So I tried to move slightly inner layers upp, but still I observe "staris" effect. I also tried with Surface color tool, but it uses interpolation and blends colors togather. In other words I need something like this: https://1.bp.blogspot.com/-5B5ueLHCBPY/T4II4IddvwI/AAAAAAAAB_w/rYwqsJNldY0/… instead of something like this: https://i.stack.imgur.com/hQmry.png -- pozdrawiam serdecznie Michał Kadlof <m.kadlof(a)cent.uw.edu.pl>

2 1

probably a bug
by Michał Kadlof 16 Dec '16

16 Dec '16

Hello, I visualize density of my data two different way with the same thresholds and observe different sizes of regions. Isn't it a bug? -- pozdrawiam serdecznie Michał Kadlof <m.kadlof(a)cent.uw.edu.pl>

2 1

Hide side chain
by Zhu, Xiaoyan 12 Dec '16

12 Dec '16

Hello, I am using UCSF Chimera to generate figures of protein crystal structures. Some crystals, if open in Chimera, automatically show the side chains of some amino acids. I want to hide those side chains in my figures. I wonder if there's an easy way to do it. Thank you! Have a good one, Xiaoyan Zhu Ph.D. Candidate Department of Biological Sciences, Marquette University

2 2

Multiple transparent objects
by Daniel Asarnow 12 Dec '16

12 Dec '16

Hi, I'm having some trouble with a movie of a morph between structures of the same protein bound to two different ligands. The basic idea is a closeup of the binding pocket, and as the morph progresses the initial ligand becomes transparent while the second one becomes opaque. The result is a similar effect to crossfade, only the morph animation is playing simultaneously. I've tried both the screen capture and ray tracing methods, with the following results: Screen capture - When a ligand is almost entirely transparent, a white silhouette is shown at the highest z-order, blocking what's behind. - There is a "pop up" effect when the ligands suddenly switch z-order at the beginning and end of my animation. - Single-layer and flat transparency options don't help. Ray tracing - Multiple transparency is fine. - Reflections and shadows are not made transparent. - E.g. when one ligand is 100% transparent, the reflections and shadows are still visible as white lines and dark spots. Maybe someone here has some insights on my transparency problems? A way to use crossfade instead would be even better. I can record movies with the two ligands separately and then crossfade them manually with ffmpeg. However, it would be nice to have this working entirely in chimera, especially for some more complicated animations I have in the works. Best, -da

2 1

.pyc session file
by Luca Pellegrini 12 Dec '16

12 Dec '16

Hi, How can I open in Chimera a saved .pyc session file? When I try, Chimera reports the following error message: Error while processing /Users/lp212/Dropbox/<filename.pyc>: .pyc files are not portable; please use .py file instead Unfortunately I don’t have the .py file, Chimera doesn’t seem to have saved it (not sure why). Thanks, Luca Luca Pellegrini, PhD Department of Biochemistry University of Cambridge 80 Tennis Court Road Cambridge CB2 1GA UK Email: lp212(a)cam.ac.uk Phone: 0044-1223-760469

3 2

Selective lighting
by Oliver Grant 12 Dec '16

12 Dec '16

Hi all, Is it possible to use the lighting command and have it apply only to particular models? I'd like my ligand to "pop" out by having a light that only shines on it, and not the protein. Oliver

2 2

Fit-to-window for procedurally generated images
by Mike Sivley 09 Dec '16

09 Dec '16

I am procedurally generating images of protein structures using a Chimera python script. I can center the protein structure within the viewing window, align the protein to emphasize certain elements, set the output size of the final image, and manually apply a scale factor to the view, but I have not found a way to ensure that the protein structure is maximized (but fully contained) within the final image. Is there a "fit-to-window" command that will automatically scale the view in this way? Thanks! Mike Sivley

3 2

How to find symmetry parameters for icosahedral molecule
by Haselbach, David 08 Dec '16

08 Dec '16

Hi everyone, sorry to ask such a basic question. I have a cryo EM density map 3.2 A of a icosahedral molecule and I already build the model of the monomer. I would like to place the remaining 59 copies and calculate the BioMT matrix for it. I started playing around with the sym command, but I couldn't get it working nicely. It's probably a problem of the coordinate system, but I couldn't figure it out. Than I thought there must be an automatic way and I tried Multifit but also this failed me miserably as the runs fail continuously with "stderr". So I was wondering is there a protocol to do this? Already thank you for any suggestions, David Dr. David Haselbach Max-Planck-Institute for biophysical Chemistry Department for structural Dynamics Am Fassberg 11 37077 Germany Tel. +495512011302

2 1

symmetry axis
by chemocev marker 08 Dec '16

08 Dec '16

Hi All I am interested to measure the symmetry axis of individual sub-unit (chain A & chain B) along with the symmetry axis of the heterodimer (AB). Each chain has 2 domains and 2 fold rotations axis, and I can measure by removing 1 chain and measure for the other. Is there way to measure the inter-domain symmetry axis of the heterodimer molecule. best Jiri

2 3

Fwd: chimera minimize option
by Jacqueline Vitali 07 Dec '16

07 Dec '16

It puckers ASN and ASP (and probably other amino acids). (I run it to relieve steric conflicts only) How do you fix this? And can you run it with the whole protein selected?

4 5

question
by A G Szabo 06 Dec '16

06 Dec '16

You have been very helpful in explaining different items that use the command line. I have a PDB file that was determined by NMR. Thus there are a large number of chains number as pdb (#0.1) Chain A etc. So I know that I can select a single chain, and build up an image with Chains A, B, etc using the #0.1 chain number. I want to visualize three Gly residues in one chain. So I learned that I can use the command ribbackbone to visualize the backbone atoms of the amino acids that I want to show. So when I use that command all the backbone atoms of the residues selected are visualized including the side chains. I also know that if I use the command ~ribbackbone all the backbone atoms are suppressed leaving only the atoms of the side chains. I know that without any residues or atoms specified the command applies to the whole model. Now as I said that I want to visualize three Gly residues in one of the chains. The Gly residue numbers are 33, 37, 38. They are on Chain #0.1 K. The information on the ribbackbone in the Chimera users guide, indicates that I should be able to do this by using parameters termed atom-spec. i.e. ribbackbone atom-spec <frameatom_spec.html> I looked elsewhere for how to designate the parameters for atom-spec and tried a few things, but I was not able to achieve the specificity of showing the backbone atoms of only the Gly residues. I don't think that I have to indicate that I want only the Gly residues in Chain K #0.1 because I can select the Gly residues from the sequence of Chain K #0.1. So after this long description of what I want to do, would you kindly inform me of how I can visualize the backbone atoms of only Gly. Thank you Arthur G. Szabo

2 8

Re: [Chimera-users] Chimera -> Povray
by Elaine Meng 05 Dec '16

05 Dec '16

Hi Lothar, If you use "File… Export Scene” (instead of Save Image) or the “export” command, there is an option to save a file for POV-Ray without rendering: <http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/export.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 5, 2016, at 3:07 PM, esserlo(a)helix.nih.gov wrote: > > Hi, > currently, I am using chimera to make a figure for publication and render > the result in povray. However the figure is a bit complex - with a gradient > background just to mention one thing- so that I have to edit the .pov file. > Is there a way to ask chimera to output the .pov file but not render it in > povray right away (that would save time as I need to edit it first!). > > Unfortunately in the "save image" menu, using povray, I did not find a > checkbox that would enable chimera to write out .pov only. > Thanks for any advice. > Lothar >

1 0

difficulty in using chimera for large PDB files
by Shruti Khare 05 Dec '16

05 Dec '16

Dear Sir/Madam, When we use chimera to visualize PDB files with ~18,000 atoms (PDB ID 5FUU) even rotating the molecule is slow in chimera. The computer has 6GB RAM. Would upgrading RAM help or we need to get a graphics card installed? Graphics card capacity should be 2GB or higher? Could you please tell me the minimum computer configuration required by chimera to handle large PDB files? Also what is chimera's maximum limit for handling PDB files given a computer with very good configuration? Thank you. -- Shruti Khare Graduate Student Prof. Raghavan Varadarajan Lab Molecular Biophysics Unit Indian Institute of Science Bangalore +91-80-2293-2612 -- This message has been scanned for viruses and dangerous content by MailScanner, and is believed to be clean.

2 1

pipes and planks
by Kyle Morris 02 Dec '16

02 Dec '16

Dear Chimera dev, I am using the pipes representation. In this case is it still possible to color/highlight a segment by residue number? Using the normal commands when displaying pipes doesn’t update the color, only on the original model. Thanks for your help! Best wishes, Kyle

2 3

Selected measure angles
by Soumya Govinda Remesh 01 Dec '16

01 Dec '16

Hello, I have just modified the loop script available here ( http://www.rbvi.ucsf.edu/chimera/docs/ProgrammersGuide/basicPrimer.html) to measure certain angles and distances. I just want to know if there is a way to export the individual values of angles and distances to a file other than parsing from the reply-log. Thank you. Best, Soumya Govinda Remesh, Ph.D. Postdoctoral Fellow at Advanced Light Source Lawrence Berkeley National Laboratory 1 Cyclotron Road MS 6R2100 Berkeley, CA 94720 (phone) 510-495-8179 (cell) 804-402-8730 sgremesh(a)lbl.gov

2 2

saving pdb with transformed coordinates
by Michael Elbaum 01 Dec '16

01 Dec '16

Dear Chimera people, I'm trying to create a series of pdb files showing a single structure rotated around a given axis. Using the GUI I can turn the molecule and save the pdb after unclicking the "Use untransformed coordinates" box. However, when I run a script using the write command I don't know how to invoke this option. The docs on the "relative" option suggest that the default is to write transformed coordinates. I've tried the following: for i in xrange(0,90,10): rc("rock " + "y " + str(i) + " 1") fn = basefn + str(i) + '.pdb' rc("write format pdb 0 " + fn) and the files end up identical although the rotation around y does appear on the screen. I had the problem in version 1.10.2 and it remains after upgrading to 1.11. thanks for your help, Michael

3 3

Unit cell
by Neha Gandhi 01 Dec '16

01 Dec '16

Hi Elaine, I have tried to look at the unit cell information of pdb Id :3HR2. It is a fiber diffraction structure of collagen with only CA coordinates. I tried using different programs like VMD to visualise unit cell, the unit cell is at the beginning of N terminus in agreement with experimental paper. Whereas in Chimera it is not the case. How can I visualize the unit cell of 3HR2 structure in Chimera and also number of molecules within one unit cell? Many thanks, Neha -- Regards, Dr. Neha S. Gandhi, Vice Chancellor's Research Fellow, Queensland University of Technology, 2 George Street, Brisbane, QLD 4000 Australia LinkedIn Research Gate

2 5

Re: [Chimera-users] failing to open saved PDBs
by Healey, Joe 30 Nov '16

30 Nov '16

Hi Elaine, OK, that is what I had expected to happen yes - and that NMR ensemble looks to be what I was anticipating. I've attached the PDB chimera created here. I'm no expert on PDB internal formatting but it at least looks to me like there are coordinates, sequence information etc etc. I had thought that I had taken care to make sure whole models were saved (not just sel atoms), but I'll replicate it to make sure. Applying the preset in any event didn't work. Maybe you can spot what's wrong with it? Joe Healey M.Sc. B.Sc. (Hons) PhD Student MOAC CDT, Senate House University of Warwick Coventry CV47AL Mob: +44 (0) 7536 042620 | Email: J.R.J.Healey(a)warwick.ac.uk Jointly working in: Waterfield Lab<http://www2.warwick.ac.uk/fac/med/research/tsm/microinfect/staff/waterfield…> (WMS Microbiology and Infection Unit) and the Gibson Lab<http://www2.warwick.ac.uk/fac/sci/chemistry/research/gibson/gibsongroup/> (Warwick Chemistry) Twitter: @JRJHealey<https://twitter.com/JRJHealey> | Website: MOAC Page<http://www2.warwick.ac.uk/fac/sci/moac/people/students/2013/joseph_healey>

2 2

more questions about visualizing the backbone.
by A G Szabo 30 Nov '16

30 Nov '16

Last week I asked a question regarding about how to visualize the backbone atoms. I have figured that out. I am working with a peptide that forms an amyloid structure. The peptide strands form a beta sheet. When I select two of the peptide units that are hydrogen bonded to form two strands in the beta sheet, the strands are shown as a coil type structure rather than a flat beta strand. How can I change the coil backbone to a beta strand structure? Second part is there a simple way to suppress the side chains but leaving the backbone atoms? I have included two peptide units that are hydrogen bonded to one another. Thanks Arthur G. Szabo

2 7