March 2008 - Chimera-users - RBVI Mailing Lists

turning off PipesAndPlanks or Axes
by Boaz Shaanan 21 Mar '08

21 Mar '08

Hi, How do I turn off PipesAndPlanks once they're depicted ? The same holds for Axes. (build 2492). Thanks, Boaz Boaz Shaanan, Ph.D. Dept. of Life Sciences Ben-Gurion University of the Negev Beer-Sheva 84105 Israel Phone: 972-8-647-2220 ; Fax: 646-1710 Skype: boaz.shaanan‎

2 1

source headers, release
by Randy Heiland 20 Mar '08

20 Mar '08

Before I attempt to rebuild a plugin, can someone tell me if the 1.2349 headers that you currently provide: http://www.cgl.ucsf.edu/ chimera/sourcecode.html will play nice with the current, 1.2470, release? thanks, Randy

2 1

problem with 17/03 build linux 64
by Jean-Didier Maréchal 19 Mar '08

19 Mar '08

Hi guys, thanks a lot for the different answers. I wanted to test Elaine's suggestion and I install the last linux 64 build. When running chimera I have: Traceback (most recent call last): File "/usr/local/chimera/share/__main__.py", line 65, in <module> value = chimeraInit.init(sys.argv) File "/usr/local/chimera/share/chimeraInit.py", line 304, in init from chimera import tkgui File "CHIMERA/share/chimera/tkgui.py", line 80, in <module> File "CHIMERA/share/chimera/printer.py", line 444, in <module> AttributeError: 'module' object has no attribute 'No' Some problem in my installation ? All the best JD El mié, 19-03-2008 a las 12:00 -0700, chimera-users-request(a)cgl.ucsf.edu escribió: > 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: Principal axis of inertia (Eric Pettersen) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Wed, 19 Mar 2008 11:47:51 -0700 > From: Eric Pettersen <pett(a)cgl.ucsf.edu> > Subject: Re: [Chimera-users] Principal axis of inertia > To: "'Chimera BB'" <chimera-users(a)cgl.ucsf.edu> > Message-ID: <DDE1C65B-AD42-48C1-9414-91625662D812(a)cgl.ucsf.edu> > Content-Type: text/plain; charset=ISO-8859-1; delsp=yes; format=flowed > > Just wondering why you weighted by atomic number rather than mass > (a.element.mass)? > > --Eric > > On Mar 19, 2008, at 11:42 AM, Tom Goddard wrote: > > > Hi JD, > > > > Attached is a script that computes the inertia axes of molecules. > > Open a molecule, then open the script. It will print the axes in > > the reply log and show an inertia ellipsoid. It is weighting the > > atoms by their atomic number, not by the average isotopic mass. > > Tested in Chimera 1.2498. Should work in older Chimera versions. > > > > Tom > > > > Jean-Didier Mar?chal wrote: > >> Hi everyone, > >> > >> I was wondering if there is a method available to calculate the > >> principal axis of inertia of a protein (or a part of a protein) in > >> chimera? > >> > >> If not what would be the easiest way to move forward? > >> > >> All the best, > >> > >> JD > >> > >> > > > > > > # > > ---------------------------------------------------------------------- > > ------- > > # Compute inertia tensor principle axes for molecule. > > # > > def inertia_ellipsoid(m): > > > > atoms = m.atoms > > n = len(atoms) > > from _multiscale import get_atom_coordinates > > xyz = get_atom_coordinates(atoms) > > anum = [a.element.number for a in atoms] > > from numpy import array, dot, outer, argsort, linalg > > wxyz = array(anum).reshape((n,1)) * xyz > > mass = sum(anum) > > > > c = wxyz.sum(axis = 0) / mass # Center of mass > > v33 = dot(xyz.transpose(), wxyz) / mass - outer(c,c) > > eval, evect = linalg.eigh(v33) > > > > # Sort by eigenvalue size. > > order = argsort(eval) > > seval = eval[order] > > sevect = evect[:,order] > > > > return sevect, seval, c > > > > # > > ---------------------------------------------------------------------- > > ------- > > # > > def ellipsoid_surface(center, axes, lengths, color = (.7,.7,.7,1)): > > > > from Icosahedron import icosahedron_triangulation > > varray, tarray = icosahedron_triangulation(subdivision_levels = 3, > > sphere_factor = 1.0) > > from numpy import dot, multiply > > es = dot(varray, axes) > > ee = multiply(es, lengths) > > ev = dot(ee, axes.transpose()) > > ev += center > > > > import _surface > > sm = _surface.SurfaceModel() > > sm.addPiece(ev, tarray, color) > > return sm > > > > # > > ---------------------------------------------------------------------- > > ------- > > # > > def show_ellipsoid(axes, d2, center, model): > > > > from math import sqrt > > d = [sqrt(e) for e in d2] > > sm = ellipsoid_surface(center, axes, d) > > sm.name = 'inertia ellipsoid for %s' % m.name > > from chimera import openModels as om > > om.add([sm], sameAs = model) > > > > # > > ---------------------------------------------------------------------- > > ------- > > # > > def print_axes(axes, d2, m): > > > > from math import sqrt > > paxes = ['\tv%d = %6.3f %6.3f %6.3f d%d = %6.3f' % > > (a+1, axes[a][0], axes[a][1], axes[a][2], a+1, sqrt(d2[a])) > > for a in range(3)] > > from chimera.replyobj import info > > info('Inertia axes for %s\n%s\n' % (m.name, '\n'.join(paxes))) > > from Accelerators.standard_accelerators import show_reply_log > > show_reply_log() > > > > # > > ---------------------------------------------------------------------- > > ------- > > # > > from chimera import openModels as om, Molecule > > mlist = om.list(modelTypes = [Molecule]) > > for m in mlist: > > axes, d2, center = inertia_ellipsoid(m) > > print_axes(axes, d2, m) > > show_ellipsoid(axes, d2, center, m) > > _______________________________________________ > > 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 59, Issue 18 > *********************************************

