- Chimera-users - RBVI Mailing Lists

Stereo angle value
by Prof. Fabrizio Briganti 02 May '06

02 May '06

Hi, I have not been able to find any reference to the value for the angle difference between Right and Left views in Stereo Images made with Chimera. Can you help me find it? Is there any way to change such angle value? Thanks Fabrizio

2 1

Paper virus models
by Thomas Goddard 02 May '06

02 May '06

I made a Chimera tool that flattens icosahedral virus models. Printing, cutting, folding, and taping produces a paper model. http://www.cgl.ucsf.edu/chimera/experimental/flatten_icosahedron/flaticos.h… The flatten icosahedron tool can be downloaded from the Chimera experimental features web page: http://www.cgl.ucsf.edu/chimera/experimental/experimental.html Tom

2 1

writing mrc files
by Eran Shacham 21 Apr '06

21 Apr '06

Hey there, I am trying to write an mrc file using the write_mrc2000_grid_data module. My function is below. My array of data is a <scitbx_array_family_flex_ext.double object at 0x58b9502c>. Which is a basically a matrix with dimensions: (243, 243, 243). When trying to write an mrc file using this array I get the following error: ************************* File "/home/eshacham/cctbx/mrc_functions.py", line 20, in write_mrc write_mrc2000_grid_data(g, path) File "CHIMERA/share/VolumeData/mrc/writemrc.py", line 44, in write_mrc2000_grid_data File "CHIMERA/share/VolumeData/arraygrid.py", line 48, in submatrix ValueError: need more than 1 value to unpack Code used *********************** def write_mrc(data): print data origin = (1.54, 2.45, 3.45) # XYZ coordinates of array index (0,0,0) step = (.25, .25, .25) # Spacing between grid points import Numeric a = Numeric.array(data, Numeric.Float32) arrays = [a] import VolumeData g = VolumeData.Array_Grid_Data(arrays, origin, step) path = '/home/eshacham/structures/result.mrc' from VolumeData.mrc import write_mrc2000_grid_data write_mrc2000_grid_data(g, path) I am assuming that it has to do with my array dimensions. Any help would be greatly appreciated. Cheers, Eran.

2 2

Re: [Chimera-users] Rotation about an axis
by Thomas Goddard 21 Apr '06

21 Apr '06

