S. pombe vms1 has no experimental molecular-function annotation (GOA carries
GO:0003674 molecular_function with evidence ND). Its only biological-process
annotation is an IBA to GO:0036503 ERAD pathway. The Vms1/ANKZF1 family's
best-characterised activity is instead cleavage of polypeptidyl-tRNA on stalled
60S ribosome-nascent-chain complexes, which in human ANKZF1 depends on a single
catalytic glutamine in the VLRF1 (Vms1-like release factor 1) domain.
Does S. pombe vms1 retain that catalytic glutamine, i.e. is it plausibly a
catalytically competent family member rather than a degenerate pseudoenzyme?
vlrf1_catalytic_conservation.py performs global BLOSUM62 pairwise alignments
(Biopython PairwiseAligner, gap open -11 / extend -1) of full-length UniProt
sequences and projects the human ANKZF1 catalytic position onto each ortholog:
| Accession | Protein | UniProt ACT_SITE |
|---|---|---|
| Q9H8Y5 | human ANKZF1 | 246 (Q246L abolishes polypeptidyl-tRNA cleavage) |
| Q04311 | S. cerevisiae Vms1 | 295 |
| O74977 | S. pombe vms1 | 249 (rule-inferred, PROSITE PRU01389) |
Human ANKZF1 is used as the reference because it is the only family member whose
catalytic residue has direct mutagenesis evidence in UniProt.
Reproduce with:
cd genes/SCHPO/vms1/vms1-bioinformatics
uv run python vlrf1_catalytic_conservation.py
FASTA inputs were downloaded from https://rest.uniprot.org/uniprotkb/<ACC>.fasta.
Reference: Q9H8Y5 (human ANKZF1), 726 aa
annotated catalytic site 246: AKRGTA[Q]GLRDAR
Q04311 (S. cerevisiae Vms1), 632 aa [alignment score 242.0]
UniProt ACT_SITE : 295 RKQGGS[Q]SAMDNA
aligned to Q9H8Y5:246 : 295 RKQGGS[Q]SAMDNA
alignment agrees with ACT_SITE annotation: yes
residue at aligned position: Q
O74977 (S. pombe vms1), 600 aa [alignment score 313.0]
UniProt ACT_SITE : 249 RKQGGS[Q]GAADNT
aligned to Q9H8Y5:246 : 249 RKQGGS[Q]GAADNT
alignment agrees with ACT_SITE annotation: yes
residue at aligned position: Q
ACT_SITE call for bothRKQGGSQ motif that isACT_SITE in both queries is the