- Chimera-users - RBVI Mailing Lists

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

suggest improvements/group viz
by Dougherty, Matthew T 31 May '16

31 May '16

1) in a group viz it would be desirable to have two cursors. A regular cursor that would operate as expected under the OS, and a second cursor that would be restricted to the view window overlay. The first cursor would be controlled by the viz operator used to adjust parameters. The second cursor would be controlled by the group through a gamepad etc, that could be used as a pointer. Possibly the volume eraser might be a quick work around. 2) volume eraser when erasing sets the density to zero. Need to be able to set the value. Zero is a valid density value, usually in the middle of the density histogram, which causes problems for volume rendering. Matthew Dougherty National Center for Macromolecular Imaging Baylor College of Medicine ================================================= =================================================

1 0

Rotamers in Chimera
by Lew Jacobson 31 May '16

31 May '16

I am using Chmiera 1.10 build 40304 on a Mac. (I am a leftover Midas user!) So far as I can tell, this version is still using the Dunbrack 2002 rotamer library. Is any subsequent version using the improved 2011 Dunbrack library? If not, can you instruct me how to incorporate that newer library into Chimera? thanks, Lew Jacobson -- Dr. Lewis Jacobson Professor of Biological Sciences Univ. of Pittsburgh Pittsburgh PA 15260 412-624-4647 Fax 412-624-4759 LJAC(a)pitt.edu

2 1

Can I use wait command over Rest service?
by Mungo Carstairs (Staff) 31 May '16

31 May '16

Hi, I am trying to devise a way to highlight locations on a structure in Chimera while mousing over the sequence in Jalview. My latest attempt tries to 'flash' the side-chain display at the residue position i.e. - hide then show atoms (if currently displayed e.g. display is in atom mode) - or show then hide atoms (if currently not displayed e.g. display is in ribbon mode) I'm sending a command like display #0:30.A;wait 20;~display #0:30.A to the Rest service. However while this seems to do what I want when used at the Chimera command line, via the Rest service it ignores the wait command, with the result that the action is too quick to see. Any suggestions please for how to get this to work, or a better alternative? Thanks, Mungo Mungo Carstairs Jalview Computational Scientist The Barton Group Division of Computational Biology School of Life Sciences University of Dundee, Dundee, Scotland, UK. www.jalview.org<http://www.jalview.org/> www.compbio.dundee.ac.uk<http://www.compbio.dundee.ac.uk/> The University of Dundee is a registered Scottish Charity, No: SC015096

3 8

broken lines
by shazia sabir 29 May '16

29 May '16

Hi dear All, Why chimera is showing broken lines (region encircled in red at one area). What can we predict from structure. I attached sequence of structure that shows amino acid is there(gey region in sequence). Regards, shazia.

2 1

Modeller awareness of other proteins
by Eiros Zamora, Juan 26 May '16

26 May '16

Dear Chimera experts, I’ve got a system of three interacting proteins, and each of them is missing some linker regions in the crystal structure. I’ve succesfully used the Modeller interface in Chimera to add these regions, but the problem I’m facing now is that it adds the residues without considering that other atoms from another protein of the system are in the way. Is there a way to tell Modeller to be sensitive to other structures nearby? Thanks for any comments, [Imperial College London] <http://www.icb-cdt.co.uk/> Juan Eiros Zamora / PhD Student je714(a)ic.ac.uk<mailto:je714@ic.ac.uk> / +44 (0) 7783018326 Imperial College London Imperial College Road SW7 2AZ | UK http://www.icb-cdt.co.uk/

2 2

Evolutional-based analysis in Chimera
by James Starlight 25 May '16

25 May '16

Dear Chimera Users! For my particular analysis I have 2 homologue membrane proteins (one is bacterial, and another one- mammalian cytocrome-C oxydase) each bound to same smaller (ligand) protein called for simplicity "the_guy". For that particular case the conservation in sequence of both protein is quite low however 1) the overall fold of both proteins especially within the membrane-embedded region as well as 2) ligand-bound (which is mostly the water exposed loops) region are quite similar. Now using some Chimera tool I need to estimate binding interfaces of two homologous proteins in terms of the conservation of the resides involved in binding of the_guy. Then I need to make more complex task- performing the same evolutional analysis but comparing 2 ensembles extracted from the MD trajectory represent simulation of binding of each of the protein with the_guy to track correlations of the statistics between the residues involved in the binding during MD as well as overall conservation for particular fold. Does the "Evolutionary tracing" might be useful here? Thanks for help! Gleb

2 1

PDB post-processing
by James Starlight 25 May '16

25 May '16

Dear Chimera users! After some post-processing of MD simulation I have extracted from the trajectory several snapshots as individual pdbs which I need to 1- to re-assign information regarding chain letters (which was lost after some operations on pdbs) for all of those structures assuming that within each of these pdbs residues are ordered correctly e.g from (1-104 which should be chain A) than from 1-100 (it should be chain B), than 1-70 (in should be chain C) etc. 2- superimpose and combine all of the pdbs within one model in nmr-like format. Thanks for help!

2 1

Long-ranged electrostatics analysis
by James Starlight 25 May '16

25 May '16

Dear Chimera users! I am analysing several snapshots etracted from the full atomistic md simulation of protein-protein assosiation. I wounder to ask what are the Chiemra plugins might be useful to study long-rang electrostatic effect mediated initial recognition of both proteins in this system besides APBS caclucation which is more focused on short-ranged electrostatics. Thanks for help! Gleb

