- Chimera-users - RBVI Mailing Lists

Save marker sets from command line or python
by Kyle Morris 30 Jan '17

30 Jan '17

Dear Chimera dev and community, I would like to save a marker set from the command line and ultimately do so in a script. Using the write command which I would use for PDB's, marker sets get saved in PDB format, of course. Is there a command to save the marker set as XML .cmm as it would from the volume tracer menu? I notice at the bottom of the following page this is a suggested development: https://www.cgl.ucsf.edu/chimera/docs/ContributedSoftware/volumepathtracer/… - supply Python code for reading and writing the XML files containing saved marker sets Is it there already and I can't find perhaps? Thanks for your help as ever! Best wishes, Kyle

3 2

on chimera unit cell function
by Smith Liu 29 Jan '17

29 Jan '17

Dear All, In Chimera, we do Tools →Higher-Order Structure→Unit Cell, then what will be the diffenrence if we clicked "Pack molecules in unit cell" in comparison with if we did not click "Pack molecules in unit cell"? Does it mean that what we observed (or what Chimera displayed) if we click "Pack molecules in unit cell" , was really part of what it displayed if we did not click "Pack molecules in unit cell" but we increased the "Number of Cells" significantly? I am looking forward to getting a reply from you. Smith

2 3

Metal Geometry Analysis - per frame
by Ilkow, Veronica 27 Jan '17

27 Jan '17

Dear All, Is there a way I could include the metal geometry analysis as a per frame script for a MD movie? Ideally I would like to have an output file which lists the best geometry and essentially prints out the co-ordination table for two atom selections per frame. I am trying to figure out how to use the python or chimera per-frame commands to specifically select an atom and then print an output file. Any hints and suggestions on how to do this will be greatly appreciated. Thank you in advance for your help. -Veronica

2 1

Re: [Chimera-users] dipoles
by Tom Goddard 25 Jan '17

25 Jan '17

Hi Vatsala, To change the arrow thickness change the line of code in the dipole.py script bildString = ".arrow %g %g %g %g %g %g .1 .2 .9" The values at the end control the shape, the “.1” is the cylinder radius, the “.2” is the arrow tip cone radius, and the “.9” is the fraction of the arrow length taken by the cylinder. This is described in the Chimera BILD format documentation here https://www.cgl.ucsf.edu/chimera/docs/UsersGuide/bild.html <https://www.cgl.ucsf.edu/chimera/docs/UsersGuide/bild.html> Tom > On Jan 25, 2017, at 10:02 AM, Sagar, Vatsala (NIH/NEI) [E] wrote: > > Hi Tom, > Do you know how to change the dipole vector thickness. I am using dipole.py program to calculate dipoles in Chimera, and it is giving me a very thin line. > Thanks a lot, I sincerely appreciate it, > Vatsala Sagar > > From: "Sagar, Vatsala (NIH/NEI) [F]" > Date: Friday, January 6, 2017 at 10:09 AM > To: Tom Goddard , "chimera-users(a)cgl.ucsf.edu <mailto:chimera-users@cgl.ucsf.edu> BB" <chimera-users(a)cgl.ucsf.edu <mailto:chimera-users@cgl.ucsf.edu>>, "Sagar, Vatsala (NIH/NEI) [F]" > Subject: Re: [Chimera-users] dipoles > > Hi Tom, > I am attaching the PQR files and raw PDB files of the crystal structure of human gamma D with two monomers rotated 90 degrees to one another. Please take a look at it and tell me what you think, > Thanks, > Vatsala > From: Tom Goddard > Date: Thursday, January 5, 2017 at 7:07 PM > To: "chimera-users(a)cgl.ucsf.edu <mailto:chimera-users@cgl.ucsf.edu> BB" <chimera-users(a)cgl.ucsf.edu <mailto:chimera-users@cgl.ucsf.edu>>, "Sagar, Vatsala (NIH/NEI) [F]" > Subject: Re: [Chimera-users] dipoles > > Hi Vatsala, > > I guess you are refering to the dipole.py script on the Chimera Python scripts web page: > > http://plato.cgl.ucsf.edu/trac/chimera/wiki/Scripts <http://plato.cgl.ucsf.edu/trac/chimera/wiki/Scripts> > > Yes it should give the same answer for a molecule and a copy of that molecule rotated by 90 degrees. Without have the structure and rotated PQR files I don’t know what is going wrong. > > Tom > > >> On Jan 4, 2017, at 11:17 AM, Sagar, Vatsala (NIH/NEI) [F] wrote: >> >> Hi, >> Do you know if there is a method to calculate dipoles in Chimera? I have generated a PQR structure in predicted protonation state and uploaded it onto Chimera. I then used a downloaded app called dimer.py to calculate the dipole moments of two of the same monomers rotated 90degrees. Technically they should both give the same dipole as the structures align 100% but according to the PDB coordinates are only rotated relative to one another. The dipole.py scipt is calculating them as two different dipoles. Can you help me with this? >> Thanks, >> Vatsala Sagar >> _______________________________________________ >> Chimera-users mailing list: Chimera-users(a)cgl.ucsf.edu <mailto:Chimera-users@cgl.ucsf.edu> >> Manage subscription: http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users <http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users> >