2 1

Un-tile
by Jim Nettles 19 Mar '08

19 Mar '08

Thanks that is a useful function I hadn't explored. However, how to you un-tile to go back to a previously aligned system? Thanks, Jim > > Message: 4 > Date: Wed, 19 Mar 2008 16:53:11 +0530 > From: "Bharat Yadav" <bharatbioinformaticshau(a)gmail.com> > Subject: [Chimera-users] Requesting to solve my poblem > To: chimera-users(a)cgl.ucsf.edu > Message-ID: > <e162fe4e0803190423k5651193amf41e192251b613c3(a)mail.gmail.com> > Content-Type: text/plain; charset=ISO-8859-1 > > Hello Sir / Madam > > Good Evening > > My name is Bharat Yadav and I am perusing Master's Degree in > Bioinformatics from Chaudhary Charan Singh Haryana Agriculrual > University - Hisar (India). I am working here on Molecular Modeling. > With a long time I am using UCSF Chimera. It is the one of wonderful > Molecular visiualization software haivng some very peculear features. > Respected Sir / Madam, In my present work I am facing a bit problem. > Hoping form you to help regarding this. My problem is as follows. I > have generated five models of protein. when I open all these modles > these were overlapping each other. I want to view all these five > models in a single panel non overlapping. I mean I want to view these > 3D models in panel as follows > ----------------------------------------- > | | | > | | | > | Model 1 | Model 2 | > | | | > ----------------------------------------- > | | | > | | | > | Model 3 | Model 4 | > | | | > ----------------------------------------- > | | | > | | | > | Model 5 | Model 6 | > | | | > ----------------------------------------- > > Hope that you now easily understand what actually I want. Please Help > me regarding this and mail me its solution or steps or commands for > this task. Please free to ask any question. > > Thanks in anticipation > > Warm regards from Bharat Yadav > > > ------------------------------ > > Message: 5 > Date: Wed, 19 Mar 2008 10:11:52 -0700 > From: Elaine Meng <meng(a)cgl.ucsf.edu> > Subject: Re: [Chimera-users] Requesting to solve my poblem > To: "Bharat Yadav" <bharatbioinformaticshau(a)gmail.com> > Cc: chimera-users(a)cgl.ucsf.edu > Message-ID: <37BB555C-E967-43BB-B6EA-5CD8A75853FD(a)cgl.ucsf.edu> > Content-Type: text/plain; charset=US-ASCII; delsp=yes; format=flowed > > Hi Bharat, > > Thanks for your kind words! I think you just need to use "tile" in > the Model Panel. Open Model Panel (Favorites... Model Panel). In > there, choose all 5 models on the left side. Then, on the right, > change the menu from "frequently used" to "infrequently used" - now > there should be a "tile" button on the right - try using that. > > 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 > >