2 1

Script to open density from Coot in Chimera.
by Oliver Clarke 23 May '16

23 May '16

Hi all, Posting this in case it is useful to anyone else who uses both programs. This script will allow Coot to open the currently displayed molecules/maps in Chimera, and adjust display settings appropriately for taking a picture of electron density. Either put the appended code in your coot.py and use from the scripting interface, or place this longer script ( https://www.dropbox.com/s/0b4bebwxw0p9x0e/oli_custom.py?dl=0) in your ~/.coot-preferences directory, which will add a new "Custom" menu including this option under the "Display" submenu (as well as a bunch of other stuff, and some extra custom key bindings). The orientation of the molecules should be correct after loading in Chimera, but the zoom and clipping may require some adjustment (I am not sure exactly how to match these precisely between the two programs). The script will only load a fragment of the maps (the same region that is displayed in Coot), but this has the advantage that it takes advantage of Coot's ability to deal with crystal symmetry, so it can be used at crystal contacts with no problems. Let me know if you run into any bugs! Cheers, Oli. Code: #Open displayed models and maps in Chimera, set map levels and view def open_in_chimera(): if find_exe("chimera"): #If chimera exists do stuff, else raise an error import subprocess pwd=os.getcwd() #Dir from which coot was launched cofr=rotation_centre() #3-membered list [x,y,z] make_directory_maybe("coot-chimera") #Put coot droppings here coot_chimera_path=pwd+"/coot-chimera/" check_path=coot_chimera_path+"chimera_launcher.cmd" #Chimera run script that will be written later check_path2=coot_chimera_path+"matrix.txt" #Orientation matrix for chimera input view_number=add_view_here("tmp") # initial_zoom=zoom_factor() #Coot's zoom factor. 100 is the initial state, zooming out gets larger, in gets smaller. # reset_view() zoom=zoom_factor() # go_to_view_number(view_number,1) # scale_factor=initial_zoom/base_zoom # chimera_scale=base_zoom/initial_zoom chimera_scale=100.0/zoom #Approx - breaks down when initial views upon loading mol in Coot or Chimera are different map_radius=get_map_radius() #We'll use this to export a map_radius sized fragment if os.path.isfile(check_path): #if chimera script exists, delete it. os.remove(check_path) if os.path.isfile(check_path2): #same for matrix.txt os.remove(check_path2) map_list=map_molecule_list() mol_list=model_molecule_list() matrix_list=view_matrix() #9-membered list of rotation matrix elements. with open(check_path2,"a") as matrix_file: #Write orientation matrix (specify first pdb as model, then use matrixcopy to apply to all others). Translation vector in last column is 0,0,0 because we will take care of that separately. matrix_file.write("Model {model_id}.0\n".format(model_id=mol_list[0])) matrix_file.write("\t {a} {b} {c} 0.0\n".format(a=matrix_list[0],b=matrix_list[1],c=matrix_list[2])) matrix_file.write("\t {a} {b} {c} 0.0\n".format(a=matrix_list[3],b=matrix_list[4],c=matrix_list[5])) matrix_file.write("\t {a} {b} {c} 0.0\n".format(a=matrix_list[6],b=matrix_list[7],c=matrix_list[8])) with open(check_path,"a") as cmd_file: #Start writing stuff to chimera launch script. for mol_id in model_molecule_list(): if mol_is_displayed(mol_id): file_name=coot_chimera_path+"mol_{mol_id}.pdb".format(mol_id=mol_id) write_pdb_file(mol_id,file_name) path_to_file=file_name model_id=mol_id cmd_file.write("open #{model_id} {mol}\n".format(model_id=model_id,mol=path_to_file)) cmd_file.write("color gold #{model_id}; color byhet #{model_id}; ca_and_sidechains #{model_id}; disp ~protein\n".format(model_id=model_id)) cmd_file.write("matrixset matrix.txt\n") cmd_file.write("matrixcopy #{mol_id0} #{mol_id}\n".format(mol_id0=mol_list[0],mol_id=model_id)) map_id0=mol_id for map_id in map_molecule_list(): if map_is_displayed(map_id): map_level=get_contour_level_absolute(map_id) if map_is_difference_map(map_id): model_id=map_id0+map_id #make sure that the assigned model number of each map is unique and does not overlap with those of pdbs. if map_id==0: model_id=model_id+1 file_name=coot_chimera_path+"diff_map_{model_id}.mrc".format(model_id=model_id) path_to_file=file_name export_map_fragment(map_id,cofr[0],cofr[1],cofr[2],map_radius,file_name) cmd_file.write("open #{model_id} {diff_map} \n".format(model_id=model_id,diff_map=path_to_file)) cmd_file.write("volume #{model_id} capfaces false style mesh meshlighting false squaremesh false color \"#08882eefa222\" step 1 ; sop cap off ; set depthCue ; set dcStart 0.2 ; set dcEnd 1 ; background solid white ; set showcofr ; cofr view ; clip on; volume #{model_id} level -{map_level} color #da1200000000 level {map_level} color #0000bda00000 \n".format(model_id=model_id,map_level=map_level)) cmd_file.write("matrixset matrix.txt\n") cmd_file.write("matrixcopy #{mol_id0} #{model_id}\n".format(mol_id0=mol_list[0],model_id=model_id)) else: model_id=map_id0+map_id if map_id==0: model_id=model_id+1 file_name=coot_chimera_path+"map_{model_id}.mrc".format(model_id=model_id) export_map_fragment(map_id,cofr[0],cofr[1],cofr[2],map_radius,file_name) path_to_file=file_name cmd_file.write("open #{model_id} {map} \n".format(model_id=model_id,map=path_to_file)) cmd_file.write("volume #{model_id} capfaces false style mesh meshlighting false squaremesh false color \"#08882eefa222\" step 1 ; sop cap off ; set depthCue ; set dcStart 0.2 ; set dcEnd 1 ; background solid white ; set showcofr ; cofr view ; clip on; volume #{model_id} level {map_level} \n".format(model_id=model_id,map_level=map_level)) cmd_file.write("matrixset matrix.txt\n") cmd_file.write("matrixcopy #{mol_id0} #{model_id}\n".format(mol_id0=mol_list[0],model_id=model_id)) cmd_file.write("~sel; cofr {x},{y},{z} coordinatesystem #{mol_id0}; ac mc; center sel; cofr sel; clip hither 5 fromCenter true; clip yon -5 fromCenter true; cofr view; ~set showcofr; del sel; scale {chimera_scale}".format(x=cofr[0],y=cofr[1],z=cofr[2],mol_id0=mol_list[0],chimera_scale=chimera_scale)) chimera_exe=find_exe("chimera") info_dialog("Opening in Chimera...\nOrientation should be right but you will probably need to \nadjust the scale and clipping to get the same view as in Coot. \n If you can't see anything, try zooming out.") subprocess.Popen(["chimera",check_path]) else: info_dialog("Sorry, you need UCSF Chimera installed and accessible from the terminal for this to work!") #This works okay, though adjustment of zoom and clipping still not ideal. #would it be possible to add a subprocess.communicate() to use Coot as a controller for Chimera?

1 0

Chimera does not start with mesa-11.2.2
by Jaime Rodríguez-Guerra Pedregal 23 May '16

23 May '16

Hi guys, just to let you know that mesa-11.2.2 is messing with Chimera startup sequence.: $ chimera --debug initializing general preferences loading Tix initializing graphics libGL error: unable to load driver: i965_dri.so libGL error: driver pointer missing libGL error: failed to load driver: i965 libGL error: unable to load driver: swrast_dri.so libGL error: failed to load driver: swrast X Error of failed request: BadValue (integer parameter out of range for operation) Major opcode of failed request: 154 (GLX) Minor opcode of failed request: 3 (X_GLXCreateContext) Value in failed request: 0x0 Serial number of failed request: 1410 Current serial number in output stream: 1411 $ `chimera --root`/bin/linuxdist /etc/os-release: NAME="Antergos Linux" VERSION="2015.11-ISO-Rolling" ID="antergos" ID_LIKE="arch" PRETTY_NAME="Antergos Linux" CPE_NAME="cpe:/o:antergosproject:antergos:2015.11" ANSI_COLOR="1;34;40" HOME_URL="https://antergos.com/" SUPPORT_URL="https://forum.antergos.com/" BUG_REPORT_URL="https://github.com/antergos" /etc/arch-release: Antergos Linux release 2015.11 (ISO-Rolling) /proc/cpuinfo: cpu: 8 GenuineIntel Intel(R) Core(TM) i7-4790 CPU @ 3.60GHz (x86_64) cpu MHz: 3884.625 cache: 8192 KB /proc/meminfo: MemTotal: 16296392 kB Graphics: 00:02.0 VGA compatible controller: Intel Corporation Xeon E3-1200 v3/4th Gen Core Processor Integrated Graphics Controller (rev 06) 01:00.0 VGA compatible controller: NVIDIA Corporation GM204 [GeForce GTX 970] (rev a1) I had to rollback to mesa-11.2.1 to make it work again. Cheers, Jaime Rodríguez-Guerra Pedregal Molecular modeling of transition metal systems group Department of Chemistry, Faculty of Sciences Edif. C · C7/153 Campus de la UAB · 08193 Bellaterra (Cerdanyola del Vallès) · Barcelona · Spain +34 93 581 28 57 www.uab.cat<http://www.uab.cat/> [http://www.uab.cat/vcard/Logo_ORCID_0.png] <http://orcid.org/0000-0001-8974-1566> [http://www.uab.cat/vcard/Logo_LinkedIn_0.png] <https://www.linkedin.com/in/jaimergp> This message may contain private or confidential information, and is addressed exclusively to its intended recipient. If you have received it in error, please notify the sender and delete it. Please bear in mind that you are not authorized to use it for any purposes. Please consider the environment before printing this email.

1 0

How to convert Coot rotation matrix to something Chimera can read?
by Oliver Clarke 19 May '16

19 May '16

Hi all, I'd like to write a little python function in coot to open the current view with Chimera, because manually replicating in Chimera views that were identified in Coot is a little painful. The file handling aspects of this seem straightforward - the one aspect I haven't quite got yet is how to set the view in Chimera to match that seen in Coot. I would like to load up all molecules/maps that are displayed in coot, in their default orientations, and then apply the view matrix from coot using matrixset to get the equivalent view. Coot outputs a 3X3 rotation matrix which seems to be the same as that generated for an equivalent view by matrixget in Chimera, except that in Chimera the matrix is 3X4, with the last column representing the translation vector. This I cannot access directly within coot, but I can probably calculate it if I know what it represents. This is probably a stupid question, but is this simply the vector from the origin to the current center of rotation? Or is it something else? Cheers, Oliver.

1 1

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

16 May '16

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 >> >>

2 5

question
by zona cat 14 May '16

14 May '16

Good afternoon, My name is Daniel Heller, and I am working with Helen Berman at Rutgers University on a PDB project. We are trying to analyze multiple DNA binding proteins simultaneously (via matchmaker) and looking to visualize the surface electrostatic potential of each structure as a tiled picture. We are deleting the nucleotide chains for all proteins when complexed with DNA. However, I am only able to get the surface electrostatic coloring for one structure (the reference structure), subsequent structures will not apply the columbic coloring after i hit apply, they just maintain the monotone default surface structure . Could you please help me with this issue. Thank you very much,Dan

2 5

Analysis of the contact maps from the MARTINI MD trajectory
by James Starlight 13 May '16

13 May '16

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

2 1

Suggestion for ChimeraX - interruptible Matchmaker?
by Oliver Clarke 12 May '16

12 May '16

Hi, Matchmaker takes a *really* long time when dealing with long chains, particularly if it has been (unwisely/accidentally) directed to find the best matching pair of chains. This is as expected, but as Chimera is unresponsive while this happens, if the wrong matchmaker parameters are accidentally selected, there is no way to recover the session except by waiting for the operation to complete or force-quitting the app. If feasible, maybe in ChimeraX it would be useful to execute this operation in a manner that allows for cancelling a structural alignment in progress? And perhaps displaying an estimate of time to completion in the taskbar? Cheers, Oli.

2 1

Averaging multiple structures
by Yehuda Halfon 11 May '16

11 May '16

Hi there, I'm using Rosetta using comparative modeling of protein, and I get 1000 new PDBs. I wanted to upload all of them together and average them into one structure and to see if the average one is the better then the last one the software creates, which needs to be the best. Thanks, Yehuda

2 4

building and positioning disordered sequence
by thereal sisterdot 11 May '16

11 May '16

Dear Chimera Users, i would like to do something in chimera, that is not exactly science but more graphics :-) i have a known structure with 5 positions that i know are crosslinking to a 100aa intrinsically disordered partner sequence. now just in order to better imagine how the disordered binding partner wraps around the structure, i would like to put 5 lysines in the open space near the known crosslinking positions in the pdb structure, and than somehow build the rest of the disordered partner sequence connecting those lysines. its not about getting a estimate of the structural arrangement of the partner. its just so i can better imagine how the partner sequence extends around the interacting residues (as the residues are not in a patch but really all over the 100 aa) or if the arrangement would be sterically possible. any ideas if and how i could achieve that in chimera ? i know that is sound more like a task for illustrator, but :-) THANKS A LOT!!! sisterdot

