- Chimera-users - RBVI Mailing Lists

center of mass for multiple models.
by Jean Didier Pie Marechal 18 Jul '08

18 Jul '08

Hi everyone, I am stuck with a script and my mistake seems to come from my poor understanding on how to treat atom coordinates with multiple models. I join to this e-mail a chimera file with a simple system to illustrate my problem. In the file, I have two indentique phenyl moleculer each of them in different models and I would like to stock the center of mass of both molecules in a list. I wrote this very simple script using the way of calcuting the center of mass you provided me recently. import chimera import BuildStructure from chimera import runCommand as run from chimera import selection, Vector from chimera import openModels, Molecule, Point #Stock the two fragments in the list FRAG run("select") frag=[] for i in chimera.selection.currentResidues(): frag.append(i) #Calculate center of mass of both entities #stocked in the list MAS mas=[] for i in range(len(frag)): sum=0 center = Vector(0, 0, 0) for a in frag[i].atoms: sum += a.element.mass center += a.element.mass * Vector(*a.coord().data()) mas.append(Point(center[0]/sum,center[1]/sum,center[2]/sum)) print i,mas[i] The results I have are 0 20.0465 70.3735 23.9063 1 20.0465 70.3735 23.9064 Basically, the coordinates of the center of mass are the same? What did I do wrong in the way I treated this problem? Thanks a lot for any answer, All the best JD Dr. Jean-Didier Maréchal Professor Lector Unitat de Química Física Departament de Química Universitat Autònoma de Barcelona Edifici C.n. 08193 Cerdanyola (Barcelona) Tel: +34.935814936 e-mail: JeanDidier.Marechal(a)uab.es ----- Missatge original ----- De: chimera-users-request(a)cgl.ucsf.edu Data: Dijous, Juliol 17, 2008 9:00 pm Assumpte: Chimera-users Digest, Vol 63, Issue 17 > 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: Possible image saving bug (Greg Couch) > 2. Possible image saving bug (Elaine Meng) > > > -------------------------------------------------------------------- > -- > > Message: 1 > Date: Thu, 17 Jul 2008 10:35:26 -0700 (PDT) > From: Greg Couch <gregc(a)cgl.ucsf.edu> > Subject: Re: [Chimera-users] Possible image saving bug > To: Walter Novak <wnovak(a)brandeis.edu> > Cc: Chimera forum <chimera-users(a)cgl.ucsf.edu> > Message-ID: <Pine.OSF.4.63.0807171001010.1939086(a)guanine.cgl.ucsf.edu> > Content-Type: TEXT/PLAIN; charset=US-ASCII; format=flowed > > On Thu, 17 Jul 2008, Walter Novak wrote: > > > Hi All, > > > > I submitted a bug report on this, but the behavior is both odd > and > > reproducible for me. Using the latest production release and > trying to save > > an image using supersampling, the image is not stitched back > together > > appropriately. It is odd, because it always happens for a certain > session > > file I have. It doesn't seem to happen if I open a new molecule > and make a > > picture. I can provide the session file if needed. First image is > using 3x3 > > and second is 6x6. > > I see from the bug report that you're using Mac OS X 10.5. We are > having > "issues" with the new X11 application in 10.5. And we have seen > that > slightly different X11 applications are installed on otherwise > identical > Macs and that reinstalling from the Mac OS X Install Disk fixed a > non-chimera bug. > > Also, I have read on the web that if you install the 10.4 X11 > application > on 10.5, then OpenGL applications will work better. Since that is > not an > option for everyone, we are spending more time to improve the Aqua > version > that is available as a Snapshot release on the chimera download > page. It > has different bugs, but the graphics is better. So you should try > it out > and see how well it works for you. > > Greg Couch > UCSF Computer Graphics Lab > > > ------------------------------ > > Message: 2 > Date: Thu, 17 Jul 2008 10:40:13 -0700 > From: Elaine Meng <meng(a)cgl.ucsf.edu> > Subject: [Chimera-users] Possible image saving bug > To: wnovak(a)brandeis.edu > Cc: Chimera BB <chimera-users(a)cgl.ucsf.edu> > Message-ID: <B15630DD-299C-4B48-A8B4-469FA8995872(a)cgl.ucsf.edu> > Content-Type: text/plain; charset=US-ASCII; format=flowed; delsp=yes > > Hi Wally, > I also had this tiling problem. We think it may be something about > > X11 in Mac OS 10.5 (possibly just 10.5.4), as the problem has not > been > observed with 10.4/X11. > > Two workarounds that were successful in my case (I don't know if > they > will fix all instances of the problem): > > (A) simply show the Command Line before saving the image. I tried > hiding it again after that, but then the problem reoccurred, albeit > in > a more subtle fashion. > > (B) save the session, restore it in the aqua version of Chimera > (still > on Mac OS 10.5.4 but not using X11) and save the image. Yeah, (A) > is > easier, but I didn't discover that one until this morning, after > receiving your mail. > > I hope these solutions are general! Sorry about the inconvenience, > 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 > > > > > ------------------------------ > > _______________________________________________ > Chimera-users mailing list > Chimera-users(a)cgl.ucsf.edu > http://www.cgl.ucsf.edu/mailman/listinfo/chimera-users > > > End of Chimera-users Digest, Vol 63, Issue 17 > ********************************************* >

2 1

center of mass for multiple models
by Jean Didier Pie Marechal 18 Jul '08

18 Jul '08

