February 2015 - Chimera-users - RBVI Mailing Lists

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