- Chimera-users - RBVI Mailing Lists

vrml lighting
by Dougherty, Matthew T 12 Feb '14

12 Feb '14

I have a multiple vrml models that vary by scale. 1, 10 & 100 I notice the lighting significantly drops off as the scale increases. When I have another model, say PDB, that is in approximate relation (ie, viewable in regards to the volume size similar to the vrml files volumes); the PDB file looks good regardless of size (eg, tRNA or ribosome). I can increase the brightness to compensate (10-20) but it washes out the PDB models when brightness goes beyond 2. is the lighting handled differently between these different types of models? Matthew Dougherty National Center for Macromolecular Imaging Baylor College of Medicine

2 1

superimposing two structures based on a specific selection
by David Martin 10 Feb '14

10 Feb '14

Chimera 1.8 on Windows I have been trying to align (superimpose) models based on a specific conserved region. Unfortunately Chimera doesn't seem to do what the documentation claims it does. I run (example): match #0:5-15@CA,20-30@CA #1: 5-15@CA,20-30@CA and get an RMSD calculation (more on this below), but no movement of the model occurs. I specify 'move true show true' as options. Still no movement but I have a transformation matrix. Unfortunately this is not in a form that can be understood by Transform Coordinates. How do I persuade my structures to superimpose on the specific atoms? The RMSD calculation appears to be broken. If I take the command rmsd #0:5-15@CA,20-30@CA #1: 5-15@CA,20-30@CA then I will get a number If I then run rmsd #0:5-15@CA,20-30@CA #1: 5-15@CA,21-31@CA I get exactly the same number. In fact, if I specify atom sets of different sizes (as long as it is not the first block) then I continue to get exactly the same number rmsd #0:5-15@CA,20-30@CA #1: 5-15@CA,20-35@CA This should fail but doesn't Have I completely misunderstood what Chimera should be doing here? ..d The University of Dundee is a registered Scottish Charity, No: SC015096

2 3

Query on introducing an external protein
by Aditya Padhi 06 Feb '14

06 Feb '14

Dear All, I want to place a protein molecule (appx. 20 kDa in size) into a viral capsid loop region which is about 7 amino acid long. What I want to do is to remove the 5 amino acids out from the loop region and place the 20 kDa protein, so that the loop region can accommodate the 20 kDa protein. I am finding it difficult as the co-ordinates of both the molecules are different. Although, this is not related to MD simulation or associated stuff, it would be grateful if I can get some suggestions. Thanking you, Regards Aditya.

2 3

Query on introducing a protein
by Aditya Padhi 04 Feb '14

04 Feb '14

Dear All, I want to place a protein molecule (appx. 20 kDa in size) into a viral capsid loop region which is about 7 amino acid long. What I want to do is to remove the 5 amino acids out from the loop region and place the 20 kDa protein, so that the loop region can accommodate the 20 kDa protein. I am finding it difficult as the co-ordinates of both the molecules are different. Although, this is not related to MD simulation or associated stuff, it would be grateful if I can get some suggestions. Thanking you, Regards Aditya.

1 0

Question about "Color Zone" and "Render by Attribute"
by Helge Paternoga 04 Feb '14

04 Feb '14

Hello, I have a question regarding the colouring of surfaces and their caps. I am trying to recreate the colouring by hydrophobicity on capped surfaces as shown in (on the right side): http://www.cgl.ucsf.edu/chimera/docs/ContributedSoftware/2dlabels/2dlabels.… So first I've opened the "Render by Attribute" dialog and assigned "kdHydrophobicity". This only coloured the surface and not the surface cap though. I am assuming this is done by "Color Zone"? Indeed this colours the surface cap but I have one problem: The colouring of the surface is different between "Render by Attribute" and "Color Zone" and I don't know why. The differences are small but I would like to know which one is more accurate (intuitively I would say the one from "Render by Attribute"). I've played around with the radius of "Color Zone" but this doesn't achieve the same colouring as "Render by Attribute". So my question is: How can I achieve the same surface colouring by "Color Zone" as in "Render by Attribute"? Does it have something to to with "bondzone" or the hierachy in which the atoms apply their respective colour in "Color Zone"? Thanks for reading! Helge

4 6

Question about opening mrc and PDB files
by Min Woo Sung 03 Feb '14

03 Feb '14

Hi, This is Min Woo Sung from Texas A&M University. I was using Chimera to see EM structure(.mrc file) and wanted to put a PDB file to fit that in EM map. I'm wondering if I can use Chimera to analyze the volume, size or domains of these two files for the purpose of analysis. I saw one PDB file seemed much bigger that EM map when I opened it in Chimera. Do I have to normalize the size of any files? If so, how can I do that? Any suggests would greatly appreciated. Best regard, Min Woo Sung

