Thanks for the explanation, but still one more thing confuse me.<div><br></div><div>Indeed weather delete waters first then measure buriedArea, or just use measure buriedArea :.d&protein :.e&protein command gave BaveSAS 349.</div>
<div><br></div><div>But when I use surfcat, I mean I delete waters first.</div><div><br></div><div>then use command </div><div><br></div><div>surfcat 1 :.d</div><div>surf 1</div><div><br></div><div>measure surface area of d, I got 2972.6, </div>
<div><br></div><div>I do the same thing for e, </div><div><br></div><div><div>surfcat 2 :.e</div><div>surf 2</div></div><div><br></div><div>I got 5504; </div><div><br></div><div><br></div><div>then I do that for d,e as one model,</div>
<div><br></div><div><div>surfcat 3 :.d-e</div><div>surf 3</div><div><br></div><div>I got 8117.6.</div><div><br></div><div>From these values I do the calculation (2972.6+5504-8117.6)/2=179.5.</div><div><br></div><div>This value does not match with BaveSAS 349.4, or BaveSES 186.2.</div>
<div><br></div><div>Do you know where the difference come from? Does chimera define different surfaces based on different methods?</div><div><br></div><div>thanks a lot </div><div><br></div><div>bobo </div><div><br><div class="gmail_quote">
On Thu, Jan 10, 2013 at 1:54 PM, Elaine Meng <span dir="ltr"><<a href="mailto:meng@cgl.ucsf.edu" target="_blank">meng@cgl.ucsf.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi Bobo,<br>
One issue to be careful about with both "surfcat" and "measure buriedArea" is to NOT include stuff like waters.  I forgot to mention that issue in my first reply.  You could either delete water first, or when using the commands take care to use specifications that exclude water, something like this:<br>

<br>
measure buriedArea :.d&protein :.e&protein<br>
<br>
(because there are waters with chain ID D and waters with chain ID E)<br>
<br>
When I did that, I got the same answer both ways, BaveSAS 349 for 4gln protein-only.<br>
<div class="im">I hope this helps,<br>
Elaine<br>
<br>
</div><div><div class="h5">On Jan 10, 2013, at 11:14 AM, Bobo Dang wrote:<br>
<br>
> Thanks a lot.<br>
><br>
> Your method B is exactly what I want to get.<br>
><br>
> As for the binding area, it seems that using different methods, I get different values.<br>
><br>
> I tried surfcat command to get the surface area of D, E and DE. and then calculate by (D+E-DE)/2, it gave me 363. but using "measure buriedArea" command it gave me 372 for BaveSAS. Do you know why these two results are not the same even though the difference is small?<br>

><br>
><br>
> thanks<br>
><br>
> bobo<br>
><br>
> On Thu, Jan 10, 2013 at 11:59 AM, Elaine Meng <<a href="mailto:meng@cgl.ucsf.edu">meng@cgl.ucsf.edu</a>> wrote:<br>
> Hi Bobo,<br>
> (A) If you mostly want the area measurement (not display), you can just use the command "measure buriedArea."  Use the Favorites menu to show the Chimera command line.  Structure 4GLN has 4 chains: D, E, F, and H ... I could see that by rainbow-coloring each chain a different color, for example with Chimera commands:<br>

><br>
> open 4gln<br>
> rainbow chain<br>
><br>
> You can see which chain is which by putting the mouse cursor over them, which will show an information "balloon."  So for example, I could measure area between chains E and F with:<br>
><br>
> measure buriedArea :.e  :.f<br>
><br>
> ...or between D+E and F+H with:<br>
><br>
> measure buriedArea :.d:.e  :.f:.h<br>
><br>
> The information goes to the Reply Log (show using Favorites menu).  For details on this command, see:<br>
> <<a href="http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/measure.html#buriedArea" target="_blank">http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/measure.html#buriedArea</a>><br>
><br>
> (B) If you care about display too, then more steps are required. By default, Chimera will make one surface that includes all the protein chains.  The "split" command will put each chain in a separate model, for example:<br>

><br>
> split<br>
> surface<br>
> rainbow model<br>
><br>
> <<a href="http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/split.html" target="_blank">http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/split.html</a>><br>
><br>
> You can see what models you have by showing the Model Panel (also in Favorites menu).<br>
><br>
> Then you still want to show just the interface.  For example, to show only surface of E within 1.5 Angstroms of surface of F:<br>
><br>
> ~surface<br>
> measure contactArea #0.2 #0.3 1.5 color orange offset 0.1<br>
><br>
> <<a href="http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/measure.html#contactArea" target="_blank">http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/measure.html#contactArea</a>><br>
><br>
> I had to look in the Model Panel to see what model numbers to use for the "measure contactArea" command.  If offset is >0, this command will make a new model, also listed in the Model Panel. To also show the whole surface of E, for example, you could use the command:<br>

><br>
> surface :.e<br>
><br>
> (C) instead of B, a completely different approach is to use the Intersurf tool to show an interface surface instead of the molecular surface of either chain.  This tool is in the menu under Tools... Surface/Binding Analysis.  In that dialog you would choose the two chains and then Apply.   You can use the histogram at the bottom to color by distance. You could also use the command "intersurf" to do the same thing.<br>

> <<a href="http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/intersurf.html" target="_blank">http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/intersurf.html</a>><br>
> <<a href="http://www.cgl.ucsf.edu/chimera/docs/ContributedSoftware/intersurf/intersurf.html" target="_blank">http://www.cgl.ucsf.edu/chimera/docs/ContributedSoftware/intersurf/intersurf.html</a>><br>
><br>
> Example Intersurf result is shown here:<br>
> <<a href="http://www.cgl.ucsf.edu/chimera/ImageGallery/entries/alphabeta/alphabeta.html" target="_blank">http://www.cgl.ucsf.edu/chimera/ImageGallery/entries/alphabeta/alphabeta.html</a>><br>
><br>
> The Intersurf results are more "chunky," not as smooth as the molecular surface(s) shown by method B above.<br>
><br>
> I hope this helps,<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>
><br>
> On Jan 10, 2013, at 8:21 AM, Bobo Dang wrote:<br>
><br>
> > Dear Chimera Users,<br>
> ><br>
> > I am trying to use Chimera to show the binding surface between two chains in one Model and measure the area of this binding surface,<br>
> ><br>
> > but I don't know how I can select the bind surface, could anybody give me a hand? the model I am using is 4GLN from PDB.<br>
> ><br>
> > thanks a lot<br>
> ><br>
> > bobo<br>
><br>
><br>
</div></div>> _______________________________________________<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" target="_blank">http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users</a><br>
<br>
</blockquote></div><br></div></div>