- Chimera-users - RBVI Mailing Lists

Chimera 1.6 officially released
by Eric Pettersen 02 Mar '12

02 Mar '12

Hello everybody, A new production release of UCSF Chimera (version 1.6) is available: <http://www.cgl.ucsf.edu/chimera/download.html> Download is free for noncommercial use. Platforms: Windows, Mac, Linux, including 64-bit versions (recommended for working with large datasets on machines with at least 4GB memory). This will be the last Chimera release supporting Mac PowerPC platforms and OS X 10.5 (Leopard). New features include a Rapid Access interface to frequently used data and tools, new ribbon spline options, improved transparency rendering, new fetch type: Uniprot protein sequence and annotations, and several options for fitting atomic structures into density maps: global search, iterative sequential fitting, symmetric fitting, etc. More details are given below; see release notes for the full list: <http://www.cgl.ucsf.edu/chimera/docs/relnotes/1.6.html> Graphics, Depiction: * new spline options for ribbon closer to backbone atom positions * new initial model colors, pastel rather than saturated * improved transparency rendering (global single-layer on by default) * 3D labels drawn on top by default * fancy backgrounds: images from files, color gradients from Palette Editor * options for cleaner pocket/tunnel surface display (see "sop" command) Structure Analysis and Modeling: * AddH allows specifying less common states of Asp, Glu, Lys, Cys * FindHBond can label H-bonds with distances * Blast Protein shows more info for PDB hits (ligands, resolution, etc.) * Build Structure improvements for joining peptides, swapping substituents * tools to change chain IDs, renumber residues * minimization: - can be conjugate gradient (previously steepest descent only) - parts of structures can be ignored * distance and angle plotting for trajectories * Cage Builder tool for generating polyhedral cage models <http://www.cgl.ucsf.edu/chimera/docs/ContributedSoftware/cagebuilder/cagebu… > * interface to MultiFit web service (simultaneous rigid fitting of multiple structures to density, developed by the Sali Lab, UCSF) * "fitmap" command: - global search with random initial placement - iterative sequential fitting - symmetric fitting * "measure" command options: - map symmetry - surface contact area - center of mass User Interface and I/O: * Rapid Access interface to frequently used data and tools shown at startup, can be hidden/shown at any time by clicking the lightning-bolt icon <http://www.cgl.ucsf.edu/chimera/1.6/docs/UsersGuide/rapid.html> * caching fetched data is turned on by default * scaling with mouse/touchpad scrolling can be turned off (Mouse prefs) * Model Panel: - group/ungroup function, such as for NMR ensembles - functions reorganized into favorites/all * new input types: SPHGEN spheres, Maestro/Glide docking output (ViewDock) * WebGL export (experimental feature) Sequence-related: * new fetch type UniProt: protein sequence and feature annotations * Multalign Viewer: - UniProt features can be mapped onto already open sequences - Region Browser improvements for handling feature annotation regions - calculation of region overall RMSD - interface to Modeller homology modeling allows including water, other HET Enjoy! On behalf of the Chimera team, --Eric Eric Pettersen UCSF Computer Graphics Lab http://www.cgl.ucsf.edu

3 3

Re: [Chimera-users] reading SMILES description into Chimera
by Elaine Meng 01 Mar '12

01 Mar '12

Hi Boaz, You can just use the command line, for example: open smiles:CC1=CC=C(C=C1)S(=O)(=O)[O-] I'm thinking it may be a mistake that SMILES shows up as a file type in the "Open File" dialog. The only other place I can think of is in Build Structure (under Tools... Structure Editing), where the Start Structure section includes a SMILES option. In that case you would just type or paste in the SMILES string in the indicated field. 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 Mar 1, 2012, at 7:08 AM, Boaz Shaanan wrote: > Hi, > > I'm trying to read a small molecule for which I have a SMILES description into Chimera. I selected 'smiles' from the input options but it's not clear to me how to proceed. Do I have to generate a text file with the smiles description in it? do I paste the smiles description in the file name line? or what? > > Thanks for your help. > > Boaz > > > > Boaz Shaanan, Ph.D. > Dept. of Life Sciences > Ben-Gurion University of the Negev > Beer-Sheva 84105 > Israel > > E-mail: bshaanan(a)bgu.ac.il > Phone: 972-8-647-2220 Skype: boaz.shaanan > Fax: 972-8-647-2992 or 972-8-646-1710 > > > > > _______________________________________________ > Chimera-users mailing list > Chimera-users(a)cgl.ucsf.edu > http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users

3 3

how to copy partial charges from a .mol2 file to another
by Thomas Evangelidis 29 Feb '12

