- Chimera-users - RBVI Mailing Lists

x,y coordinates of the Chimera screen
by Jan Kosinski 13 May '14

13 May '14

Hi, When I load a model in Chimera I can move the model on the plane of the screenwith the move command along the x,y axes. Is it possible to get the current range of the coordinates of the x and y axes of the screen, so I could automatically (without setting positions with savepos) move my models to particular points on the screen? For example, move to left upper corner in a movie... Or move the model outside the screen and create nice fly-in effect in movies. Thanks, Jan

2 4

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

A question on tunnelling chimera through ssh between a linux server and OSX client
by Miguel Ortiz Lombardía 01 May '14

01 May '14

Hi all, I understand that this is not a typical use of chimera, so I don't expect much support, but I would appreciate if someone can give me a clue to solve the issue... We have a GNU/Linux server (Debian wheezy) with some EM software installed on it. One of these programs calls at some point chimera to show something on which the user has to act upon (or just check, not sure, I'm not the user as you've guessed) This machine is accessed from client OSX machines via ssh. All permissions are set so X11 is seamlessly tunnelled between them (checked with many graphical programs) Now, for a OSX client running on 10.6.8 (snow leopard) and using Apple's X11 the whole process works without problems. But for a OSX on 10.7.5 (lion) with a parallel installation of Apple's X11 and Xquartz 2.7.5 the first small window opens and then the chimera process hangs (it can be killed via Ctrl-C, which returns the control to the program that launched chimera and things can proceed, though without the information provided by chimera) If I connect to the Linux server from the lion machine and directly execute chimera I make a different observation: the first small window opens but then chimera crashes leaving this error in the terminal: Error of failed request: GLXBadCurrentWindow Major opcode of failed request: 149 (GLX) Minor opcode of failed request: 5 (X_GLXMakeCurrent) Serial number of failed request: 1360 Current serial number in output stream: 1360 Now, I would appreciate any clues on the meaning and possible solution of this error or at least any directions on how to properly remove both Apple's X11 and Xquartz (they are supposed to be able to work side-by-side but I presume that a possible source of problems could arise from some conflict between them) AND re-install them to test them in a fresh install. Downgrading the lion machine to snow leopard does not seem to be an option... Thank you! Cheers, -- Miguel Ortiz Lombardía Architecture et Fonction des Macromolécules Biologiques (UMR7257) CNRS, Aix-Marseille Université Case 932, 163 Avenue de Luminy, 13288 Marseille cedex 9, France Tel: +33(0) 491 82 86 44 Fax: +33(0) 491 26 67 20 mailto:miguel.ortiz-lombardia@afmb.univ-mrs.fr http://www.afmb.univ-mrs.fr/Miguel-Ortiz-Lombardia

3 2

Preparation of the protein-ligand system for the amber simulation
by James Starlight 30 Apr '14

30 Apr '14

Dear Chimera Users, I wounder to know if it's possible to use combination of the amber ff with the cgenff for the both parametrization of the small ligands in complex with the proteins using amber tools built in Chimera's? Does it possible to use some scripts instead of GUI for such tasks ? I'll be thankfull for some tutorial or example Thanks for help, James

3 3

multi-channel solid rendering
by Nick Walter 30 Apr '14

30 Apr '14

Hello, I've been attempting to create a movie in Chimera using multiple data files in solid style volume rendering, and I have been running into an issue where Chimera consistently shows the last channel to be loaded on top, regardless of the spacial position of the data with respect to the other solid channels. I found this message ( http://plato.cgl.ucsf.edu/pipermail/chimera-users/2011-May/006329.html) which explains that multi-channel solid rendering was removed in 2008. I was wondering, if there are any plans to reinstate this feature, or if you were aware of any work-arounds that would work well for creating movies. Thank you for your time, Nick Walter

2 1

related with Fit in map option
by SRDA.DRAKULIC.599589 30 Apr '14

30 Apr '14

Hi, I have a question related to Fit in Map option and how to evaluate subsequently given results. I have a EM map at 19 ang. resolution of protein at yet non-characterised activity state, and there are two atomic structures of homologue of that protein, at clenched and open, active state. I would like to relate that state to one of the atomic structures, thus, both of them were fitted in my EM map and for the clenched state I got following values: cross correlation: 0.8447; average map value 9.526; and number of atoms outside the contour: 1709/7621 (or 22.42%) For the open state: cross correlation: 0.8385; average map value 9.085; and number of atoms outside the contour: 2064/8124 (or 25.42%) Are those differences big enough to say with certainty that my protein is at conformation similar to the clenched one or due to the EM map resolution this can not be said. Sincerely, Srdja Drakulic