3 2

Creating output with IDATM naming
by Nikolay Igorovich Rodionov 03 Feb '14

03 Feb '14

Hi, Does anyone know a way to rename all of the atom names in a PDB file to their IDATM types, so that a new PDB file can be generated with the IDATM naming scheme? Best, Nikolay Rodionov

2 1

Zoom problem?
by Damien Larivière 03 Feb '14

03 Feb '14

Dear all, I try to zoom in the scene, to get very close to my model shown here below. I cannot be closer than in the image below, the zooming function is at its maximum. When zooming, I get the error message : /refocus to continue scaling/ I do not know whether it is a bug or something I'm doing wrong. Thanks for your help, Damien

5 5

Need help in Coulombic surface coloring
by Haresh Tukaram More 30 Jan '14

30 Jan '14

Hi, I tried to get the Coulombic surface coloring for one of the protein that I modified by removing certain residue at the C-termini and then mutated some residues to the non-natural amino acid 5,5,5-trifluoroleucine using Swisssidechin option. I am getting this error which I am not able to understand. Disabled GPU programs because a graphics driver bug was encountered while compiling a vertex shader. #0, chain A: COMP #0, chain B: COMP #0, chain C: COMP #0, chain D: COMP #0, chain E: COMP C:\Program Files\Chimera 1.8.1\bin\mscalc.exe 1.400000 2.000000 1 MSMSLIB 1.3 started on Local PC Copyright M.F. Sanner (March 2000) Compilation flags Surface CC TFLE.pdb, category main, probe radius 1.4, vertex density 2 22 connected surface components Total solvent excluded surface area = 10393.6 component areas = 8599.46, 624.842, 225.9, 216.153, 81.7059, 54.8894, 54.839, 48.0624, 45.9525, 43.9588, 41.5667, 39.1204, 38.7549, 37.0513, 35.5418, 33.6712, 32.9167, 28.9293, 27.8262, 27.7523, 27.7319, 26.9429 Total solvent accessible surface area = 10372 component areas = 9971.67, 202.012, 92.2225, 70.7414, 13.8165, 4.12717, 3.62753, 3.24856, 2.597, 1.1866, 1.72654, 1.33137, 1.08737, 0.962375, 0.64812, 0.231422, 0.40193, 0.120068, 0.0570046, 0.0515394, 0.0674059, 0.0311705 No SEQRES records for combination (#-2147483648.-2147483648) chain A; guessing terminii instead No SEQRES records for combination (#-2147483648.-2147483648) chain B; guessing terminii instead No SEQRES records for combination (#-2147483648.-2147483648) chain C; guessing terminii instead No SEQRES records for combination (#-2147483648.-2147483648) chain D; guessing terminii instead No SEQRES records for combination (#-2147483648.-2147483648) chain E; guessing terminii instead Chain-initial residues that are actual N terminii: #-2147483648.-2147483648 MET 27.A, #-2147483648.-2147483648 MET 27.B, #-2147483648.-2147483648 MET 27.C, #-2147483648.-2147483648 MET 27.D, #-2147483648.-2147483648 MET 27.E Chain-initial residues that are not actual N terminii: Chain-final residues that are actual C terminii: #-2147483648.-2147483648 TFLE 67.A, #-2147483648.-2147483648 TFLE 67.B, #-2147483648.-2147483648 TFLE 67.C, #-2147483648.-2147483648 TFLE 67.D, #-2147483648.-2147483648 TFLE 67.E Chain-final residues that are not actual C terminii: 424 hydrogen bonds Hydrogens added Charge model: AMBER ff12SB Assigning partial charges to residue TFLE+TFLE (net charge +0) with am1-bcc method Running ANTECHAMBER command: C:/Program Files/Chimera 1.8.1/bin/amber12/bin\antechamber -ek qm_theory='AM1', -i c:\users\haresh~1.pol\appdata\local\temp\tmpipssqi\ante.in.mol2 -fi mol2 -o c:\users\haresh~1.pol\appdata\local\temp\tmpipssqi\ante.out.mol2 -fo mol2 -c bcc -nc 0 -j 5 -s 2 (TFLE+TFLE) Cannot automatically determine charges for residue TFLE+TFLE; Run AddCharge tool manually to add charges and then rerun ESP Can someone please tell me a step by step procedure to make this work. Thanks & Regards, Haresh

4 13

Re: [Chimera-users] Coloring map by curvature
by Tom Goddard 29 Jan '14

29 Jan '14

