<html>
  <head>
    <meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
  </head>
  <body>
    <div class="moz-cite-prefix">On 4/22/20 10:08 AM, Elaine Meng wrote:<br>
    </div>
    <blockquote type="cite"
      cite="mid:C4E6EC2A-3A17-41FE-9CFA-75277D6FD95F@cgl.ucsf.edu">
      <pre class="moz-quote-pre" wrap="">Hi Kenward,
Sorry for the difficulties.  Intersurf was an implementation based on a paper, and has turned out to be somewhat of a novelty item that doesn't get much use.  I suspect there are bugs, or quirks that I failed to explain in the documentation (likely not having understood it myself).  Its main distinguishing characteristic is that the "inter"surface is not the molecular surface of either set of atoms in the interface but derived from both sets of atoms.  I too am mystified by the behavior of the "prune" option in the command, as it seems that, at least in your example, all prune values (even 0) give the same thing.  I can make a bug report but we are mostly focused on ChimeraX these days.</pre>
    </blockquote>
    <p><br>
    </p>
    <p>Hi Elaine!<br>
    </p>
    <p>As is usual, thank you for the detailed reply with the several
      (detailed) suggestions!   <br>
    </p>
    <p>My goal is simply to give the students a visualization of the
      region between the antibody and the epitope itself, hoping that
      this would give them a better reference when looking at what
      appears to be a jungle of residues.  Intersurf struck me as the
      ideal, does it all at once, answer to this idea.  Another for me
      involves different reps for the two sides (thick line vs. ball and
      stick, etc.).<br>
    </p>
    <p>What appears to be working based on your suggestions is a cheap
      way to get those residues and a usable surface:</p>
    <p>split</p>
    <p>Intersurf B and K, selecting the atoms</p>
    <p>up arrow to collect all atoms of each residue, demo'd by Bonds
      > show<br>
    </p>
    <p>name the selection</p>
    <p>redo Intersurf, between B and M, reusing the surface.  Up arrow
      to gather all atoms of all residues again. Bonds > show.</p>
    <p>Trash the surface in model panel.<br>
    </p>
    <p>Select mode to append, then select the original named selection. 
      Now I have all residues for the interface.  New name for the
      selection.</p>
    <p>Intersect selection with chain B.</p>
    <p>Surf, colorize, transparency. See attached sample...<br>
    </p>
    <p><br>
    </p>
    <p>I'm unfamiliar with the buried area approach, so for quick
      results did the above despite the number of steps.  I'm learning
      all too much with our going online right now.  If there's a fast
      way to do this by command line I'm all for it, but ugh...<br>
    </p>
    <p>Don't know how well this'll work with them, but...  Still working
      on what colors would work well.<br>
    </p>
    <p>Thank you for the leading suggestions!  <br>
    </p>
    <p><br>
    </p>
    <p>Kenward</p>
    <p><br>
    </p>
    <blockquote type="cite"
      cite="mid:C4E6EC2A-3A17-41FE-9CFA-75277D6FD95F@cgl.ucsf.edu">
      <pre class="moz-quote-pre" wrap="">Depending on what you want the students to get from the exercise, there may be better alternatives, such as simply selecting atoms or residues with buried area in the interface between the two sets of atoms.  You could color those residues or their molecular surfaces.  (e.g. as in the "open book" image tutorial)
<a class="moz-txt-link-rfc2396E" href="http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/tutorials/openbook.html"><http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/tutorials/openbook.html></a>

Maybe the problem is that the default surface wraps around all the proteins together?  You can get separate ones by using "split" beforehand to make each chain a separate model (as in that tutorial above), or if you can keep the chains in the same model but use "surfcat" to wrap K and M separately from B, e.g. commands:

surfcat bee :.b
surfcat kandm :.k:.m
surf bee
surf kandm

Another command of possible interest is "measure buriedArea" which will assign buriedSESArea and buriedSASArea attribute values to the atoms.  Then one can select atoms by ranges of those values.
<a class="moz-txt-link-rfc2396E" href="http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/measure.html#buriedArea"><http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/measure.html#buriedArea></a>

measure buriedArea :.b :.k:.m
start Render by Attribute 
 (...that shows a dialog with histogram of values)
select @/buriedSASArea>2.0
(... or use the Select tab of that dialog to do the selection... whatever cutoff you like, of course)

Then you can "surf sel" or "color magenta sel" or whatever you were planning to do with the selection of interface atoms.
 
Sorry about the lame answer re Intersurf.  If there are other questions about how to accomplish various steps of the project, please don't hesitate to ask.  Best regards,
Elaine
-----
Elaine C. Meng, Ph.D.                       
UCSF Chimera(X) team
Department of Pharmaceutical Chemistry
University of California, San Francisco

</pre>
      <blockquote type="cite">
        <pre class="moz-quote-pre" wrap="">On Apr 21, 2020, at 8:49 PM, Kenward Vaughan <a class="moz-txt-link-rfc2396E" href="mailto:kay_jay@kvaughan.family"><kay_jay@kvaughan.family></a> wrote:

On 4/21/20 8:40 PM, Kenward Vaughan wrote:
...

</pre>
        <blockquote type="cite">
          <pre class="moz-quote-pre" wrap="">The command I use is

intersurf   :.b   :.k,.m   pair atom   prune 3   bias 0.5   select
true

followed by adjusting the transparency.  The result surprised me,
both before and after lighting up the atoms and bonds.  I've attached
a (hopefully small) jpg.

</pre>
        </blockquote>
        <pre class="moz-quote-pre" wrap="">
...

By the way, I just realized with further experimentation that a prune distance of 3 gives NO result ("no interface found") using the graphical interface with chain B and either K or M.  Changing this to 10 in my command still gives the same weird results...




Kenward
-- 
In a completely rational society, the best of us would aspire to be 
_teachers_ and the rest of us would have to settle for something less, 
because passing civilization along from one generation to the next 
ought to be the highest honor and the highest responsibility anyone 
could have.     - Lee Iacocca

Any fool can know. The point is to understand. - Albert Einstein

_______________________________________________
Chimera-users mailing list: <a class="moz-txt-link-abbreviated" href="mailto:Chimera-users@cgl.ucsf.edu">Chimera-users@cgl.ucsf.edu</a>
Manage subscription: <a class="moz-txt-link-freetext" href="http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users">http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users</a>
</pre>
      </blockquote>
      <pre class="moz-quote-pre" wrap="">
</pre>
    </blockquote>
    <p><br>
    </p>
    <pre class="moz-signature" cols="72">-- 
In a completely rational society, the best of us would aspire to be 
_teachers_ and the rest of us would have to settle for something less, 
because passing civilization along from one generation to the next 
ought to be the highest honor and the highest responsibility anyone 
could have.     - Lee Iacocca

Any fool can know. The point is to understand. - Albert Einstein
</pre>
  </body>
</html>