October 2016 - Chimera-users - RBVI Mailing Lists

Chimera crashes on Mac OS Sierra (10.12)
by Tom Goddard 20 Oct '16

20 Oct '16

The Chimera developers have received 8 Chimera crash reports on the new Mac OS Sierra (10.12) somehow related to scrollbars. You may want to hold off updating to Sierra if you make heavy use of Chimera. We have not yet been able to reproduce the crash. It is very likely a bug in the Mac Tk window toolkit that Chimera uses, so even when we figure out how to reproduce it it may take time to figure out a fix since the problem is not in our code. Tom

1 1

creating O2 molecule in Chimera
by Elaine Meng 20 Oct '16

20 Oct '16

> Subject: Query about chimera > Date: Thu, 20 Oct 2016 06:47:23 +0000 (UTC) > From: Shumaila Khan <shumailakhan432(a)yahoo.com> > > Hello, > I m student from Quaid-i-Azam University, Pakistan, doing my MS research project and I have some problems while working on chimera. I want some help regarding Attachment of oxygen atom to plant hemoglobin.Kindly please help me out how will I do that. kindly tell me how to draw oxygen molecule on chimera? While drawing Oxygen molecule it has only only Oxygen atom and two hydrogen atoms attached. Kindly help me out. Waiting for your kind response. Hi Shumalia, Your question was forwarded to me — Chimera questions should be sent to this address, chimera-users(a)cgl.ucsf.edu There are several different approaches you could use. One way is to use command: open smiles:oo If you hold the mouse over it you can see a balloon that says the bond length is 1.004. You could change the bond length to a better value by (1) selecting the bond (Ctrl-click on the bond so it has a green outline) (2) using Build Structure to change the length of the bond. Menu: Tools… Structure Editing… Build Structure, and in that dialog change from the “Start Structure” section to the “Adjust Bonds” section. Now you can use the slider or type in a value to change the length to what you want. <http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/editing/editing.h…> That would make the oxygen molecule a separate model from your hemoglobin. You would probably want to first move the oxygen molecule by hand to the right position, and then (optionally) combine it together with the hemoglobin to make a single model depending on what you are planning to do. If you don’t need them in the same model, there is no need to combine. <http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/mouse.html#activedef> <http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/mouse.html> <http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/combine.html> Or, another totally different approach is to find some other hemoglobin structure that already has an oxygen molecule and superimpose the two structures to get the oxygen in the right place, and then delete all the other atoms of the other structure except for the oxygen molecule. <http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/superposition.html> I hope this helps, Elaine ---------- Elaine C. Meng, Ph.D. UCSF Computer Graphics Lab (Chimera team) and Babbitt Lab Department of Pharmaceutical Chemistry University of California, San Francisco

1 0

chimera 1.11.2rc crashes
by Yanni Goudeuk 18 Oct '16

18 Oct '16

Hello Chimera developers I've got the latest chimera (1.11.2rc), installed on windows 7 (64bit). My card is Nvidia Quadro K2000 and is using the latest drivers: 369.49 - *3D OpenGL stereo* is *enabled* under "*manage 3D settings"* in *nvidia Control Panel* My problem is that chimera crashes a few minutes after I open a .pyc file. Sometimes it will crash when I add a map in "*model panel*", sometimes when I rotate the structure, sometimes when I go though the tabs under " *Favourites*" --> "*Side view*". This is happening with other .pyc files too. I can also confirm that this is happening with older versions of chimera. Please have a look at the Application error log from event viewer: The program pythonw.exe version 0.0.0.0 stopped interacting with Windows and was closed. To see if more information about the problem is available, check the problem history in the Action Center control panel. Process ID: ab0 Start Time: 01d2246bcf21dbd4 Termination Time: 46 Application Path: C:\Program Files\Chimera 1.11.2rc\bin\pythonw.exe Report Id: 3c965b43-905f-11e6-bd99-b8ca3a8bb0d4 If you have any suggestions please let me know. Many thanks

3 7

Guidance in finding Interactive Aminoacids with a ligand
by Shreyas Iyer 17 Oct '16

17 Oct '16

Respected Sir, I am unable to find the interactive aminoacids with the ligand present in a protein. Please guide me how to determine interatcive aminoacids. Awaiting a positive reply Yours Sincerely, Shreyas M Iyer

2 1

MatchMaker algorithm
by 陈柳青 15 Oct '16

15 Oct '16

Respected Sir/Madam The ucsf matchmaker algorithm is "Needleman-Wunsch" and "Smith-Waterman", but other program use SSM algorithm, such as ccp4mg, cuemol et al., can chimera add SSM algorithm? and whether these algorithm will affect results differently? sincerely Liuqing chen -- Liuqing Chen School of Life Science & Biotechnology Minhang Campus Shanghai Jiao Tong University No. 800 Dongchuan Road Shanghai 200240 P.R. China Tel: +86 21 34207247

