- Chimera-users - RBVI Mailing Lists

Find pseudo 2-fold axes relating heavy and light chain
by Jian Guan 30 Dec '13

30 Dec '13

Dear all, I am trying to find pseudo 2-fold axes relating heavy and light chains in the constant domain and variable domain respectively in Fab. I tried Define Axes, it failed. And I made artificial double heavy chain, it failed again. I can only get axis which is horizontal passing through both light chain and heavy chain. It is not what I want. Can anyone help to give some clues? Thanks Jian

3 4

Failure of energy miminization
by Ross KK Leung 30 Dec '13

30 Dec '13

In August there was a discussion entitled "Minimization and mystery nonstandard HIS atom name" and that's the exact problem supposed to have been solved last year. What I did is similar to Jason (http://www.cgl.ucsf.edu/pipermail/chimera-users/2013-August/009069.html) and unfortunately I'm unable to resolve the problem, giving that I'm using chimera 1.8.1 on a Win 7 32G RAM computer. I attach two files, one succeeded till completion while the other one failed early. What I did is simply "Structure editing", next, ok, next, .. Charges of 0.0 were assigned to the unknown atoms 1 model(s) had non-integral total charge Details in reply log No MMTK name for atom "HC1" in standard residue "ALA" Any suggestions to fix the problem is highly appreciated because Chimera is so impressive and I've been using it for many good analyses already. A pity to give up and switch to other programs at this very last stage.

2 1

Reset surface programatically from python
by Thomas Hrabe 30 Dec '13

30 Dec '13

Hi everyone, I have this line of code where I want to load a density and export a negative surface of it: import chimera as ch mod = ch.openModels.open(‘./density.em')[0] mod.surface_levels = [-1,0] mod.display = True ch.exports.doExportCommand('WebGL', 'tmp.html’) However, the surface of the density is rendered for positive values. How can I refresh the display after I specify the surface_levels? Thanks, Thomas

2 6

Re: [Chimera-users] [chimera-dev] writing in b-factor field
by Elaine Meng 30 Dec '13

30 Dec '13

On Dec 29, 2013, at 8:43 AM, MARIA VITTORIA CUBELLIS <cubellis(a)unina.it> wrote: > Dear sir/madam, > I have a file that associates to each amino acid of a protein A a given value v. > I would like to transfer this value to the b-factor field of the protein A the value v. I have a single value for each aminoacid N1,N2... I would like to transfer the value to each atom > > let say > > N1 C v1 > N1 CO v1 > N1 N v1 > .... > .... > N2 C v2 > N2 CO v2 > N2 N v2 > > and so on > > is it possible to this withe chimera? > thank you > > mvcubellis Dear M V Cubellis, Yes, you can do it, but another question is why do you want to do it? I'm guessing you want to color by those values. In Chimera you can assign the values to the atoms or residues directly as "attributes" and then do coloring to show those values. The B-factor from an input file is one atom attribute named bfactor, but you can create additional attributes with whatever names you wish, and color by any of them. The main work for you is to create an "attribute assignment file" which is mostly tab-separated columns with the atom (or residue) specifiers and values. The necessary format is described here, and there are also some example files: <http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/defineattrib/defi…> Probably the trickiest part is how to specify the residues in the file, e.g. could be something like :12.a for residue 12 in chain A. See instructions on atom/residue specification: <http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/atom_spec.html#basic> Once you have created that text file, you can read it in with Define Attribute (in menu under Tools… Structure Analysis) or command defattr. <http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/defineattrib/defi…> <http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/defattr.html> Then Render by Attribute (which may automatically appear when you create the attribute, but also is in menu under Tools… Structure Analysis) can be used to show the values with colors and/or radii: <http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/render/render.htm…> If you really wanted to write a PDB file with values in the B-factor column, you would still need to create an attribute assignment file as described above, but more specifically to assign an atom attribute named bfactor. Then when you use the File menu or write command to output a PDB file, your values will be in the B-factor column. See also the Attributes tutorial for examples of coloring by attribute value. <http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/tutorials/attributes.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

keyboard shortcuts
by Dougherty, Matthew T 30 Dec '13

30 Dec '13

I am writing a bunch of USB-HID keyboard shortcuts for some extensions I am writing. Usually one shortcut per def. Since my decoding of HID is integrally tied to the extension, I would like to simplify all my shortcuts coming from one device to one def. As of now if I have 100 button commands I have to 1) register 100 shortcuts, 2) create 100 defs, 3) do corresponding stuff in my extension. I would prefer to call one def that passes the actual keyboard string to my extension, thereby simplifying my user accelerators to just registration. How do I get the string within my keyboard shortcuts python file? thanks, Matthew Dougherty National Center for Macromolecular Imaging Baylor College of Medicine

