- Chimera-users - RBVI Mailing Lists

Calculating SAS and SES with Chimera-MSMS wrap
by Jaime Rodríguez-Guerra 07 Apr '14

07 Apr '14

Dear Chimera dev team, I am trying to use your MSMS wrapper to compute some SAS/SES values for a given set of atoms. As suggested in a previous post in this list, I've taken a look into Measure.measure.buried_area() and I think I've got it, but just wanted to make sure. Here is my current procedure: ``` import chimera, Measure, MoleculeSurface atoms = chimera.selection.currentAtoms() xyzr_data = Measure.measure.atom_xyzr(atoms) surfaces = MoleculeSurface.xyzr_surface_geometry(xyzr_data) for i, a in enumerate(atoms): print "Atom: {0}, SAS: {1}, SES: {2}".format(a.name, surfaces[3][i,0], surfaces[3][i,1]) ``` Is this OK? I don't know how to handle all that output that MSMS is giving back! Thanks in advance, Jaime.

2 1

Re: [Chimera-users] Reaction Mechanism Movie
by Eric Pettersen 07 Apr '14

07 Apr '14

On Apr 5, 2014, at 1:06 PM, Moses Adenaike wrote: > Dear Dr. Pettersen, > I am wondering how I would go about making a moive detailing a biochemical reaction. For example, I want to show a DNA polymerase reaction, essentially a ligation between the backbone of two nucleotides. How would I go about doing so? I know chimera has a make new bond function but I haven't been too successful at getting that to work. Hi Moses, It is typically better to send questions like this to the chimera-users mailing list than to me directly so that others can benefit from the answer and in case others on the list have more expertise on the topic than I do. The simplest method is to decide which frame# you want the bond to start existing and use the "Define script…" entry in MD Movie's Per-Frame menu to define a one-command script to execute at that frame. For example, if you wanted a bond between :337.A@P and :338.A@O3' to exist starting with frame 53, in the Per-Frame dialog you would set "Interpret script as" to "Chimera commands" and enter as your script: #53 bond :337.A@P :338.A@O3' There is also a hidden capability in Chimera to compute all bonds on a per-frame basis. I can tell you how to access that if it turns out you really need it. However, for the purpose of making a movie it might be nice to have a more "continuous" representation of the bond forming than just the binary on/off that the above approaches give you. I have attached a Python script that you can use the Per-Frame dialog (with "Interpret script as" set to Python, and reading in the script with the "Insert text file…" button) to have a bond form (by starting out very thin and becoming thicker) as two designated atoms approach each other. There are some values at the top of the script that you can edit to make it work best in your situation. There is a distance value for when the bond begins to exist (noBond, default 3.5) and when it reaches full size (fullBond, default 2.0). There is a two-member list of residue IDs (resIDs) for the atoms involved in the bond, and a list of atom names (atomNames). Let me know if you need more help than this. Hmmm, this all assumes you have a trajectory of the reaction. You may have to use the Morph Conformations tool to create such a trajectory if you don't have one already. --Eric Eric Pettersen UCSF Computer Graphics Lab http://www.cgl.ucsf.edu

1 0

presets-interactive 3-hydrophobicity surface- those water molecules stay!
by Ying-Chieh Sun 06 Apr '14

06 Apr '14

Hi, I did a simulation using Amber package and saved a snapshot into its pdb file, and use Chimera to open it. The pdb file contains a protein-ligand complex plus a box of about 20000 TIP3P water molecules. When I pressed presets-interactive 3-hydrophobicity surface, there are 8 water molecules stayed close to the ligand binding site. I am wondering how these 8 water molecules were determined? The reason I ask is because we are trying to find a way to determine 'crystal (or say, ordered) water molecules'. This presets-interactive 3-hydrophobicity surface seems to perform this function in some way. Thank you for your help. Ying-chieh Ying-Chieh Sun, Professor of Computational Chemistry Chemistry Department National Taiwan Normal University Taipei, Taiwan 11617 Tel: +886-2-7734-6215 Fax: +886-2-2932-4249

3 2

Resolution of .stl file of ribbon representatiom
by Gavin Whittaker 04 Apr '14

04 Apr '14

