- Chimera-users - RBVI Mailing Lists

MD setup from command line?
by Dr. techn. Thomas Mitterfellner 04 Aug '11

04 Aug '11

Hello! I'm using Chimera to view DFT optimizations and vibrational modes. I split the large xyz files into multiple small xyz files and load them with the MD tool. Now my question is: can I do the last step from the command line so I could write a script which sets the first and the last frames automatically? Also, is it possible to create the respective movie (as mpeg) in an automated fashion? Thank you, Thomas

3 5

Fwd: [3DEM]asymmetric unit segmentation
by C.J. 04 Aug '11

04 Aug '11

Dear all, I post here for the suggestion of Steven Ludtke in the 3DEM list. best regards! C.J. ---------- Forwarded message ---------- From: C.J. <biocjh(a)gmail.com> Date: 2011/7/13 Subject: [3DEM]asymmetric unit segmentation To: 3dem(a)ncmir.ucsd.edu Dear all, Any one have a experience on segmentation of asymmetric unit from icosahedron virus structure volume? I wonder whether there is any programme could process that automatically? Any comment would be appreciated! -- Regards! C.J. -- Regards! C.J.

3 3

Torsion adjust causes weird behavior
by Greg Friedland 03 Aug '11

03 Aug '11

Hi, I'm attempting to adjust a side chain torsion by hand for an asparagine residue but when I change it in the 'adjust torsion' panel, another residue (only one) 5 angstroms or so away is rotated around the bond axis as well. I made sure that the dihedral bond vector in the panel is pointing away from the backbone pointing in the direction and that only one torsion is selected. Is there a way I can change this behavior to only adjust the relevant residue? Thanks, Greg

2 1

Announcement: Release of DOCK 6.5
by Scott Brozell 01 Aug '11

01 Aug '11

We are pleased to announce the release of DOCK 6.5. DOCK is a suite of programs for molecular docking. In version 6.5 a new scoring function Descriptor Score is available. It reports van der Waals and Coulombic energies, hydrogen bonds, and per-residue energy decomposition comparisons with a reference molecule. DOCK continues to provide a variety of other scoring methods including grid based, PB/SA, and GB/SA, as well as a simple interface to the AmberTools MD engine that enables receptor flexibility and a complete MM force field. For full information on what is new in DOCK 6.5, please visit: http://dock.compbio.ucsf.edu/DOCK_6/new_in_6.5.txt Sincerely, The DOCK Team Please visit us at the DOCK Web site. http://dock.compbio.ucsf.edu

1 0

centroid trace problems
by Francesco Pietra 01 Aug '11

01 Aug '11

hello: "centroid" works fine to me in getting the trace of a moving molecule in and out of a protein, using command define centroid radius 0.2 color red :residuename when the residue is unique. In fact, the method has produced extraordinarily illustrative (I am using the words of the referee) in a paper that will appear shortly on Chemistry & Biodiversity. Now I have three diatomic molecules, only one of which is made moving around. In this case define centroid radius 0.2 color red :residuename:residuenumber for the moving residue does not work. What happens is that "red balls" for the specified residue are being accumulated at the place of origin of the residue, while the residue is leaving the protein as a short bar, not a ball. Could that be fixed in the above command? thanks francesco pietra

3 5

hide a peptide in Chimera
by C.J. 29 Jul '11

29 Jul '11

Hi,all, I have a pdb file containing several peptides. And I want to hide some of them when visualizing the file in Chimera. Anyone would like to help me to that? I'm appreciated of any suggestions. -- Regards! C.J.

2 3

Hbond in CLI + cacheDA
by Maciek Wójcikowski 29 Jul '11

29 Jul '11

Hello everyone, I'm trying to compute hbonds for quite large molecular database, so i do it in CLI. I've added cacheDA parameter which speeds up whole process at first by 10 times, although after time cache is getting bigger and bigger operations slows down as one can expect, because Chimera is caching every compound. Is there a way to limit the size of a cache, or even better to tell chimera to cache only protein donors and acceptors? ---- Pozdrawiam, | Best regards, Maciek Wójcikowski maciek(a)wojcikowski.pl

2 6

Chimera 1.5.3, nvidia quadro 4000 and stereo
by Veli-Pekka Kestilä 27 Jul '11

27 Jul '11

Hello, I found following problem when trying to use Chimera 1.5.3 with stereo in Nvidia Quadro 4000. When starting the chimera in windows I will get error: Error Initializing OpenGL: Couldn't configure togl widget When stating in OpenGL debug mode it starts without any problem. Nvidia driver is 275.65 and everything works in linux without problem. In Manage 3D settings I have put on Stereo Enable. Greetings, Veli-Pekka Kestilä -- University of Helsinki / Institute of Biotechnology / veli-pekka.kestila(a)helsinki.fi / +358 9 191 58921 / Room 4316, Biocenter 3

