May 2011 - Chimera-users - RBVI Mailing Lists

Re: [Chimera-users] Pseudobonds disappearing during morphing movie
by Elaine Meng 17 May '11

17 May '11

Hi Alan, Sorry, there isn't currently a command to apply residue class, only scaling (ribscale) and style (ribrepr): <http://www.cgl.ucsf.edu/chimera/docs/ContributedSoftware/ribbonstyle/ribbon…> Some possibilities: - save session with morph already calculated and residue class applied. Your script would then restore the session, play the trajectory, and record the movie. I agree this is not as nice, since then you are saving two files instead of one, and the session file is large. You wouldn't have to wait for the morph to recalculate, however. - use a small python script to apply the residue class. Your Chimera script could open that python file, which automatically executes it. One of the other developers would have to supply the python (beyond my skill set, sorry), and that also requires having two files, albeit much smaller than a session file. In general when playing a morph with ribbons that show secondary structure, consider whether you need to recalculate the secondary structure assignments with "ksdssp" at each frame. To execute a set of commands at each frame, one would use the "perframe" command in combination with "alias": <http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/perframe.html> <http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/alias.html> <http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/ksdssp.html> For example, your script could contain something like: alias ^redefine ksdssp #2 perframe redefine movie record coordset #2 1,21; wait 21; ~perframe I hadn't used ksdssp/perframe in my examples since (perhaps fortuitously) the secondary structure didn't change that much. If you are just displaying a tube or licorice which doesn't show secondary structure with different widths, it doesn't matter, however. You're right that the shape tube surface will not automatically morph along with the structures, but it might work to destroy/recreate it at each coordinate frame using "perframe". Interesting idea... I haven't tried it, though! I hope this helps, Elaine On May 13, 2011, at 6:10 AM, Alan Cheung wrote: > Thanks Elaine! Yes, the reported thickness was 3.333 (pixels?) and that seems to match the width of the pseudobond i see on the video. The width isn't too bad, and i can live with it. > > Thanks for the explanations and tip too. > > One last question : can I apply a saved residue class (from the ribbon style editor) via a command? I'm morphing a mixture of ribbon and side chain atoms and (for a leucine residue) the CA-CB bond is unrealistically long, but seems ok if i remove the orientation atom (O) from the residue class. I saw the "shape tube" feature in the documentation that draws a tube through CA atoms, but i guessed it's not morphable (and might also not include my precious pseudobond ;) ) > > Many thanks again for the quick replies! > > Alan > > On 12/05/2011 19:22, Elaine Meng wrote: >> Hi Alan, >> Never mind about the PDBs, I made up my own by deleting parts of the structures in my favorite morph test case. >> >> I think I answered all of your issues except for the pseudobonds completely disappearing with "morph movie nogui true" -- I did reproduce that problem, and it seems like a bug, albeit mysterious. I'll file a report about that one, but in the meanwhile I hope you will be able to manage without using that option. >> >> One minor thing is that I generally recommend adding a few more wait frames to the recording after the morph plays, to avoid any "motion blur" effect at the end you might get as a result of encoding. I didn't see much of that effect in my test case proteins using your script, however. It may be because the morph is spread out over 50 frames so there is less differences between adjacent frames, but I was also using different PDBs than you. >> Best, >> Elaine >> >> On May 12, 2011, at 9:25 AM, Elaine Meng wrote: >> >>> Hi Alan, >>> I haven't tried the whole script yet (can you send me the PDB files? I'll keep them private, i.e. only for the development group here if you send them just to me) but will try to address the low-hanging fruit in your set of questions: >>> >>> (A) lines thinner in saved movie. In general, lines such as silhouette edges and pseudobonds may turn out thinner in saved images and movies when you are doing "supersampling," which initially saves a larger image and samples it down to the final requested size. You can check if this is the issue by using File... Save Image to save a single image -- the dialog should report the maximum achievable linewidth given the requested size and amount of supersampling. The reason is that there is an upper limit to how fat your computer can draw the lines (at least via the method used by Chimera). If you have 2x2 supersampling, for example, Chimera will try to make the lines twice as fat in the initial larger image, but may run into this upper limit. Then when the image is sampled down to the final size, the lines may not be as thick as they appear in the interactive display. You could try using lower supersampling. Also make sure you are using a recent version (1.5.3 candida > te ! >>> or 1.6) because there was a lines-too-thin-in-images bug a couple of months ago. >>> >>> (B) setattr "p" vs. "g" -- this is because some attributes are for individual pseudobonds (p) and some attributes are defined for whole pseudobond groups (g). Individual pseudobonds don't have a linewidth attribute, it is set at the group level. All the pseudobonds in that group (for example, hydrogen bonds) would get the same linewidth and it is not possible to make different ones in the same group have different widths. If you are looking at the Selection Inspector or some other attributes inspector, it will show you which attributes are for individual pseudobonds and which are for groups. The balloon help will give the attribute name and allowed values. >>> <http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/modelattrib.html> >>> >>> I will try to look at your other issues and get back to you, >>> 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 May 12, 2011, at 8:54 AM, Alan Cheung wrote: >>> >>>> Hi there - i'm scripting a complex movie and am having problems with >>>> pseudobonds during morphing. Script is pasted at the end. >>>> >>>> Part of the movie contains a morph between two structures, both of >>>> which are missing the same amino acid segment. Loading the two >>>> structures show the pseudobond for the missing residues and the >>>> resulting morph also contains the same pseudobond. >>>> >>>> I've been changing the thickness of the pseudobond with : >>>> >>>> setattr g lineWidth 8 >>>> >>>> and while the two static structures and the morph in the chimera >>>> window have a thicker pseudobond, the resulting movie (encoded to >>>> .mov) does not. Am i doing something wrong? >>>> >>>> Alan >>>> >>>> p.s. incidentally, when morphing by script, if one includes "nogui >>>> true" with the "morph movie" command, pseudobonds aren't kept in the >>>> morph at all. A bug? >>>> p.p.s. can someone explain to me the difference between "setattr g" >>>> and "setattr p"? I've read the manual but i don't understand why only >>>> the former seems to work for me. >>>> >>>> ## movie script >>>> >>>> open model1.pdb >>>> open model2.pdb >>>> morph start #0 frames 50 name morph1 method linear >>>> morph interpolate #1 name morph1 >>>> morph movie name morph1 >>>> coordset #2 1,1 >>>> setattr g lineWidth 8 >>>> ~modeldisplay # >>>> >>>> movie record supersample 3 directory ~/TEMP/ >>>> modeldisplay #2 >>>> coordset #2 1,51 ; wait 51 >>>> movie stop >>>> movie encode output morph1.mov >>> >>> >>> _______________________________________________ >>> Chimera-users mailing list >>> Chimera-users(a)cgl.ucsf.edu >>> http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users >> >> > > -- > Alan Cheung > Gene Center > Ludwig-Maximilians-University > Feodor-Lynen-Str. 25 > 81377 Munich > Germany > Phone: +49-89-2180-76845 > Fax: +49-89-2180-76999 > E-mail: cheung(a)lmb.uni-muenchen.de