Hi Jordan, Aha! It turns out I wrote code that does exactly what you want last year. It is called the symmetric molecule copies command at the bottom of the Chimera experimental features web page: http://www.cgl.ucsf.edu/chimera/experimental/experimental.html It is not distributed with Chimera yet but you can download it on that page and install it with your existing Chimera. I hate to discourage you from Python programming in Chimera but I think this is just what you need. You open one copy of your molecule and then use this "sym" command to make the hexamer. Then when you move the original copy of the molecule relative to your map, the other 5 copies automatically reposition to maintain the 6-fold symmetry. To get a more detailed description look at the "details" link on the experimental features page. The hexamer symmetry comes from the BIOMT matrices in the PDB file. This sym command treats those matrices in the coordinate system of the density map (yes, a bit weird). The key here is to make sure the symmetry axis of the density map matches the symmetry axis of the BIOMT matrices. The easiest approach is probably to make the BIOMT matrices produce a hexamer about the z-axis through (0,0,0). When you open your density map in Chimera it may not have the hexamer density centered at (0,0,0). You will need to adjust its origin so the density hexamer is centered at (0,0,0). Use the volume dialog menu entry Features / Origin and Scale to adjust the position of the lower left front corner of the map. Assuming the hexamer is centered in a 128 x 128 x 100 grid with voxel size 0.8 angstroms you would set the origin to (-50.8, -50.8, -39.6) = (-63.5 * 0.8, -63.5 * 0.8, -49.5 * 0.8). The weird 63.5 = 128/2 - .5 instead of 64 = 128/2 is easiest to understand if you think of a map that is only 2 by 2 by 2 in size. Remember to press the Enter key in one of the origin fields in the volume dialog to have the new position take effect. In any case that may still be a half voxel off depending on your map, or it may be way off. Of course you have to know the position of the center of the hexamer in your density to make this work. A way to verify that you have set the correct origin is to open 2 copies of the map -- set the map origins, then use constrained move to rotate just one copy keeping the other one fixed. It is easy then to visually see if the two copies are concentric. For convenience you can save your map as an MRC file with the new origin with the volume dialog File / Save As... menu entry. Saving a Chimera session also saves the origin in the Chimera session file (the map is not modified). I can see various troubles with using the Multiscale tool. First Multiscale low resolution surfaces are not reliably synchronized to the individual PDB models. If you change a PDB model position Multiscale won't realize it and the low resolution surface will not move to the new location. If you are using the latest Chimera snapshot (1.2224, April 10, 2006) then the problems will be even more severe because I changed the way Multiscale positions PDB models. Instead of setting the pdb_model.openState.xform it changes the atom coordinates of the model (setting atom.coord = x,y,z). All the xforms of the 6 models in a hexamer would be the same in Chimera 1.2224. This change was made to allow correctly moving individual chains of multiscale models, possible with the "movement mouse mode" feature on the experimental features web page. Lastly some clarification about Chimera's coordinate systems. First there is a global Chimera coordinates system with the x axis always horizontal relative to your graphics window, the y axis always vertical, and the z axis always pointing out of the screen at you. When you rotate models with the mouse, this coordinate system does not change -- the individual matrices positioning the models are changed. The origin of this global coordinate system is not necessarily in the center of the screen. It can be anywhere including off-screen. The global coordinates for the center of the screen are set to match the center of the bounding box of the first model you open in Chimera. Besides the global coordinate system each model has a local coordinate system. With a PDB file, the atom coordinate positions are in the local coordinate system. If you open a PDB model in Chimera and rotate it, the atom.coord x,y,z values do not change -- they are expressed in the local coordinate system which is fixed to the atoms and rotates with the model. (You can however change the atom.coord positions using Python code.) The pdb_model.openState.xform 3 row by 4 column matrix gives the rotation and translation that maps the local coordinates (e.g. a PDB atom position) into Chimera global coordinates. Other types of Chimera models like density maps are just like PDB models having a local coordinate system -- the density map grid points have unchanging coordinates in the local coordinate system. More than you wanted to know. Tom > Date: Thu, 20 Apr 2006 19:51:40 -0700 > From: "Jordan Benjamin" <jordan.benjamin(a)gmail.com> > To: "Thomas Goddard" <goddard(a)cgl.ucsf.edu> > Subject: Re: [Chimera-users] Rotation about an axis > > Ah, it looks like it still works for normal models. The weird behavior > was with a model that had been generated with the Multiscale Models > tool from crystal contacts. I still wish I understood the > global-vs-local issue better though. > > > On 4/20/06, Jordan Benjamin <jordan.benjamin(a)gmail.com> wrote: > > Thank you for your help on this. I have a further question. Basically, > I'm trying to dock into electron density the crystal structure of a > protein that forms a hexamer. I would like to be able to have six > copies of the monomer structure loaded at once, and when I move or > rotate one copy, I'd like to be able to transform the other five to > match. So, I thought the ConstrainedMove code would be a good starting > point, since what I'm trying to do is the same thing it does, just > with the rotation specified by degrees instead of by dragging the > mouse. > > After you sent the email below, I modifed the function I was writing > accordingly, and it worked. However, I'm trying to use it now, and > it's not working. I'm a little lost on debugging it, because I don't > even know what Chimera's concept of a global coordinate system is, nor > do I have any way to visually display the local origin and axis for > each object. Is there a way to do these things that you know of? Do > you know of a simple way to rotate an object X degrees around a > specified axis (relative to another object if needed)? > > If these are really hard questions and the answer is just that I > should read the source code until I figure it out, that's totally > fine. I just thought I'd ask. > > At least with Chimera the code IS available, unlike the situation with > Amira where things would be undocumented and broken, and there'd be no > way to fix them... > > Thanks so much, > Jordan > > The function I was using is as follows (yes, it's brittle and ugly, > but it's all I had time for): > def update_hexamer(ref_model,model,*rest): > origin = (0,0,0) > axis = (0,0,1) > angle_deg = 60 > > mvo = chimera.Point(*origin) > vo = ref_model.openState.xform.apply(mvo) > mva = apply(chimera.Vector, axis) > va = ref_model.openState.xform.apply(mva) > > xf_original = model.openState.xform > > for i in range(5): > xf = chimera.Xform() > xf.translate(*vo) > xf.rotate(va, angle_deg+60*i) > xf.translate(-vo.x, -vo.y, -vo.z) > model.openState.globalXform(xf) > rest[i].openState.xform = model.openState.xform > model.openState.xform = xf_original >

