May 2013 - Chimera-users - RBVI Mailing Lists

22 May '13

Dear Chimera, I want to know something very simple which is making my life very difficult at the moment. I have an EM map of a virus capsid and I would like to make various distance measurements. I simply want to draw a strait line and measure it or place two points on the map and measure the distance between them. I have read all your help files on this matter and I cannot make any sense out of the instructions Can you help me out ? Thankyou Mr A.I

2 1

Fitmap in batch mode
by ingvar 22 May '13

22 May '13

Hello, I am trying to run fitmap on a set of volumes on a server without a display. I would like to collect the results of the top hits to store in a database – correlation coefficients, inclusion scores, and the transforms. The only way I found to collect the fitting results is to save the session after displaying the listFits dialog. The listFits dialog cannot be displayed when running Chimera in –nogui mode, otherwise this would have been a solution. I tried to use a virtual display, with Xvfb. I have used Xvfb with some java applications in the past. When I use Xvfb with Chimera, it appears that Chimera do not want to start up properly, the process just hangs, i.e. no cpu usage. Is there a way around this, or am I asking for an enhancement, e.g. a “writeListFile toFile” option. I am on RHEL 6.2/6.3 if that matter. Noticed that there is a local correlation option in the pipeline for Chimera 1.8, looking forward to try it out. Kind Regards, Ingvar Typical python script: from chimera import runCommand as rc rc("open emd_5500.map") rc("volume #0 level 1.34 transparency 0.5") rc("open emd_2017.map") rc("volume #1 level 0.58 transparency 0.5") rc("fitmap #1 #0 search 200 metric overlap envelope true inside 0.2 listFits false") #rc("fitmap #1 #0 search 200 metric overlap envelope true inside 0.2") rc("save session5500_2017") rc("stop now") Setup to run with Xvfb (in bash): i=0 while [ -f /tmp/.X$i-lock ] do i=$[$i+1] done Xvfb :$i & pid=$! export DISPLAY=:$i chimera fit5500_2017.py kill $pid

2 2

Pipes and planks representation
by Ernesto Arias Palomo 21 May '13

21 May '13

Dear chimera developers, I have noticed that in the last Chimera releases you have significantly improved the 'pipes and planks' representation. I find it particularly useful because, in my opinion, it makes the structures look simpler and clearer, especially when you are displaying large PDBs. I would like to suggest an improvement. It would be really useful if you could make that representation an attribute of the structure (like the 'ribbons' one) instead of a separate object. That way one could easily show, hide or color different regions of the structure without having to modify the actual PDB file. Thank you very much for your help, and congratulations for your great work. Cheers, Ernesto.

2 1

Re: [Chimera-users] RMSD i.b and RMSD u.b
by Elaine Meng 20 May '13

20 May '13

Please send Chimera questions to chimera-users(a)cgl.ucsf.edu instead of just my address. The RMSD lb (lower bound) and RMSD ub (upper bound) of a docked ligand are calculated by AutoDock Vina and are described in the manual for that program: <http://vina.scripps.edu/manual.html> (go to that page and search for "RMSD") I believe they just indicate how similar each position is to the best-scoring position found. I hope this helps, Elaine ---------- Elaine C. Meng, Ph.D. UCSF Computer Graphics Lab (Chimera team) and Babbitt Lab Department of Pharmaceutical Chemistry University of California, San Francisco On May 20, 2013, at 6:10 AM, Muhammad ALi <muh.ali741(a)yahoo.com> wrote: > AOA > > sir > when I check my result with ViewDock their are two RMSD values i.e RMSD i.b and RMSD u.b plz tell me about them.....

1 0

3D Printing our protein
by Bart Janssen 20 May '13

20 May '13

Hi My apologies for what are likely to be stupid questions. We are plant developmental biologists and we just obtained our first X-ray crystal structure for a new hormone receptor in plants. We'd really like to make a 3D model of the protein using one of the 3D printers at our institute. I've been trying to use and understand Chimera for the last week and I'm a bit out of my depth :). I can import our protein 4DNP.pdb and display it using solid surface or using the sphere representations of atoms. I can also display it with the "lid" domain hidden and vice versa. My problem has been getting files that our 3D printer can read. The solid surface model of the whole protein can be exported as a scene in .stl format and that can be printed just fine. Although I can't figure out how to adjust resolution or the smoothness of the exported file. But I'd really like to print the two domains of the protein separately. If I use Solid surface model with the lid hidden what I get is an eggshell like model which has no internal detail. If I use atoms shown as spheres instead I get the internal detail I want BUT when I try and export the scene I get a 2 Gb+ file that is unreadable by the printer. If I export .obj format instead the colleague who runs the printer says the file won't load and appears to have errors. I have a horrible feeling I am going about this entirely the wrong way but I can't find a tutorial that describes what I'm trying to achieve. I do realise chimera is not really written for the somewhat trivial use that I am attempting but it really would be helpful if we had a physical model of the protein to look at. Is there a tutorial I should follow or another part of Chimera (or another program) that I should be using instead? Any advice would be much appreciated. Cheers Bart Bart Janssen Scientist Plant Development team Breeding and genomics Plant and Food Research Private Bag 92169 Auckland 1142 New Zealand Ph +64 09 9257179 bart.janssen(a)plantandfood.co.nz The contents of this e-mail are confidential and may be subject to legal privilege. If you are not the intended recipient you must not use, disseminate, distribute or reproduce all or any part of this e-mail or attachments. If you have received this e-mail in error, please notify the sender and delete all material pertaining to this e-mail. Any opinion or views expressed in this e-mail are those of the individual sender and may not represent those of The New Zealand Institute for Plant and Food Research Limited.

5 7

showing salt bridges in pdb trajectory movie
by Albert Lyons 18 May '13