3 4

gtx-470 GeForce cards
by Francesco Pietra 26 Jul '11

26 Jul '11

With GeForce2 MX200 / Debian GNU-Linux i386 wheezy I have no problems in running movie with "define centroid". In contrast, with GeForce gtx-470 (two such cards and PhenomII 1075T (six cores)/ Debian amd64 wheezy, on running movie to get the trace of the centroid, the X server hangs and the keyboard is no more sensed. I have to restart the computer from a ssh-linked desktop. In both cases chimera 1.5.3. Chimera was started as the only application, just after booting the OS (I am saying that because, to run molecular dynamics, I have first to command # nvidia-smi -L # nvidia-smi -pm 1 Any experience to this regard? thanks francesco pietra

2 1

castp .poc file is not opened in restored session (plain text version)
by Boaz Shaanan 26 Jul '11

26 Jul '11

Hi, I save a session in which I used castp .poc file (very nice indeed!). However, after saving the session, the restored session does not open the .poc file. Is that how it's supposed to be? Cheers, Boaz P.S. sent again in plain text format since it seems the previous html version has been mangled up Boaz Shaanan, Ph.D. Dept. of Life Sciences Ben-Gurion University of the Negev Beer-Sheva 84105 Israel E-mail: bshaanan(a)bgu.ac.il Phone: 972-8-647-2220 Skype: boaz.shaanan Fax: 972-8-647-2992 or 972-8-646-1710

2 4

Query regarding Ksdssp assignment
by Arun Prasad Pandurangan 25 Jul '11

25 Jul '11

Dear Users and Developers, I would like to know whether its possible to retrieve the full secondary structural assignment using Ksdssp within Chimera. For instance, I would like to know all the strands that belong to a beta sheet. If I am right, running Ksdssp from Chimera only gives details of the strand (using isStrand) and helix (using isHelix) assignment and not about sheet information. Using the full output of an another program DSSP, I am able deduce which strands belongs to which sheets. I would like to know whether similar kind of detailed output is available from Ksdssp as well. If so how to retrieve those information within Chimera? Looking forward for the advices. Thank you. with best regards, Arun Prasad

3 4

Display of bond between CA's
by Greg Friedland 25 Jul '11

25 Jul '11

Hi, I'm trying to prepare an image from a structure in which I'd like to hide the backbone atoms except the CAs. When I do this for two neighboring residues however a long bond is drawn between the CAs. How can I hide this without hiding the CAs? Thanks, Greg

2 1

creating multimeric units
by ALEXANDER MUCHUGIA 25 Jul '11

25 Jul '11

Hi, My name is Alexander Muchugia, I am a Bioinformatics student at Rhodes University, South Africa and I am working on modelling HIV protease. It is a dimer but Its modelled as a monomer. I was trying to look for software than can convert it as a dimer. I am enquiring if it is possible to use chimera to do it and how can I do it. Regards. Alexander.

2 1

Chimera for Mac OSX Lion
by Victoria James 22 Jul '11

22 Jul '11

Hi there! I was just wondering if you know if the latest build of Chimera will run on Mac OSX Lion? If not, will you be releasing another build soon that will be able to? Many thanks and kind regards Victoria

5 4

Re: [Chimera-users] Structure -> Match in a script
by Eric Pettersen 21 Jul '11

21 Jul '11

On Jul 20, 2011, at 2:09 PM, Sierk, Michael wrote: > Eric, > > Thanks for the script! It works but I have a couple of issues still: > > 1) The second sequence in the alignment is not automatically being > associated with the second structure. (It looks like the first > sequence/structure association is done automatically, but not the > second.) I can do it manually with Structure->Associations (with > zero mismatches), and then the script works. Is there a way to > create the second association via either the Chimera command line or > using python? Is that association supposed to occur automatically > as well? AFAIK the second sequence should associate. If you send me the alignment and the structures (or tell me their PDB codes if they're standard PDB entries) I can investigate (send them off-list!). > 2) What I am trying to do is load in two molecules and 3 different > sequence alignments of the two molecules. In my Chimera script I > make 2 copies of the 2nd molecule (using "combine") and color them > differently. I would like to then do the Structure->Match process > for each of the 3 alignments. Does the "findMAV" function in the > mavMatch.py script return all instances of MAVs or just one? I know > the titles so could I perform 3 different mav.match statements, one > for each MAV instance? Well, this is a more complicated scenario than what the script I wrote really handles -- since you have multiple MAV instances and only want certain open structures to associate with each instance despite the fact that those structures are identical. So it needs two scripts! Here's the approach: In your Chimera script: 1) Open the three alignments. 2) Open the first Python script (no-auto-assoc.py) to turn off auto- association in the MAV instances. 3) Open your two structures, "reference" structure first (more specifically, in a lower model number than the matching structures) 4) Make the two copies, or just open the second structure two more times. 5) Run the second Python script (assoc-and-match.py) to associate the reference structure with the first sequence in all the alignments, and one copy of the matching structure with the second sequence in each alignment, and then match them. --Eric