1 0

Re: [Chimera-users] Sum of volume data values
by Thomas Goddard 13 Apr '06

13 Apr '06

Hi Eran, Ok, here is code to read an MRC file and print the sum of map values without using Chimera. You need Python, the Numeric Python module which handles multidimensional arrays and the VolumeData module from Chimera. Numeric is obtained from sourceforge: http://sourceforge.net/project/showfiles.php?group_id=1369&package_id=1351 The VolumeData module is in your Chimera distribution. chimera/share/VolumeData Here's an example running the script where I've set PYTHONPATH to have access to the VolumeData module: % env PYTHONPATH=/usr/local/chimera/share python2.4 volumesum2.py Sum of values of groel.mrc = 21720.8984375 Send future questions to chimera-users(a)cgl.ucsf.edu so others who have similar problems can benefit. I get that email just as fast. Tom ------- file volumesum2.py follows: # Script to print the sum of volume data values. # Does not need Chimera, just the VolumeData module and Numeric and Python. path = 'groel.mrc' # Read the MRC volume data and get the 3-D Numeric array. from VolumeData.mrc.mrc_format import MRC_Data data = MRC_Data(path) m = data.matrix() # Compute the sum of the elements of the volume matrix. import Numeric sum = Numeric.sum(m.flat) # Print the sum print 'Sum of values of', path, ' = ', sum ------- > Date: Wed, 12 Apr 2006 14:46:14 -0700 > From: Eran Shacham <eshacham(a)burnham.org> > To: Thomas Goddard <goddard(a)cgl.ucsf.edu> > Subject: Re: Sum of volume data values > > Hey Tom, > > Thanks for your reply. However, I would like to do that from within a > script I already wrote. i.e I would like to integrate a read_in method > for an MRC file, access it values and then calculate the sum. I would > like to avoid doing that by opening my file in Chimera at this stage of > the program. Can it be done by importing some of your code? > > Thanks, > Eran.

2 2

Sum of volume data values
by Thomas Goddard 12 Apr '06

12 Apr '06

Hi Eran, You can open your MRC file in Chimera, then open the following script "volumesum.py" which contains Python code to have the sum of all the volume values printed. Use menu entry File / Open to open both files. If you want to do this many times I can give you a modified version that creates a Chimera menu entry (under the Tools menu) that runs the script. Tom ------------ File volumesum.py follows: # Script to print the sum of volume data values. # Get the volume dialog object import VolumeViewer d = VolumeViewer.volume_dialog() # Get the volume data currently shown by the volume dialog dr = d.focus_region m = dr.full_matrix() # Compute the sum of the elements of the volume matrix. import Numeric sum = Numeric.sum(m.flat) # Print the sum to the reply log dialog. message = 'Sum of values of %s = %.6g' % (dr.full_name(), sum) print message # Display the reply log dialog. from chimera import dialogs dialogs.display('reply') ------------- > Date: Wed, 12 Apr 2006 13:52:26 -0700 > From: Eran Shacham <eshacham(a)burnham.org> > To: Thomas Goddard <goddard(a)cgl.ucsf.edu> > Subject: Reading mrc files > > Hi Tom, > > Can I use Chimera's code for reading in mrc files? Essentially, I would > like to read in an mrc file and access its values. It is a mask file, > and I simply need to add up it values. Any suggestions on how to do that > easily with a Chimera module? > > Cheers, > Eran

