July 2008 - Chimera-users - RBVI Mailing Lists

Distance between atoms using Python
by Tom Goddard 18 Jul '08

18 Jul '08

Hi Ed, Attached is some Python code for controlling ligand display based on distance between two atoms. Tom Edward Brignole wrote: > > Tom, > > Thanks again for your help with that atom display. Got another > MDtrajectory issue for you. I've got a ligand in each of 10 domains > in the structure and I only want to display the ligands in the frames > that have less than a certain distance between a particular atom in > the ligand and a particular atom in a neighboring domain. > > The distance command won't work for me because all it does is add a > distance to the distance dialog when I need it as a variable in the > script to evaluate with "if distance < 10: display atomspecs" > > Sorry to keep bugging you with these things but I assure you it isn't > without trying for several hours to no avail before finally resorting > to asking. > > Ed > def atom_distance(spec1, spec2): from chimera.specifier import evalSpec atom1 = evalSpec(spec1).atoms()[0] atom2 = evalSpec(spec2).atoms()[0] d = atom1.xformCoord().distance(atom2.xformCoord()) return d from chimera import runCommand if atom_distance('#0:15.B@C4', '#0:7.A@N3') >= 10: runCommand('display #0:SPM') else: runCommand('~display #0:SPM')

1 0

center of mass for multiple models.
by Jean Didier Pie Marechal 18 Jul '08

18 Jul '08

