<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=Windows-1252">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; color: rgb(0, 0, 0); font-size: 14px; font-family: Calibri, sans-serif;">
<div>
<div>Dear Alex,</div>
<div>
<div><br>
</div>
<div>Instead of calling chimera from a python script, you can simply use chimera's python to run the script and call the command line commands with the runCommand() function. This way, you can have just the python script and avoid the cmd dependency.</div>
<div><br>
</div>
<div>You can run the following script with "chimera --nogui script.py":</div>
<div><br>
</div>
<div>"""</div>
<div>
<div>import chimera</div>
<div>from chimera import runCommand as rc</div>
<div>rc('open 1BPT')</div>
<div>rc('swapaa tyr :5; swapaa phe :7')</div>
</div>
<div>"""</div>
<div><br>
</div>
<div>Best,</div>
<div>Matej</div>
<div><br>
</div>
<div><span style="font-family: Consolas; font-size: medium;">---------------------------------</span></div>
<div>
<div style="font-family: Consolas; font-size: medium;">Dr. Matej Repic</div>
<div>
<div><font face="Consolas">Ecole Polytechnique Fédérale de Lausanne</font></div>
<div><font face="Consolas">Laboratory of Computational Chemistry and Biochemistry</font></div>
<div><font face="Consolas">SB - ISIC – LCBC</font></div>
<div><span style="font-family: Consolas;">BCH 4108</span></div>
<div><font face="Consolas">CH - 1015 Lausanne</font></div>
</div>
<div style="font-family: Consolas; font-size: medium;">------------------------------------------------------</div>
</div>
</div>
</div>
<div><br>
</div>
<span id="OLK_SRC_BODY_SECTION">
<div style="font-family:Calibri; font-size:11pt; text-align:left; color:black; BORDER-BOTTOM: medium none; BORDER-LEFT: medium none; PADDING-BOTTOM: 0in; PADDING-LEFT: 0in; PADDING-RIGHT: 0in; BORDER-TOP: #b5c4df 1pt solid; BORDER-RIGHT: medium none; PADDING-TOP: 3pt">
<span style="font-weight:bold">From: </span><<a href="mailto:chimera-users-bounces@cgl.ucsf.edu">chimera-users-bounces@cgl.ucsf.edu</a>> on behalf of roz orchidea <<a href="mailto:roz.orchidea@gmail.com">roz.orchidea@gmail.com</a>><br>
<span style="font-weight:bold">Date: </span>Thursday, July 9, 2015 at 12:49<br>
<span style="font-weight:bold">To: </span>"<a href="mailto:chimera-users@cgl.ucsf.edu">chimera-users@cgl.ucsf.edu</a>" <<a href="mailto:chimera-users@cgl.ucsf.edu">chimera-users@cgl.ucsf.edu</a>><br>
<span style="font-weight:bold">Subject: </span>Re: [Chimera-users] swapaa: choosing rotamers<br>
</div>
<div><br>
</div>
<div>
<div>
<div dir="ltr">
<div>
<div>Thank you both for your answers! <br>
Eric I am not trying all rotamers, I just want to have some specific rotamers in the binding sites of my protein-DNA models. I am using a python script to call chimera and a cmd script with all the swapaa commands that I need. Is there another way?<br>
</div>
<br>
</div>
<div>Thank you for your time,<br>
</div>
<div>Alex<br>
</div>
</div>
<div class="gmail_extra"><br>
<div class="gmail_quote">2015-07-01 19:03 GMT+02:00 Eric Pettersen <span dir="ltr">
<<a href="mailto:pett@cgl.ucsf.edu" target="_blank">pett@cgl.ucsf.edu</a>></span>:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Also, the "lib" keyword to the swapaa command allows you to pick the library to use, e.g.: swapaa ASP :41.a lib dynameomics criteria 2<br>
<span class="HOEnZb"><font color="#888888"><br>
--Eric<br>
</font></span>
<div class="HOEnZb">
<div class="h5"><br>
On Jul 1, 2015, at 10:01 AM, Eric Pettersen <<a href="mailto:pett@cgl.ucsf.EDU">pett@cgl.ucsf.EDU</a>> wrote:<br>
<br>
> Hi Alex,<br>
> To get specifically get the second most probable rotamer you add the "criteria" keyword to the swapaa command with the value "2", e.g.: swapaa ASP :41.a criteria 2<br>
> You don't really say why you specifically want the second rotamer. Are you trying out all the rotamers one by one? You might want to consider doing that via the a Python script that uses the functions in the Rotamers module. That way you can find
out how many rotamers are available, which isn't really possible with just the swapaa command.<br>
> Another thing to consider is that if you know the chi angles you want, you can just swap to the new residue type and then request those specific chi values. For example to change residue 41.A to ASP and set its chis to 40 and 75:<br>
><br>
> swapaa ASP :41.a<br>
> setattr a chi1 40 :41.a<br>
> setattr a chi2 75 :41.a<br>
><br>
> --Eric<br>
><br>
> Eric Pettersen<br>
> UCSF Computer Graphics Lab<br>
> <a href="http://www.cgl.ucsf.edu" rel="noreferrer" target="_blank">
http://www.cgl.ucsf.edu</a><br>
><br>
> On Jul 1, 2015, at 9:32 AM, Elaine Meng <<a href="mailto:meng@cgl.ucsf.edu">meng@cgl.ucsf.edu</a>> wrote:<br>
><br>
>> Dear Alex,<br>
>> Not directly, the “swapaa” command only uses the ranking rules that you tell it to use. Only if the rules happened to choose that rotamer would you get it. To directly choose a specific rotamer as listed in the rotamers dialog, you would have to use that
dialog (Rotamers graphical user interface).<br>
>><br>
>> The rules are specified with the “criteria” option of swapaa, see:<br>
>> <<a href="http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/swapaa.html#criteria" rel="noreferrer" target="_blank">http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/swapaa.html#criteria</a>><br>
>><br>
>> Rotamers graphical interface:<br>
>> <<a href="http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/rotamers/framerot.html" rel="noreferrer" target="_blank">http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/rotamers/framerot.html</a>><br>
>><br>
>> Best,<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>
>> On Jul 1, 2015, at 4:57 AM, roz orchidea <<a href="mailto:roz.orchidea@gmail.com">roz.orchidea@gmail.com</a>> wrote:<br>
>><br>
>>> Dear all,<br>
>>> I would like to use swapaa and choose which rotamer I want. Foe example I want from library dynameomics the second rotamer that appears in the pop-up menu of Arg. is it possible?<br>
>>> Thank you for your time,<br>
>>> Alex<br>
>><br>
>><br>
>> _______________________________________________<br>
>> Chimera-users mailing list<br>
>> <a href="mailto:Chimera-users@cgl.ucsf.edu">Chimera-users@cgl.ucsf.edu</a><br>
>> <a href="http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users" rel="noreferrer" target="_blank">
http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users</a><br>
>><br>
><br>
><br>
><br>
<br>
</div>
</div>
</blockquote>
</div>
<br>
</div>
</div>
</div>
</span>
</body>
</html>