2 3

Join models
by shazia sabir 10 May '16

10 May '16

Hi dear all, I have to join two models in one single model. For this I have to select N and C terminal of each of the model. But I couldn't find a way how to select/pick atom from N terminal of one model and C terminal of the other model. Please can anyone help me on that? kind Regards, Shazia.

2 3

surfaces
by A G Szabo 10 May '16

10 May '16

I am trying to represent a large homodimer with 840 residues per subunit, as a surface. On a smaller protein, 350 residues I was able to get a surface representation and color different domains. When I try to produce the surface representation of the larger molecule, I receive an error message saying that I have to split it because it is too large a molecule. So how do I do that? Thanks AGS

2 1

repeat question
by A G Szabo 10 May '16

10 May '16

I am trying to make the sphere of an ion in an active site larger than what appears in the structure. I asked a similar question back in December, and on checking my notes I tried to use the command setattr ballScale 0.5 I first selected the ion using control<click). The ion was selected. I then tried the command line as I show it above. Nothing appeared to happen. Indeed I tried to do the same for an atom of a residue that I had selected using control <click> again nothing. I must have forgotten something. A message from Elaine suggested that I use setattr m ballScale 0.5 Am I supposed to select the residue and leave out the m? Thanks for your patience. Arthur Szabo

2 2

