- Chimera-users - RBVI Mailing Lists

About volume viewer in chimera 1.6.1
by Yuhong Li 30 Apr '12

30 Apr '12

Dear Chimera developers, Hi! I am a student from Sun Yat-Sen University. I recently began using the Chimera for my study, but I have a problem about how to read the displayed information in volume viewer. Here I attach a picture and inside the picture there are two places I've point out with arrows which represent two different map files in mrc fornat. Now I was wondering if you could spare some time to explain it to me that what do those different numbers mean? Hope to hear from your reply in very soon. Best Wishes, Yuhong Lee

2 2

addaa command
by Mahendra B Thapa 27 Apr '12

27 Apr '12

Dear chimera users, I also added an amino acid residue at the N-terminus of a small protein following the posting: [Chimera-users] Using addaa at the N-terminus?*[ posted date:Fri Sep 30 10:36:31 PDT 2011]* The residue is added at the n-terminus but the added residue is not connected to rest of the protein structure in the graphical display [i.e, a little gap between the added residue and the rest of the protein at N-terminus]. Is this due to wrong command ? Thanks, Mahendra Thapa Graduate student Department of Physics University of Cincinnati, OH

2 1

Per model clipping not working ?
by Damien Larivière 27 Apr '12

27 Apr '12

Dear all, I try to use the Per-Model clipping tool but this error is returned, whatever the PDB file: Model cannot currently be clipped [no valid bounding box] Do I make something wrong? I use Chimera version 1.6.1 (build 35849). My configuration is Windows 7 64bits, ATI Mobility Radeon HD5650. Many thanks for your help! Damien

2 1

loading trajectories with no box line
by Michael F. Bruist 25 Apr '12

25 Apr '12

Is there a way to load trajectories that do not have box lines? I have trajectories from which most solvent and ions have been stripped. They were written with ptraj as amber coordinate files with the "nobox" option. I cannot get these to stay in register for MD movie. I do not have this problem for trajectories for which the box line has not been removed. Mike Bruist Dr. Michael F. Bruist Associate Professor of Biochemistry Department of Chemistry & Biochemistry Griffith Hall 147 University of the Sciences in Philadelphia Philadelphia, PA 19104-4495 Phone (Office & Voice Mail): (215) 596-8530 Department Phone: -8839 Department Fax: -8543 E-mail: m.bruist @usp.edu web page: http://tonga.usp.edu/bruist

2 1

saving transformed map
by Elaine Meng 25 Apr '12

25 Apr '12

On Apr 25, 2012, at 9:04 AM, Ariane Briegel wrote: > Dear Elaine, > I have been using chimera to align several EM maps in respect to each other. > These maps are generated by imod (mrc format). I know I can save the transformed maps in chimera when I use the 'save session' command, but that only allows me to use the transformed maps inside chimera. I would like to save the individual transformed maps- either in mrc or CCP4 format. > I tried to find any info about that in the manual, then checked the old chimera list comments. I found a similar question from the year 2006 (http://plato.cgl.ucsf.edu/pipermail/chimera-users/2006-July/000864.html) but at that time your conclusion was it is not possible. > Do you know if there is a way to do this now? > Thank you so much for your answer, > Ariane Hi Ariane, The map formats don't include transformation information, so there is no way to save EXACTLY the same data, except transformed. However, you can resample the map onto a new grid, and then save the new map. The values at the new grid points are determined by interpolation, so there will be some (generally minor) data loss. <http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/savemodel.html#afterfitting> The resampling of one map onto the grid of another can be done with command "vop resample" <http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/vop.html#resample> For example: vop #1 resample onGrid #0 To save the resampling result to a file, you could use "File... Save Map as" in the Volume Viewer dialog or the command "volume" with its "save" option. 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

2 1

Re: [Chimera-users] Chimera-users Digest, Vol 108, Issue 25
by Boaz Shaanan 24 Apr '12

24 Apr '12

Hi, Elaine is correct in that newly deposited structures (since a year or so?) have to have their TLS information converted to aniso B and have those values in the appropriate place on the ATOM card. However, this was not the case for the older entries. The 1ss8 entry seems to have been refined by refmac5 from the ccp4 package. In this package there is a program called tlsanl which will do the job of converting the TLS infotmation to B-aniso. (a several seconds job) So maybe Steve could use it and thus solve his problem? Cheers, Boaz 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 ________________________________________ From: chimera-users-bounces(a)cgl.ucsf.edu [chimera-users-bounces(a)cgl.ucsf.edu] on behalf of chimera-users-request(a)cgl.ucsf.edu [chimera-users-request(a)cgl.ucsf.edu] Sent: Tuesday, April 24, 2012 7:19 PM To: chimera-users(a)cgl.ucsf.edu Subject: Chimera-users Digest, Vol 108, Issue 25 Send Chimera-users mailing list submissions to chimera-users(a)cgl.ucsf.edu To subscribe or unsubscribe via the World Wide Web, visit http://plato.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: remote session (Greg Couch) 2. Rendering TLS information in Chimera ? (Steven Ludtke) 3. Incomplete chain contain many dotted lines while displaying (Junfeng Huang) ---------------------------------------------------------------------- Message: 1 Date: Mon, 23 Apr 2012 12:05:41 -0700 From: Greg Couch <gregc(a)cgl.ucsf.edu> To: Hernando J Sosa <hernando.sosa(a)einstein.yu.edu> Cc: "chimera-users(a)cgl.ucsf.edu" <chimera-users(a)cgl.ucsf.edu> Subject: Re: [Chimera-users] remote session Message-ID: <4F95A805.1080609(a)cgl.ucsf.edu> Content-Type: text/plain; charset="iso-8859-1"; Format="flowed" The error is that the OpenGL library is looking for a NVidia-specific extention to the remote protocol and failing when it is not there. That is a driver bug -- it shouldn't matter what graphics driver you have on either of the systems. There are two things I'd do to try and fix it: (1) try updating the NVidia driver on the remote computer that is running chimera (not the computer displaying chimera), and (2) try uninstalling the NVidia driver on the remote computer, and reverting to the software implementation of OpenGL (mesa). The first option is preferable if it works, because the proprietary (NVidia/AMD) drivers have been faster and less buggy than the mesa driver for most uses. Good luck, Greg On 04/20/2012 02:25 PM, Hernando J Sosa wrote: > > Hi > > I used to be able to run chimera remotely from a machine with an X- > windows server (the one that comes with Cygwin in this case). But now > I get the following error: > > Xlib: extension "NV-GLX" missing on display "localhost:10.0". > > X Error of failed request: BadRequest (invalid request code or no > such operation) > > Major opcode of failed request: 148 (GLX) > > Minor opcode of failed request: 181 () > > Serial number of failed request: 12779 > > Current serial number in output stream: 12784 > > I t does work fine if run directly from the computer where chimera is > installed. Any idea what the problem could be? > > Thanks > > Server computer: > > Chimera v 1.5 (similar problem with 1.6) > > OS: Linux OpenSuse 11.3 64b > > Client computer: > > OS: Windows 7 64. > > Graphics: ATI Radeon HD 5450 >

1 0

electron density ven diagram
by Terry Lang 24 Apr '12

24 Apr '12

Hey Guys, I am interested in making kind of a ven diagram between two maps for electron density (see attached). What I would like to do is where the density is above a certain threshold for both maps, color the map grey, and where it is above the threshold for map A only, color red, and for map B only, blue. Plan #1 --superimpose maps (have this working) --create a third map (have this working) --do a direct < > = point-by-point comparison between map A and B and reset the values of the third map to be -1 (map A only), 0 (both maps) or 1 (map B only) --color the map using the histograms PROBLEM: The maps are not gridded the same. After some careful thought, I don't think I actually can get them gridded the same for technical reasons. Plan #2 --superimpose maps (have this working) --create copy of map A (coarser map) with map A gridding (have this working) --use interpolation function to interpolate map B (finer map) onto map A gridding --create copy of map A and do a direct < > = point-by-point comparison between map A and interpolated map B and reset the values of the fourth map to be -1 (map A only), 0 (both maps) or 1 (map B only) --color the map using the histograms QUESTION: Is there an easier way to do this with tools that already exist in Chimera? If not, can I get a little bit of info on how the data is stored in the ccp4 style maps (eg start at origin and loop over x, loop over y, loop over z or something else?) Let me know if you need more details or clarification on this rats nest. Thanks! Terry -- Paula Therese Lang Postdoctoral Scholar Alber Lab UC Berkeley/QB3

2 2

Incomplete chain contain many dotted lines while displaying
by Junfeng Huang 24 Apr '12

24 Apr '12

Hello Everyone， I have a problem while display PDB file. I extract the residues from protein-protein interface by ligplus into a new PDB file, and want to display and analysis it. As the file contain only the residues form interface, they are discontinuous. While display with pymol, it is OK. However, While using Chimera, there are many dotted lines between the residues. I know what the dotted lines are, but I do want to hide all these lines, can somebody please give me some hints? I sent two figures in the attachment. It's the first time I used mail-list, if it doesn't support figures, please check here: Interface display by pymol http://bbs.sysu.edu.cn/attach/Pictures/1335235173.JPG Interface display by Chimera http://bbs.sysu.edu.cn/attach/Pictures/1335235254.JPG Thank you very much!

2 2

Charge data
by Nikolay Igorovich Rodionov 24 Apr '12

24 Apr '12

Hi everyone, I know chimera does not produce PQR files, but I was wondering if it is possible to create a data text file with the atoms' three dimensional coordinates and charges after applying the amber forcefield? Sincerely, Nikolay Rodionov

3 5

Rendering TLS information in Chimera ?
by Steven Ludtke 24 Apr '12

24 Apr '12

Hi. A question about rendering B-factor information: Over the last decade, rather than simple anisotropic B-factors, it has become increasingly common to represent local "motion" using the translation-libration-screw (TLS) model, which provides for anisotropic representation of such motions. Is there a good way to represent this data in Chimera ? Here is an example. If you render the PDB structure 1OEL as atoms, and color the atoms by B-factor, you can immediately observe regions of high variability in the apical domain. This same data was later re-refined using a TLS model as PDB structure 1SS8. When you render this structure in Chimera in the same way, you see high B-factors only for a few specific sidechains. The overall large-scale variability is gone, presumably because the TLS parameters have absorbed it, and the "B-factor" is now only a residual B-factor. However, I can't see any obvious way of rendering the TLS information. Any tips ? ---------------------------------------------------------------------------- Steven Ludtke, Ph.D. Associate Professor Co-Director National Center For Macromolecular Imaging (ncmi.bcm.edu) Co-Director CIBR Center (www.bcm.edu/research/cibr) Dept of Biochemistry and Mol. Biol. (www.bcm.edu/biochem) Baylor College of Medicine sludtke(a)bcm.edu stevel(a)alumni.caltech.edu

3 2

remote session
by Hernando J Sosa 23 Apr '12

23 Apr '12

Hi I used to be able to run chimera remotely from a machine with an X- windows server (the one that comes with Cygwin in this case). But now I get the following error: Xlib: extension "NV-GLX" missing on display "localhost:10.0". X Error of failed request: BadRequest (invalid request code or no such operation) Major opcode of failed request: 148 (GLX) Minor opcode of failed request: 181 () Serial number of failed request: 12779 Current serial number in output stream: 12784 I t does work fine if run directly from the computer where chimera is installed. Any idea what the problem could be? Thanks Server computer: Chimera v 1.5 (similar problem with 1.6) OS: Linux OpenSuse 11.3 64b Client computer: OS: Windows 7 64. Graphics: ATI Radeon HD 5450

2 1

Displaying atoms within a distance
by Nikolay Igorovich Rodionov 21 Apr '12

21 Apr '12

Hi everyone, I am trying to display waters within a set distance from a solvated protein. Can someone please tell me how I can use the display or show commands to selectively show only the residues that are within a specific distance from a residue type? I've looked at the chimera manual and I am still confused on how these commands function. Thank you, Nikolay Rodionov

3 3

Charge calculations PDB2PQR
by Nikolay Igorovich Rodionov 20 Apr '12

20 Apr '12

Hi all, Does anyone knows any alternative charge calculating software packages to PDB2PQR. I need to calculate partial charge distribution but the local copy of PDB2PQR that I installed has some kind of glitch where it is scrapping all of the water molecules from my solvated protein while doing the analysis. I think it has something to do with the side of my PDB file since I am using about 30K atoms, at least, and everything goes perfectly when I work with less 10K atoms. So, I guess the second question is, has anyone else ever encountered this problem? Sincerely, Nikolay Rodionov

3 7

Volume Viewer Mesh
by Terry Lang 19 Apr '12

19 Apr '12

Hey Guys, I am using the Volume Viewer to look at some electron density for a figure. I'd like to use the mesh option but when I generate the figure, the mesh is too fine, even at the starting default values. Is there a way to make the mesh coarser? Thanks, Terry -- Paula Therese Lang Postdoctoral Scholar Alber Lab UC Berkeley/QB3

3 4

modeller interface
by Sadler, Evan 19 Apr '12

19 Apr '12

Hi, I'm trying to construct what I think should be a simple homology model and it's failing as below whether I use the web server or run modeller locally: Target seq: VWF_CK Template structures: 2TGI.pdb (#0)Run locally, the Modeller binary location: /Library/modeller-9.10/bin/mod9.10 Now, modeller is running locally... <type 'exceptions.UnboundLocalError'> Exception in Tk callback Function: <function command at 0x17982770> (type: <type 'function'>) Args: () Traceback (innermost last): File "/Applications/Chimera.app/Contents/Resources/lib/python2.7/site-packages/Pmw/Pmw_1_3/lib/PmwBase.py", line 1747, in __call__ return apply(self.func, args) File "/Applications/Chimera.app/Contents/Resources/share/chimera/baseDialog.py", line 449, in command getattr(s, buttonFuncName(txt))() File "/Applications/Chimera.app/Contents/Resources/share/chimera/baseDialog.py", line 700, in OK self.Apply() File "/Applications/Chimera.app/Contents/Resources/share/MultAlignViewer/ModellerDialog.py", line 579, in Apply veryFast=self.veryFastVar.get(), tempPath=self.pathTempOpt.get().strip(), **kw) File "/Applications/Chimera.app/Contents/Resources/share/MultAlignViewer/ModellerBase.py", line 247, in model if loopData: <type 'exceptions.UnboundLocalError'>: local variable 'loopData' referenced before assignment UnboundLocalError: local variable 'loopData' referenced before assignment File "/Applications/Chimera.app/Contents/Resources/share/MultAlignViewer/ModellerBase.py", line 247, in model if loopData: See reply log for Python traceback. My version of Chimera is alpha version 1.7 build 36205, running under Mac 10.6.8. Do you have any advice? Thanks, Evan Sadler -- J. Evan Sadler, M.D., Ph.D. Chief, Division of Hematology Department of Medicine Washington University 660 S. Euclid Ave., Box 8125 St. Louis, MO 63110 Tel: 314-362-9067 FAX: 314-454-3012 email: esadler(a)wustl.edu

2 1

How to save high resolution image?
by Yogendra 19 Apr '12

19 Apr '12

Dear Sir, I am a new user of Chimera software. I am using this software for academic purposes (i.e. research publications). I am trying to save the image with more then 600 dpi but every time i am ending up with an image of 96 dpi. I am using the POV-Ray method for image improvement but everything is going to the drain when taking the final print of the document. Can you please help me ho can i save a high resolution image? I am briefly describing my protocol of image generation: 1) I am making a density map of particular PDB file. 2) loading in the Chimera and adjusting the volume by volume viewer option 3) Coloring it using described colors 4) Using Pov-Ray for increasing effects on plane 5) Increasing the pixels to about 5000 in the pixel option 6) finally trying to save an image with 600 or 1200 dpi. But at the end of this, i am just getting an image of 96 dpi every time. Can you please help me out to make an image of higher dpi (600, 1200 or more if possible). It will be very much helpful for me if can do it. I am stuck at this point.I am badly waiting for your reply whether i can do it or should i look for some another software. Regards Yogendra S Rathore Research Scholar ______________________________________________________________________ सूक्ष्मजीव प्रौद्योगिकी संस्थान (वैज्ञानिक औद्योगिक अनुसंधान परिषद) Institute of Microbial Technology (A CONSTITUENT ESTABLISHMENT OF CSIR) सैक्टर 39 ए, चण्डीगढ़ / Sector 39-A, Chandigarh पिन कोड/PIN CODE :160036 दूरभाष/EPABX :0172 6665 201-202

