- Chimera-users - RBVI Mailing Lists

1.9 production release
by Eric Pettersen 13 May '14

13 May '14

Hello everybody, A new production release of UCSF Chimera (version 1.9) is available: http://www.rbvi.ucsf.edu/chimera/download.html Download is free for noncommercial use. Platforms: Windows, Mac, Linux, including 64-bit versions (recommended for working with large datasets on machines with at least 4GB memory). New features include multiple sequence alignment using Clustal Omega or MUSCLE web services, building double-helical nucleic acids, "colorkey" command, "struts" command to reinforce structures for 3D printing, "vseries" volume series playback and processing options, more efficient save/restore of coordinates in sessions, Dynameomics amino acid rotamer library, WMV2 movie output, COLLADA export. More details are given below; see release notes for the full list: http://www.rbvi.ucsf.edu/chimera/docs/relnotes/1.9.html Display, Interface: * icon for menu selection mode: Append, Intersect, Replace, Subtract * distance monitor line style and color can be saved as preferences * silhouette outlines can be controlled per-model * Axes/Planes/Centroids objects can be recolored with Actions menu * "colorkey" command implementation of Color Key tool Structure Analysis and Modeling: * new Align Chain Sequences tool for multiple sequence alignment of structure chains using a Clustal Omega or MUSCLE web service (these can also be called to realign the sequences in an existing Multalign Viewer window) * "mda" command to Blast-search and retrieve structures similar to a query sequence, lay them out left->right by N->C along the query * APBS interface omits explicit solvent by default * "coulombic" command allows specifying atoms/charges to use for ESP calculation independent of the surface to color * Build Structure tool allows generating nucleic acid double helices, specifying phi angle when joining peptides * "adjust" command to change bond angle or length * MD Movie session support for partly loaded trajectories * Morph Conformations pairs chains by ID if the sets of IDs are same * "ramachandran" command can assign probabilities as attributes * Rotamers/swapaa include Dynameomics rotamer library (Daggett group) Volume Data (density maps, ESP maps): * "vop" volume operations to take pointwise maximum or multiply maps * "vseries" command options to open and close volume time series, jump playback to a specific frame, perform various processing such as alignment, normalization, cropping, compressing, save to file * new volume input format VTK structured points (ASCII) General, I/O: * some Preferences reorganization, new categories: Labels, Selection * more efficient save/restore of coordinates in sessions * new output movie format WMV2 * new export formats COLLADA, VTK * "struts" command to add strengthening bonds for 3D printing * "leap" command for using Chimera with a Leap Motion Controller Enjoy! On behalf of the Chimera team, --Eric Eric Pettersen UCSF Computer Graphics Lab http://www.cgl.ucsf.edu

1 0

how to get translation and rotation vectors for move and turn commands
by Jan Kosinski 13 May '14

13 May '14

Hi, I have the same model loaded in Chimera twice. One model has been moved to a different position and rotated manually using mouse. How can I get x,y,z coordinates for the centers of the to models to calculate x,y,z axis definition compatible with the move command? And how to get rotation information for turn command? I tried with cofr command but it seems to give coordinates in different reference frame. I would like to avoid using savepos/reset/fly commands in this case. Python solution would be fine. Thanks in advance for your help, Jan

3 4

EM density map rescaling issue???
by Christian Geiss 13 May '14

13 May '14

Dear Chimera User list, I have a fundamental problem when loading pdb structures and EM density maps as *.em. The scale doesn't fit at all, with other words my pdb structure is tremendously bigger than the EM structure. One work-around is to resample the EM structure from current 64 px boxsize with 8.6 A pixelsize to 1 A pixelsize, giving ~ a 500 pixel volume, but then the volume gets quite big und my chimera gets really slow in handling those structures (pdb and big EM map at once)...is there any other way to solve this issue like scaling the EM structure on the fly in not such a memory consuming way...I also wonder if there is any information about the pixelsize in the EM header, because then the scaling should be done usually automatically??? Thanks in advance, Chris P.S. I couldn't find anything related to this in the chimera archive... -- *Christian Geiss, M.sc.* Universitaet Frankfurt BMLS Frangakis Group Max-von-Laue-Strasse 15 D - 60438 Frankfurt am Main Phone: +49-(0)69-79842597 Fax: +49-(0)69-79846467 geiss(a)biophysik.org

3 2

threshold level problem
by 李智海 13 May '14

13 May '14

Hello, I am a freshman for Chimera, sorry for a possibly stupid question but i seem to be confused by the difference between threshold level in volume view and the sigma value that is used to contour the density map in some articles. please help me out! Thank you very much! Zhihai Li

3 2

help with surface coloring
by Benda, Christian 09 May '14

09 May '14

Hello, sorry for a possibly stupid question but I seem to have a problem getting a molecular surface colored the same way as the atoms or ribbons of the underlying structure. Here’s what I did: - selected two domains of a larger protein separately (sequence viewer) and colored each domain differently (say green and blue, ribbons and atoms are in the respective color) - then selected the molecule in the model panel and clicked surface main - I get a surface, but it is either transparent grey or a solid color (pink) of the color well displayed in the model panel (not the domain colors I selected). - when opening the attributes of the surface, I can toggle between color source “atoms" or “model" but that only toggles between transparent gray and pink. I tried all kinds of combinations and clicks but I don’t manage to get the surface colored in the colors of the two domains. Any hint? It would be highly appreciated! Many thanks! Christian

