- Chimera-users - RBVI Mailing Lists

Suggestion on biological assembly load
by Miguel Ortiz-Lombardia 11 Apr '07

11 Apr '07

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Dear all, I have small suggestion to make that I would find helpful if implemented in chimera: I think it would be great if the "biological assembly" could be fetched by ID from the File menu, just as it goes now for the pdb coordinates file from the RCSB site. Such assemblies are kept at the PQS database at the EBI (Cambridge). The links there follow this syntax: http://pqs.ebi.ac.uk/pqs-doc/macmol/<ID>.mmol where <ID> is obviously that of the PDB structure. Despite the extension, it is a plain PDB file. The CRYST1 record (for crystallographic structures) of such coordinates files has been modified to space group P1, since the crystallographic operations of the original structure don't hold. A potential problem if these files is that they don't seem to care about having duplicated atom numbers... These assemblies are most often correct and validated by experiments. However, you might want to add a caveat or, even better, an intermediate window with a link to the place where the PQS explains how that particular oligomer was generated and why it is considered biologically relevant: http://pqs.ebi.ac.uk/pqs-bin/macmol.pl?filename=<ID> Best regards, Miguel - -- Miguel Ortiz Lombardía Centro de Investigaciones Oncológicas C/ Melchor Fernández Almagro, 3 28029 Madrid, Spain Tel. +34 912 246 900 Fax. +34 912 246 976 e-mail: molatwork(a)yahoo.es - ---------------------------------------------------------------------- Et ainsi ne pouvant faire que ce qui est juste fût fort, on a fait que ce qui est fort fût juste. Blaise Pascal, Pensées -----BEGIN PGP SIGNATURE----- Version: GnuPG v1.4.5 (Darwin) iD8DBQFGGL+gF6oOrDvhbQIRAtQ2AKCi9oIpIJffYwshk1PgB2aiG4e9rgCghjzZ Uz0yTug1VXEHXNP2MG2HsmI= =EGzy -----END PGP SIGNATURE-----

4 3

exporting a mineral structure as VRML
by Janet Iwasa 11 Apr '07

11 Apr '07

Hello, I'm trying to use Chimera to export a mineral file ( http://virtual-museum.soils.wisc.edu/soil_smectite/soil_smectite.pdb) as a ball and stick VRML structure, but I'm running into a bit of difficulty since some of the bonds don't get made into sticks -- rather, they appear as purple dotted lines in the viewer. Any suggestions on how to get these bonds to be made into geometry? thanks! Janet Iwasa

4 3

Extra Chimera marker attributes
by Tom Goddard 10 Apr '07

10 Apr '07

Hi Jeff, The Chimera path tracer markers accommodate extra attributes that are read and written to the XML marker files. There is no user interface to define new marker attributes in Chimera. A user interface to define new marker attributes does not seem hard. You could specify an attribute name and value in a command or dialog that assigns that attribute for all selected markers. For these new attributes to be of any use within Chimera more capabilities would be needed. For example some command or dialog that would allow you to select all markers with a certain attribute and value (e.g. is_filament = true). These kinds of capabilities already exist for attributes of atoms of molecules. For example, see the "setattr" and "defattr" commands, and "Render/Select by Attribute" tool, and atom specifiers that can test for specific attribute values. http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/setattr.html http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/defattr.html http://www.cgl.ucsf.edu/chimera/docs/ContributedSoftware/render/render.html http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/atom_spec.html#descri… These will work with path tracer markers since markers are implemented as atoms. But there are problems. For example "defattr" uses names of atoms and the names of the path tracer atoms are not meaningful. Probably "setattr" is more useful and you can set an attribute for all currently selected atoms. The main weakness is that these atom attributes will not be saved in marker XML files. That is probably easy to fix but I would need to investigate how atom attributes are implemented. For markers the attributes are kept in a Python dictionary marker.extra_attributes. I do not know if atom attributes on markers will be saved in Chimera session files. Would have to test that. Tom Jeff wrote: > Tom, > > how hard would it be to add other attribute(s) to the markers > generated with the volume path tracer tool, aside from > marker-to-marker connectivity? Such as the ability to enumerate > individual filaments/objects within a single marker set, or a flag > that would tell whether a set/group is composed of filaments, or > surfaces, or just a point cloud? Either of of these (enumeration, or > object type) could be stored as an integer. > > -Jeff

1 2

transparency factor
by jean-francois menetret 09 Apr '07

09 Apr '07

