- Chimera-users - RBVI Mailing Lists

new Chimera production release
by Eric Pettersen 14 May '04

14 May '04

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

1 0

by riccardodesantis＠fastwebnet.it 05 May '04

05 May '04

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

2 1

by riccardodesantis＠fastwebnet.it 05 May '04

05 May '04

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

1 0

Coloring by bfactor
by Thomas Goddard 23 Apr '04

23 Apr '04

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

2 1

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

23 Apr '04

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

1 0

install question
by Randy Heiland 23 Apr '04

23 Apr '04

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

2 1

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

23 Apr '04

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

2 1

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

22 Apr '04

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

2 1

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

13 Apr '04

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

1 0

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

09 Apr '04

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

1 0

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

08 Apr '04

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

1 0

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

08 Apr '04

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

2 1

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

02 Apr '04

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

2 1

correction re sessions
by Elaine Meng 01 Apr '04

01 Apr '04

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

1 0

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

01 Apr '04

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

1 0

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

01 Apr '04

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

1 0

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

30 Mar '04

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

1 0

Rotating molecules
by Eduardo Perozo 30 Mar '04

30 Mar '04

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

1 0

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

30 Mar '04

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

2 1

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

30 Mar '04

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

1 0

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

29 Mar '04

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

2 1

export pdb file
by Shuanghong Huo 29 Mar '04

29 Mar '04

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

1 0

question
by Wendy Ochoa 28 Mar '04

28 Mar '04

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

2 1

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

13 Mar '04

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

1 0

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

13 Mar '04

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

2 1

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

13 Mar '04

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

2 1

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

12 Mar '04

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

2 1

mac version of 1.1917 snapshot
by Eric Pettersen 11 Mar '04

11 Mar '04

Hi again, If you downloaded the Mac version of 1.1917 before about 5pm PST today, you should download it again because we updated it to fix a bug. The non-Mac versions are fine. The Mac version was not displaying atoms in stick, ball-and-stick, or sphere modes (wireframe mode was okay). Sorry for the possible inconvenience! Eric Pettersen UCSF Computer Graphics Lab pett(a)cgl.ucsf.edu

1 0

new snapshot release (1.1917)
by Eric Pettersen 10 Mar '04

10 Mar '04

Hi, We've just put a new snapshot release of Chimera up for download on the web site (http://www.cgl.ucsf.edu/chimera) Here are the highlights of the changes between this snapshot (1.1917) and the previous one (1.1892): tools: Match -> Align can create sequence alignments for multiple matched structures (not just pairs) Multalign Viewer can show/edit sequence annotations for relevant formats (RSF, Stockholm) SimpleSession -- saves label font/size and secondary structure assignments -- compiles to .pyc file: 50% less memory use, 15% less time to restore added Export Surfaces as VRML (in Utilities category) -- save surfaces made by Multiscale Models and Volume Viewer as VRML files added Browser Configuration (also in Utilities) -- configure browsers to send certain file types to Chimera (see also new Web Data preferences) FindHBond and PBReader allow line width to be set Movie has preliminary support for reading GROMOS trajectories newly available commands: bondcolor -- color bonds independent of flanking atoms bonddisplay -- set how bond display depends on flanking atom display bondrepr -- set bond representation to be wire or stick (independent of atoms) hbonds -- command equivalent of FindHBond tool rlabel -- whether residue name/number is included in the first atom label per residue Model Panel new 'render attrs...' button allows coloring by model/residue/atom attributes (e.g. B-factor) The 'render attrs...' interface is in preliminary form and will be changing. It also isn't documented yet but is pretty easy to use. Just pick an attribute of interest and a histogram of values will be displayed. By default there will be three color bars on the histogram that you can grab and move around with the mouse, or add or remove with control-click. You can use the 'Color' color well to set the color of the last bar you picked. The 'No value' color well is used in the rare case of the attribute being missing (or set to None) in the atom/residue/model. Click Apply or OK to make the coloring happen. Several parts of the interface are grayed out since they haven't been implemented yet. In the future you will be able to select by attribute value, and change atom radii or ribbon widths based on attribute values. If you are using the last production release (1.1872) instead of the previous snapshot, then you will also be getting these features when you go to 1.1917: general: additional file formats for saved images (PS, EPS) different lighting scheme in graphics window (over-the-shoulder instead of headlight) tools: added Shininess Control (in Viewing Parameters) -- adjust shininess, brightness, and the color used for shiny highlights Collaboratory section in Preferences allows designation of port number and IP address (for firewalls/NAT) Movie -- allows definition of scripts to be run at each time step -- looping can be turned off Multalign Viewer allows structure selection by arbitrary % residue conservation Volume Viewer options added to subdivide and smooth isosurfaces commands: copy command now offers supersampling and PS/EPS formats new ribbackbone command controls whether ribbon and backbone atoms are shown simultaneously Enjoy! Eric Pettersen UCSF Computer Graphics Lab pett(a)cgl.ucsf.edu http://www.cgl.ucsf.edu