Hi Joel, I don't think surface curvature will be helpful understand the differences between two EM maps. But I was curious what it would look like. So here's a Python script and an image on a simulated map at 15 Angstroms. You use the script by opening your map, selecting the surface (ctrl-click on it), then open the curvature.py script (menu File / Open…). Here's a description of the method it uses from the comments at the top of the curvature.py file. # Color selected surface pieces by mean curvature. # # Gray at the average curvature value over the surface, and blue and red # at +/- 3 standard deviations of curvature values across the surface. # # The curvature is estimated from the vertices and normals of the triangulated # surface in simple way which will show artifacts from non-isotropic meshes. # For each triangle edge it computes the normal vector rotation from one vertex # to the other divided by the edge length. The vertex mean curvature is the # mean of the curvatures computed for each edge. # And the steps to make the example image open 1grl molmap #0 15 grid 2 model #1 select #1 open ~/Desktop/curvature.py Tom On Jan 29, 2014, at 9:40 AM, Joel Meyerson wrote: > Hi Darrell, > Thanks for the tip! I hadn't considered mesh lab but that link you sent looks promising. I'll also wait to see if Tom has any further suggestions. > Joel > > > On Wed, Jan 29, 2014 at 11:44 AM, Hurt, Darrell (NIH/NIAID) [E] wrote: > Hi Joel, > > I love Chimera and use it all the time. However, a quick-and-dirty solution to your problem might be to try Meshlab. It is kind of buggy software, but can be very useful. Here's a tutorial/blog entry that I found in a quick search. At the very least, it describes what I think Tom was communicating about enhancing surface shading: > http://meshlabstuff.blogspot.com/2010/03/mean-curvature-cavity-map-zbrush-a… > > FWIW, > Darrell > > -- > Darrell Hurt, Ph.D. > Section Head, Computational Biology > Bioinformatics and Computational Biosciences Branch (BCBB) > OCICB/OSMO/OD/NIAID/NIH > > 31 Center Drive, Room 3B62B, MSC 2135 > Bethesda, MD 20892-2135 > > > Disclaimer: The information in this e-mail and any of its attachments is confidential and may contain sensitive information. It should not be used by anyone who is not the original intended recipient. If you have received this e-mail in error please inform the sender and delete it from your mailbox or any other storage devices. National Institute of Allergy and Infectious Diseases shall not accept liability for any statements made that are sender's own and not expressly made on behalf of the NIAID by one of its representatives. > > From: <Meyerson>, "Joel [F] (NIH/NCI)" > Date: Tuesday, January 28, 2014 9:08 PM > To: Tom Goddard > Cc: "chimera-users(a)cgl.ucsf.edu List" <chimera-users(a)cgl.ucsf.edu> > > Subject: Re: [Chimera-users] Coloring map by curvature > > Hi Tom, > I have cryo-EM maps for two conformations of a protein, both in the 15 Angstrom resolution range. The conformations are visibly different, but I am also interested in seeing where they differ in terms of their surface curvature, as it could have bearing on how I interpret the conformations. If there's any other info I can provide just let me know. > Thanks! > Joel > > > On Tue, Jan 28, 2014 at 7:36 PM, Tom Goddard wrote: > Hi Joel, > > No, Chimera does not compute surface curvature. It would not be too hard to make a Python script that computed it and used it to color a surface. The main trouble is defining numerically the curvature for a triangulated surface at each vertex. Why are you interested in this? Is the idea to simulate ambient occlusion lighting where surface cavities are dark and projections are brighter? Is the idea to try it on EM maps or molecular surfaces? > > Tom > > > On Jan 28, 2014, at 2:29 PM, Joel Meyerson wrote: > > > Hi, > > Is it possible to color a map based on surface curvature? > > Thanks, > > Joel > > _______________________________________________ > > Chimera-users mailing list > > Chimera-users(a)cgl.ucsf.edu > > http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users > > >

2 1

Coloring map by curvature
by Joel Meyerson 29 Jan '14

29 Jan '14

Hi, Is it possible to color a map based on surface curvature? Thanks, Joel

4 4

adding GalNac moiety on Threonine residue
by Elsa Marquez 28 Jan '14

28 Jan '14

Hi there! Do you guys know if I can add or somehow draw a carbohydrate moiety such as GalNac to a single Threonine in a protein complex? Thanks! Elsa Marquez University of Alberta

2 1

Surface Transparency Command Issues
by Joshua Carter 28 Jan '14

28 Jan '14

