- Chimera-users - RBVI Mailing Lists

Measuring the distance between residues and a surface
by Stanley, Julian A 25 Jul '19

25 Jul '19

Hello! I'm fairly new to Chimera, so if you're able to point me towards any resources that can help me answer questions like this one independently in the future, please do! Here's the question: I would like to measure the distance between the ends arginine and lysine residues and a surface, using Chimera's python modules. Typically, I would load my .pdb file, generate a surface, and then use the "measure distance" command and parse the output from the Reply Log. This approach works fine, except since "measure distance" prints to the reply log, I need to run my script in GUI mode. If possible, I'd prefer to run everything from the command line only. I thought that I might be able to do this with the distance() command, but I can't figure out how to pass a surface to that command. Does anyone have any tips/advice? In summary, what I usually do: ``` from chimera import runCommand as rc # load the file rc("open " + pdbFile) # Generate a surface rc("select protein") rc("split") rc("surface allComponents false") rc("~select") # Measure distance between arg and lys and the surface rc("select #0") rc("~select #0:@") rc("measure distance :arg@cd|:lys@ce selection multiple true show true") ``` And then I parse out the distances from the Reply Log. What I would prefer to do: It would be better to have the output to go to STDOUT when I run the script. I imagine that I'd use something like the `distance()` command, of the premade scripts (like here: http://plato.cgl.ucsf.edu/trac/chimera/attachment/wiki/Scripts/atomdist.py) but they are designed to measure the distance between atoms or residues, not between residues and a surface. Thanks for your help! Sincerely, Julian

3 2

Modeller Restraints
by Danilo Boškovic 24 Jul '19

24 Jul '19

Dear all, how do i add restraints specific to my templates that i use for comparative modelling in chimera? Specifically what are the names assigned to those templates? Do i need the names to specify the template residues that i want to restrain? Are the names same as in the tab when i choose the templates for the run? Thanks, D. B.

3 5

Re: [Chimera-users] Cannot resample to voxel size of 1
by Tom Goddard 22 Jul '19

22 Jul '19

It sounds like you just want to resample a map at a new grid spacing. Chimera does not have such a feature. The vop resample command is intended make the grids of two different maps match each other. Tom > On Jul 22, 2019, at 1:54 PM, Ryan Harlich <ryanharlich(a)hotmail.com> wrote: > > Hi Tom, > > I need the voxel size to be 1, so subdivideSurface and subdivisionLevel options won’t work. Also, vop new won’t work because I need it to have the same bounding box as the input surface. > > Best, > Ryan Harlich > > Sent from Mail <https://go.microsoft.com/fwlink/?LinkId=550986> for Windows 10 > > From: Tom Goddard <goddard(a)sonic.net> > Sent: Monday, July 22, 2019 1:49:38 PM > To: Ryan Harlich <ryanharlich(a)hotmail.com> > Cc: chimera-users(a)cgl.ucsf.edu BB <chimera-users(a)cgl.ucsf.edu> > Subject: Re: [Chimera-users] Cannot resample to voxel size of 1 > > Hi Ryan, > > Resampling only works on the grid of another map as you can see from the documentation so "vop new" allows you to make that map. If you just want to make the mesh of a surface depiction of a map finer then look at the volume command subdivideSurface and subdivisionLevel options. > > https://www.rbvi.ucsf.edu/chimerax/docs/user/commands/volume.html#dispsurf <https://eur03.safelinks.protection.outlook.com/?url=https%3A%2F%2Fwww.rbvi.…> > > Tom > > >> On Jul 22, 2019, at 12:03 PM, Ryan Harlich <ryanharlich(a)hotmail.com <mailto:ryanharlich@hotmail.com>> wrote: >> >> Hi Elaine, >> >> I think you missed the point. I am trying to take an existing surface and resample it to a grid Spacing of 1. I am not sure I see the point of making a new map. >> >> Best, >> Ryan Harlich >> >> Sent from Mail <https://eur03.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgo.micros…> for Windows 10 >> >> From: Elaine Meng <meng(a)cgl.ucsf.edu <mailto:meng@cgl.ucsf.edu>> >> Sent: Monday, July 22, 2019 11:59:06 AM >> To: Ryan Harlich <ryanharlich(a)hotmail.com <mailto:ryanharlich@hotmail.com>> >> Cc: chimera-users(a)cgl.ucsf.edu <mailto:chimera-users@cgl.ucsf.edu> BB <chimera-users(a)cgl.ucsf.edu <mailto:chimera-users@cgl.ucsf.edu>> >> Subject: Re: [Chimera-users] Cannot resample to voxel size of 1 >> >> Hi Ryan, >> You can make a new map of whatever grid spacing you want without any atoms using “vop new”: >> <https://nam05.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.rbvi.u… <https://eur03.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.rbvi.u…>> >> >> I hope this helps, >> Elaine >> ----- >> Elaine C. Meng, Ph.D. >> UCSF Chimera(X) team >> Department of Pharmaceutical Chemistry >> University of California, San Francisco >> >> > On Jul 22, 2019, at 11:33 AM, Ryan Harlich <ryanharlich(a)hotmail.com <mailto:ryanharlich@hotmail.com>> wrote: >> > >> > Hi Chimera Bug Team, >> > >> > I am unable to resample surface to a grid spacing of 1 to be 1 voxel. >> > >> > Currently I do this: >> > 'molmap #2 6 gridSpacing 1\n' >> > 'vop resample #3 onGrid #2.1\n' >> > >> > But this requires me to have a pdb structure. I would like to do this without a pdb structure. >> > >> > Another alternative is I could automatically without GUI interface make a 3 atom pdb structure with python script that is the size of the surface’s boundary box. Although there is no way to get surface’s boundary box x,y,z. >> > >> > Any help would be much appreciated. Option 1 would be a much better approach. >> > >> > Best, >> > Ryan Harlich >> >> _______________________________________________ >> Chimera-users mailing list: Chimera-users(a)cgl.ucsf.edu <mailto:Chimera-users@cgl.ucsf.edu> >> Manage subscription: http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users <https://eur03.safelinks.protection.outlook.com/?url=http%3A%2F%2Fplato.cgl.…>