4 6

Methanesulfinamide Stability and Synthesis
by Nancy 19 Apr '12

19 Apr '12

Hello All, I would greatly appreciate it if anyone has information on the stability of methanesulfinamide (SMILES: C[S](=O)N([H])[H]), as well as it's stability in water. Thank you.

3 2

Command line "Render by attribute"
by rswett＠chem.wayne.edu 18 Apr '12

18 Apr '12

So, I"m wondering if there's a way to call "Render by Attribute" command line for movie making purposes. I have several Defined Attributes and i'd like to be able to recolor surfaces by the various attributes at different points in a movie. What would the format for that type of command be? Thanks for any help you guys have! ~Rebecca

2 1

Check Marks in 2D Labels
by Paul Cernek 17 Apr '12

17 Apr '12

Hello, I'm having trouble inserting a check mark (Unicode page here: http://www.fileformat.info/info/unicode/char/2713/index.htm) into a 2D label, using the Chimera Daily Build from April 10 2012 on Windows Vista. I know that you can reference a special character via its Unicode escape sequence (in this case \u2713). Indeed, when I specify such text from the command line, or when I copy/paste the check mark into the text box of the 2D Labels window, it displays exactly what I want IN THE 2D LABELS WINDOW. However, in the main Chimera model display window the check mark is NOT properly depicted in the label. Instead, I get that empty square you get whenever a program encounters an unrecognized character. Other more conventional symbols (such as the Greek letter 'alpha') display normally. Why would the "2D Labels" window display the check mark properly, while the main display window does not? Any assistance in this matter is very much appreciated. Thanks, -Paul

3 2

APBS maps
by Rebecca Swett 12 Apr '12

12 Apr '12

Hey guys, i'm working with some electrostatic potential maps and am interested in doing a difference map. Two nearly homologous proteins with nearly identical surface shapes, and radically different electrostatics. I"m thinking a difference map would illustrate the differences well, but in looking at the APBS output I have no idea how to do that manually. Does anyone here have an idea about how to do that in chimera?

4 7

Re: [Chimera-users] [Dock-fans] Help with antechamber
by Trent E. Balius 10 Apr '12

10 Apr '12

Hi Brad, I am cc'ing this to the chimera list. I think that you should write a python script to run DOCK prep on the command line. See the following: http://www.cgl.ucsf.edu/pipermail/chimera-dev/2011/000784.html For other example scripts see the following: http://www.cgl.ucsf.edu/pipermail/chimera-users/2009-February/003539.html You should not have to specify the AMBERHOME path. This is done internally to chimera. as an aside, It is unnecessary to run antichamber with sudo. I hope this helps, Trent E. Balius Graduate Student, Rizzo Group, Department of Applied Mathematics and Statistics, Stony Brook University. Office: Math Tower 3-129, Phone: (631) 632-8519 URL: http://www.ams.sunysb.edu/~tbalius/ -----Original Message----- From: Brad Ridder <clustro(a)gmail.com> To: dock-fans <dock-fans(a)docking.org> Sent: Wed, Feb 29, 2012 4:38 pm Subject: [Dock-fans] Help with antechamber Hi Dockfans, My name is Brad; I am a grad student in the dept. of chemical engineering at Purdue University. We are experimenting with DOCK for various drug discovery work. At the moment, I am trying to do execute all of the "dock prep" commands through the Terminal, since it is not practical to do large volumes of Dock Prep by hand through the CHIMERA GUI. Corresponding to what is in the DOCK tutorial, I am on "Receptor and Ligand Preparation." I got reduce to work, in order to add hydrogens. Now, I am trying to get antechamber to work, in order to add the charges, but it is not working. Here is the command I am running: sudo ./antechamber -i $DESKTOP/3GLR.pdb -fi pdb -o $DESKTOP/3GLRcharged.mol2 -fo mol2 Which returns the error: Warning: $AMBERHOME and $DOCK_HOME enviornment strings are not set, use "mopac.sh" in the work directory Cannot open CONNECT.TPL , exit The atom number exceeds the MAXATOM, reallocate memory I looked around in the dock6 directory, and found a folder called "parameters/antechamber", which contains CONNECT.TPL. I figured this is what should be the $AMBERHOME directory. For DOCK_HOME, I figured it should be the dock6 folder. Here is what I have in my .bash_profile file in my home folder. export AMBERHOME='/home/.../DOCK6.5/dock6/parameters/antechamber' export DOCK_HOME='/home/.../DOCK6.5/dock6' I also tried using the "amber11" folder that comes with CHIMERA as AMBERHOME, but to no effect. I don't actually have AMBER installed on my computer, but since CHIMERA can add charges without actually requiring the purchase of AMBER, I do not think that is the problem. I did not know what to do with "mopac.sh." I ran the script, but it says it needs input files. What am I doing wrong that is preventing antechamber from executing correctly? -- Bradley James Ridder Chakrabarti Group Graduate Student School of Chemical Engineering Purdue University _______________________________________________ Dock-fans mailing list Dock-fans(a)docking.org http://mailman.docking.org/mailman/listinfo/dock-fans

4 7

vector of axis
by Gwo-Yu Chuang 09 Apr '12

09 Apr '12

Dear all, I am trying to figure of the principal inertia axis of a protein domain using Chimera. Is there any way to get the vector of the axis (the values), besides just showing it on the screen? Many thanks Yours Sincerely, Paul

2 1

Morph with Cofactors
by Marco Jost 06 Apr '12

06 Apr '12

Hi, I am trying to morph two conformations of a protein that includes multiple cofactors. While the protein structure is morphed just fine, the cofactors are not included in the morph (or in the output PDB). I numbered the cofactors identically in both files and everything else should be consistent as well. Does anybody know how to include cofactors in the morphing and have them be written into the output PDB? Thanks! Marco Jost Ph. D. Student, Department of Chemistry Massachusetts Institute of Technology