1 0

problem with a script calculating distance between residues
by Damien Larivière 20 Jul '11

20 Jul '11

Dear all, The script below has been validated for me by the Chimera team sometimes ago. In fact this script does not work each time I run it: When reaching the pdb dockSL0007.pdb, or dockSL0010.pdb, or dockSL0020.pdb,... Chimera can remained fixed on the current PDB file as if it faces problem to open it. I'm forced to close Chimera which does not respond anymore. The exit code is 805306369. I attached dockSL0009.pdb and dockSL0010.pdb. What can I do to report this bug (which may be due to something external to Chimera; my Chimera version is 1.5.2 Build 32411; my plateform is Windows 7 64 bits) since I have to brutally stop Chimera, and know more about this problem? Many thanks for your help Damien from chimera import runCommand as rc #from chimera.tkgui import saveReplyLog for i in range(1, 101): rc(r"open C:\Program Files (x86)\Hex 6.3\examples\New Hex run_SL_250211\FirstRun_S_Ldimer\dockSL%04d.pdb" % i) rc("distance :246.C@CB :297.A@CB") rc("distance :246.C@CB :297.B@CB") rc("close all") #rc("saveReplyLog('replyLog%04d.txt')" % i)

3 2

Structure -> Match in a script
by Sierk, Michael 20 Jul '11

20 Jul '11

All - I would like to write a command file that will automatically read in two structures and a corresponding alignment file, then perform the superposition (and coloring, etc.). I can read in the pdb and pir files using the command line, and I can do the superposition in the Multalign Viewer using Structure->Match, but I cannot figure out how to do the superposition based on the sequence using the command line. I think if there was a way to select the residues that are aligned in the sequence alignment I could use the "match" command, but I don't see a way to do that. Is there an attribute analogous to mavPercentConserved that corresponds to "aligned"? Or is there an automatically generated region based on the sequence alignment? Thanks, Mike Sierk +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ Michael Sierk Email: michael.sierk<AT>email<dot>stvincent<dot>edu Assistant Professor of Bioinformatics Phone: (724) 805-2367 Saint Vincent College Office: Physics 205 ----------------------------------------------------------------------------------------------------------------------------------------- "There is something fascinating about science. One gets such wholesome returns of conjecture out of such a trifling investment of fact." -- Mark Twain. +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

2 1

Re: [Chimera-users] Listing hydrogen bond geometries
by Elaine Meng 18 Jul '11

18 Jul '11

Hi Greg, For FindHBond, Chimera "imagines" the possible locations of hydrogens given the atom types (element, hybridization state) of the heavy atoms. <http://www.cgl.ucsf.edu/chimera/docs/ContributedSoftware/findhbond/findhbon…> "When there are no explicit hydrogens on a potential donor atom, FindHBond will use the donor atom type to infer their presence. For functional groups where the hydrogen positions are well-determined by the positions of the nonhydrogen atoms, the H-bonds detected in the presence and absence of explicit hydrogen atoms should be virtually identical. For functional groups where the hydrogen position is not well-determined (such as rotatable hydroxyls), the presence of an explicit hydrogen atom will limit the detected H-bonds to those appropriate for the observed position, instead of all of its possible positions." If you want explicit hydrogens, you could use Chimera's AddH tool (under Tools... Structure Editing) or addh command. By default, it will look at the surroundings and try to place hydrogens to make hydrogen bonds. However, it does not exhaustively optimize the H-bonding network of the whole structure. The goal is to determine reasonable positions based on immediate surroundings in a fairly short time. However, if you add explicit hydrogens and then use FindHBond, there will be fewer possibilities listed, as explained above. <http://www.cgl.ucsf.edu/chimera/docs/ContributedSoftware/addh/addh.html> For the next level of computational intensity in adding hydrogens, we recommend the "Reduce" program from the Richardson lab. From the bottom of the AddH page, "When a more intensive approach is desired, the program Reduce (developed by the Richardson Laboratory) is a good alternative. Reduce places hydrogens to optimize local H-bonding networks and avoid steric overlaps, while flipping certain sidechains 180 degrees as deemed appropriate to fulfill these criteria. Asparagine and glutamine sidechains may be flipped to switch their terminal N and O atoms, and the imidazole ring of histidine may be flipped to switch N and C identities. The protonation state of histidine is adjusted based on the local environment. Reduce is available free at http://kinemage.biochem.duke.edu: • for on-line use as part of the MolProbity service, on a file uploaded or chosen by PDB code (individual sidechain flips can be accepted or rejected) • for download (from the Software section) to run on most platforms and is described in: Asparagine and glutamine: using hydrogen atom contacts in the choice of side-chain amide orientation. Word JM, Lovell SC, Richardson JS, Richardson DC. J Mol Biol. 1999 Jan 29;285(4):1735-47. " I hope this helps, Elaine On Jul 18, 2011, at 2:16 PM, Greg Friedland wrote: > Hi Elain, > Thanks very much for you response and sorry for the delay in getting back to you. > > I am trying to determine the potential hydrogen bonding potential (i.e. the ability to form these bonds) from hydrogen-less structures (they are from the PDB). I have been trying to model in the hydrogens but it is not trivial since there are many hydroxyls, which can rotate about their dihedrals. (if you have suggestions on how to model these hydrogens it'd be great to know about them). > > So until I can get a good model with the hydrogens attached I am wondering how Chimera is able to determine potential hbonds from heavy atom positions only. Any insight you can provide would be appreciated. > > Cheers, > Greg

