PPI benchmark evaluation

PPI Benchmark Studio · engine 0.1.0 (framework 1.0.0) · 2026-10-14T17:12:25+0000 · seed 42 · preset standard

Shortcut-dominated — Most of this benchmark's apparent signal is reproduced by network degree, and little survives when that shortcut is removed.

Readiness score 25/100. Heuristic summary of the diagnostics below (100 − 15 per concern − 5 per caution); not a statistical quantity.

Diagnostics

pass

Data integrity

24 rows removed during cleaning (0.2%)

24 duplicate pairs (0 reversed), 0 self-pairs and 0 pairs with conflicting labels were removed; 11,164 canonical pairs remain.

For model building: Reversed duplicates (A–B and B–A) placed on both sides of a random split leak labels directly. Train on the canonical pair list, not the raw file.

fail

Degree bias in supplied negatives

Endpoint degree SMD 1.16 between positives and supplied negatives

Positive pairs join proteins with mean positive-network degree 6.72, supplied negatives 3.34. Degree-matched N2 brings the SMD to 0.01.

For model building: A model can tell the classes apart by how well-studied the proteins are, without learning anything about interaction. Compare against degree-matched negatives before trusting a score.

Constructed negatives are putative (unobserved) pairs, not experimentally established non-interactions.

fail

Degree-only shortcut

Degree alone reaches ROC-AUC 0.962 on the reference labels

A logistic model on four endpoint-degree summaries, with no sequence information, scores 0.962 on N0 and 0.502 on degree-matched N2 negatives.

For model building: Most of the benchmark's label separation is reproducible from network degree. Headline scores on these labels mix interaction signal with annotation density.

Degree is computed from the same incomplete positive graph that supplies the labels; it is a retrospective audit feature, not an input available for a protein with no known interactions.

pass

Length-only shortcut

Sequence length alone: ROC-AUC 0.540

Endpoint sequence-length summaries used as the only features.

For model building: Sequence length does not separate the labels.

fail

Performance on unseen proteins

ROC-AUC falls by 0.109 from seen (R1) to unseen proteins (R3)

AAC+CTD MLP with the training set held fixed: R0 0.701, R1 0.697, R3 0.587.

For model building: Random pair splits overstate performance on new proteins. Report the protein-disjoint R3 split if the model is meant for proteins without known interactions.

R3 test sets are small; see the uncertainty chapter.

fail

Signal after degree matching

Close to chance (0.517) with degree-matched negatives on unseen proteins

Best representation AAC+CTD, MLP trained and tested with N2 negatives at R3: ROC-AUC 0.517 (N0: 0.587). Component bootstrap 95% CI 0.400–0.960; above chance under both corrected tests: no.

For model building: With degree shortcuts removed and proteins unseen, there is little sequence signal for a model to learn from this dataset as constructed.

A non-significant test does not establish equivalence to chance; few connected components limit what can be concluded.

fail

Independent test units

R3 test graphs contain 5–19 connected components

Test pairs that share proteins are not independent. Uncertainty is computed by resampling whole connected components of the test-pair graph.

For model building: Confidence intervals are wide and tests have little power; differences between models of a few hundredths are not interpretable.

A non-significant test does not establish equivalence to chance; few connected components limit what can be concluded.

Recommendations

Negative-set balance

negative_setn_negativespositive_meannegative_meansmdks_statisticlength_smd
N055836.7243.3421.1600.4800.043
N155816.7244.4520.7650.3050.026
N255816.7246.6950.0110.005-0.009
N355816.7245.5730.3900.1550.005
N555816.7246.6830.0150.007-0.001
N2-exact55816.7246.7240.0000.0000.011
N2-nearest55816.7246.7240.0000.0000.011

Splits and leakage gates

regimetrain_pairstest_pairstrain_proteinstest_proteinsprotein_overlaptest_componentspassed
R08931223324972080208097True
R157031426199715611561178True
R2570335701997216416672True
R357034621997420019True

Protein-novelty curve

representationregimemodeltest_pairsroc_aucpr_aucmcc
AAC+CTDR0lr22330.6400.6150.212
AAC+CTDR0rf22330.6840.6550.291
AAC+CTDR0mlp22330.7010.6800.310
AAC+CTDR1lr14260.6600.6530.252
AAC+CTDR1rf14260.7010.6850.321
AAC+CTDR1mlp14260.6970.6840.305
AAC+CTDR2lr35700.6060.5690.151
AAC+CTDR2rf35700.6520.6030.218
AAC+CTDR2mlp35700.6260.5870.191
AAC+CTDR3lr4620.5580.5220.109
AAC+CTDR3rf4620.5870.5230.136
AAC+CTDR3mlp4620.5870.5490.129

Topology-only and length-only controls

featureevaluated_onroc_aucpr_aucmcc
degreeN00.9620.9650.779
degreeN10.8130.7710.486
degreeN20.5020.5050.019
degreeN30.7000.6500.333
degreeN50.5010.5020.029
degreeN2-exact0.4990.5000.026
degreeN2-nearest0.4990.5000.026
lengthN00.5400.5320.064
lengthN10.5200.5170.028
lengthN20.5060.5090.019
lengthN30.5090.5130.020
lengthN50.5090.5120.012
lengthN2-exact0.5160.5150.024
lengthN2-nearest0.5160.5150.024

Consistent-negative revalidation

representationregimenegative_settest_pairsroc_aucpr_aucmcc
AAC+CTDR0N022340.7200.6910.329
AAC+CTDR3N04620.5870.5490.129
AAC+CTDR0N222340.4550.464-0.061
AAC+CTDR3N24460.5170.528-0.054
AAC+CTDR0N522340.4460.464-0.074
AAC+CTDR3N54450.4620.477-0.015

Component-aware uncertainty (R3)

representationnegative_setn_testn_componentsroc_aucboot_lowboot_highjack_lowjack_highboot_p_bhjack_p_bh
AAC+CTDN0462190.5870.5360.8440.4030.7720.0840.728
AAC+CTDN244650.5170.4000.9600.1370.8970.1050.906
AAC+CTDN5445100.4620.4270.5780.3430.5810.4800.728

Methods

The dataset (2,497 proteins, 11,188 pairs) was audited and canonicalized to 11,164 unique undirected pairs (5,581 positive, 5,583 negative). Following the protein-novelty and negative-sampling control framework, a random pair split (R0) and a matched protein-novelty curve (R1 both proteins seen, R2 one unseen, R3 both unseen; shared training set of 5,703 pairs) were built from one partition of identical-sequence groups, with automated leakage gates. Putative negative sets (N0, N1, N2, N3, N5, N2-exact, N2-nearest) were constructed by one-sided substitution matched on positive-network degree and/or sequence length deciles. AAC+CTD descriptors (167-d) with symmetric pair fusion were used with lr, rf, mlp models. Degree-only and length-only classifiers were trained on the reference negatives and evaluated on each alternative set, and an MLP was re-evaluated with each negative definition used consistently for training and testing at R0 and R3. Uncertainty was quantified with a connected-component bootstrap (1000 replicates) and a delete-one-component jackknife, with Benjamini–Hochberg adjustment. Seed 42.