2 2

Data of Slice display
by Jason 19 Mar '08

19 Mar '08

Hi all I am using slice display to plot the profile of density map. That is the function that I want to exam the profile of virus density map without spherical average. Where can I get the data of each dot of profile, so I can plot profile at Excel for further use. Thank you for your help Kind regards Jason Chen -- Open WebMail Project (http://openwebmail.org)

2 1

Average slab thickness
by Tom Goddard 19 Mar '08

19 Mar '08

Hi Jeff, The Chimera mask command code can be adapted to calculate the average thickness of a slab defined by two surfaces. I've attached the adapted Python code. You select the surfaces then open the script and it reports the mean thickness, standard deviation of thickness, min/max thickness and number of thickness samples used. The script contains a projection axis which should be perpendicular to the slab and a grid step which defines the grid spacing. The thickness is measured on a grid perpendicular to the axis and only grid points where the projection axis intersects the surfaces at least twice are used. You have to edit the script (at the bottom) to set the projection axis and grid step to appropriate values. Those values could be automatically set using the Python code of the inertia ellipsoid script I gave you earlier. Tested in Chimera version 1.2498 on a spherical surface. Tom Jeff wrote: > hi Tom, > > I'd like to compute the 'average distance' between two surfaces to > calculate the width of a cellular compartment (have attached an image, > the two surfaces delineating the compartment are colored blue). > > one way that has been suggested to me is to calculate the volume > between the two surfaces, and then just divide by a mean path length > along one of the surfaces. I know that chimera has a volume > calculation tool for an isosurface, and was wondering if I could use > the surfaces you see to do a volume calculation (and, if so, how I > would go about capping those surfaces to create a closed surface). If > able to calculate the volume that way, I could just use the path > length of a contour (which I have) that generated one of the two > surfaces to derive the "average width". > # ----------------------------------------------------------------------------- # Script to calculate average thickness of a slab defined by surfaces. # Adapted from mask command code for Jeff Triffo, March 18, 2008. # # Adapted version of masked_volume() from depthmask.py. # # Only uses first two surface layers along projection axis. # def average_thickness(projection_axis = (0,0,1), grid_step = 1.0, surfaces = None): # Use selected surfaces if none are passed in. if surfaces is None: from Surface import selected_surface_pieces plist = selected_surface_pieces() from Mask.depthmask import surface_geometry surfaces = surface_geometry(plist, tf = None, pad = None) # Determine transform from vertex coordinates to depth array indices from Mask.depthmask import orthonormal_frame fx,fy,fz = orthonormal_frame(projection_axis) from numpy import array, float32, intc, zeros tf = array(((fx[0], fx[1], fx[2], 0), (fy[0], fy[1], fy[2], 0), (fz[0], fz[1], fz[2], 0)), float32) / grid_step # Transform vertices to depth array coordinates. zsurf = [] tcount = 0 for vertices, triangles in surfaces: varray = vertices.copy() from Mask.depthmask import apply_transform apply_transform(tf, varray) zsurf.append((varray, triangles)) tcount += len(triangles) if tcount == 0: return 0,0,0,0,0 from Mask.depthmask import bounding_box vmin, vmax = bounding_box(zsurf) voffset = -(vmin - 0.5) tf[:,3] += voffset from _contour import shift_vertices for varray, triangles in zsurf: shift_vertices(varray, voffset) from math import ceil, floor dxsize = int(ceil(vmax[0] - vmin[0] + 1)) dysize = int(ceil(vmax[1] - vmin[1] + 1)) # Create depth arrays depth = zeros((dysize,dxsize), float32) tnum = zeros((dysize,dxsize), intc) depth2 = zeros((dysize,dxsize), float32) tnum2 = zeros((dysize,dxsize), intc) # Copy volume to masked volume at masked depth intervals. max_depth = 1e37 depth.fill(max_depth) tnum.fill(-1) from Mask.depthmask import surfaces_z_depth any = surfaces_z_depth(zsurf, depth, tnum, None, None) if not any: return 0,0,0,0,0 depth2.fill(max_depth) tnum2.fill(-1) surfaces_z_depth(zsurf, depth2, tnum2, depth, tnum) # Find thicknesses at grid points that intercept surfaces twice. dlist = [] for i in range(dysize): for j in range(dxsize): if tnum[i,j] != -1 and tnum2[i,j] != -1: d = (depth2[i,j] - depth[i,j]) * grid_step dlist.append(d) if len(dlist) == 0: return 0,0,0,0,0 from numpy import mean, std return mean(dlist), std(dlist), min(dlist), max(dlist), len(dlist) # ----------------------------------------------------------------------------- # projection_axis = (0,0,1) grid_step = 1.0 m, sd, mn, mx, n = average_thickness(projection_axis, grid_step) print ('Slab thickness along axis %s sampled at grid spacing %.3g' % (str(projection_axis), grid_step)) print 'Mean %.5g, sdev %.5g, min %.5g, max %.5g, samples %d' % (m,sd,mn,mx,n) from Accelerators.standard_accelerators import show_reply_log show_reply_log()

2 1

Requesting to solve my poblem
by Bharat Yadav 19 Mar '08

19 Mar '08

Hello Sir / Madam Good Evening My name is Bharat Yadav and I am perusing Master's Degree in Bioinformatics from Chaudhary Charan Singh Haryana Agriculrual University - Hisar (India). I am working here on Molecular Modeling. With a long time I am using UCSF Chimera. It is the one of wonderful Molecular visiualization software haivng some very peculear features. Respected Sir / Madam, In my present work I am facing a bit problem. Hoping form you to help regarding this. My problem is as follows. I have generated five models of protein. when I open all these modles these were overlapping each other. I want to view all these five models in a single panel non overlapping. I mean I want to view these 3D models in panel as follows ----------------------------------------- | | | | | | | Model 1 | Model 2 | | | | ----------------------------------------- | | | | | | | Model 3 | Model 4 | | | | ----------------------------------------- | | | | | | | Model 5 | Model 6 | | | | ----------------------------------------- Hope that you now easily understand what actually I want. Please Help me regarding this and mail me its solution or steps or commands for this task. Please free to ask any question. Thanks in anticipation Warm regards from Bharat Yadav

2 1

interface residues
by Francesco Malatesta 17 Mar '08

17 Mar '08

Hello, could anyone help me out with this problem. I would like to exclusively colour interface residues in a dimeric protein represented as spheres or surface. Is this possible? Thank you very much Francesco Malatesta

2 1

preferences - povray - keep input files doesn't
by Klaas Decanniere 17 Mar '08

17 Mar '08

For at least some versions of Chimera, the "keep pov-ray input files" option in the Favorites - Preferences- Pov-ray tab does not work ( .pov files are deleted after rendering during animation) This is true e.g for build 2491 2008/02/11 (Platform: linux2 windowing system x11). Is there a way to force this option? I would like to do some post processing on the .pov files before rendering. An option to generate an animation series .pov files without doing the actual rendering would be even better - how can I do that? Thanks for the help and kind regards, Klaas Decanniere

4 4

segmentation fault - surface rendering
by Michael DiMattia 13 Mar '08

13 Mar '08

Hello, I am trying to surface render a few different proteins within the same Chimera session and I keep encountering a segmentation faullt. I imagine this is due to a memory limit being reached, although I do not know if it is a Chimera parameter that can be increased/overriden. Instead of having the different macromolecules in the same pdb, I have actually split it up into many pdbs and opened them one at a time, surface rendering one at a time. This worked for awhile but it seems I can't even get that far now. Any advice or tips on what I can do? I realize surface rendering is very computationally intensive but I would love to get this to work! Thanks Mike DiMattia

3 2

Chimera-users March 2008