1 0

Use of Chimera for Trends Plant Sci March 2004...
by John F. Allen 24 Feb '04

24 Feb '04

....for article and issue cover illustration. Final proofs attached. Acknowledged. Great program. Will send reprint(s) and final pdf on request. Many thanks John F. Allen Professor of Plant Cell Biology Plant Biochemistry Center for Chemistry and Chemical Engineering Lund University Box 124 SE-221 00 Lund Sweden Tel: +46 46 222 7788 (Office); +46 46 222 4073 (Laboratory) Fax: +46 46 222 4009 (Direct); +46 46 222 4116 (Departmental) E-mail: john.allen(a)plantbio.lu.se Plant Cell Biology: http://plantcell.lu.se Publications: http://plantcell.lu.se:591/john/

1 0

CHIMERA help - coloring density maps
by Padmaja Natarajan 02 Feb '04

02 Feb '04

Hi, I have just started using CHIMERA to look at volume files in the CCP4 format for virus particles. I would like to color the density as a color ramp continuously changing between ANY two specified colors depending on the radial distance from the center of the virus particle (as shown in the attached figure). How do I achieve this with CHIMERA? thank you very much for your help in advance. -Padma

2 1

Chimera dead after NVidia driver install...
by Patrick Ladam 12 Jan '04

12 Jan '04

Hi all, I had chimera 1.1872 working OK on my MDK 9.2 Linux box until I installed the NVidia driver for my GeFOrce2 card. What happens now is that chimera launches but stops during start-up with an error box indicating this message: "Togl: Couldn't get visual" This sounds like an OpenGL problem to me but I am far from being a specialist on this. Does anyone have a clue how I can debug/fix this ? Bye -- ----------------------------------------------------------------- | Patrick LADAM | | | Laboratoire CSSB | THE BIG BANG THEORY: | | UFR SMBH | | | 74 rue Marcel Cachin | In the begining there was | | 93017 Bobigny CEDEX | nothing at all. | | e-mail: | | | ladam(a)smbh.smbh.univ-paris13.fr | Then, it exploded... | | Tel: 01 48 38 77 26 / 76 85 | | | Fax: 01 48 38 77 77 | | -----------------------------------------------------------------

1 0

hi
by amutha ramaswamy 08 Jan '04

08 Jan '04

Dear Sir It is possible to make movies with PDB structurers in chimera. this would really be helpfull in visualing the conformatinal changes. I know chimera has the option for movie with trajectory files. but i wish to kow the option with PDB files only. Looking forward for your kind reply Sincerely Amutha Dr Amutha Ramaswamy Postdoctoral Fellow, Department of Biomedical Informatics 3184 Graves Hall, 333 West 10th Avenue The Ohio state University, Columbus, OHIO 43210 Phone: 614-688-3236 (Office); 614-291-0071(Residence) Fax: 614-292-7659 (Office) --------------------------------- Do you Yahoo!? Yahoo! Hotjobs: Enter the "Signing Bonus" Sweepstakes