1 0

Volume series tool
by Thomas Goddard 12 Apr '06

12 Apr '06

I wrote a Chimera tool to display sequences of volume data sets. They are displayed one at a time like a movie and a slider lets you choose any particular volume. It is intended for looking at 3-D optical microscopy time series, or a sequence of single particle density map reconstructions. Available on the Chimera experimental features web page: http://www.cgl.ucsf.edu/chimera/experimental/experimental.html Tom

1 0

Rotation about an axis
by Jordan Benjamin 11 Apr '06

11 Apr '06

I'm trying to rotate a model about a specified axis and origin (preferably in global coordinates). I would like the rotation to be applied on top of whatever transformations have been applied already. I've looked at the documentation, including the programmer's FAQ and the source code for ConstrainedMove. I couldn't figure out how to do this, so I wrote the following function based on the ConstrainedMove code to do the rotation relative to another model: def rotate_around_model_axis(ref_model, model, origin, axis, angle_deg): mvo = apply(chimera.Vector, origin) vo = ref_model.openState.xform.apply(mvo) mva = apply(chimera.Vector, axis) va = ref_model.openState.xform.apply(mva) xf = chimera.Xform() xf.translate(vo) xf.rotate(va, angle_deg) xf.translate(-vo) model.openState.globalXform(xf) But, this doesn't seem to work as I'd intended. Is there a way to do this? Also, I'm a little confused about the coordinate system. It seems like the concept of a global coordinate system may not apply in Chimera the way I'm used to (like in Amira for example) since the camera is fixed. Is this true, or am I totally lost? Thanks in advance for any help you can provide. Jordan

3 2

Newbie Questions
by Shafrir, Yinon (NIH/NCI) [F] 11 Apr '06

11 Apr '06

Hi there! I have two questions regarding using Chimera. How do you break a molecule to different models using the command line? I have a protein composed of 4 chains. It loads as one piece and designated model 0 I then select a single chain and I want to designate it as a new model. How do I do that using the command line (so I can write a script)? Also, when I select the chain and use the write command in this manner "write selected su1.pdb" the save dialog box pops up (even though I designated a filename) When I try "write selected 0 su1.pdb" nothing happens. My other question is: How can you force chimera to assign secondary type display to residues without editing the PDB file first? Say I want to show 30 residues that are not detected by ksdssp as beta sheets. Still, I want to represent them as beta sheets. Do I need to edit the pdb file? Best regards Yinon Shafrir

4 3

Axis display
by Jordan Benjamin 10 Apr '06

10 Apr '06

Is there a convenient way to display the axes of the coordinate system (and preferably the origin)? Forgive me if this is addressed in the documentation, but I wasn't able to find it. Thanks, Jordan

2 1

Collaboratory session
by Eduardo Sanz Garcia 10 Apr '06

10 Apr '06

I have been trying to set up a collaboratory session, but I could not see the model. I was able to chat and to see the pointer of the other person, anything else. Any idea? How can I know if it is a problem caused by the firewall?

2 1

Move atoms with mouse
by Thomas Goddard 07 Apr '06

07 Apr '06

I put a tool on the Chimera experimental features web page that enhances the standard rotation and translation mouse modes (left and middle mouse buttons) to enable moving helices, strands, turns, PDB chains, any selected set of atoms, or multiscale model chains. http://www.cgl.ucsf.edu/chimera/experimental/experimental.html The molecular pieces are translated or rotated by dragging them with the mouse. It requires Chimera version 1.2206 or newer, and that is only available for Macintosh computers right now. Tom

1 0

Newbie Questions
by Shafrir, Yinon (NIH/NCI) [F] 07 Apr '06