Is it possible to change the transparency factor of a volume at the command line ? Jean-Fran�ois M�n�tret, PhD Boston University School of Medicine Physiology and Biophysics Department 700 Albany Street W315 Boston, MA 02118 Email: menetret(a)bu.edu Mailing address: 715 Albany Street

3 3

Fetch EMDB map by id
by Thomas Goddard 09 Apr '07

09 Apr '07

Hi Cathy, It is a good idea to allow fetching EMDB maps by id in the current File / Fetch by ID dialog. Their are two issues related to the large size of some maps that I would like to solve to make that work well. First I would like Chimera not to freeze while the fetching is done. That is what currently happens when fetching files. Instead it could show a progress dialog with a cancel button (like a web browser) and open the map when it finally arrives. Second I would like the map to be saved on your local machine so that you don't have to fetch the same large file more than once. Perhaps the way that would work is the the first time you fetch any file it asks you where to save this and future fetched files, and there will be an option not to save them. This would be remembered in your Chimera preferences file. When you subsequently used the fetch dialog it could first look in the directory where you save fetched files and use the copy there if available. Thanks for the suggestion. It is on my requested features list. Tom

2 1

Re: [Chimera-users] hkcage command
by Thomas Goddard 09 Apr '07

09 Apr '07

Hi Padma, The code for the hkcage command had to be changed for Chimera version 1.2348 and newer because we switched from using the Numeric Python array package to its not completely compatible successor called NumPy. The new code for the hkcage command is on the experimental features web page http://www.cgl.ucsf.edu/chimera/experimental/experimental.html Use that new code with Chimera 1.2363. I have not had time to try flattening an icosahedral map. I made a small step toward another request of yours, plotting radial density. I tried a Python 2-d plotting library called matplotlib that could be used for making radial density plots. It is a pretty big package so we have not decided whether to include it in Chimera. Tom

1 0

GIF animations in Chimera
by Tom Goddard 09 Apr '07

09 Apr '07

Anyone interested in GIF animation output from Chimera? Their main appeal is that they play directly in most web browsers (like banner ads). It was easy to add to the Chimera movie recorder using ffmpeg which is included in the Chimera distribution. But ffmpeg uses a fixed color palette causing severe color banding and does not offer LZW compression resulting in large file sizes. That option looks unusable. I tried another approach using a separate program called "gifsicle". That gave decent coloring and compression but it would require adding gifsicle to the Chimera distribution. I am reluctant to do that (entails future maintenance) unless there is strong interest in making GIF animations. Tom

2 1

[Fwd: Re: GIF animations in Chimera]
by Thomas Goddard 09 Apr '07

09 Apr '07

Miguel uses ImageMagick to make GIF animations from Chimera movie frames. Tom -------- Original Message -------- Subject: Re: [Chimera-users] GIF animations in Chimera Date: Sat, 07 Apr 2007 09:09:43 +0200 From: Miguel Ortiz-Lombardia To: Tom Goddard Hi Tom, I make animated GIFs with Chimera by simply saving the frames and using convert (from ImageMagick). I run an extremely simple script to do so and the results are up to my needs. Cheers, Miguel En/na Tom Goddard ha escrit: > Anyone interested in GIF animation output from Chimera? Their main > appeal is that they play directly in most web browsers (like banner ads). > > It was easy to add to the Chimera movie recorder using ffmpeg which is > included in the Chimera distribution. But ffmpeg uses a fixed color > palette causing severe color banding and does not offer LZW compression > resulting in large file sizes. That option looks unusable. I tried > another approach using a separate program called "gifsicle". That gave > decent coloring and compression but it would require adding gifsicle to > the Chimera distribution. I am reluctant to do that (entails future > maintenance) unless there is strong interest in making GIF animations. > > Tom > > _______________________________________________ > Chimera-users mailing list > Chimera-users(a)cgl.ucsf.edu > http://www.cgl.ucsf.edu/mailman/listinfo/chimera-users >

1 0

Re: [Chimera-users] [chimera-dev] Building a multimer model
by Elaine Meng 07 Apr '07

07 Apr '07

