<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<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);">
Hi Tom,</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);">
Yes, "kappa" is the stiffness. Default in the current release is 5, (as specified in the documentation via "usage isolde restrain") which I agree is very soft for this particular purpose - after playing around with it a bit I find "isolde restrain distances
#1 distance 5 kappa 100 fallOff -3" works quite well. The extra "adjustForConfidence true" argument coming in the new ISOLDE release (not yet available - here in a few weeks) will help further - because it uses AlphaFold's own confidence on residue-residue
distances, it can extend the distance restraints over a longer range (e.g. the default 8 Angstroms) while ignoring spurious chance contacts between residues. Will play around a little more and share the result.</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);">
-- Tristan</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>From:</b> Tom Goddard <goddard@sonic.net><br>
<b>Sent:</b> 25 October 2021 20:01<br>
<b>To:</b> chimerax-users@cgl.ucsf.edu <chimerax-users@cgl.ucsf.edu>; Roden Deng Luo <deng.luo@kaust.edu.sa>; Tristan Croll <tic20@cam.ac.uk><br>
<b>Subject:</b> Re: [chimerax-users] Refine predicted structure by pulling</font>
<div> </div>
</div>
<div>Hi Roden, Tristan,<br>
<br>
I tried Tristan's suggestion for straightening the long AlphaFold helix model. It didn't work well for me because the restraints were too soft, so if I pulled on the first 160 residues it hardly moved the attached 2000 and instead just distorted the linker.
Tristan is there a way to make the restraints much stiffer? I see the "isolde restrain distance" command has a "kappa" option and "falloff" option but I don't know if kappa is a stiffness or what magnitude it usually uses.<br>
<br>
I think a better approach is to just rigidly rotate and translate large straight sections by hand using the "Move atoms" mouse mode in the ChimeraX 1.3 toolbar tab called Right Mouse. I tried this and have attached an image of about 10 minutes fiddling by
hand. I hover the mouse at a hinge residue to get the residue number and then select everything on one side of that hinge "select :1182-2419" then use the right mouse button in the "move atoms" mode, holding the shift key down to do rotations, and no shift
key to do translations. I zoom in on the hinge to try to adjust so the bond is not stretched and to avoid clashes. It would probably take 1 hour to do a good job. With my 10 minute test I then saved an mmCIF file and tried in ChimeraX 1.2.5 to minimize
it with ISOLDE, but it said there were severe clashes.<br>
<br>
Tom<br>
<br>
<br>
<br>
<div class="x_moz-cite-prefix">On 10/24/2021 1:49 AM, Tristan Croll via ChimeraX-users wrote:<br>
</div>
<blockquote type="cite"><style type="text/css" style="display:none">
<!--
p
{margin-top:0;
margin-bottom:0}
-->
</style>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
Hi Roden,</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)">
This should be reasonably straightforward with the ISOLDE plugin (<a href="https://cxtoolshed.rbvi.ucsf.edu/apps/chimeraxisolde" id="LPNoLPOWALinkPreview">https://cxtoolshed.rbvi.ucsf.edu/apps/chimeraxisolde</a>). You'll also want to be working on a machine
with a good GPU to make it feasible for a model this size. It'll be easier still with the new release coming in a few weeks, which has some added tools specifically for working with AlphaFold models. But with care you could do it now:</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)">
<ul>
<li><span>first, check carefully for entanglement between distant domains (this is an issue for AlphaFold with these very long, flexible proteins - if distant parts of the protein have no coevolutionary signal they don't really "know" about each other and can
end up trying to occupy the same space). Looks like your particular example is free of that problem.</span></li><li><span>Start ISOLDE. Then, assuming your protein is model #1, do the following commands:</span></li><ul>
<li><span>isolde restrain distances #1 distanceCutoff 5 </span></li><li><span>isolde restrain torsions #1 angleRange 60</span></li><li>isolde sim start #1</li></ul>
<li>Then, you'll want to use the "tug selection" mouse mode (third from bottom right of the ISOLDE panel). Make a selection for a suitable region to tug (e.g. "sel #1-100@CA" to tug the N-terminus), and tugging forces applied by right-click-and-drag will be
spread over that selection. Keep a close eye on the distance restraints - if they start to become too strained (lots of purple appearing), stop tugging and let the model settle.</li><li>If you have unwanted distance restraints between domains that should be distant from each other, you can selectively release restraints on selected atoms with "isolde release distances sel".</li><li>You might find it easier for a job like this to reduce the view to a C-alpha trace - just "hide ~@CA" - if you do this after you've started a simulation, it'll automatically revert to the all-atom view once you stop.</li></ul>
<div>Happy modelling!</div>
<div><br>
</div>
<div>-- Tristan</div>
</div>
<br>
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="x_divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>From:</b> ChimeraX-users
<a class="x_moz-txt-link-rfc2396E" href="mailto:chimerax-users-bounces@cgl.ucsf.edu">
<chimerax-users-bounces@cgl.ucsf.edu></a> on behalf of Roden Deng Luo via ChimeraX-users
<a class="x_moz-txt-link-rfc2396E" href="mailto:chimerax-users@cgl.ucsf.edu"><chimerax-users@cgl.ucsf.edu></a><br>
<b>Sent:</b> 23 October 2021 16:30<br>
<b>To:</b> <a class="x_moz-txt-link-abbreviated" href="mailto:chimerax-users@cgl.ucsf.edu">
chimerax-users@cgl.ucsf.edu</a> <a class="x_moz-txt-link-rfc2396E" href="mailto:chimerax-users@cgl.ucsf.edu">
<chimerax-users@cgl.ucsf.edu></a><br>
<b>Subject:</b> [chimerax-users] Refine predicted structure by pulling</font>
<div> </div>
</div>
<div>
<div dir="ltr">Dear ChimeraX Users,<br>
<div><br>
</div>
<div>There is this AlphaFold prediction: <a href="https://alphafold.ebi.ac.uk/entry/P02549">https://alphafold.ebi.ac.uk/entry/P02549</a>. And we know this protein should be more or less in a rod-like shape rather than a globular shape:
<a href="https://en.wikipedia.org/wiki/Spectrin">https://en.wikipedia.org/wiki/Spectrin</a> and related literature such as
<a href="https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3218374/">https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3218374/</a>. I wonder is there a way that we can pull the two ends of the predicted structure while maintaining the secondary structures, mostly alpha-helices
in this case, inside ChimeraX? </div>
<div><br>
</div>
<div>I am quite new to ChimeraX. I understand this is more fitted to an MD simulation. But I got to know there are options to fit a PDB to a density map. So I wondered if this can be achieved in ChimeraX. </div>
<div><br>
</div>
<div>Best,</div>
<div>Roden</div>
</div>
<br>
<div>
<hr>
</div>
<font size="1" face="Arial">This message and its contents, including attachments are intended solely for the original recipient. If you are not the intended recipient or have received this message in error, please notify me immediately and delete this message
from your computer system. Any unauthorized use or distribution is prohibited. Please consider the environment before printing this email.</font></div>
<br>
<fieldset class="x_mimeAttachmentHeader"></fieldset>
<pre class="x_moz-quote-pre">_______________________________________________
ChimeraX-users mailing list
<a class="x_moz-txt-link-abbreviated" href="mailto:ChimeraX-users@cgl.ucsf.edu">ChimeraX-users@cgl.ucsf.edu</a>
Manage subscription:
<a class="x_moz-txt-link-freetext" href="https://www.rbvi.ucsf.edu/mailman/listinfo/chimerax-users">https://www.rbvi.ucsf.edu/mailman/listinfo/chimerax-users</a>
</pre>
</blockquote>
<br>
</div>
</body>
</html>