2 1

MD trajectory analysis
by George Tzotzos 13 Oct '16

13 Oct '16

My task is to identify ligand interactions with protein residues throughout an MD trajectory. I thought of something along the lines sel ligand z<3.5 #for vdW interactions, etc. I don’t think this can be done using the MD movie interface but I may be wrong. I then tried the following: 1. sel protein/hold selection steady 2. sel ligand 3. calculate occupancy (MD Movie/Analysis) 4. Chimera/select/clear selection 5. turned the Volume Viewer to mesh/level2. What I obtained is a mesh grid enclosing solid surfaces (see attached snapshot). I’m not able to get rid of these surfaces? I think that the problem ma be in step 2 above (sel :ligand). My ligand is N,N-Diethyl-meta-toluamide. Should I select the aromatic ring part separately of the amide part, run the above process and then repeat selecting the amide part? I also tried MD movie/Analysis/Residue interaction network, first selecting the ligand. The results I obtained are not satisfactory. Is there a better way of dealing with my problem? Thanks in advance for your kind help George

2 5

Re: [Chimera-users] question
by Jaynes, Jesse 12 Oct '16

12 Oct '16

Hi Folks: I am having difficulty in converting C-terminal carboxyl to C-terminal amide on a beta sheet of 16 amino acids in length. When I try to get the complete amino acids displayed, it doesn't work. It just displays the main peptide cylinder. Thanks for your help. Jesse Jaynes On Thu, Oct 6, 2016 at 6:11 PM, Jaynes, Jesse <jjaynes(a)mytu.tuskegee.edu> wrote: > Hi Folks: > > I am interested in generating a 23 amino acid peptide with several > ornithines present. Is it possible? If not can I construct it with > glutamine as a replacement for ornithine and then just take the carbonyl > group and replace the oxygen with 2 hydrogens. > > I just renewed my registration. > > Thank You, Jesse Jaynes > > On Thu, Oct 6, 2016 at 3:22 PM, Jaynes, Jesse <jjaynes(a)mytu.tuskegee.edu> > wrote: > >> Hi Folks: >> >> I am interested in generating a 23 amino acid peptide with several >> ornithines present. Is it possible? If not can I construct it with >> glutamine as a replacement for ornithine and then just take the carbonyl >> group and replace the oxygen with 2 hydrogens. >> >> Thank You, Jesse Jaynes >> > >

2 3

Creating a dummy atom on a defined position
by Subha Kalyaanamoorthy 11 Oct '16

11 Oct '16

Hi There, Wondering if there is a way to create a dummy atom on a specified position(xyz coords) using a chimera command line? Thanks, Subha -- Dr. Subha Kalyaanamoorthy Post Doctoral Fellow Faculty of Pharmacy and Pharmaceutical Sciences University of Alberta Edmonton, Canada.

3 3

Computing 'basic' structure characteristics
by Healey, Joe 11 Oct '16

11 Oct '16