2 1

Re: save selected atoms in a pdb file?
by Elaine Meng 08 Jan '04

08 Jan '04

Hi Sabuj, You can make surfaces that surround a set of atoms without making a separate PDB file. For example, if you have those atoms selected already, just use the command: surfcat mygroup sel where "mygroup" could be whatever you choose to name this set of atoms and "sel" means what is currently selected. Then, the command surface mygroup will show surface surrounding that group. There are several ways you could select the atoms of interest, and it sounds like you already figured out that part; the following detailed example is just one possibility: open pdb:2gbp sel :glc z<4.5 ~sel :glc ~sel solvent surfcat nearby sel ~sel surf nearby Here I opened a glucose-binding protein structure, selected all residues with any atom within 4.5 A of the ligand glucose (glc), then subtracted glucose itself and waters from the selection. Then I made the selected atoms a surface category with "surfcat," unselected them, and showed the corresponding surface. If you really want to save the selection to a file, it can be done on unix systems. Chimera can send information on what atoms are displayed to the system, and then the system command "cat" can send the information to a file. In Chimera, if the atoms of interest are selected: show sel pdbrun nouser cat > ~homedir/myfile.pdb (where "homedir" is your username) saves the displayed atoms to a PDB file named myfile.pdb in your home directory. It is necessary to use show (which undisplays everything else) since pdbrun includes all atoms that are displayed. The coordinates are transformed, so will be different from the original coordinates if you have moved the view around. If you want the untransformed coords, use the "reset" command before the "pdbrun" command. This process doesn't work on Windows as there is no cat command. I hope this helps, Elaine =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-= Elaine C. Meng, Ph.D. Computer Graphics Lab and Babbitt Lab meng(a)cgl.ucsf.edu http://www.cgl.ucsf.edu/home/meng/index.html =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=

1 0

save selected atoms in a pdb file?
by Pattanayek, Sabuj 08 Jan '04

08 Jan '04

Hi, Is it possible to save only the currently selected set of atoms/bonds or the atoms/bonds currently shown in the display into a pdb file? This is would be useful when doing surfaces for a particular location (e.g. atoms surrounding a central cavity), since the command "surface selected" seems to leave gaps in inaccessible areas (behind the atoms surrounding the cavity). Thanks, Sabuj Pattanayek

1 0

[Chimera-users] shininess/brightness/specular color of surfnet surfaces
by Pattanayek, Sabuj 05 Jan '04

05 Jan '04

Hi, Is it possible to alter the shininess/brightness/specular color of Surfnet surfaces? Also is it possible to save the surfaces so they don't have to be recreated every time? Thanks, ----------------------------------------------------------------- Sabuj Pattanayek Martin Egli Lab Vanderbilt University Dept. Biochemistry http://linegli.mc.vanderbilt.edu

3 2

Old version for Sparky
by laurent volpon 05 Jan '04

05 Jan '04

Hello, I would like to download an old version of Chimera (1.1442) because I have to use it with the same Python version (2.1) as Sparky. I tried at http://www.cgl.ucsf.edu/chimera/distrib/oldreleases.html but I had this: "Error: The URL you used is malformed". Please can you help me to find this version, or is there a way to use the version I have (1.1872) with the version 2.1 of Python (instead of 2.2). Thank you by advance for your answer. Sincerely Dr Laurent VOLPON ______________________________ Dr Laurent Volpon 3575 Avenue du Parc, Suite 5315 P.O. Box #85 Montréal, Qc, H2X 4A3 tel.: (514) 847-9494 ______________________________ --------------------------------- Lèche-vitrine ou lèche-écran ? Yahoo! Magasinage.

3 2

vdw viewer
by Terry Downing 23 Dec '03

23 Dec '03