opengl error with quadro k620
by Krishnan, Raman 10 May '16

10 May '16

[cid:image001.jpg@01D1A9F2.DB70AF30] Get this error on running chimera -debug-opengl option. Using quadrok620 nvidia card on windows 7.0 thanks Director Structural Biology, Biocryst Pharmaceuticals.

1 0

rainbow color scheme
by Shreeramesh 09 May '16

09 May '16

Dear Sir, I guess blue represents N-terminal of protein in rainbow coloring and would like know the indication of other colors in rainbow scheme of chimera for protein structure. Thanks & regards *Ramesh M* Email: shreeramesh(a)gmail.com Mobile: (91) + 9646469938 (27) + (0) 844556164

2 1

I get this when i tried to start chimera 1.11rc on windows 7. any suggestions?
by Krishnan, Raman 09 May '16

09 May '16

[cid:image001.jpg@01D1A55D.8C415B10] krish Biocryst Pharmaceuticals.

2 1

Re: [Chimera-users] [chimera-dev] suggestions for better STL export - remove internal, reduce atom vertices
by Tom Goddard 09 May '16

09 May '16

Hi Jakob, We discussed the excessive STL vertices in atoms and bonds last week and I put a fix into the daily build and release candidate. http://plato.cgl.ucsf.edu/pipermail/chimera-users/2016-April/012134.html Apparently you are not using those versions since your images show the high vertex density. I’ve attached an image using the current Chimera release candidate that shows reasonable vertex density. Unfortunately eliminating internal structures and stitching together meshes for atoms, bonds and ribbons at the seams is very difficult. Yes $10000 CAD software packages can do it but we don’t have the resources to solve that problem. Since every 3d printer has to deal with this problem I expect the printer software to be able to handle it, but perhaps it is too hard a problem for the 3d printer companies as well. Perhaps there is free software that can eliminate internal structures and stitch mesh intersecting edges together, I have not researched it — I’d start looking at MeshLab. Tom > On Apr 15, 2016, at 4:33 AM, Jakob Suckale wrote: > > Dear Chimera developers, > > Chimera has the great feature of STL export, very useful for 3D printing models for research and teaching. May I make two suggestion to improve this feature even further. > > > 1 - Remove invisible internal and therefore unnecessary structures. > > Chimera generates a lot of internal structures that complicate printing and are pointless for STLs. This includes inside surfaces at the end of ribbon cartoons (red below). > <Screen Shot 2016-04-15 at 13.24.44.png> > as well as copious amounts of invisible internal surfaces for wire and ball and stick models. > <Screen Shot 2016-04-15 at 13.27.36.png> > If these were automatically removed, STLs would not only be much smaller but also much easier to print without remashing. > > > 2 - Reduce the number of vertices for atom models. > > Ribbon conversion to vertices is already good, but as can bee seen above atom models are converted using far too many vertices inflating the resulting STL files. I'm aware that this has been discussed and listed previously as a bug but the problem remains and makes for difficult to use STLs and bad 3D prints. > > Best regards, > > > Jakob Suckale, PhD > Lecturer in Biochemistry > University of Tübingen > +49-7071-29-73363 > _______________________________________________ > Chimera-dev mailing list > Chimera-dev(a)cgl.ucsf.edu > http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-dev