Hi again Chimera Team, I have 2 questions, one specific and one quite general. If you see the emails below, I was speaking with Jaime about computing basic structure characteristics (pKa, net surface charge, proportions/predominance of 2ary structure etc). The general question I have is what suggestions might you have for functions within chimera capable of computing any or all of these features? (ideally commandline based so that I can run this over many PDB structures). The specific question that pertains to the below is that Jaime mentioned each residue has an int value that corresponds to a particular structure designation, ints from -1 to 6. Could you let me know what the correspondence is between these values and actual structure? Joe Healey M.Sc. B.Sc. (Hons) PhD Student MOAC CDT, Senate House University of Warwick Coventry CV47AL Mob: +44 (0) 7536 042620 | Email: J.R.J.Healey(a)warwick.ac.uk Jointly working in: Waterfield Lab<http://www2.warwick.ac.uk/fac/med/research/tsm/microinfect/staff/waterfield…> (WMS Microbiology and Infection Unit) and the Gibson Lab<http://www2.warwick.ac.uk/fac/sci/chemistry/research/gibson/gibsongroup/> (Warwick Chemistry) Twitter: @JRJHealey<https://twitter.com/JRJHealey> | Website: MOAC Page<http://www2.warwick.ac.uk/fac/sci/moac/people/students/2013/joseph_healey> ________________________________ From: Jaime Rodríguez-Guerra <jaime.rodriguezguerra(a)uab.cat> Sent: 08 October 2016 09:55 To: Healey, Joe Subject: Re: [Chimera-users] Recursive structure matching and acquisition of descriptive numbers Hi again Joe! No problem at all; I'll be pleased to help you. For charges, I think Chimera can calculate them as part of the Surface/Binding analysis> Dock Prep wizard. You also have some information about 2nd structure inside chimera.Residue objects, specifically ssId attribute. However, these are int values, and I don't really know which ss is each (ranging from -1 to 6, as I have observed). For pKa, you can use PropKa (https://github.com/jensengroup/propka-3.1), and then I guess you'll be able to get the pI somewhere around the code. I've been designing a GUI to use this in Chimera, but it's not yet ready... Sorry about that! I don't know if you need anything else... but I think the Python community will have you covered one way or another! BioPython.PDB might be useful too (http://biopython.org/wiki/The_Biopython_Structural_Bioinformatics_FAQ). Cheers, Jaime. 2016-10-06 15:57 GMT+02:00 Healey, Joe <J.R.J.Healey(a)warwick.ac.uk<mailto:J.R.J.Healey@warwick.ac.uk>>: Hi again Jaime, Hope you don't mind but I had a couple of questions for you. You may have figured out from all the previous correspondance that I have a load of protein structure simulations. The RMSD script you helped me with was to sort of 'broadly' gauge whether the simulations looked realistic or not. For some of them, I'd like to build up a chart of their tertiary structure properties based on the simulations. Things like charge/hydrophobicity, isoelectric point, 2ndary structure dominance and so on. Could I ask you for any advice you might about what tools are available to analyse PDBs? I'm sure chimera is probably capable of it too. I've mostly found webservers at the moment for specific tasks, but I'd like to do most of the analysis locally. My group is full of bioinformaticians, but I'm the only one who really looks at proteins so I need the advice of a computational chemist I think! Cheers, Joe Joe Healey M.Sc. B.Sc. (Hons) PhD Student MOAC CDT, Senate House University of Warwick Coventry CV47AL Mob: +44 (0) 7536 042620<tel:%2B44%20%280%29%207536%20042620> | Email: J.R.J.Healey(a)warwick.ac.uk<mailto:J.R.J.Healey@warwick.ac.uk> Jointly working in: Waterfield Lab<http://www2.warwick.ac.uk/fac/med/research/tsm/microinfect/staff/waterfield…> (WMS Microbiology and Infection Unit) and the Gibson Lab<http://www2.warwick.ac.uk/fac/sci/chemistry/research/gibson/gibsongroup/> (Warwick Chemistry) Twitter: @JRJHealey<https://twitter.com/JRJHealey> | Website: MOAC Page<http://www2.warwick.ac.uk/fac/sci/moac/people/students/2013/joseph_healey> ________________________________ From: Healey, Joe Sent: 23 September 2016 11:51:40 To: Jaime Rodríguez-Guerra Cc: Eric Pettersen; chimera-users(a)cgl.ucsf.edu<mailto:chimera-users@cgl.ucsf.edu> Subject: Re: [Chimera-users] Recursive structure matching and acquisition of descriptive numbers Not at all! Much appreciated, as I think I mentioned, this is my first go at writing anything remotely complex in python (and not a bad turn out for about a week I don't think...) A lot of the mess comes from copying functional bits from StackOverflow in the interests of just getting it working first, and making it pretty second as I'm sure you can appreciate (you'd hate to see my LaTeX preamble)! Not needed to close the file with 'with' is useful to know (takes some of the headache out of knowing where to place the close). I'll certainly make the tidy ups you've suggested, for my own learning. Thanks again, Joe :) Joe Healey M.Sc. B.Sc. (Hons) PhD Student MOAC CDT, Senate House University of Warwick Coventry CV47AL Mob: +44 (0) 7536 042620<tel:%2B44%20%280%29%207536%20042620> | Email: J.R.J.Healey(a)warwick.ac.uk<mailto:J.R.J.Healey@warwick.ac.uk> Jointly working in: Waterfield Lab<http://www2.warwick.ac.uk/fac/med/research/tsm/microinfect/staff/waterfield…> (WMS Microbiology and Infection Unit) and the Gibson Lab<http://www2.warwick.ac.uk/fac/sci/chemistry/research/gibson/gibsongroup/> (Warwick Chemistry) Twitter: @JRJHealey<https://twitter.com/JRJHealey> | Website: MOAC Page<http://www2.warwick.ac.uk/fac/sci/moac/people/students/2013/joseph_healey> ________________________________ From: Jaime Rodríguez-Guerra <jaime.rodriguezguerra(a)uab.cat<mailto:jaime.rodriguezguerra@uab.cat>> Sent: 23 September 2016 10:56:38 To: Healey, Joe Cc: Eric Pettersen; Jaime Rodríguez-Guerra; chimera-users(a)cgl.ucsf.edu<mailto:chimera-users@cgl.ucsf.edu> Subject: Re: [Chimera-users] Recursive structure matching and acquisition of descriptive numbers Hi Joe! You're welcome! Glad to see you got it working. I've reviewed your code and would want to point out some details that might be helpful to you in the future, though. These are mostly style-related things, but those will contribute to a cleaner reading experience! 1. You don't need backslashes if you're inside parenthesis. So, in all those parse_argument lines, you can safely delete them. Related to this, you can wrap imports with parenthesis, so that backslashes are not needed. 2. Tuple unpacking does not need parenthesis in you're dealing with 1D-tuples. Only needed if you are dealing with more than one dimension (ie, (animal, fruit), sport = [['tiger', 'pear'], 'soccer'] 3. Since you are opening files with the "with" context manager, you don't need to manually close the file later. It's closed automatically as soon as you leave that with block. 4. Lines 166-170 could be replaced with a (cleaner, but not necessarily more performant) glob.glob() call. Check it out to see if it satisfies your requirements. 5. While you are at it, take a look at the PEP8 docs and the Google Python style guide. You don't need to follow all the rules, but they are really helpful in bringing consistency to your own style! That's it! Hope you don't mind these pieces of advice :) Cheers, Jaime. 2016-09-23 11:27 GMT+02:00 Healey, Joe <J.R.J.Healey(a)warwick.ac.uk<mailto:J.R.J.Healey@warwick.ac.uk>>: > That was exactly what I was after thanks! I figured the object orientation > should provide that somewhere but was barking up slightly the wrong tree! > > > Thank you all for all your help - definitely wouldn't have been able to do > it without you! > > > If it's of any interest, or it helps for any future questions where you > might want to refer back to code segments, I've put the > more-or-less-finished script in this paste: > > > http://pastebin.com/xJNSWJGq > > > Joe Healey > > > M.Sc. B.Sc. (Hons) > PhD Student > MOAC CDT, Senate House > University of Warwick > Coventry > CV47AL > Mob: +44 (0) 7536 042620<tel:%2B44%20%280%29%207536%20042620> | Email: J.R.J.Healey(a)warwick.ac.uk<mailto:J.R.J.Healey@warwick.ac.uk> > > Jointly working in: > Waterfield Lab (WMS Microbiology and Infection Unit) > and the Gibson Lab (Warwick Chemistry) > > Twitter: @JRJHealey | Website: MOAC Page > ________________________________ > From: Eric Pettersen <pett(a)cgl.ucsf.edu<mailto:pett@cgl.ucsf.edu>> > Sent: 22 September 2016 20:54:35 > To: Jaime Rodríguez-Guerra > Cc: chimera-users(a)cgl.ucsf.edu<mailto:chimera-users@cgl.ucsf.edu>; Healey, Joe > Subject: Re: [Chimera-users] Recursive structure matching and acquisition of > descriptive numbers > > Thanks Jaime — exactly right. I realized that the “atoms1, atoms2” in my > original example didn’t really give any indication of which atoms were > which, which is why in my later examples I changed to "simAtoms, refAtoms”, > so it’s actually “atoms2” that are the reference atoms. > > —Eric > >> On Sep 22, 2016, at 11:45 AM, Jaime Rodríguez-Guerra >> <jaime.rodriguezguerra(a)uab.cat<mailto:jaime.rodriguezguerra@uab.cat>> wrote: >> >> Hi! >> >> Since you are getting lists of atoms back, you only need to do some >> attribute access (rather than method calling)! >> >> Every chimera.Atom object has an attribute called molecule: a >> reference to its parent chimera.Molecule. chimera.Molecule objects >> have an attribute called 'name', but sometimes this is not very >> informative (depends on your input), so you need to resort to the pdb >> filename, stored in the first element of the openedAs tuple. >> >> In code, this looks like this: >> >> for atoms1, atoms2, rmsd, fullRmsd in match(CP_BEST, [ref, sims], >> defaults[MATRIX], "nw", defaults[GAP_OPEN], defaults[GAP_EXTEND]): >> molecule1 = atoms1[0].molecule # any atom from the list will do >> molecule2 = atoms2[0].molecule >> print molecule1.name<http://molecule1.name>, "\t", molecule2.name<http://molecule2.name>, "\t", rmsd # >> molecule1.openedAs[0] will also work here >> >> I don't know if atoms1 comes consistently from the reference molecule, >> but I'd guess it does. >> >> Hope it helps! >> _______________________________________________ >> Chimera-users mailing list: Chimera-users(a)cgl.ucsf.edu<mailto:Chimera-users@cgl.ucsf.edu> >> Manage subscription: >> http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users >

3 7

Chimera using problem
by Chonghua Ren 11 Oct '16

11 Oct '16

Dear there, I am learning how to use Chimera now, I like it so much and find it is so awesome. I have one question, if I want to input a DNA sequence like ID "1d86", how can I transfer the DNA sequence to PDB format? Another question, how can I know protein and DNA interaction using this software? Thank you so much! Best wishes, Chonghua Ren

2 1

Chimera-users October 2016