Article ID: 467 - Last Modified: March 15, 2013
I want to calculate the surface area of a protein binding site, defined as any residue 6 Å from the ligand, and also the volume of the binding site. How do I do this?
For the surface area, the best option is to use the Binding Surface Area Analysis panel, which you can open with Scripts → Project Table → Binding Surface Area Analysis.
Another option is to use use Workspace → Surface → Molecular Surface (Suite 2010 on) or Display → Surface → Molecular Surface (earlier releases), with the following ASL:
not (mol.entry 2) and fillres within 6 (mol.entry 2)
assuming your ligand is molecule #2 in the complex structure. Make sure that "Surface Context" is set to "Molecule"; also if your ligand is covalently bound or has zero-order bonds to the receptor (and thus not considered a separate molecule) you might want to break the bonds first using Delete → Bonds or Delete → Zero-Order Bonds from the Delete toolbar button. The area of the surface (Connolly) is reported in the Manage Surfaces panel in Å2.
The Script Center has scripts for calculating the solvent-accessible surface area of a ligand in a receptor (
atomic_sasa.py), and the change in surface area on ligand binding (
Probably the best option for the volume is to use Applications → SiteMap. You can pick the ligand to specify which site to consider and how far from the ligand to extend the buffer region. The SiteMap job reports the site volume in the Project Table and also reports the surface area for the surface map in the "Manage Surfaces" panel.
Keywords: surface area, volume , measurement, Manage surfaces
#251: How do I calculate the binding site volume and report the result in Å3 in Maestro?
Type the words or phrases on which you would like to search, or click here to view a list of all
Knowledge Base articles