Hi Benjamin, There is probably more than one way to do this in Chimera, but I'd open the protein of interest six times and match it to the six monomers of the related protein. Let's call your protein "A" and the similar protein "B" (for which a hexamer is available or can be easily constructed from matrices in the PDB file). Since you didn't mention the specifics, I didn't try the whole process outlined below (not having a good example to work with), but have used the tools in other contexts. This is kind of a brute force approach. Someone else may know of a more elegant way... Making the hexamer of B: If the PDB file doesn't already have a hexamer, it should have matrix information that allows it to be constructed. If BIOMT, use Multiscale Models (under Tools... Higher-Order Structure). Initially, this doesn't copy all the atoms but just makes surfaces for computational expediency. Since you will need copies of all the atoms for matching, using the Multiscale Models GUI: select "All" chains, change their style to anything other than surface (e.g. Show... Ribbons). You can see the model IDs of the copies of B in the Model Panel (under Favorites). If SMTRY, CRYST1, or MTRIX you could use the Unit Cell tool (under Tools... Higher-Order Structure) instead of Multiscale Models. Relevant manual pages: http://www.cgl.ucsf.edu/chimera/docs/ContributedSoftware/multiscale/ framemulti.html http://www.cgl.ucsf.edu/chimera/docs/ContributedSoftware/unitcell/ unitcell.html Matching the 6 copies of A to the monomers of B: Options for matching are the "match" command or MatchMaker (GUI under Tools... Structure Comparison, also implemented as a command, "matchmaker"). If you wanted to name specific atoms to use for fitting, you'd use the "match" command. That can be tedious, and I'd try MatchMaker first - if there is some sequence similarity, even fairly low, it usually works quite well. By default, it iterates the fit so that only residues that match well in space are ultimately used to superimpose the proteins. That means the parts that are different should not interfere with a good fit of the parts that are the same. You can adjust several parameters, but for example if B is model 0, chains A-F and copies of A are models 1-6, these commands would use the matchmaker defaults: matchmaker #0:.a #1 pair ss matchmaker #0:.b #2 pair ss matchmaker #0:.c #3 pair ss matchmaker #0:.d #4 pair ss matchmaker #0:.e #5 pair ss matchmaker #0:.f #6 pair ss (pair ss just means to use the specified chains rather than trying all chain combinations between the models) You could specify other parameters with various keywords, for example the stringency to use when iterating the fit, how much weight secondary structure should have vs. the residue types in the sequence, etc. Relevant manual pages: http://www.cgl.ucsf.edu/chimera/docs/ContributedSoftware/matchmaker/ matchmaker.html http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/matchmaker.html http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/match.html I hope this helps, Elaine ----- Elaine C. Meng, Ph.D. meng(a)cgl.ucsf.edu UCSF Computer Graphics Lab and Babbitt Lab Department of Pharmaceutical Chemistry University of California, San Francisco http://www.cgl.ucsf.edu/home/meng/index.html On Apr 5, 2007, at 12:10 AM, Benjamin G.S. Horsman wrote: > Hello, > > My name is Benjamin Horsman and I am a student at the University of > Toronto. I have what should be a basic (possibly naive) question > regarding protein modelling that I was hoping you might be able to > answer for me or point my in the right direction. I have been > studying the structure of a bacterial protein complex, trying to > understand how it interacts with several known binding partners. The > complex is known to be a hexameric homomultimer. Recently the > structure of the subunit protein was solved. However, only the > monomeric form of the protein was determined, and the corresponding > PDB file does not contain coordinates for the multimer (i.e., no > BIOMT1, BIOMT2, BIOMT3 matrix values). I would like to assemble a > model of the hexameric complex for visual inspection and for some > docking studies. What would be the best way to accomplish this? > > There are several structures in the PDB of hexameric protein complexes > composed of subunits with very close structural alignment to my > protein of interest, differing in two regions of the protein that are > outside of the multimer binding regions. I believe that my protein > complex posses very similar or nearly identical geometry to these > known complexes. Is there any way to basically insert my protein into > this geometry? Can this be done in Chimera? If not, is there any > program or web server that allows this to be done? Or are you aware > of a better way to get this done? > > Thank you for your time. Any help is greatly appreciated. > > > Benjamin Horsman > _______________________________________________ > Chimera-dev mailing list > Chimera-dev(a)cgl.ucsf.edu > http://www.cgl.ucsf.edu/mailman/listinfo/chimera-dev

2 1

Re: [Chimera-users] viewing & comparison of multiple overlay structures
by Roberta King 04 Apr '07

04 Apr '07