29 Feb '12

Dear Chimera community, Is there any straight forward way to copy the partial charges from one .mol2 file to another? Thanks in advance. Thomas -- ====================================================================== Thomas Evangelidis PhD student Biomedical Research Foundation, Academy of Athens 4 Soranou Ephessiou , 115 27 Athens, Greece email: tevang(a)bioacademy.gr tevang3(a)gmail.com website: https://sites.google.com/site/thomasevangelidishomepage/

2 3

partial conservation coloring
by Eric Pettersen 28 Feb '12

28 Feb '12

Jill Voreis wrote: > I'm trying to use chimera to show conservation maps. Using pdb files > that I fetch from the database using the pdb id number and > alignments created in clustal that contain the sequence of the pdb > file, I select render by conservation. For some reason, the > conservation colors only go onto my model when it's in ribbon form. > When I try to show the surface representation or use spheres so that > I can use commans to show the SAS, the colors of my model do not > match the conservation model but rather are stuck on some sort of > coloring scheme related to the compounds that are there. The > majority of the map looks brown, like the default color before > rendering by conservation. If you could help me, I'd really > appreciate it! Hi Jill, There are a few possibilities here, though it's impossible to be certain what the problem without access to a concrete example. My best guess is that much of the surface region of your protein is poorly conserved and those areas have columns where the rows are more than 50% gap characters. For such columns no conservation score is assigned. That means that when you use Render By Attribute with default settings, those atoms / surface parts will retain their original colors. If in Render By Attribute you use the "No-value color" color well to assign a color to such regions (perhaps gray) you will be able to tell if this is what is happening. Two other possibilities are that you don't have the "Color atom" and "Color surfaces" check boxes checked [though I image you would have noticed this!] or that your protein has multiple chains and that only one is associated with the alignment and only that one is getting colored. If none of these are "it", if you can provide me with a specific example so I can reproduce the behavior then I will investigate. --Eric Eric Pettersen UCSF Computer Graphics Lab http://www.cgl.ucsf.edu

1 0

Not crossfade but coinciding bond command...
by Kevin D. Walker 28 Feb '12

28 Feb '12

I temporarily disabled a bond command and found this to be the cause of the selected bond problem I am having. I don't think crossfade has anything to do with it.

1 0

Animate a Reacton Mechanism?
by Kevin D. Walker 28 Feb '12

28 Feb '12

I constructed a movie for the first part of a reaction mechanism (substrate docking). I would like the animated mechanism to continue as follows. The substrate donates a proton to a tyrosine residue of the enzyme, forming an intermediate. The intermediate rotates 180 (I have accomplished this part) and the proton rebounds to the intermediate forming a new product, which is released. I can animate the product release -- it's the reverse of substrate docking. The part of the animation I cannot envisage is how to smoothly show or morph the bond transitions (bond breaking and making) during the reaction progression. IAny advice on whether Chimera or other ancillary software can do this is appreciated.

4 6

using CalcAttr from a Python script
by Roger L. Chang 27 Feb '12

27 Feb '12

Hi, 1) I need to get residue.areaSES for two separate selections (all polar amino acids and all nonpolar amino acids). I have successfully generated a list of residue.areaSES using the GUI simply by selecting the residues of interest and then using the attribute calculator typing "residue.areaSES" in the formula box of the Attribute Calculator, which outputs the list to the Reply Log. This is the data that I need, but I need to get it for >2000 PDB structures and would like to automate this analysis through a Python script. I am somewhat new to Python and Chimera and don't yet know all the ins and outs yet. I dug through the code and found a method "getResidueAttr" in the "ResidueDataWrapper" class from "CalcAttr". This seems to be what I need, but I don't know the proper syntax to call it from my Python script. Could someone please provide an example? 2) I would like to sum the residue.areaSES values once I have them. I tried "sum(residue.areaSES)" in the formula box of the Attribute Calculator in the GUI, but this returns the error "Formula evaluation error: expecting vector argument in function". I'm not sure again on the proper syntax to use "sum()" in the Attribute Calculator. If I knew how to generate/save the residue.areaSES data from the command line, I could just parse the saved file and sum the values in Python. How can Attribute Calculator formulas be used from the command line, and if they can't, how can I save the residue.areaSES data to file from a Python script? Thanks in advance, Roger -- Roger L. Chang Graduate Student Bioinformatics and Systems Biology Graduate Program University of California, San Diego

3 2

Coloring molmap surfaces
by Tom Goddard 24 Feb '12

24 Feb '12