3 3

Alter both clipping planes at same time?
by Oliver Clarke 05 May '16

05 May '16

Hi, is there any way to alter both global clipping planes at the same time? This is handy to move the current viewable section, without altering its thickness. At the moment I use an alias for this: # move both clip planes hither or yon. # Usage: zslab x # Move hither and yon clip planes towards (positive) or away (negative) alias ^zslab clip hither $1; clip yon $1 This works fine, except when the clip planes are close together - in that case, I get an error stating that I'm trying to move the near clip plane behind the far clip plane. It would be nice to have an option to move both planes at once - e.g. maybe changing the syntax of clip so either hither, yon, or both can be specified in the same command would be possible? Cheers, Oliver.

2 1

ImportError: No module named Chimera
by Eiros Zamora, Juan 04 May '16

04 May '16

Hi, I’m running a per-frame script to display the time in a MD movie. I used to run it without problems until I downloaded the new Mac OS X 64-bit 1.11 version. (I erased the old app, also tried to re-download the old version with no success). More specifically, I load the MD movie, and in the command line I write: 2dlabels create timer text "0 ns" color white ypos .9 xpos .9 size 50 then I open a per-frame script with the following commands: from Chimera import runCommand runCommand("2dlabels change timer text '%.1f ns'" % (mdInfo['frame']*0.02)) and get the error: Problem in per-frame script: ImportError: No module named Chimera I’ve googled it and found this reply: [Chimera-users] How to install chimera as a python module http://www.cgl.ucsf.edu/pipermail/chimera-users/2014-July/010120.html I’ve tried setting the $CHIMERA_PATH environment variable to /Applications/Chimera.app/Contents/Resources/share but I still run into the same problem. Now, I don’t want to run the script directly from the command line, as it’s a per-frame script I have to run after I load the MD movie and create the 2dlabel object. Anything obvious I’m missing here? Cheers Juan

3 2

how to move the model and save as mrc
by rongjian li 04 May '16

04 May '16

Hi, I am trying to move the model/protein in some MRC file and save it as a new MRC file. I know how to do it in Windows by clicking the mouse, but I have lots of such operations expected. So I am wondering how could I realize this in Linux such that it can be done automatically. Thank you very much. Best, Rongjian

2 1

Multiple Docking Using Chimera
by SOBAHY, TURKI MOHAMMED 02 May '16

02 May '16

Hello, I have a protein of interest that I need to perform molecular docking on by multiple ligands (~1000,000). The protein is in PDB file format, while ligands are in a SDF file format. Initially, the ligands were divided into multiple SDF files (mini) for hardware sake. Then, a python script was developed to handle the docking of each mini SDF file against the protein of interest. This works nicely, however docking scores (tables) have no ligands IDs associated to distinguish each score. So, basically I have scoring tables without any key IDs!! Kindly, I want to know if you have encountered the same problem and what approach you took to solve this. Thanks,? ______________________________ Turki Sobahy Research Associate Research Center (J-04) King Faisal Specialist Hospital & Research Center (Gen.Org.)-Jeddah Branch PO Box 40047, Jeddah 21499, Saudi Arabia Tel# +966 667 7777 ext. 66736| MCD# 40716 Mobile# +966 505837027 email tsobahy(a)kfshrc.edu.sa<mailto:ESan-Diego@kfshrc.edu.sa>

2 1

Solicitud Información
by luz Martinez 02 May '16

02 May '16

Cordial Saludo Soy estudiante de licenciatura en Química de la UNIVERSIDAD PEDAGOGICA NACIONAL. Colombia, Cordialmente solicito información acerca del modelo matemático general que el programa UCSF Chimera utiliza, lo anterior con el fin llevar a cabo una investigación en el espacio académico de Química Computacional que curso. Agradezco la atención y pronta respuesta Luz Angela Martinez M. Estudiante Decimo semestre

2 1

