- Chimera-users - RBVI Mailing Lists

building molecules
by Ben Schuler 31 Aug '04

31 Aug '04

Dear Chimera developers, It seems to me that the manual formation of bonds and building of any type of molecules (not only proteins and nucleic acids), which was a strength of MIDAS, is not supported by Chimera. I would very much welcome an implementation of these features to make Chimera a real replacement for MIDAS. Ben __________________________________________________ Prof. Ben Schuler, PhD Biochemisches Institut, Universität Zürich Winterthurerstr. 190, CH-8057 Zürich Tel: +41 44 635 5535, Fax: +41 44 635 5907 Web: www.biochem.unizh.ch/schuler Email: schuler(a)bioc.unizh.ch

1 0

chimera --send with multiple Chimeras
by Thomas Goddard 22 Jul '04

22 Jul '04

Hi Sandeep, If you have more than one Chimera running on a machine then "chimera --send" sends the file to the most recent one that is accepting connections. Each time you change the Favorites/Preferences/Web Data option to accept connections, that Chimera becomes the most recent. The way this works is there is a list of running Chimera instances that are accepting connections in file (on unix /tmp/chim_webinfo-<username> where <username> is your login name). This file lists the Chimera process id and port number it is listening on. The last line indicates the most recent Chimera that is accepting connections. I wouldn't recommend relying on info in this file since the implementation may change in the future. Tom

3 3

suggest OBJ export option
by Steven McQuinn 19 Jul '04

19 Jul '04

I'm really happy with Chimera, but vexed that every attempt to export a high resolution mesh via the VRML utility ends in a python failure. Have not tried exporting molecular models, only isosurfaces with 250,000+ polygons. What I use instead for export is VMD, which has an .stl export option, providing smaller file sizes than VRML. VMD is buggy on my XP system, requiring several attempts to export correctly. I use Chimera to get my isosurface settings right, then plug these values into VMD. SUGGESTION: If you provided an OBJ mesh export option, nearly every 3d application could utilize the product. The Wavefront OBJ format is the defacto universal exchange format, fairly small file sizes, retaining textures. Should be easy for you geniuses to program. Also, the stereo camera option never "takes." The walleye and crosseye cameras work OK but not the straight stereo pair, which I have to compose from the right and left images. Would like to send the development team some work samples. How do I do that? Please reply to my email address, as I don't read the list. Steven McQuinn _________________________________________________________________ MSN Toolbar provides one-click access to Hotmail from any Web page FREE download! http://toolbar.msn.click-url.com/go/onm00200413ave/direct/01/

2 1

Re: [Chimera-users] PDB-file trajectories
by Eric Pettersen 16 Jul '04

16 Jul '04

I have PDB-file trajectory support working in Movie now. It's only in our internal development version, but I could make a custom snapshot release for anyone who wants to use it immediately. Otherwise, it will be in the next regular release, which will probably be in a month or two. Eric Pettersen UCSF Computer Graphics Lab pett(a)cgl.ucsf.edu http://www.cgl.ucsf.edu > On Thursday, January 8, 2004, at 06:52 AM, amutha ramaswamy wrote: > >> It is possible to make movies with PDB structurers in chimera. this >> would really be helpfull in visualing the conformatinal changes. I >> know chimera has the option for movie with trajectory files. but i >> wish to kow the option with PDB files only. > > Hi Amutha, > I assume you're talking about the situation where you have one PDB > file for each "frame" of a trajectory (or a multi-MODEL PDB file of > the trajectory), rather than a starting and ending PDB file that you > want to "morph" between. Either way, unfortunately, we don't currently > support what you want. However, now that you have pointed out the > usefulness for plain-PDB support in Movie, I can probably add it in > the not-too-distant future. But given my other priorities (including > support for attribute coloring that you mentioned earlier!) the > "not-too-distant future" might still be a couple of months. I'll send > mail when there is at least a preliminary version available.

1 0

How to send files to Chimera
by Thomas Goddard 02 Jul '04

02 Jul '04

Hi Sandeep, To open files in an already running Chimera from your program you can use the chimera --send option. For example, chimera --send myfile.mrc You need to tell Chimera to accept "web data" for this to work. You do this using Chimera menu entry Favorites/Preferences, category Web Data, set "Accept web data" to true. Press the Save button if you want this setting to apply to all future Chimera sessions (saved in ~/.chimera/preferences). If you don't set this preference, the --send command-line option will just start another copy of Chimera. This is described in the Chimera manual where command-line options are described: http://www.cgl.ucsf.edu/chimera/1.1951/docs/UsersGuide/options.html You can do fancier things, like send Python code (*.py files) or Chimera commands (*.cmd or *.com files) to Chimera to show your data in a certain way. You would still use the --send option to send the Python file. These remote control features were added to Chimera so that a web browser link can cause commands and files to be sent to Chimera and this is described in the Chimera manual: http://www.cgl.ucsf.edu/chimera/1.1951/docs/ContributedSoftware/webdata/chi… The hard part here is that volume data display options like threshold levels do not have simple Chimera commands to set them. Instead you have to dig around in the VolumeViewer Python code (in the share directory) in order to see what Python functions set these things. This is not easy to figure out. I can help you if you have specific questions about the Python functions for volume data display. Tom > Date: Fri, 2 Jul 2004 10:01:40 -0600 > From: Sandeep N Seshadri <sandeep(a)nmsu.edu> > To: Thomas Goddard <goddard(a)cgl.ucsf.edu> > Subject: Re: [Chimera-users] Re: [chimera-dev] regd error when opening a .mrc file > > Hi Tom, > sorry for being very confusing . i am trying to open up files generated from my > program . as of now all i do is chimera test1.mrc > > and next time i generate another file i would open it as chimera test2.mrc > > so it would open a new chimera window and show the next .mrc file . what i would > like to do is when i open the second file test2.mrc i would like to see it open > in the already opened chimera window for test1.mrc. > > do you think i can do this ? Please let me know > > Sandeep Seshadri

1 0

[Chimera-users] Re: [chimera-dev] regd error when opening a .mrc file
by Thomas Goddard 01 Jul '04

01 Jul '04

Hi Sandeep, This error indicates that the mrc file data values are all zero. The error occured when making the histogram that is shown in the Volume Viewer dialog, when it was trying to scale it to fit the width of the dialog. Of course, our code should not generate an error in this case. I will fix it for the next release. If you will be opening mrc files that may be all zero you can change chimera/share/VolumeViewer/histogram.py line 26 hscale = float(cheight) / max_height to if max_height == 0: hscale = 1 else: hscale = float(cheight) / max_height Then you will not get any error and the histogram will display correctly. Thanks for pointing out this problem. I am not familiar with pdb2mrc so I do not know why it is generating a volume consisting of all zeros. Tom > Date: Wed, 30 Jun 2004 21:33:12 -0600 > From: Sandeep N Seshadri <sandeep(a)nmsu.edu> > To: chimera-dev(a)cgl.ucsf.edu > Subject: [chimera-dev] regd error when opening a .mrc file > > i am getting the following error when i try to open a .mrc file . i generated a > .pdb file and use a pdb2mrc to generate a .mrc file . i tried opening the .pdb > file it works fine but when i try to open the .mrc file it gives the error. > > if anyone has any idea on why this error occurs please let me know > > thanks > sandeep > > > Error while processing test-helix.mrc: > ZeroDivisionError: float division > (see reply log for Python traceback info) > Traceback (most recent call last): > File "/home/sandeep/chimera/share/__main__.py", line 56, in ? > value = chimeraInit.init(sys.argv) > File "/home/sandeep/chimera/share/chimeraInit.py", line 275, in init > prefixableType=1) > File "/home/sandeep/chimera/share/chimera/__init__.py", line 886, in open > models = func(filename, *args, **kw) > File "/home/sandeep/chimera/share/VolumeViewer/ChimeraExtension.py", line 30, > in open_cb > d.add_data_sets([g]) > File "/home/sandeep/chimera/share/VolumeViewer/__init__.py", line 175, in > add_data_sets > self.redisplay_cb() > File "/home/sandeep/chimera/share/VolumeViewer/__init__.py", line 459, in > redisplay_cb > self.show_using_dialog_settings(self.focus_region) > File "/home/sandeep/chimera/share/VolumeViewer/__init__.py", line 380, in > show_using_dialog_settings > tp.update_histograms(data_region) > File "/home/sandeep/chimera/share/VolumeViewer/__init__.py", line 1142, in > update_histograms > cw.update_histogram(data_region, read_matrix, self.dialog.message) > File "/home/sandeep/chimera/share/VolumeViewer/__init__.py", line 1465, in > update_histogram > h.show_data(Numeric.log(s.counts + 1)) > File "/home/sandeep/chimera/share/VolumeViewer/histogram.py", line 26, in > show_data > hscale = float(cheight) / max_height > ZeroDivisionError: float division