2 1

display of improper peptide bond
by Tatiana Domitrovic 30 Dec '13

30 Dec '13

Hi, I am trying to get rid of an improper bond that chimera is displaying between two nearby segments that are actually cleaved. The bond is distorted, creating long CACB bonds if B-spline mode is active. I tried to turn off the auto-chaining, but didn’t work. Because PYmol can properly display these segments as independent polypetides, I was wondering if there is other adjustment that I can try before manually deleting the bond. The PDB is 1ohf and the bond is between residue 570 and 571. Cheers,

5 7

Failure of energy miminization
by Ross KK Leung 29 Dec '13

29 Dec '13

In August there was a discussion entitled "Minimization and mystery nonstandard HIS atom name" and that's the exact problem supposed to have been solved last year. What I did is similar to Jason (http://www.cgl.ucsf.edu/pipermail/chimera-users/2013-August/009069.html) and unfortunately I'm unable to resolve the problem, giving that I'm using chimera 1.8.1 on a Win 7 32G RAM computer. I attach two files, one succeeded till completion while the other one failed early. What I did is simply "Structure editing", next, ok, next, .. Charges of 0.0 were assigned to the unknown atoms 1 model(s) had non-integral total charge Details in reply log No MMTK name for atom "HC1" in standard residue "ALA" Any suggestions to fix the problem is highly appreciated because Chimera is so impressive and I've been using it for many good analyses already. A pity to give up and switch to other programs at this very last stage.

1 0

Sir, i donot want to recieve antomately letter from chimera, can you remove my email adress from your email list..
by 一波 28 Dec '13

28 Dec '13

2 1

How to get central cross section which pass through two asigned cross axis
by Jian Guan 27 Dec '13

27 Dec '13

Dear all, I have a density map of virus-Fab complex. I want to make a central cross section which pass through the Fab pseudo-symmetry 2-fold axis. I have measure the angle of Fab relative to the capsid. Hope someone can help me. Thank you so much. Sincerely yours, Jian -----Original Message----- From: chimera-users-bounces(a)cgl.ucsf.edu [mailto:chimera-users-bounces@cgl.ucsf.edu] On Behalf Of chimera-users-request(a)cgl.ucsf.edu Sent: 2013年12月24日 15:00 To: chimera-users(a)cgl.ucsf.edu Subject: Chimera-users Digest, Vol 128, Issue 25 Send Chimera-users mailing list submissions to chimera-users(a)cgl.ucsf.edu To subscribe or unsubscribe via the World Wide Web, visit http://plato.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: saving coords from script (Eric Pettersen) ---------------------------------------------------------------------- Message: 1 Date: Mon, 23 Dec 2013 15:59:16 -0800 From: Eric Pettersen <pett(a)cgl.ucsf.edu> To: Tetsuro Fujisawa <fujisawa(a)gifu-u.ac.jp> Cc: "chimera-users(a)cgl.ucsf.edu Mailing List" <chimera-users(a)cgl.ucsf.edu> Subject: Re: [Chimera-users] saving coords from script Message-ID: <04B15E0C-844B-4594-8AE6-760588744D85(a)cgl.ucsf.edu> Content-Type: text/plain; charset="windows-1252" On Dec 20, 2013, at 1:58 AM, Tetsuro Fujisawa wrote: > Hi, > > I found some strange behavior in saving pdb file by script file. > I tried to modify the orientation of pdb file ?tst.pdb? (#1) in reference to axis coordinate (#0), then save the atomic coordinates in pdb file. > When I typed in the commands through either command line window or IDLE, transformed coordinates were saved. It worked fine. However, the exact identical commands (inside the brancket) or lines (rc(?..?)) were run by script file shown in below, the untransformed coordinates were saved. This is also the case for ?Midas.write method?. > I wondered this might be a bug. > > Regards, > Tetsuro > > Sample.py: > ---------------------------------------------- > from chimera import runCommand as rc > > rc('open axis.bild') > rc('open tst.pdb') > rc('rock y 60 1 center #1 coord #0 models #1') # orientation of > tst.pdb was modified > rc('focus') > rc('write relative #0 #1 tstA.pdb') > # write current coordinates to tstA.pdb rc('open tstA.pdb') # opened > tstA has original untransformed orientation > ----------------------------------------------- Hi Tetsuro, The issue is that unlike using the command line or IDLE or running a chimera-command script, no frames get drawn during the execution of a Python script unless explicitly requested and therefore the "rock" command in your script doesn't actually do anything until after the entire script runs, and by that time your coordinate file has already been written. To get the rock command to complete before proceeding with the rest of the script you need to append "; wait" or "; wait 1" (since you are rocking for 1 frame) to the end of your rock command. --Eric Eric Pettersen UCSF Computer Graphics Lab http://www.cgl.ucsf.edu