2 1

Re: [Chimera-users] [Open Babel] UCSF Dock mol2 files
by Scott Brozell 06 Apr '12

06 Apr '12

Hi, This is a good suggestion. But before implementation, 1) do the DOCK developers or users have any comments ? 2) do the Chimera developers or users have any comments ? In particular, does Chimera use ########## as the lexeme for the DOCK token ? It might be useful to have a consistent lexeme across programs. Or maybe the Babel folks have comments on what would be a more general tag ? thanks, scott On Tue, Apr 03, 2012 at 05:40:46PM +0200, Maciek Wójcikowski wrote: > Here you have two example files > http://dl.dropbox.com/u/6846289/1afk_ligand_conformers.mol2 > http://dl.dropbox.com/u/6846289/3h78_ligand_conformers.mol2 > > What OB should do (in my mind) is to take "Name:" value as target molecule > name and pass all 3 grid scores (total, vdw and ele) somehow (i guess > comment is a right place, since there is not much choice), so that it can > output it in simmilar manner. > > Maybe it would be easier to have an option "d" (aka. dock) which would copy > ALL the ###### data, and some other one, such as "s" (aka. sufix) to take > suffixed name into account. Disadvantage of above is that it wouldn't work > for other formats. > ---- > Pozdrawiam, | Best regards, > Maciek Wójcikowski > maciek(a)wojcikowski.pl > > > W dniu 3 kwietnia 2012 17:23 użytkownik Geoffrey Hutchison < > geoff(a)geoffhutchison.net> napisał: > > > > I have a thought about mol2 files form Dock program, which contain > > > additional information in commentary (###). Those information is being > > lost > > > during openbabel operations, which is a pity. I'd like to propose a > > in/out > > > option for mol2 for Dock'ed mol2's to read/write compatible information > > from > > > those files. Especially conformer ID is important for further anylisis. > > > Simmilar approach is used with pdbqt's contained in Autodock's .dlg > > files. > > > > > > If you can give us a selection of example files and a list of what > > information you'd like retained, it's fairly easy to patch the parsers. > > > > Cheers, > > -Geoff

3 4

Selecting submodels when they are density maps
by Javier Velazquez 05 Apr '12

05 Apr '12

Hi, I am interested in generating a map from a PDB with 2 chains, and color each chain differently. I can do it very easily using the graphical interface, but I couldn't figure out how to do it from the command line. Here are my commands: molmap #0:.A 7 molmap #0:.B 7 And now adding color is simple: volume #0 color red But this adds color to both maps. Is there a way to specify a submap, so I can add color separately. I am looking for something like this: volume #0:.A color red which I tried and didn't work. Thank you for your help. Javier

2 1

Merge the protein sequences of two different chains
by Hsih-Te Yang 05 Apr '12

05 Apr '12

Dear all: I am a new beginner of Chemira software package. By using "Multalign Viewer", I would like to merge the protein sequences of two different chains (200~300 a.a.). After that, I want to map some epitopes (8 ~ 20 a.a.) to this merged sequence. How can I do? Thanks for any kind help further. Hsih-Te

2 4

Solvate models
by Nikolay Igorovich Rodionov 05 Apr '12

05 Apr '12

Hi all, I was wondering if some could explain to me the differences between the water models available under the solvate option. * POL3BOX - POL3 water model * SPCBOX - SPC/E water model * SPCFWBOX - SPC/Fw water model * QSPCFWBOX - qSPC/Fw water model * TIP3PBOX - TIP3P water model * TIP3PFBOX - TIP3P/F water model * TIP4PBOX - TIP4P water model * TIP4PEWBOX - TIP4P/Ew water model Thank you all in advance, Nikolay Rodionov

2 1

Coloring distance labels
by George Tzotzos 04 Apr '12

04 Apr '12

I'm somewhat lost in the manual. I'm measured some specific distances which are clearly shown. Changing to Presets: Publication 3 the distances are shown in yellow in white background. I'd like to change the color of the distances. I'm pretty sure there's a quick solution for this. I simply can't find it. Any help would be much appreciated. Regards George

2 2

Coloring distance labels
by George Tzotzos 04 Apr '12

04 Apr '12

I'm somewhat lost in the manual. I'm measured some specific distances which are clearly shown. Changing to Presets: Publication 3 the distances are shown in yellow in white background. I'd like to change the color of the distances. I'm pretty sure there's a quick solution for this. I simply can't find it. Any help would be much appreciated. Regards George

1 0

color by RMSD
by Thomas Evangelidis 04 Apr '12

04 Apr '12

Dear UCSF Chimera users and developers, Is it possible to color a pair of structures representing different stages of the thermodynamic cycle of a protein by the RMSD, without aligning them? thanks in advance, Thomas

3 6

Re: [Chimera-users] [Chimera] #10638: heme not included in morph trajectory
by Elaine Meng 04 Apr '12

04 Apr '12

Hi Francesco, There is some progress... in the newest daily build (Apr 4) I can run the commands below to successfully show a morph of 1hho and 2hhb that includes heme and iron. 1hho only has two chains, so it is necessary to delete C,D from 2hhb before morphing. To include heme and iron, it is necessary to delete all the other het residues that may be different between the structures. I hope this helps, Elaine Here between "----" are contents of Chimera command file, which you could save as plain text file named something.cmd and open in Chimera to execute the commands. Lines starting with # are just comments/explanation and don't do anything. ------------- open 2hhb open 1hho # # delete 2hhb chains C,D # delete #0:.c-d # # delete het residues that may be different between the structures so that # the HEM residues will be included in the morph (if you want to see those # residues, open another copy of structure not used for morphing) # del :hoh,po4,oxy # # superimpose the structures; I used "mm" with defaults but it could be # done with options, or MatchMaker GUI, or "match" command, # or manually with mouse # mm #0 #1 # # create morph trajectory which will be model #2 # morph start #0 frames 25 morph interpolate #1 frames 25 morph movie ~modeldisp #0,1 --------------------- Then I tried to get fancy and make a tetramer of 1hho so that the morph could be from tetramer to tetramer. However, then only the heme/iron from A and B were included, not from C and D. I don't know why, because they are all in the input structures!! Here is that script in case it might be useful: --------------------- open 2hhb open 1hho # # make 1hho tetramer; #1 is original, #2 additional copy # sym #1 # # combine #1,2 dimers into tetramer model #3, close #1,2 # combine #1,2 name 1hho-tetra newchainids true close true # # delete het residues that may be different between the structures so that # the HEM residues will be included in the morph (if you want to see those # residues, open another copy of structure not used for morphing) # del :hoh,po4,oxy # # superimpose #0 and #3 # mm #0 #3 # # create morph trajectory which will be model #1 # morph start #0 frames 25 morph interpolate #3 frames 25 morph movie ~modeldisp #0,3 ---------------------

1 0

occupancy field
by Darjan Andrejc 03 Apr '12

03 Apr '12

Hey, Is it possible to change occupancy field to 1.0 in Chimera? I need a pdb file of my protein for docking and the file needs to have an occupancy field 1.0. Thanks for the help Darjan Andrejc, Mpharm PhD student, Department of Biotechnology National institute of Chemistry, Ljubljana, Slovenia

2 2

Minimize with fixed atoms
by Hurt, Darrell (NIH/NIAID) [E] 02 Apr '12

02 Apr '12

Hi there, I'm trying to minimize a protein after mutating a residue with the Rotamers panel. I first do a Check Clashes and select the clashing atoms. Then I expand the selection using Zone to <5 A. I use these selected atoms as the non-fixed atoms and the unselected atoms as the fixed ones. When I then run the minimization with the appropriate settings, I get an RMSD of 0.000 A for each step, even though the step size is 0.02 A. Nothing moves. Same thing if I have no fixed atoms. Am I in a local minima? Is that why nothing moves? The energy is very high (from the Reply Log): Potential energy: 222274199.525130, Gradient norm: 833984762.552943 Updated 6029 atoms. RMSD: 0.000000 Finished 1 of 5 steepest descent minimization steps Any ideas? I do have some hetero atoms... Thanks, Darrell Darrell Hurt, Ph.D. Section Head, Computational Biology Bioinformatics and Computational Biosciences Branch (BCBB) OCICB/OSMO/OD/NIAID/NIH 31 Center Drive, Room 3B62B, MSC 2135 Bethesda, MD 20892-2135 Office 301-402-0095 Mobile 301-758-3559 http://bioinformatics.niaid.nih.gov (Within NIH) http://exon.niaid.nih.gov (Public) Disclaimer: The information in this e-mail and any of its attachments is confidential and may contain sensitive information. It should not be used by anyone who is not the original intended recipient. If you have received this e-mail in error please inform the sender and delete it from your mailbox or any other storage devices. National Institute of Allergy and Infectious Diseases shall not accept liability for any statements made that are sender's own and not expressly made on behalf of the NIAID by one of its representatives.

3 5

Re: [Chimera-users] morphing
by Elaine Meng 31 Mar '12

31 Mar '12

Hello Francesco! First make sure you are hiding the input conformations, so that you are seeing only the morph (for example, you could show Model Panel from Favorites menu, then uncheck the S(hown) boxes for everything except the morph results). I say that because if you see something not moving, it might be the original structure. As explained here, the morph should contain all of the atoms in the protein chain that are included in all of the input structures. <http://www.cgl.ucsf.edu/chimera/docs/ContributedSoftware/morph/morph.html#p…> Since your input structures are the same protein, if each one includes the sidechain atoms, the morph should contain all the same sidechains and they should move from one position to the other during the morph. They may not be displayed automatically, however. Try displaying all atoms (again better to hide the other structures first) to make sure of what is there. About the ion: even if the ion is in both structures, it may not be in the morph if you are using an older version of Chimera. However, there have been some fixes, and I believe if you are using Chimera 1.6 or newer, the ion should be included. The fix I'm thinking of was 6 months ago and is described here: <http://plato.cgl.ucsf.edu/trac/chimera/ticket/9886> If you try Chimera 1.6 and the suggestions above but you believe there is a bug, you could use Help... Report a Bug in the menu and describe how to reproduce the problem, attaching session file if necessary. 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 Mar 31, 2012, at 5:24 AM, Francesco Malatesta wrote: > Hello from Rome! > I'm trying to morph 2 different conformations of the same protein. I get a very very nice movie, but I can't get amino acid side chains or heteroatoms (heme iron) to move also. > Is there any way to do this? > thank you for help > Francesco

2 1

How to save VIPER transformation matrix
by snoze pa 30 Mar '12

30 Mar '12

Dear Chimera Users, Is it possible to save the coordinates after the VIPER transformation (xyz-2 fold symmetry) in Chimera. I can create a multiscale biological molecule in chimera but it is not saving the viper matrix in pdb file. Please note that this pdb is a model not the crystal structure and I am trying to save the viper matrix for this pdb. Thank you, S

4 5

Chimera is not loading properly
by Andrea Edwards 27 Mar '12

27 Mar '12

I have not used Chimera in awhile. I just downloaded the new version for Mac-OSX-64-bit and when I open it, none of the buttons on the gui are visible. There are ovals and things where there should be buttons, but not words on them. How can I fix this? attached is a screenshot. Thanks, Andrea

