<html><head><meta http-equiv="Content-Type" content="text/html; charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">Great!<div class=""><br class=""></div><div class="">—Eric<br class=""><div><br class=""><blockquote type="cite" class=""><div class="">On Sep 9, 2019, at 6:31 AM, Kamalesh Damodaran <<a href="mailto:kamaleshdams@gmail.com" class="">kamaleshdams@gmail.com</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div dir="ltr" class=""><div class="">Hi Eric,</div><div class=""><br class=""></div><div class="">Thank you so much for your guidance and the python script. I got everything what I was looking for. <br class=""></div><div class=""><br class=""></div><div class="">Thanks and regards,</div><div class="">Kamalesh D.<br class=""></div></div><br class=""><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, 4 Sep 2019 at 04:54, Eric Pettersen <<a href="mailto:pett@cgl.ucsf.edu" class="">pett@cgl.ucsf.edu</a>> wrote:<br class=""></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div style="overflow-wrap: break-word;" class="">Hi Kamalesh,<div class=""><span class="gmail-m_-8897730003222345391Apple-tab-span" style="white-space:pre-wrap">     </span>If you are familiar with Python, there is a script attached to this previous mailing-list message:  <a href="http://www.cgl.ucsf.edu/pipermail/chimera-users/2007-December/002090.html" target="_blank" class="">[Chimera-users] Writing alignment and RMSD out</a> that would be easy to modify to do what you want.  You would just remove the parts the parts that used Multalign Viewer to write the sequence alignment.  The RMSD that the MatchMaker.match() function returns is a carbon-alpha RMSD, so if you needed a full-backbone RMSD that would be more work.  Let me know and I can offer guidance on that.  Also, match() returns the overall RMSD of the match.  To get the per-residue RMSD you would loop over the atom lists returned by the match() function and for corresponding pairs get the distance (== RMSD for two atoms) with:  atom1.xformCoord().distance(atom2.xformCoord()).  For output purposes you can get a string representation of an atom’s residue with str(atom.residue).</div><div class=""><span class="gmail-m_-8897730003222345391Apple-tab-span" style="white-space:pre-wrap">     </span>Also, unfortunately the attachment to the message has a “.bin” suffix despite being from a Python file that had a “.py” suffix, so if that’s a problem let me know and I can send you the .py file directly.</div><div class=""><br class=""></div><div class="">—Eric<br class=""><div class="">
<div style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none;" class=""><span class="gmail-m_-8897730003222345391Apple-tab-span" style="white-space:pre-wrap"><br class="gmail-m_-8897730003222345391Apple-interchange-newline">       </span>Eric Pettersen</div><div style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none;" class=""><span class="gmail-m_-8897730003222345391Apple-tab-span" style="white-space:pre-wrap">  </span>UCSF Computer Graphics Lab</div><div style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none;" class=""><br class=""></div><br class="gmail-m_-8897730003222345391Apple-interchange-newline">
</div>
<div class=""><br class=""><blockquote type="cite" class=""><div class="">On Aug 31, 2019, at 7:20 AM, Kamalesh Damodaran <<a href="mailto:kamaleshdams@gmail.com" target="_blank" class="">kamaleshdams@gmail.com</a>> wrote:</div><br class="gmail-m_-8897730003222345391Apple-interchange-newline"><div class=""><div dir="ltr" class=""><div class="">Dear sir/madam,</div><div class=""><br class=""></div><div class="">I would like to superimpose the two structures and compute the residue wise RMSD values. With Gui support in chimera, I opened the two structures and used Matchmaker to superimpose two structures. In Matchmaker, I checked the box 'Show pairwise alignment(s)' and click Apply, subsequently it displays MultAlignViewer dialogue box, where --> Headers --> RMSD : backbone --> Save --> filenamegiven. This process would give me residue wise RMSD score only for backbone atoms in .hdr format. <br class=""></div><div class=""><br class=""></div><div class="">I wants to write python script to do the same, so that I will compare large number of structures, and my script is able to perform till the procedure of superimposition using matchmaker and display MultAlignViewer dialogue box. I don't know the way to store the RMSD : backbone and RMSD : calpha in .hdr format using python script through chimera.</div><div class=""><br class=""></div><div class="">Can some one please help me in this. The reason why I specifically want to do in this approach is, because among the two structures, some residues of one structure is missing. Since this approach will compute RMSD only for the matching residues, I am interested in this. <br class=""></div><div class=""><br class=""></div><div class="">Hope my question is not confusing. Thanks in advance.</div><div class=""><br class=""></div><div class="">Best regards,</div><div class="">Kamalesh D.<br class=""></div></div>
_______________________________________________<br class="">Chimera-users mailing list: <a href="mailto:Chimera-users@cgl.ucsf.edu" target="_blank" class="">Chimera-users@cgl.ucsf.edu</a><br class="">Manage subscription: <a href="http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users" target="_blank" class="">http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users</a><br class=""></div></blockquote></div><br class=""></div></div></blockquote></div>
</div></blockquote></div><br class=""></div></body></html>