1 0

on Multiscale Models
by Smith Liu 23 Jan '17

23 Jan '17

Dear Al, In the "Mutiscale Models" Windows , will you please tell me what does it mean for the item "Range"? I am looking forward to getting a reply from you. Smith

2 1

Python script for reading volume values at atom positions
by Kyle Morris 19 Jan '17

19 Jan '17

Dear chimera dev and users, Can anyone confirm with the latest chimera version (i’m running v1.11.2 build 41376) if this python script for measuring map values at atoms works for them? The script was reported previously on the users list here: http://plato.cgl.ucsf.edu/pipermail/chimera-users/2007-May/001587.html For me I get a python import error: from VolumeViewer import Data_Region, full_region, Rendering_Options Thanks for your insight! Best wishes, Kyle

3 6

beta sheet to loop
by chemocev marker 19 Jan '17

19 Jan '17

Hi I am looking at some region of the protein in chimera and it looks like a loop but the from the procheck analysis it shows a beta-hair pin like structure. In the electron density, its shows a hydrogen bonding between the main chain atom. Why is this difference and how to overcome on it. best Jiri

2 1

forget last question
by A G Szabo 18 Jan '17

18 Jan '17

Forget the last question regarding visualizing nucleotides as balls and sticks. Figured it out. ags

1 0

quick question
by A G Szabo 18 Jan '17

18 Jan '17

I am looking at a protein polynucleotide structure. How do I turn the rectangles that represent the individual nucleotide into balls and sticks? Thanks Arthur G. Szabo

1 0

Autodock vina + modeller + ATOM syntax error
by Zirpel, Bastian 18 Jan '17

18 Jan '17

Hi everyone, I'm using Chimera 1.11.2 on windows. I am experience following error after my workflow. I fetched a known pdb model from the servers and generated several structural models of similar sequences using the integrated MODELLER script. Since my pdb model contains a covalently attached cofactor in the active site, I let MODELLER include non-water HETATM into my structural models as well. Afterwards I tried to do substrate/product binding studies with the known pdb model as well as with the generated structural models using the AutoDockVina script. However, when I run the script with the known pdb model, everything works fine. But when I'm running the script with the structural models, I receive following error msg: "Parse error on line 8553 in file "receptor.pdbqt": ATOM syntax incorrect: "Ho" is not a valid AutoDock type. Note that AutoDock atom types are case-sensitive." It should be noted that this error does not occur if I won't include non-water HETATM by MODELLER. Makes sense, since having a look at the generated pdbqt file, the FAD lines are generating the error. " HETATM 8554 HO3' FAD 525 62.339 60.023 -29.788 1.00 0.00 0.238 Ho HETATM 8555 HO2' FAD 525 66.741 58.106 -28.406 1.00 0.00 0.397 Ho HETATM 8556 HO3B FAD 525 60.634 66.984 -32.274 1.00 0.00 0.221 Ho HETATM 8557 HO4' FAD 525 65.182 62.428 -29.602 1.00 0.00 0.217 Ho HETATM 8558 1H6A FAD 525 68.751 67.392 -38.565 1.00 0.00 0.157 HD HETATM 8559 2H6A FAD 525 68.156 65.783 -38.183 1.00 0.00 0.157 HD HETATM 8560 HO2B FAD 525 62.183 68.217 -34.408 1.00 0.00 0.253 Ho " I found that AutoDockVina needs a "HD" instead of a "Ho" entry at these positions? Is this correct? So my question is, why is AutoDock vina generating a pdbqt file in the first place which it cannot read properly in the next step? How can I prevent this? I cannot just change the pdbqt entries as AutoDockVina will overwrite it again. Thanks in advance for any help! Bastian

2 1