2 2

Is it possible to undo an action?
by Lesheng Kong 01 Jul '04

01 Jul '04

Hi, Is it possible to undo an action in Chimera? I notice that some actions in Chimera can be undid, such as "undo selection" or "undo move". Can every action be undid in Chimera? For example, I colored a selection zone in a different color. But I want to undo this action and recover each atom's previous color. Thank you very much! With best regards, Lesheng

2 2

Re: Chimera questions.
by Thomas Goddard 30 Jun '04

30 Jun '04

Hi Matt, Support for DCD trajectories is in our current Chimera development version, but is not in the most recent release, Chimera 1.1951. It is described in the development version of the Chimera manual: http://www.cgl.ucsf.edu/chimera/docs/ContributedSoftware/movie/movie.html If you like, I can give you the new Python code for DCD trajectories. I have not yet tested whether it will work in Chimera 1951, but I can test it. Otherwise, you can get it in our next release which is at least a few months away. To reposition a PDB model using Python code see the Programmer FAQ: http://www.cgl.ucsf.edu/chimera/docs/ProgrammersGuide/faq.html#q4 Tom > Date: Wed, 30 Jun 2004 16:53:00 +0000 > From: Matthew Baker <mlbaker(a)gmail.com> > To: goddard(a)cgl.ucsf.edu > Subject: Chimera questions. > > Tom, > > Two quick questions: > 1. Is there a reader/display for DCD or netCDF trajectory files in > Chimera? I see that you can read AMBER trajectory files in the > Molecular Dynamics module. > > 2. How do you update a PDB file loaded into Chimera? I have a PDB file > and I want to apply a transformation to it and have the updated > coordinates applied in the viewer. > > Thanks, > > Matt

1 0

Filled in/triangulated nucleotide base ring planes
by Pattanayek, Sabuj 28 Jun '04

28 Jun '04

Hi, I was wondering if I missed this feature in Chimera somewhere...a way to fill in nucleotide base rings/planes, this is what I mean: http://linegli.mc.vanderbilt.edu/~sabuj/1GQU_ribbons.png (using ribbons v3.2) http://linegli.mc.vanderbilt.edu/~sabuj/1GQU_vmd.png (using vmd 1.8.2) If you've ever used Ribbons you know it's a real pain to setup. I actually ended up using Chimera's selection methods to create separate pdb files for inputs into Ribbon's scripts, and it doesn't do multi-colored surfaces/transparency very well. VMD as you can see is filling in the ribose sugars (which I don't want) and is failing to fill in some of the bases. I tried recompiling it's code with some modifications but failed due to some header problems. Anyways I was wondering if this feature might be added, if it already isn't there, with the ability to choose colors for the base rings, or attribute certain colors to groups of bases, etc. Thanks, Sabuj

2 1

Nice postscript output from chimera.
by J B Procter 26 Jun '04

26 Jun '04

Does anyone know of a way of generating postscript from chimera via a vector rendering of the structure, rather than simply embedding the rendered image ? What would be even more fun is if one could generate a molscript diagram from a chimera view. j.

2 2

Re: Save options
by Eric Pettersen 24 Jun '04

24 Jun '04

On Jun 24, 2004, at 2:30 PM, Karl N. Kirschner wrote: > Hi, > We are looking at chimera to use in our department here at Hamilton. > I am currently evaluating it and I was wondering is there a way to save > a structure under a different format other than filename.py (ie PDB)? > From your website it appears that one should be able to do this. We > are > using the chimera-1.1951-linux.exe and chimera-1.1917-linux_x86_64.exe > versions of the program. Maybe it is a bug - or maybe I am not just > seeing something. > > Thanks, > Karl Hi Karl, You can write out structures in PDB format either using the "write PDB" button of the model panel (Favorites..Model Panel in the menus), or by using the "write" command in the command line ("help write" will provide details of using the write command). Both of these methods save the entirety of the models specified, displayed or not. To save parts of models, one can either delete the unwanted portions (Actions..Atoms/Bonds..delete in the menus or using the "delete" command), or (on UNIX-based systems) use the pdbrun command to output only the displayed portions ("help pdbrun" for details). Please feel free to ask about any other things that may not be obvious -- it helps us to refine the usability of Chimera's interface. I should note that usage-style questions like these should probably go to the chimera-users mailing list rather than chimera-bugs though. Thanks for giving Chimera a try! Eric Pettersen UCSF Computer Graphics Lab

1 0

[Chimera-users] Chimera 1.1951 on OSX 10.3.4 fails to do much at all...
by Grant Jacobs 20 Jun '04

20 Jun '04

re: Chimera 1.1951 on OSX 10.3.4 fails to do much at all... I'm a first time user obviously! :-) Trying to open Side View or any of the Tools:Viewing Preferences fails. Many commands via the command interface likewise fail, e.g. surface #0 (no module msms), (no module webrowser). AFAICS from the error messages, the errors appear to be either Python modules missing or Tcl/TK errors. My download is the correct size. Q: Does Chimera require a Tcl/TK installation? - which doesn't show these up and they aren't hidden in one of the sbin's (I'm assuming that they are named tcl and tk - I've never had to work with them directly myself). If so is this tested on start up? There appears to be nothing about this in the docs AFAICS. The OSX installation instructions only say to drag the app. to the disk, nothing about other requirements. So naturally that's all I've done so far! If it does require more stuff, would someone be so kind as to tell me exactly what and where to get it? Two suggestions: It'd be good to have a simple test of what modules are present to "verify" an installation. (e.g. Tools:test:modules) It'd be good to have a test of the presence of Tcl/TK & let the user know it its not around. Grant E.g. of traceback: ImportError: No module named Lighting.controller (see reply log for Python traceback info) Traceback (most recent call last): File "/Applications/Science/Chimera-1.1951.app/Contents/Resources/share/__main__.py", line 56, in ? File "/private/automount/usr/local/src/staff/chimera-build/Darwin-X11/install/Chimera.app/Contents/Resources/share/chimeraInit.py", line 301, in init File "tkgui.py", line 2891, in eventLoop File "/usr/local/src/staff/chimera-build/Darwin-X11/foreign/Python-2.3.2/lib/python2.3/lib-tk/Tkinter.py", line 965, in mainloop File "/private/automount/usr/local/src/staff/chimera-build/Darwin-X11/install/Chimera.app/Contents/Resources/lib/python2.3/site-packages/Pmw/Pmw_1_2/lib/PmwBase.py", line 1751, in __call__ File "/private/automount/usr/local/src/staff/chimera-build/Darwin-X11/install/Chimera.app/Contents/Resources/lib/python2.3/site-packages/Pmw/Pmw_1_2/lib/PmwBase.py", line 1814, in _reporterror File "/private/automount/usr/local/src/staff/chimera-build/Darwin-X11/install/Chimera.app/Contents/Resources/lib/python2.3/site-packages/Pmw/Pmw_1_2/lib/PmwBase.py", line 1377, in displayerror File "/private/automount/usr/local/src/staff/chimera-build/Darwin-X11/install/Chimera.app/Contents/Resources/lib/python2.3/site-packages/Pmw/Pmw_1_2/lib/PmwBase.py", line 1747, in __call__ File "base.py", line 434, in activate File "base.py", line 80, in menuActivate File "ChimeraExtension.py", line 20, in activate File "base.py", line 102, in module ImportError: No module named Lighting.controller AttributeError: ViewerDialog instance has no attribute 'nb' (see reply log for Python traceback info) Traceback (most recent call last): File "/Applications/Science/Chimera-1.1951.app/Contents/Resources/share/__main__.py", line 56, in ? File "/private/automount/usr/local/src/staff/chimera-build/Darwin-X11/install/Chimera.app/Contents/Resources/share/chimeraInit.py", line 301, in init File "tkgui.py", line 2891, in eventLoop File "/usr/local/src/staff/chimera-build/Darwin-X11/foreign/Python-2.3.2/lib/python2.3/lib-tk/Tkinter.py", line 965, in mainloop File "/private/automount/usr/local/src/staff/chimera-build/Darwin-X11/install/Chimera.app/Contents/Resources/lib/python2.3/site-packages/Pmw/Pmw_1_2/lib/PmwBase.py", line 1751, in __call__ File "/private/automount/usr/local/src/staff/chimera-build/Darwin-X11/install/Chimera.app/Contents/Resources/lib/python2.3/site-packages/Pmw/Pmw_1_2/lib/PmwBase.py", line 1814, in _reporterror File "/private/automount/usr/local/src/staff/chimera-build/Darwin-X11/install/Chimera.app/Contents/Resources/lib/python2.3/site-packages/Pmw/Pmw_1_2/lib/PmwBase.py", line 1377, in displayerror File "/private/automount/usr/local/src/staff/chimera-build/Darwin-X11/install/Chimera.app/Contents/Resources/lib/python2.3/site-packages/Pmw/Pmw_1_2/lib/PmwBase.py", line 1747, in __call__ File "base.py", line 495, in activate File "base.py", line 80, in menuActivate File "StdTools.py", line 81, in activate File "StdTools.py", line 206, in raiseViewingTab AttributeError: ViewerDialog instance has no attribute 'nb' -- ------------------------------------------------------------------- Grant Jacobs Ph.D. BioinfoTools ph. +64 3 478 0095 (office, after 10am) PO Box 6129, or +64 25 601 5917 (mobile) Dunedin, gjacobs(a)bioinfotools.com NEW ZEALAND. Bioinformatics tools: deriving knowledge from biological data Bioinformatics tools - software development - consulting - training Check out the website for more details: http://www.bioinfotools.com The information contained in this mail message is confidential and may be legally privileged. Readers of this message who are not the intended recipient are hereby notified that any use, dissemination, distribution or reproduction of this message is prohibited. If you have received this message in error please notify the sender immed- iately and destroy the original message. This applies also to any attached documents.