> Greetings Tom > > How can color the molmap surface of certain atoms/residues. > > Many Thanks, Vijay > Hi Vijay, Use the Color Zone dialog to color a molmap surface to match nearby atom colors (menu Tools / Volume Data / Color Zone). You select the atoms you want to use for color (e.g. menu Select / Select All) choose a distance range from those atoms to color in the Color Zone dialog and press the Color button. In Chimera 1.6 (release candidate available on the web) or daily builds there is a command version of color zone. For example scolor #0 zone @ca range 5 says to color molmap surface #0 to match C-alpha atom colors within a distance of 5 Angstroms. To avoid jagged edges in the coloring you should use the molmap gridSpacing option to specify a grid size about 1/5 to 1/10 the map resolution. The default value is 1/3 of map resolution making a fairly coarse surface mesh. Colors are only assigned to the mesh vertices, so a fine mesh is needed for smooth boundaries between different color regions. Tom

1 0

batch peptide builder
by Xiao-Ping Zhang 23 Feb '12

23 Feb '12

Dear Chimera developers, I would like to generate some peptides from a protein sequence. The peptide builder in Chimera is a very nice application. I wonder if there is a stand alone script which I can use to generate series peptides from the command line. Thank you for your help. Xiao-Ping Zhang

2 1

Drawing arrows without using .bild file
by Victor Muñoz 21 Feb '12

21 Feb '12

Dear all, I'm trying to make a script that draws an arrow on every atom and every frame during a trajectory using the MD movie. Right now, the script creates a .bild file to draw the arrows, but it is not a very good approach since it is creating a new file for every frame of the trajectory. I was wondering if there's a way to do the same thing without having to create the .bild file. I already checked the "shape" command but you can´t use it to draw an arrow. Maybe I should make a combination between a cylinder and a rectangle to draw an arrow? Is there any other easy way to do that? Thanks in advance. Kind regards, --

3 2

bug
by Eifler, Jay Q. (UMKC-Student) 21 Feb '12

21 Feb '12

Whew, too much work to get you guys a message. The bug report form was lengthy. Anyway, can't join models in your latest 1.6 version but 1.5 whatever is fine. Should have never updated.

3 2

(no subject)
by Valia Dimitrova 21 Feb '12

21 Feb '12

Hi, I am PhD student, and I am working with program Gaussian, I want ask you some questions. I want to make maps with electrostatic potential, I saw that with this program may be make very beautiful pictures. What type of files I have to have, so I can do these pictures? I have ".fchk" files and ".cube" files. May I use them to make this? I will be very appreciated if some one can help me. thank you in advance Valia Nikolova

2 1

how to find
by Bishwajit Das 18 Feb '12

18 Feb '12

suppose i want so select ARG 188 and TRP 144 and also want to see the distance between them ....how i do this in chimera..please help me. -- BISHWAJIT DAS KOLKATA,West Bengal INDIA* * --

2 3

bug report
by Longfei Wang 16 Feb '12

16 Feb '12

Hi there, For nogui Chimera: Calls with script in nogui mode create the following traceback. terry@daisy[14]% chimera --nogui --script "coord_for_cluster_ATP_ bothChiFlip.py" Traceback (most recent call last): File "/nfs/turing2/labs/sbgrid/programs/i386-linux/chimera/1.6rc/share/__main__.py", line 69, in <module> value = chimeraInit.init(sys.argv) File "/nfs/turing2/labs/sbgrid/programs/i386-linux/chimera/1.6rc/share/chimeraInit.py", line 618, in init chimera.initializeGraphics() File "/nfs/turing2/labs/sbgrid/programs/i386-linux/chimera/1. 6rc/share/chimera/__init__.py", line 2114, in initializeGraphics bgprefs.initialize() File "/nfs/turing2/labs/sbgrid/programs/i386-linux/chimera/1.6rc/share/chimera/bgprefs.py", line 211, in initialize onDisplay=_bgPrefDisp, onHide=_bgPrefHide) File "/nfs/turing2/labs/sbgrid/programs/i386-linux/chimera/1. 6rc/share/chimera/preferences/base.py", line 567, in register self._loadCategory(category, cat, inherit, convert) File "/nfs/turing2/labs/sbgrid/programs/i386-linux/chimera/1. 6rc/share/chimera/preferences/base.py", line 482, in _loadCategory cat.load(self._loadDict[name]) File "/nfs/turing2/labs/sbgrid/programs/i386-linux/chimera/1. 6rc/share/chimera/preferences/base.py", line 249, in load opt.set(value, asSaved=True, fromPref=True) File "/nfs/turing2/labs/sbgrid/programs/i386-linux/chimera/1. 6rc/share/chimera/preferences/base.py", line 64, in set self._callback(self) File "/nfs/turing2/labs/sbgrid/programs/i386-linux/chimera/1.6rc/share/chimera/bgprefs.py", line 53, in _backgroundColorCB Midas.background(color=color, setMethod=False) File "/nfs/turing2/labs/sbgrid/programs/i386-linux/chimera/1.6rc/share/Midas/__init__.py", line 966, in background method = chimera.viewer.Solid AttributeError: '_chimera.NoGuiViewer' object has no attribute 'Solid' Best Longfei -- *Longfei Wang * Postdoctoral Fellow Harvard Medical School 250 Longwood Ave Boston, MA 02115

