March 2011 - Chimera-users - RBVI Mailing Lists

matchmaker output
by Tiziano Tuccinardi 18 Mar '11

18 Mar '11

Dear all, I am using MatchMaker for aligning about 100 structures, and it works very well!! Is it possible to automatically save the pairwise alignment without opening the MultAlignViewer? Due to the large number of structures, I cannot use the Match-Align tool but it would be enough to save the pairwise alignments Thanks in advance Tiziano -- +++++++++++++++++++++++++++++++++++++ Tiziano Tuccinardi, PhD Dip. Scienze Farmaceutiche Via Bonanno 6, 56126 PISA Tel ++39 050 2219595 Fax ++39 050 2219605 E-mail tuccinardi(a)farm.unipi.it <mailto:tuccinardi@farm.unipi.it> http://www.mmvsl.farm.unipi.it/staff/tiziano-tuccinardi/ Adjunct Assistant Professor SbarroInstitute for Cancer Research and Molecular Medicine Centerfor Biotechnology TempleUniversity PhiladelphiaPA 19122 (USA) +++++++++++++++++++++++++++++++++++++

2 3

Movie frames are blank
by Gökhan Tolun 17 Mar '11

17 Mar '11

Dear all, I am trying to create a very simple movie, but my movies turn out blank. I am using the Spin Movie example script on http://www.cgl.ucsf.edu/chimera/tutorials/movies08/moviemaking.html My movie seems to render, and the movie file is created, but when I open the file and play it, its just a blank (all black) video. When I examined the captured frames, theyre also all blank black frames. If I choose Save image, it saves it just fine, so I guess the image capturing algorithm for movies is different the one for saving the images. This is on Windows 7 64-bit, Chimera 1.5.2, Intel Mobile 965 express Chipset GPU (its a TabletPC, hence the GPU). Any suggestions to fix this would be appreciated. Thank you. _______________________________________ Gökhan Tolun, M.S., PhD. Post-doctoral Research Associate Lineberger Comprehensive Cancer Center University of North Carolina 450 West Dr., CB#7295, Room 11-119 Chapel Hill, NC, USA 27599

2 1

Average Backbone B-factor
by Simon Sham 17 Mar '11

17 Mar '11

Hi, I need to read a pdb file and calculate the average backbone (CA, C, & N) B-factor for each residue and compare it with the result some MD. Is there anyway I can do it in Chimera? Thanks, Simon

2 1

Re: [Chimera-users] error: from DockPrep import prep
by wascm 16 Mar '11

16 Mar '11

Dear sir: Thanks for the detailed reply! When the following code was added to the start of sciprt: import chimeraInit chimeraInit.init(nogui=True, silent=True) The interpreter print the following confused error: --------------------------------------------- Traceback (most recent call last): File "./xx.py", line 5, in <module> chimeraInit.init(nogui=True, silent=True) TypeError: init() takes at least 1 argument (2 given) ----------------------------------------------- it's really easy and convenient to running python script with chimera. The reason for i want to run python script independently is that i want add some custom command option. However, if i run python script with chimera, it's hard for me to add some custom command option and handle it in the script though getopt module. For example, when i run some command like this: chimera –nogui test.py -t test test.py ################### #!/usr/bin/python t=None try: opts, args = getopt.getopt(sys.argv[1:], "t:") except getopt.GetoptError, err: pass finally: for arg_k, arg_v in opts: if (arg_k =="-t"): t=arg_v #################################### This will raise running exception! is there a way for me to add the coustom option to python script and run the script with chimera? the exception message as follow: ------------------------------------------------- Executing test.py... Traceback (most recent call last): File "/home/wascm/chimera/share/chimeraInit.py", line 588, in init chimera.openModels.open(a, prefixableType=1) File "/home/wascm/chimera/share/chimera/__init__.py", line 1537, in open models = func(filename, *args, **kw) File "/home/wascm/chimera/share/chimera/__init__.py", line 972, in _openPython loadFunc(sandboxName, fileName, f) File "test.py", line 10, in <module> for arg_k, arg_v in opts: NameError: name 'opts' is not defined Error while processing test.py: NameError: name 'opts' is not defined File "test.py", line 10, in <module> for arg_k, arg_v in opts: See reply log for Python traceback. Traceback (most recent call last): File "/home/wascm/chimera/share/chimeraInit.py", line 588, in init chimera.openModels.open(a, prefixableType=1) File "/home/wascm/chimera/share/chimera/__init__.py", line 1503, in open raise ValueError, "Unknown model type" ValueError: Unknown model type Error while processing -t: ValueError: Unknown model type File "/home/wascm/chimera/share/chimera/__init__.py", line 1503, in open raise ValueError, "Unknown model type" See reply log for Python traceback. Traceback (most recent call last): File "/home/wascm/chimera/share/chimeraInit.py", line 588, in init chimera.openModels.open(a, prefixableType=1) File "/home/wascm/chimera/share/chimera/__init__.py", line 1503, in open raise ValueError, "Unknown model type" ValueError: Unknown model type Error while processing test: ValueError: Unknown model type File "/home/wascm/chimera/share/chimera/__init__.py", line 1503, in open raise ValueError, "Unknown model type" See reply log for Python traceback. ---------------------------------------------------------- At 2011-03-15 09:02:09，"Greg Couch" <gregc(a)cgl.ucsf.edu> wrote: (This is really achimera-dev(a)cgl.ucsf.edu mailing list question.) My guess is that your problem is that chimeraInit.init() was never called, so the chimera internals were never set up. So in the beginning of your Python script, add: import chimeraInit chimeraInit.init(nogui=True, silent=True) before using chimera's modules. But, in general, you would be better off running those python scripts with chimera. For example: % chimera --nogui --silent script.py That would set up the the appropriate environment for chimera, execute the python script, and exit. And you would be protected against any future changes in how chimera is initialized. If you want to use other Python modules with chimera, set your CHIMERAPATH environment variable with a list of directories that hold the Python modules (same syntax as PYTHONPATH, but just for chimera). If you want to install other Python packages into chimera's Python that have a setup.py file, use: % chimera --nogui --silent setup.py install It is also possible to make easy_install work with chimera by temporarily putting chimera's bin directory on your path and running the appropriate egg (seehttp://pypi.python.org/pypi/setuptools for details). Right now, chimera is using Python 2.7. HTH, Greg