18 May '13

Hi all, I am having a problem trying to show salt bridges between selected atoms in a movie of my pdb format trajectory. I am using the following commands in the Per-Frame script, to display H-bonds between the residues but with a less strong distance criterion (to bring the distance of interaction to 6 A) and angle criterion (to make the angle 0deg and so make all angles available): select :1,6,7,10,16,18 findhbond selRestrict :1@HZ1,HZ2,HZ3 :6@2HH2,2HH1,HE,2HH1,1HH1 :7@O3P,O2P,O1P :10@HZ1,HZ2,HZ3 :16@HZ1,HZ2,HZ3 :18@2HH2,2HH1,HE,2HH1,1HH1 linewidth 2 reveal true relax true distSlop 2.97 angleSlop 145 When I apply this script I get a Chimera error saying: Error while sourcing select :1,6,7,10,16,18 findhbond selRestrict :1@HZ1,HZ2,HZ3 :6@2HH2,2HH1,HE,2HH1,1HH1 :7@O3P,O2P,O1P :10@HZ1,HZ2,HZ3 :16@HZ1,HZ2,HZ3 :18@2HH2,2HH1,HE,2HH1,1HH1 linewidth 2 reveal true relax true distSlop 2.97 angleSlop 145 ,line 1: "findhbond selRestrict :1@HZ1,HZ2,HZ3 :6@2HH2,2HH1,HE,2HH1,1HH1 :7@O3P,O2P,O1P :10@HZ1,HZ2,HZ3 :16@HZ1,HZ2,HZ3 :18@2HH2,2HH1,HE,2HH1,1HH1 linewidth 2 reveal true relax true distSlop 2.97 angleSlop 145" Non-alphanumeric character in keyword ':6@2HH2,2HH1,HE,2HH1,1HH1' I have tried tweaking it different ways but just get other error messages, I can't seem to write it in the correct format. Can anyone see a problem with this script and how I might fix it? Cheers Albie

2 1

showing salt bridges in pdb trajectory movie
by Albert Lyons 18 May '13

18 May '13

Hi all, I am having a problem trying to show salt bridges between selected atoms in a movie of my pdb format trajectory. I am using the following commands in the Per-Frame script, to display H-bonds between the residues but with a less strong distance criterion (to bring the distance of interaction to 6 A) and angle criterion (to make the angle 0deg and so make all angles available): select :1,6,7,10,16,18 findhbond selRestrict :1@HZ1,HZ2,HZ3 :6@2HH2,2HH1,HE,2HH1,1HH1 :7@O3P,O2P,O1P :10@HZ1,HZ2,HZ3 :16@HZ1,HZ2,HZ3 :18@2HH2,2HH1,HE,2HH1,1HH1 linewidth 2 reveal true relax true distSlop 2.97 angleSlop 145 When I apply this script I get a Chimera error saying: Error while sourcing select :1,6,7,10,16,18 findhbond selRestrict :1@HZ1,HZ2,HZ3 :6@2HH2,2HH1,HE,2HH1,1HH1 :7@O3P,O2P,O1P :10@HZ1,HZ2,HZ3 :16@HZ1,HZ2,HZ3 :18@2HH2,2HH1,HE,2HH1,1HH1 linewidth 2 reveal true relax true distSlop 2.97 angleSlop 145 ,line 1: "findhbond selRestrict :1@HZ1,HZ2,HZ3 :6@2HH2,2HH1,HE,2HH1,1HH1 :7@O3P,O2P,O1P :10@HZ1,HZ2,HZ3 :16@HZ1,HZ2,HZ3 :18@2HH2,2HH1,HE,2HH1,1HH1 linewidth 2 reveal true relax true distSlop 2.97 angleSlop 145" Non-alphanumeric character in keyword ':6@2HH2,2HH1,HE,2HH1,1HH1' I have tried tweaking it different ways but just get other error messages, I can't seem to write it in the correct format. Can anyone see a problem with this script and how I might fix it? Cheers Albie

1 1

Re: [Chimera-users] Conserved amino acids
by Rebecca Swett 17 May '13

17 May '13

Use match-> align to get a structural superposition and sequence alignment, then set your alignment type to AL2CO. Go to render by attribute and then hit the select by attribute tab. Select the range from 0-whatever your maximum is. Should be around 4 if you've got any conserved tryptophans. That will select all invariant residues. Jean-Paul Boissel <boissel(a)uni-mainz.de> wrote: >Hi, > >By homology modelling, I got putative structures for five closely related >amino acid transporters. Is there any script that will allow the quick >selection and visualisation of the invariant residues on the structures? >Thanks, > >Jean-Paul > > > > > >_______________________________________________ >Chimera-users mailing list >Chimera-users(a)cgl.ucsf.edu >http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users

2 1

Conserved amino acids
by Jean-Paul Boissel 17 May '13

17 May '13

Hi, By homology modelling, I got putative structures for five closely related amino acid transporters. Is there any script that will allow the quick selection and visualisation of the invariant residues on the structures? Thanks, Jean-Paul

1 0

Carbon color of 1 chain
by Markus Heller 16 May '13

16 May '13

Hello, I don't seem to be able to figure out how to color a scene by element, but with different carbon colors for different chains. I can't figure out how to get an expression like "color green, chain C C" Or something like that ... How to go about this? Thanks and Cheers Markus -- Markus Heller, Ph.D. NMR Scientist CDRD - The Centre for Drug Research and Development 2405 Wesbrook Mall, Fourth Floor | Vancouver, BC V6T 1Z3 | Main: (604) 827-1147 Direct: (604) 827-1122 | F: (604) 827-1299 | E: mheller(a)cdrd.ca<mailto:kessabar@cdrd.ca> | www.cdrd.ca

4 7