- Chimera-users - RBVI Mailing Lists

failing to open saved PDBs
by Healey, Joe 28 Nov '16

28 Nov '16

Hi (again) Chimera Team! Encountered a problem where chimera seems to fail to open PDBs that it created. I saved multiple PDBs to a single file as I wanted to open them and have them behave like a single structure. I started a new session with another PDB open that I want to compare the multiple PDBs to, went to open the multiple PDBs from before that I saved as a single file and nothing appears. Chimera claims to have opened the files successfully and I get no errors etc. Am I maybe misunderstanding the point/use of the multiPDB file form? 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>

2 1

Fit in map equivalent for surfaces?
by Healey, Joe 28 Nov '16

28 Nov '16

Hi Chimera Team, Is it possible/does a function exist to fit one structure inside the calculated surface of another? To be clear, i don't mean inside a density map, I mean an MSMS surface calculated from a PDB inside chimera itself. I have 2 proteins that are sequentially very different, but I suspect fulfilling the same role. They don't overlay via MatchMaker well, but they look to have quite similar space occupation, so I'd like to see how well they fit in one anothers 'map'. Any alternative suggestions appreciated! Thanks, Joe 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>

2 3

angle between subunits
by chemocev marker 26 Nov '16

26 Nov '16

Hi All I am interested to measure the angle between chain A and chain B of heterodimer, and both sequence identity is around 14%. I can measure the angle by selecting 3 atoms on chain A and chain B but this is bias because, if I change the residues then the angle will be different. I can also measure by structure measurement options and defining the axes and all the create axes for all the helices in the chain A and chain B and measure the angle between two helices on chain A and chain B. I am not sure if its a right way to claim the angle between chain A and chain B in one PDB file. I also try to use the "match :.A at CA <http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users> :.B at CA <http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users> move false show true" But it gives the rotation angle. As I am comparing the different homologous structureswhere chain B shows more identity but the each structure have a different state giving a unique angle between chains. What can be best option to measure this??? best Jiri

2 1

simple question
by A G Szabo 25 Nov '16

25 Nov '16

I should be able to find this but haven't looked too hard. I am interested in a salt bridge between the C-terminal carboxylate and a Lys residue in the protein. I know how to visualize the Lys side chain, and the C-terminal residue. How do I visualize the carboxylate terminus? I tried backbone only actions, but nothing worked. Thanks AGSzabo

2 5

Sequence tool
by Hernando J Sosa 23 Nov '16

23 Nov '16

Dear all Is there a way using the sequence (or any other) tool in chimera to save a text file with sequence and corresponding secondary structure assignments as displayed in the sequence tool. Thanks Hernando

2 2

usage of mcopy
by Jiying Jia 23 Nov '16

23 Nov '16

Hi, I want to use python script to adjust the bond angle randomly and only accept the case if the random new angle is larger than old one. I use "mcopy" to recover the old angle if the new is smaller. But I encounter a problem: if I use "mcopy #old #new setting a" or "x" or "y", the after "adjust angle anglenew" command, the angle adjusted is not the anglenew, but a different one. I don't know why. I hope anyone can solve my problem.. Below is the python script and init.mol2. PYTHON SCRIPT: import os import time import chimera import random from chimera import runCommand as rc from chimera.tkgui import saveReplyLog import MoleculeTransform from MoleculeTransform import transform_atom_coordinates import numpy as np import sys import re #open initial conformation init.mol2 opened = chimera.openModels.open("init.mol2")# opened as model 0 in chimera chimera.openModels.open("init.mol2") # opened as model 1 in chimera, used to store old conformation in monte carlo simulation rho = opened[0] print "angleold anglenew angle_after_adjust" for pdbi in range(0,20): random_atom = 1 an0 = rho.atoms[random_atom-1] an1 = rho.atoms[random_atom] an2 = rho.atoms[random_atom+1] angleold = chimera.angle(an0.xformCoord(),an1.xformCoord(),an2.xformCoord()) anglenew = random.random()*180 rc("adjust angle %f %s %s %s"%(anglenew,an0.oslIdent(),an1.oslIdent(),an2.oslIdent())) angleafter = chimera.angle(an0.xformCoord(),an1.xformCoord(),an2.xformCoord()) print str(angleold)+" "+str(anglenew)+" "+str(angleafter) if anglenew>angleold: rc("mcopy #0 #1 settings x") else: rc("mcopy #1 #0 settings x") rc("close all") rc("stop now") INIT.MOL2: @<TRIPOS>MOLECULE init.mol2 3 2 3 0 0 NUCLEIC_ACID NO_CHARGES @<TRIPOS>ATOM 1 C 1.7324 -0.5083 0.8663 C.3 1 DA 0.0000 2 C -0.6180 0.4940 -0.0060 C.3 2 DA 0.0000 3 C 1.2450 0.9460 1.8950 C.3 3 DA 0.0000 @<TRIPOS>BOND 1 1 2 1 2 2 3 1 @<TRIPOS>SUBSTRUCTURE 1 DA 1 RESIDUE 4 A DA 1 ROOT 2 DA 2 RESIDUE 4 A DA 2 3 DA 3 RESIDUE 4 A DA 1

