<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=Windows-1252">
</head>
<body>
<div dir="ltr">
<div dir="ltr" data-ogsc="" style="">
<div dir="ltr">
<div dir="ltr" data-ogsc="" style="">
<div dir="ltr">
<div dir="ltr" data-ogsc="" style="">
<div>
<meta content="text/html; charset=Windows-1252">
</div>
<div>
<div>
<div style="direction:ltr">Hi Elaine,</div>
<div style="direction:ltr" dir="ltr"></div>
</div>
<div>Yes I tried but it always gave me error message saying unequal numbers of atoms chosen when I used n,ca,c,o. However they are the same protein so shouldn’t have this problem. Now I tried “all” and it works so I think I will just use this value. Thank you
so much!</div>
<div dir="ltr">Best,</div>
<div dir="ltr">Bixia</div>
<div class="ms-outlook-ios-signature"></div>
</div>
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size: 11pt;" data-ogsc=""><b>From:</b> Elaine Meng <meng@cgl.ucsf.edu><br>
<b>Sent:</b> Friday, September 18, 2020 8:44:22 AM<br>
<b>To:</b> Zhang, Bixia <bixia.zhang@wsu.edu><br>
<b>Cc:</b> Chimera <chimera-users@cgl.ucsf.edu><br>
<b>Subject:</b> Re: [Chimera-users] Backbone RMSD</font>
<div> </div>
</div>
<div style="word-wrap:break-word; line-break:after-white-space">Hi Bixia,
<div class="">Did you read my previous reply? In that previous message, also shown below, I said to use the "match" command instead of matchmaker to specify whatever atoms you want (could be all atoms or all backbone, for example) for fitting/RMSD -- or if
you want RMSD without fitting (without moving the structures from their current positions) use the "rmsd" command. There are links to the command help pages and tutorial examples in that message. Those commands would give a single overall RMSD value for
the atoms that you specified.</div>
<div class=""><br class="">
</div>
<div class="">If you don't mean a single RMSD value but instead that you want to see it at each position of a sequence alignment, then use matchmaker option to show the pairwise sequence alignment. The sequence alignment window shows above the sequences a
histogram "RMSD: ca" but you can use the Headers menu of that window to show "RMSD: full" or "RMSD: backbone" or all three if you want, as shown in the example image. They might not be very different from each other, however, as in my example.</div>
<div class=""><br class="">
</div>
<div class="">This header stuff is explained in the help for the sequence viewer (Multalign Viewer):</div>
<div class=""><<a href="https://urldefense.com/v3/__http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/multalignviewer/framemav.html__;!!JmPEgBY0HMszNaDT!99VUICUgB4sdYAGMQoM-Z4t1ndgRXPIRZ3V2DcLU2xBQszvFD0kLpSeAbUFgHzwUjeA$" class="" data-ogsc="" style="">http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/multalignviewer/framemav.html</a>></div>
<div class=""><<a href="https://urldefense.com/v3/__http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/multalignviewer/multalignviewer.html*headers__;Iw!!JmPEgBY0HMszNaDT!99VUICUgB4sdYAGMQoM-Z4t1ndgRXPIRZ3V2DcLU2xBQszvFD0kLpSeAbUFgVqHKdJc$" class="" data-ogsc="" style="">http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/multalignviewer/multalignviewer.html#headers</a>></div>
<div class=""><br class="">
</div>
<div class="">Best,</div>
<div class="">Elaine<br class="">
<div class="">-----<br class="">
Elaine C. Meng, Ph.D. <br class="">
UCSF Chimera(X) team<br class="">
Department of Pharmaceutical Chemistry<br class="">
University of California, San Francisco<br class="">
<br class="">
</div>
<div class=""><img id="x_B05822C2-CDA1-4643-8AFC-1ABFAEC548FB" width="640" height="378" src="cid:CEBD7333-C1DC-46E8-A0CC-FE9A7FDCA451@gateway.sonic.net" class=""></div>
<blockquote type="cite" class="">On Sep 17, 2020, at 9:05 PM, Zhang, Bixia <<a href="mailto:bixia.zhang@wsu.edu" class="" data-ogsc="" style="">bixia.zhang@wsu.edu</a>> wrote:<br class="">
<br class="">
Hi Elaine,<br class="">
I was trying to figure out the backbone RMSD yesterday and I just realized that, as you said, by matchmaker I already got RMSD for ca. So now I am wondering if there’s any way to get RMSD for the full protein (including side chain). I am investigating the solution
affect on structure, so the change is very tiny, if I can see the RMSD for full protein and backbone I might get some ideas on where it changes. Thank you!<br class="">
Best,<br class="">
Bixia <br class="">
<br class="">
From: Elaine Meng <<a href="mailto:meng@cgl.ucsf.edu" class="" data-ogsc="" style="">meng@cgl.ucsf.edu</a>><br class="">
Sent: Wednesday, September 16, 2020 12:24<br class="">
To: Zhang, Bixia<br class="">
Cc: chimera-users@cgl.ucsf.edu<br class="">
Subject: Re: [Chimera-users] Backbone RMSD<br class="">
<br class="">
Hi Bixia,<br class="">
Matchmaker always uses only one atom per residue, CA for amino acid residues.<br class="">
<https://urldefense.com/v3/__http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/superposition.html__;!!JmPEgBY0HMszNaDT!5O_WiKmxcIblQ26zQKLIYV7Eg7QfEMgarJp3JpDss8ggtnmYbFSbMFNmzLk6Mw2g-s0$ %3E<br class="">
<https://urldefense.com/v3/__http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/matchmaker/matchmaker.html__;!!JmPEgBY0HMszNaDT!5O_WiKmxcIblQ26zQKLIYV7Eg7QfEMgarJp3JpDss8ggtnmYbFSbMFNmzLk6o96l658$ %3E<br class="">
<br class="">
However, you can do fitting with the "match" command on specific atoms, or "rmsd" command just to measure the current RMSD of those specific atoms without fitting. The harder part is in that case you have to specify all the atoms in the command and there
must be equal numbers of atoms from the two structures, which will be paired in order. This can be difficult if the structures have different numbers of residues and/or different residue numbering. <br class="">
<https://urldefense.com/v3/__http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/match.html__;!!JmPEgBY0HMszNaDT!5O_WiKmxcIblQ26zQKLIYV7Eg7QfEMgarJp3JpDss8ggtnmYbFSbMFNmzLk6YjD9ZH8$ %3E<br class="">
<https://urldefense.com/v3/__http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/rmsd.html__;!!JmPEgBY0HMszNaDT!5O_WiKmxcIblQ26zQKLIYV7Eg7QfEMgarJp3JpDss8ggtnmYbFSbMFNmzLk63QOas4I$ %3E<br class="">
<br class="">
There are examples in the links above and also in tutorials, e.g.<br class="">
<https://urldefense.com/v3/__http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/tutorials/images.html*matching__;Iw!!JmPEgBY0HMszNaDT!5O_WiKmxcIblQ26zQKLIYV7Eg7QfEMgarJp3JpDss8ggtnmYbFSbMFNmzLk6nOITpWA$ %3E<br class="">
<br class="">
If your two structures were #0 and #1 and they had same numbering, for example, you could match using backbone atoms of residues 10-300 from each with:<br class="">
<br class="">
match #0:10-300@n,ca,c,o #1:10-300@n,ca,c,o<br class="">
<br class="">
It will be harder if they have different numbering, and/or you wish to omit loop areas that are not well superimposed.<br class="">
<br class="">
Matchmaker vs. match is discussed here:<br class="">
<https://urldefense.com/v3/__http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/superposition.html__;!!JmPEgBY0HMszNaDT!5O_WiKmxcIblQ26zQKLIYV7Eg7QfEMgarJp3JpDss8ggtnmYbFSbMFNmzLk6Mw2g-s0$ %3E<br class="">
<br class="">
I hope this helps,<br class="">
ELaine<br class="">
-----<br class="">
Elaine C. Meng, Ph.D. <br class="">
UCSF Chimera(X) team<br class="">
Department of Pharmaceutical Chemistry<br class="">
University of California, San Francisco<br class="">
<br class="">
> On Sep 16, 2020, at 12:06 PM, Zhang, Bixia <bixia.zhang@wsu.edu> wrote:<br class="">
> <br class="">
> Hi,<br class="">
> I am trying to get the backbone RMSD between two structures. However I tried selecting backbone and use matchmaker, it gives me same value with the whole protein. So there must be something wrong. Could you please teach me a command that can do this? Thank
you so much for helping!<br class="">
> Best,<br class="">
> Bixia <br class="">
</blockquote>
<br class="">
</div>
</div>
<!--{"start":[1,1,0,12],"end":[1,1,0,12]}--></div>
<!--{"start":[0,1,1,0,12],"end":[0,1,1,0,12]}--></div>
<!--{"start":[0,0,1,1,0,70],"end":[0,0,1,1,0,70]}--></div>
<!--{"start":[0,0,0,1,1,0,70],"end":[0,0,0,1,1,0,70]}--></div>
<!--{"start":[0,0,0,0,1,1,0,80],"end":[0,0,0,0,1,1,0,80]}--></div>
<!--{"start":[0,0,0,0,0,1,1,0,80],"end":[0,0,0,0,0,1,1,0,80]}--></div>
</body>
</html>