3 2

PDB model merging
by Nikolay Igorovich Rodionov 27 Mar '12

27 Mar '12

Hi all, I'm working on a project that requires me to merge models created through PDB BIOMT records. I have been using the copy&combine feature in the model panel to merge all of my models, and then writing a new PDB of only the combined model. I've noticed that this method deletes a large amount of information about the structure, including secondary structure records. Is there a way to merge models without losing such data? If it is relevant, I am merging 90+ models.

3 5

Trajectory analysis: hbond
by George Tzotzos 26 Mar '12

26 Mar '12

I'm wondering if there's a way to determine hbond formation during a trajectory in a similar way to the way the Plot distances tool in MD Movie/Analysis works. Or dumb values to a file. Many thanks in advance for any suggestions Regards George

3 4

Re: [Chimera-users] Electrostatic Potential using APBS .dx file
by Tom Goddard 26 Mar '12

26 Mar '12

Hi Soca, Sounds like you are doing everything right. That makes me suspect that something is wrong with your APBS file -- like all the potential values are zero. To see what the distribution of the potential values are, open the APBS file and then use menu Tools / Volume Data / Volume Viewer, and click on the histogram (empty rectangle in the middle of the volume dialog). That will show a histogram of the potential values. You can move the vertical bar on the histogram to see what the potential value is at that position. Check if -5 to 5 is a reasonable range. Tom > Hello, > > I have been attempting to generate a electrostatic surface map in > chimera using an external potential file generated from APBS (.dx). > > I understand that Chimera has an internal surface calculation protocol > via Coulombic Surface Coloring, however, I would like to use the > Electrostatic Surface Coloring. > > However, when I open my molecule (PDB), show surfaces, and choose its > MSMS main surface to color by electrostatic potential using the .dx > file, the surface map is not the red, white, blue that I have set in > my parameterts under Surface color. Everything turns to grey. My other > parameters include: > > red -5, white 0, blue 5. Colors 3. Palette Red-Blue. Surface offset 1.4. > > I feel as though I am missing something trivia that would allow me to > generate the surface image I want, however I am lost after may > searches for advice through the forums and tutorials. > > Please let me know what I am missing. > > Thanks > > Soca > > -- > Ardian Soca Wibowo, Ph.D. > Postdoctoral Associate > Dann Research Group > Department of Chemistry > Indiana University

1 1

Dimethyl Sulfoxide Electromagnetic Absorbance
by Nancy 25 Mar '12

25 Mar '12

Hi All, I would greatly appreciate it if anyone could provide a reference or information on the electromagnetic absorbance characteristics of dimethyl sulfoxide, especially in the region of 800-1400nm. Thank you, Nancy

1 1

RMSD plot
by lecan 23 Mar '12

23 Mar '12

Hi everybody. I am interested in plotting RMSD values vs time from molecular dynamics. Is there a way to extract the RMSD values from a the RMSD map in the analysis tool of MD movie. Or maybe exists another way to do this with Chimera. Thank you in advance Luis L.

2 1

how to measure the volume of the active site
by Thomas Evangelidis 23 Mar '12

23 Mar '12

Dear UCSF Chimera developers, I am interested in measuring the volume of the active site of a kinase. I have searched through the documentation and the only thing I found was how to measure the enclosed volume of a surface object. Is it possible to do what I want with UCSF Chimera? Could you please point me towards the relevant commands? Thanks in advance for any help. Thomas -- ====================================================================== Thomas Evangelidis PhD student Biomedical Research Foundation, Academy of Athens 4 Soranou Ephessiou , 115 27 Athens, Greece email: tevang(a)bioacademy.gr tevang3(a)gmail.com website: https://sites.google.com/site/thomasevangelidishomepage/

2 1

Multi-PDB fitting question
by Tom Goddard 21 Mar '12

21 Mar '12

> Dear Tom, > > ... > Currently we are trying to use a lot of PDBs (hundreds!) to fit one > density map and trying to find out the best fitted PDB, and we found > that fitting the PDBs one by one manually was very time-consuming and > low-efficiency. So we are wondering that if there is any advanced way > or batch file that can speedup this process? > > Thank you so much and we are looking forward your reply. > > Best Regards, > Lei > Hi Lei, Ok, I added a new "eachModel" option to the Chimera fitmap command so that one command can independently fit many molecules in the same map. For example, fit #1-50 #0 resolution 13.5 eachModel true listFits true will fit PDB models with id numbers 1 to 50 into map with id #0 and list the fits in a dialog. Each of the fits is done independently (no account for clashes between molecules). Details on each fit are given in the Reply Log dialog (menu Favorites / Reply Log). This will be in tonight's Chimera daily builds (run at midnight) if they succeed: http://www.cgl.ucsf.edu/chimera/download.html#daily Here is the fitmap command documentation which will be updated later this week or next week. http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/fitmap.html You will probably find that you can't open hundreds of PDB molecules in Chimera because it takes too much memory. But you can fit them in a few batches if this happens. Tom

1 0

disulfides: naive question
by George Tzotzos 21 Mar '12

21 Mar '12

Apologies for the earlier message. I need to clarify. Loading the topology and trajectory files onto Chimera and selecting functional group/disulfides shows the bonds throughout the trajectory. The problem I have is that I don't see the bonds in pdb files that I've generated from the topology files. Thanks again George

2 4

disulfides: naive question
by George Tzotzos 20 Mar '12

20 Mar '12

I've built a number of disulfide bonds in Amber topology files. From these, I've generated pdb files. When I load them in Chimera and do Select/Chemistry/Functional group/disulfed, I fail to see the selected residues and/or the bonds between them. Obviously, I'm doing something wrong. Any help would be greatly appreciated Regards George

1 0

depicting secondary structure elements
by Jordi Querol Audi 20 Mar '12

20 Mar '12

Dear Chimera developers, I would like to show a fitting of a pdb into my EM reconstruction but instead of displaying helices and b-sheets I would prefer to use cylinders or tubes. Is there any way to display a pdb this way? Thank you so much in advance, Jordi

3 2

smoothMD.py
by Anne Bowen 19 Mar '12

19 Mar '12