Hi everyone, I am stuck with a script and my mistake seems to come from my poor understanding on how to treat atom coordinates with multiple models. I join to this e-mail a chimera file with a simple system to illustrate my problem. In the file, I have two indentique phenyl moleculer each of them in different models and I would like to stock the center of mass of both molecules in a list. I wrote this very simple script using the way of calcuting the center of mass you provided me recently. import chimera import BuildStructure from chimera import runCommand as run from chimera import selection, Vector from chimera import openModels, Molecule, Point #Stock the two fragments in the list FRAG run("select") frag=[] for i in chimera.selection.currentResidues(): frag.append(i) #Calculate center of mass of both entities #stocked in the list MAS mas=[] for i in range(len(frag)): sum=0 center = Vector(0, 0, 0) for a in frag[i].atoms: sum += a.element.mass center += a.element.mass * Vector(*a.coord().data()) mas.append(Point(center[0]/sum,center[1]/sum,center[2]/sum)) print i,mas[i] The results I have are 0 20.0465 70.3735 23.9063 1 20.0465 70.3735 23.9064 Basically, the coordinates of the center of mass are the same? What did I do wrong in the way I treated this problem? Thanks a lot for any answer, All the best JD Dr. Jean-Didier Maréchal Professor Lector Unitat de Química Física Departament de Química Universitat Autònoma de Barcelona Edifici C.n. 08193 Cerdanyola (Barcelona) Tel: +34.935814936 e-mail: JeanDidier.Marechal(a)uab.es ----- Missatge original ----- De: chimera-users-request(a)cgl.ucsf.edu Data: Dijous, Juliol 17, 2008 9:00 pm Assumpte: Chimera-users Digest, Vol 63, Issue 17 > Send Chimera-users mailing list submissions to > chimera-users(a)cgl.ucsf.edu > > To subscribe or unsubscribe via the World Wide Web, visit > http://www.cgl.ucsf.edu/mailman/listinfo/chimera-users > or, via email, send a message with subject or body 'help' to > chimera-users-request(a)cgl.ucsf.edu > > You can reach the person managing the list at > chimera-users-owner(a)cgl.ucsf.edu > > When replying, please edit your Subject line so it is more specific > than "Re: Contents of Chimera-users digest..." > > > Today's Topics: > > 1. Re: Possible image saving bug (Greg Couch) > 2. Possible image saving bug (Elaine Meng) > > > -------------------------------------------------------------------- > -- > > Message: 1 > Date: Thu, 17 Jul 2008 10:35:26 -0700 (PDT) > From: Greg Couch <gregc(a)cgl.ucsf.edu> > Subject: Re: [Chimera-users] Possible image saving bug > To: Walter Novak <wnovak(a)brandeis.edu> > Cc: Chimera forum <chimera-users(a)cgl.ucsf.edu> > Message-ID: <Pine.OSF.4.63.0807171001010.1939086(a)guanine.cgl.ucsf.edu> > Content-Type: TEXT/PLAIN; charset=US-ASCII; format=flowed > > On Thu, 17 Jul 2008, Walter Novak wrote: > > > Hi All, > > > > I submitted a bug report on this, but the behavior is both odd > and > > reproducible for me. Using the latest production release and > trying to save > > an image using supersampling, the image is not stitched back > together > > appropriately. It is odd, because it always happens for a certain > session > > file I have. It doesn't seem to happen if I open a new molecule > and make a > > picture. I can provide the session file if needed. First image is > using 3x3 > > and second is 6x6. > > I see from the bug report that you're using Mac OS X 10.5. We are > having > "issues" with the new X11 application in 10.5. And we have seen > that > slightly different X11 applications are installed on otherwise > identical > Macs and that reinstalling from the Mac OS X Install Disk fixed a > non-chimera bug. > > Also, I have read on the web that if you install the 10.4 X11 > application > on 10.5, then OpenGL applications will work better. Since that is > not an > option for everyone, we are spending more time to improve the Aqua > version > that is available as a Snapshot release on the chimera download > page. It > has different bugs, but the graphics is better. So you should try > it out > and see how well it works for you. > > Greg Couch > UCSF Computer Graphics Lab > > > ------------------------------ > > Message: 2 > Date: Thu, 17 Jul 2008 10:40:13 -0700 > From: Elaine Meng <meng(a)cgl.ucsf.edu> > Subject: [Chimera-users] Possible image saving bug > To: wnovak(a)brandeis.edu > Cc: Chimera BB <chimera-users(a)cgl.ucsf.edu> > Message-ID: <B15630DD-299C-4B48-A8B4-469FA8995872(a)cgl.ucsf.edu> > Content-Type: text/plain; charset=US-ASCII; format=flowed; delsp=yes > > Hi Wally, > I also had this tiling problem. We think it may be something about > > X11 in Mac OS 10.5 (possibly just 10.5.4), as the problem has not > been > observed with 10.4/X11. > > Two workarounds that were successful in my case (I don't know if > they > will fix all instances of the problem): > > (A) simply show the Command Line before saving the image. I tried > hiding it again after that, but then the problem reoccurred, albeit > in > a more subtle fashion. > > (B) save the session, restore it in the aqua version of Chimera > (still > on Mac OS 10.5.4 but not using X11) and save the image. Yeah, (A) > is > easier, but I didn't discover that one until this morning, after > receiving your mail. > > I hope these solutions are general! Sorry about the inconvenience, > 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 > > > > > ------------------------------ > > _______________________________________________ > Chimera-users mailing list > Chimera-users(a)cgl.ucsf.edu > http://www.cgl.ucsf.edu/mailman/listinfo/chimera-users > > > End of Chimera-users Digest, Vol 63, Issue 17 > ********************************************* >

2 1

center of mass for multiple models
by Jean Didier Pie Marechal 18 Jul '08

18 Jul '08

