February 2008 - Chimera-users - RBVI Mailing Lists

Making an omnimax POVray movie
by Kristina Furse 26 Feb '08

26 Feb '08

Hi- We're interested in using Chimera to prepare POVray files for a movie to be shown on a hemispherical dome. The first part of the movie will be a 360 degree rotation of a DNA/drug complex showing water occupancy around the drug, and the second half will come from an AMBER MD trajectory. We've figured out the content and rendering. The two remaining issues are that the camera must be changed from "perspective" to "omnimax" (fisheye), and the extensive rendering time that will be required for the large system at very high resolution (4096x4096). Here is our plan: 1. Use a command-line script to rotate the system and export the scene to x3d files. 2. Use per-frame scripting in MD movie to export each trajectory frame to x3d files. 3. Use x3d2pov to translate into *.pov files. 4. Change the camera from perspective to omnimax in each .pov file. 5. Render the individual images in parallel on our cluster using povray (different version than the one that ships with Chimera). 6. Compress to jpeg. 7. Assemble movie as .mov (format compatible with projection system in the dome). Does anyone have a better idea? I'm assuming that I can't change the camera to omnimax from within Chimera? Thanks for any suggestions! Kristina -- Kristina Furse Postdoctoral Research Associate 262 Stepan Chemistry Hall Notre Dame, IN 46556 (574)631-3904

2 2

[chimera-users] --nogui
by 吴超 25 Feb '08

25 Feb '08

Hi, folks! I came across with a problem when using chimera --nogui file.cmd > data.log the file.cmd is as follows: open /home/lynn/project/hivrt_pdb/holo_without_mutation/1C1C.pdb findclash :.het@C1 overlap -1 hbond 0 ignoreIntraRes true log true findclash :.het@C2 overlap -1 hbond 0 ignoreIntraRes true log true findclash :.het@C3 overlap -1 hbond 0 ignoreIntraRes true log true findclash :.het@C4 overlap -1 hbond 0 ignoreIntraRes true log true findclash :.het@C5 overlap -1 hbond 0 ignoreIntraRes true log true findclash :.het@C6 overlap -1 hbond 0 ignoreIntraRes true log true findclash :.het@S overlap -1 hbond 0 ignoreIntraRes true log true findclash :.het@C7 overlap -1 hbond 0 ignoreIntraRes true log true findclash :.het@C12 overlap -1 hbond 0 ignoreIntraRes true log true findclash :.het@CF overlap -1 hbond 0 ignoreIntraRes true log true findclash :.het@CH overlap -1 hbond 0 ignoreIntraRes true log true findclash :.het@CI overlap -1 hbond 0 ignoreIntraRes true log true findclash :.het@C11 overlap -1 hbond 0 ignoreIntraRes true log true findclash :.het@O11 overlap -1 hbond 0 ignoreIntraRes true log true findclash :.het@N8 overlap -1 hbond 0 ignoreIntraRes true log true findclash :.het@C9 overlap -1 hbond 0 ignoreIntraRes true log true findclash :.het@O9 overlap -1 hbond 0 ignoreIntraRes true log true findclash :.het@N10 overlap -1 hbond 0 ignoreIntraRes true log true findclash :.het@CA overlap -1 hbond 0 ignoreIntraRes true log true findclash :.het@OB overlap -1 hbond 0 ignoreIntraRes true log true findclash :.het@CC overlap -1 hbond 0 ignoreIntraRes true log true findclash :.het@CD overlap -1 hbond 0 ignoreIntraRes true log true If i use the graphic version and read the cmd file, it worked quite well and I can find all I need in the reply log. But when it came to the nogui version, the data file only contains the result of the first two lines of the cmd file. That is open /home/lynn/project/hivrt_pdb/holo_without_mutation/1C1C.pdb findclash :.het@C1 overlap -1 hbond 0 ignoreIntraRes true log true Results are: Opening 1C1C.pdb... 1 model opened Opened 1C1C.pdb containing 1 model, 7997 atoms, and 1096 residues Allowed overlap: -1 H-bond overlap reduction: 0 Ignore overlap between atoms separated by 4 bonds or less Ignore intra-residue overlaps 5 overlaps atom1 atom2 overlap distance #0:999.het@C1 #0:100.A@CD1 -0.352 4.112 #0:999.het@C1 #0:100.A@CD2 -0.643 4.403 #0:999.het@C1 #0:188.A@CD2 -0.723 4.363 #0:999.het@C1 #0:188.A@CB -0.883 4.643 #0:999.het@C1 #0:188.A@CG -0.989 4.479 Error in the shell is as follows: $ /usr/local/chimera/bin/chimera --nogui /home/lynn/cmd/1C1C.cmd > /home/lynn/project/data/distance/1C1C.dist Traceback (most recent call last): File "/usr/local/chimera/share/__main__.py", line 65, in ? value = chimeraInit.init(sys.argv) File "/usr/local/chimera/share/chimeraInit.py", line 358, in init replyobj.reportException( File "CHIMERA/share/chimera/replyobj.py", line 610, in reportException File "CHIMERA/share/chimera/replyobj.py", line 438, in nogui_message ValueError: I/O operation on closed file $ Could anyone give me hand? I am really frustrated! Thanks in advance... -- Chao Wu Currently: National Institute of Biological Sciences, Beijing Huang Niu's Lab http://www.nibs.ac.cn/ Tel: 80726688-8572 Address: Beijing ZhongGuanCun Life Science Park, NO.7 Science Park East Road, NIBS B427 Postcode: 102206 Life Science Honors Program, China Agricultural University PO.1055, NO.2 YuanMingYuan west road, Beingjing, 100094 P.R. China Tel: 86-10-62812484 E-mail: lynnwuchao(a)gmail.com