2 2

labeling residues
by Daniel Gurnon 16 May '11

16 May '11

Hi all, A student of mine wants to make a wire sculpture of a protein, and I'm thinking of ways Chimera could make the task easier. Ideally, I would create a script that would display any protein in wire-form, with residue positions clearly and unobtrusively labeled. I'd like to use the command line to display a backbone as a licorice ribbon with residue numbers and 1-letter codes at the alpha carbon of each residue. so starting out, I would use: ribscale licorice; ribbackbone; rlabel @ca; But instead of "ALA 24.B", I want "A 24". Though I can get what I want using the Actions menu, I can't find a way to create custom rlabels via the command line. Is there a way to do this? If possible, I would also like to position each label near the spot on the ribbon closest to the alpha carbon (the "primary atom"). I can see how to do this with the molecule attributes panel, but I'm having difficulty setting the placement through the command line. And while I'm at it...is there a way to color a spot on the ribbon closest to the alpha carbon a different color than the rest of the ribbon? Thanks for the help! Dan -- ____________________________ Daniel Gurnon, Ph. D. Associate Professor of Chemistry DePauw University Greencastle, IN 46135 p: 765-658-6279 e: danielgurnon(a)depauw.edu

3 5

Can eraser position be defined by coordinates?
by Kelly Dryden 16 May '11

16 May '11

Hi, I would like to place the origin of the eraser to correspond with the center of my map but haven't found a way to define its position based on coordinates, only by moving it with the mouse. Is there a line command or other menu option that allows me to set the position more precisely? Thanks, Kelly

2 1

rigid body rotation of subunit
by Dhiraj Srivastava 13 May '11

13 May '11