2 1

Re: [Chimera-users] radius of gyration
by Eric Pettersen 15 Feb '12

15 Feb '12

On Feb 14, 2012, at 6:07 PM, Daniel Gurnon wrote: > Hi Eric, > The radius of gyration interest stems from a little project I'm > starting with an undergrad. I'm interested in the radius of gyration > because I'm curious about how accurate our pictures of biomolecules > are, at least the ones that come from crystal structures. I'd like > to compare predictions of size-related behavior derived from the pdb > with measurements of size in solution, using gel filtration, SAXS, > or light scattering. That said, I don't know if I want center of > geometry or principal axis! Any advice? Hi Dan, I'm no expert, but it sounds like you want the radius of gyration from the center of geometry. Nonetheless, I'm cc'ing this to the chimera-users in case any others have thoughts. In an old posting to chimera-users, Elaine Meng provided a Fortran program to compute the radius of gyration: http://plato.cgl.ucsf.edu/pipermail/chimera-users/2009-March/003648.html I've attached a Python script that will print the center of geometry and radius of gyration from that center in Chimera's reply log. Run the script simply by opening it with File->Open. Depending on whether you feel nearby waters and ions will contribute to the measurements you are trying to model, you may or may not want to delete water and ions before running the script. --Eric Eric Pettersen UCSF Computer Graphics Lab http://www.cgl.ucsf.edu

2 2

radius of gyration
by Daniel Gurnon 14 Feb '12

14 Feb '12

Hi all, Can Chimera determine (or has a script been written to determine) the radius of gyration for a pdb structure? Dan -- ____________________________ Daniel Gurnon, Ph. D. Associate Professor of Chemistry and Biochemistry DePauw University Greencastle, IN 46135 p: 765-658-6279 e: danielgurnon(a)depauw.edu

3 2

Re: [Chimera-users] possible Chimera bug?
by Ofeliia Mada 14 Feb '12

14 Feb '12

