Regenerate with uv run python sigma5_structure.py in this directory. Structures are
downloaded from RCSB and sequences from UniProt at run time; expected sequence lengths
are asserted, and no result below is hardcoded.
| chain | protein | modelled residues | UniProt span | numbering offset | fraction matching after offset |
|---|---|---|---|---|---|
| A | zeta5 (AP5Z1, mouse) | 671 | 1-787 | -1 | 0.9985 |
| B | beta5 (AP5B1, human) | 551 | 7-631 | 0 | 1.0 |
| C | sigma5 (AP5S1, human) -- the target | 182 | 1-199 | 0 | 1.0 |
| D | SPG11 WD40-hairpin (spatacsin, human) | 402 | 22-521 | 0 | 1.0 |
| E | mu5 (AP5M1, mouse) | 468 | 2-490 | 0 | 1.0 |
| S | sigma2 (AP2S1, mouse) in 2JKR | 142 | 1-142 | 0 | 1.0 |
Every chain indexes its UniProt sequence directly except zeta, whose deposited
numbering is shifted by one (an extra N-terminal residue in the construct); the offset
is found by scanning rather than assumed, and all zeta residue numbers below are
reported in UniProt coordinates.
Sequence lengths verified: {'Q9NUS5': 200, 'Q3U829': 807, 'Q2VPB7': 878, 'Q96JI7': 2443, 'Q8BJ63': 490, 'P62743': 142}
Free solvent-accessible surface of the isolated sigma-5 chain: 11230.0 A^2.
| partner subunit | sigma-5 residues in contact | buried sigma-5 surface (A^2) | closest approach (A) |
|---|---|---|---|
| zeta5 (AP5Z1, mouse) | 60 | 2621.8 | 2.16 |
| beta5 (AP5B1, human) | 0 | 0.0 | - |
| SPG11 WD40-hairpin (spatacsin, human) | 20 | 755.4 | 1.52 |
| mu5 (AP5M1, mouse) | 0 | 0.0 | - |
The zeta N-terminal segment tested is UniProt residues 1-40. It contacts 9 sigma-5 residues (L24, L120, L152, L153, A154, S156, T157, L160, L161); the zeta residues involved are [1, 2, 3, 5, 6, 9, 10].
The site is defined from 2JKR: the sigma2 residues within 4.5 A of the bound CD4 dileucine peptide (modelled sequence MSQIKRLLS). That gives 19 residues: N9, R10, A11, R15, Y62, A63, G64, L65, F67, H85, V88, E89, N92, N97, V98, C99, E100, L101, L103.
Superposition method: Bio.PDB.cealign.CEAligner (Combinatorial Extension, CA atoms, sequence-independent), reference 2JKR chain S (sigma2).
| superposed chain | RMSD (A) | reference CA positions covered within 4 A (of 142) |
|---|---|---|
| sigma-5 (8YAB chain C) | 3.66 | 135 |
| positive control: second sigma2 copy (2JKR chain I) | 0.03 | 142 |
| negative control: beta5 solenoid (8YAB chain B) | 6.82 | 113 |
Read the RMSD column, not the coverage column: beta5 is a 551-residue solenoid, so its
CA atoms are dense enough to fall within 4 A of most reference positions whatever the
fold. RMSD separates the three cases cleanly, and sigma-5 sits between the identical
control and the unrelated fold, which is what a divergent homologue should do.
| sigma2 pocket residue | nearest sigma-5 CA | CA-CA distance (A) | structurally equivalent | same residue | that sigma-5 position contacts the zeta N-terminus |
|---|---|---|---|---|---|
| N9 | T9 | 1.3 | True | False | False |
| R10 | L10 | 2.71 | True | False | False |
| A11 | C21 | 2.0 | True | False | False |
| R15 | Y25 | 1.81 | True | False | False |
| Y62 | V117 | 1.62 | True | False | False |
| A63 | L118 | 1.86 | True | False | False |
| G64 | S119 | 1.79 | True | False | False |
| L65 | L120 | 1.2 | True | True | True |
| F67 | F122 | 0.41 | True | True | False |
| H85 | R140 | 0.88 | True | False | False |
| V88 | T143 | 1.64 | True | False | False |
| E89 | R144 | 1.9 | True | False | False |
| N92 | D148 | 1.25 | True | False | False |
| N97 | R151 | 3.12 | True | False | False |
| V98 | R151 | 1.92 | True | False | False |
| C99 | L152 | 2.18 | True | False | True |
| E100 | S156 | 4.07 | False | False | True |
| L101 | S156 | 2.18 | True | False | True |
| L103 | L159 | 1.99 | True | True | False |
Summary: of 19 sigma2 residues that contact the dileucine peptide, 18 have a structurally equivalent position in sigma-5 (CA within 4 A after superposition), 3 of those carry the identical residue, and 3 of them are in contact with the zeta N-terminus in the assembled AP-5 core.
The dileucine site is defined on P62743 (mouse Ap2s1, the sigma2 chain of 2JKR), while the functional
weight of these positions comes from work on the human protein, P53680 (AP2S1).
Same length: True. Identical over the full chain: True. All peptide-contacting positions identical between the two: True.
| paralog | length | % identity over aligned columns | alignment score | shuffled mean +- sd | z |
|---|---|---|---|---|---|
| AP2S1 (sigma2) | 142 | 20.7 | -55.0 | -81.1 +- 13.2 | 1.98 |
| AP1S1 (sigma1A) | 158 | 27.3 | -32.0 | -68.9 +- 11.9 | 3.09 |
| AP1S2 (sigma1B) | 157 | 24.1 | -51.0 | -71.1 +- 15.2 | 1.32 |
| AP1S3 (sigma1C) | 154 | 21.8 | -40.0 | -72.8 +- 14.9 | 2.21 |
| AP3S1 (sigma3A) | 193 | 23.0 | -49.0 | -69.6 +- 15.0 | 1.37 |
| AP3S2 (sigma3B) | 193 | 24.3 | -60.0 | -67.4 +- 14.2 | 0.52 |
| AP4S1 (sigma4) | 144 | 15.3 | -63.0 | -78.3 +- 11.7 | 1.3 |
Every pairwise alignment of sigma-5 against a human sigma paralog scores negative and
sits within a few standard deviations of shuffled sequence, so the correspondence used
in Q2 could not have been obtained from pairwise sequence alignment. That is consistent
with Hirst et al. 2011 (PMID:22022230), who needed profile methods (HHpred, probability
96.8% against sigma2) rather than BLAST to recognise this protein as a sigma subunit --
and it is why the mapping above is done on structures.
Sigma-5 is the zeta-adaptin partner: it buries 2621.8 A^2 against the zeta trunk over 60 residues, and touches neither beta5 nor mu5. Structurally it occupies the position that sigma1/sigma2/sigma3 occupy in AP-1/AP-2/AP-3.
Sigma-5 is also part of the SPG11 binding surface, burying 755.4 A^2 against the SPG11 WD40-hairpin over 20 residues. This is a role the other AP sigma subunits do not have, and it is the structural counterpart of the pull-down result that SPG11 associates with the zeta/sigma-5 subcomplex and is needed to assemble the AP-5 heterotetramer.
The fold that carries the acidic dileucine site is retained (18 of 19 peptide-contacting positions have a structural equivalent) but the chemistry is not: only 3 of those positions carry the same residue. The hydrophobic pocket linings and the basic patch are the positions that change, and in the assembled AP-5 core 3 of the equivalent positions are covered by the zeta N-terminus. Retention of the sigma fold is therefore not evidence that sigma-5 reads dileucine sorting signals, and the structure gives a positive reason to think it does not.
This is a structural argument, not a binding assay. No one has tested AP-5 against a dileucine motif, so the claim defended here is the negative one: nothing in the structure supports annotating cargo-signal recognition to sigma-5.