April 2008 - Chimera-users - RBVI Mailing Lists

using EMAN/AIRS/SSEHunter
by Peter Hildebrand 10 Dec '08

10 Dec '08

Hallo there, we have installed the AIRS (EMAN 1.8) Package for Chimera (build 2470), but have difficulties to get it run. steps to reproduce: launch Chimera load .mrc file calling tools/AIRS/SSEHunter there is no file to select in the 'MRC file' selection box. If we provide parameters for Angstrom, Resolution and threshold and press 'Apply' we get this: ================================================================= UnboundLocalError Exception in Tk callback Function: <function command at 0x89f67d4> (type: <type 'function'>) Args: () Traceback (innermost last): File "CHIMERA/lib/python2.4/site-packages/Pmw/Pmw_1_2/lib/PmwBase.py", line 1747, in __call__ None File "CHIMERA/share/chimera/baseDialog.py", line 238, in command None File "/usr/local/EMAN/chimeraext/Airs/ssehunter.py", line 89, in Apply targetPath=airsIO.getPath(target) File "/usr/local/EMAN/chimeraext/Airs/airsIO.py", line 62, in getPath return path UnboundLocalError: local variable 'path' referenced before assignment ================================================================= Are we using this module in the wrong way (wrong parameters e.g.) or is it a bug (where?). Thanks in advance, jan and peter

3 3

Drawing a circle or square
by Elaine Meng 12 May '08

12 May '08

Hi Pradeep, Other people have recommended Adobe Illustrator and Inkscape for such purposes: http://www.cgl.ucsf.edu/pipermail/chimera-users/2008-April/002579.html However, you CAN draw a rectangle in the BILD format and read that into Chimera. For example, put the following 7 lines in a text file named square.bild and open it with File... Open: .scale 1.5 .color red .move 0 0 0 .draw 5 0 0 .draw 5 5 0 .draw 0 5 0 .draw 0 0 0 I also attach said file below, for your convenience. You can edit the coordinates, scale, and color. Once you open it in Chimera, it is a separate model, and you can move it independently using the checkboxes below the Command Line or the "Active" checkboxes in the Model Panel. Advantage: since it does exist in 3D, you can place it between the surface and the nucleic acid. Disadvantage: cannot make curves/circles, may be hard to position/orient, cannot control line thickness. The thickness issue may be addressed by instead making the edges cylinders and the corners balls. This can also be done in the simple BILD format, but it is sort of a hack to make a 3D object when you really want a 2D object. BILD format description: http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/bild.html I hope this helps, Elaine ----- Elaine C. Meng, Ph.D. meng(a)cgl.ucsf.edu UCSF Computer Graphics Lab and Babbitt Lab Department of Pharmaceutical Chemistry University of California, San Francisco http://www.cgl.ucsf.edu/home/meng/index.html On Apr 30, 2008, at 11:57 AM, Pallan, Pradeep S wrote: > Hi Elaine, > I am wondering whether one can draw a circle or a square > to indicate a site/s of a protein? Attached is a tiff file that has a > square and a circle. I want the regions to be pointed out. > > So what I did: > 1) Generated the image in chimera > 2) added a circle and a square in microsoft paint (Doesn't look > good!). > > I was going through the chimera tutorials, but didn't find the > feature. > Any suggestions for adding a circle or square? In my case, the > surface image is at the back, and nucleic acid in the front. > I would like to draw the circle and square and then load the nucleic > acids, so that will look better. In other words the square or > circle in > between > the surface image (at the back) and the nucleic acid (in front). > > If editing the tiff file using Adobe photoshop (instead of MS paint) > it does a better job, thick lines etc, but it leaves patches on the > nucleic > acid. > I appreciate your help. > Thanks, > Pradeep

2 1

pretty nice article on PDB website about Chimera
by Eric Pettersen 30 Apr '08

30 Apr '08