Hello, I'm having trouble getting the surface transparency command to work correctly. I'm currently making a movie about the formation of a protein complex and I want the subunits to be shown as ribbons while they dock and then after they dock the electrostatic surface to appear. I have the script written and the surface generated but when I use either the surftransp or transparency command nothing happens. The command is correct, if I use it to change the transparency of a ribbon or atom it does so, but if I try to make the surface appear after docking nothing changes. In addition, if I go to the Actions menu and change the surface transparency from there it also works. If you could help me figure out why the command won't work that would be extremely helpful. Thanks,Josh

2 1

Re: [Chimera-users] Overlay structures in Chimera (UNCLASSIFIED)
by Elaine Meng 28 Jan '14

28 Jan '14

On Jan 28, 2014, at 9:45 AM, "Keasey, Sarah L CTR USARMY MEDCOM USAMRIID (US)" <sarah.l.keasey.ctr(a)mail.mil> wrote: > Elaine, > > I have a new question regarding overlaying structures using Chimera. > > I have opened two structures in Chimera. I am trying to move them so that they are overlaid. It appears as though I should be able to inactivate one of the structures in order to freely move the other, but this doesn't work for me and I'm not sure where I am going wrong. When I inactivate one structure, I can freely rotate the other, but not move it from one place to another. I have watched the youtube video here: http://www.youtube.com/watch?v=6VGYo1pRRZ8 which deals with merging two structures, and also watched the "Sequential Fitting" video on the "Videos" page of the UCSF Chimera site. > > Any suggestions for what I am doing wrong? I apologize for such elementary questions, I am a new Chimera user. > > Thanks in advance, > Sarah Keasey Hi Sarah, (we generally ask everyone to send questions to chimera-users(a)cgl.ucsf.edu, CC'd here, so that everybody can benefit, unless the data are private) Those videos aren't that relevant to your situation. More relevant is the manual page discussing the various ways to superimpose structures. The Chimera User's Guide can be searched or browsed from the Help menu, or you can see the same thing on our website, e.g.: <http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/superposition.html> If using the mouse to move one structure, what you did is exactly correct: inactivate the other structure. The mouse should continue to work exactly the same for the movable structure as when everything was movable (both rotation and translation). The default control for XY-translation is middle mouse button, if you have one. For touchpad or mouse with fewer buttons, see here: <http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/mouse.html> However, more often one doesn't use the tedious approach of trying to do it by hand, and instead uses one of the other superposition methods. The easiest is usually MatchMaker (in the menu under Tools… Structure Comparison), so I would try that first. <http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/matchmaker/matchm…> To help you become more familiar with Chimera, there are several tutorials in the User's Guide (again also available from Help menu), ranging from "getting started" to more advanced. <http://www.rbvi.ucsf.edu/chimera/current/docs/UsersGuide/frametut.html> 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

1 0

Multialign Viewer Header Definition File (UNCLASSIFIED)
by Keasey, Sarah L CTR USARMY MEDCOM USAMRIID (US) 24 Jan '14

24 Jan '14

Classification: UNCLASSIFIED Caveats: NONE Hello, I would like to generate a Header Definition File to load in to Multialign Viewer. I am unclear on the format of the file though, I keep getting errors when I try to load the file in and I think it has something to do with the way I am formatting the file. Would you be able to provide an example Header Definition file? Thank you in advance, Sarah Keasey Graduate Student University of Maryland, Baltimore County ORISE participant, USAMRIID, Fort Detrick MD Classification: UNCLASSIFIED Caveats: NONE

2 3

Segger : command line
by agnel-praveen.joseph＠stfc.ac.uk 24 Jan '14

24 Jan '14

Hi, I would like to automate watershed segmentation (Segger) on a set of volumes. Given that I have the details on the volume threshold, number of smoothing steps and minimum voxels in each region, is there a way to automate segmentation with different parameter inputs? Many Thanks. -- Scanned by iCritical.

2 2

Chimera Issues
by Roy Assis 21 Jan '14

21 Jan '14

Hello ! For some time i've been using chimera, but for some reason a week ago, when I try to lunch chimera, I get an error message (attached). Please help Roy Assis

2 1

movie recode
by Wenda Wang 21 Jan '14

21 Jan '14