Hello, I would like to smooth an MD trajectory using the script " smoothMD.py" mentioned here (http://plato.cgl.ucsf.edu/pipermail/chimera-users/2011-September/006797.html) I can't access the attachment. Could someone send it to me? Thanks! Anne

2 1

Chimera Modeller interface trouble
by Karel Berka 16 Mar '12

16 Mar '12

Dear all, I have a strange problem with Chimera 1.6.1 Modeller (9.10) interface - when I tried to prepare a model based on template 3DDQ (CDK2 kinase), it showed that Modeller locally is running: 0%, but it stayed on this value. When I checked Task list, I have found that although it is still "running", the result was "progress: finished" with no protein model produced. I dwelwed more and I later found that Modeller script was prepared, but Modeller was unable to start due to the misspelling between sequence used for the template (probably from SEQRES) and sequence calculated from template structure. When I deleted residues that Modeller was unable to find in the structure, I was able to build the models. As for the web homology modelling interface - I was also unable to get model from Opal server, probably for the same reason. What to repair in Chimera to get Modeller interface running? because otherwise it seems like homology modelling is not working at all, which is not true. Best wishes Karel "Krápník" Berka P.S.: I have tried both Win64 and Debian64 binaries of Chimera 1.6.1. **************************************************************** RNDr. Karel Berka, Ph.D. Palacký University in Olomouc Faculty of Science Department of Physical Chemistry tř. 17. listopadu 1192/12 771 46 Olomouc tel: +420-585634769 fax: +420-585634769 e-mail: karel.berka(a)upol.cz ****************************************************************

2 1

Calculating partial charges of a charged ligand
by Christos Deligkaris 15 Mar '12

15 Mar '12

Dear Chimera users, I have a charged ligand (formal charge of +2) and I have tried using Chimera to calculate the partial charges of that molecule. Even though Chimera predicts correctly that the net charge is +2, after calculating the partial charges (using the Gasteiger method) I see that the sum of all partial charges is zero. Changing the net charge of the molecule seems to have no effect on the calculation of the partial charges. It also seems that the Gasteiger website does not offer the service of calculating the partial charges any more... ( http://www2.chemie.uni-erlangen.de/services/petra/smiles.html) Do you have any suggestions on how to calculate the Gasteiger charges so that the net charge on my molecule is +2? Thank you for your help, -- Christos

2 1

Match structures and density maps
by Sam Meyer 15 Mar '12

15 Mar '12

Good morning, I have two density maps, which I want to align for comparison. They are associated to 2 different PDB structures that I can use as references. I want to match the two PDB structures, and use the rotation/translation matrix of this alignment to move the density map as well. Is that possible ? Thanks Sam Meyer

2 1

1.6.1 release out
by Eric Pettersen 14 Mar '12

14 Mar '12

Hi all, We have put out a 1.6.1 release which, along with some minor fixes, corrects these problems: - Graphics would be corrupted or chimera would crash on recent versions of Linux, e.g., Ubuntu 11.10 and Fedora 16 - Planes (of Axes/Planes/Centroids tool) failed to start - Inappropriate "file types" (e.g. PDB ID) included in the file-type list of the main File→Open dialog - Very old sessions (pre-2006) would not restore on 64-bit Windows Please download 1.6.1 if these problems are affecting you. --Eric Eric Pettersen UCSF Computer Graphics Lab http://www.cgl.ucsf.edu

1 0

plz ans the ques
by rashmi gupta 12 Mar '12

12 Mar '12

Sir, I am working on insilico metabolomic and genomic study of edible cyanobacteria. What plan should i do to accomplish my goal. -- Rashmi research scholar School of Biomedical Engineering IT BHU

1 0

resolution for volume fitting
by Bala subramanian 07 Mar '12

07 Mar '12

Friends, I am a newbie to fitting molecular structures on a volume map. I would appreciate if someone could suggest me some guidelines to choose an optimal resolution of the map before fitting. Thanks, Bala -- C. Balasubramanian

3 4

Increasing the mesh thickness in Chimera Volume Viewer
by Pradeep Pallan 07 Mar '12

07 Mar '12

Hi, I wanted to know whether there is an option to increase the mesh thickness in Chimera? I have loaded the map in Volume Viewer, but do not know how to increase the thickness of the mesh. Thanks for your help, Pradeep -- -------------------------------------------------------------------

2 2

Linux and Chimera 1.6
by Greg Couch 07 Mar '12

07 Mar '12

If you have are running a recent version of Linux, for example, Ubuntu 11.10 or Fedora 16, then you should wait for Chimera 1.6.1 before updating or use a daily build. -- Greg

2 1

Re: [Chimera-users] [Dock-fans] Help with antechamber
by Scott Brozell 07 Mar '12

07 Mar '12

Hi, On Tue, Mar 06, 2012 at 08:04:51PM -0500, Brad Ridder wrote: > Thanks for the responses, but I am still stuck on antechamber. > > We are holding off on this python scripting; that was the plan before, but > I want to make sure I can get the individual components working properly > before I go to the effort of writing a python script to do dock prepping. > > I cannot get either the DOCK-supplied version of antechamber or the > CHIMERA-supplied version of antechamber to work correctly. ==== dock > First, I tried the DOCK distribution again. $AMBERHOME was not defined. I > used: > > *$DOCK_HOME/bin/antechamber -i $PDBFILES/3GLR.pdb -fi pdb -o > $DESKTOP/3GLRcharged.mol2 -fo mol2 -c bcc* > > This returns the error: > * > Total number of electrons: 0; net charge: 0 > > Running: $DOCK_HOME/bin/mopac.sh > Error: unable to find mopac charges in mopac.out > * It is possible that antechamber is working and that mopac really fails. Did you look in mopac.out ? Are you really trying to run antechamber on 3glr, which is a protein ? antechamber is not the way to assign charges to a protein. Read the AmberTools manual http://ambermd.org/doc11/AmberTools.pdf Here is one way to test your dock's antechamber: cd dock6/install/test/amber_score_181l make ../../../bin/antechamber -s 2 -i lig.1.mol2 -fi mol2 -o lig.1.gaff.mol2 -fo mol2 -nc 0 -c bcc This last command should emit: ... Running: $DOCK_HOME/bin/mopac.sh Running: .../dock6/bin/am1bcc -i ... ... > I then repeated this later with $AMBERHOME defined to $CHIMERA/bin/amber11, > which gives the new error message: > > *Total number of electrons: 0; net charge: 0 > > Running: $AMBERHOME/bin/divcon > sh: > /home/abhay/Desktop/rcg-research/personal-work/CHIMERA/bin/amber11/bin/divcon: > not found > Error: cannot run "$AMBERHOME/bin/divcon" of bcc() in charge.c properly, > exit* This approach - using dock's antechamber with amberhome defined to be inside chimera - is not a good one. However, you might get it to work: It is very likely that sqm, not divcon, is the qm program that should be invoked. So get help to look for the option to force sqm usage: antechamber -h ... -df am1-bcc flag, 2 - use sqm(default); 0 - use mopac ... is probably what you'll find; if so then add -df 2 to your antechamber command line. > Now, a few questions/reflections: > > 1. As per the first error message, why is mopac.out causing trouble? Is > this a file I am supposed to make? mopac.out is the output of the mopac qm program which antechamber invokes and then processes to make am1-bcc charges. > 2. As per the second error message, I have no idea what "divcon" is, or > what it is needed for. There is no divcon file that comes with the CHIMERA > amber11 directory, nor is there any divcon file in the AmberTools amber11 > directory I have on my Desktop. divcon is another qm program. ==== chimera > I then tried the version of antechamber with the CHIMERA distribution. > First I typed: > > *$CHIMERA/bin/amber11/bin/antechamber* > > And antechamber said "AMBERHOME is not set!" > > So, I set $AMBERHOME to be $CHIMERA/bin/amber11. Then I used > > *$CHIMERA/bin/amber11/bin/antechamber -i $PDBFILES/3GLR.pdb -fi pdb -o > $DESKTOP/3GLRcharged.mol2 -fo mol2 -c bcc* > > The Terminal then spams: > > *Warning: detected more than 10 Residue sequence numbers; > this may be a large multiple residue PDB file; > large multiple residue PDB files are not supported. > Continuing, but problems may be encountered.* > > And then finally stops with > > *Segmentation fault* This warning was absent in early versions of antechamber. And it likely right on target. mopac (above) probably failed because you are really trying to run antechamber on 3glr. antechamber is not the way to assign charges to a protein. > So, in short, I am utterly bamboozled. Obviously, antechamber can handle > "large multiple residue PDB files", since the Dock Prep tutorial used the > AM1-BCC charge model on the large 1ABE.pdb protein, and it worked fine. Why > is it not working for me? No, the receptor charges were assigned using residue matching. See step 1 of http://dock.compbio.ucsf.edu/DOCK_6/tutorials/struct_prep/prepping_molecule… and the chimera dockprep docs: http://www.cgl.ucsf.edu/chimera/current/docs/UsersGuide/framecontrib.html "Charges for standard residues ... are taken from Amber (details)." Read the force field chapter of the AmberTools manual and references therein for even more details. scott > On Wed, Feb 29, 2012 at 10:30 PM, Scott Brozell <sbrozell(a)rci.rutgers.edu>wrote: > > Here's another similar ligand prep script: > > http://mailman.docking.org/pipermail/dock-fans/2011-December/002591.html > > While it is a good pedagogical exercise to prep once manually, > > scripting is the only practical approach for lots of ligands. > > > > As far as running antechamber, be careful and be specific about which > > one you run: > > If you are manually running Dock's antechamber then define DOCK_HOME. > > Do not define AMBERHOME. > > See > > http://dock.compbio.ucsf.edu/DOCK_6/faq.htm#Error%20from%20amberize_ligand > > ; > > > > If you are manually running Chimera's antechamber then DOCK_HOME > > is irrelevant but AMBERHOME may be relevant. > > ( I have not manually run Chimera's antechamber, so my advice may be > > wrong. ) > > > > Also note that Dock's antechamber is circa 2006 with some key subsequent > > bugfixes. Chimera's antechamber is much more recent. > > Dock's antechamber uses mopac for qm; > > Chimera's antechamber uses sqm for qm; > > (http://ambermd.org/#AmberTools) > > I have found mopac to be more robust than sqm. > > Your mileage may vary. > > > > scott > > > > On Wed, Feb 29, 2012 at 04:08:22PM -0800, Eric Pettersen wrote: > > > This previous posting to the dock-fans list seems to be more directly > > > pertinent: > > > > > > http://mailman.docking.org/pipermail/dock-fans/2007-May/001043.html > > > > > > On Feb 29, 2012, at 2:47 PM, Trent E. Balius wrote: > > > >I am cc'ing this to the chimera list. > > > > > > > >I think that you should write a python script to run DOCK prep on > > > >the command line. See the following: > > > > > > > >http://www.cgl.ucsf.edu/pipermail/chimera-dev/2011/000784.html > > > > > > > >For other example scripts see the following: > > > > > > > > > > http://www.cgl.ucsf.edu/pipermail/chimera-users/2009-February/003539.html > > > > > > > >You should not have to specify the AMBERHOME path. This is done > > > >internally to chimera. > > > > > > > >as an aside, It is unnecessary to run antichamber with sudo. > > > > > > > >-----Original Message----- > > > >From: Brad Ridder <clustro(a)gmail.com> > > > >To: dock-fans <dock-fans(a)docking.org> > > > > > > > >My name is Brad; I am a grad student in the dept. of chemical > > > >engineering at Purdue University. > > > > > > > >We are experimenting with DOCK for various drug discovery work. At > > > >the moment, I am trying to do execute all of the "dock prep" > > > >commands through the Terminal, since it is not practical to do large > > > >volumes of Dock Prep by hand through the CHIMERA GUI. > > > > > > > >Corresponding to what is in the DOCK tutorial, I am on "Receptor and > > > >Ligand Preparation." I got reduce to work, in order to add hydrogens. > > > > > > > >Now, I am trying to get antechamber to work, in order to add the > > > >charges, but it is not working. Here is the command I am running: > > > > > > > >sudo ./antechamber -i $DESKTOP/3GLR.pdb -fi pdb -o $DESKTOP/ > > > >3GLRcharged.mol2 -fo mol2 > > > > > > > >Which returns the error: > > > > > > > >Warning: $AMBERHOME and $DOCK_HOME enviornment strings are not set, > > > >use "mopac.sh" in the work directory > > > > > > > >Cannot open CONNECT.TPL , exit > > > >The atom number exceeds the MAXATOM, reallocate memory > > > > > > > >I looked around in the dock6 directory, and found a folder called > > > >"parameters/antechamber", which contains CONNECT.TPL. I figured this > > > >is what should be the $AMBERHOME directory. > > > >For DOCK_HOME, I figured it should be the dock6 folder. > > > > > > > >Here is what I have in my .bash_profile file in my home folder. > > > > > > > >export AMBERHOME='/home/.../DOCK6.5/dock6/parameters/antechamber' > > > >export DOCK_HOME='/home/.../DOCK6.5/dock6' > > > > > > > >I also tried using the "amber11" folder that comes with CHIMERA as > > > >AMBERHOME, but to no effect. > > > > > > > >I don't actually have AMBER installed on my computer, but since > > > >CHIMERA can add charges without actually requiring the purchase of > > > >AMBER, I do not think that is the problem. > > > > > > > >I did not know what to do with "mopac.sh." I ran the script, but it > > > >says it needs input files. > > > > > > > >What am I doing wrong that is preventing antechamber from executing > > > >correctly?

2 1

Surface charge density!
by IONEL , ALINA ELENA 05 Mar '12

05 Mar '12

Hi, My name is Alina Ionel and I am PhD student. I have a pseudoatomic model (in pdb format) for a capsid protein. I did a representation of the electrostatic potential in Chimera with the option “Coulombic Surface Coloring” for a hexamer, for example. The unit of this potential are in Kcal/mol*e and the surface is colored in a range of values -10 to 10 Kcal/mol*e. I want to know the surface charge density for this hexamer (in C/m2). I need the value. Can you help me, please? Thank you, Alina

3 3

Atom distances
by Kyle Morris 05 Mar '12

05 Mar '12

Hi, Is there a way in chimera to search a PDB (or any structure file..) for a particular distance and report back in text or visually between what atoms this occurs, i.e. query whether a structure contains atoms spaced at between 5.5-5.6 Å and which ones? Thanks in advance for your help, appreciate it! Kyle

3 2

how to show backbone carbonyl oxygen
by james09 pruza 03 Mar '12

03 Mar '12

Hi All, I am using latest Chimera v1.6. I want to show the backbone carbonyl oxygen of one of the residue. The backbone is in ribbon representation. Also, how to make dotted line representation of the missing residues joining two domains of the same chain. Thanks. James

2 5

How to show the electron density of just selected residues? Thanks
by Pingyang Liu 02 Mar '12

02 Mar '12

Hi, How to show the electron density of just selected residues? Thanks I followed the Chimera tutorial, but still the electron density of the whole protein was displayed. Another question: what does the convexity of a protein structure tell? I did not find much information about convexity online. Thanks a lot. Best, Andy --

2 2

segmentation fault
by nauton lionel 02 Mar '12

02 Mar '12

Hello I send this Email, because I can't used chimera-1.6 for linux 64 bits ( 2-29-2012 version) it crash with this comment "segmentation fault", everytime I tried to load pdb or mol2 files. and I don't have this problem with chimera-1.6rc for linux-32 bits or chimera-1.6( 2-29-2012 version) and chimera-1.5.3 for linux 64bits works well on the same PC I used fedora 16... regards

2 1

Chimera 1.6 officially released
by Eric Pettersen 02 Mar '12

02 Mar '12

Hello everybody, A new production release of UCSF Chimera (version 1.6) is available: <http://www.cgl.ucsf.edu/chimera/download.html> Download is free for noncommercial use. Platforms: Windows, Mac, Linux, including 64-bit versions (recommended for working with large datasets on machines with at least 4GB memory). This will be the last Chimera release supporting Mac PowerPC platforms and OS X 10.5 (Leopard). New features include a Rapid Access interface to frequently used data and tools, new ribbon spline options, improved transparency rendering, new fetch type: Uniprot protein sequence and annotations, and several options for fitting atomic structures into density maps: global search, iterative sequential fitting, symmetric fitting, etc. More details are given below; see release notes for the full list: <http://www.cgl.ucsf.edu/chimera/docs/relnotes/1.6.html> Graphics, Depiction: * new spline options for ribbon closer to backbone atom positions * new initial model colors, pastel rather than saturated * improved transparency rendering (global single-layer on by default) * 3D labels drawn on top by default * fancy backgrounds: images from files, color gradients from Palette Editor * options for cleaner pocket/tunnel surface display (see "sop" command) Structure Analysis and Modeling: * AddH allows specifying less common states of Asp, Glu, Lys, Cys * FindHBond can label H-bonds with distances * Blast Protein shows more info for PDB hits (ligands, resolution, etc.) * Build Structure improvements for joining peptides, swapping substituents * tools to change chain IDs, renumber residues * minimization: - can be conjugate gradient (previously steepest descent only) - parts of structures can be ignored * distance and angle plotting for trajectories * Cage Builder tool for generating polyhedral cage models <http://www.cgl.ucsf.edu/chimera/docs/ContributedSoftware/cagebuilder/cagebu… > * interface to MultiFit web service (simultaneous rigid fitting of multiple structures to density, developed by the Sali Lab, UCSF) * "fitmap" command: - global search with random initial placement - iterative sequential fitting - symmetric fitting * "measure" command options: - map symmetry - surface contact area - center of mass User Interface and I/O: * Rapid Access interface to frequently used data and tools shown at startup, can be hidden/shown at any time by clicking the lightning-bolt icon <http://www.cgl.ucsf.edu/chimera/1.6/docs/UsersGuide/rapid.html> * caching fetched data is turned on by default * scaling with mouse/touchpad scrolling can be turned off (Mouse prefs) * Model Panel: - group/ungroup function, such as for NMR ensembles - functions reorganized into favorites/all * new input types: SPHGEN spheres, Maestro/Glide docking output (ViewDock) * WebGL export (experimental feature) Sequence-related: * new fetch type UniProt: protein sequence and feature annotations * Multalign Viewer: - UniProt features can be mapped onto already open sequences - Region Browser improvements for handling feature annotation regions - calculation of region overall RMSD - interface to Modeller homology modeling allows including water, other HET Enjoy! On behalf of the Chimera team, --Eric Eric Pettersen UCSF Computer Graphics Lab http://www.cgl.ucsf.edu

3 3

Re: [Chimera-users] reading SMILES description into Chimera
by Elaine Meng 01 Mar '12

01 Mar '12

Hi Boaz, You can just use the command line, for example: open smiles:CC1=CC=C(C=C1)S(=O)(=O)[O-] I'm thinking it may be a mistake that SMILES shows up as a file type in the "Open File" dialog. The only other place I can think of is in Build Structure (under Tools... Structure Editing), where the Start Structure section includes a SMILES option. In that case you would just type or paste in the SMILES string in the indicated field. 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 Mar 1, 2012, at 7:08 AM, Boaz Shaanan wrote: > Hi, > > I'm trying to read a small molecule for which I have a SMILES description into Chimera. I selected 'smiles' from the input options but it's not clear to me how to proceed. Do I have to generate a text file with the smiles description in it? do I paste the smiles description in the file name line? or what? > > Thanks for your help. > > Boaz > > > > 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 > > > > > _______________________________________________ > Chimera-users mailing list > Chimera-users(a)cgl.ucsf.edu > http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users

3 3

how to copy partial charges from a .mol2 file to another
by Thomas Evangelidis 29 Feb '12

29 Feb '12

Dear Chimera community, Is there any straight forward way to copy the partial charges from one .mol2 file to another? Thanks in advance. Thomas -- ====================================================================== Thomas Evangelidis PhD student Biomedical Research Foundation, Academy of Athens 4 Soranou Ephessiou , 115 27 Athens, Greece email: tevang(a)bioacademy.gr tevang3(a)gmail.com website: https://sites.google.com/site/thomasevangelidishomepage/

2 3

partial conservation coloring
by Eric Pettersen 28 Feb '12

28 Feb '12

Jill Voreis wrote: > I'm trying to use chimera to show conservation maps. Using pdb files > that I fetch from the database using the pdb id number and > alignments created in clustal that contain the sequence of the pdb > file, I select render by conservation. For some reason, the > conservation colors only go onto my model when it's in ribbon form. > When I try to show the surface representation or use spheres so that > I can use commans to show the SAS, the colors of my model do not > match the conservation model but rather are stuck on some sort of > coloring scheme related to the compounds that are there. The > majority of the map looks brown, like the default color before > rendering by conservation. If you could help me, I'd really > appreciate it! Hi Jill, There are a few possibilities here, though it's impossible to be certain what the problem without access to a concrete example. My best guess is that much of the surface region of your protein is poorly conserved and those areas have columns where the rows are more than 50% gap characters. For such columns no conservation score is assigned. That means that when you use Render By Attribute with default settings, those atoms / surface parts will retain their original colors. If in Render By Attribute you use the "No-value color" color well to assign a color to such regions (perhaps gray) you will be able to tell if this is what is happening. Two other possibilities are that you don't have the "Color atom" and "Color surfaces" check boxes checked [though I image you would have noticed this!] or that your protein has multiple chains and that only one is associated with the alignment and only that one is getting colored. If none of these are "it", if you can provide me with a specific example so I can reproduce the behavior then I will investigate. --Eric Eric Pettersen UCSF Computer Graphics Lab http://www.cgl.ucsf.edu

1 0

Not crossfade but coinciding bond command...
by Kevin D. Walker 28 Feb '12

28 Feb '12

I temporarily disabled a bond command and found this to be the cause of the selected bond problem I am having. I don't think crossfade has anything to do with it.

1 0

Animate a Reacton Mechanism?
by Kevin D. Walker 28 Feb '12

28 Feb '12

I constructed a movie for the first part of a reaction mechanism (substrate docking). I would like the animated mechanism to continue as follows. The substrate donates a proton to a tyrosine residue of the enzyme, forming an intermediate. The intermediate rotates 180 (I have accomplished this part) and the proton rebounds to the intermediate forming a new product, which is released. I can animate the product release -- it's the reverse of substrate docking. The part of the animation I cannot envisage is how to smoothly show or morph the bond transitions (bond breaking and making) during the reaction progression. IAny advice on whether Chimera or other ancillary software can do this is appreciated.

4 6

using CalcAttr from a Python script
by Roger L. Chang 27 Feb '12

27 Feb '12

Hi, 1) I need to get residue.areaSES for two separate selections (all polar amino acids and all nonpolar amino acids). I have successfully generated a list of residue.areaSES using the GUI simply by selecting the residues of interest and then using the attribute calculator typing "residue.areaSES" in the formula box of the Attribute Calculator, which outputs the list to the Reply Log. This is the data that I need, but I need to get it for >2000 PDB structures and would like to automate this analysis through a Python script. I am somewhat new to Python and Chimera and don't yet know all the ins and outs yet. I dug through the code and found a method "getResidueAttr" in the "ResidueDataWrapper" class from "CalcAttr". This seems to be what I need, but I don't know the proper syntax to call it from my Python script. Could someone please provide an example? 2) I would like to sum the residue.areaSES values once I have them. I tried "sum(residue.areaSES)" in the formula box of the Attribute Calculator in the GUI, but this returns the error "Formula evaluation error: expecting vector argument in function". I'm not sure again on the proper syntax to use "sum()" in the Attribute Calculator. If I knew how to generate/save the residue.areaSES data from the command line, I could just parse the saved file and sum the values in Python. How can Attribute Calculator formulas be used from the command line, and if they can't, how can I save the residue.areaSES data to file from a Python script? Thanks in advance, Roger -- Roger L. Chang Graduate Student Bioinformatics and Systems Biology Graduate Program University of California, San Diego