http://www.rcsb.org/pdb/general_information/news_publications/ newsletters/2008q1/education_corner.html

2 1

Initial window size
by Bernard Heymann 30 Apr '08

30 Apr '08

Where is the initial graphics window size set? I did not see any preference for that. I would like to change that to have Chimera open with a bigger window. Thanks Bernard Heymann Rm 1515, 50 South Dr, NIH, Bethesda, 20892-8025 Tel 301-451-8241, Fax 301-480-7629

2 1

Can I re-save the mol2 file of compounds sorted by grid score?
by �� 29 Apr '08

29 Apr '08

Hi there, I'm a novice in Chimera. I have done the ViewDock tutorial, but now I have a question. I know, in the ViewDock, the listing can be sorted by any column(e.g. Grid Score). And I expect to know if I could re-save the mol2 file of the compounds sorted by grid score. Any answer would be appreciated! Best regards, Yolanda Guo Northeast Normal University --------------------------------- Be a better friend, newshound, and know-it-all with Yahoo! Mobile. Try it now.

2 1

Calling Chimera from another program
by Tom Goddard 29 Apr '08

29 Apr '08

Hi Brian, I talked to Niels Volkmann on Friday about how PyCoan can use Chimera for displaying maps. There are two existing ways for programs to communicate with Chimera. Here are some info about those. Since PyCoan uses Python it would be possible for it to run under the same Python interpreter as Chimera. Current Chimera daily builds use a standard Python 2.5 distribution with a bunch of additional modules (numpy, pil, pytables, ...) included. Although this approach gives the most power for passing data between PyCoan and Chimera, it requires the two programs to use exactly the same versions of third party Python modules and that can cause problems. Also it allows unlimited access to Chimera internals and use of that functionality would likely make PyCoan/Chimera communication break as Chimera code evolves. There are two ways to communicate with Chimera from a separate process. One involves sending Chimera commands (e.g. "open mymap.mrc") to Chimera's stdin. You have to start Chimera with the following option to get it to read stdin "chimera --start ReadStdin". The other method involves sending Python code to Chimera using a socket. The second approach works by running "chimera --send somecode.py" which will connect to an existing Chimera on the local machine or start one if none is running. Both methods require Chimera and the controlling process to run on the same machine. The second method does not require that limitation but we impose it because there is no authentication method -- so if you allow remote connections, anyone could run arbitrary Python code on your machine (unlikely to be abused, but we're not going to distribute that capability). The first approach is limited to Chimera commands and many Chimera capabilities available in the graphical user interface have no equivalent command. All the capabilities are available if Python code is used. Sending Chimera commands returns the text output of the command although this rarely contains useful information. Sending Python code with "chimera --send xxx.py" returns nothing. So the communication is essentially uni-directional. The advantage of these communication methods over just start Chimera with a Python script as an argument is that you can maintain a connection to a single instance of Chimera and update the visualization as needed. Both the stdin and socket methods need map data to be written to a separate file if you want to open that. Niels was interested in avoiding the use of temporary files and suggested using shared memory. Python has some shared memory support but I recommend against that. The performance is unlikely to be dramatically better. If you write a temporary file the operating system is going to cache it in memory rather than writing it to disk and Chimera will then read it from memory. If you delete the temporary file soon it may never be written to disk. You could force this behavior using shared memory but that is just a more complex way to achieve the same optimization that the operating system is already providing. In any case I would not pursue that until you test the performance of using files. I've attached an example Python script that starts Chimera and sends commands to open a volume and set its contour level and color. I tested it with a Chimera daily build (1.2502). I recommend using the daily builds instead of the November 2007 production release since many volume data improvements have been made. Here is Chimera documentation on the ReadStdin and "--send" communication methods. http://www.cgl.ucsf.edu/chimera/docs/ContributedSoftware/readstdin/readstdi… http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/options.html Tom # ----------------------------------------------------------------------------- # Example program that sends commands to Chimera. # class Chimera_IO: def __init__(self): args = ['chimera', '--start', 'ReadStdin'] self.prog = run_program(args) def send_command(self, command): p = self.prog p.stdin.write(command + '\n') p.stdin.flush() # Collect output. lines = [] while True: line = p.stdout.readline() if line == '' or line == 'END\n': break lines.append(line) output = ''.join(lines) return output # ----------------------------------------------------------------------------- # def run_program(args): # Hide console window on Windows. import sys if sys.platform == 'win32': from subprocess import STARTUPINFO, STARTF_USESHOWWINDOW spinfo = STARTUPINFO() spinfo.dwFlags |= STARTF_USESHOWWINDOW else: spinfo = None from subprocess import Popen, PIPE try: p = Popen(args, stdin = PIPE, stdout = PIPE, startupinfo = spinfo) except: print 'Starting "%s" failed.\n' % ' '.join(args) return None return p # ----------------------------------------------------------------------------- # c = Chimera_IO() print c.send_command('open 0 /tmp/sfv.mrc') print c.send_command('volume #0 level 2500 color tan')