1 0

Cannot resample to voxel size of 1
by Ryan Harlich 22 Jul '19

22 Jul '19

Hi Chimera Bug Team, I am unable to resample surface to a grid spacing of 1 to be 1 voxel. Currently I do this: 'molmap #2 6 gridSpacing 1\n' 'vop resample #3 onGrid #2.1\n' But this requires me to have a pdb structure. I would like to do this without a pdb structure. Another alternative is I could automatically without GUI interface make a 3 atom pdb structure with python script that is the size of the surface’s boundary box. Although there is no way to get surface’s boundary box x,y,z. Any help would be much appreciated. Option 1 would be a much better approach. Best, Ryan Harlich Sent from Mail <https://go.microsoft.com/fwlink/?LinkId=550986> for Windows 10

3 3

same volume data values over the entire volume
by Dieter Blaas 22 Jul '19

22 Jul '19

Hi, Trying to color a volume according to "volume data" I get the same numerical values for all colors and a speckled appearance. When using "volume data gradient norm" I get the expected colors, i.e. strong local changes of the volume data are red and spots where there is essentially no change are in blue. How is this possible? If I understand correctly, "volume data gradient norm" colors according to the first derivative of "volume data" or is this wrong? Do I need something like a normalization? Thanks a lot for explanations, best, Dieter ------------------------------------------------------------------------ Dieter Blaas, Max F. Perutz Laboratories Medical University of Vienna, Inst. Med. Biochem., Vienna Biocenter (VBC), Dr. Bohr Gasse 9/3, A-1030 Vienna, Austria, Tel: 0043 1 4277 61630, Fax: 0043 1 4277 9616, e-mail: dieter.blaas(a)meduniwien.ac.at ------------------------------------------------------------------------

4 4

Fitting the structure into the map after changing origin
by Ahmad Khalifa 22 Jul '19

22 Jul '19

Is there a shortcut way of transforming the model to fit the map after I change, center or reset the origin? Best regards.

3 2

RMSD Chimera
by Ishita Gupta 22 Jul '19

22 Jul '19

Hi, I have a few questions in mind regarding rmsd calculation in chimera: 1. Does the matchmaker tool calculate c alpha? 2. Do we consider RMSD of pruned atoms or the total atoms while reporting? 3. If we use structural comparison match align tool then rmsd is based on multiple sequence alignment? And it also gives C alpha? 4. If we calculate rmsd from multi align viewer by clicking on the structure and then match. What these rmsd values denote? Please elaborate on these points. Thanks Ishita

2 1

Measurement similarity protein structure
by Nail Besli 18 Jul '19

18 Jul '19

Hi I superimposed two proteins in chimera with winch algorithms( blossom62) and p250. The result is the same. Why is the result same for both two proteins, although I tried two different algorithm settings? I am using the chimera's MatchMaker tool. The only thing different is the sequence alignment score. I attached the result file. kind regards

3 7

what laptop configuration should I get for using Chimera and ChimeraX with very large volumes?
by Ilona Nudelman 18 Jul '19

18 Jul '19

Hi, I'm looking for an upgrade to my PC laptop (Windows) and wanted to ask if anyone has any advice on the configuration needed to run Chimera and ChimeraX with very large volumes and lots of heavy pdb's docked in. I am looking for a mobile workstation type laptop but I don't want to spend too much money - maybe 2,500-3,500$. Within these confines what memory, GPU, processor etc. would you recommend? Thanks! Ilona

3 2

how to obtain information regarding complexes
by Quadir, Farhan (MU-Student) 18 Jul '19

18 Jul '19

Hello, I am relatively new to chimera. I am stuck with a problem of how to perform the following from a pdb file which contains multiple chains: 1. I would like to find regions which are 5.0 Angstroms or less between selected chains. 2. After obtaining that I would like to get from (1) a list of the residues and their positions in the chain that are interacting with each other inter-chainwise. Regards, Farhan Quadir PhD Student and Graduate Research Assistant Bioinformatics, Data Mining and Machine Learning (DBM) Lab 304 Engineering Building 1 Electrical Engineering and Computer Science (EECS) Department University of Missouri, Columbia, MO 65201, USA Email: fqg7h(a)mail.missouri.edu Phone (Cell): 573-846-4047

2 1