Hi everyone, I am stuck with a script and my mistake seems to come from my poor understanding on how to treat atom coordinates with multiple models. I join to this e-mail a chimera file with a simple system to illustrate my problem. In the file, I have two indentique phenyl moleculer each of them in different models and I would like to stock the center of mass of both molecules in a list. I wrote this very simple script using the way of calcuting the center of mass you provided me recently. import chimera import BuildStructure from chimera import runCommand as run from chimera import selection, Vector from chimera import openModels, Molecule, Point #Stock the two fragments in the list FRAG run("select") frag=[] for i in chimera.selection.currentResidues(): frag.append(i) #Calculate center of mass of both entities #stocked in the list MAS mas=[] for i in range(len(frag)): sum=0 center = Vector(0, 0, 0) for a in frag[i].atoms: sum += a.element.mass center += a.element.mass * Vector(*a.coord().data()) mas.append(Point(center[0]/sum,center[1]/sum,center[2]/sum)) print i,mas[i] The results I have are 0 20.0465 70.3735 23.9063 1 20.0465 70.3735 23.9064 Basically, the coordinates of the center of mass are the same? What did I do wrong in the way I treated this problem? Thanks a lot for any answer, All the best JD Dr. Jean-Didier Maréchal Professor Lector Unitat de Química Física Departament de Química Universitat Autònoma de Barcelona Edifici C.n. 08193 Cerdanyola (Barcelona) Tel: +34.935814936 e-mail: JeanDidier.Marechal(a)uab.es ----- Missatge original ----- De: chimera-users-request(a)cgl.ucsf.edu Data: Dijous, Juliol 17, 2008 9:00 pm Assumpte: Chimera-users Digest, Vol 63, Issue 17 > 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: Possible image saving bug (Greg Couch) > 2. Possible image saving bug (Elaine Meng) > > > -------------------------------------------------------------------- > -- > > Message: 1 > Date: Thu, 17 Jul 2008 10:35:26 -0700 (PDT) > From: Greg Couch <gregc(a)cgl.ucsf.edu> > Subject: Re: [Chimera-users] Possible image saving bug > To: Walter Novak <wnovak(a)brandeis.edu> > Cc: Chimera forum <chimera-users(a)cgl.ucsf.edu> > Message-ID: <Pine.OSF.4.63.0807171001010.1939086(a)guanine.cgl.ucsf.edu> > Content-Type: TEXT/PLAIN; charset=US-ASCII; format=flowed > > On Thu, 17 Jul 2008, Walter Novak wrote: > > > Hi All, > > > > I submitted a bug report on this, but the behavior is both odd > and > > reproducible for me. Using the latest production release and > trying to save > > an image using supersampling, the image is not stitched back > together > > appropriately. It is odd, because it always happens for a certain > session > > file I have. It doesn't seem to happen if I open a new molecule > and make a > > picture. I can provide the session file if needed. First image is > using 3x3 > > and second is 6x6. > > I see from the bug report that you're using Mac OS X 10.5. We are > having > "issues" with the new X11 application in 10.5. And we have seen > that > slightly different X11 applications are installed on otherwise > identical > Macs and that reinstalling from the Mac OS X Install Disk fixed a > non-chimera bug. > > Also, I have read on the web that if you install the 10.4 X11 > application > on 10.5, then OpenGL applications will work better. Since that is > not an > option for everyone, we are spending more time to improve the Aqua > version > that is available as a Snapshot release on the chimera download > page. It > has different bugs, but the graphics is better. So you should try > it out > and see how well it works for you. > > Greg Couch > UCSF Computer Graphics Lab > > > ------------------------------ > > Message: 2 > Date: Thu, 17 Jul 2008 10:40:13 -0700 > From: Elaine Meng <meng(a)cgl.ucsf.edu> > Subject: [Chimera-users] Possible image saving bug > To: wnovak(a)brandeis.edu > Cc: Chimera BB <chimera-users(a)cgl.ucsf.edu> > Message-ID: <B15630DD-299C-4B48-A8B4-469FA8995872(a)cgl.ucsf.edu> > Content-Type: text/plain; charset=US-ASCII; format=flowed; delsp=yes > > Hi Wally, > I also had this tiling problem. We think it may be something about > > X11 in Mac OS 10.5 (possibly just 10.5.4), as the problem has not > been > observed with 10.4/X11. > > Two workarounds that were successful in my case (I don't know if > they > will fix all instances of the problem): > > (A) simply show the Command Line before saving the image. I tried > hiding it again after that, but then the problem reoccurred, albeit > in > a more subtle fashion. > > (B) save the session, restore it in the aqua version of Chimera > (still > on Mac OS 10.5.4 but not using X11) and save the image. Yeah, (A) > is > easier, but I didn't discover that one until this morning, after > receiving your mail. > > I hope these solutions are general! Sorry about the inconvenience, > 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 > > > > > ------------------------------ > > _______________________________________________ > Chimera-users mailing list > Chimera-users(a)cgl.ucsf.edu > http://www.cgl.ucsf.edu/mailman/listinfo/chimera-users > > > End of Chimera-users Digest, Vol 63, Issue 17 > ********************************************* >

1 0

Possible image saving bug
by Elaine Meng 17 Jul '08

17 Jul '08