3 2

Coloring molmap surfaces
by Tom Goddard 24 Feb '12

24 Feb '12

> Greetings Tom > > How can color the molmap surface of certain atoms/residues. > > Many Thanks, Vijay > Hi Vijay, Use the Color Zone dialog to color a molmap surface to match nearby atom colors (menu Tools / Volume Data / Color Zone). You select the atoms you want to use for color (e.g. menu Select / Select All) choose a distance range from those atoms to color in the Color Zone dialog and press the Color button. In Chimera 1.6 (release candidate available on the web) or daily builds there is a command version of color zone. For example scolor #0 zone @ca range 5 says to color molmap surface #0 to match C-alpha atom colors within a distance of 5 Angstroms. To avoid jagged edges in the coloring you should use the molmap gridSpacing option to specify a grid size about 1/5 to 1/10 the map resolution. The default value is 1/3 of map resolution making a fairly coarse surface mesh. Colors are only assigned to the mesh vertices, so a fine mesh is needed for smooth boundaries between different color regions. Tom

1 0

batch peptide builder
by Xiao-Ping Zhang 23 Feb '12

23 Feb '12

Dear Chimera developers, I would like to generate some peptides from a protein sequence. The peptide builder in Chimera is a very nice application. I wonder if there is a stand alone script which I can use to generate series peptides from the command line. Thank you for your help. Xiao-Ping Zhang

2 1

Drawing arrows without using .bild file
by Victor Muñoz 21 Feb '12

21 Feb '12

Dear all, I'm trying to make a script that draws an arrow on every atom and every frame during a trajectory using the MD movie. Right now, the script creates a .bild file to draw the arrows, but it is not a very good approach since it is creating a new file for every frame of the trajectory. I was wondering if there's a way to do the same thing without having to create the .bild file. I already checked the "shape" command but you can´t use it to draw an arrow. Maybe I should make a combination between a cylinder and a rectangle to draw an arrow? Is there any other easy way to do that? Thanks in advance. Kind regards, --

3 2

bug
by Eifler, Jay Q. (UMKC-Student) 21 Feb '12

21 Feb '12

Whew, too much work to get you guys a message. The bug report form was lengthy. Anyway, can't join models in your latest 1.6 version but 1.5 whatever is fine. Should have never updated.

3 2

(no subject)
by Valia Dimitrova 21 Feb '12

21 Feb '12

Hi, I am PhD student, and I am working with program Gaussian, I want ask you some questions. I want to make maps with electrostatic potential, I saw that with this program may be make very beautiful pictures. What type of files I have to have, so I can do these pictures? I have ".fchk" files and ".cube" files. May I use them to make this? I will be very appreciated if some one can help me. thank you in advance Valia Nikolova

2 1

how to find
by Bishwajit Das 18 Feb '12

18 Feb '12

suppose i want so select ARG 188 and TRP 144 and also want to see the distance between them ....how i do this in chimera..please help me. -- BISHWAJIT DAS KOLKATA,West Bengal INDIA* * --

2 3

bug report
by Longfei Wang 16 Feb '12

16 Feb '12