4 5

Re: chimera questions
by Thomas Goddard 11 Jun '04

11 Jun '04

Hi Wen, There is not a Chimera command to print the transformation of one model relative to another. But you can get it with some Python code. Below is some Python code to add two keyboard shortcuts that print the transform of model 0 relative to model 1 or model 1 relative to model 0. The model numbers are shown in the Model Panel dialog (Tools/Inspectors or Favorites menu), 0 being the first model you opened. To use the keyboard shortcuts copy the Python code below to a file, say accel.py. Show the accelerator dialog, Tools/Keyboard/Accelerator List, and enter the path to the accelerator file and press Load. Also click the "Accelerators on" switch. After loading two models you can type x0 or x1 in the Chimera graphics window and it will print out the transform to the Chimera reply log (Favorites/Reply Log). The accelerator file is remembered so you only need to do this setup once. To have the accelerators be on whenever you start Chimera, set Accelerators On to "auto start" in Favorites/Preferences/Category Tools. The transformation matrix is a 3 by 4 matrix. The first 3 columns give the rotation, and the 4th column is a translation. The rotation is applied first, then the translation. Tom def print_transform(model_id_1, model_id_2): '''Print transformation matrix of model 1 relative to model 2.''' import chimera m1 = chimera.openModels.list(id = model_id_1)[0] # returns a list m2 = chimera.openModels.list(id = model_id_2)[0] xf1 = m1.openState.xform xf2 = m2.openState.xform xf2.invert() xf1.multiply(xf2) print xf1 def register_accelerators(): from Accelerators import standard_accelerators, add_accelerator standard_accelerators.register_accelerators() add_accelerator('x0', 'Print model 0 transform relative to model 1', lambda: print_transform(0,1)) add_accelerator('x1', 'Print model 1 transform relative to model 0', lambda: print_transform(1,0))

1 0

Re: chimera questions
by Thomas Goddard 10 Jun '04

10 Jun '04

Hi Wen You can write out coordinates of a PDB model relative to another model such volume you hand docked to, using Model Panel (under Favorites menu or Tools/Inspectors menu). Show the Model Panel dialog, click on the PDB model in the list, press the "Write PDB" button in the right column of the Model Panel. Then a Save file dialog comes up and at the bottom is a checkbutton that says "Save relative to model chosen in model panel when 'OK' is hit". You then click on the volume listed in the Model Panel that the PDB coordinates should be relative to. Then press OK in the Save dialog. A bit crazy. We'll probably improve this confusing interface in a future release. Tom > From: Wen Jiang <wjiang(a)bcm.tmc.edu> > Subject: Re: chimera questions > Date: Thu, 10 Jun 2004 17:06:28 -0500 > To: Thomas Goddard <goddard(a)cgl.ucsf.edu> > > Tom, > Is it possible to do manual fit of PDB to a density map in Chimera then > save the transformed PDB coordinates to a new PDB file? If it is not > possible, how to read out the transformations applied to the PDB model? > Thanks, > > Wen

1 0

Re: save session/grasp surfaces
by Thomas Goddard 07 Jun '04

07 Jun '04

Hi Sabuj, We'll consider adding DelphiViewer session saving. It is a fair bit of work and we have been thinking of replacing DelphiViewer's surface coloring capabilities with something that works with volume data other than electrostatic potential maps. I'll discuss it with other Chimera developer's and get back to you. Tom > Date: Mon, 7 Jun 2004 11:49:00 -0500 > From: "Pattanayek, Sabuj" <sabuj.pattanayek(a)Vanderbilt.Edu> > To: goddard(a)cgl.ucsf.edu > Subject: save session/grasp surfaces > > Hi, > > Method 1 would be what I am looking for (reload potential file, grasp surface > file, recalculates colors according to saved preferences, and places the surface > where it was before saving). I accidentally deleted your message, but I think > that was the gist of choice #1. > > Thanks, > Sabuj Pattanayek

1 0

[Chimera-users] v1.1951 save session not reopening grasp surfaces that were there before i saved session
by Pattanayek, Sabuj 04 Jun '04

04 Jun '04

Hi, I had to resort to using grasp surfaces because MSMS doesn't behave "nicely" when dealing with multiple free floating residues that come directly from a PDB file. Although grasp surfaces are lower in visual quality than MSMS, grasp rendered all the areas I wanted whereas the MSMS plugin could not. So, after I did save session and re-opened the session the grasp surfaces did not reappear. Is there anyway to add functionality so that the grasp surfaces, colors, and the corresponding locations can be recalled when I load the saved session file? This saves time especially when there are many grasp surface components each with different colors. One more minor detail/UI glitch I noticed?: Color Editor -> click opacity : the "A" row along with slider appears -> unclick opacity : the sliders disappear but the "A" along with row space remains ...don't know if that was done on purpose. Thanks for the great program!

2 1

Re: [Chimera-users] non-GUI use of Chimera
by Charles Moad 25 May '04

25 May '04

Hi, I am working with Randy on this, and I think we have a web service type solution that will work using your suggestions. The one problem that remains is killing chimera after a request. The webdata facility is targeted to run on a client's machine where a user quits chimera, but we will be running it on the server. I tried adding the 'stop' command at the end of the chimerax file, but I received an error from chimera and it continued to run. Same story for 'sys.exit()' Running the chimera_send script a second time only opens a new instance of chimera which also remains open. Is there a way to attach chimera_send to an existing process of chimera, or have chimera quit at the end of a chimera_send call? Thanks, Charlie > -----Original Message----- > From: chimera-users-bounces(a)cgl.ucsf.edu > [mailto:chimera-users-bounces@cgl.ucsf.edu] On Behalf Of Thomas Goddard > Sent: Monday, May 24, 2004 6:20 PM > To: heiland(a)indiana.edju > Cc: chimera-users(a)cgl.ucsf.edu > Subject: [Chimera-users] non-GUI use of Chimera > > Hi Randy, > > If I understand correctly you want to be able to type some Chimera > Python > code or Chimera commands to a web browser, those get sent to the web > server > which sends them to Chimera running on the web server machine. Then you > capture the Chimera image and send it back to the web browser. > > I believe this can be done using the Chimera "Web Data" facility. > This > allows programs to send commands to Chimera. The intended use is for > displaying models in Chimera when the user clicks a link in a web > browser. > It works by having the *.chimerax file associated with the link lauch > helper application <chimera-install-location>/bin/chimera_send. > The same mechanism allows any program to send commands to Chimera. > The Chimera User's Guide describes this under Tools/Web Browser > Configuration. > Here's a URL: > > > http://www.cgl.ucsf.edu/chimera/1.1951/docs/ContributedSoftware/webdata/ > chimerax.html > > What you might do is have your web server create an appropriate > *.chimerax that executes the user's commands and then prints the > Chimera graphics window to a file. The web server than sends the > image (jpg) to the web client. There is a bunch of somewhat tricky > glue to make this work. As previous posts mentioned, Chimera will > need to pop up windows on some screen to render an image. For the web > server to know when Chimera has finished making an image, the > *.chimerax file might end in a command that writes to some secondary > status file after the image file is written. That will tell the web > server that the image is ready. > > We haven't tried this. And there are other pitfalls, like you are > likely to run into trouble if the web server starts 2 instances of > Chimera > and then tries to dispatch commands with the chimera_send program. > > Please let me know if you get this to work. Feel free to send > additional > questions. > > Tom > _______________________________________________ > Chimera-users mailing list > Chimera-users(a)cgl.ucsf.edu > http://www.cgl.ucsf.edu/mailman/listinfo/chimera-users > > > >

3 5

vrml 2.0 export of pdb representation
by Pattanayek, Sabuj 25 May '04