Hey CGL, I would like to be able to visualize several different vdw parameter sets on the same molecule. Is there a way to visualize just the vdw surface? Can I change the parameter files easily? As always, thanks in advance for your help. Sincerely, Terry Downing -- P. Therese Downing Chemistry and Chemical Biology Graduate Program Kuntz and James Labs, UCSF Voice: (415)476-3986 Fax: (415)502-1411

1 0

ribbon+atoms
by Eric Pettersen 23 Dec '03

23 Dec '03

Hi, the 1892 snapshot release we just put out allows control over whether backbone atoms can be displayed with ribbons (on a per-model basis). There is a "ribbackbone" command, an entry in the New Molecules preference category, and a control in the attributes panel (reached via the "attributes..." button in the model panel). Eric Pettersen UCSF Computer Graphics Lab pett(a)cgl.ucsf.edu

1 0

interactive volume segmentation
by Juha Huiskonen 18 Dec '03

18 Dec '03

Hello, currently it is possible to select only box-shaped subregions in Volume Viewer. I would be especially interested in segmenting electron density maps or their isosurfaces to different subregions with any shape interactively. Is anyone planning to add a feature like this? Regards, Juha Huiskonen

2 3

Re: [Chimera-users] Chimera export to VRML?
by Thomas Goddard 17 Dec '03

17 Dec '03

Hi Morgan, I've written an extension to export Chimera multi-scale and volume surfaces as VRML. I wasn't able to export the molecular (MSMS) surfaces because the code does not provide access to the surface coloring. I've submitted a request to enhance the molecular surface code to make surface colors readable. I'm not sure if the VRML surface exporting will make it into the coming snapshot release. If it doesn't I will provide it to you by email. I'll let you know when we have the snapshot ready, probably later this week. Tom

1 0

Re: [Chimera-users] Chimera export to VRML?
by Thomas Goddard 16 Dec '03

16 Dec '03

Hi Morgan, We hope to have a snapshot Chimera release later this week. I'll try to get VRML output of multi-scale and msms surfaces in that release. Snapshot releases aren't as well tested as production releases, but are intended for trying out new features like this. Tom PS: Your email does not go to chimera-users unless you explicitly CC it to that list. We like to post to that list so others can benefit from it. I've included your email below and CC'ed this to chimera-users. Let me know if you want to keep your email private. > From: Morgan Ryan <mryan1(a)nc.rr.com> > Subject: Re: [Chimera-users] Chimera export to VRML? > Date: Tue, 16 Dec 2003 14:35:27 -0500 > > Thanks, Thomas. (Am I broadcasting to the list? I'm not sure if I'm > talking to a list or just you and Eric Petterson...) > > > For a large structure like a > > ribosome the VRML file ...might have a million triangles, and the file > > will > > at least be 30 megabytes. > > The 3D programs I use have fairly refined mesh management tools for > reducing complexity selectively, although being able to choose the mesh > density directly through MSMS would be pleasing. > > > You can choose the resolution of the multi-scale surfaces. We don't > > currently have VRML output for the multi-scale surfaces but that > > could be added fairly easily. I could look into that. > > The advantage of VRML output (or 3dmf, or dxf, or some other > interchange format) is that you could open the model in Maya, Studio > Max, Strata, Cinema 4D, any of the 3D animation programs, and make > startling movies wiith Hollywood special effects--light sabers! > Atomic-resolution recreations of enzyme reaction mechanisms! (The > latter is what I'm working on...) > > It would be a wonderful extension of Chimera if VRML output were added. > Morgan

1 0

Chimera export to VRML?
by Morgan Ryan 16 Dec '03

16 Dec '03

VRML seems to be a native internal format for (the magnificent) Chimera. Is it possible to export VRML models from Chimera? Or do you worldly folks know of a mechanism for converting msms output into vrml files? Thanks! Morgan Ryan, enthusiastic amateur ____________________ Morgan Ryan Development editor W H Freeman 209 Adams Way Chapel Hill NC 27516 phone: 919 942 4890 fax: 651 846 5265 mryan(a)whfreeman.com

3 3