Hi everyone, I am stuck with a script and my mistake seems to come from my poor understanding on how to treat atom coordinates with multiple models. I join to this e-mail a chimera file with a simple system to illustrate my problem. In the file, I have two indentique phenyl moleculer each of them in different models and I would like to stock the center of mass of both molecules in a list. I wrote this very simple script using the way of calcuting the center of mass you provided me recently. import chimera import BuildStructure from chimera import runCommand as run from chimera import selection, Vector from chimera import openModels, Molecule, Point #Stock the two fragments in the list FRAG run("select") frag=[] for i in chimera.selection.currentResidues(): frag.append(i) #Calculate center of mass of both entities #stocked in the list MAS mas=[] for i in range(len(frag)): sum=0 center = Vector(0, 0, 0) for a in frag[i].atoms: sum += a.element.mass center += a.element.mass * Vector(*a.coord().data()) mas.append(Point(center[0]/sum,center[1]/sum,center[2]/sum)) print i,mas[i] The results I have are 0 20.0465 70.3735 23.9063 1 20.0465 70.3735 23.9064 Basically, the coordinates of the center of mass are the same? What did I do wrong in the way I treated this problem? Thanks a lot for any answer, All the best JD Dr. Jean-Didier Maréchal Professor Lector Unitat de Química Física Departament de Química Universitat Autònoma de Barcelona Edifici C.n. 08193 Cerdanyola (Barcelona) Tel: +34.935814936 e-mail: JeanDidier.Marechal(a)uab.es ----- Missatge original ----- De: chimera-users-request(a)cgl.ucsf.edu Data: Dijous, Juliol 17, 2008 9:00 pm Assumpte: Chimera-users Digest, Vol 63, Issue 17 > Send Chimera-users mailing list submissions to > chimera-users(a)cgl.ucsf.edu > > To subscribe or unsubscribe via the World Wide Web, visit > http://www.cgl.ucsf.edu/mailman/listinfo/chimera-users > or, via email, send a message with subject or body 'help' to > chimera-users-request(a)cgl.ucsf.edu > > You can reach the person managing the list at > chimera-users-owner(a)cgl.ucsf.edu > > When replying, please edit your Subject line so it is more specific > than "Re: Contents of Chimera-users digest..." > > > Today's Topics: > > 1. Re: Possible image saving bug (Greg Couch) > 2. Possible image saving bug (Elaine Meng) > > > -------------------------------------------------------------------- > -- > > Message: 1 > Date: Thu, 17 Jul 2008 10:35:26 -0700 (PDT) > From: Greg Couch <gregc(a)cgl.ucsf.edu> > Subject: Re: [Chimera-users] Possible image saving bug > To: Walter Novak <wnovak(a)brandeis.edu> > Cc: Chimera forum <chimera-users(a)cgl.ucsf.edu> > Message-ID: <Pine.OSF.4.63.0807171001010.1939086(a)guanine.cgl.ucsf.edu> > Content-Type: TEXT/PLAIN; charset=US-ASCII; format=flowed > > On Thu, 17 Jul 2008, Walter Novak wrote: > > > Hi All, > > > > I submitted a bug report on this, but the behavior is both odd > and > > reproducible for me. Using the latest production release and > trying to save > > an image using supersampling, the image is not stitched back > together > > appropriately. It is odd, because it always happens for a certain > session > > file I have. It doesn't seem to happen if I open a new molecule > and make a > > picture. I can provide the session file if needed. First image is > using 3x3 > > and second is 6x6. > > I see from the bug report that you're using Mac OS X 10.5. We are > having > "issues" with the new X11 application in 10.5. And we have seen > that > slightly different X11 applications are installed on otherwise > identical > Macs and that reinstalling from the Mac OS X Install Disk fixed a > non-chimera bug. > > Also, I have read on the web that if you install the 10.4 X11 > application > on 10.5, then OpenGL applications will work better. Since that is > not an > option for everyone, we are spending more time to improve the Aqua > version > that is available as a Snapshot release on the chimera download > page. It > has different bugs, but the graphics is better. So you should try > it out > and see how well it works for you. > > Greg Couch > UCSF Computer Graphics Lab > > > ------------------------------ > > Message: 2 > Date: Thu, 17 Jul 2008 10:40:13 -0700 > From: Elaine Meng <meng(a)cgl.ucsf.edu> > Subject: [Chimera-users] Possible image saving bug > To: wnovak(a)brandeis.edu > Cc: Chimera BB <chimera-users(a)cgl.ucsf.edu> > Message-ID: <B15630DD-299C-4B48-A8B4-469FA8995872(a)cgl.ucsf.edu> > Content-Type: text/plain; charset=US-ASCII; format=flowed; delsp=yes > > Hi Wally, > I also had this tiling problem. We think it may be something about > > X11 in Mac OS 10.5 (possibly just 10.5.4), as the problem has not > been > observed with 10.4/X11. > > Two workarounds that were successful in my case (I don't know if > they > will fix all instances of the problem): > > (A) simply show the Command Line before saving the image. I tried > hiding it again after that, but then the problem reoccurred, albeit > in > a more subtle fashion. > > (B) save the session, restore it in the aqua version of Chimera > (still > on Mac OS 10.5.4 but not using X11) and save the image. Yeah, (A) > is > easier, but I didn't discover that one until this morning, after > receiving your mail. > > I hope these solutions are general! Sorry about the inconvenience, > 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 > > > > > ------------------------------ > > _______________________________________________ > Chimera-users mailing list > Chimera-users(a)cgl.ucsf.edu > http://www.cgl.ucsf.edu/mailman/listinfo/chimera-users > > > End of Chimera-users Digest, Vol 63, Issue 17 > ********************************************* >

