[Chimera-users] Percent Identity for specified parts of a structure
meng at cgl.ucsf.edu
Tue Nov 8 11:07:49 PST 2011
There isn't a command, but there is an option in the menu of Multalign Viewer (the Chimera tool that shows sequences and sequence alignments). Currently it can only use the entire sequence alignment, but with Chimera you could save and re-open sub-alignments containing just the helix of interest. It is a little tedious, but at least you wouldn't have to duplicate your structure files.
You can simply open a sequence alignment file in Chimera,
or If you don't yet have a sequence alignment, it can be created in Chimera:
If you want Chimera to do the fitting/superposition and if the sequences are reasonably easy to align (similar at least in the parts you care about), you could use MatchMaker with "Show pairwise alignments" turned on.
If the sequences are not easy to align, and/or you already have the proteins superimposed the way you want and don't want to change that, use Match->Align to create a sequence alignment from the pre-existing superposition.
These tools are under Tools... Structure Comparison in the menu.
Then, when you have the overall sequence alignment, you can highlight the part of interest by dragging a box with the mouse, and save just that part with "File... Save as" in the sequence window menu, turning on the option to "Restrict save to active region." Then open that newly saved shorter alignment in Chimera, and use its menu item "Tools... Percent Identity."
You can show the secondary structure assignments on the overall alignment (may help in dragging a box in the right place) with sequence window menu item "Structure... Secondary Structure... show actual"
I hope this helps,
Elaine C. Meng, Ph.D.
UCSF Computer Graphics Lab (Chimera team) and Babbitt Lab
Department of Pharmaceutical Chemistry
University of California, San Francisco
On Nov 8, 2011, at 7:20 AM, Victoria James wrote:
> Hi there,
> I don't know if this command already exists, but I can't seem to find it in the Chimera tools or manual.
> I'd like to be able to calculate the percent identity of specified parts of 2 matched sequences - in my case the individual transmembrane helices of the structures.
> I wonder if a command already exists or if there is a way of doing this in Chimera without having to make new .pdb files of each of the areas of the proteins that I want to compare?
> I think it would be a useful command if it does not already exist.
More information about the Chimera-users