3 2

radial clipping
by Kelly Dryden 16 Mar '11

16 Mar '11

Hi, I would like to radially clip a density map. There are nice ways to do radial coloring, but I haven't found anything that allows a similar cropping (preferably with capping). Is there an easy way to do this? I do not want to erase the density, but rather be able to see the surface at varying radii, just as you can do with planes. (At the moment, I can only run v1.4). Thank you, Kelly

3 5

Re: [Chimera-users] [chimera-dev] error: from DockPrep import prep
by wascm 15 Mar '11

15 Mar '11

Dear sir: Your suggestions do work well! The problems have been resolved well by your help. Thank you very much! Best Wishes Tao Cui From:Eric Pettersen Date:2011-03-16 05:12 Subject:Re: [chimera-dev] [Chimera-users] error: from DockPrep import prep To:wascm Cc:Greg Couch,chimera-dev BB,UCSF Chimera Mailing List On Mar 15, 2011, at 5:37 AM, wascm wrote: Dear sir: Thanks for the detailed reply! When the following code was added to the start of sciprt: import chimeraInit chimeraInit.init(nogui=True, silent=True) The interpreter print the following confused error: --------------------------------------------- Traceback (most recent call last): File "./xx.py", line 5, in <module> chimeraInit.init(nogui=True, silent=True) TypeError: init() takes at least 1 argument (2 given) There was an oversight in the code Greg provided: chimeraInit.init takes a mandatory "argv" argument, which is exactly like sys.argv. ----------------------------------------------- it's really easy and convenient to running python script with chimera. The reason for i want to run python script independently is that i want add some custom command option. However, if i run python script with chimera, it's hard for me to add some custom command option and handle it in the script though getopt module. For example, when i run some command like this: chimera –nogui test.py -t test test.py ################### #!/usr/bin/python t=None try: opts, args = getopt.getopt(sys.argv[1:], "t:") except getopt.GetoptError, err: pass finally: for arg_k, arg_v in opts: if (arg_k =="-t"): t=arg_v #################################### This will raise running exception! is there a way for me to add the coustom option to python script and run the script with chimera? Yes. Use a --script argument, like this: chimera --nogui --script "test.py -t test" --Eric

1 0

I am having a problem
by Burt Goldberg 15 Mar '11

15 Mar '11

Dear, Chimera, as I opened the program it prompted me to re-register. I tried 3 time but got a Error message PKCS7. HELP!!!! Burt -- Burt Goldberg, PhD Professor of Chemistry & Biochemistry Chemistry Department Affiliated Research Professor Medical Parasitology New York University School of Medicine New York University

2 1

Re: [Chimera-users] Selection through bond and fragment naming
by Jean Didier Pie Marechal 14 Mar '11

14 Mar '11