Hi there, For nogui Chimera: Calls with script in nogui mode create the following traceback. terry@daisy[14]% chimera --nogui --script "coord_for_cluster_ATP_ bothChiFlip.py" Traceback (most recent call last): File "/nfs/turing2/labs/sbgrid/programs/i386-linux/chimera/1.6rc/share/__main__.py", line 69, in <module> value = chimeraInit.init(sys.argv) File "/nfs/turing2/labs/sbgrid/programs/i386-linux/chimera/1.6rc/share/chimeraInit.py", line 618, in init chimera.initializeGraphics() File "/nfs/turing2/labs/sbgrid/programs/i386-linux/chimera/1. 6rc/share/chimera/__init__.py", line 2114, in initializeGraphics bgprefs.initialize() File "/nfs/turing2/labs/sbgrid/programs/i386-linux/chimera/1.6rc/share/chimera/bgprefs.py", line 211, in initialize onDisplay=_bgPrefDisp, onHide=_bgPrefHide) File "/nfs/turing2/labs/sbgrid/programs/i386-linux/chimera/1. 6rc/share/chimera/preferences/base.py", line 567, in register self._loadCategory(category, cat, inherit, convert) File "/nfs/turing2/labs/sbgrid/programs/i386-linux/chimera/1. 6rc/share/chimera/preferences/base.py", line 482, in _loadCategory cat.load(self._loadDict[name]) File "/nfs/turing2/labs/sbgrid/programs/i386-linux/chimera/1. 6rc/share/chimera/preferences/base.py", line 249, in load opt.set(value, asSaved=True, fromPref=True) File "/nfs/turing2/labs/sbgrid/programs/i386-linux/chimera/1. 6rc/share/chimera/preferences/base.py", line 64, in set self._callback(self) File "/nfs/turing2/labs/sbgrid/programs/i386-linux/chimera/1.6rc/share/chimera/bgprefs.py", line 53, in _backgroundColorCB Midas.background(color=color, setMethod=False) File "/nfs/turing2/labs/sbgrid/programs/i386-linux/chimera/1.6rc/share/Midas/__init__.py", line 966, in background method = chimera.viewer.Solid AttributeError: '_chimera.NoGuiViewer' object has no attribute 'Solid' Best Longfei -- *Longfei Wang * Postdoctoral Fellow Harvard Medical School 250 Longwood Ave Boston, MA 02115

2 1

Re: [Chimera-users] radius of gyration
by Eric Pettersen 15 Feb '12

15 Feb '12

On Feb 14, 2012, at 6:07 PM, Daniel Gurnon wrote: > Hi Eric, > The radius of gyration interest stems from a little project I'm > starting with an undergrad. I'm interested in the radius of gyration > because I'm curious about how accurate our pictures of biomolecules > are, at least the ones that come from crystal structures. I'd like > to compare predictions of size-related behavior derived from the pdb > with measurements of size in solution, using gel filtration, SAXS, > or light scattering. That said, I don't know if I want center of > geometry or principal axis! Any advice? Hi Dan, I'm no expert, but it sounds like you want the radius of gyration from the center of geometry. Nonetheless, I'm cc'ing this to the chimera-users in case any others have thoughts. In an old posting to chimera-users, Elaine Meng provided a Fortran program to compute the radius of gyration: http://plato.cgl.ucsf.edu/pipermail/chimera-users/2009-March/003648.html I've attached a Python script that will print the center of geometry and radius of gyration from that center in Chimera's reply log. Run the script simply by opening it with File->Open. Depending on whether you feel nearby waters and ions will contribute to the measurements you are trying to model, you may or may not want to delete water and ions before running the script. --Eric Eric Pettersen UCSF Computer Graphics Lab http://www.cgl.ucsf.edu

2 2

radius of gyration
by Daniel Gurnon 14 Feb '12

14 Feb '12

Hi all, Can Chimera determine (or has a script been written to determine) the radius of gyration for a pdb structure? Dan -- ____________________________ Daniel Gurnon, Ph. D. Associate Professor of Chemistry and Biochemistry DePauw University Greencastle, IN 46135 p: 765-658-6279 e: danielgurnon(a)depauw.edu

3 2

Re: [Chimera-users] possible Chimera bug?
by Ofeliia Mada 14 Feb '12

14 Feb '12