2 1

NOT ABLE TO VIEW SURFACE USING PHYRE MODEL DATA
by nyyount＠cox.net 25 Feb '08

25 Feb '08

Hi, I have not been able to view my Phyre pdb files using Chimera�s surface view. I get the error code -1073741785 and am immediately shut out of the program. I am not a computer person � yet. Are there any quick fixes that I can add to my Phyre generated PDB files to be able to view them in �surface� view? I�m very pleased with Chimera and have used it for all of the molecular graphics for our next publication. I would like to be able to view my Phyre models using the same software. Thank you! Nannette Y. Yount, Ph.D. Harbor-UCLA

3 2

save manipulated electron density maps in spider
by Jan-Christian Schuette 25 Feb '08

25 Feb '08

Dear Chimera Developers, I quite often have to manipulate electron density maps in the spider format. When using chimeras volume eraser function, it doesn't allow me to save the new map back to spider or another easy-to-convert file format (such as brix or ccp4). I would be great if you could add such a feature in the next version of the software. Best regards, Jan-Christian Schuette

3 2

Re: [Chimera-users] distance
by Elaine Meng 22 Feb '08

22 Feb '08

Hi! Eric has kindly provided a python script (attached at the bottom) that will write all the current distance measurements to the Reply Log. Just open the python file with Chimera and it will happen. Alternatively, depending on your situation, you could use "findclash," as discussed in these chimera-users messages a few months ago: http://www.cgl.ucsf.edu/pipermail/chimera-users/2007-November/ 002061.html http://www.cgl.ucsf.edu/pipermail/chimera-users/2007-November/ 002069.html http://www.cgl.ucsf.edu/pipermail/chimera-users/2007-November/ 002070.html I hope this helps, Elaine ----- Elaine C. Meng, Ph.D. meng(a)cgl.ucsf.edu UCSF Computer Graphics Lab and Babbitt Lab Department of Pharmaceutical Chemistry University of California, San Francisco http://www.cgl.ucsf.edu/home/meng/index.html On Feb 21, 2008, at 5:24 PM, 吴超 wrote: > Thanks, Elaine! > The problem is that I want to automatically calculate a number of > distances in non-gui interface, and sent the date to a file, just > like the logfile. > So is there a chance that I can achieve it using chimera? > > On Fri, Feb 22, 2008 at 1:00 AM, Elaine Meng <meng(a)cgl.ucsf.edu> > wrote: > > The distances are listed in a separate dialog - choose "Tools... > Structure Analysis... Distances" from the menu to open that dialog. > You can adjust the number of decimal places in the measurement, and > there is a "Save" button if you want to save distances or other > measurements to a text file.

