June 2016 - Chimera-users - RBVI Mailing Lists

Cap color precedence when global clip planes slice through multiple volumes?
by Oliver Clarke 06 Jun '16

06 Jun '16

Hi, When slicing through multiple superimposed volumes with caps shown, the color of the resulting cap seems to be a bit of a patchwork - sometimes dominated by one volume, sometimes by another, sometimes a mix of the two It would be good if there were a way to control the precedence of the caps, to make a cleaner and more interpretable picture - it is often good to be able to see the outline of both overlayed volumes. I know there is one way to do this, by using per model clipping planes and changing the offsets such that one per model clip plane is just slightly in front of the other, but it might be helpful in ChimeraX to have more granular control over this. Cheers, Oli

2 2

Cut asymmetric unit from volume
by Morris, Kyle 03 Jun '16

03 Jun '16

Hi Chimera team and users, Is there a way to extract the asymmetric unit from a volume which has been 'symmetrised’ from, say for instance, Relion? In this case what are the conventions used by chimera to define the bounds of the sub-volume to extract? Thanks for your help! Best wishes, Kyle

2 3

Request for support
by Mendoza Wilson, Ana Maria 03 Jun '16

03 Jun '16

Hello, How I can select an isodensity value in the Chimera software (1.11rc) for calculating molecular surface and molecular volume from one cube data file or fchk file previously generated in Gaussian 09? Thanks in advance -- Dra. Ana María Mendoza Wilson Investigador Titular Centro de Investigación en Alimentación y Desarrollo, A.C. Coordinación de Tecnología de Alimentos de Origen Vegetal Carretera a La Victoria Km. 0.6 Hermosillo, Sonora, México C.P. 83000 Tel +52 (662) 289 24 00 Ext. 354

2 3

FW: Chimera Log File
by Williams, Joanne 03 Jun '16

03 Jun '16

From: Erin Coppola <erincoppola8(a)gmail.com<mailto:erincoppola8@gmail.com>> Date: Friday, June 3, 2016 8:10 AM To: "chimera(a)cgl.ucsf.edu<mailto:chimera@cgl.ucsf.edu>" <chimera(a)cgl.ucsf.edu<mailto:chimera@cgl.ucsf.edu>> Subject: Chimera Log File Hi, When using Chimera, is there a way of tracking the script command versions of GUI actions performed?

2 1

Creating maps from models for map/model FSC calculations
by Ardan Patwardhan 03 Jun '16

03 Jun '16

Dear all We are using Chimera to align PDB models to maps, then creating maps from the models using the following command (as an example): molmap #0 5 gridSpacing 1.3 This should create a map with the voxel spacing of 1.3A and a gaussian atom spread function with the width related to 5A resolution (https://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/molmap.html). Our problem is a) the map that is created can have an arbitrary size - we need it to be cubic and have the same (known) dimensions as the second map against which it will be compared b) we have tried padding it in other programs like bsoft but then we sometimes run into problems of centring the model. I was wondering if there was a way to get a centred map of a certain size from a model directly from Chimera? Perhaps there are certain parameters that can be used in molmap? You thoughts and advice on this issue are much appreciated! Many thanks and best wishes Ardan Patwardhan Coordinator - cellular structure & PDBe production Protein Data Bank in Europe (PDBe) European Bioinformatics Institute (EMBL-EBI) European Molecular Biology Laboratory Wellcome Trust Genome Campus Hinxton, Cambridge CB10 1SD Tel: +44 1223 492649

3 7

mask
by Dougherty, Matthew T 02 Jun '16

02 Jun '16

I have a volume whose density values are -127 thru 127. Zero is a valid density & approximately the peak of the Guassian distribution. I would like to use the command mask. It sets the unmasked density values to zero. How can I set the unmasked density values to -128 using the mask command? Will I need to hack the code? Matthew Dougherty National Center for Macromolecular Imaging Baylor College of Medicine ================================================= =================================================

2 1

finding residues with large effects on local electrostatic potential
by Elaine Meng 01 Jun '16

01 Jun '16

