<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body dir="ltr">
<div class="elementToProof"><span style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);">Hi Elaine,
<br>
Thank you a lot for this clarification, my group and I will reconsider the parameters to use for the rotamers selection, because they are indeed really interesting.
<br>
</span></div>
<div class="elementToProof"><span style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);">One thing though, and it's just a semantic question out of curiosity, but how would you
 write the option to run swapaa on the first rotamer depending on the probability only ? It would be something comprising "criteria", "1" for the rank and "p" for the parameter concerned, but I don't know how to take the 2 last into account at the same time.<br>
Thank you again for your time, and best wishes for 2023 !</span></div>
<div class="elementToProof"><span style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);">Robert Noé<br>
</span></div>
<div id="appendonsend"></div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>De :</b> Elaine Meng <meng@cgl.ucsf.edu><br>
<b>Envoyé :</b> lundi 2 janvier 2023 18:48<br>
<b>À :</b> Noe Robert <noe.robert@etu.univ-cotedazur.fr><br>
<b>Cc :</b> chimera-users@cgl.ucsf.edu <chimera-users@cgl.ucsf.edu><br>
<b>Objet :</b> Re: [Chimera-users] Question concerning rotamer choosing</font>
<div> </div>
</div>
<div class="BodyFragment"><font size="2"><span style="font-size:11pt;">
<div class="PlainText">Hi Noé Robert,<br>
I can't answer the python part, but wanted to make sure that you understand the meaning of the probability: it is the probability based on backbone conformation only.  It is NOT the total probability for the specific position, because it does not look at sidechain
 interactions.  So the one with the highest probability might be a very poor choice for the specific position (might have high clashes, low H-bonds). 
<br>
<br>
It is a common misconception that the highest probability is the best one to choose.<br>
<br>
Rotamers mutagenesis is implemented as the "swapaa" command where the default is to take account first clashes, then H-bonds, the probability (again because that value is only based on the backbone conformation).   See the "swapaa" command "criteria" option:<br>
<<a href="https://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/swapaa.html#options">https://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/swapaa.html#options</a>><br>
<br>
So even though you may still want to use python for various reasons like looping, you should consider carefully how you choose the rotamer, and possibly use the default settings of the "swapaa" command.  Even if you decide you only want to use probability,
 then you can still use the "swapaa" command option with "criteria p" as explained in the help link above.<br>
<br>
I hope this clarifies the situation,<br>
Elaine<br>
-----<br>
Elaine C. Meng, Ph.D.                       <br>
UCSF Chimera(X) team<br>
Department of Pharmaceutical Chemistry<br>
University of California, San Francisco<br>
<br>
> On Jan 2, 2023, at 6:26 AM, Noe Robert via Chimera-users <chimera-users@cgl.ucsf.edu> wrote:<br>
> <br>
> Hi, <br>
> I'm currently trying to write a python script which automatizes mutageneses on Chimera 1.16, and calculates clashes inducted by the change of amino acid through a local minimization (which explains why I'm using Chimera instead of ChimeraX).
<br>
> I would like to choose which rotamer to use depending on their probability from the rotamer library Dunbrack (for instance, when probability is >0.1), which is visible on the "Rotamers graphical user interface". Is there a way to extract these probability
 values using a script, in order to run it on a certain number of rotamers ? That way it would be quite automated.<br>
> On Jul 1, 2015, Eric Pettersen advised to modify the "Python script that uses the functions in the Rotamers module" (here), but I don't which one it is.<br>
> <br>
> Thank you,<br>
> <br>
> Noé Robert,<br>
> Fifth-year student in bioengineering, minor in Bioinformatics and Modelling, at engineering school Polytech Nice Sophia, France<br>
> Phone : +33 7 81 00 52 02<br>
> "let's limit the use of our resources and lighten our mails"<br>
> <br>
> <br>
> _______________________________________________<br>
> Chimera-users mailing list: Chimera-users@cgl.ucsf.edu<br>
> Manage subscription: <a href="https://www.rbvi.ucsf.edu/mailman/listinfo/chimera-users">
https://www.rbvi.ucsf.edu/mailman/listinfo/chimera-users</a><br>
<br>
</div>
</span></font></div>
</body>
</html>