Script: active_site_check.py (run with uv run active_site_check.py); raw output in
active_site_check.out. Sequences and feature positions are fetched live from UniProt.
Method: Biopython global pairwise alignment, BLOSUM62, gap open -10, extend -0.5.
Pairwise alignment at 25-35% identity is approximate; residue correspondences below are
a screening result, not a structure-based assignment.
| Reference | Reference position | CG5611 residue |
|---|---|---|
| rat ECHS1 P14604 (enoyl-CoA hydratase) | G141 (oxyanion hole, BINDING) | G148 |
| rat ECHS1 P14604 | E144 (catalytic Glu) | E151 |
| rat ECHS1 P14604 | E164 (SITE, important for catalysis) | N171 |
| CamK Q93TU6 | E124 (ACT_SITE nucleophile, by similarity) | E151 |
| CamK Q93TU6 | H145 (binding) | F169 |
| CamK Q93TU6 | D154 (binding) | D180 |
Identity over aligned columns: ECHS1 34.9%, CamK 24.7%, BamA 30.6%, Srlp 54.9%.
CG5611 retains the oxyanion-hole glycine and the first catalytic glutamate of the
crotonase fold, but in this alignment the second catalytic glutamate of canonical
enoyl-CoA hydratase (ECHS1 E164) corresponds to asparagine. UniProt also reports no
PROSITE enoyl-CoA hydratase/isomerase signature for Q9VB17. Together with the chemically
heterogeneous SF1 membership, this does not support assigning enoyl-CoA hydratase
activity. The exact reaction remains unresolved; the closest family relative, Srlp, has
only localisation (mitochondrial matrix, IDA) and a hypoxia-response phenotype in GO.