Hi Deniz, No need to apologize, that’s why we have an address for questions. :-) I am guessing you were coloring by electrostatic potential and want to figure out which residues are having large effects in some places. The question is not that simple because electrostatic interactions are long-range, meaning that the potential at a point is based on many atoms in many residues, including those farther away. However, the largest contributions at a surface point will generally be from the residues with net positive or negative charges that are not too far away. In other words, if you see a blue patch, you could then try to figure out positively charged residues are near it, or if a red patch, which negatively charged residues are near it. The problem in Chimera that you may be having is that when you put the mouse over the surface, the surface is blocking the atoms, so it does not tell you what residue is under that surface. One way to try to figure it out is to temporarily hide the surface (Favorites… Model Panel, uncheck the “Shown” button for the surface), then click back into the main window *without moving the structure* and put the mouse over the same area to see what residues are there, using the pop-up balloons with residue information. You would probably want to show all the sidechain atoms first (e.g. command: display protein), and maybe select residues that are known to be positive (or negative) to help guide yourself to the right place. Command: select :lys,arg,his (… the ones likely to contribute to positive patches…) Command: select :asp,glu (… the ones likely to contribute to negative patches…) You can keep turning the surface on and off with the checkbox in the Model panel. This is interactive guessing and may not be that easy. Just keep in mind that you can’t say that some patch is all from one or two residues, only that they have significant contributions. 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 Jun 1, 2016, at 3:44 PM, Deniz Simsek Buck <simsekd(a)stanford.edu> wrote: > Hi Elaine, > I am sorry for bothering you. But, I need urgent help. > If I tell chimera to find the surface charge, negative and positive, and if I want to know which residues are responsible for this charge, what do I do? > > thank you so much, > deniz

1 0

How to build symmetry in Chimera
by Ana Ramos 01 Jun '16

01 Jun '16

To whom it may concern, I am working with a protein with P61 symmetry, and I would like to build symmetry mates; however in the pdb, the symmetry is not stated. How can I build in chimera symmetry mates with a P61 symmetry? Best, Ana

3 2

Chimera 1.11 - Combining Transparancies of Ribbon and Surface
by Falko Jähnert 01 Jun '16

01 Jun '16

To whom it may concern, since I'm using Chimera a lot lately I enjoy the usage of this program and can't await the release of Chimera X. Now I like to ask if it is possible to combine the transparency of the ribbon and the surface representation using Chimera so that you can visualize both parts in varying transparency levels. My problem with this is if I use "transparency 80,r [Tag]" the part of my protein gets transparent (as requested) whilst activating the transparency function of the surface results in complete hiding of the previously chosen amino acids. Is there a workaround available? Thanks in advance for your help. Best regards, Falko Jähnert. _____ Martin-Luther-University Halle-Wittenberg Institute for Biochemistry and Biotechnology Departement for Naturstoffbiochemie Kurt-Mothes-Str. 3 D-06110 Halle/Saale _____ Tel. +49 (0)345 - 55 24812 Mail. <mailto:falko.jaehnert@biochemtech.uni-halle.de> falko.jaehnert(a)biochemtech.uni-halle.de

2 1

Re: [Chimera-users] Analysis of the contact maps from the MARTINI MD trajectory
by James Starlight 01 Jun '16

01 Jun '16