2 3

Re: [Chimera-users] Display a spider volume data in a specific orientation
by Elaine Meng 08 May '14

08 May '14

Hi Zongli, This is not my area of expertise (our expert is currently out of town) -- at first I thought you could use "Transform Molecule Coordinates" (under Tools... Movement), but now I see that it works only on atomic structures rather than maps: <http://www.cgl.ucsf.edu/chimera/docs/ContributedSoftware/transform/transfor… > It only works on molecular structures because it actually changes the atomic coordinates, without changing the model transformation. I've figured out one workaround, but it is both crazy and slightly tedious: (1) open your map (say as model 0). (2) open any atomic-resolution structure, it doesn't matter what, and open it twice -- now you have three models. Don't move them separately from the map (yet). Let's say the atomic structures are models 1 and 2. (3) use "Transform Molecule Coordinates" with desired Euler angles and shift on one copy of the atomic structure (say model 1). (4) use the match command to superimpose the other copy of the atomic structure on the one that was transformed: match #2 #1 (the command can use the whole models since they should have exactly the same atoms) (5) use the matrixcopy command to apply the transformation matrix of the second copy of the atomic structure to the map: matrixcopy 2 0 You might think "why do I need two copies of the atomic structure? can't I just transform one copy and then use matrixcopy?" ... in fact that's what I tried first, but after it didn't work, I realized it was because the Transform tool works by changing the coordinate values and not the stored matrix of the model. It is necessary to match the second copy to generate a corresponding matrix, and then matrixcopy can be used. I'm CC-ing the chimera-users(a)cgl.ucsf.edu list in case other Chimera users have suggestions. Our expert Tom G should be back next week and he might also have ideas. In general, it's better to send Chimera questions to that mailing list instead of just me, because I could be away from work or (like now) lacking a good answer! 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 Apr 7, 2010, at 11:51 AM, Li, Zongli wrote: > Hi Elaine, > I have a spider volume data (or mrc density map) from which I got a > serious 2D projections. I have the Euler angles for each projection. > Is there a way in Chimera that I can use those Euler angles to > display the 3D spider volume data (or mrc density map) in a specific > orientation that corresponding to that 2D projection? > Thank you very much in advance. > Best, > Zongli > > **************************************** > Zongli Li, Ph.D. > > Howard Hughes Medical Institute > Department of Cell Biology > Harvard Medical School > C Building, Room 420B > 240 Longwood Avenue > Boston, MA 02115 > > Tel: 617-432-4095 (Office) > 617-432-4708 (EM Suite) > 617-432-4091 (Computer room) > Fax: 617-432-1144 > Email: zongli_li(a)hms.harvard.edu > *************************************** >

3 5

Enquiry: Match Maker
by Anna Hernandez 06 May '14

06 May '14

Hi, I have a question regarding structural alignments with Match Maker. I would like to know if there is any way to obtain a list of all possible alignment solutions (or a subset) tested by Match Maker, instead of only the "best" one, so that I can further assess the different options manually. Specifically I am trying to find an homologous substructure/motif from my template structure within a query structure. The pattern to be searched is rather small as compared to the whole protein and it can be aligned in multiple positions. I would like to assess each of the possibilities more carefully based on the biological knowledge of the protein. Thank you in advance. Best wishes, Anna Hernandez

2 1

Re: [Chimera-users] how to save distance to a file
by Eric Pettersen 05 May '14

05 May '14

On May 5, 2014, at 4:07 PM, chang YG <changyg2008(a)gmail.com> wrote: > Hi Eric, > > Thank you very much for the quick response. I can't wait to try the script you provided. > > For the RMSD, I have two conformers for the same sequence. If I encounter any problem, I will give it my best shot before posting questions. In that case it should be easy: from chimera import openModels, Molecule m1, m2 = openModels.list(modelTypes=[Molecule]) import Midas rmsd = Midas.rmsd(m1.atoms, m2.atoms) etc. --Eric > > Yonggang > > > On Mon, May 5, 2014 at 3:59 PM, Eric Pettersen <pett(a)cgl.ucsf.edu> wrote: > On May 5, 2014, at 3:57 PM, Eric Pettersen <pett(a)cgl.ucsf.edu> wrote: > >> print>>f, a1, a2, diet > > Stupid auto-correct spelling. "diet" should be "dist". > > --Eric > >

1 1

how to save distance to a file
by chang YG 05 May '14

05 May '14

Hi, Thanks for taking time to answer the questions. Could you advise me of how to save distances using a python script. Here is what I have: from chimera import runCommand as rc rc("open 1t4y") rc("distance :33@H :34@H") #how to save distance. Another related question: how to save "rmsd" of two structures using a python script? Thanks, Yonggang

2 2

patching ACE NME
by Francesco Pietra 05 May '14

05 May '14

Hello; In parameterizing a new metal cluster, I would like to patch the various amino acids (which stand isolated from one another) with what are called by AMBER ACE and NME. I would like to do that within the whole block, not for each amino acid isolated from the block. Also, I would like to avoid tleap, which gets into troubles here. I could do that by modifying structure, but a problem arises of residue names, i.e, the original residue takes the name ACE or NME, or viceversa the capping groups take the name of the amino acid being capped. Is there a specific tool with chimera? thanks francesco pietra

2 2