SASA calculator

What is SASA?

Solvent Accessible Surface Area (SASA) measures how much of a protein's surface is exposed to the surrounding solvent. This property reveals which residues are on the protein's exterior versus buried in its hydrophobic core.

SASA is fundamental to understanding protein folding, stability, and function. Exposed hydrophobic residues often indicate binding sites or regions that may aggregate. Changes in SASA between conformational states can quantify domain movements or ligand-induced structural changes.

For a comprehensive analysis of your structure, combine SASA with other structural tools like the Ramachandran Plot for backbone geometry, DSSP for secondary structure assignment, or PDBsum for structural summaries.

How does the Shrake-Rupley algorithm work?

The Shrake-Rupley algorithm calculates SASA by computationally "rolling" a probe sphere (representing a water molecule) over the protein surface. First introduced in 1973, it remains the standard method for SASA calculation.

The rolling ball concept

Each atom is represented as a sphere with its van der Waals radius. The algorithm expands these radii by the probe radius (typically 1.4 Å for water) to create an accessible sphere. Points distributed on this expanded sphere are tested for overlap with neighboring atoms.

ProteinIQ uses 92 test points per atom. For structures with alternate atom locations, it selects one coherent residue conformer using occupancy and deterministic tie-breaking before calculating SASA.

The accessible surface area for each atom equals the fraction of test points not buried by neighbors, multiplied by the sphere's surface area:

$$\text{SASA}_{\text{atom}} = \frac{N_{\text{accessible}}}{N_{\text{total}}} \times 4\pi r^2$$

where $r$ is the van der Waals radius plus the probe radius, and $N$ is the number of test points.

Van der Waals radii

The algorithm uses standard van der Waals radii for each element. Carbon atoms have a radius of 1.7 Å, nitrogen 1.55 Å, oxygen 1.52 Å, and sulfur 1.8 Å. These radii define the physical size of each atom.

Relative accessibility

For residue-level output, the calculator also reports relative accessibility—the percentage of a residue's surface that is exposed compared to its maximum possible exposure. This is calculated as:

$$\text{Relative accessibility} = \frac{\text{SASA}_{\text{residue}}}{\text{SASA}_{\text{max}}} \times 100\%$$

Maximum SASA values use the theoretical ALLOWED Gly-X-Gly reference scale from Tien et al. (2013). Residues with relative accessibility below 20% are typically considered buried, while those above 50% are surface-exposed. The reference scale is tied to the 1.4 Å rolling-probe definition and is not defined reliably for chain-terminal residues, so the calculator leaves relative accessibility blank for larger probes and chain termini while still reporting absolute SASA.

Input requirements

PDB structure

Upload one or more PDB files containing your protein structure. The calculator processes protein ATOM records and excludes hydrogen and HETATM records. You can also fetch structures directly from the RCSB PDB using their 4-character IDs.

For a multi-model PDB, the calculator analyzes the first coordinate model. Submit individual models as separate files when you need model-by-model or ensemble comparisons.

Settings

• Output level: Choose the granularity of results. Structure returns a single total SASA value. Chain breaks down SASA by each polypeptide chain. Residue provides per-residue accessibility, which is most useful for identifying surface-exposed positions.

• Probe radius: The radius of the virtual solvent sphere. The standard value of 1.4 Å represents a water molecule and supports relative accessibility. Larger probes (1.8–2.0 Å) can model bulkier solvents or identify only the most accessible regions, but return absolute SASA only.

Understanding the results

Structure-level output

At the structure level, you get the total SASA in Ų along with counts of chains and residues. Typical globular proteins have SASA values ranging from a few thousand to tens of thousands of Ų, depending on size.

Chain-level output

Each chain is listed separately with its total SASA and residue count. This is useful for comparing surface exposure between subunits or identifying which chains contribute most to the complex's surface.

Residue-level output

The most detailed view shows:

• Residue: The three-letter amino acid code
• Chain: PDB chain identifier
• Position: Sequence position within the chain
• Insertion Code: PDB insertion code when present
• SASA (Ų): Absolute surface area for that residue
• Rel. Accessibility (%): Percentage of maximum possible exposure
• RSA Status: Whether RSA was calculated, omitted for a terminal residue, or unavailable for a residue outside the standard reference scale

Residues with high relative accessibility are good candidates for surface mutations or chemical modifications. Those with unexpectedly low accessibility despite being charged (Lys, Arg, Glu, Asp) may indicate buried salt bridges.

Use cases

SASA analysis is valuable when designing mutations—you can confirm a target residue is surface-exposed before introducing modifications. It also helps identify potential binding interfaces, which often show intermediate accessibility values.

Comparing matched structures can reveal changes in protein surface exposure, but this calculator does not directly calculate buried interface area. A buried-surface calculation requires consistent SASA calculations for the isolated partners and their complex; ligand atoms stored as HETATM are not included here. Monitoring SASA changes across molecular dynamics trajectories reveals conformational dynamics.

Limitations

The Shrake-Rupley algorithm treats each submitted model as a static structure. For flexibility, submit individual conformations separately or use molecular-dynamics analysis. The calculation excludes hydrogen atoms and all HETATM records.

Relative accessibility values assume the Tien et al. 1.4 Å Gly-X-Gly reference state, which may not perfectly represent every local environment. Non-terminal values exceeding 100% remain visible rather than being clipped. RSA is unavailable with larger probes and for chain-terminal or unsupported non-standard residues.