2 3

How to place a small volumed object in place of markers in marker set?
by BRANDEN WAYMON WEST 23 Nov '16

23 Nov '16

How do I place an object with small volume in the locations of markers in a marker set? Essentially, I would like to replace the sphere markers with copies of my own object. This is my first post in this mailing list. I very much appreciate all of your help. Thank you. Branden West

2 1

how to find clashes/contacts for residue of our interest?
by Seera Suryanarayana 23 Nov '16

23 Nov '16

Dear chimera users, I am new to chimera and I looked into the manual. But I did not get perfectly that how to use find clashes command. I would like to know clashes/contacts of residue of my interest. But I am not able to select the residue in graphical interface under select tool. Kindly tell me how to find clashes/contacts ASP32 in a protein. Thanks in advance Surya Graduate student India.

2 1

Rotating marker sets
by Kacher, Joshua P 22 Nov '16

22 Nov '16

Hello, Is there a way to rotate marker sets? I have a 3D crystal system that is defined by a hexagonal coordinate system. I can load in a marker set in the form of a hexagon to act as the coordinate axes, but I need to then rotate the marker set to align it with my system. So far, I've only been able to see how to move one marker at a time, not the entire marker set. Thanks, Josh -- Josh Kacher Assistant professor Georgia Institute of Technology KacherLab.gatech.edu

2 1

Render by Attribute - Radii and individual radius for each HETATM
by Kuda Ondřej 22 Nov '16

22 Nov '16

Hello, I have a path (chain of heteroatoms with defined radii, see below) forming a pore in a protein. I would like to create an envelope (surface view) to illustrate the 3D path within the protein. My problem is that Chimera uses default vdw value and not the values I defined. Therefore, the "pore" has a constant diameter. According to https://www.cgl.ucsf.edu/chimera/docs/ContributedSoftware/render/render.htm… , I tried "Render by Attribute - Occupancy - Radii". I assumed that the VDW values are rendered for each HETATM according to the input, but it seems that the VDW values are sorted (smallest > largest) and then applied on the chain starting from #1. Therefore, it creates kind of a funnel structure. Could you please suggest how to create this tunnel with variable width? Thank you. Best, Ondrej REMARK ATOM NAM RES TUNID X Y Z Distnm RadiusA HETATM 1 X TUN T 1 -45.770 -35.888 -4.154 0.00 2.01 HETATM 2 X TUN T 1 -45.766 -35.893 -4.132 0.02 2.03 HETATM 3 X TUN T 1 -45.762 -35.897 -4.111 0.04 2.11 HETATM 4 X TUN T 1 -45.759 -35.901 -4.091 0.07 2.16 HETATM 5 X TUN T 1 -45.756 -35.903 -4.073 0.08 2.14 HETATM 6 X TUN T 1 -45.755 -35.903 -4.057 0.10 2.17 HETATM 7 X TUN T 1 -45.756 -35.901 -4.045 0.11 2.10 HETATM 8 X TUN T 1 -45.758 -35.896 -4.037 0.12 2.01 HETATM 9 X TUN T 1 -45.763 -35.888 -4.032 0.13 2.01 HETATM 10 X TUN T 1 -45.770 -35.876 -4.033 0.15 2.00 HETATM 11 X TUN T 1 -45.779 -35.861 -4.037 0.16 1.99 HETATM 12 X TUN T 1 -45.787 -35.843 -4.039 0.18 1.97 HETATM 13 X TUN T 1 -45.794 -35.822 -4.034 0.21 1.95 HETATM 14 X TUN T 1 -45.796 -35.800 -4.017 0.23 1.93 HETATM 15 X TUN T 1 -45.792 -35.777 -3.984 0.27 1.92 ________________________________ This email was scanned by Bitdefender Upozornění: Není-li v této zprávě výslovně uvedeno jinak, má tato E-mailová zpráva nebo její přílohy pouze informativní charakter. Tato zpráva ani její přílohy v žádném ohledu Fyziologický ústav AV ČR, v. v. i. k ničemu nezavazují. Text této zprávy nebo jejích příloh není návrhem na uzavření smlouvy, ani přijetím případného návrhu na uzavření smlouvy, ani jiným právním jednáním směřujícím k uzavření jakékoliv smlouvy a nezakládá předsmluvní odpovědnost Fyziologického ústavu AV ČR, v. v. i. Disclaimer: If not expressly stated otherwise, this e-mail message (including any attached files) is intended purely for informational purposes and does not represent a binding agreement on the part of Institute of Physiology CAS. The text of this message and its attachments cannot be considered as a proposal to conclude a contract, nor the acceptance of a proposal to conclude a contract, nor any other legal act leading to concluding any contract; nor does it create any pre-contractual liability on the part of Institute of Physiology CAS.

2 3