June 2012 - Chimera-users - RBVI Mailing Lists

Re: [Chimera-users] Ellipsoids for 800 pdbs
by Elaine Meng 14 Jun '12

14 Jun '12

On Jun 14, 2012, at 4:54 AM, mselvaraj(a)mbu.iisc.ernet.in wrote: > Dear Meng, > I need to calculate ellipsoid axis for some 800 pdbs that I have in a > directory using measur inertia command. I was trying to load one by one and > calcualte, but I was huge. I read about python script for running in > batches. Is there any chimera script for doing this?. I would like to get > your help in doing this. Thanks in advance. > > Yours sincerely, > M. Selvaraj, > Graduate student, > Molecular Biophysics Unit, > Indian Institute of Science, > Bangalore 560 012. Dear M. Selvaraj, The Chimera Programmer's Guide includes a primer on how to write a python script to process multiple structures: <http://www.cgl.ucsf.edu/chimera/docs/ProgrammersGuide/basicPrimer.html> It is better to send Chimera questions to chimera-users(a)cgl.ucsf.edu (or if advanced programming issues, to chimera-dev(a)cgl.ucsf.edu) instead of to me directly. 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

1 0

addH: Weird results
by George Tzotzos 11 Jun '12

11 Jun '12

Hi everybody, For several months I have been saving docking results from Autodock, adding hydrogens and submitting them to Amber for molecular dynamics. Chimera worked perfectly well. A typical result is attached as <file1.pdb>. I repeated the process with a different pose of the same ligand and got very weird results. In fact the atoms were numbered differently and some atoms were not numbered at all, or rather they were given ***. I submitted some of the original ligands, that is the ligands to which Chimera added Hs correctly. A typical example is attached <docking.result_8_2.pdb>. The result of <file1.pdb> was not reproduced. The errors described above were repeated (see <file2.pdb>). Trying to renumber the atoms according to <file1.pdb> did not help. I am a bit baffled by this. Any help will be most welcome. Best regards George PS. The submitted ligand is bisabolol.

2 3

TER characters between residues in exported PDB ?
by Marek Maly 07 Jun '12

07 Jun '12

Hello all, I would like to know if there is any way how to export from Chimera molecular structure in PDB format but with TER string between all residues in the file. This would be very helpful if one prepares in Chimera PDB inputs (containing also non protein or non DNA/RNA molecules) for Amber because Amber (Leap) automatically tries to create bonds between all consecutive residues unless they are separated by TER string. Thank you very much for any useful information. Best wishes, Marek -- Tato zpráva byla vytvořena převratným poštovním klientem Opery: http://www.opera.com/mail/

3 11

error in add chain ID
by Hamzeh, Rahimi 07 Jun '12

07 Jun '12

Dear Eric use I want rename_pdb_for_chimera.py for adding chain ID to pdb file but i got following erro. buksanbio@hamzeh:/media/6C166CAD166C7A4A/kit/zdock/apc10-apc3$ python rename_pdb_for_chimera.py 1 avgapc10final.pdb File "rename_pdb_for_chimera.py", line 100 there's a chain identifier in the "TER" line... (the condition len ^ IndentationError: expected an indented block Please let me know about my wrong best wish

2 1

Re: [Chimera-users] undisplay label in distance monitor
by Bala subramanian 07 Jun '12

07 Jun '12

Friends, > I have written a script to draw lines between atoms (using distance > monitor). I dnt want the distance label to be displayed with the pseudo > bond. Kindly write me the way to remove the label. My script is below. > > model=chimera.openModels.open('myavg.pdb') > res=model[0].residues > MAT=loadtxt('mymatrix',dtype=float) > chimera.runCommand('focus') > for x in range(len(res)): > for y in range(len(res)): > if x <> y : > value=MAT[x,y] > > if 0.5 < value > 0.75: > > b=distanceMonitor.newPseudoBond(res[x].atomsMap['CA'][0],res[y].atomsMap['CA'][0]) > b.drawMode=1 > b.radius=0.05 > #b.label=None I tried keeping this value as None and empty > string but it doesnt help. > b.color=getColorByName('red') > > else: pass > > else: continue > > Thanks, Bala -- C. Balasubramanian

2 3

index out of range
by Francesco Pietra 06 Jun '12

06 Jun '12