1 0

Possible image saving bug
by Elaine Meng 17 Jul '08

17 Jul '08

Hi Wally, I also had this tiling problem. We think it may be something about X11 in Mac OS 10.5 (possibly just 10.5.4), as the problem has not been observed with 10.4/X11. Two workarounds that were successful in my case (I don't know if they will fix all instances of the problem): (A) simply show the Command Line before saving the image. I tried hiding it again after that, but then the problem reoccurred, albeit in a more subtle fashion. (B) save the session, restore it in the aqua version of Chimera (still on Mac OS 10.5.4 but not using X11) and save the image. Yeah, (A) is easier, but I didn't discover that one until this morning, after receiving your mail. I hope these solutions are general! Sorry about the inconvenience, 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

1 0

Possible image saving bug
by Walter Novak 17 Jul '08

17 Jul '08

Hi All, I submitted a bug report on this, but the behavior is both odd and reproducible for me. Using the latest production release and trying to save an image using supersampling, the image is not stitched back together appropriately. It is odd, because it always happens for a certain session file I have. It doesn't seem to happen if I open a new molecule and make a picture. I can provide the session file if needed. First image is using 3x3 and second is 6x6. Also, when using pov-ray, the balls are completely devoid of shininess, regardless of my shininess settings (I set it to 128). That is the last image. Thanks, Wally Walter R.P. Novak, Ph.D. Postdoctoral Fellow Rosenstiel Basic Medical Research Center Brandeis University 415 South St. MS 029 Waltham, MA 02454-9110 Phone: (781) 736-4944 Fax: (781) 736-2405

2 1

Coloring a binding site
by Tom Goddard 17 Jul '08

17 Jul '08

Hi Andy, To color a binding site I would select the ligand (DNA in your case) then select all the nearby atoms using Select / Zone..., then use Actions / Color to color those atoms. Tom Anindito Sen wrote: > Another question - I need to color a portion of the crystal structure > differently (ribbon diagram) in order to indicate the binding sites of > the DNA with the protein molecule and how do I do that. > > Thanks and Best Wishes > Andy

2 2

possible bug
by Ben Keshet 16 Jul '08

16 Jul '08

Hi, when I un-select atoms, I had some weird long lines crossing my structures across (picture attached). I restarted the program how ever the same exact thing happened again. The two files that were opened worked perfectly in the past and before I un-selected the atom (ctrl+left click in the "air"). Other wise - great software. Thank you. Ben

2 1

Re: [Chimera-users] History feature for command line between sessions
by Eric Pettersen 15 Jul '08

15 Jul '08

> On Apr 17, 2008, at 8:52 AM, Boaz Shaanan wrote: > > > The history of commands in the command line mode is retained > > within a session and is very convenient but I was wondering whether > > there is a feature that will allow you to recall commands from > > previous sessions (such as in O for example - it generates an .o- > > history file which is being constantly updated and is very > > convenient). If not, is it possible to request the addition of such > > a feature ? There are many ways to implement this of course and I'm > > sure you'll find the best way, if it's not already there lurking in > > Chimera. > > Hi Boaz, > Good idea -- and shouldn't be too hard to implement. I was thinking > of remembering the last 100 commands. Will that be enough you think? > > --Eric This is implemented now. It should show up in the next daily build that works. Only remembers 60 by default, but that number is changeable in Preferences. --Eric Eric Pettersen UCSF Computer Graphics Lab http://www.cgl.ucsf.edu

1 0

Can't type in Chimera
by Thomas Goddard 15 Jul '08

15 Jul '08

Hi Andy, I don't know how Chimera could get in a state where it will not accept any typed input. Do you recall the last thing you did in Chimera before it stopped accepting key input? If I were stuck in this situation and could not type a file name to save the Chimera session I would create an empty file with the session file name I wanted using a text editor. Then choose that file with the mouse in the Chimera session saving dialog. Desperate measures. If you can reproduce the problem please tell me how and I will fix Chimera. Tom Anindito Sen wrote: > HI Tom > > I have installed the latest version of Chimera on my LINUX box Fedora > Core 6.0. I faced a strange problem- while carrying out a docking of a > crystal structure in a density map I was using the command line for the > differnt movements and controls. After some time I am unable to type > anything on the command line. Infact I am even unable to save the > session as I cannot type the "NAME" of the session. I checked my keyboad > and typed on the terminal and on text editiors and it worked normally. > > Waiting for your response. > > Best Wishes > Andy > > > Dr. Anindito Sen (Ph.D) > Research Associate , Dept. of Biochemistry and Molecular Genetics > University of Virginia > Box 800733 Charlottesville, VA 22908

1 0

Re: [Chimera-users] reg. ribbon reprsentation
by Elaine Meng 15 Jul '08

15 Jul '08

>> >> Greetings. I have already posted a similar question on the forum >> previously and hence i am writing to your email. My problem is when >> i open a protein pdb, i get wonderful ribbon representation when i >> opt for it. But when i give a protein pdb with only a CA atoms/ >> CA,C,N atoms, i dnt get the ribbon representation. I understand >> that ribbon is constructed by considering CA,C,N,O and OXT atoms. >> What i did was that i simply removed O and OXT atoms by going >> through Tools--> Depiction --> Ribbon style Editor. It may be a >> wrong idea to do that. But i would appreciate if you could write me >> some way to draw the ribbon for protein having only CA or CA,C,N >> atoms IF POSSIBLE. >> >> Regards, >> Bala >> Hi Bala, If you have N-CA-C atoms, it actually works pretty well if you open Ribbon Style Editor, go to the "Residue Class" tab, choose "amino acid" but then change the "Orientation atom" to C instead of O, and Apply. (You still need to show ribbon if you didn't already do that.) You can save and name the new residue class for re-use in some later Chimera session if you want. The resulting ribbon looks the same to me as the standard ribbon made from N,CA,C,O (I tested this by writing a PDB of only the N,CA,C atoms from a complete structure, reading it back in, and experimenting with Ribbon Style Editor). http://www.cgl.ucsf.edu/chimera/docs/ContributedSoftware/ribbonstyle/ribbon… If you have CA atoms only, then there is no way to figure out the plane of the ribbon - it does not work to set both the "Guide atom" and "Orientation atom" to CA. Unfortunately this prevents you from displaying "licorice" ribbons even though they don't use an orientation plane. My only solution in this case is to use some program to build the additional atoms from the CA trace, for example this web server: http://www.ebi.ac.uk/Tools/maxsprout/index.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

1 0