Thank you and a suggestion
by Mohamed Elrobh 01 May '16

01 May '16

Dear UCSF Chimera, Just wish to say big thank you for your program. If you kindly allowed me to ask if there is a possibility to make a live webinar tutorial on Mohamed Elrobh, Ph.D., Assistant Professor, Room 2A 66 B5, Department of Biochemistry, College of Science, King Saud University, PO Box 2455, Riyadh 11451, Saudi Arabia Phone +966 114675944, Fax +966 114675791, Mobile +966 502072937 Email melrobh(a)ksu.edu.sa<mailto:melrobh@ksu.edu.sa> Site http://fac.ksu.edu.sa/melrobh<http://faculty.ksu.edu.sa/elrobh> ________________________________ Disclaimer: This communication is intended for the above named person and is confidential and / or legally privileged. Any opinion(s) expressed in this communication are not necessarily those of KSU (King Saud University). If it has come to you in error you must take no action based upon it, nor must you print it, copy it, forward it, or show it to anyone. Please delete and destroy the e-mail and any attachments and inform the sender immediately. Thank you. KSU is not responsible for the political, religious, racial or partisan opinion in any correspondence conducted by its domain users. Therefore, any such opinion expressed, whether explicitly or implicitly, in any said correspondence is not to be interpreted as that of KSU. KSU may monitor all incoming and outgoing e-mails in line with KSU business practice. Although KSU has taken steps to ensure that e-mails and attachments are free from any virus, we advise that, in keeping with best business practice, the recipient must ensure they are actually virus free.

1 0

Manually adjusting alignment
by Visvaldas K. 30 Apr '16

30 Apr '16

Dear All, I wasn't able to adjust alignment in the Alignment window, and search did not find an answer (because I get great number of false negatives)...I tried shift-arrow, ctrl-shif-arrow and many times but only once the region moved "by accident", and I wasn't able to repeat this feat. To simplify the task, let's say I have in the alignment window: seqA NQDF.RFRDseqB MKTFIEHKD I want to move "RF" by one residue to the left so that the top part becomes NQDFRF.RD What would be the sequence of keys, clicks/drags and so on to do that (also, how is the active region which is ready to move signified: green, green in a dashed frame, or perhaps white dashed frame)? Sincerely yours, Vis Kairys

2 3

need help joining existing pdb files with new sequence
by John Hess 28 Apr '16

28 Apr '16

HI folks, I have asked Elaine a couple questions about this project of mine. She’s been very helpful, but I’m such a beginner, I feel guilty wasting her time and the rest of the folks reading the entire chimera list as I ask how to connect this and move that. I would like to take at least 2 existing pdb files and link them to some additional peptides sequence that would be between them. The resulting pdb would describe a structure that looks like this: > protein 1 -> peptide linker -> protein 2 I’m dealing mainly with alpha helices (in coiled-coils), and I know how I want them positioned. I joined a helical peptide sequence (made using Build Structure) to one pdb file but the orientation was way off and I’m wondering how to “grab” the entire helical peptide region and move it where I want it (sorry, in the past, I have tried to modify structures with Foldit and have had a collaborator who used Rosetta to try and model my data.). Is anyone else on this list from UC Davis and willing to get together for some tutoring? Does anyone not at Davis have time to answer my questions and help? If we work well by email and I make progress, I think that this could be a $ for time or co-authorship situation. thanks, John Hess jfhess(a)ucdavis.edu

2 1

Labeling and coloring chains
by Elaine Meng 28 Apr '16

28 Apr '16

> On Apr 28, 2016, at 8:15 AM, Deniz Simsek Buck <simsekd(a)stanford.edu> wrote: > > Hi Elaine, > Thank you for your reply. > > I already changed the chain IDs to different names. I am working with ribosome structure: 44 chains are large subunit and then I want to color them the same color. > If I already changed chain A to L23, lets say, how do I choose them? > > Also, is there a way to choose different colors than the available ones and also how do I make the label the same color as the chain? > > Thank you very much, > Deniz > Hi Deniz, I do not know how you changed the chain names, so it is hard to say. In PDB format there is only one column allotted to chain ID in the ATOM lines, <http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/tutorials/pdbintro.html#co…> ...so it can only be one letter or number, as far as I know. So it would depend on what Chimera thinks are the chain IDs…. for example, it may have only read the single character in column 22 as the chain ID. You could probably tell by putting the mouse over some atom in the chain and see what information is in the pop-up balloon. But maybe you used some other format that allows more characters of chain ID. Nevertheless, it might not be possible to use the range syntax I showed in my previous message (first example included below). You could try it, but if ranges didn’t work, you’d have to specify all the IDs (again, all of what Chimera thinks are the IDs) in one long command, or in multiple separate coloring commands. > color red,r :.A-D > (meaning: color red ribbons-only in chains A-D in all models) > > color green,a #0:.B,.D > (meaning: color green atoms-only in chains B and D of model #0) Instead of changing chain IDs I would recommend using “alias” to give new names to the set of atoms and then use these new names in the coloring commands, e.g. alias L23 :.X (meaning assign the name L23 to chain X) color hot pink,r L23 … of course you could assign the name to a collection of multiple chains if you wanted, e.g. alias L23 :.X,.Y,.H You could also use alias with your existing changed chain IDs providing you can figure out what Chimera thinks they are. As for colors, you can specify an infinite number of colors as RGBs or hexadecimal or first define new color names and use those names, if the built-in set of named colors doesn’t have what you need. Color specification is described in the same “color” manpage I mentioned before. <http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/color.html> Also described in that page are how to specify coloring labels, e.g. color orange,r,l (meaning coloring ribbons and labels only) Again I’ve CC’d the recommended address for chimera questions, chimera-users(a)cgl.ucsf.edu 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 1