Dear List, I am trying to use UCSF chimera to make a movie after TEM reconstruction using IMOD. But an error message always pops out, even I use the script in the examples. Belows shows the error messages in the Reply Log. Can anyone take a look and give me some clue where the problem is? Thank you. ----------------------------- C:\Program Files\Chimera 1.8.1\bin\ffmpeg.exe -r 25 -i c:\users\337\appdata\local\temp\chimovie_zUrR-%05d.ppm -vf crop=(in_w/2)*2:(in_h/2)*2:0:0 -y -vcodec libx264 -f mp4 -vb 200.0k C:\Users\Public\Desktop\rock.mp4 An error occurred during encoding. See Reply Log for details. Error during MPEG encoding: ----------------------------- Exit value: 1 Error message: ----------------------------- C:\Program Files\Chimera 1.8.1\bin\ffmpeg.exe -r 25 -i c:\users\337\appdata\local\temp\chimovie_zUrR-%05d.ppm -y -vcodec libvpx -f webm -vb 200.0k C:\Users\Public\Desktop\rock.webm An error occurred during encoding. See Reply Log for details. Error during MPEG encoding: ----------------------------- Exit value: 1 Error message: ----------------------------- C:\Program Files\Chimera 1.8.1\bin\ffmpeg.exe -r 25 -i c:\users\337\appdata\local\temp\chimovie_zUrR-%05d.ppm -y -vcodec libtheora -f ogg -vb 200.0k C:\Users\Public\Desktop\rock.ogv An error occurred during encoding. See Reply Log for details. Error during MPEG encoding: ----------------------------- Exit value: 1 Error message: Thank you very much, Wenda Wenda Wang, PhD candidate Department of Materials Sciences and Engineering, Drexel University, Philadelphia, PA 215-6883348

2 1

Morphing
by James Naismith 20 Jan '14

20 Jan '14

I am experimenting with the morphing feature. I know I am being very stupid but I cannot get the analysis to work. I am using a MAC OSX latest version. I have created my morph which is running nicely. What I would like to do is analyse the morph for VDW clashes. (That is in going from structure 1 to structure 2, what residues would have to move). First problem is when I run cluster (under analysis in MD movie) then generate residue interaction for cluster, up pop a new dialogue with warnings but seems to ask about calculations When I click OK or apply I get a new box residue interaction contact parameters This always generates an error “No atoms designated for clash detection’ I seem only to have four options themselves all others atoms (default) other atoms in same model second set of designated atoms Any help gratefully received (Ps I have installed the latest cytoscope but it does not have the structure.viz plugin; I keep getting errors trying to install older version) best Jim Jim Naismith The University of St Andrews is a charity registered in Scotland : No SC013532

2 2

Re: [Chimera-users] [Chimera] #12778: failure assigning charges to FAD (was: Chimera bug report submission)
by Eric Pettersen 17 Jan '14

17 Jan '14

Hi Will, The underlying program that Chimera uses to compute charges, antechamber/sqm from AmberTools, frequently fails for molecules with regions of concentrated negative charge, such as the central phosphate connector of FAD. It not really something I can do much about. Nonetheless, the failure not only depends on the charge density but also on the specific atomic geometry. So for instance, I can compute the charges for the FAD in 4kpu, though it does take a long time (about 10 minutes). This means that there is kind of a work around. Compute the charges in a structure that works and then transfer them to your structure. You do this by bringing up the Render by Attribute tool after successfully computing charges and using it's File->Save Attributes menu entry to save the charges to a file. You should select the FAD and then restrict the save to selected atoms. Edit the resulting file and change the residue name and chain ID (in this case 401.A) to FAD so that it can be applied to all FADs in your structure. Then in the Chimera with your structure open, read in the charge-attribute file with the Define Attributes tool. To save you a lot of the effort here, I've attached the edited charge-attribute file I saved from 4kpu. --Eric Eric Pettersen UCSF Computer Graphics Lab http://www.cgl.ucsf.edu

1 0

how to display coordinate axes
by Michael H Peters 17 Jan '14

17 Jan '14

Hi: I have a simple question - How do I display the lab or fixed frame coordinate axes? I am using the windows 32 bit UCSF Chimera on my desktop. In rasmol the command is "axes on". I cannot find anything similar in Chimera. Thanks for your help. An axes display is critical if you need to know the spatial orientation wrt the lab frame. Best wishes, Michael Dr. Michael H. Peters mpeters(a)vcu.edu 804-828-7790 www.engineering.vcu.edu/proteinengineering

3 2

pyroglutamate distorted in 3PFJ
by Craig Dana 17 Jan '14

17 Jan '14

I'm looking at the N-terminal pyroglutamate residue in PDB: 3PFJ. It is clearly distorted, but when I open the same PDB file in VMD, the pyroglutamate is not distorted.

3 2

Failed to start chimera
by jupiter 17 Jan '14

17 Jan '14

Hi, I installed Chimera 1.7 on CentOS 6.4. I can run any opengl applications, but failed to run chimera , it got following error message. Appreciate any tips to fix this issue. X Error of failed request: BadMatch (invalid parameter attributes) Major opcode of failed request: 72 (X_PutImage) Serial number of failed request: 34 Current serial number in output stream: 37

1 1

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