25 May '04

Hi, Is there anyway to export a pdb model into vrml once I have it just the way I want it (ribbon & atom sizes, colors, etc)? I tried a pdb2vrml converter from http://www.pc.chemie.tu-darmstadt.de/research/vrml/pdb2vrml_right.shtml but freewrl crashes when loading the output (I think pdb2vrml might be outputting vrml 1.0), it's also not very versatile as far as what I want to display into vrml from the pdb. Thanks

3 2

non-GUI use of Chimera
by Thomas Goddard 24 May '04

24 May '04

Hi Randy, If I understand correctly you want to be able to type some Chimera Python code or Chimera commands to a web browser, those get sent to the web server which sends them to Chimera running on the web server machine. Then you capture the Chimera image and send it back to the web browser. I believe this can be done using the Chimera "Web Data" facility. This allows programs to send commands to Chimera. The intended use is for displaying models in Chimera when the user clicks a link in a web browser. It works by having the *.chimerax file associated with the link lauch helper application <chimera-install-location>/bin/chimera_send. The same mechanism allows any program to send commands to Chimera. The Chimera User's Guide describes this under Tools/Web Browser Configuration. Here's a URL: http://www.cgl.ucsf.edu/chimera/1.1951/docs/ContributedSoftware/webdata/chi… What you might do is have your web server create an appropriate *.chimerax that executes the user's commands and then prints the Chimera graphics window to a file. The web server than sends the image (jpg) to the web client. There is a bunch of somewhat tricky glue to make this work. As previous posts mentioned, Chimera will need to pop up windows on some screen to render an image. For the web server to know when Chimera has finished making an image, the *.chimerax file might end in a command that writes to some secondary status file after the image file is written. That will tell the web server that the image is ready. We haven't tried this. And there are other pitfalls, like you are likely to run into trouble if the web server starts 2 instances of Chimera and then tries to dispatch commands with the chimera_send program. Please let me know if you get this to work. Feel free to send additional questions. Tom

1 0

Re: [Chimera-users] vrml 2.0 export of pdb representation
by Pattanayek, Sabuj 22 May '04

22 May '04

imho, 3DS is currently more widely supported than X3D by expensive, cheap, and free 3d modellers. Most 3D format translators can convert from 3DS to other types (DXF, iv, etc). However it looks like from web3d.org that X3D is much more versatile and feature rich than old 3ds... Quoting Greg Couch <gregc(a)cgl.ucsf.edu>: > We currently only support VRML export of volume isosurfaces and multiscale > surfaces. So the answer is no, not yet. > > Exporting to a 3D graphics file format is definately something we would > like to do. However, it is unclear what format would be best to generate. > VRML 2.0 will never support transparency like chimera does. VRML's > successor, X3D, looks promising, but there aren't many users yet. Other > possiblities include the Renderman format, popular with animation studios; > the .3DS format, popular on PCs; and many others. > > So, a question for everyone reading this, if you had to pick one format, > would it be VRML97? And if so, why? And if not, why not? We are leaning > torwards outputing X3D, but would appreciate comments from people who > would actually be using the output. > > Greg Couch > UCSF Computer Graphics Lab > gregc(a)cgl.ucsf.edu > > > On Fri, 21 May 2004, Pattanayek, Sabuj wrote: > > > Is there anyway to export a pdb model into vrml once I have it just the way > I > > want it (ribbon & atom sizes, colors, etc)? I tried a pdb2vrml converter > from > > http://www.pc.chemie.tu-darmstadt.de/research/vrml/pdb2vrml_right.shtml > but > > freewrl crashes when loading the output (I think pdb2vrml might be > outputting > > vrml 1.0), it's also not very versatile as far as what I want to display > into > > vrml from the pdb. >

1 0

Color blending on ribbons
by Emma Hill 20 May '04

20 May '04

Hi, I wondered if there is any easy way to blend colors from one residue to the next on a ribbon chain (as in the rainbow facility.) I want to be able to precisely define where this takes place (ie whenever I switch from one color to another) - more specifically I am coloring different regions of a superimposed pair to show where is aligned and where is not, and the drastic switch from one color to the next would definitely look better more blended. Thanks in advance Emma -- Emma Hill Brenner Research Lab U.C. Berkeley

2 1

non-GUI use of Chimera
by Randy Heiland 20 May '04

20 May '04

I am evaluating Chimera as a tool for providing (pseudo-interactive) vis functionality from a web server and therefore would like to use it in a non-gui mode. I saw the earlier post/reply to this list (May 4) on this topic, and Eric indicated that one way to do this was to put all your cmds in a file then run "chimera --nogui <cmd-file>". Obviously this isn't ideal for my purposes. So, I'm wondering what the alternatives are? I'm a newbie to Chimera, but I've read how it was based on Midas, and so I decided to see how far I could get in reading in a pdb file and... I didn't get very far: (after install Chimera (Linux, beta 1 build 1951), plus dependent Python pkgs) %python >>> import chimera >>> import Midas >>> f=Midas.open('foo.pdb') and I get an error about 'missing default material' and realize I'm doing the wrong thing anyway. Basic question: can I, from the Python interpreter, use Chimera/Midas to read in a pdb file and do basic Chimera functionality? For now, I don't even care if it assumes X is running and pops up windows. Thanks, Randy

2 3

rainbow command with sub-models
by Jarrod Smith 19 May '04

19 May '04

Hello all, Is it possible to add a new "level" to the rainbow command such that each submodel of a multi-model PDB file (MODEL/ENDMDL records) can be visited separately by the rainbow command? I have some pdb files with several tens to hundreds of submodels in them. Separating them by chain ID (so I can use the rainbow chains command) has proven to be impractical for other reasons, so this would be useful in such a case. Thanks for the great software. Jarrod Smith -- Jarrod A. Smith, Ph.D. Asst. Director, Center for Structural Biology Research Asst. Professor, Biochemistry Vanderbilt University 203 processes: 140 sleeping, 63 running

2 4

new Chimera production release
by Eric Pettersen 14 May '04

14 May '04

Hi all, A new Chimera production release (version 1.1951) is now available for download. I've summarized the major additions since the last snapshot (1.1917) below. A list of all the changes since the last production release (1.1872) is in the release notes: http://www.cgl.ucsf.edu/chimera/docs/relnotes/1.1951.html Eric Pettersen UCSF Computer Graphics Lab pett(a)cgl.ucsf.edu http://www.cgl.ucsf.edu Changes since 1.1917: General changes --------------- Newly opened models are shown larger (as well as models after a 'window' or 'focus' command) Selections can be written out as parsable text files (Actions->Write... or from the Selection Inspector) New residue attributes: kdHydrophobicity (Kyte-Doolittle hydropathy, amino acids only) mavPercentConserved (% conserved in a sequence alignment open in Multalign Viewer) above for use with Render/Select by attribute tools Double-picking a bond brings up a menu with Rotate Bond and Select Bonded (the latter is new and selects the atoms); double-picking a pseudobond bring up a menu with just Select Bonded; "select bonded" also now in Pseudobond Panel As much as possible, can handle PDB files with 100000+ atoms (e.g. AMBER output). Some PDB records (e.g. CONECT records) have no room for 6-digit atom serial numbers, so only can be used for the first 99,999 atoms. A new tutorial describing preparing images for publication Select by Attribute Value... in Selection menu and Model Panel New Tools --------- Lighting (Viewing Parameters): change and save lighting parameters Color Secondary Structure (Graphics): color peptides/proteins by secondary structure Render Attributes (Graphics): color or otherwise change atom, residue, or model displays based on attribute values; was in 1.1917 release but several changes have been made and documentation added Demo (its own category; entries are individual demos): create/replay demos in Chimera; example demo included (Cyclooxygenase) Tool Changes ------------ Adjust Torsions now allows a bond to be rotated with mouse motions in the graphics window Movie has File->Save PDB for writing out either the current frame or all frames that have been viewed Multalign Viewer: ClustalX-style residue letter coloring added and made the default; can save alignment depiction as an EPS file; can hide Consensus and/or Conservation lines RibbonStyleEditor now has two panels in which scalings and cross-sections can be defined and named; named cross-section are analogous to the built-in cross-sections (flat, edged, round) and can be accessed from the Actions menu or with ribrepr; nucleic acids are now treated separately in the scalings panel New "Graphics" category; Per-Model Clipping, PipesAndPlanks, ResProp, and Ribbon Style Editor moved into it SimpleSession saves lighting, ribbon scaling and cross-section(s), B-factor and occupancy values Web Data asks for confirmation before accepting certain files; associated option in Preferences Command Changes --------------- alias: $1, $2, $3, ... may be used to indicate the first, second, third... arguments of the alias conic: works on Windows, and the -s flag now works on Linux (besides Mac, SGI) write: has "trajectory" keyword to write all frames that have been viewed (loaded) [analogous option added to Model Panel's "write PDB" dialog] ribbackbone: now applies to nucleic acids also