1 0

crystal contacts to show molecules
by Miguel Ortiz-Lombardía 28 Apr '08

28 Apr '08

Dear all, Would it be possible to modify the code of the Crystal Contact tool so it can produce, as an option, the molecules that are in contact with the model as full atom models instead of spheres? I would be happy to do so, but then I would appreciate if you can point me to the relevant code in chimera to do so. Thank you! Miguel -- http://www.pangea.org/mol/spip.php?rubrique2 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Je suis de la mauvaise herbe, Braves gens, braves gens, Je pousse en liberté Dans les jardins mal fréquentés! Georges Brassens

2 1

post-processing of Chimera images
by Jonathan Hilmer 28 Apr '08

28 Apr '08

Francesco: Especially since you're preparing images for publication, I would recommend you use software such as Inkscape (a excellent free alternative to Adobe Illustrator) to handle the composition of figures. It's a vector-based program that's much better suited to handling geometric shapes as compared to a raster program like Photoshop or Gimp. The final post in the following discussion has a good overview with screenshots showing how Inkscape can be used: http://inkscape-forum.andreas-s.net/topic/146897 Here are some of the screenshots: http://jo.irisson.free.fr/dropbox/inkscape/poster_inkscape_style2.png http://jo.irisson.free.fr/dropbox/inkscape/poster_inkscape_outline.png I've got an example figure that you can see that uses some of the basic features with images, shapes, figures, and text: ftp://guest:@bothnerlab01.chemistry.montana.edu/inkscape_example.zip The best part of using Inkscape is that it's totally resolution-independent, so you can export graphics for publication at whatever resolution or format is required. There are only two disadvantages to Inkscape that I have found so far. First, the PDF export option works fine but treats text as shapes, so the resulting files are huge if there is substantial amounts of text: ~500MB for a poster with about 50% text. Second, it doesn't embed images within the saved files the way Powerpoint does: it just saves links in the manner of OpenOffice. I rather prefer that behavior, but you have to be prepared for it if you're moving files around. Regardless, once you give it a try you'll never go back to using Powerpoint or Photoshop for making figures or posters. Jonathan On Mon, Apr 28, 2008 at 1:00 PM, <chimera-users-request(a)cgl.ucsf.edu> wrote: > Send Chimera-users mailing list submissions to > chimera-users(a)cgl.ucsf.edu > > To subscribe or unsubscribe via the World Wide Web, visit > http://www.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. Re: Numbering the amino acid sequence (Jeff Speir) > 2. Re: Chimera - Matchmaker tool (Elaine Meng) > 3. Re: Numbering the amino acid sequence (Eric Pettersen) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Mon, 28 Apr 2008 09:07:02 -0700 > From: Jeff Speir <speir(a)scripps.edu> > Subject: Re: [Chimera-users] Numbering the amino acid sequence > To: chiendarret(a)yahoo.com > Cc: chimera <chimera-users(a)cgl.ucsf.edu> > Message-ID: <B2AD9487-6EEF-4AC4-8060-51788F8E86B2(a)scripps.edu> > Content-Type: text/plain; charset=US-ASCII; delsp=yes; format=flowed > > Outside of Chimera, there are several ways to do this in image > editing programs (Gimp, Photoshop, etc.). Use layers and: 1) draw > arrows; or 2) place the number on top of a box just big enough to fit > the text (opaque to semi-transparent); or 3) use effects on the text, > such as adding an outer glow effect in Photoshop, which adds a halo > around the text so that its outline stands out from the image, making > it readable. > > I would advise not going crazy with colors, this just adds to the > complexity of the image. > > Cheers, > > Jeff > > > > On Apr 27, 2008, at 10:49 AM, Francesco Pietra wrote: > > > A referee is asking me to number the amino acid sequence involved > > in docking a ligand with a protein. Not satisfied by having a side > > table with a guideline. > > > > Is it possible to add that numbering in a least-intruding way? The > > figure exported in standard graphic formats from Chimera is already > > extremely crowded. Therefore, is any possibility of adding the > > numbering outside the protein with some kind of pointers? I can't > > do with colors alone: I already used different colors for different > > helices. > > > > Thanks > > francesco pietra