2 1

saving coords from script
by Tetsuro Fujisawa 23 Dec '13

23 Dec '13

The following bug report has been submitted: Platform: windows Chimera Version: 1.8.1 Description Hi, I found some strange behavior in saving pdb file by script file. I tried to modify the orientation of pdb file ‘tst.pdb’ (#1) in reference to axis coordinate (#0), then save the atomic coordinates in pdb file. When I typed in the commands through either command line window or IDLE, transformed coordinates were saved. It worked fine. However, the exact identical commands (inside the brancket) or lines (rc(’..’)) were run by script file shown in below, the untransformed coordinates were saved. This is also the case for “Midas.write method”. I wondered this might be a bug. Regards, Tetsuro Sample.py: ---------------------------------------------- from chimera import runCommand as rc rc('open axis.bild') rc('open tst.pdb') rc('rock y 60 1 center #1 coord #0 models #1') # orientation of tst.pdb was modified rc('focus') rc('write relative #0 #1 tstA.pdb') # write current coordinates to tstA.pdb rc('open tstA.pdb') # opened tstA has original untransformed orientation -----------------------------------------------

2 1

Color by ramachandran probability? Or molprobity overlap?
by Oliver Clarke 19 Dec '13

19 Dec '13

Hi all, Quick question, just wondering if there is any way to color a structure by Ramachandran probability? Or by steric overlap scores for each atom as generated by mol probity? Either of these would be very handy as a way of getting a quick graphical overview of trouble spots in a structure. It seems like Chimera “knows about” Ramachandran plot data, (as one can generate such plots from the model panel) but I can’t find any convenient way to display this data on the structure. Many thanks, Oliver.

2 3

Running Segger from a script
by Tom Goddard 19 Dec '13

19 Dec '13

Hi Rakesh, There is not Chimera command to run Segger, only a dialog. So you would have to use Python to run it and write out map files for each masked region if you wanted to script it. Segger uses several Chimera features such as Gaussian filtering, and watershed calculations that are in written in C++ and compiled against the Chimera Python. While you can get those to run in another Python if it has an identical version (e.g. 2.7.3), setting up the environment variables to find all the needed libraries is a difficult process. So I suggest using the Chimera Python. You can run Chimera scripts without the graphical interface by starting it from a shell "chimera --nogui myscript.py". The Segger Python functions are not documented so you would need to study the (complex and messy) code found in chimera/share/Segger or on Mac in Chimera.app/Contents/Resources/share/Segger Tom On Dec 18, 2013, at 5:54 PM, rakesh ramachandram wrote: > Hi, > > I would like to know if I can run Segger in a script in an automated fashion and then obtain the individual segmented regions as separate density maps. Is there a way to run chimera scripts from command line directly without using the chimera's python interpreter but using the default python interpreter after importing its libraries. > > Regards > Rakesh >

1 0

question
by wanxiaobo＠nibs.ac.cn 17 Dec '13

17 Dec '13

Dear chimera team / users, I was looking for a way to save the content of Axis/plane/centroid information, not only the radii via command line, because I want to analyze the structure movement from MD trajectory. It is not convenient to save the information manually for each frame. Thanks Xiaobo

2 1

Command line export to webgl
by Thomas Hrabe 17 Dec '13

17 Dec '13

Dear Chimera Team, I was curious if there was a functionality to export chimera scenes to webgl through the command line, without having to open the chimera window. Is there a set of python script commands that would do it? Thank you for your support, Thomas

2 1

Attribute Calculator by command line?
by Hurt, Darrell (NIH/NIAID) [E] 16 Dec '13

16 Dec '13

Hi Chimera friends, I'm trying to calculate convexity as demonstrated in the Attributes tutorial: http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/tutorials/attributes.html I can do it with the graphical tools without a problem, but when I try to do it via the command line, I run into trouble. I cannot find any documentation on running the Attributes Calculator on the command line. I've tried using "setattr": setattr a convexity atom.areaSAS/atom.areaSES This command generates a new atom attribute called "convexity" for every atom, but the value of this attribute for every atom is literally "atom.areaSAS/atom.areaSES" instead of the calculated value. Is there a way to calculate a new attribute value on the command line within Chimera? Sorry to be so demanding… Thanks, 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 Office: 301-402-0095 Mobile: 301-758-3559 Web: BCBB Home Page<http://www.niaid.nih.gov/about/organization/odoffices/omo/ocicb/Pages/bcbb.…> Twitter: @niaidbioit<https://twitter.com/niaidbioit> 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.

3 4

Setting selection mode command line
by rswett＠chem.wayne.edu 16 Dec '13

16 Dec '13

Hi all, I'm interested in putting together a script that will open several large pdb files, select some very specific pieces of those files and save the output to another file. The problem is that they are large files and unwieldy in a graphics window. I have no problem generating the commands to select the regions I'm interested in, but I still need to be able to switch selection modes in between commands. Is there a way to change the selection mode at the command line without navigating the selection menu? Rebecca

3 5

Export options for x3d2stl
by Hurt, Darrell (NIH/NIAID) [E] 16 Dec '13

16 Dec '13

Hi Chimera friends, Here's the second question of the day. Thanks for your patience with me. When I export the scene as an STL (for 3D printing), I think what happens is that Chimera outputs first an X3D file and then runs the utility program "x3d2stl" to create the STL. The X3D file is about 5 MB, but the STL is about 100 MB. I understand that much of this size inflation is because the "primitives" in X3D have to be converted to mesh objects for STL. However, when I look at the STL file, the meshes derived from the primitives are finer than I really need. I find that I can "get away with" a coarser meshification by using x3d2stl on the command line: x3d2stl -c -r 2 -o sticks.stl < sticks.x3d If I want to export the scene via the command line within Blender, is there any way to issue options to x3d2stl? Or should I just run x3d2stl outside of Blender? Thanks, 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 Office: 301-402-0095 Mobile: 301-758-3559 Web: BCBB Home Page<http://www.niaid.nih.gov/about/organization/odoffices/omo/ocicb/Pages/bcbb.…> Twitter: @niaidbioit<https://twitter.com/niaidbioit> 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.

2 2

import binary format .raw file
by Carol J Hirschmugl 16 Dec '13

16 Dec '13

Hello, I am new to chimera, but have seen its capabilities for rendering 3D images of complex data sets, and think that this will be an appropriate program for excellent image analysis for my 3D tomography data sets. My data sets are in general binary format, (.raw file), 32 float, big endian, x fastest order for an x,y,z data set of 128 x 128 x 128 voxels of grey scale information. I am using AVIZO and FIJI at present. Avizo is expensive, and fiji is slow and reduces the data to 8 bit before rendering it, which I find unacceptable. TO move to Chimera - I found the web page in the Users Manual that suggested ADDING A FILE TYPE But I am lost by the description of how to proceed since I do not have prior knowledge of python and how to do the mapping to CHIMERA (x,y,z) that is suggested. I gather I need to write a python code - or numpy code maybe. Can you help or provide more context or a self-tutorial link so I can fill in the details? Thank you, Carol

2 1

Local execution of pdb2pqr and APBS
by Hurt, Darrell (NIH/NIAID) [E] 16 Dec '13

16 Dec '13

Hi Chimera friends, Sorry to bother you again with so many questions. I'm trying to use PDB2PQR and APBS in a script to calculate an EP map. I know how to specify a local executable location for these external applications if I use the graphics interface, but the documentation that I can find does not indicate how to do this via the command line. Is it possible? How do I do it? Thanks! 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 Office: 301-402-0095 Mobile: 301-758-3559 Web: BCBB Home Page<http://www.niaid.nih.gov/about/organization/odoffices/omo/ocicb/Pages/bcbb.…> Twitter: @niaidbioit<https://twitter.com/niaidbioit> 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.

2 2

changing the rgb for stereo images
by Bala subramanian 13 Dec '13

13 Dec '13

Friends, I have a glass with a red-cyan color. When i use the Chimera stereo mode (red-cyan) under view controls. I still see some ghost images. I would like to know if there is any possibility to add a features something like a box where i can give the rgb (the light frequencies) values of the images to match my glass. Thank you, Bala -- C. Balasubramanian

2 1

Dynamic scolor custom cmap?
by Hurt, Darrell (NIH/NIAID) [E] 13 Dec '13

13 Dec '13

Hi Chimera friends, I have a volumetric potential file loaded that I want to map onto a molecular surface. For my application, I want the color mapping to use a custom palette of at least three colors that spans the full range of surface values. I know that I can do this by simply specifying the custom colors in the Surface Color panel and then clicking "Set full range of surface values". That works for me. However, I want to do this via command line. The following line works for me (I've omitted "cmapRange full" because that is the default setting): scolor #0 volume #1 offset 1.4 cmap cyanmaroon My problem become if I want to set a custom palette: scolor #0 volume #1 offset 1.4 cmap -10,yellow:0,white:10,green This line also works, but for values lower or higher than –10 or +10, the color ramp "plateaus". And this line only works because the extents of this particular mapping happen to go from –29 to +17. Is it possible to obtain the values for the full range and specify them for a custom palette? Or is it possible to assign a custom palette and call it like the default palettes? I'm happy to clarify if needed. Thanks! 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 Office: 301-402-0095 Mobile: 301-758-3559 Web: BCBB Home Page<http://www.niaid.nih.gov/about/organization/odoffices/omo/ocicb/Pages/bcbb.…> Twitter: @niaidbioit<https://twitter.com/niaidbioit> 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.

2 2

Modeling with cryo-EM
by Aurore Vaitinadapoule 13 Dec '13

13 Dec '13

Dear all, I'm would like tu use this following script to model with cryo-EM: from modeller import * log.verbose() env = environ() struct='./M60_9.pdb' map='new_map.mrc' resolution=10.2 box_size=48 # how to define the box size ?? apix=1.88 x=90; y=63; z=23 #origin steps=20 # Read in cryo-EM density map den = density(env, file=map, em_density_format='MRC', voxel_size=apix, resolution=resolution, em_map_size=box_size, density_type='GAUSS', px=x,py=y,pz=z) # Fit the PDB file into the map by MC simulated annealing den.grid_search(em_density_format='MRC', num_structures=1, em_pdb_name=struct, chains_num=[1], start_type='CENTER', number_of_steps=steps, angular_step_size=30., temperature=100., best_docked_models=1, translate_type='RANDOM', em_fit_output_file='modem.log') But i got this error message: _modeller.ModellerError: read_mrc_h_E> Density map must be cubic, and of dimension em_map_size (48): actual dimensions (40, 40, 36) What is exactly the box size? And how to define it? Thank's Aurore -- Aurore G. VAÏTINADAPOULÉ aurore <a.vaitinadapoule(a)inserm.fr>.vaitinadapoule(a)inserm.fr<a.vaitinadapoule(a)inserm.fr> a.vaitinadapoule(a)gmail.com PhD Student in Bioinformatics INSERM UMR-S665, DSIMB Team, Dynamique des Structures et des Interactions des Macromolécules Biologiques Mobile phone: +33(0) 778 18 52 07

2 1

chimera slow
by Daniel Murin 12 Dec '13

12 Dec '13

Hi, I have been trying to use the Unit Cell function on Chimera. It works, but the graphic is so large that moving it around is nearly impossible because the program is running so slowly. It is rendered in ribbon. Is there a way to lower the quality in order to manipulate the graphic faster? Thanks! Charles Daniel Murin Ph.D. Candidate - Ward/Ollmann Saphire Labs Integrative Structural and Computational Biology Immunology and Microbial Science Kellogg School of Science and Technology The Scripps Research Institute 10550 North Torrey Pines Road, CimBio 105 La Jolla, CA 92037 Office: 858.784.3326 Lab: 858.784.7170 Fax: 858.784.8218 email: cmurin(a)scripps.edu

2 1

(no subject)
by Danilo Di Maio 12 Dec '13

12 Dec '13

danilo.dimaio(a)gmail.com

1 0

Launch external program from Chimera python script?
by Hurt, Darrell (NIH/NIAID) [E] 11 Dec '13

11 Dec '13

Hi Chimera friends, I am thinking of writing a script for Chimera that allows the user to map lipophilicity to the surface of protein. There are a few algorithms out there, but the easiest that I've found to use is "pyMLP": http://code.google.com/p/pymlp/ I see two ways to do this: (1) Write a wrapper script (shell, Python, etc.) to first run pyMLP to generate a DX file and then launch Chimera with its own script to generate the surface and do the mapping. (2) Write a Chimera Python script to launch pyMLP to get the DX file, generate the surface, and do the mapping. For my purposes, it might be easier in the long run to choose option 2, but that seems a little trickier starting out. Can you advise me on the best way to proceed? Also, I know that you aren't looking to add new functionality to Chimera v1, but the generation of a lipophilicity map might be a nice feature to add to Chimera at some point and pyMLP seems relatively easy to package with Chimera… Thanks, 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 Office: 301-402-0095 Mobile: 301-758-3559 Web: BCBB Home Page<http://www.niaid.nih.gov/about/organization/odoffices/omo/ocicb/Pages/bcbb.…> Twitter: @niaidbioit<https://twitter.com/niaidbioit> 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.

2 3