Three-dimensional structural model of the VP7 trimer
by Dr Omer Naseer 28 Apr '16

28 Apr '16

Hi I am s student of PhD and working on the Rota virus A, I want to created "three-dimensional structural model of the VP7 trimer". by using "UCSF Chimera", can you please tell me How to do this? because I tried a lot but unable to produce a structure like yours. Regards Dr, Omer Naseer Ph.D Scholar University of Veterinary & Animal Sciences, Lahore, Pakistan. +92 321 640 89 56

2 1

The question about fitting
by Zhongjun Hu 28 Apr '16

28 Apr '16

Hi, everyone: Once I fit one map to the reference map, there will be some fitting result. I want to use this fitting result to align the map to the reference map using another software, for example, EMAN. How do I this? It seems that I am not clear about the fitting result, could anyone give some idea? or, is there other protocol/software to do this? Thank you! Best, Zhongjun

2 2

electron density map
by 陈柳青 28 Apr '16

28 Apr '16

Hi! I display the 2Fo-Fc map in chimera, but i wonder what the meaning of step 1, step2,etc in volume viewer. The other question is what the relationship between the level in chimera volume viewer and Coot map level ? Best regards -- 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

Re: [Chimera-users] electron density map
by Oliver Clarke 28 Apr '16

28 Apr '16