1 0

by riccardodesantis＠fastwebnet.it 05 May '04

05 May '04

Hallo, my name is Riccardo De Santis and I'm a new user of CHIMERA. I'd like to know how run the program from the command line. I've tryed with the command NOGUI but the program doesn't start. Another question is: how can I use matchmaker (in the graphic menu) to match the molecules from the command line? Which is the istruction to type? Thank you in advance, Riccardo De Santis.

2 1

by riccardodesantis＠fastwebnet.it 05 May '04

05 May '04

Hallo, my name is Riccardo De Santis, and I'm a new user of CHIMERA. I'd like to know how to use the program in command line mode. I noticed that there is a NOGUI option, I tryed it but unsuccessfully. Cuold you help me , please? Thank you in advance, Riccardo De Santis.

1 0

Coloring by bfactor
by Thomas Goddard 23 Apr '04

23 Apr '04

Hi Padma, We have developed a Chimera tool to color by bfactor or other attributes but it is not yet released. We expect to have a new release in the middle of May. Here are a couple approaches you could use now. You can use the color command like the following to color all atoms with bfactor between 10 and 15 green. color green @/bfactor>=10 and bfactor<=15 You could put a series of these commmands in a file to color several intervals different colors to get a rainbow effect. Use the "open" command to run the file of coloring commands. A more flexible route is to use some Python code. Putting the code shown below in a file bcolor.py and then opening this file will make a rainbow coloring by bfactor for the selected atoms. A fancier version of this could have a dialog popup to specify bfactor range and colors. If you make a molecular surface with Actions/Surface/Show the default coloring of this surface is to use the atom colors. So coloring the atoms by any of the above means will allow you to get the same coloring on a surface. Tom ----- code bcolor.py follows: # ----------------------------------------------------------------------------- # def color_atoms_by_bfactor(atoms, b0, rgb0, b1, rgb1): for a in atoms: if hasattr(a, 'bfactor'): b = a.bfactor if b < b0: f = 0 elif b < b1: f = (b - b0) / (b1 - b0) else: f = 1 r,g,b = ((1-f)*rgb0[0] + f*rgb1[0], (1-f)*rgb0[1] + f*rgb1[1], (1-f)*rgb0[2] + f*rgb1[2]) import chimera a.color = chimera.MaterialColor(r,g,b) # ----------------------------------------------------------------------------- # from chimera import selection atoms = selection.currentAtoms() # currently selected atoms b0 = 5 rgb0 = (0,0,1) # blue b1 = 20 rgb1 = (1,0,0) # red color_atoms_by_bfactor(atoms, b0, rgb0, b1, rgb1) > To: Thomas Goddard <goddard(a)cgl.ucsf.edu> > Subject: One more Chimera question! > From: Padmaja Natarajan <nataraja(a)scripps.edu> > Date: Fri, 23 Apr 2004 13:31:38 -0700 > > Hi Tom, > > Is it possible to make a molecular surface based on PDB coordinates > and be able to view this surface colored (as a spectrum going from blue > for smallest to red for largest value), > depending on values stored in say B-factor column of the PDB file. > (as shown in figure below created using AVS) > > thanks again > regards > Padma

2 1

Re: [Chimera-users] CHIMERA help - coloring density maps
by Thomas Goddard 23 Apr '04

23 Apr '04

Hi Padma, There is a Chimera command called "longbond" described in the Chimera User's Guide / Basic Functions / Commands / longbond http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/longbond.html that removes long bonds. There is currently no radial masking tool. But using the Isosurface Colorizer extension from EMAN you can set the color to transparent black outside or inside a certain radius. That will make it invisible. For example you can make the color blue at radius 10 nm, and transparent black at 10.01 nm. You will also want to turn on the Volume Viewer / Options / Dim transparent surface mesh switch to get rid of the specular highlights on the transparent black part of the surface. This trick for masking does not "cap" the isosurface at the radius where it is cut. So you see the inside of the isosurface. Turning off two sided lighting in Volume Viewer / Options will make the interior of the surface dark. To cut and cap the isosurface at a specific radius you really want to the display to act as if the volume data were zero outside a specified radius. I don't have any way to do that, although it probably would not be too hard to hack something into surface.py that copied the 3-d matrix and made it zero outside a specified radius. By the way, the modification I described to Chimera surface.py to make the Isosurface Colorizer from EMAN work is not needed in Chimera 1917. That change was included in Chimera 1917. I you are looking at models of virus particles you might be interested in the Chimera Multiscale extension. Here are some example images: http://www.cgl.ucsf.edu/home/goddard/viper/ The images at the bottom of that page show the RNA in satellite tobacco mosaic virus. Tom > From: Padmaja Natarajan <nataraja(a)scripps.edu> > Subject: Re: [Chimera-users] CHIMERA help - coloring density maps > Date: Fri, 23 Apr 2004 13:08:02 -0700 > To: Thomas Goddard <goddard(a)cgl.ucsf.edu> > > Hi Tom, > > Thank you very much for the details you had given me earlier in > making the radial coloring work for surfaces. > I , with the help of my colleague Jeff who made the changes in the > chimera (version 1917) script, am > able to produce radial coloring for my density map! > > Can ask more questions?? > > 1) I have coordinates for an icosahedron, but when I read it into > chimera, > several un-wanted connections are made (probably because of the way the > PDB file is written - for use in the graphics program O). I do not see > any distance cut off > criterion command that would help (probably) remove some of those > connections! > > 2) density map: I would like a radial masking tool, that can be used to > choose > minimum and maximum values for radius and be able to see the map only > within this > radius range. For example, this helps in looking at the inner (RNA) > density in a virus structure. > > thanks again! > regards > Padma > ______________________________ > Padmaja Natarajan, Ph.D > The Scripps Research Institute > La Jolla, California, USA.

1 0

install question
by Randy Heiland 23 Apr '04

23 Apr '04

Greetings from a newbie to Chimera, Is there a searchable archive? (I found the archive, but no associated search engine) I'd really like to be able to install Chimera without blowing away my current Python installation - is this possible? (I'll be using both Windows and Linux) Thanks, Randy

2 1

Viewing orbitals and charge distributions
by João Otávio 23 Apr '04

23 Apr '04

Hello All! I am a CHIMERA newbie and I would like to know how could I create a module in order to read orbitals and/or charge distributions from molecular orbitals programs such as GAUSSIAN, GAMESS, etc. Thanks in advance, .................................................................João -- ___________________________________________________________ Sign-up for Ads Free at Mail.com http://promo.mail.com/adsfreejump.htm

2 1

save session error in Chimera 1892
by Juha Huiskonen 22 Apr '04

22 Apr '04

Hi, when trying to save session in chimera 1892 I get an error: Error processing trigger "save session": ValueError: Can't convert type <type 'bool'> (see reply log for Python traceback info) Reply log: Traceback (most recent call last): File "/usr/local/chimera1892/share/__main__.py", line 56, in ? value = chimeraInit.init(sys.argv) File "/usr/local/chimera/share/chimeraInit.py", line 285, in init tkgui.eventLoop() File "/usr/local/chimera/share/chimera/tkgui.py", line 2809, in eventLoop app.mainloop() File "/usr/local/chimera/lib/python2.3/lib-tk/Tkinter.py", line 965, in mainloop self.tk.mainloop(n) File "/usr/local/chimera/lib/python2.3/site-packages/Pmw/Pmw_1_2/lib/PmwBase.py", line 1747, in __call__ return apply(self.func, args) File "/usr/local/chimera/share/chimera/baseDialog.py", line 409, in OK self.Apply() File "/usr/local/chimera/share/OpenSave/__init__.py", line 141, in Apply self.command(1, self) File "/usr/local/chimera/share/SimpleSession/gui.py", line 24, in _saveCB saveSession(paths[0]) File "/usr/local/chimera/share/SimpleSession/save.py", line 269, in saveSession chimera.triggers.activateTrigger(SAVE_SESSION, buf) File "/usr/local/chimera/share/chimera/triggerSet.py", line 214, in activateTrigger self._triggerDict[name].activate(data) File "/usr/local/chimera/share/chimera/triggerSet.py", line 127, in activate ret = handler.invoke(self._name, triggerData) File "/usr/local/chimera/share/chimera/triggerSet.py", line 82, in invoke self._funcData, triggerData) File "/usr/local/chimera/share/VolumeViewer/__init__.py", line 138, in save_session_cb session.save_volume_state(self, file) File "/usr/local/chimera/share/VolumeViewer/session.py", line 22, in save_volume_state t = objecttree.instance_tree_to_basic_tree(s) File "/usr/local/chimera/share/SessionUtil/objecttree.py", line 69, in instance_tree_to_basic_tree d[k] = it2bt(getattr(itree, attr)) File "/usr/local/chimera/share/SessionUtil/objecttree.py", line 37, in it2bt return instance_tree_to_basic_tree(it, c, a) File "/usr/local/chimera/share/SessionUtil/objecttree.py", line 72, in instance_tree_to_basic_tree raise ValueError, "Can't convert type " + str(t) ValueError: Can't convert type <type 'bool'> -- The session file is written, but it doesn't contain the files I have read into chimera. Cheers, Juha Huiskonen