Thank you Eric for your very detailed answers to my several questions. It will take me a few days to work through the suggestions. If I have any more troubles, I will contact you next week. Thanks very much, Roberta King ******************************** Roberta S. King, Ph.D. Department of Biomedical and Pharmaceutical Sciences College of Pharmacy University of Rhode Island http://www.uri.edu/pharmacy/faculty/king _____ From: Eric Pettersen [mailto:pett@cgl.ucsf.edu] Sent: Tuesday, April 03, 2007 7:57 PM To: Roberta King Cc: Chimera BB Subject: Re: [Chimera-users] viewing & comparison of multiple overlay structures On Apr 3, 2007, at 1:47 PM, Roberta King wrote: First, After overlaying multiple protein (pdb) structures with ligands using MultiAlign based on alignment, how does one easily view only the overlaid ligands? [Note: I have up to 36 structures overlaid, each with 0-2 ligands.] Hi Roberta, Using the command line (Favorites->Command Line), typing "~disp; show ligand" will hide all structures and then show only the ligands. If every structure had a ligand you would only need "show ligand" since the show command hides parts of structures that aren't to be shown, but if no part of a structure matches the show criteria, the structure isn't changed -- so you need the "~disp" also to guarantee that ligand-less structures get hidden. This is also doable with the menus by Select->Structure->ligand followed by Select->Invert Selection (or shift-right-arrow instead to invert the selection in all models) and then Actions->Atoms/Bonds->hide. From there, you will probably want to select the ligand clump by control dragging the mouse across the clump (there will be a green box outline as you drag the mouse) and then use Actions->Focus to zoom in on the clump and make it the center of rotation. Secondly, can chimera then calculate a solvent-accessible closed surface for the protein cavity surrounding each ligand? Alternatively, I could calculate the cavity first, then overlay based on alignment. The goal is to view the overlaid cavities for visual comparison. It depends. If the cavity is actually completely enclosed then Chimera won't depict it. If there is some channel to the binding pocket, then it should be depicted. I say "should" because the surfacing library we are currently using (MSMS) is known to have problems when there are narrow channels in a surface. We are working on a replacement library that should solve both these problems, but it's not ready yet. Providing that MSMS doesn't fail, your pockets should be depicted. It is also sometimes possible to work around the MSMS problem by adding hydrogens to a structure or by slightly changing VDW radii (with the vdwdefine command) -- but you have to do this before any surfaces fail because after a failure surfacing won't work again until you restart Chimera. You create a surface with Actions->Surface->show or the "surf" command. If with experience you know that surfacing certain models is problematic, you can restrict the surfacing to particular models by selecting them before using the menus or providing the model numbers as arguments to the surf command (e.g. "surf #0,3-7,9"). Thirdly, can chimera calculate electrostatic charge on amino acids surrounding a binding cavity and reflect that charge onto the closed cavity surface? Can this surface be included in overlay view? You can use the Add Charge tool (in the Structure Editing category) to add partial charges to all atoms. You can then use the Render By Attribute tool (Structure Analysis category) to color the surface and/or atoms based on the partial charges. For more precise electrostatics, you can use an external tool such as APBS or Delphi to generate an electrostatic grid, and open that in Chimera and then use the Electrostatic Surface Coloring tool (Surface/Binding Analysis category) to color the surface with the grid data. You can certainly show the surface(s) with the ligands overlaid. You would probably only want to show a few surfaces at most simultaneously for understandability (see next answer). Given the number of structures, I would much appreciate any advice on how to make this manageable. If Chimera cannot do the calculations, can it import multiple structures saved by another capable software, then overlay them using MultiAlign (including the surfaces)? Depending on the speed of your computer, its amount of memory and its graphics card, working with 36 structures and their surfaces at once may or may not be satisfactory (and unless you have a real high-end machine my guess would be not). You will want to use the Model Panel (under Favorites) to control which surfaces/structures are being shown -- particularly which surfaces since I imagine that looking at more than a few would be incredibly confusing. Another thing that will help is to delete parts of structures that you aren't interested in. One trick here is to use MultAlign Viewer to select one residue in each structure (create a region by mouse dragging that covers one column in the alignment -- those residues will be selected), then select the corresponding chains by clicking into the main window and hitting the up-arrow key, then right-arrow to invert the selection, and finally Actions->Atoms/Bonds->delete to delete all atoms that are not part of chains in the alignment. Hopefully your ligands will have the same chain ID as the chain in the alignment and therefore will not be deleted. You might want to look at your ligands after you've inverted the selection and make sure none are selected for the following deletion. If some are selected, deselect them with "~select ligand". As a quasi-extreme measure you could prune away parts of the structure more than a certain distance from the ligands. Select the ligand clump and use Select->Zone... to select residues further than a certain distance from the clump and then delete the selected residues. --Eric Eric Pettersen UCSF Computer Graphics Lab pett(a)cgl.ucsf.edu http://www.cgl.ucsf.edu

1 0