Hi I have the structure of two homologous protein. These proteins are dimeric. when I am superimposing chain A of one protein on to the chain A of other, it looks like chain B of one protein is related to chain B of other protein by rigid body rotation. How can I calculate the angle of rotation in chimera? Is there any way, I can draw a vector along beta sheet? Thank you. Dhiraj

3 2

shape factor
by maznar 13 May '11

13 May '11

Hello, We are using Chimera to classify viruses. We use pdb files from viper. We want a way to obtain exactly shape factor and/or compare with an icosahedron. We have problem do it. We open pdb file, draw icosahedron using "hkcage" (all models actives) and try to compare that two figures using "fit in map", or similar but doesn't work. What is the best and exactly way to do it? Thank you.

2 3

square of density volume
by 仁田亮 12 May '11

12 May '11

Hi chimera users, I want to calculate the square value for the volume of the electron density because I want to calculate the variance map. Do anyone know whether we can calculate the square value of the volume by chimera? And if we can, could you teach me how to calculate? Thanks for your time and help, Ryo

5 12

Pseudobonds disappearing during morphing movie
by Alan Cheung 12 May '11

12 May '11

Hi there - i'm scripting a complex movie and am having problems with pseudobonds during morphing. Script is pasted at the end. Part of the movie contains a morph between two structures, both of which are missing the same amino acid segment. Loading the two structures show the pseudobond for the missing residues and the resulting morph also contains the same pseudobond. I've been changing the thickness of the pseudobond with : setattr g lineWidth 8 and while the two static structures and the morph in the chimera window have a thicker pseudobond, the resulting movie (encoded to .mov) does not. Am i doing something wrong? Alan p.s. incidentally, when morphing by script, if one includes "nogui true" with the "morph movie" command, pseudobonds aren't kept in the morph at all. A bug? p.p.s. can someone explain to me the difference between "setattr g" and "setattr p"? I've read the manual but i don't understand why only the former seems to work for me. ## movie script open model1.pdb open model2.pdb morph start #0 frames 50 name morph1 method linear morph interpolate #1 name morph1 morph movie name morph1 coordset #2 1,1 setattr g lineWidth 8 ~modeldisplay # movie record supersample 3 directory ~/TEMP/ modeldisplay #2 coordset #2 1,51 ; wait 51 movie stop movie encode output morph1.mov

2 2

How to show an icosahedral cage
by Tom Goddard 11 May '11

11 May '11

Hi Yen-Chywan, The Chimera "hkcage" command makes a cage of pentagons and hexagons to show virus symmetry. For example, open emdbID:1058 hkcage 1 1 radius 170 orient 222r linewidth 4 color blue The hkcage command has options for the lattice parameters for your virus, its radius, what icosahedral coordinate system is being used. Here's the documentation. http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/hkcage.html http://www.cgl.ucsf.edu/chimera/tutorials/volumetour/volumetour.html#icosca… The meshmol command can make the edges of the cage cylinders instead of lines. http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/meshmol.html The HHMI virus viewer Chimera tool is not part of the standard Chimera distribution. It was made for a CD distributed by the Howard Hughes Medical Institute http://wwwtest.cgl.ucsf.edu/Outreach/ViralOutbreak/index.html http://www.hhmi.org/catalog/main?action=product&itemId=366 We don't distribute it separately. Tom > Dear Tom, > > I attend Scripps cryoEM workshop before and meet you there. > > Attach please find a screen capture of a virus. But there is a > icosahedra cage. My question is how to display that cage with my virus > volume? > > I did find it in HHMI virus viewer demo movie. But where can download > this viewer? > > Thanks > > Yen-Chywan >

1 1

constructing a dimer
by uzma saqib 11 May '11

11 May '11

Hello, I want to totally freeze the movement of one protein (which is a homology model) and open another molecule of the same protein, which I want to keep movable, so that I can make a dimer out of this at my chosen place on the first protein molecule. Plz help . thanks uzma

2 1

System Compatiblility
by Satterlee, James D 11 May '11

11 May '11

Hello, I've not yet downloaded or used Chimera, but I am contemplating it. I am looking for Windows System compatibility information for the various newer versions of Chimera. To me it wasn't obvious on the Chimera home site. I do not use Window's 7 and I'm wondering which version might be most trouble-free for a Windows XP version with all service packs installed. Thanks, Jim Satterlee James D. Satterlee Professor of Chemistry Washington State University Pullman, WA 99164-4630 hemeteam(a)wsu.edu 509.335.8620

2 1