2 1

Re: [Chimera-users] chimera and hydrophobic potentials
by Thomas Goddard 13 Apr '04

13 Apr '04

Hi Miguel, If GRID computes the hydrophobic potential on a surface, it seems like a bad strategy to save it in a volume format. If you interpolate onto a slightly different surface, you will get values interpolated with zeros in the volume because the hydrophobicity wasn't computed at most grid vertices. It seems it would be better to render the hydrophobicity values with exactly the GRID surface. Can GRID output the triangulated surface with hydrophobicity values? I looked at the MSMS library. It does not look like it provides a triangulated version of the solvent accessible service needed by Chimera. It can give you the surface area contributed by each atom to the solvent accessible surface. Tom > Date: Sun, 11 Apr 2004 22:39:40 +0100 > From: Miguel Ortiz <mol(a)ysbl.york.ac.uk> > To: Thomas Goddard <goddard(a)cgl.ucsf.edu> > Subject: Re: [Chimera-users] chimera and hydrophobic potentials > > Dear Tom, > > I've tested your cns2delphi script and, as expected, it works wonderful- > ly. The problem with colouring the surface using these maps actually > comes from the fact that for the hydrophobic probe, grid only calculates > the potential on what it's called the "solvent accessible surface" and > not on the whole volume (as delphi does for electrostatics and grid for > other probes). Being this the case, the best solution would be to > render in chimera this surface (at present, msms renders the "solvent > excluded surface" defined as you said in your email) > > I will try to find out how msms renders its surfaces, because the only > difference between "solvent accessible " and "solvent excluded" is that > the first is the "surface" defined by the center of the probe, while > the second (msms one) is the contact surface between the surface of the > probe and that of the protein (as defined by the Van der Waals radii of > its atoms). And it seems to me that it shouldn't be difficult to convert > from one to the other. Or perhaps I'm wrong and it's more tricky... > > Thank you again for all your help. It has been not only useful but > also instructive. > > Best regards, > > Miguel

1 0

Re: [Chimera-users] chimera and hydrophobic potentials
by Thomas Goddard 09 Apr '04

09 Apr '04

Hi Miguel, I wrote a conversion program cns2delphi.py that takes the k2f output and produces a delphi file and tried it on the data you sent. The delphi file is read properly by Chimera DelphiViewer and can be used to color an MSMS surface. A problem with the data you sent is that the hydrophobicity map is zero over the entire MSMS surface. There are blobs of hydrophobicity a few angstroms away from the surface but they don't touch the surface. I observed this by showing the map with the Chimera volume viewer (just use File/Open...). You can change the probe radius for generating the MSMS surface by going to Favorites/Model Panel, selecting the MSMS surface, clicking the Attributes button, changing the probe radius setting, and pressing Enter. (You need to set the surface color to a solid color before changing the radius or an error message is generated.) Unfortunately that does not push the molecular surface out uniformly -- only the surface in the cavities is pushed further from the molecule. Apparently the probe radius is the radius of a ball rolled over the molecule with standard atom radii. The MSMS surface is the surface of the region that cannot be touched by the ball. So a bigger ball does not change the surface location at protrusions. Maybe a better approach is to specify larger atom radii. Below is the cns2delphi.py Python program. It requires the Numeric Python module which is not a standard part of Python. The Python that comes with Chimera includes Numeric. You can convert k2f output to delphi format as follows. % /usr/local/chimera/bin/python2.3 cns2delphi.py 1ltlA.k2f 1ltlA.phi This creates the 1ltlA.phi file. To use the Chimera DelphiViewer extension, start Chimera, open 1ltlA.pdb, make a molecular surface with Actions/Surface/Show, start Extensions/Delphi/DelphiViewer, enter Potential File 1ltlA.phi, click the Browse button next to MSMS Surface in the Molecular surface pane so the surface is listed, then press Apply in the Molecular surface pane. This will color the surface. You can also just use File/Open... on 1ltlA.phi to see an isosurface of the hydrophobicity map. You can drag the vertical bar on the histogram shown in the volume viewer dialog to change the isosurface threshold level. Tom ----- cns2delphi.py follows. Preserve indentation since it is important in Python. #!/usr/bin/env python # ---------------------------------------------------------------------------- # Convert CNS volume to Delphi format. # # ---------------------------------------------------------------------------- # Returns matrix, xyz_origin, and xyz_step # def read_cns_file(path): f = open(path, 'r') f.readline() # First line is blank ntitle_line = f.readline() # integer number of comment lines try: ntitle = int(ntitle_line.split()[0]) except ValueError: raise SyntaxError, ('Invalid CNS comment line count on line 2: %s' % ntitle_line[:80]) for t in range(ntitle): f.readline() extent_line = f.readline() extent = split_fields(extent_line, 8, 9) try: na, amin, amax, nb, bmin, bmax, nc, cmin, cmax = map(int, extent) except ValueError: raise SyntaxError, 'Invalid CNS grid size line: %s' % extent_line[:80] cell_line = f.readline() cell = split_fields(cell_line, 12, 6) try: cell_params = map(float, cell) except ValueError: raise SyntaxError, ('Invalid CNS cell parameters line: %s' % cell_line[:80]) cell_size = cell_params[0:3] cell_angles = cell_params[3:6] axis_order_line = f.readline() # Must be ZYX asize, bsize, csize = (amax - amin + 1, bmax - bmin + 1, cmax - cmin + 1) if asize <= 0 or bsize <= 0 or csize <=0: raise SyntaxError, ('Bad CNS grid size (%d,%d,%d)' % (asize, bsize, csize)) import Numeric data = Numeric.zeros((csize, bsize, asize), Numeric.Float32) for c in range(csize): read_floats(f, 1) # c section number values = read_floats(f, asize * bsize) array = Numeric.array(values, Numeric.Float32) plane = Numeric.reshape(array, (bsize, asize)) data[c,:,:] = plane f.readline() # footer - int value must be -9999 f.readline() # density average and standard dev f.close() xyz_step = (cell_size[0]/na, cell_size[1]/nb, cell_size[2]/nc) xyz_origin = map(lambda gmin, step: gmin * step, (amin, bmin, cmin), xyz_step) return data, xyz_origin, xyz_step # ----------------------------------------------------------------------------- # Read ascii float values on as many lines as needed to get count values. # def read_floats(f, count): values = [0.0] * count c = 0 while c < count: line = f.readline() if line == '': raise FileFormatError('File is shorter than expected') fields = split_fields(line, 12, 6) if c + len(fields) > count: fields = fields[:count-c] for v in map(float, fields): values[c] = v c = c + 1 return values # ----------------------------------------------------------------------------- # def split_fields(line, field_size, max_fields): fields = [] for k in range(0, len(line), field_size): f = line[k:k+field_size].strip() if f: fields.append(f) else: break return fields[:max_fields] # ---------------------------------------------------------------------------- # def write_delphi_file(path, matrix, xyz_center, scale): zsize, ysize, xsize = matrix.shape if xsize != ysize or xsize != zsize: raise ValueError, "Matrix size (%d,%d,%d) has to be the same in all 3 dimensions" % (zsize, ysize, xsize) size = xsize uplbl = delphi_record('Converted from ' + path) morelabels = delphi_record('%70s' % '') import Numeric data_flat = Numeric.ravel(matrix) data = delphi_record(floats_as_string(data_flat)) botlbl = delphi_record('') params = delphi_record(floats_as_string((scale,) + tuple(xyz_center)) + integer_as_string(size)) f = open(path, 'wb') f.write(uplbl + morelabels + data + botlbl + params) f.close() # ---------------------------------------------------------------------------- # def delphi_record(data): size_str = integer_as_string(len(data)) r = size_str + data + size_str return r # ---------------------------------------------------------------------------- # def integer_as_string(i): import Numeric return Numeric.array((i,), Numeric.Int32).tostring() # ---------------------------------------------------------------------------- # def floats_as_string(f): import Numeric return Numeric.array(f, Numeric.Float32).tostring() # ---------------------------------------------------------------------------- # def cns_to_delphi_file(cns_path, delphi_path): data, xyz_origin, xyz_step = read_cns_file(cns_path) if xyz_step[0] != xyz_step[1] or xyz_step[0] != xyz_step[2]: raise Value_Error, 'Step sizes (%.2g, %.2g, %.2g) are not the same along all axes' % xyz_step scale = 1.0 / xyz_step[0] size = list(data.shape) size.reverse() xyz_center = map(lambda o, st, sz: o + st * ((sz-1)/2), xyz_origin, xyz_step, size) write_delphi_file(delphi_path, data, xyz_center, scale) # ---------------------------------------------------------------------------- # if __name__ == '__main__': import sys if len(sys.argv) != 3: sys.stderr.write('Syntax: cns2delphi <cnsfile> <delphifile>\n') sys.exit(1) cns_path = sys.argv[1] delphi_path = sys.argv[2] cns_to_delphi_file(cns_path, delphi_path)