1 0

Listing hydrogen bond geometries
by Greg Friedland 15 Jul '11

15 Jul '11

Hello, I would like to list the hydrogen bond geometries between two models. I have found the feature to do this in the FindHBond panel but it only lists the distances not the angles. This information appears to be used in identifying the hydrogen bonds however, so is it possible to extract this information as well? If not, can you suggest another tool that could do this? Thanks, Greg

3 2

(no subject)
by efratk＠techunix.technion.ac.il 14 Jul '11

14 Jul '11

Hello, I am a Phd student at the Technion, I use the program- DMS and I would like to know what is the numeration of the columns of the first six fields in the output format of the dms. Thanks in advance, Efrat

3 4

(no subject)
by efratk＠techunix.technion.ac.il 14 Jul '11

14 Jul '11

Hello, I am a Phd student at the Technion, I use the program- DMS and I would like to know what is the numeration of the columns of the first six fields in the output format of the dms. Thanks in advance, Efrat

1 0

saving coordinates from a session and re-opening them in another one
by Boaz Shaanan 14 Jul '11

14 Jul '11

Hi I save selected & displayed coordinates (using save pdb) from a session, in order to incorporate them into another session. It’s a selected zone around the active site. When I open them in the new session they look strange. Also, the saved coordinates don’t always contain the full set of displayed coordinates (I’ve done this several times on different structures and it seems particularly but not just on pdb’s with aniso card). I’m attaching the session whose coordinates I saved + the saved coordinates (I hope that’s OK with the BB policy). Please help me figure out the problem. Thanks, Boaz Boaz Shaanan, Ph.D. Dept. of Life Sciences Ben-Gurion University of the Negev Beer-Sheva 84105 Israel E-mail:bshaanan@bgu.ac.il Phone: 972-8-647-2220 ; Fax: 646-1710 Skype: boaz.shaanan

2 5

DMS output format
by Efrat Kligun 13 Jul '11

13 Jul '11

Hello, I am a Phd student at the Technion, I use the program- DMS and I would like to know what is the numeration of the columns of the first six fields in the output format of the dms. Thanks in advance, Efrat

2 1

3D chimera
by Bigley, Andrew 11 Jul '11

11 Jul '11

Dear Chimera Developers, The recent advances in 3D technology is affording some unique opportunities to people interested. I have recently seen advertisements for new glasses free 3D smart phones. Can you adapt chimera for such a platform and hence convert a smart phone into a portable 3D protein viewer that would allow manipulation of the view through typical touch screen commands? Andrew N. Bigley, Ph.D.

2 1

Curved helical filaments
by Tom Goddard 11 Jul '11

11 Jul '11

Hi Andy, The "vop unbend" command can straighten the density from a curved filament. The tricky part is you need a center-line path. You could define that by placing 3 markers. http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/vop.html#unbend Bending a straight filament map to follow a curve is not available. This is the inverse of the vop unbend operation. You can do a fake version of this by making a helical curve around the straight density map and then using vop unbend. That won't center the filament on the helical path, but it will create something similar by straightening the helix which I think will twist the density. Tom > Hi Tom- > > Not sure if there is an option for this - > > 1. can we bend an electron density map of a filament/tube by an small > angle ? > > 2. Can we generate an helix using a density map of a filament. > An example will be that of the bacterial flagella - it has the local > helical arrangement of the flagellin molecule and then thAT > ARRANGEMENT extends over several micrometers like a "sinusoidal wave" > with different helical arrangement to the filament axis. > > Best > Andy > > > Dr. Anindito Sen (Ph.D) >

1 0