3 2

Saving ribbons as PDBs?
by Maja Divjak 29 Apr '14

29 Apr '14

Hello Chimera Gurus, Is there a way to save ribbon structures as PDBs for import into Maya? I've found only the surface structures are saved in the PDB, even though I am working with the ribbon structure in chimera initially. If there is a way to do this, I would be mightily pleased! I used to export chimera domains as vrml files and then convert them to maya files and both the ribbon and surface structures would be included. You can reply to my work email at divjak.m(a)wehi.edu.au if you prefer. Thank you so much and best wishes, Maja Divjak, PhD Biomedical Animator Walter and Eliza Hall Institute of Medical Research

2 1

select zone & rename script
by Michael Garton 29 Apr '14

29 Apr '14

Hi all, Can anyone show me a good way to select all residues that are <2Å from zinc atoms, then rename any selected CYS and HIS, to CYM and HIN respectively ... I'm trying to write a script to do this on ~300 files. Thanks!! Mike Here's what I have so far import os from chimera import runCommand as rc from chimera import replyobj # gather the names of .pdb files in the folder file_names = [fn for fn in os.listdir(".") if fn.endswith(".pdb")] # loop through the files, opening, processing, and closing each in turn for fn in file_names: replyobj.status("Processing " + fn) # show what file we're working on rc("open " + fn) rc("select :ZNB zr<2") rc("writesel "reslist + fn") rc("close all") rc("stop now")

2 1

Question About Using Chimera In Python Scripts with Other Features
by John Dougherty 28 Apr '14

28 Apr '14

Hello, I am not very computer savvy and very new to programming. I am using Windows 7, Python 2.7.3 32 bit I would like to know if it is possible to use features from Chimera in a python script. For example: from chimera import * from otherprogram import * otherprogram stuff runCommand('blah blah blah') Chimera does some things the other program does not and vice versa. My life would be made much easier if I could just integrate them into one script. I do not understand the instructions about this matter which are posted on the FAQs on the website. Sincerely, -- John J. Dougherty III Biomathematics Major Rutgers, New Brunswick (P) 856 981 9491

2 1

Re: [Chimera-users] (no subject)
by Elaine Meng 28 Apr '14

28 Apr '14