Dear all, thanks a lot for your answers. I seem to have a problem though. When running the script, one way or another I always ends up with something like: a.residue.removeAtom(a) AttributeError: 'NoneType' object has no attribute 'removeAtom' AttributeError: 'NoneType' object has no attribute 'removeAtom' File "M:\REUNIONS\dibuixos\Manipulate_structure\newres.py", line 20, in <module> a.residue.removeAtom(a) See reply log for Python traceback. I think I got the origin of the problem but can't figure out how to solve it. What I observed is that, before running the script if I do the following after selecting my atoms: >>> for a in at: print a.residue I have the following list printed UNK 1 UNK 1 UNK 1 UNK 1 UNK 1 UNK 1 UNK 1 UNK 1 UNK 1 UNK 1 UNK 1 UNK 1 UNK 1 UNK 1 UNK 1 UNK 1 UNK 1 UNK 1 UNK 1 UNK 1 but doing the same AFTER having run the script gives : None UNK 1 UNK 1 UNK 1 UNK 1 UNK 1 ... It seems like the residue type of the first atom disappeared. Any clue? all the best, 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 ----- Missatge original ----- De: Eric Pettersen <pett(a)cgl.ucsf.edu> Data: Divendres, Març 11, 2011 8:29 pm Assumpte: Re: [Chimera-users] Selection through bond and fragment naming > On Mar 11, 2011, at 7:42 AM, Jean Didier Pie Marechal wrote: > > > Hi, > > > > I have a system from which I only have xyz coordinates. What I > would > > like to do is to define myself residue attribute on different > > fragments of the molecules. Could you tell me how (or if there is > a > > way to ): > > > > 1. expand a selection to bonding atoms (i.e. I click on one atom > and > > expand selection through one bond after each actions) This would > be > > of the easiest way for me to select the fragments I want to. > > Hi JD, > My only idea as to how to do this is with a short Python script: > > from chimera.selection import currentAtoms, addCurrent > addAtoms = set() > for a in currentAtoms(): > addAtoms.update(a.neighbors) > addCurrent(addAtoms) > > You could then use the open command to run it: > > open <script folder location (e.g. ~marechal/scripts)>/expand.py > > which you could short to "xp" with the alias command: > > alias ^xp open blah-blah/expand.py > > > > > 2. set a residue attribute to the given selected fragment. I > think I > > remember how to do so in scripting, but what is the command line > to > > do so? > > So, do you really want to set a residue-level attribute, or do you > want to organize the selected atoms into a residue? Setting a > residue- > level attribute for the atoms is easy: > > setattr r myAttrName my-attr-val sel > > Organizing the atoms into a residue again requires Python, and not > as > simple a script as for your first question, e.g.: > > # script designed to run under "runscript" command, therefore has > an > "arguments" variable available... > try: > rname, rnum = arguments > except ValueError: > raise chimera.UserError("must supply residue name and number > arguments") > > try: > rnum = int(rnum) > except ValueError: > raise chimera.UserError("residue number (second arg) must be an > integer") > > atoms = chimera.selection.currentAtoms() > if not atoms: > raise chimera.UserError("No atoms selected") > > # make a new residue > newRes = atoms[0].molecule.newResidue(rname, "A", rnum, ' ') > > # withdraw atoms from their current residues > for a in atoms: > a.residue.removeAtom(a) > if len(a.residue.atoms) == 0: > a.molecule.deleteResidue(a.residue) > > # add atoms to new residue > for a in atoms: > newRes.addAtom(a) > > So you would run this script with the command "runscript blah- > blah/ > newres.py TYR 1", which could again be shorted using aliasing: > > alias ^nr runscript blah-blah/newres.py $1 $2 > > --Eric > > > Eric Pettersen > UCSF Computer Graphics Lab > http://www.cgl.ucsf.edu > > >

2 1

registration and user group
by SteffenVfx 14 Mar '11

14 Mar '11

Hi there. I've been trying to register with you all for some time, but it always comes back with a registration error. I'd like to enroll in your mailing list and user group forum. can you pass along my information? thanks casey *Steffen VFX ~the science of storytelling~* *| Casey Steffen M.S. | President/Chief Creative Officer | 211 6th Ave. #2a | Brooklyn, NY 11215 | **Web:**www.steffenvfx.com | **Twitter:**SteffenVfx<http://twitter.com/steffenvfx> | **Blog:**http://blog.steffenvfx.com | **Phone:347-875-0423 |*

2 1

bond selection
by Elisabeth Ortega 11 Mar '11

11 Mar '11

Dear Chimera Dev. Team, I have a little question for you. How can I select a bond using the Command Line (or the IDLE)? Thanks in advance. Elisabeth Ortega-Carrasco Universitat Autònoma de Barcelona

3 2