help
by amutha ramaswamy 16 Dec '03

16 Dec '03

Dear sir I am using the Chimera molecular visualising software. I am not able to color the PDB structure according to the B-factors. could you plase let me know how do i color accroding to the bfactor. i would greatly apprecuate your help in this regard sincerely amutha Dr Amutha Ramaswamy Postdoctoral Fellow, Department of Biomedical Informatics 3184 Graves Hall, 333 West 10th Avenue The Ohio state University, Columbus, OHIO 43210 Phone: 614-688-3236 (Office); 614-291-0071(Residence) Fax: 614-292-7659 (Office) --------------------------------- Do you Yahoo!? New Yahoo! Photos - easier uploading and sharing

2 2

Re: Gaussians
by Thomas Goddard 11 Dec '03

11 Dec '03

Hi Lars, I'm not sure what you mean by visualizing some Gaussians. If you mean you want to make a volume which contains some specified Gaussians then I can give you a separate Python program that will do that. I already wrote it for another user. It is included at the end of this email. It outputs a netcdf format volume file. Tom % ./spots.py Make a volume array that is the sum of specified Gaussians. The first line of the input file has the grid cell size. Each subsequent line has 5 fields giving x, y, z, width, height. A volume grid is made that includes each Gaussian out to 3 widths. The resulting volume is written in netcdf format. Syntax: spots.py <infile> <outfile> % cat spotfile .5 15.49100 7.05000 17.16500 0.38 1.00 14.94600 7.55000 15.88800 0.35 1.00 16.01100 7.60400 14.81400 0.34 1.00 ... % ./spots.py spotfile gspottest.nc The spots.py code follows. Note that it contains a path to your Chimera installation because it needs Chimera modules to write the netcdf volume. #!/usr/bin/env python # ----------------------------------------------------------------------------- # Make a volume array that is the sum of specified Gaussians. # The resulting volume is written in netcdf format. # # Syntax: spots.py <infile> <outfile> # # The input file has a grid cell size on the first line followed by one # line per gaussian with 5 fields specifying x, y, z, width, height. # import os.path import sys # Need Chimera path to use VolumeData and Numeric modules CHIMERA_PATH = '/usr/local/chimera' sys.path.insert(0, os.path.join(CHIMERA_PATH, 'lib')) sys.path.insert(0, os.path.join(CHIMERA_PATH, 'share')) import Numeric # ----------------------------------------------------------------------------- # def make_spot_grid(path): step_size, spot_parameters = read_spots(path) range = 3 # in spot widths xyz_min, xyz_max = bounding_box(spot_parameters, range) xyz_origin = xyz_min xyz_step = (step_size, step_size, step_size) xyz_size = map(lambda a,b: a - b, xyz_max, xyz_min) import math grid_size = map(lambda s, step, math=math: 1 + int(math.ceil(s/step)), xyz_size, xyz_step) a = spot_array(spot_parameters, range, xyz_origin, step_size, grid_size) from VolumeData import Array_Grid_Data g = Array_Grid_Data(a, xyz_origin, xyz_step) return g # ----------------------------------------------------------------------------- # def read_spots(path): f = open(path, 'r') lines = f.readlines() f.close() step_size = float(lines[0]) spot_params = [] for line in lines[1:]: params = map(float, line.split()) spot_params.append(params) return step_size, spot_params # ----------------------------------------------------------------------------- # def bounding_box(spot_parameters, range): xmin = ymin = zmin = None xmax = ymax = zmax = None for x, y, z, w, h in spot_parameters: pad = range * w if xmin == None or x - pad < xmin: xmin = x - pad if ymin == None or y - pad < ymin: ymin = y - pad if zmin == None or z - pad < zmin: zmin = z - pad if xmax == None or x + pad > xmax: xmax = x + pad if ymax == None or y + pad > ymax: ymax = y + pad if zmax == None or z + pad > zmax: zmax = z + pad return (xmin, ymin, zmin), (xmax, ymax, zmax) # ----------------------------------------------------------------------------- # def spot_array(spot_parameters, range, xyz_origin, step_size, grid_size): shape = list(grid_size) shape.reverse() # array indices are z, y, x. array = Numeric.zeros(shape, Numeric.Float) for x, y, z, width, height in spot_parameters: add_spot((x,y,z), width, height, range, xyz_origin, step_size, array) return array # ----------------------------------------------------------------------------- # def add_spot(xyz, width, height, srange, xyz_origin, step_size, array): ijk = map(lambda a,b,s=step_size: (a - b)/s, xyz, xyz_origin) w = width / step_size w2 = w * w r = srange * width / step_size import math imin, jmin, kmin = map(lambda l, r=r, math=math: int(math.floor(l-r)), ijk) imax, jmax, kmax = map(lambda l, r=r, math=math: int(math.ceil(l+r)), ijk) for k in range(kmin, kmax+1): for j in range(jmin, jmax+1): for i in range(imin, imax+1): di, dj, dk = i-ijk[0], j-ijk[1], k-ijk[2] s2 = di*di + dj*dj + dk*dk e = height * math.exp(-s2/(2*w2)) array[k,j,i] += e # ----------------------------------------------------------------------------- # def bfactor_backbone_spots(path, molecule = None): if molecule == None: import chimera molecule = chimera.openModels.list()[0] out = open(path, 'w') out.write('.5\n') # step size import math for a in molecule.atoms: if a.name in ('C', 'N', 'CA') and hasattr(a, 'bfactor'): x, y, z = a.coord().xyz.data() w = math.sqrt(a.bfactor / (8 * math.pi**2)) h = 1 out.write('%8.5f %8.5f %8.5f %6.2f %6.2f\n' % (x,y,z,w,h)) out.close() # ----------------------------------------------------------------------------- # def syntax(): msg = ('Make a volume array that is the sum of specified Gaussians.\n' + 'The first line of the input file has the grid cell size.\n' + 'Each subsequent line has 5 fields giving x, y, z, width, height.\n' + 'A volume grid is made that includes each Gaussian out to 3 widths.\n' + 'The resulting volume is written in netcdf format.\n' + 'Syntax: spots.py <infile> <outfile>\n') sys.stderr.write(msg) sys.exit(1) # ----------------------------------------------------------------------------- # def run_command(): if len(sys.argv) != 3: syntax() inpath = sys.argv[1] outpath = sys.argv[2] grid = make_spot_grid(inpath) from VolumeData.netcdf.netcdf_grid import write_grid_as_netcdf write_grid_as_netcdf(grid, outpath) # ----------------------------------------------------------------------------- # if __name__ == '__main__': run_command()