1 0

Re: [Chimera-users] chimera and hydrophobic potentials
by Thomas Goddard 08 Apr '04

08 Apr '04

Hi Miguel, Send me some example k2a output and I'll send you a little Python program that will convert it into a Delphi file readable by Chimera DelphiViewer. Tom > Date: Thu, 8 Apr 2004 19:33:31 +0100 > From: Miguel Ortiz =?iso-8859-15?Q?Lombard=EDa?= <mol(a)ysbl.york.ac.uk> > Subject: Re: [Chimera-users] chimera and hydrophobic potentials > > Hi Tom, > > Thank you very much for your comments. > Part of the problem is that I don't trust the "delphi" maps that > GRID produces. Actually, I've been able to open some of them with > chimera. Other programs ablo to open such maps (like vmd) can't > open any of them. > > Yes, I've used k2a, and it works: it produces an ascii map that > looks sensible. My problem here is that my programming skills are > very reduced! I could work it out on a "model" program, but not > from scratch. Have you something like that? > > Thank you again, your answer has been very useful. > > Best regards, > > Miguel > > -- > Miguel Ortiz Lombardía > email: mol(a)ysbl.york.ac.uk > http://www.ysbl.york.ac.uk/~mol > ----------------------------------------------------------------------- > Our business in life is not to succeed, but > to continue to fail in good spirits. R.L. Stevenson >

1 0

chimera and hydrophobic potentials
by Miguel Ortiz Lombardía 08 Apr '04

08 Apr '04

