January 2014 - Chimera-users - RBVI Mailing Lists

Re: [Chimera-users] rna and set shadows commands
by Elaine Meng 16 Jan '14

16 Jan '14

On Jan 9, 2014, at 12:07 PM, "pythonchemia Gazeta.pl" <pythonchemia(a)gazeta.pl> wrote: > > Hello Elaine C. Meng > > Thank you very much for your reply and sorry that I answer only now. > When it comes to command rna path and rna model your tips helped me. :) > Although I do not know how to define the path for the command operations rna duplex. : ( > > Unfortunately for me the system WinXP / Chimera 1.8 - 32bit shadows set command does not work. : ( > Chimera returns this message: > > - Unable to turn on shadows (shadows not supported) > > After opening the Chimera Reply log returns this message: > > Disabled programs because the GPU and graphics driver bug > you encountered while compiling a vertex shader. > -------------------------------------------------- ----------- > > Maybe my graphics card is inappropriate? > > The User's Guide about the Build Structure writes that the helical DNA / RNA (not yet Implemented ) > and this option is disabled. > > Sincerely, Janusz Grabowski Hi Janusz, The message suggests your computer can't do the shadows. It might help to update your graphics driver, but I don't know for certain. However, other "set" options should work, for example: set bgColor white As mentioned in the "rna duplex" description, the path file is created using the command "rna path": <http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/rna.html#duplex> <http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/rna.html#path> Also, as mentioned in an earlier message you have to get a daily build (version 1.9), also available from our download page, to use the "helical DNA/RNA" option in Build Structure. <http://plato.cgl.ucsf.edu/pipermail/chimera-users/2013-December/009497.html> It is recommended to send chimera questions to chimera-users(a)cgl.ucsf.edu instead of to me directly. Thanks, Elaine

2 1

Re: [Chimera-users] batch writing of attributes
by Eric Pettersen 14 Jan '14

14 Jan '14

Hi Jay, We thought as much really. Well, you are going to have to use a Python script. The good news is that it's a really simple Python script. Put the following in a file whose name ends in ".py": from chimera import openModels, Molecule for m in openModels.list(modelTypes=[Molecule]): from OpenSave import osOpen f = osOpen("~/surfData/%s.data" % m.name, "w") for a in m.atoms: print>>f, a, getattr(a, "areaSES", 0.0) f.close() The indentation matters, so make sure to preserve it. Then once you open your model and surface it, you can run the above script simply by opening it with the "open" command. It will create a file with a list of atoms and corresponding surface areas. The file will be named <model name>.data and it will be in the "surfData" folder in your home directory, so make sure that folder exists beforehand or modify the script to have it go to another location. --Eric Eric Pettersen UCSF Computer Graphics Lab http://www.cgl.ucsf.edu On Jan 14, 2014, at 5:56 AM, "Eifler, Jay Q. (UMKC-Student)" wrote: > Actually I'm doing ab-initio quantum mechanical calculations and need > to analyze on a atom and residue basis. I've already been able to > save the atom and residue results but can't figure out how to automate > it since you have got to click buttons in the render/select attribute. The > Chimera notes seem to imply it should be possible to save attributes to > file. However, areaSAS is supposed to be an attribute but doesn't seem > to belong to the model attributes. > From: Eric Pettersen > Sent: Monday, January 13, 2014 4:05 PM > To: Eifler, Jay Q. (UMKC-Student) > Cc: <Chimera-users(a)cgl.ucsf.edu> List > Subject: Re: [Chimera-users] batch writing of attributes > > On Jan 11, 2014, at 1:58 PM, Jay Eifler wrote: > >> >> I'm trying to automate saving the areaSAS for a bunch >> of models. I've got no problem loading the models, etc. >> The problem is simply how do I use command line >> arguments to save the areaSAS. I've tried writesel. This >> works all fine in Render/Select by attribute simply >> because it has a feature to save the areaSAS to a file. >> Thanks. > > Hi Jay, > Do you really need the values on a per-atom basis? If you just need the grand total, it's reported in the reply log so you could just open and surface your structures one by one via a script (one by one to avoid possibly running out of memory) and then afterward save the reply log contents to a file manually (there's a "Save" button on the reply log). > > --Eric > > Eric Pettersen > UCSF Computer Graphics Lab > http://www.cgl.ucsf.edu

1 0

Problem with HEADER keyword in PDBs
by Thomas Hrabe 14 Jan '14

14 Jan '14

Hi everyone, I believe I found a problem with the >HEADER< keyword in PDB files. Some programs may create pdbs with the header keyword specifying the origin or a name of a chain or other metadata with the coordinates below. PyMol understands this keyword and colors chains that have identical chain ID (see example attached) in multiple colors while chimera interprets all as one chain. Thanks, Thomas

3 2

MD trajectory play by commandline
by inhee park 14 Jan '14

14 Jan '14

Dear Chimera developers and users: I am using Chimera (production version 1.8 (build 38824) 2013-06-07 21:09:20 UTC). I wanted to process MD trajectories to extract their detailed H-bond information. To do so, I would like to process by script with Hbond.py file. $ chimera --script Hbond.py #---[Hbond.py] from chimera import runCommand # pr trajectory movie with open("metafile", 'w') as m: meta_input = """amber mypdb.prmtop pr_dry_netcdf.crd """ print >> m, meta_input # runCommand("open movie:metafile") runCommand("alias ^myhb hbond relax true savefile hbond.$1") runCommand("perframe myhb") runCommand("coordset #0 1,1000,10; wait 100") #----------------------- As I have many MD trajectories to be processed, I was looking for a commandline in lieu of actually clicking buttons in MD movie GUI to play. I thought that "coordset #0 1,10000,10; wait 100" may operate like automatically starting the MD movie to load frame 1 to 1000 at every 100 interval. But resulting saved H-bond files contain same information all about the first frame (because graphically displayed was the first frame without any update unless clicking the button in MD movie GUI). Is there a way to achieve an automatic play via commandline? Thank you, inhee park

4 4

RMSD per residue NMR structure
by cantini＠cerm.unifi.it 14 Jan '14

14 Jan '14

Dear users Does anyone knows how to calculate through CHIMERA the backbone rmsd per residue of a NMR ensembles of protein structures ? thank a lot Francesca Cantini ---------------------------------------------------------------- This message was sent using IMP, the Internet Messaging Program.

2 1

batch writing of attributes
by Jay Eifler 13 Jan '14

13 Jan '14

I'm trying to automate saving the areaSAS for a bunch of models. I've got no problem loading the models, etc. The problem is simply how do I use command line arguments to save the areaSAS. I've tried writesel. This works all fine in Render/Select by attribute simply because it has a feature to save the areaSAS to a file. Thanks.

2 1

Computing solvent accessible area without internal cavities
by Tom Goddard 13 Jan '14

13 Jan '14

Hi Numan, Chimera uses the MSMS program to compute solvent accessible surface area, and MSMS has an option to only compute for the outside surface component without internal cavities. So I think this is just what you want. Although you could script it in Chimera you might find it easier to write a script that uses the standalone program directly from Michel Sanner http://mgl.scripps.edu/people/sanner/html/msms_home.html The one trouble you may run into is the calculation failing on some molecules due to numerical instabilities. The older version of the MSMS library we use in Chimera often fails to compute the surfaces and areas for molecules larger than 10,000 atoms, and sometimes fails for much smaller molecules. The standalone program from Michel failed much less often in tests I did a few years ago, probably because it used some newer more reliable code. Tom On Jan 13, 2014, at 11:13 AM, "Oezguen, Numan" wrote: > Dear Tom, > > I am facing a problem that you have most probably already have solved and hoping you could help me out. I need to identify atoms that are solvent accessible. For this I used the program Getarea from Werner Braun’s group in Galveston Texas. It works fine, however it lists also internal cavity surface atoms as solvent exposed. Technically it is correct… but I am interested only in those surface areas that are not internal i.e. surface areas that are accessible to molecules approaching it from outside. > > Do you know of any commercial or free (preferably) program that can do this. Ideally it should be possible to call it from scripts (shell or Perl). > > Thank you very much for you time and help, > Numan > > > Numan Oezguen, Dr. rer. nat. > Instructor of Pathology & Immunology > Texas Children’s Microbiome Center > Feigin Center, Suite 830 > 1102 Bates Avenue > Houston, TX 77030

1 0

Registration ?
by Jean-Paul Boissel 13 Jan '14

13 Jan '14

Just downloaded the new dally build. At opening got the message: Registration file "/Users/JPB/.chimera/registration" has expired. Do we have to register again for 2014 Thanks > > Dr. Jean-Paul Boissel, Team Leader > Solute Transport Through Biological Membranes/ Ellen I. Closs' group > Department of Pharmacology > University Medical Center of the Johannes Gutenberg-University > Obere Zahlbacher Strasse 67, D-55101 Mainz > Tel.: +49 6131-17 9189 > FAX: +49 6131-17 9329 > e-mail: boissel(a)uni-mainz.de > http://www.unimedizin-mainz.de/pharmakologie/research/solute-transport-thro… > -biological-membranes.html?L=1 > This e-mail contains confidential information and may also be legally > privileged. If you are not the intended recipient or received this mail by > mistake, please inform the sender immediately and delete the mail. Any > un-authorized use, copying, disclosure or forwarding of this mail and of > information contained therein is prohibited.

2 1

displaying electrostatic field
by Bala subramanian 10 Jan '14

10 Jan '14

Friends, Kindly see the following link (Fig 22, left panel). http://www.bioscience.org/2009/v14/af/3398/fulltext.php?bframe=figures.htm Has anyone tried making the display of electric field in Chimera. Is there any way i can make the same in chimera. Thanks, Bala -- C. Balasubramanian

2 1

running wild with Chimera :-)
by Eric Pettersen 09 Jan '14

09 Jan '14

New Year Cards in 3D

1 0

Chimera-users January 2014