Hi Joshua, The linewidth issue is mentioned in the supersampling documentation, but you wouldn't be expected to know that your problem was related to supersampling in the first place! Most people report encounter it as the lines being too thin rather than being completely invisible. <http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/print.html#supersampling> If you use the Image Save dialog, it reports near the bottom of the dialog the "effective maximum linewidth" given the specified level of supersampling and image pixel dimensions. Elaine On Apr 28, 2014, at 3:43 PM, Eric Pettersen <pett(a)cgl.ucsf.edu> wrote: > Hi Joshua, > I'm glad to hear it worked! And of course you're welcome. Don't hesitate to send other problems or concerns to the list, and I'm sure you're right that others will benefit from having the info Googlable from the list archive. > > --Eric > > On Apr 28, 2014, at 3:39 PM, Joshua Broyde <jb3401(a)c2b2.columbia.edu> wrote: > >> Dear Eric, >> Thank you so much. The first thing I did is change the supersampling from 3x3 to 1x1 and it worked! Thanks so much! >> >> Sincerel, >> Joshua Broyde >> p.s. feel free to share this email with others who have this problem. When I googled "dashed lines not saving chimera".Thanks again! >> >> ----- Original Message ----- >> From: "Eric Pettersen" <pett(a)cgl.ucsf.edu> >> To: "Joshua Broyde" <jb3401(a)c2b2.columbia.edu> >> Cc: chimera-users(a)cgl.ucsf.edu >> Sent: Monday, April 28, 2014 2:06:04 PM >> Subject: Re: [Chimera-users] (no subject) >> >> Hi Joshua, >> The first thing you should do is make sure you have the latest driver installed for your graphics card. Based on the info you provided, I believe the latest release for your card was 4/24/14 -- only 4 days ago! Go to the AMD drivers site here: >> >> >> Download Drivers >> >> >> and input the values for your system to get to the latest drivers. >> If that doesn't help, my best guess as to the issue is that the card/driver can't show dashed lines beyond a certain thickness and that the supersampling done while saving the image requires lines thicker than the card supports (each "tile" of the supersampled image is essentially a zoomed-in depiction of that part of the image and is later sampled down to its final size -- which reduces aliasing artifacts). Therefore the two things to try are to reduce the amount of supersampling as you save the image, or reduce the thickness of the dashed lines before you save the image. >> There is a "Supersample" control on the Save-Image dialog, and there is a "supersample" keyword to the "copy" command, if that's what you're using. >> You can reduce the line width of the hydrogen bonds by bringing up the PseudoBond Panel (in the General Controls tool category) and double clicking on the "hydrogen bonds" group (assuming you've made some hydrogen bonds). The "line width" control is in the lower part of the inspector that comes up. >> I hope one of these things work! :-) >> >> >> --Eric >> >> >> Eric Pettersen >> UCSF Computer Graphics Lab >> http://www.cgl.ucsf.edu >> >> >> >> >> On Apr 27, 2014, at 3:24 PM, Joshua Broyde < jb3401(a)c2b2.columbia.edu > wrote: >> >> >> Hi, >> >> >> I had a chimera-related problem. When saving images from chimera, chimera is not saving the dashed hbonds (I am NOT using POV-Ray), in the images, even though they appear on the screen. Does this have to do with the graphics driver, or is is there something internally in chimera that is likely causing this? >> I am using the AMD Radeon HD 6450 graphics driver on ubuntu 12.04. When I tried this on another computer with ubuntu 12.10 (using Gallium .04 graphics driver) , it worked fine and the image came out perfectly. This problem happens whether I use chimera 1.8, 1.8.1, or any release after that. >> >> Sincerely, >> Joshua Broyde >> _______________________________________________ >> Chimera-users mailing list >> Chimera-users(a)cgl.ucsf.edu >> http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users >> >> > >

1 0

(no subject)
by Joshua Broyde 28 Apr '14

28 Apr '14

Hi, > I had a chimera-related problem. When saving images from chimera, chimera is not saving the dashed hbonds (I am NOT using POV-Ray), in the images, even though they appear on the screen. Does this have to do with the graphics driver, or is is there something internally in chimera that is likely causing this? > I am using the AMD Radeon HD 6450 graphics driver on ubuntu 12.04. When I tried this on another computer with ubuntu 12.10 (using Gallium .04 graphics driver) , it worked fine and the image came out perfectly. This problem happens whether I use chimera 1.8, 1.8.1, or any release after that. > > Sincerely, > Joshua Broyde

2 1

Print selected
by m.maletta＠nki.nl 28 Apr '14

28 Apr '14

Dear all, I would like to select all the atoms of the chain B that are less than 3.5A far away from chain A. And I would also like to print out all the distances that Chimera finds. I manage to do the first part with this command, but I do not manage to have a list of all the distances found. select #0:.A za<3.5 & #0:.B Is there a way to get the list? Best and thank you in advance Max __________________________________ Massimiliano Maletta post-doc in NKI Netherlands Cancer Institute -Antoni van Leeuwenhoek Hospital Plesmanlaan 121 B6 1066CX Amsterdam group of Peter Peters www.nki.nl/research/peters/ | www.necen.nl m.maletta(a)nki.nl Tel : +31 20 512 2022 Fax : +31 20 512 2029

2 2

Multiscale models color
by Hernando J Sosa 26 Apr '14

26 Apr '14

Hi, How do you control the color of the surfaces generated by the multi-scale models command "Make models" so that the user can specify the color associated with each chain? I tried to follow the documentation on how to do this but still I don't get it right as I didn't manage to change the colors assigned by default. Thanks Hernando

2 2

1.9 release candidate
by Eric Pettersen 24 Apr '14

24 Apr '14

Hi everybody, There is a release candidate of Chimera 1.9 on the download page. If you can, please try it out and let us know if you have any problems or questions. The differences between 1.9 and 1.8 are covered in the release notes. Thanks! --Eric Eric Pettersen UCSF Computer Graphics Lab http://www.cgl.ucsf.edu

1 0

Export Pov-Ray file with a python script.
by jiserte＠unq.edu.ar 23 Apr '14

23 Apr '14

I'm writing a small python script to generate Pov-Ray with chimera. The script loads a pdb file containing a protein and a ligand, modifies the protein view to ribbon and the ligand to surface. Then exports the .pov file. However the image created by Pov-Ray shows only the surface of ligand. The image is shown correctly in the chimera GUI. If I run the commands one by one in the chimera command line, the .pov file generated is ok. The script is: ----------------------------------------- import chimera from chimera import runCommand runCommand('windowsize 800 600') runCommand('open ../../kinase.domain.2/A.4K9Y.pdb') runCommand('ribbon #0') runCommand('surf ligand') runCommand('export test.pov') ------------------------------------------ After running the script the ReplyLog contains: ----------------------------------------- A.4K9Y.pdb opened /opt/UCSF/Chimera64-1.7/bin/mscalc 1.400000 2.000000 1 MSMSLIB 1.3 started on hamal Copyright M.F. Sanner (March 2000) Compilation flags Surface A.4K9Y.pdb, category ligand, probe radius 1.4, vertex density 2 1 connected surface components Total solvent excluded surface area = 436.845 Total solvent accessible surface area = 772.871 ----------------------------------------- Javier Iserte

3 2

Transparency of selected atoms
by Chinaza Egbuta 23 Apr '14

23 Apr '14

Hi, Is there a way to change the transparency of selected atoms in a model? I am trying to visualize what is behind the plane of my model without using the surface capping function. Thanks! CE Chinaza Egbuta Ph.D. Candidate Department of Pharmacology & SUNY Upstate Cancer Research Institute SUNY Upstate Medical University Syracuse, NY 13210

2 1

Re: [Chimera-users] [chimera-dev] read file by python script and color pdb file with script
by Eric Pettersen 22 Apr '14

22 Apr '14

Hi Sariya, Your "check" file contains this (one) line: 67 23 which means that the residue range you generate is ":67-23" which contains nothing. You want "check" to contain "23 67" so that your residue range becomes ":23-67". --Eric Eric Pettersen UCSF Computer Graphics Lab http://www.cgl.ucsf.edu On Apr 21, 2014, at 8:33 PM, Sanjeev Sariya <s.sariya_work(a)ymail.com> wrote: > > I am having 2 coordinates in check file. I want to take/read them and paint my pdb file loaded in chimera. > I have chimera in my ~/bin/ > System specs: ubuntu 64 bit, 13.10, 4gb ram > > Below is the command I have used for script-chimera communication > > $ chimera64-1.8.1 --script "y.py check" -- Q6UXT6.pdb > > > I am storing values in listu and want to use them for coloring my file loaded. > > rc("color blue :"+listu[0]+"-"+listu[1]) > > This works well, if hard coded > rc ("color red :55-86") > > I do not get any error. The particular protein length is not colored. > > attached 3 files. > Kindly suggest, how do I proceed. > > > -- > Regards > Sariya > > <Q6UXT6.pdb><y.py><check>_______________________________________________ > Chimera-dev mailing list > Chimera-dev(a)cgl.ucsf.edu > http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-dev

1 0

read file by python script and color pdb file with script
by Sanjeev Sariya 22 Apr '14

22 Apr '14

Hello Chimera team, I am having 2 coordinates in check file. I want to take/read them and paint my pdb file loaded in chimera. I have chimera in my ~/bin/ System specs: ubuntu 64 bit, 13.10, 4gb ram Below is the command I have used for script-chimera communication $ chimera64-1.8.1 --script "y.py check" -- Q6UXT6.pdb I am storing values in listu and want to use them for coloring my file loaded. rc("color blue :"+listu[0]+"-"+listu[1]) This works well, if hard coded rc ("color red :55-86") I do not get any error. The particular protein length is not colored. attached 3 files. Kindly suggest, how do I proceed. PS: I had sent the same mail to chimera-dev(a)cgl.ucsf.edu <chimera-dev(a)cgl.ucsf.edu>; Sorry, I didn't find this mailing list before. -- Regards Sariya

1 1

Structure editing: add charge
by George Tzotzos 21 Apr '14

21 Apr '14

Hi everybody, I’m trying to add to a ligand intended for an MD simulation using Amber. The ligand is 7-octenoic acid. The charge added to the ligand is non-integer 0.999. Given that he unperturbed charge on the protein is +9.000, the charge on the complex is +9.999. Amber does not accept non-integer charges, so it adds on 9 Na+ ions. Is there a way to force the charge on the ligand to be +1.000? Apologies if this question should be addressed to the Amber list, but I used chimera to add hydrogens and charges (AM1-BCC). Regards George

3 4

Problem with charges when writing mol2 file with a python script.
by jiserte＠unq.edu.ar 21 Apr '14

21 Apr '14

I'm trying to write a simple python script to open a mol2 file with a ligand, add hydrogens, add charges and write it into another mol2. The script is this: ----------------------------------------- import chimera #1 from chimera import runCommand #2 runCommand('open ligands/lig_in.mol2') #3 runCommand('addh') #4 runCommand('addcharge') #5 runCommand('write format mol2 #0 lig_out.mol2') #6 ------------------------------------------ The result mol2 file does no contain any charges. However, if a run another script that is exactly like this but does not contain the last line (#6), and then execute the command 'write format mol2 #0 lig_out.mol2' manually, the mol2 file contains charges. I don't understand why there is a difference with both methods. best regards. javier

2 1

how to calculate the core diameter of a trimer (the individual subunits of this trimer is a helix)
by ashok rout 21 Apr '14

21 Apr '14

Hi, I have a transmembrane domain of a protein (20 amino acid long, single helix). The NMR derived structure comes out as a coiled-coil trimer. So can you please assists me how to calculate the core diameter of this trimer. Regards, Ashok -- thanks & regards, Ashok

3 2

Chimera-crashing-ubuntu
by Prasanth Kumar 18 Apr '14

18 Apr '14

Dear Sir/Madam, I have been using chimera for past 3 months. In the past two days, when i load my saved sessions (.py files), chimera is getting hanged and my ubuntu 12.04 OS is also getting hanged. So please help me on how to fix the problem. Chimera has been a really useful tool in my research. Waiting for your reply. Thank you, Yours Truly, *PRASANTH*, Undergraduate Department, Indian Institute of Science, Bangalore-560012. Contact:- prasanth.311994(a)gmail.com, prasanth(a)ug.iisc.in.

2 1

Rendering into 3d
by Rebecca Swett 18 Apr '14

18 Apr '14

Hey All, I'm looking at rendering some movies output from chimera as a 3d movie for viewing with nvidia 3d vision pro glasses. From everything I've read as long as the left/right eye frames are sequential, making the movie with ffmpeg should work fine. The problem I'm having is that ffmpeg has no way to alternate between 0-left.ppm 0-right.ppm 1-left.ppm 1-right.ppm, which is how chimera saves images when recording movies in stereo mode. If those same images were named 1.ppm 2.ppm 3.ppm 4.ppm, there would be no problem. Is there a quick fix that can be inserted to change the naming of the recorded frames? Any other suggestions? Rebecca Swett

2 1

morphs
by Moses Adenaike 17 Apr '14

17 Apr '14

Hello Chimera users, Is there a way to make a morph between two PDBs that don't have the same amount of atoms or chains. For example, I want to morph between an enzyme bound to DNA and one without DNA. Moses

2 1

call file-> open gui popup from script ?
by Katrina Lexa 17 Apr '14

17 Apr '14

Hello All, I’m wondering if there is a way to call the file->open popup within chimera when running a python script? I know that it’s possible to run the script on the command line and read the desired arguments in that way, and obviously the files I want could be coded in and then changed every time I want to look at something new, but I’m not seeing an obvious option from the mailing list or in the runCommand list for writing in a mid-script call for the user to select a file. Is such a thing possible? Thanks for your help, -- Katrina

3 3

Doubts regarding building a structure using Build Structure option
by Garg, Ayush A. 16 Apr '14

16 Apr '14

Dear Developer I was working with with UCSF Chimera. I was trying to reproduce structure 1V4S - Glucokinase structure using the Build Structure option under the tab Structure Editing in Tools. I have the sequence of the 455 amino acids and also know the sequence which is alpha helix, beta helix but the sequence of connectors between these structures is not a secondary structure. I was wondering how to reproduce those using your structure builder option. I have attached a word file with the sequence and the image of 1V4S in UCSF Chimera. I hope that will help you understand my doubt better. Thank and Regards Ayush Arpit Garg Graduate Student Ohio State University

2 1

Working with several proteins
by Navya Shilpa Josyula 16 Apr '14

16 Apr '14

Hi, I am working on a set of 400 proteins and for all these proteins I have specific list of amino acid residues. I wanted to know if there is any available script or a command to get the following 1. atomic areaSAS 2. atomic B-factor values 3. CASTp identified pocket residues and 4. CASTp identified pocket volumes for all the residues of these 400 proteins to be written in a single file. Kindly help me at your earliest. Thank you in advance for your help. Regards, Navya

2 1

Relevant interacting amino acids between molecules
by Ogunniyan, Abisola 16 Apr '14

16 Apr '14

Hi, My name is Abisola Abisoye-Ogunniyan, a first year Biology Graduate student in Tuskegee University. I got to learn about the use of Chimera in my Biochemistry Class. I have been having difficulty in trying to determine the relevant interacting amino acids between molecules. Please, I will really appreciate it if this is made clearer to me. Thank You. Best Regards, Abisola Abisoye-Ogunniyan.

2 2

Structure editing: addH
by George Tzotzos 11 Apr '14

11 Apr '14

Hi everybody, I’ve been trying to add hydrogens to a ligand (7-octenoic acid) obtained from docking. Hydrogens are added but not on the carboxyl group. I’m attaching the ligand.pdb for easy reference. I’d appreciate any suggestions regarding this behavior Best regards George

3 2

Request for guidance
by MUHAMMAD RAMZAN MANWAR 08 Apr '14

08 Apr '14

Dear Dr. Ferrin, Hope this email finds you the best of your health. First of all, I did appreciate the effort that You and Your team did in the form of UCSF Chimera. Being the part of academic/government university, I have used Chimera tool for protein modeling and cited in my two articles: one published in JCB-Wiley and other in Hundawi in 2014. Current email I am sending to request for guidance about how can I substitute an amino acid (wild type to mutant residue) within protein 3D model by using Chimera. I will be thankful for your guidance, if possible. With regards Hussain MR PACER-HD, KAU, KSA

2 1

Surprising "measure symmetry" result
by ingvar 08 Apr '14

08 Apr '14

Hello, I fetched entry EMD-2463 from EMDB and issued: measure symmetry #0 nMax 24 and got the result "Symmetry emd_2463.map: C13, center 64 64 37" Increasing the number of sampling points does not change the result, but if I increase the correlation criteria a little, to say 0.992, I get the expected C12 symmetry. Also if I change the contour level to 0.05, the author recommended level, I get C12 again. Looking at the volume at several different contour levels, it is difficult to see how it could be considered C13, and also how it would go from C13 to C12. Kind Regards, Ingvar Lagerstedt -- Ingvar Lagerstedt European Bioinformatics Institute (EMBL-EBI) European Molecular Biology Laboratory Wellcome Trust Genome Campus Hinxton Cambridge CB10 1SD Tel: +44 (0)1223 492533

2 3

Docking ligands and protein using Chimera
by Sook Wong 08 Apr '14

08 Apr '14

Good morning, I am currently doing research with a professor for protein-ligand structure and would like to dock the receptor molecule (pdb file of the protein) and ligands (pdb file of the sugars). Since I have both the structure files from Chimera, how would I predict and assume which orientation bound to each other forming a stable complex? Looking forward to hear from you. Thank you so much! Sincerely, Sook Y Wong Chemical and Biomolecular Engineering '14 A. James Clark School of Engineering University of Maryland, College Park, MD 20742 swong18(a)terpmail.umd.edu

2 1

Calculating SAS and SES with Chimera-MSMS wrap
by Jaime Rodríguez-Guerra 08 Apr '14

08 Apr '14

Dear Chimera dev team, I am trying to use your MSMS wrapper to compute some SAS/SES values for a given set of atoms. As suggested in a previous post in this list, I've taken a look into Measure.measure.buried_area() and I think I've got it, but just wanted to make sure. Here is my current procedure: ``` import chimera, Measure, MoleculeSurface atoms = chimera.selection.currentAtoms() xyzr_data = Measure.measure.atom_xyzr(atoms) surfaces = MoleculeSurface.xyzr_surface_geometry(xyzr_data) for i, a in enumerate(atoms): print "Atom: {0}, SAS: {1}, SES: {2}".format(a.name, surfaces[3][i,0], surfaces[3][i,1]) ``` Is this OK? I don't know how to handle all that output that MSMS is giving back! Thanks in advance, Jaime.

2 1

Re: [Chimera-users] Reaction Mechanism Movie
by Eric Pettersen 07 Apr '14

07 Apr '14

On Apr 5, 2014, at 1:06 PM, Moses Adenaike wrote: > Dear Dr. Pettersen, > I am wondering how I would go about making a moive detailing a biochemical reaction. For example, I want to show a DNA polymerase reaction, essentially a ligation between the backbone of two nucleotides. How would I go about doing so? I know chimera has a make new bond function but I haven't been too successful at getting that to work. Hi Moses, It is typically better to send questions like this to the chimera-users mailing list than to me directly so that others can benefit from the answer and in case others on the list have more expertise on the topic than I do. The simplest method is to decide which frame# you want the bond to start existing and use the "Define script…" entry in MD Movie's Per-Frame menu to define a one-command script to execute at that frame. For example, if you wanted a bond between :337.A@P and :338.A@O3' to exist starting with frame 53, in the Per-Frame dialog you would set "Interpret script as" to "Chimera commands" and enter as your script: #53 bond :337.A@P :338.A@O3' There is also a hidden capability in Chimera to compute all bonds on a per-frame basis. I can tell you how to access that if it turns out you really need it. However, for the purpose of making a movie it might be nice to have a more "continuous" representation of the bond forming than just the binary on/off that the above approaches give you. I have attached a Python script that you can use the Per-Frame dialog (with "Interpret script as" set to Python, and reading in the script with the "Insert text file…" button) to have a bond form (by starting out very thin and becoming thicker) as two designated atoms approach each other. There are some values at the top of the script that you can edit to make it work best in your situation. There is a distance value for when the bond begins to exist (noBond, default 3.5) and when it reaches full size (fullBond, default 2.0). There is a two-member list of residue IDs (resIDs) for the atoms involved in the bond, and a list of atom names (atomNames). Let me know if you need more help than this. Hmmm, this all assumes you have a trajectory of the reaction. You may have to use the Morph Conformations tool to create such a trajectory if you don't have one already. --Eric Eric Pettersen UCSF Computer Graphics Lab http://www.cgl.ucsf.edu

1 0

presets-interactive 3-hydrophobicity surface- those water molecules stay!
by Ying-Chieh Sun 06 Apr '14

06 Apr '14

Hi, I did a simulation using Amber package and saved a snapshot into its pdb file, and use Chimera to open it. The pdb file contains a protein-ligand complex plus a box of about 20000 TIP3P water molecules. When I pressed presets-interactive 3-hydrophobicity surface, there are 8 water molecules stayed close to the ligand binding site. I am wondering how these 8 water molecules were determined? The reason I ask is because we are trying to find a way to determine 'crystal (or say, ordered) water molecules'. This presets-interactive 3-hydrophobicity surface seems to perform this function in some way. Thank you for your help. Ying-chieh Ying-Chieh Sun, Professor of Computational Chemistry Chemistry Department National Taiwan Normal University Taipei, Taiwan 11617 Tel: +886-2-7734-6215 Fax: +886-2-2932-4249

3 2

Resolution of .stl file of ribbon representatiom
by Gavin Whittaker 04 Apr '14

04 Apr '14

Dear all, I'm in the process of 3d printing a large protein in ribbon format; I've generated the .stl file, but the resolution of the mesh that Chimera creates means that the resulting file is too large to use easily. I've managed to use other software to reduce the resolution of the mesh, so I /can/ now get around the problem, but can anyone tell me if it's possible to reduce the resolution of the mesh in Chimera to save me a step for future models? I know how to reduce the resolution of a surface mesh, but files of ribbon representations remain obstinately large. One method I've used that gives limited results is to simplify the cross-sectional form of the ribbon, but I wonder, am I missing something obvious and more effective? thanks! Gavin

3 7

query
by Bhaskar Bagchi 04 Apr '14

04 Apr '14

Respected Sir I want to calculate surface area of different atoms in a molecule. Is it possible in chimera? Please help me. Thanking you Yours faithfully, Bhaskar Bagchi Research Scholar Raiganj University College Email: bhaskarbagchi(a)gmail.com

2 1

ellipsoid main axis alignment ?
by Hugo J Bohorquez 03 Apr '14

03 Apr '14

Hi, I am writing a script in python for generating a series of ellipses linking the nuclei a molecule. Every ellipsoid should be positioned at a specific location *P* = (X,Y,Z) and should be aligned with the orthonormal vectors describing its main directions-*A*, *B*, *C*. However, when I create an ellipsoid with shape_command it gets located correctly at *P* but it is aligned with the global axis. I can't find a general way to reorient an object so its original axes point along the vectors *A*, *B*, and *C*. I will appreciate any suggestion how to do it. -- Hugo J Bohórquez <https://sites.google.com/site/hugojbohorquez/home>, PhD *Investigador* Grupo de Biofísica Fundación Instituto de Inmunología de Colombia TEL. +57 +1 -3158920 Ext. 122 CEL. +57 315-6169069 *The Pauli energy grows exponentially with the electronic localisation <http://dx.doi.org/10.1080/00268976.2013.836605>*

3 3

resampling maps
by Jason Kaelber 03 Apr '14

03 Apr '14

Hi all, I decided to do a rough alignment of density maps manually and resample onto a common coordinate system to feed the output mrc files into my fine aligner. But, when I vop resample #1 onGrid #0, I get, TypeError: Cannot cast ufunc add output from dtype('float32') to dtype('uint8') with casting rule 'same_kind' What might I be doing wrong to get this error? As far as I know, all the maps should be uint8 since they were most recently operated on by David Mastronarde's newstack command which outputs them as unsigned bytes. Thanks, Jason <type 'exceptions.TypeError'> Exception in Tk callback Function: <bound method MidasUI.processCommand of <Midas.midas_ui.MidasUI instance at 0x276b2d88>> (type: <type 'instancemethod'>) Args: (<Tkinter.Event instance at 0x276a2ea8>,) Event type: KeyPress (type num: 2) Traceback (innermost last): File "/opt/UCSF/Chimera64-1.8.1/lib/python2.7/site-packages/Pmw/Pmw_1_3/lib/PmwBase.py", line 1747, in __call__ return apply(self.func, args) File "/opt/UCSF/Chimera64-1.8.1/share/Midas/midas_ui.py", line 283, in processCommand midas_text.makeCommand(cmdText) File "/opt/UCSF/Chimera64-1.8.1/share/Midas/midas_text.py", line 69, in makeCommand f(c, args) File "/opt/UCSF/Chimera64-1.8.1/share/VolumeFilter/ChimeraExtension.py", line 24, in vop_cmd vop_command(cmdname, args) File "/opt/UCSF/Chimera64-1.8.1/share/VolumeFilter/vopcommand.py", line 91, in vop_command doExtensionFunc(func, fargs, specInfo = spec) File "/opt/UCSF/Chimera64-1.8.1/share/Midas/midas_text.py", line 435, in doExtensionFunc extFunc(*tuple(processedArgs), **kw) File "/opt/UCSF/Chimera64-1.8.1/share/VolumeFilter/vopcommand.py", line 553, in resample_op boundingGrid, False, None, modelId) File "/opt/UCSF/Chimera64-1.8.1/share/VolumeFilter/vopcommand.py", line 168, in add_operation scale = scale[i]) File "/opt/UCSF/Chimera64-1.8.1/share/VolumeViewer/volume.py", line 1535, in add_interpolated_values m[:,:,:] += values <type 'exceptions.TypeError'>: Cannot cast ufunc add output from dtype('float32') to dtype('uint8') with casting rule 'same_kind' ================================================ Event contents: char: delta: 36 height: ?? keycode: 36 keysym: Return keysym_num: 65293 num: ?? send_event: False serial: 39103 state: 16 time: 139539396 type: 2 widget: .59791048.661335648.802317088.802317376.802317448 width: ?? x: 159 x_root: 524 y: 9 y_root: 988 TypeError: Cannot cast ufunc add output from dtype('float32') to dtype('uint8') with casting rule 'same_kind' File "/opt/UCSF/Chimera64-1.8.1/share/VolumeViewer/volume.py", line 1535, in add_interpolated_values m[:,:,:] += values See reply log for Python traceback. -- National Center for Macromolecular Imaging, N420 Baylor College of Medicine, MS BCM125 1 Baylor Plaza Houston, TX, USA

2 3

Is there a way to save animation editing process?
by Li, Leon 03 Apr '14

03 Apr '14

Dear all, I'm try to save all my editing of animation, it seems that 'save session' can't save 'scenes', 'timelines', etc. I wonder is there a way to save all the process, so that I can re-open it to edit. Thanks, Leon

2 2

combining different chains
by faeze mohajer 01 Apr '14

01 Apr '14

good day I have modeled heavy and light chains of myosin now I want to combine them through chimera, would you plz help me on that your faithfully *FAEZE SADAT MOHAJERPhD Student * *Faculty of Bioscience and Medical Engineering (FBME)* *Universiti Teknologi Malaysia (UTM)* Buzzdock ads <http://www.buzzdock.com/setting/V2hpbG9raWk=> Search ResultsFaculty of Biosciences & Medical Engineering<http://fbme.utm.my/>

2 1