1 0

B-Factor
by Luis Cuello 09 Dec '03

09 Dec '03

Hello, My name is Luis Cuello and I am a graduate student at the University of Virginia..I started using chimera recently but I have not been able to colored my files by B-factors values, Is there any way to do that with chimera software??? If you could help me in this matter I would appreciated..Sincerely Luis Cuello University of Virginia Molecular Physiology 1(434)9245473 e-mail:lgc4r@virginia.edu

2 1

Re: [Chimera-users] delphi potenital map viewing
by Thomas Goddard 05 Dec '03

05 Dec '03

Hi Klaas, The Delphi Viewer is the only tool that can color a molecular surface according to the potential value. Can you email me the Python traceback you get when the "_cubes error: ..." occurs? This shows exactly where in the code the error happened. You find it by looking at Favorites/Reply Log. Having the traceback alone might be enough for me to figure out and correct the problem. If you could put the phi file on the web so I could download it I would almost certainly be able to figure out the problem. The Delphi Viewer is supposed to be able handle maps of any size. It will be much slower than volume viewer, but it should not generate an error with a large map. Tom > Date: Fri, 05 Dec 2003 21:10:20 +0100 > From: Klaas Decanniere <Klaas.Decanniere(a)pandora.be> > Reply-To: Klaas.Decanniere(a)pandora.be > To: Thomas Goddard <goddard(a)cgl.ucsf.edu> > Subject: Re: [Chimera-users] delphi potenital map viewing > > Thomas Goddard wrote: > > >Hi Klaas, > > > > The size of the Delphi potential map Chimera can display is only > >limited by the size of your computer memory. > > > > If you open a large map it does not get shown immediately if the map > >is larger than 1 Mvoxels (~100^3). > > > > > > When I tried using the "delphi viewer" option, I get _cubes error: grid > index out of range > Using the volume viewer I can indeed load everything. > Is there a way to tie the volume viewer thing to a molecular surface? > With a colour gradient showing different levels of potential? (As in the > Molecular surface option in the Delphiviewer) > Or maybe I can save the subsampled map somehow and then read it into the > delphiviewer thing? > > Thanks again for your help. > > Klaas > > P.S.: I like the program more each day. > P.P.S.: If you ever have a day with nothing to do: the automatic > fetching of pdb files doesn't work for me as I'm behind a proxy server.