07 Apr '06

Hi there! I have two questions regarding using Chimera. How do you break a molecule to different models using the command line? I have a protein composed of 4 chains. It loads as one piece and designated model 0 I then select a single chain and I want to designate it as a new model. How do I do that using the command line (so I can write a script)? Also, when I select the chain and use the write command in this manner "write selected su1.pdb" the save dialog box pops up (even though I designated a filename) When I try "write selected 0 su1.pdb" nothing happens. My other question is: How can you force chimera to assign secondary type display to residues without editing the PDB file first? Say I want to show 30 residues that are not detected by ksdssp as beta sheets. Still, I want to represent them as beta sheets. Do I need to edit the pdb file? Best regards Yinon Shafrir

1 0

Re: [Chimera-users] Read PDB file
by Elaine Meng 07 Apr '06

07 Apr '06

Hi Eduardo, Aha! You can see all the header information from the PDB file (the lines that come before the coordinates). How to do that: open the Model Panel (Favorites... Model Panel) and for the model you are interested in, click "attributes..." Click the "PDB Headers..." button at the bottom of the resulting attributes panel. Currently the header lines are shown in a different order than in the original PDB file (they are sorted alphabetically by type, but wthin each type such as REMARK the order is not changed), so you will have to scroll down to the SOURCE record to see the organism. Best, Elaine On Apr 7, 2006, at 11:29 AM, Eduardo Sanz Garcia wrote: > My real situation: > 1. I have File... Fletch by ID using 2ypi > 2. I didn't know what source organism they are using to get the > structure. > > Option a: go to www.pdb.org and obtain the information. > Option from chimera?: display the text file 2ypi.pdb to get the > necessary information? > > I was wondering if there is a command or function to display (the > text) of a PDB file. > ----- 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

Read PDB file
by Eduardo Sanz Garcia 07 Apr '06

07 Apr '06

Hi, Is it possible to read the PDB file of the molecule directly from chimera? Thank you Eduardo Sanz

2 1

Name selection
by Eduardo Sanz Garcia 06 Apr '06

06 Apr '06

Hi everybody, Is there any way to give a name to a selection from the command line? Same as: *Select... Name Selection...* <cid:part1.01070107.06060601@terra.es> but from the command line. Thank you.

2 1

further levels of detail in low-res surface of a single chain Multiscale model
by Klinkhardt＠em.uni-frankfurt.de 05 Apr '06

05 Apr '06

Hello, I would like to divide/define the low-res structure of a multiscale model which contains only one chain. I need further levels of detail in this single chain (preferentially all single nucleotide bases as small podgy shapes and the backbone of phosphate-sugars as an irregular rope) so I can can animate the resulating smaller units later. In your online hept (documentation) there is written that a python script like #bluetongue.py# can create a multiscale model with further levels of (individual?) structure. I have opend the python script but unfortunately I am no programmer so I am not able to modify this script to my requirements. Is there a #more general# script which I can use for my chain or any other workaround?? Thanks a lot in advance, best regards Ute Klinkhardt

2 1

question about CHIMERA
by Guy Nimrod 04 Apr '06

04 Apr '06

Hello, I am a PhD student in the department of Biochemistry of the Tel-Aviv University and often use your great application, CHIMERA. Is it possible to view on chimera the selected part only? Regarding this, is it possible to split a model to its units (so each chain/ligand) will be considered as a model? Thanks Guy Nimrod http://ashtoret.tau.ac.il/~nimrodg

4 5

Please consider adding PyTables to Chimera
by David E. Konerding 04 Apr '06

04 Apr '06

