[Chimera-users] mmaker question
meng at cgl.ucsf.edu
Thu Aug 21 11:00:44 PDT 2014
Not sure I understand the question. I'll try to answer some possible interpretations, but if I didn't answer your question, feel free to write back!
(A) To use your own sequence alignment to define the residue pairing for superposition and RMSD calculation, just open your alignment, make sure both structures are associated, and then from the Multalign Viewer menu (the menu of the window showing the sequence alignment) choose Structure…Match. There are some additional options in the resulting GUI. You could use the pairwise alignment from a previous mmaker command with "show true", if that helps.
(B) To just calculate the RMSD of structures in their current positions WITHOUT moving them, again starting with an alignment in Multalign Viewer with both structures associated, use the mouse to draw a region (which could be a single box in the alignment, or multiple disjoint boxes if you use Shift-drag) then from the Multalign Viewer menu choose Info…Region Browser. The resulting dialog will show the RMSD for the region you defined with the mouse.
(C ) The above use the Multalign Viewer graphical interface, not commands, and use only CA atoms. There is also a "rmsd" command to calculate RMSD between 2 sets of atoms in their current positions without moving (it does NOT calculate the best-fit RMSD, just uses the current positions). It would be harder to use because now you have to specify all the atoms in the command, e.g. if the residue numbering is different in the two structures, you would have to list the residue ranges. You could specify any atoms common to both structures, not just CA atoms.
Elaine C. Meng, Ph.D.
UCSF Computer Graphics Lab (Chimera team) and Babbitt Lab
Department of Pharmaceutical Chemistry
University of California, San Francisco
On Aug 21, 2014, at 8:49 AM, Joshua Broyde <jb3401 at c2b2.columbia.edu> wrote:
> Hi again,
> I was wondering if the mmaker command had an option (or if there is some other way within chimera) of reporting the RMSD (as calculated by mmaker) of an alignment without actually superimposing one structure on the other. The reason that I am doing this that for a group of structures I can see with mmaker that they superimpose very well, but I want to show that when you compute the RMSD of those same structures when each is bound to a partner they are not quite as well aligned.
> In essence, what I want mmaker to do is not really realign the structures for me, but to calculated the RMSD of an alignment that I provide it.
> I looked at the mmaker command page but I couldn't seem to find an option for this.
> Thanks again, and Chimera is a wonderful piece of software!
> Joshua Broyde
> Chimera-users mailing list
> Chimera-users at cgl.ucsf.edu
More information about the Chimera-users