TkkyuluR X,=)RSISXSN:-P:(:=(=(=(=(=(:O(, G Sent from my Windows Phone esdcqppwevemeem

1 0

fly command
by rswett＠chem.wayne.edu 14 Feb '12

14 Feb '12

Is anyone else having issues with the fly command lately? Seems like since RC 1.5ish flying to a savepos ignores the # of frames it's set to fly over. i.e. fly 50 in (where in is a saved position) does the fly in one frame, appearing the same as fly in or fly 1 in Additionally fly in 50 has the same effect. Has anyone else run into this? I have a few saved sessions i'd like to animate, and the older versions are having none of it. Suggestions? ~Rebecca

2 1

Re: [Chimera-users] meaning of the ANISOU PDB parameters
by Eric Pettersen 13 Feb '12

13 Feb '12

On Feb 10, 2012, at 1:55 PM, Julian Chen wrote: > hi Eric, > > First of all, hello, I was in Stroud's lab for many years, until > 2003, when I left SF. Hi Julien! > I am in a rather interesting situation where we have been dealing > with neutron and X-ray refinement, and we have been able to do > anisotropic refinement on selected hydrogen atoms, from having high > resolution neutron data. I've been using Chimera. > > I am still somewhat confused as to the meaning of the six ANISOU > terms, I understand that the first three describe the three > orthogonal axes of the ellipsoid, and that the last three somehow > relate this to the coordinate space. The six terms actually describe a 3x3 matrix that is symmetric about the diagonal. The eigenvectors and eigenvalues of that matrix are the atomic displacement axes and the mean squares of the displacements respectively. The first thing to know is that the values shown in the ANISOU records are scaled by 10**4. The other thing to know is that the six numbers correspond to these positions in the matrix: 1,1 2,2 3,3 1,2 1,3 2,3. Since the eigenvalues are mean squares, one typically works with the square roots of the eigenvalues. > Specifically, I want to look at the ANISOU for the H/D atoms, and > want to try to correlate it to its bonded neighbor (the C or the N > or the O). If possible, I'd like to extract some vector for the > 'directionality' of the aniso atoms and relate them to one another. > I want to see if the vibrational motions of the hydrogen are > correlated to the vibrational motions of its bonded neighbor. Is > there anything in CHIMERA that could do something like this, or some > appropriate mathematical formula that could output the right sort of > information? The ellipsoids option in the program gives a > qualitative description, but I'd like to get something a bit more > quantitative in nature. Depending on how comfortable you are with programming in Python, you could either write your own external code based on the info above, or write something that runs in Chimera. The latter does have the advantage of you not having to write PDB-parsing code an so on. In Chimera, atoms with corresponding ANISOU records have an 'anisoU' attribute which is a 3x3 matrix (the matrix is often referred to as 'U' in the literature, hence ANISOU record). This code snippet (basically copied from the Chimera's Aniso module) will extract the eigenvectors/values for an atom and take the root of the values: from numpy.linalg import svd ignore, lengths, axes = svd(atom.anisoU) from numpy import sqrt lengths2 = sqrt(lengths) Let me know if you need more info or help with scripting. --Eric Eric Pettersen UCSF Computer Graphics Lab http://www.cgl.ucsf.edu

2 2

Visualising AutoLigand affinity points
by George Tzotzos 09 Feb '12

09 Feb '12

Hi, I'm attaching a jpg image of ligand affinity map fill points generated by AutoLigand. My question is whether there's a way to depict these fill points as "disconnected" spheres. Actions/spheres generates something like I'd like to see but the size of the spheres is too large. Regards George

2 2

pi-cation interaction
by Bishwajit Das 09 Feb '12

09 Feb '12

how to find cation pi-interaction in a docked molecule.......plz help me -- BISHWAJIT DAS KOLKATA,West Bengal INDIA* *

2 1

How to select
by Bishwajit Das 09 Feb '12

09 Feb '12

suppose i want so select ARG 188 and TRP 144 and also want to see the distance between them ....how i do this in chimera..please help me. -- BISHWAJIT DAS KOLKATA,West Bengal INDIA* *

1 0

Re: [Chimera-users] Rotamer copying- Specific Residues only
by Eric Pettersen 08 Feb '12

08 Feb '12

On Feb 8, 2012, at 8:27 AM, Ahir Pushpanath wrote: > Eric, although the script is highly useful for what I previously > asked, I have come to the realisation that simply copying all the > apo chi angles to the holo form may not be ideal for me in this > particular situation. I noted that of around 200 residues modelled, > around 90 of them are modelled different to the apo, the rest are > modelled similar to the apo form. Of these 90, only around half are > truly changed in the holo form(literature evidence from homolog), > therefore I do require a script that is able to accept an input of > which residues to copy the chi angles for, I think this will be > useful for other users who want to copy chi angles from apo-holo > forms and other applications, and probably could be included in > later build versions of chimera, under structure editing. > > Do you think you could help me with this? I am only learning python > now, so my scripting is not at a good level for these types of things. I've attached a modified script that will only copy the chi angles for residues that are selected in the lowered-numbered model. It's up to you to select them somehow. :-) > On another note, do you have any experience with using chimera to > model apo-holo transitions any other way without using Modeller? I don't. Maybe others on the list will have suggestions... ? > Also, does chimera have the ability to do docking of ligands, i.e. > given you know the apo crystal structure, can you tell it to dock a > ligand and change conformation of protein,(similar to autodock)? Currently Chimera cannot dock ligands. If our next 5 years of funding gets approved by the NIH, we have proposed interfacing Chimera to an Autodock Vina web service that would allow docking to be carried out directly from Chimera. Right now all you can do is prepare the input for / analyze the output of external docking programs in Chimera --Eric Eric Pettersen UCSF Computer Graphics Lab http://www.cgl.ucsf.edu

1 0

how to find hydrophobic interaction
by Bishwajit Das 08 Feb '12

08 Feb '12

how to find hydrophobic interaction in a docked protein structure......please help me.... -- BISHWAJIT DAS +91-9433237920 KOLKATA,West Bengal INDIA* *

2 1

how to setup a default settings for a new pdb model?
by Xing Zhang 07 Feb '12

07 Feb '12

Hi, I have many different pdb models that need to be displayed at exactly the same settings, including color, background, stick scale and so on, so I want to set all these settings for the first model, and then apply these settings to a newly opened pdb model. Is there any way to save and load a default setting for a new pdb model? Thanks much. Xing

2 2