Hello, Recently I have been converting most of my scientific data to the HDF5 format. HDF5 is a powerful, high performance data format library. I found an absolutely wonderful open source code, PyTables, which provides an HDF5 interface. It adds enormous value to HDF5 by making it much easier to use, and exposing a lot of the complex functionality that HDF5 supports (such as variable-length records). http://www.pytables.org/moin I would really appreciate it if you guys would consider adding PyTables. I would definitely convert a lot of my existing extensions and processing tools to use PyTables- for example, when I analyze trajectories to find hydrogen bond data and use that to render long movies, PyTables makes it a lot easier (I already used the Chimera headers to build my own copy of PyTables on Linux, but it would be a lot more painful to do on Windows). PyTables depends on numarray (although it also provides Numeric Python support). I'd like to see numarray in Chimera as well (although I could live with just Numeric support) Dave

3 5

trigonometry
by jean-francois menetret 03 Apr '06

03 Apr '06

What is Chimera's trigonometry convention ? How can get angles out of the rotation matrix (saved with "matrixget") ? Jean-Fran�ois M�n�tret, PhD Boston University School of Medicine Physiology and Biophysics Department 700 Albany Street W315 Boston, MA 02118 Email: menetret(a)bu.edu Mailing address: 715 Albany Street

3 4

Different colors for individual chains within a model
by Karpen, Saul J. 31 Mar '06

31 Mar '06

Dear Chimera group, Simple question from a new user--How do I distinctly color ribbon representations of different chains within the same model? Thanks, Saul J. Karpen, M.D., Ph.D. Baylor College of Medicine

3 2

Re: [Chimera-users] Question about background transparency
by Weldon, Scott 28 Mar '06

28 Mar '06

Thanks to all for the help. I have transparent images now. Thanks, Scott Weldon, M.A. Supervisor, Medical Illustrator Division of Cardiothoracic Surgery Michael E. DeBakey Department of Surgery Baylor College of Medicine One Baylor Plaza, BCM 390 Houston, TX 77030 TEL: (832) 355-9941 FAX: (832) 355-9948 E-mail: sweldon(a)bcm.edu -----Original Message----- From: Thomas Goddard [mailto:goddard@cgl.ucsf.edu] Sent: Tuesday, March 28, 2006 11:43 AM To: Weldon, Scott Cc: chimera-users(a)cgl.ucsf.edu Subject: Re: [Chimera-users] Question about background transparency Hi Scott, You need 2 dashes in the background opacity option --bgopacity According to the Chimera manual that option should occur before input PDB files in the Chimera startup command. Here are the relevant pages of the Chimera User's Guide: http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/startup.html http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/options.html Tom

1 0

Question about background transparency
by Weldon, Scott 28 Mar '06

28 Mar '06

Hello, I'm not having any luck using command line to change the background to transparent. I've tried "-bgopacity" and "bgopacity" and both say it is not a valid command. Would you please provide the exact syntax that should be entered to achieve this? Does this need to be entered before the PDB file is opened? Thank you for your help, Scott Weldon, M.A. Supervisor, Medical Illustrator Division of Cardiothoracic Surgery Michael E. DeBakey Department of Surgery Baylor College of Medicine One Baylor Plaza, BCM 390 Houston, TX 77030 TEL: (832) 355-9941 FAX: (832) 355-9948 E-mail: sweldon(a)bcm.edu

3 2

ses volume
by Charles Spink 21 Mar '06

21 Mar '06

Is it possible to get a calculated volume for the solvent excluded area for a protein or dna? I know the volume for a solvent accessible surface is possible through "measure volume" command. What about SES volumes? Thanks, Charlie Spink

3 2

superimposing proteins on co-factor
by Boaz Shaanan 08 Mar '06

08 Mar '06

Hi Chimera prople, If I want to superimpose proteins on their co-factors (or part of them), is ensembleMatch the way to go or are there other ways to do this ? I'd appreciate your advice. Regards, Boaz -- ************************************************************* Boaz Shaanan, Ph.D. Associate Prof. Dept. of Life Sciences Ben-Gurion University of the Negev Beer-Sheva 84105 Israel Phone: 972-8-647-2220 e-mail:bshaanan@bgu.ac.il Fax: 972-8-647-2992 or 972-8-646-1710 *************************************************************

3 2