[Chimera-users] Calculating SAS and SES with Chimera-MSMS wrap

[Jaime Rodríguez-Guerra]

Dear Chimera dev team,

I am trying to use your MSMS wrapper to compute some SAS/SES values for a
given set of atoms. As suggested in a previous post in this list, I've
taken a look into Measure.measure.buried_area() and I think I've got it,
but just wanted to make sure.

Here is my current procedure:

```
import chimera, Measure, MoleculeSurface
atoms = chimera.selection.currentAtoms()
xyzr_data = Measure.measure.atom_xyzr(atoms)
surfaces = MoleculeSurface.xyzr_surface_geometry(xyzr_data)

for i, a in enumerate(atoms):
    print "Atom: {0}, SAS: {1}, SES: {2}".format(a.name, surfaces[3][i,0],
surfaces[3][i,1])

```

Is this OK? I don't know how to handle all that output that MSMS is giving
back!

Thanks in advance,
Jaime.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://plato.cgl.ucsf.edu/pipermail/chimera-users/attachments/20140407/0aa79867/attachment.html>