Running movie for prmtop/dcd proved extremely slow at amd64, finally showing the protein in a water box, albeit offering very limited analysis. Error message: index out of range, error processing "monitor changes". The log is 16MB!, below the startr and end of the log: MDTools 1.2 Exp (2005-08-20) by James Phillips. Trajectory.DCD.MDToolsMarch97.md.help() for more info. HomoCoord 1.3 Exp (2005-08-20) - Constants 1.1 Exp (2004-05-17) - AtomGroup 1.4 Exp (2007-02-07) - local 1.1 Exp (2004-05-17) Unknown OSTYPE ??? set. Viewing and plotting are disabled. - Trans 1.2 Exp (2005-08-20) - DCD 1.10 Exp (2007-02-07) - Data 1.2 Exp (2005-08-20) using memory mapping REMARKS FILENAME=./heat-01.dcd CREATED BY NAMD REMARKS DATE: 06/03/12 CREATED BY USER: francesco Model 0 (heat-01.dcd) appears to be a protein without secondary structure assignments. Automatically computing assignments using 'ksdssp' and parameter values: energy cutoff -0.5 minimum helix length 3 minimum strand length 3 Use command 'help ksdssp' for more information. Traceback (most recent call last): File "/usr/local/chimera_Nov3.2011/share/chimera/triggerSet.py", line 83, in invoke self._funcData, triggerData) File "/usr/local/chimera_Nov3.2011/share/chimera/Sequence.py", line 1161, in _cacheCheck if len(kids) == molecule.rootForAtom(kids[0], IndexError: list index out of range Error processing trigger "monitor changes": .................................... .................................... Traceback (most recent call last): File "/usr/local/chimera_Nov3.2011/share/chimera/triggerSet.py", line 83, in invoke self._funcData, triggerData) File "/usr/local/chimera_Nov3.2011/share/chimera/Sequence.py", line 1161, in _cacheCheck if len(kids) == molecule.rootForAtom(kids[0], IndexError: list index out of range Error processing trigger "monitor changes": IndexError: list index out of range File "/usr/local/chimera_Nov3.2011/share/chimera/Sequence.py", line 1161, in _cacheCheck if len(kids) == molecule.rootForAtom(kids[0], See reply log for Python traceback. Thanks for the attention francesco pietra

2 2

transparency for surface in intersurf
by Jean-Didier Maréchal 06 Jun '12

06 Jun '12

Hi guys, is there a way to change the transparency of intersurf objects? I'd like to have the surface not hiden all the object in behing one given view. thanks a lot, cheers jd Dr. Jean-Didier Maréchal Lecturer The Computational Biotechnological Chemistry Team 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 www: http://asklipio.qf.uab.es

2 2

fitmap
by Hernando J Sosa 06 Jun '12

06 Jun '12

HI, What is the correct syntax to run the fitmap command?. I have a pdb structure named C1 and a density map (spider format) named MV (names as appear in the Model panel). If I try in the command line: fitmap C1 MV I get the error: fitmap invalid volume argument "MV" What am I doing wrong? Thanks

3 5

Fwd: ESP calculations
by Eric Pettersen 05 Jun '12

05 Jun '12

Begin forwarded message: > From: Nikolay Igorovich Rodionov <nirodion(a)syr.edu> > Date: May 21, 2012 12:43:34 PM PDT > To: Eric Pettersen <pett(a)cgl.ucsf.edu> > Subject: RE: [Chimera-users] ESP calculations > > Thank you both, > > Can either of you tell me how the information in a volume data set > is organized? I am looking at a .dx file that I created using APBS > which I would assume to be very similar to the volumetric data sets > created by Chimera. What does each column of data represent? Is > there any way for me to sort the ESP data based on position using > Excel or Matlab? I am only concerned with ESP data for one region of > molecule. If you could point me in the direction of some good > documentation on this file type it would be great. I don't know, so I'm forwarding this out to chimera-users land, where someone might have an answer... --Eric > > Thanks again, > Nikolay Rodionov > > -----Original Message----- > From: Eric Pettersen [mailto:pett@cgl.ucsf.edu] > Sent: Monday, May 21, 2012 2:57 PM > To: chimera-users(a)cgl.ucsf.edu BB > Cc: Nikolay Igorovich Rodionov > Subject: Re: [Chimera-users] ESP calculations > > Also, if you compute a volume, you can use the "Values at Atom > Positions" tool to, um, get the values at atom positions (say, a > ligand), if that's relevant to what you're doing. > > --Eric > > On May 21, 2012, at 9:09 AM, Elaine Meng wrote: > >> Hi Nikolay, >> You can generate a grid of values with the "Compute grid" option of >> Coulombic Surface Coloring (or corresponding option of the >> "coulombic" >> command). The grid is a volume data set, and when it is generated >> the >> Volume Viewer and Surface Color dialogs will automatically appear. >> >> <http://www.cgl.ucsf.edu/chimera/docs/ContributedSoftware/coulombic/ >> co >> ulombic.html >>> >> <http://www.cgl.ucsf.edu/chimera/docs/ContributedSoftware/coulombic/ >> co >> ulombic.html#volume >>> >> >> You can save the grid to a file with the FIle menu of Volume Viewer. >> <http://www.cgl.ucsf.edu/chimera/docs/ContributedSoftware/ >> volumeviewer >> /framevolumeviewer.html >>> >> >> As for the values at each point of the surface, I don't know of a way >> to write these all to a file, but in Surface Color if you click the >> Options button and then turn on the option to "Report value at mouse >> position" it will report the values in the status line at the bottom >> of the Chimera window when you click into the graphics window and >> then >> mouse over the surface. I just noticed that to make this work it is >> necessary to recolor the surface first by clicking the Color button >> on >> this dialog, which may be a bug, but that's the workaround! >> <http://www.cgl.ucsf.edu/chimera/docs/ContributedSoftware/surfcolor/ >> su >> rfcolor.html >>> >> >> If you have the grid file saved you can reopen it in Chimera later >> and >> use it in Surface Color without rerunning the ESP calculation. >> >> However, keep in mind that Coulombic ESP is more for qualitative >> purposes of visualization or making comparisons between related >> structures than for absolute quantitative accuracy. For example, you >> could easily make all the ESP values half as large simply by doubling >> the dielectric constant used in their calculation, as shown in the >> equation in the Coulombic Surface Coloring manual page. >> >> 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 May 21, 2012, at 8:05 AM, Nikolay Igorovich Rodionov wrote: >> >>> Hi all, >>> I was wondering if there was a way to get an output data file >>> regarding ESP coloring and calculations based on Coulomb's law. The >>> visualization is great but I would also like to quantify the data. >>> Thank you, >>> Nikolay Rodionov >> >> >> _______________________________________________ >> Chimera-users mailing list >> Chimera-users(a)cgl.ucsf.edu >> http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users >> > > > >

3 2

ESP calculations
by Nikolay Igorovich Rodionov 05 Jun '12

05 Jun '12

Hi all, I was wondering if there was a way to get an output data file regarding ESP coloring and calculations based on Coulomb's law. The visualization is great but I would also like to quantify the data. Thank you, Nikolay Rodionov

4 3

Chimera-users June 2012