Hi Wally, I also had this tiling problem. We think it may be something about X11 in Mac OS 10.5 (possibly just 10.5.4), as the problem has not been observed with 10.4/X11. Two workarounds that were successful in my case (I don't know if they will fix all instances of the problem): (A) simply show the Command Line before saving the image. I tried hiding it again after that, but then the problem reoccurred, albeit in a more subtle fashion. (B) save the session, restore it in the aqua version of Chimera (still on Mac OS 10.5.4 but not using X11) and save the image. Yeah, (A) is easier, but I didn't discover that one until this morning, after receiving your mail. I hope these solutions are general! Sorry about the inconvenience, 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

Possible image saving bug
by Walter Novak 17 Jul '08

17 Jul '08

Hi All, I submitted a bug report on this, but the behavior is both odd and reproducible for me. Using the latest production release and trying to save an image using supersampling, the image is not stitched back together appropriately. It is odd, because it always happens for a certain session file I have. It doesn't seem to happen if I open a new molecule and make a picture. I can provide the session file if needed. First image is using 3x3 and second is 6x6. Also, when using pov-ray, the balls are completely devoid of shininess, regardless of my shininess settings (I set it to 128). That is the last image. Thanks, Wally Walter R.P. Novak, Ph.D. Postdoctoral Fellow Rosenstiel Basic Medical Research Center Brandeis University 415 South St. MS 029 Waltham, MA 02454-9110 Phone: (781) 736-4944 Fax: (781) 736-2405

2 1

Coloring a binding site
by Tom Goddard 16 Jul '08

16 Jul '08

Hi Andy, To color a binding site I would select the ligand (DNA in your case) then select all the nearby atoms using Select / Zone..., then use Actions / Color to color those atoms. Tom Anindito Sen wrote: > Another question - I need to color a portion of the crystal structure > differently (ribbon diagram) in order to indicate the binding sites of > the DNA with the protein molecule and how do I do that. > > Thanks and Best Wishes > Andy

2 2

possible bug
by Ben Keshet 16 Jul '08

16 Jul '08

Hi, when I un-select atoms, I had some weird long lines crossing my structures across (picture attached). I restarted the program how ever the same exact thing happened again. The two files that were opened worked perfectly in the past and before I un-selected the atom (ctrl+left click in the "air"). Other wise - great software. Thank you. Ben

2 1

Re: [Chimera-users] History feature for command line between sessions
by Eric Pettersen 15 Jul '08

15 Jul '08

> On Apr 17, 2008, at 8:52 AM, Boaz Shaanan wrote: > > > The history of commands in the command line mode is retained > > within a session and is very convenient but I was wondering whether > > there is a feature that will allow you to recall commands from > > previous sessions (such as in O for example - it generates an .o- > > history file which is being constantly updated and is very > > convenient). If not, is it possible to request the addition of such > > a feature ? There are many ways to implement this of course and I'm > > sure you'll find the best way, if it's not already there lurking in > > Chimera. > > Hi Boaz, > Good idea -- and shouldn't be too hard to implement. I was thinking > of remembering the last 100 commands. Will that be enough you think? > > --Eric This is implemented now. It should show up in the next daily build that works. Only remembers 60 by default, but that number is changeable in Preferences. --Eric Eric Pettersen UCSF Computer Graphics Lab http://www.cgl.ucsf.edu

1 0

Can't type in Chimera
by Thomas Goddard 15 Jul '08

15 Jul '08

Hi Andy, I don't know how Chimera could get in a state where it will not accept any typed input. Do you recall the last thing you did in Chimera before it stopped accepting key input? If I were stuck in this situation and could not type a file name to save the Chimera session I would create an empty file with the session file name I wanted using a text editor. Then choose that file with the mouse in the Chimera session saving dialog. Desperate measures. If you can reproduce the problem please tell me how and I will fix Chimera. Tom Anindito Sen wrote: > HI Tom > > I have installed the latest version of Chimera on my LINUX box Fedora > Core 6.0. I faced a strange problem- while carrying out a docking of a > crystal structure in a density map I was using the command line for the > differnt movements and controls. After some time I am unable to type > anything on the command line. Infact I am even unable to save the > session as I cannot type the "NAME" of the session. I checked my keyboad > and typed on the terminal and on text editiors and it worked normally. > > Waiting for your response. > > Best Wishes > Andy > > > Dr. Anindito Sen (Ph.D) > Research Associate , Dept. of Biochemistry and Molecular Genetics > University of Virginia > Box 800733 Charlottesville, VA 22908

1 0

Re: [Chimera-users] reg. ribbon reprsentation
by Elaine Meng 15 Jul '08

15 Jul '08

>> >> Greetings. I have already posted a similar question on the forum >> previously and hence i am writing to your email. My problem is when >> i open a protein pdb, i get wonderful ribbon representation when i >> opt for it. But when i give a protein pdb with only a CA atoms/ >> CA,C,N atoms, i dnt get the ribbon representation. I understand >> that ribbon is constructed by considering CA,C,N,O and OXT atoms. >> What i did was that i simply removed O and OXT atoms by going >> through Tools--> Depiction --> Ribbon style Editor. It may be a >> wrong idea to do that. But i would appreciate if you could write me >> some way to draw the ribbon for protein having only CA or CA,C,N >> atoms IF POSSIBLE. >> >> Regards, >> Bala >> Hi Bala, If you have N-CA-C atoms, it actually works pretty well if you open Ribbon Style Editor, go to the "Residue Class" tab, choose "amino acid" but then change the "Orientation atom" to C instead of O, and Apply. (You still need to show ribbon if you didn't already do that.) You can save and name the new residue class for re-use in some later Chimera session if you want. The resulting ribbon looks the same to me as the standard ribbon made from N,CA,C,O (I tested this by writing a PDB of only the N,CA,C atoms from a complete structure, reading it back in, and experimenting with Ribbon Style Editor). http://www.cgl.ucsf.edu/chimera/docs/ContributedSoftware/ribbonstyle/ribbon… If you have CA atoms only, then there is no way to figure out the plane of the ribbon - it does not work to set both the "Guide atom" and "Orientation atom" to CA. Unfortunately this prevents you from displaying "licorice" ribbons even though they don't use an orientation plane. My only solution in this case is to use some program to build the additional atoms from the CA trace, for example this web server: http://www.ebi.ac.uk/Tools/maxsprout/index.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

1 0

Re: [Chimera-users] email address
by Khanei Mkumeni 14 Jul '08

14 Jul '08

khmlukundo(a)gmail.com

1 0

SpaceNavigator with Linux
by D. Greening 14 Jul '08

14 Jul '08

Hi, I just downloaded Chimera 1.2540 and was pleased to see that the SpaceNavigator 3-D mouse is now supported for Windows and Mac. Will those of us who use Linux be able to use our SpaceNavigators with Chimera any time soon? Thanks. D. Greening Memorial University of Newfoundland

2 1

drawing molecules
by Kathleen.Zocher＠ernaehrung.uni-giessen.de 14 Jul '08

14 Jul '08

Dear Sir or Madam, I`ve just one question concerning some special tools of Chimera. Is it possible to draw small molecules (ligands) within chimera? I can not find some information for this problem in the manual. For further docking experiments I just need a 3D-structure which is minimized (in a mol2 file). Many thanks in advance! Sincerely yours Kathleen Zocher ---------------------------------------------------------------- This message was sent using IMP, the Internet Messaging Program.

2 1

Question about bond length modification
by Benjamin G.S. Horsman 14 Jul '08

14 Jul '08

Hello, I have what appears to be simple problem but is I can't find the solution in previous mailing list queries. My problem is that I am unable to modify the lengths of any bonds in a PDB molecule. I am rather frustrated with myself over such a stupid problem, having tried dozens of times over dozens molecules to get it right but to no avail. Any possible solutions from the community? I am using MacpyMol to generate PDB structures and Chimera v1 build 1470. Thanks for any help, Ben

2 1

Re: [Chimera-users] [Chimera-announce] Chimera production release
by Elaine Meng 13 Jul '08

13 Jul '08

Hello Yao, Don't worry, the new versions still include the ability to fit map in map. It is simply that two tools were combined into a new tool with a slightly different name. Changes in this production release relative to the previous one include: "Fit Model in Map" and "Fit Map in Map" (Volume Data) — * merged into new "Fit in Map" tool * can update fit values continuously * can simulate map from atomic coordinates, report map-map correlation current documentation: http://www.cgl.ucsf.edu/chimera/current/docs/ContributedSoftware/ fitmaps/fitmaps.html Best, 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 Jul 12, 2008, at 11:49 PM, Yao Cong wrote: > Hello Chimera developer, > > Appreciate your great efforts to make Chimera an perfect tool for > Cryo-EM study. > > I and many of my colleges notice that in your latest versions, you > omit the "fit map into map" module, which is very convenient in > many cases. Is there a specific reason for that? And we all > consider it will be very helpful if you could add this function > back in your next version, maybe daily release? Otherwise, we still > keep the old version to have that function and meanwhile install > your newer version to enjoy the newer functions. > > Thanks, > Yao > ----------------------- > Yao Cong, Ph.D. >

1 0

Creating abridged multimer
by Francesco Pietra 12 Jul '08

12 Jul '08

I am posting again after having once again subscribed In order to create an abridged multimer made of non covalently bound subunits (from available pdb of subunit) which is the preferable route? (1) Load the monomer, delete substructures, and then assemble an abridged multimer? Are the headers furnished with the monomer OK for the abridged subunit? (2) Build the unabridged multimer and proceed to simplify each subunit the same way? Thanks francesco pietra

2 1

distance measurements on surfaces
by Daniel Hautzinger 10 Jul '08

10 Jul '08

Dear all, from the Volume Visualization Tour, I discovered the topography command that now allows me to dsiplay atomic force microscopy data (basically bitmap images with grey scales representing height). In that way I got some very nice images. I also found out how to scale the voxel size, so the scale of the structures should be correct. To test this, and to be able to add some length markers, I wanted to measure distances on the surface generated by the topography command. But while that seems possible for volume data (by setting markers), I had no success with the surface representation of my data nor the original image (which is handled as kind of a volume, but a single slice) when it came to setting markers or measuring distances in another way. Did I miss anything, or is that just not possible? I am using build number 2509, because of the topography feature. Setting markers is described as working by setting the marker in the first maximum in the line of sight at the click position. For surfaces, you could just take the position of the visible pixel clicked. Or maybe it is possible to convert surfaces to volumes? Thank you for your help Daniel Hautzinger

2 1

changing selection color
by Michael DiMattia 10 Jul '08

10 Jul '08

Hi all, I have a quick question regarding a figure I am trying to make. I have a multi-chain protein loaded with each subunit surface-rendered individually. We wanted to find a way to highlight the outline of one of the subunits in the 3d surface rendering and selecting the surface of the chain works really well for this in Chimera. The familiar green outline highlights the subunit well, even though I realize this is probably not the primary purpose for such highlighting. I imagine it is just to let the user know what is selected. But if i were to use the selection highlighting to show such an outline in a figure, i was wondering if it is possible to change the color of the selection outline from green to any other color or if it is possible to change the width of the outline. I realize it's a long shot, but thought I'd ask anyway in case anyone has tried this before. Many thanks Mike DiMattia

2 1

Chimera production release
by Elaine Meng 10 Jul '08

10 Jul '08

Hi everybody, A new production release, version 1.2540, is now available from: http://www.cgl.ucsf.edu/chimera/download.html#production Many improvements and new features have been added since the previous production release (Nov 2007), as detailed in the release notes: http://www.cgl.ucsf.edu/chimera/docs/relnotes/1.2540.html For the time-challenged, here are just a few of the reasons you may want to use the new version: 1) new fetch types: - comparative models from ModBase - protein pocket data from CASTp database - small molecule structures from Pub3D - small molecule structures built from SMILES strings 2) surface improvements: - increased reliability - surface areas automatically calculated - Actions menu, other tools work on more types of surfaces 3) appearance: - default lighting changed to improve POV-Ray results - larger default scale - helix ribbon insides can be colored separately - 2D labels can include Greek letters and other symbols 4) new tools: Axes, EnsembleCluster, Solvate, Structure Diagram 5) new commands: export, mclip, molmap, perframe, ribinsidecolor, shape, solvate, topography, vop See the release notes for details. The production release includes versions for Windows, Mac/X11, Linux, SGI IRIX, and HP Tru64. The Mac Aqua version has been improved significantly, but for now remains a snapshot download: http://www.cgl.ucsf.edu/chimera/download.html#snapshots Enjoy! The Chimera Team

1 0

Multiscale models and surface coloring
by David Sept 09 Jul '08

09 Jul '08

I need to create a series of images of large protein complexes, and since it is tedious to do this through the menus, I am trying to script this for Chimera. I want to create low-res, EM looking surfaces and color the individual protein chains different colors. So far, I can read in multiple structures and create the surfaces, but I can't figure out how to change the coloring. Here's what I have: ---- test.py ------- from chimera import openModels m = openModels.open('model-1.pdb')[0] import MultiScale d = MultiScale.show_multiscale_model_dialog() d.make_multimers([m]) chimera.viewer.viewAll() ---- test.py ------- The above python script work and makes very nice blobs. I can color the entire structure (say red) by adding something like: d.select_all_chains_cb() d.change_color_cb((1, 0, 0, 1)) but I can't figure out how to select individual chains and color them (not sure how to use select_chains). This also gets more complicated when I load in two structures at the same time and they have common chain letters. Since I found a way to get this (almost) working in python, I've been sticking with it, but if this can easily be done directly using a chimera command file, I'm happy to do that. Any help is appreciated. Thanks, Dave ---------------------------------------------------- David Sept, Ph.D. Associate Professor Biomedical Engineering and Center for Computational Biology Washington University, Campus Box 1097, St. Louis, MO 63130 Phone: (314) 935-8837 Fax: (314) 935-7448 dsept(a)biomed.wustl.edu http://habanero.ibme.wustl.edu

3 3

Coloring different set of marker with different radius
by Bui Khanh Huy 08 Jul '08

08 Jul '08

Hi, I tried to color my EM electron density and thought that if I can color different maker sets with different radii, that would make thing much easierMy idea is to color one set using radius a, then select another set, color with radius b. But then the color from the first set disappeared. . Is there any way to do that? Best regards, Huy

2 1

Mac Chimera in-window menu bars
by Thomas Goddard 07 Jul '08

07 Jul '08

Hi Kevin, The Chimera Aqua version has a menu bar at the top of the screen as normal Mac applications do and also at the top of the graphics window. I agree there should be a preference setting to eliminate the in-window menu bar. Here's the reason for this weirdness. There are 13 Chimera dialogs that have their own menu bars and menus. Use of these involves switching between actions in the main window and the dialog repeatedly. With menus for the active window only shown at the top of the screen this interaction becomes very tedious. For this reason I added menu bars within each Chimera window that has its own set of menus. Those in-window menus don't work very well. If you post one menu then put the pointer over an adjacent menu it will not pop up. An even worse atrocity in Aqua Chimera is that it implements "focus follows mouse". As the pointer is moved over the various Chimera windows the one underneath becomes active. This again was added because common Chimera usage involves switching frequently between windows and the normal Mac "click to focus" behavior is inconvenient for this. We would like instead to use "click through" where the first mouse click on a window both activates that window and presses the button or activates the entry field that was clicked on. I have not been able to do this in Aqua Tk, the windowing toolkit we use. The problem we're facing is that the Chimera user interface with its many windows and frequent switching between windows is poorly suited to the Mac GUI conventions. The in-window menus and focus-follows-mouse mode are a first stab at a compromise solution -- still much to be desired. Tom Kevin Walzer wrote: The latest build of Chimera from 7/3/08 on OS X/Aqua (10.5.4) includes menubars at the top of the screen (the standard OS X behavior) *and* on the main toplevel window itself, as on Windows and X11. This behavior is simply weird--two menus? It can only b e regarded as a bug, not a feature. What purpose does it serve? The correct Mac behavior can be restored by commenting out lines 1470-1476 of tkgui.py. If some of your Mac users are requesting this kind of behavior in menubars, at the very least the non- standard window-based menubar should be a preferences item that can be toggled off.

1 0

Tile Structure default view
by Eduardo Sanz Garcia 07 Jul '08

07 Jul '08

I would like to request an extra feature for chimera, if it is not there already. I normally open multiple structures in chimera and then I use the "Tile Structure" button to visualize them in an array. This is nice, for example, to visualize the evolution of a cryo-EM 3D reconstruction. I would like that chimera could open multiple structures in an array view by default. Could a "Default view" button be added to the Tile Structure window, so that the "Border scale" value setting is saved and then the next time I open several structures they would be displayed in array view? Thank you very much.

3 2

renumbering pdb
by bustik＠gmx.net 05 Jul '08

05 Jul '08

dear chimera-users, i am just starting to use chimera, so this will be a really simple question to you (although i didn't find anything about this in the manual...): is there a simple command to change residue numbers in a pdb file (i built up a homology model and want to elevate all residue numbers in the model the same amount)? thanks in advance sebastian -- Sebastian Kruggel Eppendorfer Weg 17d D 20259 Hamburg fon: +49 (0)40 4192 9579 mobil: +49 (0)151 1789 2232 e-mail: bustik(a)gmx.net

2 2

graphics performance: amira v. chimera
by Matthew Dougherty 02 Jul '08

02 Jul '08

I loaded a MRC file, isosurfaced it, gave it a transparency around 0.48, and generated a vrml in chimera. I loaded the vrml into amira. initial comparison: I reloaded the vrml into chimera, I turned off the fog on chimera, moving the object around the refresh rate seems to be around 8 fps. On amira it seems to be around 24fps. results: The graphics on Chimera seems to be better during the rotation of transparent objects: full rotation in chimera is flawless, full rotation in amira results in flipping of the object (vertices in back quantum-ly jump to the front, front vertices go to back, looks like a z buffer error) periodically at 180 degrees, making it useless for many camera positions. Any comments or suggestions?

4 3

OS X Chimera
by Christopher Vavricka 02 Jul '08

02 Jul '08

I have OS X 10.4.11 and I am trying to run Chimera. When I start Chimera.app it tells me to start X windowing system (Apple X server). X11 is not located in my applications/utlilties folder. How do I get this application? Is it on my installation disc for OS X? Regards, Chris Vavricka

3 2

(no subject)
by Levy, Hazel C. 02 Jul '08

02 Jul '08

Hi, I am a postdoctoral researcher and an enthusiastic Chimera user. We are trying to finish up a manuscript to submit tomorrow to JSB and we have been frantically trying to figure out how to convert from Chimera's "contour level" to the more conventional "sigma level" for density map viewing at different noise levels. We understand that the higher the contour level the smaller the map volume displayed (less noisy), but we are not quit sure how to report this as a sigma level. Each map that we are comparing has a different contour range assigned, I am wondering if sigma level corresponds to some proportion of the maximum volume. Please let us know how to do this conversion or please point us in the direction toward the person that can help us. Thanks you so much...this information will settle a long-standing mystery for us. Hazel Levy

2 1

Mutate by D- res...
by Yasser Almeida Hernández 30 Jun '08

30 Jun '08

Can Chimera mutate by D-residues? If this is possible, how i can do it? Thanks... -- Yasser Almeida Hernández Estudiante, 4to Año Bioquímica Facultad de Biología Universidad de la Habana

3 2

Dummy atoms
by Jean Didier Pie Marechal 30 Jun '08

30 Jun '08

Dear all, I think the question was already posted but I can't find it. I would like to make an atom I placed as an helium to a Dummy: What is the corresponding atom type in Chimera? All the best JD Dr. Jean-Didier Maréchal Professor Lector Unitat de Química Física Departament de Química Universitat Autònoma de Barcelona Edifici C.n. 08193 Cerdanyola (Barcelona) Tel: +34.935814936 e-mail: JeanDidier.Marechal(a)uab.es ----- Missatge original ----- De: chimera-users-request(a)cgl.ucsf.edu Data: Dimarts, Juny 24, 2008 9:00 pm Assumpte: Chimera-users Digest, Vol 62, Issue 21 > 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: how could u know the order of atoms? (Eric Pettersen) > > > -------------------------------------------------------------------- > -- > > Message: 1 > Date: Mon, 23 Jun 2008 14:22:31 -0700 > From: Eric Pettersen <pett(a)cgl.ucsf.edu> > Subject: Re: [Chimera-users] how could u know the order of atoms? > To: Kamal Al Nasr <k_alnasr(a)yahoo.com> > Cc: chimera-users(a)cgl.ucsf.edu > Message-ID: <8901DE4E-D117-4B18-968B-D3095AD735D4(a)cgl.ucsf.edu> > Content-Type: text/plain; charset="us-ascii" > > On Jun 19, 2008, at 10:08 AM, Kamal Al Nasr wrote: > > > Dear Sir/Madam, > > > > I am using Chimera viewer..and it is a great viewer ever... > > Thanks! I'm glad you like it. > > > I kindly have a question... how could you know (from PDB file) > the > > coorect bonding between atoms? I mean that how could you know > that > > this atom is connected with CA not N nor C.....do u have all > these > > bonds stored already in a database? or you use atom positions to > > determine that? > > For standard amino and nucleic acids, Chimera has templates > containing the correct connectivity. For non-standard residues it > bases the connectivity on inter-atomic distances, with a few added > heuristics (e.g. hydrogens can only have one bond). > > The templates are accessible in the Python layer with the > chimera.restmplFindResidue() function. You can look in the > Rotamers > extension for examples of its use. The distance-based connectivity > > function is also available for an entire molecule as > chimera.connectMolecule(). It takes the molecule as its only > argument and adds bonds to it. > > --Eric > > Eric Pettersen > UCSF Computer Graphics Lab > http://www.cgl.ucsf.edu > >

2 1

"Floating exception" error message
by Dong-Yan Tan 27 Jun '08

27 Jun '08

Dear all, We have just installed chimera-1.2470-linux.exe in our Linux computer (SUSE LINUX 9.0 ). We keep getting error message "Floating exeptions" when we are just trying to open a PDB file and the program quits automatically. Has anyone encountered such a problem or have clues about where to look for solutions? Thanks a lot. Dong-Yan -- Dong-Yan Tan Laboratory of Dr. Hernando Sosa Dept. of Physiology and Biophysics Albert Einstein College of Medicine Bronx, NY 10461 Tel: 718-430-3451 Email: dtan(a)aecom.yu.edu -- Dong-Yan Tan Laboratory of Dr. Hernando Sosa Dept. of Physiology and Biophysics Albert Einstein College of Medicine Bronx, NY 10461 Tel: 718-430-3451 Email: dtan(a)aecom.yu.edu

3 3

"Floating exception" error message
by Dong-Yan Tan 27 Jun '08

27 Jun '08

Dear all, We have just installed chimera-1.2470-linux.exe in our Linux computer (SUSE LINUX 9.0 ). We keep getting error message "Floating exeptions" when we are just trying to open a PDB file and the program quits automatically. Has anyone encountered such a problem or have clues about where to look for solutions? Thanks a lot. Dong-Yan -- Dong-Yan Tan Laboratory of Dr. Hernando Sosa Dept. of Physiology and Biophysics Albert Einstein College of Medicine Bronx, NY 10461 Tel: 718-430-3451 Email: dtan(a)aecom.yu.edu -- Dong-Yan Tan Laboratory of Dr. Hernando Sosa Dept. of Physiology and Biophysics Albert Einstein College of Medicine Bronx, NY 10461 Tel: 718-430-3451 Email: dtan(a)aecom.yu.edu

2 1

(no subject)
by Dong-Yan Tan 27 Jun '08

27 Jun '08

Dear all, We have just installed chimera-1.2470-linux.exe in our Linux computer (SUSE LINUX 9.0 ). We keep getting error message "Floating exeptions" when we are just trying to open a PDB file and the program quits automatically. Has anyone encountered such a problem or have clues about where to look for solutions? Thanks a lot. Dong-Yan -- Dong-Yan Tan Laboratory of Dr. Hernando Sosa Dept. of Physiology and Biophysics Albert Einstein College of Medicine Bronx, NY 10461 Tel: 718-430-3451 Email: dtan(a)aecom.yu.edu

1 0

hydrophobic pocket
by bala 27 Jun '08

27 Jun '08

Dear Chimerian's, How do i detect hydrophobic pockets in a protein using chimera ?. Thanks in advance, Bala

3 2

Re: [Chimera-users] Question about avi files in chimera
by Tom Goddard 27 Jun '08

27 Jun '08

Hi Jaap, The AVI files produced by Chimera movie recorder use the MSMPEG-4v2 (Microsoft MPEG4 version 2) video codec and the encoding is done by the open source ffmpeg program. The file size depends on the movie recorder "bit rate" setting (default 2000 Kbits/sec = 250 Kbytes/sec = 15 Mbytes/min). Changing this setting allows you to trade off file size versus image quality. http://www.cgl.ucsf.edu/chimera/docs/ContributedSoftware/recorder/recorder.… Tom Jaap Brink wrote: > > Hi Tom, > > > > I was wondering given the size difference in avi files that chimera > generates versus other program (Video Edit Magic for instance), which > video codec is being used? I searched the list of embedded software > packages in chimera and found the library of libavcodecs. > > > > Thanks in advance, > > > > Jaap >

1 0

upcoming production release
by Greg Couch 26 Jun '08

26 Jun '08

Hello, We are currently getting ready to put out another chimera production release. At this point, we consider the daily builds to be release candidates, so please download one and try it out. Not all of the daily builds will be production releases, in particular, the Mac OS X Aqua version is not quite ready, and will be a snapshot release. And the Linux off-screen OpenGL version is experimental and will probably stay as a daily build. Greg Couch UCSF Computer Graphics Lab Daily builds: <http://www.cgl.ucsf.edu/chimera/alpha-downloads.html> New features: <http://www.cgl.ucsf.edu/chimera/docs/relnotes/snapshot.html>

1 0

Re: [Chimera-users] how to draw lines in Chimera
by Elaine Meng 24 Jun '08

24 Jun '08

On Jun 24, 2008, at 4:24 PM, Qing Xie wrote: > Hi, > I'm trying to draw a bunch of axes from center (0,0,0). I could do > this > in 'O' by calling this macro: > > begin_object symax > colour blue > stick_xyz 139.05, -0.00, -16.25, -139.05, 0.00, 16.25, 0.25 > stick_xyz -111.40, 43.12, 73.02, 111.40, -43.12, -73.02, 0.25 > stick_xyz -107.48, 80.74, 39.10, 107.48, -80.74, -39.10, 0.25 > stick_xyz -66.72, 119.11, -31.00, 66.72, -119.11, 31.00, 0.25 > stick_xyz 50.41, 73.78, -107.77, -50.41, -73.78, 107.77, 0.25 > stick_xyz 100.65, 0.20, -97.31, -100.65, -0.20, 97.31, 0.25 > stick_xyz -68.88, 113.16, 45.26, 68.88, -113.16, -45.26, 0.25 > stick_xyz -39.02, 80.98, -107.34, 39.02, -80.98, 107.34, 0.25 > stick_xyz 10.15, 0.20, -139.64, -10.15, -0.20, 139.64, 0.25 > stick_xyz 50.47, -73.42, -108.00, -50.47, 73.42, 108.00, 0.25 > end_object symax > > Since I had to use some features in Chimera, do you know what > commands I > should combine to do similar things above? > Thank you! > Sincerely, > Qing Hi Qing, You can make a BILD file defining the objects, for example containing the following: .color medium blue .cylinder 139.05 -0.00 -16.25 -139.05 0.00 16.25 1 .cylinder -111.40 43.12 73.02 111.40 -43.12 -73.02 1 .cylinder -107.48 80.74 39.10 107.48 -80.74 -39.10 1 .cylinder -66.72 119.11 -31.00 66.72 -119.11 31.00 1 .cylinder 50.41 73.78 -107.77 -50.41 -73.78 107.77 1 .cylinder 100.65 0.20 -97.31 -100.65 -0.20 97.31 1 .cylinder -68.88 113.16 45.26 68.88 -113.16 -45.26 1 .cylinder -39.02 80.98 -107.34 39.02 -80.98 107.34 1 .cylinder 10.15 0.20 -139.64 -10.15 -0.20 139.64 1 .cylinder 50.47 -73.42 -108.00 -50.47 73.42 108.00 1 and open it in Chimera to create the objects. I also attached this below as the file sticks.bild to be opened directly in Chimera. I didn't think 0.25 looked that good, but you can edit the file to whatever radius (and color, etc.) that you want. The BILD language also includes many other objects: http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/bild.html Best, 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

Re: [Chimera-users] find levels and morph conformations
by Thomas Goddard 24 Jun '08

24 Jun '08

Hi Ed, The Morph Conformations tool works with submodels. It will not work between two models or submodels that have different numbers of chains. The error message you get means that you are giving it models with different numbers of chains. Tom Edward Brignole wrote: > Tom, > > ... I'd like to use morph conformations between two pdbs each > composed of the same submodels/chains. I've tried to do this through the > dialog but it displays an error telling me "cannot interpolate the > models have different numbers of chains". The only way I can see to do > this is to write the pdb as a mol2 so that all atoms are in a single > model without submodels, open the mol2, and run morph conformations. > Is there another method or is this the simplest? > > Thanks, > Ed

3 2

selecting/displaying zones
by Elaine Meng 24 Jun '08

24 Jun '08

On Jun 24, 2008, at 11:56 AM, Thomas C. Marlovits wrote: > > Am using Chimera and would like to do the following: > - would first like to focus on specific side chain residues (for > example lysines) > - at a specific atom of either one or all lysines present in my > sample, I would like to draw a region of interest (which should be > a sphere with a given distance ... then make it rather transparent > like soap bobbles ...) ... the idea of this is to highlight areas > within a given distance away from this atom > - then I would like to read out all other side chains that are > within such an area ( in a table or so) ... and eventually > highlight them with specific colors ... > > Elaine, is there an easy way to do this ? > Tx, > -Thomas > ---------------------------------------------------------------------- > --------- > > Thomas C. Marlovits, PhD, M.A.S. > IMP-IMBA Research Center > - Research Institute of Molecular Pathology (IMP) > - Institute of Molecular Biotechnology, Austrian Academy of > Sciences (IMBA) Hi Thomas, You can definitely (1) select the atoms or residues within a specified cutoff of some other atoms and write a list of what is selected to the Reply Log or to a file. Slightly trickier is (2) displaying a "bubble." (1) The following command selects all atoms within 4.5 angstroms of atom NZ in lysine 92, chain A (example structure 2gbp): select :92.a@nz za<4.5 to select the whole residues with any atom in the zone, use "z" or "zr" instead of "za", and to use all the atoms of that lysine, just omit "@nz" from the command. http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/ atom_spec.html#zones This selection process can also be done by Ctrl-clicking to select the atom and then using Select... Zone. You can write out a list of the selection with Actions... Write List (click Log if you want to send it to the Reply Log instead). The default is a list of all the atoms individually, but you can set it to "residues" instead to just get one line of output per residue. The command to do that same thing is "writesel" http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/writesel.html (2) The "bubble." I imagine something like Volume Eraser would be nice (check it out: Tools... Volume Data... Volume Eraser). However, that does not let you specify particular coordinates for the center. You have to move it around by hand. Volume Eraser really just makes a big atom and gives it a transparent color. You could add an atom in Build Structure (under Tools... Structure Editing), color it a transparent color, and change its radius with "vdwdefine" http://www.cgl.ucsf.edu/chimera/docs/ContributedSoftware/editing/ editing.html#build http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/vdwdefine.html It doesn't let you enter specific coordinates for the atom. However, if you add the atom as a "new model" you can match it to different locations on the first. For example, if the new atom is model 1 and the structure of interest (say 2gbp) is model 0: rep sphere #1 colordef tpink 1 .6 .6 .6 color tpink #1 vdwdef 3.0 #1 vdwdef +1 #1 vdwdef +1 #1 match #1 #0:147.a@nz select #1 za<5 & #0 [... etc. ...] The radius would just be for display purposes. You'd still have to specify the zone distance when selecting. A couple of failed ideas: (a) You can make a spherical surface of specified radius with the command "shape" http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/shape.html but you cannot directly specify coordinates for the center or easily modify radius after it has been created. (b) Create a spherical object in BILD format. http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/bild.html Disadvantages are inability to make it transparent and having to delete/recreate the object just to change its radius. Best, 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

how could u know the order of atoms?
by Kamal Al Nasr 23 Jun '08

23 Jun '08

Dear Sir/Madam, I am using Chimera viewer..and it is a great viewer ever... I kindly have a question... how could you know (from PDB file) the coorect bonding between atoms? I mean that how could you know that this atom is connected with CA not N nor C.....do u have all these bonds stored already in a database? or you use atom positions to determine that? I appreciate your time.... Kamal

2 1

connectivity
by Jean Didier Pie Marechal 18 Jun '08

18 Jun '08

Dear all, I would like to use chimera to generate inputs for gaussian QM/MM calculations (ONIOM). One of the element I need is the connectivity of the atoms. I looked around the programmer's guide but couldn't find the right module that takes this in charge. what is the way to have access the connectivity of a given molecule? All the best, JD Dr. Jean-Didier Maréchal Professor Lector Unitat de Química Física Departament de Química Universitat Autònoma de Barcelona Edifici C.n. 08193 Cerdanyola (Barcelona) Tel: +34.935814936 e-mail: JeanDidier.Marechal(a)uab.es ----- Missatge original ----- De: chimera-users-request(a)cgl.ucsf.edu Data: Dimecres, Juny 18, 2008 9:00 pm Assumpte: Chimera-users Digest, Vol 62, Issue 18 > 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. Find contour level to enclose a specified volume (Thomas > Goddard) 2. Re: Find contour level to enclose a specified volume > (Tom Goddard) > > > -------------------------------------------------------------------- > -- > > Message: 1 > Date: Tue, 17 Jun 2008 15:44:57 -0700 > From: Thomas Goddard <goddard(a)cgl.ucsf.edu> > Subject: [Chimera-users] Find contour level to enclose a specified > volume > To: "'Chimera BB'" <chimera-users(a)cgl.ucsf.edu> > Message-ID: <48583E69.5090608(a)cgl.ucsf.edu> > Content-Type: text/plain; charset="iso-8859-1" > > Ed Brignole asked about setting a volume contouring level so that > the > resulting surface encloses a desired volume (cubic Angstroms). > Attached > is a Python script that does this in two ways. The first way is > fast > but approximate finding the level such that the number of enclosed > grid > points times the voxel volume equals the requested volume. > Estimating > the enclosed volume of the surface this way can give a much > different > result than actually computing analytically the volume within the > surface. The second approach computes surfaces at different > contour > levels and calculates analytically the volume. It uses bisection > to > find the correct contour level. This approach can be slow, about > 30 > seconds for a 256^3 volume (sfv.mrc virus map) with ~2.5e6 surface > triangles. Most of that time (80%) is actually computing the > enclosed > volume. That should be very fast but the code is checking for > holes in > the surface in a rather slow way. Both methods use the step size > (subsampling) set in the volume dialog. > > Tom >

2 1

Find contour level to enclose a specified volume
by Thomas Goddard 18 Jun '08

18 Jun '08

Ed Brignole asked about setting a volume contouring level so that the resulting surface encloses a desired volume (cubic Angstroms). Attached is a Python script that does this in two ways. The first way is fast but approximate finding the level such that the number of enclosed grid points times the voxel volume equals the requested volume. Estimating the enclosed volume of the surface this way can give a much different result than actually computing analytically the volume within the surface. The second approach computes surfaces at different contour levels and calculates analytically the volume. It uses bisection to find the correct contour level. This approach can be slow, about 30 seconds for a 256^3 volume (sfv.mrc virus map) with ~2.5e6 surface triangles. Most of that time (80%) is actually computing the enclosed volume. That should be very fast but the code is checking for holes in the surface in a rather slow way. Both methods use the step size (subsampling) set in the volume dialog. Tom # ----------------------------------------------------------------------------- # Find threshold level so that a contour surface encloses a specified volume. # # Method: Estimates how many grid points should lie within # surface then finds the level having that many grid points with higher values. # This method can produce a level where the enclosed volume is very # different from the requested volume. # def level_enclosing_grid_points(v, vol): (x0,y0,z0), (x1,y1,z1) = v.xyz_bounds() tvol = (x1-x0)*(y1-y0)*(z1-z0) r = 1.0 - vol/tvol s = v.matrix_value_statistics() level = s.rank_data_value(r) return level # ----------------------------------------------------------------------------- # Find threshold level so that a contour surface encloses a specified volume. # # Method: Computes surfaces for different levels and resulting volumes # using bisection to determine the correct level. # def level_by_bisection(v, vol, iterations = 10): m = v.matrix() l0, l1 = m.min(), m.max() ijk_step = v.region[2] from numpy import product voxel_volume = product(ijk_step) * product(v.data.step) mvol = vol / voxel_volume return find_level(m, l0, l1, mvol, iterations) # ----------------------------------------------------------------------------- # def find_level(m, l0, l1, mvol, iterations): l = 0.5*(l0+l1) if iterations <= 1: return l from _contour import surface varray, tarray = surface(m, l, cap_faces = True) from _surface import enclosed_volume vl, hole_count = enclosed_volume(varray, tarray) varray = tarray = None if vl == mvol: return l if vl > mvol: l0 = l else: l1 = l return find_level(m, l0, l1, mvol, iterations - 1) # ----------------------------------------------------------------------------- # vol = 3.227e6 # sfv.mrc, threshold level 1000. from chimera import openModels from VolumeViewer import volume_list for v in volume_list(): level1 = level_enclosing_grid_points(v, vol) level2 = level_by_bisection(v, vol) print v.name, level1, level2

2 1

3D surfaces from MacMolPlt or Gamess
by David Chenoweth 13 Jun '08

13 Jun '08

Dear Chimera team, I am having a bit of trouble reading a file into Chimera. Currently I am using Gamess to run DFT calculations and I can open these files in MacMolPlt which allows me to display orbitals or plot electrostatic potential surfaces. I can export the 3D surface data however I can't get Chimera to read the 3D grid files. Can you help me? I attached the file format info below and I can send you example files if you need. I am really trying to use Chimera for all my graphics needs wether it be viewing surfaces or orbitals from DFT calculations or looking at electron density maps from crystal structure data. If I could convert the 3D surface data files below into ccp4 or cns formated files that should also work for Chimera I just haven't found a way to do this yet. Thanks in advance, Dave ********************************************** David M. Chenoweth California Institute of Technology Division of Chemistry and Chemical Engineering Mail Code: 164-30 1200 California Boulevard, 91125 Pasadena California, USA Phone: 626-395-6074 Email: dchen(a)caltech.edu ********************************************** File format data info: http://www.scl.ameslab.gov/~brett/MacMolPlt 2D and 3D Surfaces (Import and Export) MacMolPlt implements a simple format for importing and exporting 2D and 3D surface data. To export a surface first create a surface and make sure that the surface is updated so that the grid data is available for export. Then choose the Export menu command and the Surface export format. MacMolPlt will then export the surface grid data to the file you specify (this is useful to check the input format as well). To import a 2D or 3D surface create a new General 2D or 3D surface. From the surface dialog you can then choose the file containing the grid data. The format for the 2D and 3D files are similar but there are a couple of differences. Floating point numbers may be given in a variety of formats and the grid data may contain any number of values per line. The units for the coordinate system should normally be angstroms. For 2D planes the format is: 1Line Label 30 # grid points -4.165882 -4.306789 0.839681 //Origin of the 2D grid 0.297019 0.000000 0.000000 //X inc vector 0.000000 0.297019 0.000000 //Y inc vector -0.000018 -0.000021... (# of grid points squared grid values) For 3D surfaces: 1 Line label 35 12 41 //nx ny nz -4.687880 -1.587537 -5.059790 //Origin of the 3D grid 0.280144 0.288642 0.275207 //x increment, y inc, z inc/ grid(x(y(z))) 0.000000 0.000000... (nx*ny*nz grid values)

2 2

3D surfaces from MacMolPlt or Gamess
by David Chenoweth 12 Jun '08

12 Jun '08

Dear Chimera team, I am having a bit of trouble reading a file into Chimera. Currently I am using Gamess to run DFT calculations and I can open these files in MacMolPlt which allows me to display orbitals or plot electrostatic potential surfaces. I can export the 3D surface data however I can't get Chimera to read the 3D grid files. Can you help me? I attached the file format info below with an exported density file and electrostatic potential map. I am really trying to use Chimera for all my graphics needs wether it be viewing surfaces or orbitals from DFT calculations or looking at electron density maps from crystal structure data. If I could convert the 3D surface data files below into ccp4 or cns formated files that should also work for Chimera I just haven't found a way to do this yet. Thanks in advance, Dave ********************************************** David M. Chenoweth California Institute of Technology Division of Chemistry and Chemical Engineering Mail Code: 164-30 1200 California Boulevard, 91125 Pasadena California, USA Phone: 626-395-6074 Email: dchen(a)caltech.edu ********************************************** File format data info: http://www.scl.ameslab.gov/~brett/MacMolPlt 2D and 3D Surfaces (Import and Export) MacMolPlt implements a simple format for importing and exporting 2D and 3D surface data. To export a surface first create a surface and make sure that the surface is updated so that the grid data is available for export. Then choose the Export menu command and the Surface export format. MacMolPlt will then export the surface grid data to the file you specify (this is useful to check the input format as well). To import a 2D or 3D surface create a new General 2D or 3D surface. From the surface dialog you can then choose the file containing the grid data. The format for the 2D and 3D files are similar but there are a couple of differences. Floating point numbers may be given in a variety of formats and the grid data may contain any number of values per line. The units for the coordinate system should normally be angstroms. For 2D planes the format is: 1Line Label 30 # grid points -4.165882 -4.306789 0.839681 //Origin of the 2D grid 0.297019 0.000000 0.000000 //X inc vector 0.000000 0.297019 0.000000 //Y inc vector -0.000018 -0.000021... (# of grid points squared grid values) For 3D surfaces: 1 Line label 35 12 41 //nx ny nz -4.687880 -1.587537 -5.059790 //Origin of the 3D grid 0.280144 0.288642 0.275207 //x increment, y inc, z inc/ grid(x(y(z))) 0.000000 0.000000... (nx*ny*nz grid values)

1 0

Swapaa from Java Runtime.exec
by ProtEUS 11 Jun '08

11 Jun '08

Hello all! I am sorry if this is a dump question, but is it possible for swapaa to run fromJava server with Runtime.exec()? Thank you very much!! Andreas, Greece.

2 3

translation + rotation of a whole scene
by Nils Becker 11 Jun '08

11 Jun '08

Hi list, I have a question about rotating and translating a bunch of models together while preserving their spatial relations. In particular, I would like to apply a matrix almost as in the matrixset command. THe difference is that I do not want to apply it individually to each model, but to one chosen models, while having the rest move together with the chosen model. Is that feasible in chimera? cheers, Nils

3 4

saving 2 maps for difference
by Gabriel Lander 09 Jun '08

09 Jun '08

I have 2 density maps, and use the "fit in map" function to best align the 2 of them. I'd like to use this fitting to create a difference map, but the original maps are not the same size. Does anyone know how to save both maps to the same grid in Chimera? For now I've been using the "subregion selection" tool in combination with "resampling voxel size", but if someone has a more direct method it'd be much appreciated. Thanks, -gabe

2 1

<< Plotting Gaussian MO >>
by Hatuey Hack 09 Jun '08

09 Jun '08

Hi all, I have a cube file with the Molecular Orbitals calculated with Gaussian. When I load the cub file into Chimera, I am only seen half of the orbital surface. I would like to know what should I do to get the whole molecular orbital surface colored (as usually) with only to colors (each on for one half). Regards, Hatuey

3 3

gaussian orbitals, Chimera-users Digest, Vol 62, Issue 11
by Jonathan Hilmer 09 Jun '08

09 Jun '08

Hatuey, If I remember correctly, you need to open the density map twice, then contour at the same level for both maps: one at the positive value and the other negative. Each map is colored independently, so they can be depicted however you want. Jonathan On Mon, Jun 9, 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. << Plotting Gaussian MO >> (Hatuey Hack) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Mon, 9 Jun 2008 10:24:29 -0300 > From: "Hatuey Hack" <hatuey.hack(a)gmail.com> > Subject: [Chimera-users] << Plotting Gaussian MO >> > To: chimera-users(a)cgl.ucsf.edu > Message-ID: > <612c95fa0806090624l5395ab11k6c9b235379e1a4fe(a)mail.gmail.com> > Content-Type: text/plain; charset=ISO-8859-1 > > Hi all, > > I have a cube file with the Molecular Orbitals calculated with Gaussian. > > When I load the cub file into Chimera, I am only seen half of the > orbital surface. > > I would like to know what should I do to get the whole molecular > orbital surface colored (as usually) with only to colors (each on for > one half). > > Regards, > > Hatuey > > > ------------------------------ > > _______________________________________________ > Chimera-users mailing list > Chimera-users(a)cgl.ucsf.edu > http://www.cgl.ucsf.edu/mailman/listinfo/chimera-users > > > End of Chimera-users Digest, Vol 62, Issue 11 > ********************************************* >

1 0

Molecular voodoo
by Tom Goddard 06 Jun '08

06 Jun '08

While debugging a bit of Chimera code I stumbled upon molecular voodoo. Stick pins in any evil molecule. Ask your doctor if this treatment may be right for you. A snake toxin (PDB 1a0m) has been fully cursed in the attached image. I attach the voodoo doll Chimera script for your amusement. You open voodoodoll.py in Chimera and then ctrl-right-mouse-button lets you place pins. Works in Chimera production release 1.2470 and newer. Tom def mouse_down_cb(viewer, event): from VolumeViewer import slice xyz_in, xyz_out = slice.clip_plane_points(event.x, event.y) from chimera import selection selection.clearCurrent() from VolumePath import place_marker place_marker(xyz_in) place_marker(xyz_out) from chimera import mousemodes as m m.addFunction('hide triangle', (mouse_down_cb, None, None)) m.setButtonFunction('3', ['Ctrl'], 'hide triangle')

1 0

smoothening the edge of molecular surface
by Maurice Ho 06 Jun '08

06 Jun '08

Is there any way to smoothen the edge of molecular surface when only part of it would be displayed? It's quite often that a couple of triangles sticking out from the edges which makes the image a bit absurd. Please kindly advice. rgds, Maurice Ho

3 2

Installing Chimera in Ubunty Hardy
by David 05 Jun '08

05 Jun '08

Hi, Im a long time user of Chimera, but now I changed XP for Ubuntu for 6 gb of ram that couldn't be used in XP 32 bits. I need help installing Linux version of Chimera, I have read and goggled web for instructions, and none of them seems to work, I run .exe in terminal (knowing linux doesn't use .exe i find this confusing) and nothing happens, I try that Chmod command and running again and nothing too... Im a really beginner here in linux, also browsed Ubuntu forums for help and at least till today have got no answer. All help would be greatly appreciated. sincerely, David Bulnes Abundis

2 1