<html><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">In addition to the methods Elaine describes, you could alternatively do the following.<div><br></div><div>1) match the B chains</div><div>2) open the Sequence tool for one of the A chains</div><div>3) use the Structure->Associations menu item to also associate the other chain A with the sequence</div><div>4) show the RMSD header</div><div>5) proceed as Elaine described in yesterday's mail</div><div><br></div><div>Note that in order to do this you will have to download the 1.7 daily build since a bug got fixed yesterday that prevented RMSD attributes from being assigned when only a single sequence is shown in an "alignment". Also, I hope you aren't using Linux since the daily build didn't successfully build last night for Linux, so you would have to wait until tomorrow's build for a Linux version...</div><div><br></div><div>--Eric</div><div><br><div><span class="Apple-style-span" style="border-collapse: separate; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 16px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; -webkit-text-decorations-in-effect: none; text-indent: 0px; -webkit-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font face="Helvetica" size="5" style="font: normal normal normal 16px/normal Helvetica; "> Eric Pettersen</font></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font face="Helvetica" size="5" style="font: normal normal normal 16px/normal Helvetica; "> UCSF Computer Graphics Lab</font></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font face="Helvetica" size="5" style="font: normal normal normal 16px/normal Helvetica; "> <a href="http://www.cgl.ucsf.edu">http://www.cgl.ucsf.edu</a></font></div><div><font face="Helvetica" size="5" style="font: normal normal normal 16px/normal Helvetica; "><br></font></div></div></span></div><div><div>On Apr 4, 2012, at 10:17 AM, Elaine Meng wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div>Hi Thomas,<br>Absolutely you could, but it is not as simple... There is an "rmsd" command that measures RMSD of whatever atom sets you specify, in their current positions (without fitting). So after you superimposed the structures using chain B only (Matchmaker GUI and mmaker command both have options to use specific chains for superposition), then you would need to run a series of commands something like:<br><br>rmsd #0:1.a #1:1.a<br>... calculate RMSD between all atoms of model 0 residue 1 in chain A vs. model 1 residue 1 in chain A<br><br>"rmsd" instructions:<br><<a href="http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/rmsd.html">http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/rmsd.html</a>><br><br>The two residues would have to be the same type (contain the same number of atoms). If the numbering is different in the two structures, you'd have to adjust for that, e.g.<br><br>rmsd #0:1.a #1:9.a<br><br>Then you'd have to write all the results from the Reply Log and process them into the appropriate attribute definition file format, and then read that file in to assign the RMSD values as attributes, as described for method (B) in my previous response.<br><br>Actually there is one new spiffy way to loop through a residue range with commands (Chimera 1.7, so daily builds), for example if you wanted to measure RMSD of model 0 chain A residues 1-85 vs. model 1 chain A residues 6-90:<br><br>perframe "rmsd #0:$1.a #1:$2.a" range 1,85 range 6,90<br><br>If there are gaps in the sequence alignment of the two proteins, you'd need a separate perframe command for each contiguous nongapped segment of alignment, since their numbering offset would change. Also if there are partial residues (due to missing density) or residues with alternate locations, there will be errors since the numbers of atoms in the two residues wouldn't be the same. You'd have to adjust your script to skip missing or partial residues or change to using only backbone atoms if the sidechain atoms are missing. I'd check for and delete alternate locations beforehand (for example, the Dock Prep tool has an option to do that, or you could delete all alternate locations B with "delete @.b"). "perframe" instructions:<br><<a href="http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/perframe.html">http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/perframe.html</a>><br><br>For example, I can superimpose 2gbp and 2fw0, remove alternate locations B, and then calculate the all-atom RMSDs for residues 3-300 vs 3-300. I would only have to use one substitution variable and range since these two structures have the same numbering:<br><br>open 2gbp<br>open 2fw0<br>mm #0 #1<br>delete @.b<br>perframe "rmsd #0:$1.a #1:$1.a" range 3,300<br><br>The reformatting of results into attribute assignment format would probably be the most tedious step, and possibly could be done with your own scripts. For example, the Reply Log would contain something like this<br><br>Executing rmsd ['#0:299.a', '#1:299.a']<br>RMSD between 8 atom pairs is 6.970 angstroms<br><br>... which you would want to reformat to lines like this in the assignment file:<br>[tab]#0:299.a[tab]6.970<br>[tab]#1:299.a[tab]6.970<br><br><<a href="http://www.cgl.ucsf.edu/chimera/docs/ContributedSoftware/defineattrib/defineattrib.html#attrfile">http://www.cgl.ucsf.edu/chimera/docs/ContributedSoftware/defineattrib/defineattrib.html#attrfile</a>><br><br>I hope this helps,<br>Elaine<br>-----<br>Elaine C. Meng, Ph.D. <br>UCSF Computer Graphics Lab (Chimera team) and Babbitt Lab<br>Department of Pharmaceutical Chemistry<br>University of California, San Francisco<br><br><br>On Apr 3, 2012, at 2:36 PM, Thomas Evangelidis wrote:<br><br><blockquote type="cite">Hi Elaine,<br></blockquote><blockquote type="cite">Many thanks for the detailed information! I have managed to use matchmaker and save the RMSD values as you described in method (A), but the problem is that MultialignViewer measures only the RMSD of the CA atoms that were used in the alignment. In my case I want to align chains B of the two structures but measure the RMSD between chains A. Is it possible to do that?<br></blockquote><blockquote type="cite">Thanks in advance for any help.<br></blockquote><blockquote type="cite">Thomas<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">On 3 April 2012 19:36, Elaine Meng <<a href="mailto:meng@cgl.ucsf.edu">meng@cgl.ucsf.edu</a>> wrote:<br></blockquote><blockquote type="cite">Hi Thomas,<br></blockquote><blockquote type="cite">Two general possibilities:<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">(A) If the RMSD is to be calculated by Chimera, you would need to superimpose the proteins to calculate the RMSD, then color by that value using Render by Attribute or command rangecolor. However, after the coloring step you could move the proteins apart. The RMSD calculation in Chimera uses CA atoms, so if you only have two proteins, the per-residue "RMSD" boils down to simply the CA-CA distance between the two proteins in the superimposed state.<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><<a href="http://www.cgl.ucsf.edu/chimera/docs/ContributedSoftware/multalignviewer/multalignviewer.html#headers">http://www.cgl.ucsf.edu/chimera/docs/ContributedSoftware/multalignviewer/multalignviewer.html#headers</a>><br></blockquote><blockquote type="cite"><<a href="http://www.cgl.ucsf.edu/chimera/docs/ContributedSoftware/render/render.html#render">http://www.cgl.ucsf.edu/chimera/docs/ContributedSoftware/render/render.html#render</a>><br></blockquote><blockquote type="cite"><<a href="http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/rangecolor.html">http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/rangecolor.html</a>><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">(B) Alternatively, if the RMSD is calculated by some other program you could create an "attribute definition" file that lists the values, read the file into Chimera to assign the attribute values to residues (or atoms), and then color by your custom attribute the same way as in (A). The attribute file format is quite simple, see description and examples:<br></blockquote><blockquote type="cite"><<a href="http://www.cgl.ucsf.edu/chimera/docs/ContributedSoftware/defineattrib/defineattrib.html#attrfile">http://www.cgl.ucsf.edu/chimera/docs/ContributedSoftware/defineattrib/defineattrib.html#attrfile</a>><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">However, method (A) is probably easiest if the CA-CA distance is an acceptable measure. Example using MatchMaker command (mm) to superimpose and show sequence alignment:<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">open 2gbp<br></blockquote><blockquote type="cite">open 2fw0<br></blockquote><blockquote type="cite">mm #0 #1 show true<br></blockquote><blockquote type="cite"> ... sequence alignment should appear, including RMSD histogram above the sequences ...<br></blockquote><blockquote type="cite"><<a href="http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/matchmaker.html">http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/matchmaker.html</a>><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">To define and apply the coloring, you could use Render by Attribute (see the attribute of residues named mavRMSD) or rangecolor, for example:<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">rangecolor mavRMSD min purple mid red max yellow<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Then you can un-superimpose the structures, for example, with command: reset, or manually with the mouse.<br></blockquote><blockquote type="cite"><<a href="http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/mouse.html#activedef">http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/mouse.html#activedef</a>><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">I hope this helps,<br></blockquote><blockquote type="cite">Elaine<br></blockquote><blockquote type="cite">----------<br></blockquote><blockquote type="cite">Elaine C. Meng, Ph.D.<br></blockquote><blockquote type="cite">UCSF Computer Graphics Lab (Chimera team) and Babbitt Lab<br></blockquote><blockquote type="cite">Department of Pharmaceutical Chemistry<br></blockquote><blockquote type="cite">University of California, San Francisco<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">On Apr 3, 2012, at 5:46 AM, Thomas Evangelidis wrote:<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><blockquote type="cite">Dear UCSF Chimera users and developers,<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">Is it possible to color a pair of structures representing different stages of the thermodynamic cycle of a protein by the RMSD, without aligning them?<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">thanks in advance,<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">Thomas<br></blockquote></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">-- <br></blockquote><blockquote type="cite">======================================================================<br></blockquote><blockquote type="cite">Thomas Evangelidis<br></blockquote><blockquote type="cite">PhD student<br></blockquote><blockquote type="cite">Biomedical Research Foundation, Academy of Athens<br></blockquote><blockquote type="cite">4 Soranou Ephessiou , 115 27 Athens, Greece<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">email: <a href="mailto:tevang@bioacademy.gr">tevang@bioacademy.gr</a><br></blockquote><blockquote type="cite"> <span class="Apple-tab-span" style="white-space:pre"> </span><a href="mailto:tevang3@gmail.com">tevang3@gmail.com</a><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">website: <a href="https://sites.google.com/site/thomasevangelidishomepage/">https://sites.google.com/site/thomasevangelidishomepage/</a><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">_______________________________________________<br></blockquote><blockquote type="cite">Chimera-users mailing list<br></blockquote><blockquote type="cite"><a href="mailto:Chimera-users@cgl.ucsf.edu">Chimera-users@cgl.ucsf.edu</a><br></blockquote><blockquote type="cite"><a href="http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users">http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users</a><br></blockquote><br><br>_______________________________________________<br>Chimera-users mailing list<br><a href="mailto:Chimera-users@cgl.ucsf.edu">Chimera-users@cgl.ucsf.edu</a><br>http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users<br></div></blockquote></div><br></div><div apple-content-edited="true"><span class="Apple-style-span" style="border-collapse: separate; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; -webkit-text-decorations-in-effect: none; text-indent: 0px; -webkit-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><br class="Apple-interchange-newline"></div></span> </div><br></body></html>