1 0

delphi potenital map viewing
by Klaas Decanniere 04 Dec '03

04 Dec '03

Hello, Thinking that more is better I ended up with Delphi potential maps on a rather large grid which I would like to visualize with chimera (gsize=213) Chimera seems to think the grid is to large. Is there an easy way to have chimera accept bigger grids (speed is not really an issue)? Alternatively, what is the largest grid that chimera can accept? Thanks for your help. Klaas Decanniere

2 1

ribbon+atoms
by Elaine Meng 25 Nov '03

25 Nov '03

Dear Soren, Thanks for the compliments - we're trying! You're right that in recent versions of Chimera , showing ribbon turns off the backbone atoms. There is an easy workaround: open two copies of the same molecule, use one copy to show ribbon and the other copy to show atoms. The reason for the behavior is when the ribbon is drawn as a smooth curve, it does not exactly pass through the backbone atoms. Often one wants to show just the ribbon and certain side chains, and the current behavior is to extend the side chain to meet the ribbon (even when the backbone atoms are not on the ribbon) so that they are not floating in space. However, we agree that sometimes one will want to display both the ribbon and the backbone atoms, and in the future there will a preference to switch between: (a) show ribbon, not backbone atoms, and extend the side chains to meet the ribbon (b) show ribbon and whatever backbone atoms are displayed by the user, do not extend side chains to meet the ribbon Until this switchable preference becomes available, we hope that the suggested workaround will let you do what you want. Sincerely, Elaine =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-= Elaine C. Meng, Ph.D. Babbitt Lab and Computer Graphics Lab meng(a)cgl.ucsf.edu http://www.cgl.ucsf.edu/home/meng/index.html =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=

1 0

ribbon+atoms
by Soren M. Kristensen 25 Nov '03

25 Nov '03

Dear Chimera developers You have createdt an absolutely brilliant program. I have a small question/suggestion about displaying atoms and ribbons simultaneously. When displaying a number of backbone atoms and then putting on a ribbon, then some of the atoms dissapear in release 1700 and onwards. Is there a way to avoid this behaviour? It basically means that I am stuck with the pre-1700 version, but maybe there is some parameter I can tweak to make it work in the new versions? I would be grateful for your help. I reproted this as a bug some 6 months ago, but was told that it was supposed to be like this. I measure experimental values for selected backbone atoms and display them as colorcoded spheres on top of a ribbon to indicate what is measured different places in the structure. If you have suggestions to acieve this behaviour in another way I would be grateful. Thank you in anticipation, Søren -- Søren M. Kristensen, Ph.D. Phone: (+45) 3532 0301 University of Copenhagen FAX: (+45) 3535 0609 Department of Chemistry E-mail: smk(a)kl5.ki.ku.dk H. C. Ørsted Institutet http://kl5.ki.ku.dk/~smk Universitetsparken 5 Room: C508 DK-2100 København Ø Denmark

1 0