pySCA Now Available on Vecura: Uncover Protein Coevolution and Functional Sectors

What is pySCA?

pySCA (Statistical Coupling Analysis) is a Python 3 toolkit that extracts evolutionary constraints from protein family multiple sequence alignments (MSAs). Starting from aligned sequences, the tool filters and reweights sequences to account for phylogenetic redundancy, computes per-position conservation scores, and calculates a coevolution matrix that reveals how pairs of residues co-vary across the protein family. The method then decomposes this matrix using Independent Component Analysis to identify "sectors" — groups of physically connected, co-evolving residues that typically encode specific functions like catalysis, allostery, or binding specificity.

It helps users identify which residues in a protein family evolved together as functional units, revealing the evolutionary architecture underlying protein function. It is especially useful for mapping allosteric communication pathways, identifying specificity-determining positions, and guiding protein engineering experiments by pinpointing evolutionarily coupled residue networks.

What can users do with pySCA on Vecura?

With pySCA on Vecura, users can:

• Run the complete SCA pipeline from raw MSA to identified sectors in a single workflow, including sequence filtering, conservation analysis, and coevolution matrix computation

• Identify protein sectors — groups of co-evolving residues that form functional or allosteric units — with automatic selection of significant eigenmodes

• Map sector positions onto 3D protein structures using PDB references to visualize co-evolving networks in structural context

• Export comprehensive results including conservation scores, coevolution matrices, eigenvalue spectra, and sector assignments for downstream analysis

What the output means

The output provides conservation scores per position (relative entropy Di), the full SCA coevolution matrix showing residue-residue coupling strengths, eigenvalue spectra for determining significant modes, and identified sectors with their constituent positions (optionally mapped to PDB residue labels). Each sector represents a group of residues that co-evolved as a unit, typically underlying a specific biological function.

This output should be used to support scientific decision making. It does not replace experimental validation. Sector predictions identify candidate functional networks that require experimental confirmation through mutagenesis, functional assays, or structural studies.

Why this matters

Proteins are not uniform entities — they contain functionally specialized regions that evolved together to perform specific tasks. Traditional conservation analysis identifies important individual residues, but misses the collective evolutionary units that encode complex functions like allostery or multi-site catalysis. Statistical Coupling Analysis, pioneered by the Ranganathan lab and formalized in landmark studies (Halabi et al., Cell 2009; Rivoire et al., PLoS Comput Biol 2016), provides a data-driven framework to decompose proteins into these evolutionary sectors.

By revealing which residues co-evolved as functional modules, pySCA enables researchers to understand the evolutionary logic of protein families, predict allosteric pathways, identify residues critical for specific functions, and design targeted mutations for protein engineering. This approach has transformed our understanding of how evolution builds complex protein functions through coordinated changes in distributed residue networks.

• Developed by: Ranganathan Lab (University of Texas Southwestern Medical Center)

• Source: Official GitHub Repository | Documentation

• Reference: Halabi et al., Cell 2009 (https://doi.org/10.1016/j.cell.2009.07.038); Rivoire et al., PLoS Comput Biol 2016 (https://doi.org/10.1371/journal.pcbi.1004817)