<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=Windows-1252">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body dir="ltr">
<div style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
Steven,</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
rmsdAlign is not the same command - it is from my SEQCROW bundle. rmsdAlign does not return the RMSD value (yet). It also uses the entire model, not subsets of models. I will likely change this at some point.</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
When I use "rmsd" on the command line, it uses the built-in rmsd command. I need to type out "rmsdA" at least to use rmsdAlign. If "rmsd" defaults to rmsdAlign for you, I can rename rmsdAlign in a quick patch.</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
Best,</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
Tony<br>
</div>
<div id="appendonsend"></div>
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="divRplyFwdMsg" dir="ltr"><font style="font-size:11pt" face="Calibri, sans-serif" color="#000000"><b>From:</b> ChimeraX-users <chimerax-users-bounces@cgl.ucsf.edu> on behalf of Steven Truong via ChimeraX-users <chimerax-users@cgl.ucsf.edu><br>
<b>Sent:</b> Tuesday, October 12, 2021 5:39 PM<br>
<b>To:</b> Eric Pettersen <pett@cgl.ucsf.edu><br>
<b>Cc:</b> ChimeraX Users Help <chimerax-users@cgl.ucsf.edu><br>
<b>Subject:</b> Re: [chimerax-users] Use values printed from Log inside a script (e.g. RMSD)</font>
<div> </div>
</div>
<div class="" style="word-wrap:break-word; line-break:after-white-space"><font color="BA0C2F">[EXTERNAL SENDER - PROCEED CAUTIOUSLY]</font><br>
<br>
<div>Hi Eric, Tony, et al.,
<div class=""><br class="">
</div>
<div class="">Thank you for answering my question. For now, I’ll be doing what Tony suggested re: assigning a value with the RMSD command. Currently, I’m having issues with using “rmsd” by itself (it seems like rmsdAlign is the only valid command). I’m assuming
this is the same command, just named differently.</div>
<div class=""><br class="">
</div>
<div class="">Re: your inquiry, the protein does have the same sequences. However, the chain ID is critical, since I am attempting to “stitch” fragments of proteins at different residues. For example, if I align a fragment to the wrong chain, it will have
a bond hanging on the other side of the protein (not good!). This stitching can be done manually (e.g. if it’s not the right chain or residue ID, I can simply move the fragment), but I am trying to automate this process via Python scripts. Additionally,
I have been doing the bond-forming on Chimera (with the command “bond”). I’m not sure if there’s a better way of going about this, but my idea currently is to 1. Identify which chains match up, 2. Align fragments from 1 protein to the other (which may have
different chain ID identifiers), 3. Make bonds on Chimera.</div>
<div class=""><br class="">
</div>
<div class="">Many thanks,</div>
<div class="">Steven <br class="">
<div class="">
<div><br class="">
<blockquote type="cite" class="">
<div class="">On Oct 12, 2021, at 9:31 PM, Eric Pettersen <<a href="mailto:pett@cgl.ucsf.edu" class="">pett@cgl.ucsf.edu</a>> wrote:</div>
<br class="x_Apple-interchange-newline">
<div class="">
<div class="" style="word-wrap:break-word; line-break:after-white-space">Hi Steven,
<div class=""><br class="">
<div class="">
<blockquote type="cite" class="">
<div class="">On Oct 12, 2021, at 12:39 PM, Steven Truong via ChimeraX-users <<a href="mailto:chimerax-users@cgl.ucsf.edu" class="">chimerax-users@cgl.ucsf.edu</a>> wrote:</div>
<br class="x_Apple-interchange-newline">
<div class="">
<div class="">Dear ChimeraX Admins,<br class="">
<br class="">
I was hoping to use the “rmsd” command in ChimeraX to measure distances between models within a .py script.
<br class="">
<br class="">
Is there a way to use the output of “rmsd” (or generally any other output in Log) within the script? For example, if rmsd > 2 Angstroms, then do X.
<br class="">
</div>
</div>
</blockquote>
<div class=""><br class="">
</div>
Tony answered this.</div>
<div class=""><br class="">
<blockquote type="cite" class="">
<div class="">
<div class="">Ultimately, I am attempting to align chains between models. Because my protein is asymmetric, chain A on Model #1 might actually be chain B on Model #2. I was hoping to use “rmsd” to map one chain to the other. Is it a better idea to do this in
another program and re-write chain IDs in the .PDB file? If so, what other program would be better at making RMSD calculations between .PDB chains?<br class="">
</div>
</div>
</blockquote>
<div class=""><br class="">
</div>
Do these chains have the same sequence? Also, I'm a little uncertain what "map one chain to the other" means exactly. At any rate, the
<a href="https://www.cgl.ucsf.edu/chimerax/docs/user/commands/matchmaker.html" class="">
matchmaker</a> command can be used to superimpose one chain onto another and report the backbone (CA or C4'/P) RMSD for the match. By default matchmaker iterates its match, pruning off poorly matching regions and re-matching to get a good fitting "core" match,
and will report RMSD values for both the core region and across the entire sequence. You can turn off the iteration if desired with matchmaker's "cutoffDistance none" option.</div>
<div class=""><br class="">
</div>
<div class="">
<div class="">--Eric</div>
<div class=""><br class="">
</div>
<div class=""><span class="x_Apple-tab-span" style="white-space:pre"></span>Eric Pettersen</div>
<div class=""><span class="x_Apple-tab-span" style="white-space:pre"></span>UCSF Computer Graphics Lab</div>
<div class=""><br class="">
</div>
</div>
<div class="">
<blockquote type="cite" class="">
<div class="">
<div class=""><br class="">
Thank you,<br class="">
Steven Truong<br class="">
Cambridge University<br class="">
<br class="">
_______________________________________________<br class="">
ChimeraX-users mailing list<br class="">
<a href="mailto:ChimeraX-users@cgl.ucsf.edu" class="">ChimeraX-users@cgl.ucsf.edu</a><br class="">
Manage subscription:<br class="">
<a href="https://www.rbvi.ucsf.edu/mailman/listinfo/chimerax-users" class="">https://www.rbvi.ucsf.edu/mailman/listinfo/chimerax-users</a><br class="">
<br class="">
</div>
</div>
</blockquote>
</div>
<br class="">
</div>
</div>
</div>
</blockquote>
</div>
<br class="">
</div>
</div>
</div>
</div>
</body>
</html>