Dear all, I'm in the process of 3d printing a large protein in ribbon format; I've generated the .stl file, but the resolution of the mesh that Chimera creates means that the resulting file is too large to use easily. I've managed to use other software to reduce the resolution of the mesh, so I /can/ now get around the problem, but can anyone tell me if it's possible to reduce the resolution of the mesh in Chimera to save me a step for future models? I know how to reduce the resolution of a surface mesh, but files of ribbon representations remain obstinately large. One method I've used that gives limited results is to simplify the cross-sectional form of the ribbon, but I wonder, am I missing something obvious and more effective? thanks! Gavin

3 7

query
by Bhaskar Bagchi 04 Apr '14

04 Apr '14

Respected Sir I want to calculate surface area of different atoms in a molecule. Is it possible in chimera? Please help me. Thanking you Yours faithfully, Bhaskar Bagchi Research Scholar Raiganj University College Email: bhaskarbagchi(a)gmail.com

2 1

ellipsoid main axis alignment ?
by Hugo J Bohorquez 03 Apr '14

03 Apr '14

Hi, I am writing a script in python for generating a series of ellipses linking the nuclei a molecule. Every ellipsoid should be positioned at a specific location *P* = (X,Y,Z) and should be aligned with the orthonormal vectors describing its main directions-*A*, *B*, *C*. However, when I create an ellipsoid with shape_command it gets located correctly at *P* but it is aligned with the global axis. I can't find a general way to reorient an object so its original axes point along the vectors *A*, *B*, and *C*. I will appreciate any suggestion how to do it. -- Hugo J Bohórquez <https://sites.google.com/site/hugojbohorquez/home>, PhD *Investigador* Grupo de Biofísica Fundación Instituto de Inmunología de Colombia TEL. +57 +1 -3158920 Ext. 122 CEL. +57 315-6169069 *The Pauli energy grows exponentially with the electronic localisation <http://dx.doi.org/10.1080/00268976.2013.836605>*

3 3

resampling maps
by Jason Kaelber 03 Apr '14

03 Apr '14

Hi all, I decided to do a rough alignment of density maps manually and resample onto a common coordinate system to feed the output mrc files into my fine aligner. But, when I vop resample #1 onGrid #0, I get, TypeError: Cannot cast ufunc add output from dtype('float32') to dtype('uint8') with casting rule 'same_kind' What might I be doing wrong to get this error? As far as I know, all the maps should be uint8 since they were most recently operated on by David Mastronarde's newstack command which outputs them as unsigned bytes. Thanks, Jason <type 'exceptions.TypeError'> Exception in Tk callback Function: <bound method MidasUI.processCommand of <Midas.midas_ui.MidasUI instance at 0x276b2d88>> (type: <type 'instancemethod'>) Args: (<Tkinter.Event instance at 0x276a2ea8>,) Event type: KeyPress (type num: 2) Traceback (innermost last): File "/opt/UCSF/Chimera64-1.8.1/lib/python2.7/site-packages/Pmw/Pmw_1_3/lib/PmwBase.py", line 1747, in __call__ return apply(self.func, args) File "/opt/UCSF/Chimera64-1.8.1/share/Midas/midas_ui.py", line 283, in processCommand midas_text.makeCommand(cmdText) File "/opt/UCSF/Chimera64-1.8.1/share/Midas/midas_text.py", line 69, in makeCommand f(c, args) File "/opt/UCSF/Chimera64-1.8.1/share/VolumeFilter/ChimeraExtension.py", line 24, in vop_cmd vop_command(cmdname, args) File "/opt/UCSF/Chimera64-1.8.1/share/VolumeFilter/vopcommand.py", line 91, in vop_command doExtensionFunc(func, fargs, specInfo = spec) File "/opt/UCSF/Chimera64-1.8.1/share/Midas/midas_text.py", line 435, in doExtensionFunc extFunc(*tuple(processedArgs), **kw) File "/opt/UCSF/Chimera64-1.8.1/share/VolumeFilter/vopcommand.py", line 553, in resample_op boundingGrid, False, None, modelId) File "/opt/UCSF/Chimera64-1.8.1/share/VolumeFilter/vopcommand.py", line 168, in add_operation scale = scale[i]) File "/opt/UCSF/Chimera64-1.8.1/share/VolumeViewer/volume.py", line 1535, in add_interpolated_values m[:,:,:] += values <type 'exceptions.TypeError'>: Cannot cast ufunc add output from dtype('float32') to dtype('uint8') with casting rule 'same_kind' ================================================ Event contents: char: delta: 36 height: ?? keycode: 36 keysym: Return keysym_num: 65293 num: ?? send_event: False serial: 39103 state: 16 time: 139539396 type: 2 widget: .59791048.661335648.802317088.802317376.802317448 width: ?? x: 159 x_root: 524 y: 9 y_root: 988 TypeError: Cannot cast ufunc add output from dtype('float32') to dtype('uint8') with casting rule 'same_kind' File "/opt/UCSF/Chimera64-1.8.1/share/VolumeViewer/volume.py", line 1535, in add_interpolated_values m[:,:,:] += values See reply log for Python traceback. -- National Center for Macromolecular Imaging, N420 Baylor College of Medicine, MS BCM125 1 Baylor Plaza Houston, TX, USA