3 2

Driving matchAtoms from the IDLE
by Jean-Didier Maréchal 22 Feb '08

22 Feb '08

Hi Guys, I'd like to superpose two inorganic complexes using two differently "ordered" atom sets. My compound has quite similar substructures but I am sure at all what is the best "ordered" selection of atom pairs to align. Therefore, I'd like to make a "extensive" search (with an heuristic approach to define...) passing by the idle. I am unfortunately missing something before getting to this point! Lets say that I just want to align the nitrogen atoms of the my two structures. I have: import chimera import numpy from chimera import match from chimera import runCommand runCommand("select #1@N=") an1=chimera.selection.currentAtoms() runCommand("select #0@N=) an2=chimera.selection.currentAtoms() chimera.match.matchAtoms(an1,an2) And then it's when I am already stuck! the resulting prompt does not provide mistake an gives: chimera.match.matchAtoms(an1,an2) (<_chimera.Xform object at 0x9463ef8>, 0.034256915079603358) But how this results can be send to the graphical window? I guess I have to use the Xform object, but I can not clearly understand how to do this. Thanks for any help! JD On Thu, 2008-02-21 at 12:00 -0800, chimera-users-request(a)cgl.ucsf.edu wrote: > Send Chimera-users mailing list submissions to > chimera-users(a)cgl.ucsf.edu > > To subscribe or unsubscribe via the World Wide Web, visit > http://www.cgl.ucsf.edu/mailman/listinfo/chimera-users > or, via email, send a message with subject or body 'help' to > chimera-users-request(a)cgl.ucsf.edu > > You can reach the person managing the list at > chimera-users-owner(a)cgl.ucsf.edu > > When replying, please edit your Subject line so it is more specific > than "Re: Contents of Chimera-users digest..." > > > Today's Topics: > > 1. Re: distance (Elaine Meng) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Thu, 21 Feb 2008 09:00:37 -0800 > From: Elaine Meng <meng(a)cgl.ucsf.edu> > Subject: Re: [Chimera-users] distance > To: "??" <lynnwuchao(a)gmail.com> > Cc: chimera-users(a)cgl.ucsf.edu > Message-ID: <0E85C396-8C7D-4182-A029-5C0155245C3A(a)cgl.ucsf.edu> > Content-Type: text/plain; charset=UTF-8; delsp=yes; format=flowed > > Hi! > > The distances are listed in a separate dialog - choose "Tools... > Structure Analysis... Distances" from the menu to open that dialog. > You can adjust the number of decimal places in the measurement, and > there is a "Save" button if you want to save distances or other > measurements to a text file. > > I hope this helps, > Elaine > ----- > Elaine C. Meng, Ph.D. meng(a)cgl.ucsf.edu > UCSF Computer Graphics Lab and Babbitt Lab > Department of Pharmaceutical Chemistry > University of California, San Francisco > http://www.cgl.ucsf.edu/home/meng/index.html > > > > > On Feb 20, 2008, at 11:38 PM, ?? wrote: > > > Hi, folks, > > I select two atoms and measured their distance using the command > > distance and can see the exact distance value in the graphic > > interface, but when I go the log file i can not find it. > > Can anybody give me a hand and tell me why can I find it? Thanks.. > > _______________________________________________ > > Chimera-users mailing list > > Chimera-users(a)cgl.ucsf.edu > > http://www.cgl.ucsf.edu/mailman/listinfo/chimera-users > > > > ------------------------------ > > _______________________________________________ > 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 58, Issue 24 > *********************************************

2 1

placing a helium atom at a given position.
by Jean-Didier Maréchal 22 Feb '08

22 Feb '08

Hi again, that follows an e-mail exchange with Eric a while ago. I am trying to place a given atom (for the moment, helium is fine) to given x,y,z coordinates. You told me that the new BuildStructure module could do that, unfortunately after having imported the module import BuildStructure when I do BuildStructure.placeHelium(res,model=0,position=13.2909,75.4474,25.4248) I have SyntaxError: non-keyword arg after keyword arg I tried different argumentations of this but I can't have it working, if you can tell me how to do this, that would be perfect. All the best, JD some time ago Eric On Thu, 2008-02-21 at 12:00 -0800, chimera-users-request(a)cgl.ucsf.edu wrote: > Send Chimera-users mailing list submissions to > chimera-users(a)cgl.ucsf.edu > > To subscribe or unsubscribe via the World Wide Web, visit > http://www.cgl.ucsf.edu/mailman/listinfo/chimera-users > or, via email, send a message with subject or body 'help' to > chimera-users-request(a)cgl.ucsf.edu > > You can reach the person managing the list at > chimera-users-owner(a)cgl.ucsf.edu > > When replying, please edit your Subject line so it is more specific > than "Re: Contents of Chimera-users digest..." > > > Today's Topics: > > 1. Re: distance (Elaine Meng) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Thu, 21 Feb 2008 09:00:37 -0800 > From: Elaine Meng <meng(a)cgl.ucsf.edu> > Subject: Re: [Chimera-users] distance > To: "??" <lynnwuchao(a)gmail.com> > Cc: chimera-users(a)cgl.ucsf.edu > Message-ID: <0E85C396-8C7D-4182-A029-5C0155245C3A(a)cgl.ucsf.edu> > Content-Type: text/plain; charset=UTF-8; delsp=yes; format=flowed > > Hi! > > The distances are listed in a separate dialog - choose "Tools... > Structure Analysis... Distances" from the menu to open that dialog. > You can adjust the number of decimal places in the measurement, and > there is a "Save" button if you want to save distances or other > measurements to a text file. > > I hope this helps, > Elaine > ----- > Elaine C. Meng, Ph.D. meng(a)cgl.ucsf.edu > UCSF Computer Graphics Lab and Babbitt Lab > Department of Pharmaceutical Chemistry > University of California, San Francisco > http://www.cgl.ucsf.edu/home/meng/index.html > > > > > On Feb 20, 2008, at 11:38 PM, ?? wrote: > > > Hi, folks, > > I select two atoms and measured their distance using the command > > distance and can see the exact distance value in the graphic > > interface, but when I go the log file i can not find it. > > Can anybody give me a hand and tell me why can I find it? Thanks.. > > _______________________________________________ > > Chimera-users mailing list > > Chimera-users(a)cgl.ucsf.edu > > http://www.cgl.ucsf.edu/mailman/listinfo/chimera-users > > > > ------------------------------ > > _______________________________________________ > 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 58, Issue 24 > *********************************************

2 1

managing bond order
by Victor Kenyon 21 Feb '08

21 Feb '08

hi all, typically when i import a paired back active site of a protein, chimera wants to connect atoms that should not be connected. is there a way to select bonds that violate a certain bond length so i can remove them? or a possible work around? thanks for your help. victor Victor A. Kenyon kenyon(a)chemistry.ucsc.edu 831.252.7667.

2 1

distance
by 吴超 21 Feb '08

21 Feb '08

Hi, folks, I select two atoms and measured their distance using the command distance and can see the exact distance value in the graphic interface, but when I go the log file i can not find it. Can anybody give me a hand and tell me why can I find it? Thanks..

2 1

modifying presets
by glenn millhauser 20 Feb '08

20 Feb '08

Hi Folks, Is there a way to edit the Presets to include additional protein features? Specifically, I work with quite a few disulfide rich proteins and it would be great to select the Preset Interactive 1 and have it automatically draw Cys side chains with disulfide links. Many thanks, glenn Glenn L. Millhauser Department of Chemistry & Biochemistry UC Santa Cruz Santa Cruz, CA 95064 831 459 2176 voice 831 566 3337 cell 831 459 2935 fax http://chemistry.ucsc.edu/~glennm http://www.chemistry.ucsc.edu/faculty/millhauser.html

2 4

Chimera-users February 2008