Yep thanks so much Eric! In my case the contacts have been started to be recognized both from only after choosing the vdw cutoff > -4.0 A and decreasing statistical (weight discard edge threshold) factor to 0.3 (here I am analyzing big trajectory consisted of 7 independent MD trajectories for the same system merged together where both proteins are present both in bound and unbound states during the progression of simulation). Eventually I obtained LYS 25.A vdw SER 712.C 0.39322 CYS 14.A vdw TRP 1586.I 0.427119 HIS 33.A vdw ions 0.335593 HIS 33.A vdw LYS 1577.I 0.444068 LYS 27.A vdw TRP 1586.I 0.369492 GLU 838.C vdw PHE 46.A 0.444068 GLN 12.A vdw TRP 1586.I 0.501695 ASP 1580.I vdw LYS 7.A 0.318644 PHE 46.A vdw TRP 1586.I 0.491525 GLU 732.C vdw LYS 79.A 0.372881 PHE 46.A vdw solvent 0.501695 CYS 1584.I vdw HIS 26.A 0.616949 HIS 26.A vdw SER 712.C 0.423729 ASP 1580.I vdw LYS 8.A 0.305085 TYR 1542.I vdw VAL 11.A 0.311864 LYS 27.A vdw solvent 0.40678 THR 28.A vdw TYR 731.C 0.379661 PHE 46.A vdw TYR 1587.I 0.484746 Also may I ask within this output the first column should correspond to first selection, right? here the contacts was defined as between atoms from chain A which is the smaller protein and the rest of the system consisted of atoms of receptor. Thanks again! Gleb 2016-05-17 22:37 GMT+02:00 Eric Pettersen <pett(a)cgl.ucsf.edu>: > Well, if I were you I would go to a trajectory frame that I think should have contacts, and then bring up the Find Clashes/Contacts tool (in the Structure Analysis category of Tools) and play around with the criteria to see what yields contacts. You might also change the atoms to sphere mode (Actions->Atoms/Bonds->sphere) to see if they have any VDW overlap or are close. You might also measure some distances (control-click an atom, then control-shift-double-click a second atom and then pick “Show Distance” from the popup menu). Without actual access to your data, I can’t offer any more specific suggestions. > > —Eric > >> On May 17, 2016, at 2:08 AM, James Starlight <jmsstarlight(a)gmail.com> wrote: >> >> Thanks Eric! >> >> it quite strange that no interactions are found- because in my >> trajectory the complex between both proteins are emerged very soon >> after beginning of the simulation and is quite stable until ed of it. >> Seems like the Chimera does not recognize martini's CG atoms properly. >> Are there any suggestions to fix it e.g by the variations using any of >> options - e.g vdw settings etc? >> >> Gleb >> >> 2016-05-17 1:32 GMT+02:00 Eric Pettersen <pett(a)cgl.ucsf.edu>: >>> The clustering is based on least-squares-fit RMSDs between corresponding atoms, and therefore global translations and rotations are irrelevant. >>> >>> —Eric >>> >>>> On May 16, 2016, at 9:08 AM, James Starlight <jmsstarlight(a)gmail.com> wrote: >>>> >>>> update :-) >>>> >>>> in addition to the previous question which is still actual I am >>>> interesting in the clusterization of binding interfaces established >>>> during MD e.g based on the position of one protein regarding another. >>>> >>>> For my case I have performed such clusterization from MD movie plugin >>>> (which was based on current selection which was the atoms of smaller >>>> protein) obtaining alot of overlapped clusters in case of smallest >>>> step size and afew clusters increasing step size twisly. Is it >>>> possible here to increase accuracy of the clusterization and to remove >>>> some degrees of freedom which are not important for my case ? E.g to >>>> pre-process trajectory using PCA before clusterization to discard >>>> rotation translation of the system will be good option? >>>> >>>> >>>> >>>> James >>>> >>>> 2016-05-16 15:20 GMT+02:00 James Starlight <jmsstarlight(a)gmail.com>: >>>>> update : >>>>> >>>>> already playing with the plugin MD movie and "residue iteration >>>>> network within the trajectory" tool >>>>> just selecting atoms from the chain A as 1st selection >>>>> and other atoms from other chains as 2nd selection >>>>> >>>>> running calculation (trying to varry cut-off corresponded to vdw overlap radii) >>>>> >>>>> obtained error >>>>> AttributeError: 'ResInteractionDialog' object has no attribute 'markers' >>>>> >>>>> File "/opt/UCSF/Chimera64-1.10.2/share/Movie/ResInteraction.py", >>>>> line 535, in Apply >>>>> markerSettings = [(m['xy'], m['rgba']) for m in self.markers] >>>>> >>>>> See reply log for Python traceback. >>>>> >>>>> Does it means that somethig special should be specified for the >>>>> MARTINI CG atoms? >>>>> >>>>> Thanks ! >>>>> >>>>> Gleb >>>>> >>>>> 2016-05-16 12:13 GMT+02:00 James Starlight <jmsstarlight(a)gmail.com>: >>>>>> something relevant in my model: >>>>>> >>>>>> this is martini CG model consisted of big membrane protein and small >>>>>> water soluble protein. In my particular trajectory I have only atoms >>>>>> of both proteins recognized by chimera as individual chains (e.g small >>>>>> protein is chain A and the biggest is consisted of chains B-N). >>>>>> >>>>>> Gleb >>>>>> >>>>>> 2016-05-16 11:58 GMT+02:00 James Starlight <jmsstarlight(a)gmail.com>: >>>>>>> update: >>>>>>> >>>>>>> with the up-to-date chimera Trajectory is loaded fine >>>>>>> now the question is related to the analysis of the binding interface. >>>>>>> is it possible to make contact maps plots based on the 2 selections >>>>>>> from the input trajectory assuming that I am using martini CG model? >>>>>>> >>>>>>> Gleb >>>>>>> >>>>>>> 2016-05-16 10:39 GMT+02:00 James Starlight <jmsstarlight(a)gmail.com>: >>>>>>>> OK will try to update Chimera today ! >>>>>>>> >>>>>>>> What chimera's tools could be useful for the contact map based >>>>>>>> analysis of md trajectories related to my case ? E.g to determine >>>>>>>> binding interface between 2 proteins from long md simulation. >>>>>>>> >>>>>>>> Gleb >>>>>>>> >>>>>>>> 2016-05-13 19:52 GMT+02:00 Eric Pettersen <pett(a)cgl.ucsf.edu>: >>>>>>>>> Hi Gleb, >>>>>>>>> If you are using Gromacs 5 trajectories, you have to be using at least >>>>>>>>> version 1.10.2 of Chimera. If you are using that version (or later), please >>>>>>>>> use “Report a Bug” in Chimera’s Help menu to file a bug report and attach >>>>>>>>> your topology file to the bug report. Thanks! >>>>>>>>> >>>>>>>>> —Eric >>>>>>>>> >>>>>>>>> Eric Pettersen >>>>>>>>> UCSF Computer Graphics Lab >>>>>>>>> >>>>>>>>> >>>>>>>>> On May 13, 2016, at 1:46 AM, James Starlight <jmsstarlight(a)gmail.com> wrote: >>>>>>>>> >>>>>>>>> Dear Chimera users! >>>>>>>>> >>>>>>>>> >>>>>>>>> I am in charge with the analysis of protein-protein association during >>>>>>>>> long molecular dynamic simulation where my system was parametrized >>>>>>>>> using MARTINI CG force field. In particularly I am interesting to >>>>>>>>> find residues on one of the protein which are crustal for the binding >>>>>>>>> interface established during this MD. >>>>>>>>> For that purpose I am trying to use Chimera to load trajectory and >>>>>>>>> corresponded tpr file using MD movie plugin and than to >>>>>>>>> map contact maps produced by Gromacs onto the 3D structure using Chimera. >>>>>>>>> The problem that Chimera does not recognize properly the trajectory >>>>>>>>> and topology. Briefly I have removed all solvent from both files >>>>>>>>> before loading them to the Chimera using editconf and gmx convert-tpr >>>>>>>>> obtaining eventually trajectory and topology with the same number of >>>>>>>>> atoms. >>>>>>>>> >>>>>>>>> >>>>>>>>> Than when I load it to Chimera that is the error which MD movie sent me: >>>>>>>>> >>>>>>>>> VERSION 5.0.2 >>>>>>>>> using floats >>>>>>>>> version 100, generation 26 >>>>>>>>> 8 atoms >>>>>>>>> Traceback (most recent call last): >>>>>>>>> File "CHIMERA/lib/python2.5/site-packages/Pmw/Pmw_1_3/lib/PmwBase.py", >>>>>>>>> line 1747, in __call__ >>>>>>>>> File "CHIMERA/share/chimera/baseDialog.py", line 328, in command >>>>>>>>> File "CHIMERA/share/chimera/baseDialog.py", line 543, in OK >>>>>>>>> File "CHIMERA/share/Trajectory/EnsembleLoader.py", line 92, in Apply >>>>>>>>> File "CHIMERA/share/Trajectory/formats/Gromacs/__init__.py", line >>>>>>>>> 56, in loadEnsemble >>>>>>>>> File "CHIMERA/share/Trajectory/formats/Gromacs/__init__.py", line >>>>>>>>> 67, in loadEnsemble >>>>>>>>> File "CHIMERA/share/Trajectory/formats/Gromacs/Gromacs.py", line 25, >>>>>>>>> in __init__ >>>>>>>>> File "CHIMERA/share/Trajectory/formats/Gromacs/Gromacs.py", line 86, >>>>>>>>> in __init__ >>>>>>>>> File "CHIMERA/share/Trajectory/formats/Gromacs/Gromacs.py", line >>>>>>>>> 345, in _readTopology >>>>>>>>> File "CHIMERA/share/Trajectory/formats/Gromacs/Gromacs.py", line >>>>>>>>> 338, in _readString >>>>>>>>> File "CHIMERA/lib/python2.5/xdrlib.py", line 197, in unpack_fstring >>>>>>>>> raise EOFError >>>>>>>>> EOFError >>>>>>>>> EOFError >>>>>>>>> >>>>>>>>> File "CHIMERA/lib/python2.5/xdrlib.py", line 197, in unpack_fstring >>>>>>>>> raise EOFError >>>>>>>>> >>>>>>>>> See reply log for Python traceback. >>>>>>>>> >>>>>>>>> Will be thankful for any suggestions! >>>>>>>>> >>>>>>>>> Gleb >>>>>>>>> _______________________________________________ >>>>>>>>> Chimera-users mailing list: Chimera-users(a)cgl.ucsf.edu >>>>>>>>> Manage subscription: >>>>>>>>> http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users >>>>>>>>> >>>>>>>>> >>>> >>>> _______________________________________________ >>>> Chimera-users mailing list: Chimera-users(a)cgl.ucsf.edu >>>> Manage subscription: http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users >>>> >>> >> >

3 14

Chimera-users June 2016