2 3

Is there a way to save animation editing process?
by Li, Leon 03 Apr '14

03 Apr '14

Dear all, I'm try to save all my editing of animation, it seems that 'save session' can't save 'scenes', 'timelines', etc. I wonder is there a way to save all the process, so that I can re-open it to edit. Thanks, Leon

2 2

combining different chains
by faeze mohajer 01 Apr '14

01 Apr '14

good day I have modeled heavy and light chains of myosin now I want to combine them through chimera, would you plz help me on that your faithfully *FAEZE SADAT MOHAJERPhD Student * *Faculty of Bioscience and Medical Engineering (FBME)* *Universiti Teknologi Malaysia (UTM)* Buzzdock ads <http://www.buzzdock.com/setting/V2hpbG9raWk=> Search ResultsFaculty of Biosciences & Medical Engineering<http://fbme.utm.my/>

2 1

Re: [Chimera-users] Chimera Shutting Down
by Tom Goddard 31 Mar '14

31 Mar '14

Hi Joshua, Glad to hear that the new graphics driver fixed the Chimera side view crash on your Ubuntu 12.04 system. Maybe this will help some other Chimera users. Tom On Mar 31, 2014, at 12:33 PM, Joshua Broyde wrote: > Dear Tom, > I installed a new graphics driver, so now it works. In case anybody else has this same problem, I installed the Gallium Radeon driver and it worked fine. I followed step 2 from this link: > http://arstechnica.com/civis/viewtopic.php?f=16&t=1137372 > > Thanks for the tip. > Sincerely, > Joshua Broyde > ----- Original Message ----- > From: "Tom Goddard" <goddard(a)sonic.net> > To: "Joshua Broyde" <jb3401(a)c2b2.columbia.edu> > Cc: chimera-users(a)cgl.ucsf.edu > Sent: Thursday, March 27, 2014 12:58:50 AM > Subject: Re: [Chimera-users] Chimera Shutting Down > > The crash is due to a buggy graphics driver and the only fix is likely to be installing a working graphics driver. Sometimes some trick like starting the side view before opening models as described in the post you referenced avoids the problem, but usually there is no way to avoid a graphics driver crash other than replace the driver with one that doesn’t have the bug. > > Tom > > > On Mar 26, 2014, at 6:05 AM, Joshua Broyde wrote: > >> Hi, >> I am using the Chimera Linux 64-bit (chimera-alpha-linux_x86_64.bin), and Chimera shuts down when I try to open the Side View option. I tried the workaround mentioned here: http://plato.cgl.ucsf.edu/pipermail/chimera-users/2013-April/008765.html but it is not working for me. >> >> OS: Ubuntu 12.04 >> Version of Chimera: Linux 64-bit Daily Build (Mar 24, 2014) >> >> >> For what it's worth, on another computer where I am running Ubuntu 12.10, with Chimera Linux 64-bit Daily Build (July 10, 2013), I did not have this problem. >> Is there another workaround that I could do to open the side view? >> >> Thanks! >> >> Sincerely, >> Joshua Broyde >> _______________________________________________ >> Chimera-users mailing list >> Chimera-users(a)cgl.ucsf.edu >> http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users >> >

1 0