Dear Liuqing, Step size refers to whether every point in the map is displayed (step 1), or every second point (step 2), or every fourth point (step 4) and so on - basically, the coarseness of the grid. Loading a map at a larger step size will have less detail, but will be faster to load and manipulate. See here (https://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/volume.html#sampling) for details. The level in Chimera is the same as in coot in absolute terms (e-/Å^3) - it is just that in Coot, by default, maps are adjusted in steps of r.m.s. (often called “sigma”), whereas the default behaviour in Chimera is to deal directly with absolute electron density levels. You can mimic the Coot behaviour, however, by adjusting the volume level on the command line with the option sdlevel, for example “volume #1 sdlevel 1” will set the volume with id #1 to 2xr.m.s. See here (https://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/volume.html#level) for details. For convenience (really due to my forgetfulness), I use an alias to do this, as follows: "alias ^map_sig volume $1 sdlevel $2" so then one can use the command “map_sig #1 5” to adjust map #1 to 5xr.m.s., which should match the corresponding level in coot. If you want to mimic how maps look in coot, adjusting some of the mesh options is handy, as follows: "alias ^cootmode volume $1 capfaces false style mesh meshlighting false squaremesh false color "#08882eefa222" step 1; sop cap off; set depthCue; set dcStart 0.2; set dcEnd 1; background solid white; cofron” Where “cofron" is an alias that mimics the centering and navigation behaviour of coot: "alias ^cofron set showcofr; cofr view; clip on” Cheers, Oli.

1 0

Re: [Chimera-users] Viewing one biological unit of a model
by Derrick Lau 28 Apr '16

28 Apr '16

Hi Elaine, Very clear instructions. Thank you very much. Best regards, Derrick Lau PhD student University of New South Wales - Molecular Machines Group - (02) 9385 8639 *Cellular and Molecular Biology* 2016-04-26 8:09 GMT+10:00 Elaine Meng <meng(a)cgl.ucsf.edu>: > Hi Derrick, > That PDB entry is an “NMR ensemble” including 100 different conformations > of the same thing. You can hide or delete any of those 100. To first see > what you have, show Model Panel from the Favorites menu. Initially all the > conformations are collapsed into a single line with ID “0.1 … 0.100” or > something like that. If you scroll down to click the “group/ungroup” > function in the right side of the Model Panel, that entry will expand out > (ungroup) to 100 entries. > > You could show any one by clicking to highlight only that line on the > left-hand side of the Model Panel, then clicking the “show only” button on > the right-hand side, for example. You can also toggle the Shown checkboxes > individually, or highlight multiple lines and “show only” those or “close” > (delete) them leaving only the others. > > Or you don’t have to use the Model Panel at all. For example, if the > structure was opened as model #0 and you wanted to hide all except the > first conformation, command: > > ~modeldisplay #0.2-100 > - OR - > ~modeldisplay #0.2-end > > …or even delete them > > close #0.2-end > > It might take a few seconds to execute the command. > > Note that a single NMR conformation is not necessarily a biological unit. > It is just one copy of the coordinate set. > 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 Apr 24, 2016, at 8:00 AM, Derrick Lau < > derrick.lau(a)student.unsw.edu.au> wrote: > > > > Dear Chimera, > > > > Hi. My name is Derrick. I have recently started learning to use Chimera > to look at the capsid structure of HIV. I would like some help in the > interface. One particular crystal model, the HIV-dimer model PDB 2M8L, > strangely loads a bunch of ribbon structures. Would you know the way of > switching off other biological unit so I can see one pair of dimer? > > > > I didn't have this problem with other capsid structures (like 3P05). > Your help is appreciated. Thank you in advance. > > > > Best regards, > > > > Derrick Lau > >

1 0

Re: [Chimera-users] [chimera-dev] Labeling and coloring chains
by Elaine Meng 27 Apr '16

27 Apr '16

Hi Deniz, (1) It is easy to specify multiple chains in the command line. Show the command line from the Favorites menu. Example commands: color red,r :.A-D (meaning: color red ribbons-only in chains A-D in all models) color green,a #0:.B,.D (meaning: color green atoms-only in chains B and D of model #0) See manpages for color command and for command-line atom-residue-chain specification: <http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/color.html> <http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/frameatom_spec.html> … and you may want to try a “getting started” tutorial: <http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/frametut.html> <https://www.rbvi.ucsf.edu/Outreach/Tutorials/GettingStarted.html> (2) For labeling the chains there are 2 general choices. A. For 3D labels that move along with the structure, there aren’t exactly chain labels, but you can choose to label one residue in each chain and then only show the chain ID and not the other residue information. For example, command “select :10” will select residue 10 in every chain that has a residue 10. Then you could use menu: Actions… Label… residue… custom, and in that dialog specify you want labels of the form: %(chain)s , Apply. Disadvantage is that they may be on top of things and hard to see. Advantage is they move around with the structure and will stay near that residue when you rotate and translate. B. you can arbitrarily add 2D labels that are just in the XY plane and do not move with the structure, containing whatever text and symbols you want. In menu, Tools… Depiction… 2D Labels. Click the Help button on that dialog to see its help. You can type in each label manually and move it around in the XY plane. You would want to get the final view for a figure or something close to it, before bothering to add 2D labels. There are examples of 2D labels in these tutorials: <http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/tutorials/bfactor.html> <http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/tutorials/density.html> C. Actually there is a 3rd related choice, just putting the mouse over an atom or residue will pop up temporary information that includes its chain ID, and sometimes that will suffice. As you can see, there are lots of tutorials and manual pages to explore. For future reference, chimera-users(a)cgl.ucsf.edu (CC’d here) is the address for Chimera usage questions. The chimera-dev list is more for programming-related questions. 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 Apr 27, 2016, at 12:48 PM, Deniz Simsek Buck <simsekd(a)stanford.edu> wrote: > Hi, > I am very new to Chimera. If you have time, I would really appreciate your answers. > I am working with a structure that contains 80 chains. > 1. I would like to color 40 chains at once the same color, rather than manually doing it one by one. > How do I do this? > 2. Also, I would like to add labels to the individual chains I choose on the structure. > How do I do this? > Thank you, > Deniz

1 0

Define plane or planes matched to hither or yon clip planes
by Oliver Clarke 27 Apr '16

27 Apr '16

Hi, I often encounter cases where I am showing a slab of density or model (or both), clipped using the global clip planes, and would like to show the viewer a thumbnail illustrating the position of these clip planes from an orientation orthogonal to the viewing direction (similar to what is shown in side view), to give context. Is there any way to define either two planes, one matched to the current yon and one matched to the global hither, or one thick plane with the thickness defined such that one side matches hither and the other yon? Alternatively I guess I can do this using the model-specific clip planes, but that's a little cumbersome to do when many models/maps are involved. Cheers, Oli

3 4

how to change map origin in the map file
by 许靖蔚 27 Apr '16

27 Apr '16

Dear all: I recently have to match two maps and calculate the CC value between them. The origin of these two maps are different and the pixel value are also not the same. I could open these two maps in chimera and align them by map fitting. However, after I matched one of them to the other and saved the map, the re-opened map was still at the same position as the previous! So do you know how to same the map after fitting? Best wishes.

2 1

HBond on every frames of a large trajectory
by Jérémie KNOOPS [531802] 26 Apr '16

26 Apr '16

Hi, What is the easiest and fastest way to get the number of H bonds on every frame of a large trajectory (>10 000 frames) and print it to a file? The hbonds command does not seem to work with the perframe command and displaying the whole trajectory to run a per-frame script is quite slow. Jérémie Knoops PhD Student Laboratory for Chemistry of Novel Materials<http://morris.umh.ac.be/default.aspx> University of Mons<http://portail.umons.ac.be/en2/> 20, Place du Parc B-7000 Mons (Belgium) Tel : +32(0)65 37 38 65<tel:%2B32%280%2965%2037%2038%2065>

3 5

Save as GIF (animated or otherwise)
by Oliver Clarke 26 Apr '16

26 Apr '16

Hi, It would be very useful if Chimera (prob ChimeraX) was able to save either single frames, or, preferably, record movies in GIF format. This is not so much out of a desire to make animated gifs (although these can be handy for presentations and use on the web), but more because it makes it very easy to create labeled movies in vector graphics programs (one can edit a single frame in a vector graphics program and then insert it and alter the duration in for example GIMP or similar), which can then be easily converted back to mp4. Otherwise adding labels and custom annotations to movies is a pain for a novice like me, though if any one has any suggestions for recommended software I'd love to hear them! I've tried using iMovie and found it incredibly frustrating for the kinds of simple things like to do - adding labels, boxes, arrows etc to highlight specific domains during a morph, for example. At the moment I make a movie in chimera and convert it to a GIF using videoblend, edit in affinity designer and GIMP and then convert back to mp4 using online tools, but a more direct route would definitely be useful. Cheers, Oliver.

6 10

need help to plot electrostatic surface charge distribution
by ashok rout 26 Apr '16

26 Apr '16

Hi, I have NMR structure of a protein. I want to plot electrostatic surface charge distribution. Can you please write me step by step commands how I can do it. I want to remove hydrogens first and then want to plot electrostatic surface charge distribution. Ashok -- thanks & regards, Ashok

2 1