[Chimera-users] Calculating ASA for side chain

Elaine Meng meng at cgl.ucsf.edu
Fri Jun 27 09:46:11 PDT 2014


Dear Saad,
There is no need to modify: even though only the total is reported in the Reply Log, the surface area values for each atom and residue are also calculated and stored as attributes, and you can add up the values for any set of atoms that you like.

When you show a molecular surface, 
"Analytical solvent-excluded and solvent-accessible surface areas per atom and residue are assigned as attributes named areaSES and areaSAS, respectively."
<http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/representation.html#surfaces>
<http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/defineattrib/defineattrib.html#attribdef>

You can sum the atomic values over any set of atoms of interest using the Attribute Calculator tool:
<http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/calculator/calculator.html>

For example, after showing the surface (you can hide it so you can see the atoms if they are displayed):
(A) select the atoms of interest, e.g. command: select :35.A & without CA/C1'
... would select side chain atoms of residue 35 chain A, not including CA atom (use "with" instead of "without" to include CA atom)
(B) use Attribute Calculator  (in menu under Tools... Structure Analysis) to calculate a new attribute of *residues* named whatever you want, e.g. sidechainSAS, with formula:
sum(atom.areaSAS)
...and also use the option to "Restrict formula domain to current selection" and "Open Render/Select by Attribute"

The Render/Select by Attribute tool will show the result, the sum of areaSAS over the selected atoms.  Since this would define a new residue attribute with a different value for each residue, you could select the side chains of many sidechains at once (e.g. command: select protein & without CA/C1') and this process would give a separate total for each residue.  To find out the values of your new residue attribute, you could write them out as a list using Render/Select by Attribute's menu: File... Save Attributes, or show them with labels, see main Chimera menu: Actions... Label... residue... custom, if using my example name, with Residue label: %(sidechainSAS).2f

The list I wrote for my test structure looked like this:

attribute: sidechainSAS
recipient: residues
        :1.A    130.75730800628662
        :2.A    76.83372020721436
        :3.A    131.37474918365479
        :4.A    119.24635934829712
        :5.A    24.345003843307495
        :6.A    74.95613121986389
(... etc. ...)

See also the discussion on calculating normalized surface areas (however, the normalizing values are only available for SES, not SAS):
<http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/surfnorm.html>

I hope this helps,
Elaine
-------
Elaine C. Meng, Ph.D.                       
UCSF Computer Graphics Lab (Chimera team) and Babbitt Lab
Department of Pharmaceutical Chemistry
University of California, San Francisco

P.S. Recently some of my emails have been going out twice, and I don't know why.  Apologies, and I hope the problem will not continue!

On Jun 26, 2014, at 11:32 PM, saadraza at bs.qau.edu.pk wrote:

> Respected Admin
> I would like to know if there is any way to modify the calculation of
> Accessible Surface Area (ASA) to calculate only the ASA for the side chain
> or a particular group of atoms in chimera. If not possible could someone
> guide me to some other tool. I will be obliged if someone sort this
> problem for me
> Regards
> Saad Raza




More information about the Chimera-users mailing list