Hi, I would like to display the hydrophobic potential of a protein onto its MSMS surface. One program that calculates such potentials is GRID (from Molecular Discovery: http://www.moldiscovery.com/ ) GRID can produce the map in different formats. One of them is like delphi (grasp mode) map. This has the drawback that the map has to be 65x65x65 grid points and the grid points separation has to be an integer starting with 1 Ansgtroem. Hence, it's not appropriate for medium-big proteins. The GRID format hasn't those limitations. I can __sometimes__ read in chimera (throught the DelphiViewer) the delphi maps produced by GRID. Sometimes they are unreadable (and this seems to me more a problem with GRID than with chimera: something wrong in the header) But even when they are readable I'm not sure they are covering the region they are supposed to be covering. Does anyone know if chimera can read any other of the GRID map formats: GRID (this would be the optimal) Chem-X Frodo Insight Molcad Quanta Sybyl Otherwise: does anybody how to the kind of representation I'm looking for in any other way? Thank you for any advice! Miguel -- Miguel Ortiz Lombardía email: mol(a)ysbl.york.ac.uk http://www.ysbl.york.ac.uk/~mol ----------------------------------------------------------------------- Our business in life is not to succeed, but to continue to fail in good spirits. R.L. Stevenson

2 1

compiling minrms on mac os10.2.8
by Chris Adams 02 Apr '04

02 Apr '04

Hi: I am also trying to get minrms on my OS10.2.8 machine. I followed the commands you recommended to Howard to get the "make" command to work. However, my ls output was not at all like what you said would happen in the message. I get this: ls -1 /usr/bin/*make /usr/bin/automake and am unable to compile the minrms program. I am not very unix savvy but am learning, perhaps you can help? thank you, Chris Adams UW-Biochemistry Madison, WI 53704

2 1

correction re sessions
by Elaine Meng 01 Apr '04

01 Apr '04

Hello, in my previous message about saving parts of PDB files, I included some text from older replies to users, and this part about sessions is no longer true: >The command "save" (or File... Save Session As... ) doesn't save coordinates, >only how they are positioned/colored/displayed, etc. In fact, restarting >a saved session will use the original input file(s) for coordinates. Actually in the last few versions, the session files do include coordinates and do not look for the original input files. Sorry about that, Elaine

1 0

writing out part of a PDB
by Elaine Meng 01 Apr '04

01 Apr '04

Hi Shuanghong, Question (1), surface problem: I don't think the problem where nothing happens when you type the command "surface" has been reported. Sometimes there is a delay while the surface is being calculated, but it shouldn't be terribly long. Problems have been reported in some cases, but they generally cause some sort of error message. To register a problem more formally, you could send mail to chimera-bugs(a)cgl.ucsf.edu and be sure to include what version of Chimera and what kind of computer you are using, and what structure(s) cause the problem. Question (2), writing out part of a PDB file: If you are on a **non-Windows** system, you can use the command "pdbrun" to save just the atoms that are displayed. You could display only the atoms you want by first selecting what you want to measure the distance to. Then, use Select...Zone... to select all the atoms in the region you want to write out. Then the command "show sel" will undisplay the rest of the stuff. (Or, instead of selecting you could do it all in one command, for example show :212.a za<4.5 for all atoms within 4.5 A of any atom in residue 212 of chain A. Specifying zones in the command line is described in the manual page on atom specification, which is listed in the Basic Functions index.) Next, use pdbrun (see man page using the command "help pdbrun"), for example: pdbrun nouser cat > myzone.pdb Pdbrun literally sends the information for only the displayed atoms in their displayed positions to another command, in this case "cat" to send the information to the file myzone.pdb. If you want the coordinates to be the same as in the original file, first use the command "reset." The command "save" (or File... Save Session As... ) doesn't save coordinates, only how they are positioned/colored/displayed, etc. In fact, restarting a saved session will use the original input file(s) for coordinates. The command "write" (same as the function "write PDB" in the Model Panel) saves coordinates into a file. However, these will save the entire coordinates whether displayed or not. If you are on **Windows** the only way I can think of is to actually delete the atoms/bonds you don't want to write out ("delete" command or Actions... Atoms/Bonds... delete). For example, when only what you want to write out is displayed, delete @/!display actually deletes the other stuff. Subsequently, a file created with "write" would only contain the atoms/bonds that had not been deleted. I hope this helps, Elaine =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-= Elaine C. Meng, Ph.D. Computer Graphics Lab and Babbitt Lab meng(a)cgl.ucsf.edu http://www.cgl.ucsf.edu/home/meng/index.html =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=

1 0

crashing when calculating larger surfaces
by Elaine Meng 01 Apr '04

01 Apr '04

>From eperozo(a)virginia.edu Wed Mar 31 15:15:53 2004 >Sorry to bother you, but I'm having terrible problems with the generation of >surfaces. The program crashes as soon as I request a surface to be built (a >tetramer with ~250 residues/monomer). Interestingly, I don't seem to have >problems with a smaller protein (120 residues). When it crashes it gives me >exit code - 1073741819. Tried two versions of chimera in two different >computers. Same result. >Is this a known bug? >Thanks in advance. >Eduardo Hi Eduardo, It is a known problem that the surface routine is a memory hog. If you are on a unix-type machine it might help to type "unlimit" in the shell before starting Chimera. Otherwise, it might just be trying to use more memory than your machines have. We have been using code developed elsewhere, MSMS, for these surfaces. Due to memory/robustness problems, we are generating some code in-house which will eventually replace MSMS. However, this is still in progress, and has not been incorporated into any Chimera release. You could try deleting extra monomers before using surface or instead show the low-resolution surface generated by the Multiscale Models tool (although it builds multimers with matrix info, even if there is no matrix info, it will still show the low-res surface for the structure you have open). Tools...Multiscale... Multiscale Models, click the "Make models" button. Then if you select All chains (button near the top) you can change the Resolution value and then click Resurface to see different levels of detail. There is a man page for this tool in the Tools section of the User's Guide if you want more info. If you want to register the problem more formally, please send a message to chimera-bugs(a)cgl.ucsf.edu including version, platform, and an example or few of which structures (e.g. their PDB IDs) cause the problem. Best, Elaine

1 0

Re: [Chimera-users] Rotating molecules
by Elaine Meng 30 Mar '04

30 Mar '04

Hi Eduardo, yes, the command "move" does translations by a specified amount and the command "turn" rotates by a specified amount. Type help move in the command line to see the man page for move, and help turn for the man page for turn. You might also want to take a look at the commands "rock" and "roll" Hope this helps, Elaine =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-= Elaine C. Meng, Ph.D. Computer Graphics Lab and Babbitt Lab meng(a)cgl.ucsf.edu http://www.cgl.ucsf.edu/home/meng/index.html =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=

1 0

Rotating molecules
by Eduardo Perozo 30 Mar '04

30 Mar '04

Hello, Is there a way to rotate/translate a molecule by a specific amount? Clearly I can do "freehand" rotations with the mouse but I will like to define rotations and translation numerically. Thanks! Eduardo

1 0

Re: [Chimera-users] export pdb file
by Elaine Meng 30 Mar '04

30 Mar '04

Hi Shuanghong, there are two ways to export a PDB file. One is by choosing the "write PDB" function from the Model Panel. The other is by using the "write" command. Each way has an option to save the transformed coordinates or the coordinates relative to the untransformed coordinates of a particular model (take a look at the User's Guide for detailed descriptions). Hope this helps! Elaine =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-= Elaine C. Meng, Ph.D. Computer Graphics Lab and Babbitt Lab meng(a)cgl.ucsf.edu http://www.cgl.ucsf.edu/home/meng/index.html =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=

2 1

Re: [Chimera-users] v1.1917 OpenGL misconfiguration detected.
by Greg Couch 30 Mar '04

30 Mar '04

Got it. The /usr/X11R6/lib/ligGL.so.1 file should have been removed as part of the NVidia driver installation. We just confirmed that with the Red Hat 9 and the NVIDIA-Linux-x86-1.0-5336-pkg1.run driver distribution. You had version 53.28, so maybe it was a bug in that driver? We had previous installed NVidia version 43.63 using the same mechanism, so it's probably not a problem with the driver, but something went wrong during the installation. You should uninstall and install that driver again (or just upgrade to the 53.36). Good luck, Greg On Mon, 29 Mar 2004, Pattanayek, Sabuj wrote: > Here is the output > > -rwxr-xr-x 1 root root 367092 Jan 3 12:17 /usr/lib/libGL.so.1 > -rwxr-xr-x 1 root root 264888 Feb 12 13:22 /usr/X11R6/lib/libGL.so.1 > > I am using RH9. > > Quoting Greg Couch <gregc(a)cgl.ucsf.edu>: > > > Chimera checks /usr/lib/libGL.so.1 and /usr/X11R6/lib/libGL.so.1 to see > > if they point to the same file or not. If they do not, then that indicates > > that there was an OpenGL driver installation error -- sometimes resulting > > in chimera not using the hardware accelerated OpenGL driver. > > > > Could you send the output from: > > > > ls -Ll /usr/lib/libGL.so.1 /usr/X11R6/lib/libGL.so.1 > > > > I'd like to confirm that they are different. And which version of Linux > > are your running -- Red Hat 8.0? Mandrake? Suse? Debian? We'd like to > > try to recreate the problem here. > > > > Greg Couch > > UCSF Computer Graphics Lab > > gregc(a)cgl.ucsf.edu > > > > On Mon, 29 Mar 2004, Pattanayek, Sabuj wrote: > > > > > Finally got around to installing the latest development linux version. At > > > startup I'm getting the following messages: > > > > > > OpenGL misconfiguration detected. > > > Please reinstall graphics driver. > > > > > > v1.1872 & v1.1892 still work on the machine with no problems. > > > > > > The video card on the machine is nvidia quadro 750 xgl running this driver > > from > > > nvidia's ftp site: > > > > > > NVIDIA-Linux-x86-1.0-5328-pkg1.run > > > > > > I know this is not the latest version of the driver, would upgrading it > > make a > > > difference? > > >

1 0

v1.1917 OpenGL misconfiguration detected.
by Pattanayek, Sabuj 29 Mar '04

29 Mar '04

Hi, Finally got around to installing the latest development linux version. At startup I'm getting the following messages: OpenGL misconfiguration detected. Please reinstall graphics driver. v1.1872 & v1.1892 still work on the machine with no problems. The video card on the machine is nvidia quadro 750 xgl running this driver from nvidia's ftp site: NVIDIA-Linux-x86-1.0-5328-pkg1.run I know this is not the latest version of the driver, would upgrading it make a difference? Thanks, ----------------------------------------------------------------- Sabuj Pattanayek Martin Egli Lab Vanderbilt University Dept. Biochemistry http://linegli.mc.vanderbilt.edu

2 1

export pdb file
by Shuanghong Huo 29 Mar '04

29 Mar '04

Hi, I remember that Midus can export pdb file, but I cannot find this function in Chimera. Thanks, Shuanghong

1 0

question
by Wendy Ochoa 28 Mar '04

28 Mar '04

Could you please tell me how can I define the secondary structure of the protein i am displaying with chimera. I know the secondary structure of my pdb, so I would like to draw a picture of the protein showing the B sheets and alpha helix in different colors. Thansk a lot!!!

2 1

32 bit / 64 bit Linux mixup
by Thomas Goddard 13 Mar '04

13 Mar '04

The 32 bit Linux Chimera 1917 was accidentally replaced with the 64 bit Linux version yesterday (Thursday 11 AM PST, March 11). The 32 bit version has been restored Friday 7 PM PST, March 12. Tom

1 0

Re: [Chimera-users] chimera-1.1917-linux.exe: Exec format error. Wrong Architecture.
by Pattanayek, Sabuj 13 Mar '04

13 Mar '04

[root@eglimp]# uname -iorsp Linux 2.4.22-openmosix2smp athlon i386 GNU/Linux it seems the link to chimera-1.1917-linux.exe on the download page is actually the link to the x86_64 version: [root@eglimp]# head chimera-1.1917-linux.exe ELF> @@@@"@@@@@�@@@@d�d�PP��``P`@@ P�d � � �||Q�/lib64/ld-linux-x86-64.so.2GNU etc.. and [root@eglimp]# head chimera-1.1872-linux.exe ELF44 (#�#��x▒�@ /lib/ld-linux.so.2GNUCVN>LU@SB<5QF/4)1K$.TG!J▒= Can you please put the correct link to the x86_32 linux version on the download page? Thanks, Sabuj Pattanayek Quoting Greg Couch <gregc(a)cgl.ucsf.edu>: > The linux distribution is only for the Intel/AMD x86 architecture. Are > you trying in install on something else? For example, and Alpha or > PowerPC? > > - Greg > > On Fri, 12 Mar 2004, Pattanayek, Sabuj wrote: > > > I'm getting this error when I try to run the installation > > > > -rwx------ 1 root root 26226504 Mar 12 12:02 > chimera-1.1917-linux.exe > > [root@eglimp]# chimera-1.1917-linux.exe > > chimera-1.1917-linux.exe: Exec format error. Wrong Architecture. > > [root@eglimp]# ./chimera-1.1917-linux.exe > > ./chimera-1.1917-linux.exe: Exec format error. Wrong Architecture. > > [root@eglimp]# md5sum chimera-1.1917-linux.exe > > 67e955ce87d000916039d94c7d6d0a56 chimera-1.1917-linux.exe > > [root@eglimp]# uname -rs > > Linux 2.4.22-openmosix2smp > > > > web page: > > > > linux chimera-1.1917-linux.exe > > Size: 26226504 bytes > > MD5: 67e955ce87d000916039d94c7d6d0a56 >

2 1

save session stripping TER statements from PDB's
by Pattanayek, Sabuj 13 Mar '04

13 Mar '04

In version 1.1892 (and older possibly), file > save session seems to be stripping TER statements from PDB files. The TER statement is helpful for specifically telling Chimera and other PDB viewers not to connect consecutive residues in a PDB (for example ATP/water molecules at the end of a PDB that are not part of the main protein chain). Re-inserting TER's manually into the .py file can become cumbersome especially when one PDB insertion in the .py is 27000+ columns long. Would it be possible to fix this in future versions? Thanks in advance, Sabuj Pattanayek

2 1

[Chimera-users] chimera-1.1917-linux.exe: Exec format error. Wrong Architecture.
by Pattanayek, Sabuj 12 Mar '04

12 Mar '04

I'm getting this error when I try to run the installation -rwx------ 1 root root 26226504 Mar 12 12:02 chimera-1.1917-linux.exe [root@eglimp]# chimera-1.1917-linux.exe chimera-1.1917-linux.exe: Exec format error. Wrong Architecture. [root@eglimp]# ./chimera-1.1917-linux.exe ./chimera-1.1917-linux.exe: Exec format error. Wrong Architecture. [root@eglimp]# md5sum chimera-1.1917-linux.exe 67e955ce87d000916039d94c7d6d0a56 chimera-1.1917-linux.exe [root@eglimp]# uname -rs Linux 2.4.22-openmosix2smp web page: linux chimera-1.1917-linux.exe Size: 26226504 bytes MD5: 67e955ce87d000916039d94c7d6d0a56 Thanks in advance Sabuj Pattanayek

2 1