TkkyuluR X,=)RSISXSN:-P:(:=(=(=(=(=(:O(, G Sent from my Windows Phone esdcqppwevemeem

1 0

fly command
by rswett＠chem.wayne.edu 14 Feb '12

14 Feb '12

Is anyone else having issues with the fly command lately? Seems like since RC 1.5ish flying to a savepos ignores the # of frames it's set to fly over. i.e. fly 50 in (where in is a saved position) does the fly in one frame, appearing the same as fly in or fly 1 in Additionally fly in 50 has the same effect. Has anyone else run into this? I have a few saved sessions i'd like to animate, and the older versions are having none of it. Suggestions? ~Rebecca

2 1

Re: [Chimera-users] meaning of the ANISOU PDB parameters
by Eric Pettersen 13 Feb '12

13 Feb '12

On Feb 10, 2012, at 1:55 PM, Julian Chen wrote: > hi Eric, > > First of all, hello, I was in Stroud's lab for many years, until > 2003, when I left SF. Hi Julien! > I am in a rather interesting situation where we have been dealing > with neutron and X-ray refinement, and we have been able to do > anisotropic refinement on selected hydrogen atoms, from having high > resolution neutron data. I've been using Chimera. > > I am still somewhat confused as to the meaning of the six ANISOU > terms, I understand that the first three describe the three > orthogonal axes of the ellipsoid, and that the last three somehow > relate this to the coordinate space. The six terms actually describe a 3x3 matrix that is symmetric about the diagonal. The eigenvectors and eigenvalues of that matrix are the atomic displacement axes and the mean squares of the displacements respectively. The first thing to know is that the values shown in the ANISOU records are scaled by 10**4. The other thing to know is that the six numbers correspond to these positions in the matrix: 1,1 2,2 3,3 1,2 1,3 2,3. Since the eigenvalues are mean squares, one typically works with the square roots of the eigenvalues. > Specifically, I want to look at the ANISOU for the H/D atoms, and > want to try to correlate it to its bonded neighbor (the C or the N > or the O). If possible, I'd like to extract some vector for the > 'directionality' of the aniso atoms and relate them to one another. > I want to see if the vibrational motions of the hydrogen are > correlated to the vibrational motions of its bonded neighbor. Is > there anything in CHIMERA that could do something like this, or some > appropriate mathematical formula that could output the right sort of > information? The ellipsoids option in the program gives a > qualitative description, but I'd like to get something a bit more > quantitative in nature. Depending on how comfortable you are with programming in Python, you could either write your own external code based on the info above, or write something that runs in Chimera. The latter does have the advantage of you not having to write PDB-parsing code an so on. In Chimera, atoms with corresponding ANISOU records have an 'anisoU' attribute which is a 3x3 matrix (the matrix is often referred to as 'U' in the literature, hence ANISOU record). This code snippet (basically copied from the Chimera's Aniso module) will extract the eigenvectors/values for an atom and take the root of the values: from numpy.linalg import svd ignore, lengths, axes = svd(atom.anisoU) from numpy import sqrt lengths2 = sqrt(lengths) Let me know if you need more info or help with scripting. --Eric Eric Pettersen UCSF Computer Graphics Lab http://www.cgl.ucsf.edu

2 2

Visualising AutoLigand affinity points
by George Tzotzos 09 Feb '12

09 Feb '12

Hi, I'm attaching a jpg image of ligand affinity map fill points generated by AutoLigand. My question is whether there's a way to depict these fill points as "disconnected" spheres. Actions/spheres generates something like I'd like to see but the size of the spheres is too large. Regards George

2 2

pi-cation interaction
by Bishwajit Das 09 Feb '12

09 Feb '12

how to find cation pi-interaction in a docked molecule.......plz help me -- BISHWAJIT DAS KOLKATA,West Bengal INDIA* *

2 1

How to select
by Bishwajit Das 09 Feb '12

09 Feb '12

suppose i want so select ARG 188 and TRP 144 and also want to see the distance between them ....how i do this in chimera..please help me. -- BISHWAJIT DAS KOLKATA,West Bengal INDIA* *

1 0

Re: [Chimera-users] Rotamer copying- Specific Residues only
by Eric Pettersen 08 Feb '12

08 Feb '12

On Feb 8, 2012, at 8:27 AM, Ahir Pushpanath wrote: > Eric, although the script is highly useful for what I previously > asked, I have come to the realisation that simply copying all the > apo chi angles to the holo form may not be ideal for me in this > particular situation. I noted that of around 200 residues modelled, > around 90 of them are modelled different to the apo, the rest are > modelled similar to the apo form. Of these 90, only around half are > truly changed in the holo form(literature evidence from homolog), > therefore I do require a script that is able to accept an input of > which residues to copy the chi angles for, I think this will be > useful for other users who want to copy chi angles from apo-holo > forms and other applications, and probably could be included in > later build versions of chimera, under structure editing. > > Do you think you could help me with this? I am only learning python > now, so my scripting is not at a good level for these types of things. I've attached a modified script that will only copy the chi angles for residues that are selected in the lowered-numbered model. It's up to you to select them somehow. :-) > On another note, do you have any experience with using chimera to > model apo-holo transitions any other way without using Modeller? I don't. Maybe others on the list will have suggestions... ? > Also, does chimera have the ability to do docking of ligands, i.e. > given you know the apo crystal structure, can you tell it to dock a > ligand and change conformation of protein,(similar to autodock)? Currently Chimera cannot dock ligands. If our next 5 years of funding gets approved by the NIH, we have proposed interfacing Chimera to an Autodock Vina web service that would allow docking to be carried out directly from Chimera. Right now all you can do is prepare the input for / analyze the output of external docking programs in Chimera --Eric Eric Pettersen UCSF Computer Graphics Lab http://www.cgl.ucsf.edu

1 0

how to find hydrophobic interaction
by Bishwajit Das 08 Feb '12

08 Feb '12

how to find hydrophobic interaction in a docked protein structure......please help me.... -- BISHWAJIT DAS +91-9433237920 KOLKATA,West Bengal INDIA* *

2 1

how to setup a default settings for a new pdb model?
by Xing Zhang 07 Feb '12

07 Feb '12

Hi, I have many different pdb models that need to be displayed at exactly the same settings, including color, background, stick scale and so on, so I want to set all these settings for the first model, and then apply these settings to a newly opened pdb model. Is there any way to save and load a default setting for a new pdb model? Thanks much. Xing

2 2

Rotamer copying
by Ahir Pushpanath 07 Feb '12

07 Feb '12

Hello there, I have a question regarding MODELLER in Chimera and the ability of chimera to adjust the model using some inbuilt features quickly(i.e. automatically, instead of interactively) Firstly, I am trying to model the "holo" form of an enzyme with known apo-form structure. I did this usinig the two independent domains of the apo form as a template (in a bid to get the domains accurate) and for the holo-form closure I used a very similar homologue that has been solved in the holo form. The results were good and I got a decent model that had closed up like the homolog. However, I noticed that in the model a lot of the sidechains have not been modelled with the right rotamer. This is odd, as I gave it the two independent domains as a template, so it had all the information it needed to get that perfectly right. Either way, this means, I need to manually "copy" the rotamers from the apo form onto the "holo model", which is time consuming. I know from the homolog that only two or three residues in the active site change rotamer upon apo-holo conformational change. So, I just want an easy way to edit the rotamers in the holo-model to copy the rotamers of the apo-form for all residues. Is there an easy way to do this, is it possible to this with a script? I have no idea how to script, mind you! Ahir

2 3

associating multiple structures and sequences
by Uppal, Karan 07 Feb '12

07 Feb '12

Hi, I am writing a custom extension based on MultiAlign Viewer that highlights regions in a multiple sequence alignment based on numerical values stored in an external file. I modified the RegionBrowser.py script, and added a function that works just like "showSS". It works fine for one structure and sequence combination. However, I am trying to associate multiple copies of the same protein structure and its sequence. So, the structures are identical and so are the sequences in the sequence alignment, but each copy represents different information. My problem is that the sequences in the alignment viewer dialog box are highlighted in a random fashion. And, I want the sequences to be highlighted in an ordered fashion starting from first sequence to last sequence since we want to be able to map the highlighted regions back to the external file. Is there a way to control the order in which the sequences are highlighted in the alignment viewer dialog? Or is there a way to get the model # associated with a particular sequence and using the model # to highlight the alignment viewer? I would highly appreciate your help. Thank you, Karan Uppal Application Developer Emory University ________________________________ This e-mail message (including any attachments) is for the sole use of the intended recipient(s) and may contain confidential and privileged information. If the reader of this message is not the intended recipient, you are hereby notified that any dissemination, distribution or copying of this message (including any attachments) is strictly prohibited. If you have received this message in error, please contact the sender by reply e-mail message and destroy all copies of the original message (including attachments).

2 2

minimization with nonstandard atoms
by Jon N. Rumbley 07 Feb '12

07 Feb '12

I was hoping to minimize a homology model containing two heme groups. To understand how I might do this I have been starting with the known structure 1a6m. As suggested in the documentation I begin with Dock Prep and calculate charges etc. To my surprise Chimera alpha version 1.6 (build 34702) actually assigned charges to the heme. The problem is with the minimization itself. Within Minimize Structure I turn off all of the Dock Prep steps and simply try to minimize the structure but get an error,"Element Fe (atom #0:154.A@FE) is not currently supported". I was assuming the charge calc was going to be the hard part but now I am stuck at this subsequent step. Is there something I can do to get the heme atoms recognized or is this simply not possible? Thank you, Jon Rumbley Dr. Jon N. Rumbley Assistant Professor Dept. of Pharmacy Practice and Pharmaceutical Sciences University of Minnesota-Duluth 127 Life Sciences 1110 Kirby Dr. Duluth, MN 55812 Phone: 218-726-6047 email: jrumbley(a)d.umn.edu

3 2

Manipulating representation during MD movie
by Jean Didier Pie Marechal 07 Feb '12

07 Feb '12

Hi guys, I have a MD movie I generated that I want to record. However, I would need to ask on the graphical representation at given frames (i.e. frame 60) so that the system changes (i.e. change the representation en stick or add a given label). Is there a way to do so directedly? Or should I go for a command file? In this case, how would I mix the MD movie with my own succession of events? 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: chimera-users-request(a)cgl.ucsf.edu Data: Dilluns, Febrer 6, 2012 9:00 pm Assumpte: Chimera-users Digest, Vol 106, Issue 7 > Send Chimera-users mailing list submissions to > chimera-users(a)cgl.ucsf.edu > > To subscribe or unsubscribe via the World Wide Web, visit > http://plato.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. Rotamer copying (Ahir Pushpanath) > 2. Re: Rotamer copying (Eric Pettersen) > > > -------------------------------------------------------------------- > -- > > Message: 1 > Date: Mon, 6 Feb 2012 18:51:22 +0000 > From: Ahir Pushpanath <ahir29(a)googlemail.com> > To: chimera-users(a)cgl.ucsf.edu > Subject: [Chimera-users] Rotamer copying > Message-ID: > <CAMVKkxyEoXbnQUo7hxa+m2-u5wLCvMJn7DDY=hzZPJEbC_OvSw(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Hello there, I have a question regarding MODELLER in Chimera and the > ability of chimera to adjust the model using some inbuilt features > quickly(i.e. automatically, instead of interactively) > > Firstly, I am trying to model the "holo" form of an enzyme with known > apo-form structure. I did this usinig the two independent domains > of the > apo form as a template (in a bid to get the domains accurate) and > for the > holo-form closure I used a very similar homologue that has been > solved in > the holo form. The results were good and I got a decent model that had > closed up like the homolog. > > However, I noticed that in the model a lot of the sidechains have > not been > modelled with the right rotamer. This is odd, as I gave it the two > independent domains as a template, so it had all the information it > neededto get that perfectly right. > > Either way, this means, I need to manually "copy" the rotamers from > the apo > form onto the "holo model", which is time consuming. I know from the > homolog that only two or three residues in the active site change > rotamerupon apo-holo conformational change. So, I just want an easy > way to edit > the rotamers in the holo-model to copy the rotamers of the apo-form > for all > residues. > > Is there an easy way to do this, is it possible to this with a > script? I > have no idea how to script, mind you! > > Ahir >

2 1

how to find salt bridge
by Bishwajit Das 04 Feb '12

04 Feb '12

how to find salt bridge in a docked structure........plz help me -- * BISHWAJIT DAS* *+91-9433237920* KOLKATA

2 1

Coloring electrostatic potential
by George Tzotzos 03 Feb '12

03 Feb '12

Hi everybody, I'm using SiteHound for ligand binding identification sites. I'm comparing the highest affinity sites for two superimposable receptors. I'm attaching an example where the sites favourable to protein 1 are shown red, whereas those for protein two blue. I've opened in Chimera the file that contains the red and blue clusters (difference.pdb) as well as the corresponding potential file (potential.dx). I used the Volume Viewer to different affinity areas but had .... no luck. Any help would be greatly appreciated Best regards George

2 3

export saved session
by Gabor PAPAI 02 Feb '12

02 Feb '12

Dear Chimera team, Is it possible to export a session with all of its files to a different location - for e.g. to a USB flash drive for portability reason? I'd like to open a session saved to a server on an offline computer (during travel). I know that I can do it manually just I'm searching for quick solution. Thank you in advance for the help! Best regards, Gabor --- Gabor Papai IGBMC Integrated Structural Biology 1, rue Laurent Fries, Bp 10142 67404 Illkirch, France Tel: +33-3-88655748

3 2

selection based on ring names
by Bala subramanian 02 Feb '12

02 Feb '12

Friends, I have two ligands with a common benzoxadiazole ring. I want to do the rmsd comparison only on the ring, but the names of the atoms forming the benzoxadiazole ring are different in both the ligands. Is there any way to use ring name as a string the compare the rmsd. ie. rmsd of benzoxadiazole of ligand 1 and benzoxadiazole ring of ligand 2 -- C. Balasubramanian

2 1

stereo with zalman monitor
by Jose Miguel Mancheno 01 Feb '12

01 Feb '12

Dear colleagues, I am trying to see stereo images with chimera (candidate version 1.6 for mac) using a Zalman monitor ZM240W which works OK with COOT, which is connected to an iMac. As far as I know, I have to select (within the camera mode) the option: row stereo, right eye even (or odd). However, none of these options provided stereo images. If I select sequential stereo (or reverse sequential stereo) I see the image flickering, as expected, but this is not synchronized with the zalman glasses. Any suggestion about how to proceed to see in stereo mode? Thank you very much in advance. José Miguel --

2 1

Selecting equivalent surfaces
by Ankush Aggarwal 01 Feb '12

01 Feb '12

Hi, I am opening up a viperDB file of a T=3 virus in Chimera. I want to select the 60 equivalent copies of a protein, so that I can change their color together. I am able to select selecting each of them separately. Is there a smarter way to do it? Also, after I rotate the view, is there a command to reset it to the default. Right now, I am "undoing" the move, but it can be a pain if I have moved the model a few times. Thanks, Ankush.

2 2

Notebook
by Miguel Ortiz Lombardía 01 Feb '12

01 Feb '12

Dear all, Was the notebook feature implemented in Chimera? I found a message from some time ago but nothing else? Any plans to include a session notebook? I would ask for something really simple, just apt for text notes. Currently, I create text files with the same name than sessions but with .txt extension, however I think it would be useful to have this built-in in Chimera. Best regards, -- Miguel Architecture et Fonction des Macromolécules Biologiques (UMR7257) CNRS, Aix-Marseille Université Case 932, 163 Avenue de Luminy, 13288 Marseille cedex 9, France Tel: +33(0) 491 82 55 93 Fax: +33(0) 491 26 67 20 mailto:miguel.ortiz-lombardia@afmb.univ-mrs.fr http://www.afmb.univ-mrs.fr/Miguel-Ortiz-Lombardia

2 3

Re: [Chimera-users] HIV ligands
by Eric Pettersen 31 Jan '12

31 Jan '12

On Jan 30, 2012, at 5:27 PM, Forbes J. Burkowski wrote: > Hi: > > I think I have found a bug or least a small inconsistency... > > If you fetch 3NDW and do the menu invocation: > > Select... Residue... > > The list in the submenu offerings includes RIT. > > In a command line, I can use: > > sel :RIT > > to get the same selection. > > Now, repeat the above invocations for 2FDE. > > The Select... Residue... leads us to believe that the ligand can be > accessed as 385. But when I use: > > sel :385 > > I do not get a selection. > > This alternate strategy to do the selection is important for me > because I > am getting the same problem when I put this command into the > runCommand > string within a script. > > The same problem seems to exist for 3SM2 and 3SO9. Both have the > ligand > designated by a number so I have a suspicion that one must do > something > different in the sel command when a number is used to designate the > ligand > instead of a character triplet (such as RIT). > > Is this a bug or am I missing something in the syntax of the sel > command? Hi Forbes, Pedantically speaking, the single-colon (e.g ":532") syntax takes a residue sequence number afterward, whereas the lesser-known double- colon syntax (e.g. "::ala") takes the residue name. Now, as a convenience, if you type something that is not an integer after the single colon, then it assumes you mean the residue name and switches over to interpreting it like the double colon does. But in your case you have character strings that you always want interpreted as residues names regardless of whether they look like integers, so you should just use "::" all the time. --Eric Eric Pettersen UCSF Computer Graphics Lab http://www.cgl.ucsf.edu

1 0

Re: [Chimera-users] how2 select named selections using command-line [chimera1.5.3]
by Eric Pettersen 31 Jan '12

31 Jan '12

On Jan 31, 2012, at 2:34 PM, frederick lee wrote: > thanks VERY much eric and elaine. > > last question: how do i find out about the alias definitions, if > any, stored > in a session from within chimera instead of looking inside the > python script > (of that stored session). Hi Fred, "alias" with no arguments will list all current aliases in the reply log. --Eric

1 0

how2 select named selections using command-line [chimera1.5.3]
by frederick lee 31 Jan '12

31 Jan '12

dear chimera users, i can defined named selections and they work under the select menu. however, i canNOT use the command line to select them; e.g., i have 4 named selections 0a,0b,1a,1b dis 0a | 1a works sel 1a NOT working sel :1a NOT working (said selection cleared) i like to be able to select multiple named selections and further apply filters on them; e.g., "select 1a and 1b and atoms=side chain atoms" how do i accomplish this? as always thank you very much. fred

3 3