- Chimera-users - RBVI Mailing Lists

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