1 0

Numbering the amino acid sequence
by Francesco Pietra 28 Apr '08

28 Apr '08

A referee is asking me to number the amino acid sequence involved in docking a ligand with a protein. Not satisfied by having a side table with a guideline. Is it possible to add that numbering in a least-intruding way? The figure exported in standard graphic formats from Chimera is already extremely crowded. Therefore, is any possibility of adding the numbering outside the protein with some kind of pointers? I can't do with colors alone: I already used different colors for different helices. Thanks francesco pietra ____________________________________________________________________________________ Be a better friend, newshound, and know-it-all with Yahoo! Mobile. Try it now. http://mobile.yahoo.com/;_ylt=Ahu06i62sR8HDtDypao8Wcj9tAcJ

3 2

Re: [Chimera-users] Chimera - Matchmaker tool
by Elaine Meng 28 Apr '08

28 Apr '08

On Apr 27, 2008, at 10:16 AM, Sujata wrote: > Hi, > I am using Chimera for awhile now to analysis modeled structures > and I love it, way better for visualizing than a lot of other > softwares and very user freindly. > I do have a question about Matchmaker- is there a way that I can > get the total number of atom-pairs in the 'matchlist' it generates > in Matchmaker? Or in other words, I want to know the final > percentage of atom-pairs that Chimera gives me after pruning the > long atom pairs that do not exceed 2 Angstroms to compute the RMSD. > > Thank you very much. > Sujata Sovani > Hi Sujata, Thanks, tell your friends! 8-) The final match RMSD and number of alpha-carbon pairs are reported in the Reply Log (under Favorites). They are also shown in the status line at the bottom of the Chimera window, but only temporarily. If you also show the sequence alignment from MatchMaker, these residue pairs used in the final fit iteration are shown in boxes with orange- red outline and light orange interior. I am not sure which percentage you wanted; one possibility is to divide by the total number of residues in a structure. You can also find out how many pairs were pruned by doing another match with the same parameters except without iteration. That does not change the sequence alignment, only the fitting. If you already showed the sequence alignment from matchmaker, you wouldn't have to rerun matchmaker, only choose Structure... Match from the sequence alignment window and apply it with all its options turned off. I hope this helps, Elaine ----- Elaine C. Meng, Ph.D. meng(a)cgl.ucsf.edu UCSF Computer Graphics Lab and Babbitt Lab Department of Pharmaceutical Chemistry University of California, San Francisco http://www.cgl.ucsf.edu/home/meng/index.html

1 0

Chimera-users April 2008