Beyond Action Entropy: Quotient-Space Exploration for Genome-Scale Metabolic Model Repair
Organizations: Shanghai Jiao Tong University · Shanghai AI Laboratory
Abstract
Repairing scientific models from functional observations differs fundamentally from supervised prediction: feedback may certify a solution without revealing which structural correction is responsible. We study this setting for genome-scale metabolic model (GEM) repair, where multiple reaction edits can explain the same phenotypes and many apparently distinct edits correspond to the same biological mechanism. This many-to-one structure creates a hidden failure mode for conventional exploration: diversity in the output space need not translate into diversity of scientific hypotheses. We introduce QuotientPO, which collapses equivalent repairs into canonical mechanisms and optimizes exploration directly over the resulting quotient space. To make quotient exploration informative under finite rollouts, we derive a kernelized Rényi estimator that resolves graded crowding among distinct repair cores beyond coarse exact-match counts. On 2,212 held-out GEMs, QuotientPO improves Success@32 from 17.93% to 20.10% (+12.1% relative) while consistently increasing distinct successful-core discovery under the same sampling budget. These results establish quotient-space exploration as a principled approach to mechanism-level discovery under verifier-induced equivalence.
Figures & tables
| Method | Greedy | Success@1 | Success@8 | Success@32 | Improved (mean@32) | Success Cores@32 |
| MILP (Oracle) | 15.01 | - | - | - | 25.63 | 0.15 |
| SFT-Hard | 1.72 | 0.51 | 1.96 | 4.49 | 1.64 | 0.09 |
| SFT-Soft | 2.65 | 2.40 | 5.27 | 7.52 | 9.92 | 0.16 |
| Vanilla GRPO | 7.97 | 7.43 | 13.79 | 17.93 | 30.21 | 0.50 |
| MaxEnt | 7.29 | 5.97 | 12.69 | 17.12 | 29.35 | 0.48 |
| GFlowNet | 3.20 | 3.23 | 8.35 | 13.10 | 14.91 | 0.34 |
| Training (200→864) | Successful Cores | |||
| Method | Exact ties (%) | Resolved (%) | Collision | Mean |
| Vanilla | 8.19 8.59 | 60.55 84.86 | 0.5769 | 0.2408 |
| Shannon | 8.60 9.72 | 59.23 80.20 | 0.5687 | 0.2431 |
| Exact R2 | 7.34 10.23 | 59.71 78.20 | 0.5285 | 0.2861 |
| Kernel R2 | 8.57 7.99 | 63.81 72.78 | 0.5176 | 0.3008 |
| Param. | Setting | S@32 |
| – | Default | 20.10 |
| 0.125 | 19.41 | |
| 0.5 | 18.57 | |
| 0.01 | 19.05 | |
| 0.04 | 20.00 |
Appendix figures & tables19 assets
Supplementary material from the paper’s appendix.
Appendix
| Assay category | BacDive families |
| Carbon source | glucose, maltose, sucrose, glycerol, acetate, lactate, pyruvate, succinate, citrate, glutamate |
| Nitrogen source | ammonium, urea, nitrate, nitrite |
| Sulfur source | sulfate |
| Phosphorus source | phosphate |
| Anchor | Best | Best step | Last | Last ADD |
| (no anchor) | 7.91 | 15 | 5.80 | 7367.57 |
| 8.46 | 745 | 7.64 | 23.05 | |
| 864 | 10.98 | |||
| 8.00 | 715 | 7.27 | 6.05 |
| Configuration | Hit@1 | dev-select loss |
| Pretrained + ESM2 | 0.207 | 3.248 |
| Random + ESM2 | 0.114 | 4.524 |
| Pretrained, no ESM2 | 0.143 | 3.922 |
| Parameter | Setting |
| Actor initialization | Soft-SFT |
| Total / trainable parameters | 92.85M / 59.46M |
| Training / validation GEMs | 9,196 / 1,087 |
| RL seeds | 42, 43, 44 |
| Training devices | 8 GPU ranks |
| GEMs per update | 32 |
| Method | Update (s) | FBA drain (s) | Overhead |
| Vanilla GRPO | 21.323 | 2.229 | – |
| QuotientPO | 23.218 | 4.186 | +8.9% |
| Template | Method | Success@1 | Success@8 | Success@32 | Greedy | Improved | Success Cores@32 |
| Default ( ) | |||||||
| Vanilla GRPO | 7.37 0.58 | 12.72 0.79 | 16.79 1.08 | 7.84 0.36 | 28.83 1.08 | 0.516 0.003 | |
| MaxEnt | 5.69 0.96 | 11.81 1.58 | 15.93 1.28 | 7.32 1.52 | 28.04 2.33 | 0.462 0.088 | |
| GFlowNet | 1.77 0.22 | 6.07 0.54 | 10.38 0.33 | 2.54 0.30 | 11.48 0.70 | 0.265 0.017 | |
| QuotientPO | 7.89 0.50 | 14.25 0.22 | 18.75 0.82 | 7.80 0.33 | 30.01 1.19 | 0.582 0.006 | |
| Gram-positive ( ) | |||||||
| Paired GEMs (42/43/44) | Vanilla GRPO | QuotientPO | |
| 2 | 295/269/323 | ||
| 4 | 247/228/258 | ||
| 8 | 209/186/206 | ||
| 16 | 147/123/130 |
| Exact-credit ties (%) | Kernel-resolved ties (%) | |||||||
| Method | 0 | 200 | 400 | 864 | 0 | 200 | 400 | 864 |
| Vanilla GRPO | 5.67 | 8.19 | 8.31 | 8.59 | 30.61 | 60.55 | 57.50 | 84.86 |
| Discrete Shannon | 5.67 | 8.60 | 9.99 | 9.72 | 30.61 | 59.23 | 59.77 | 80.20 |
| Exact Rényi-2 | 5.67 | 7.34 | 9.06 | 10.23 | 30.61 | 59.71 | 63.71 | 78.20 |
| Kernel Rényi-2 | 5.67 | 8.57 | 8.60 | 7.99 | 30.61 | 63.81 | 57.23 | 72.78 |
| Method | Success@32 | Success Cores@32 |
| Vanilla GRPO | 17.93 1.54 | 0.50 0.01 |
| Discrete Shannon | 18.61 1.09 | 0.52 0.01 |
| Exact Rényi-2 | 18.04 1.08 | 0.54 0.03 |
| Kernel Rényi-2 | 20.10 0.52 | 0.57 0.02 |
| Method | Success@32 |
| Vanilla | 15.75 |
| MaxEnt | 14.29 |
| GFlowNet | 11.57 |
| QuotientPO | 17.21 |
| Repair alternatives & encoded chemistry | Functional verification & policy diagnosis |
| Case 1 (C1): Mycobacterium terrae Target: Urea utilization Focal pipeline ( ): : vs : Yield ( ): : : | |
| Shared deletions (3): -HEX1 , -MALTabc , -NARK Pathway 1: Uptake blockade -UREAt Pathway 2: Carboxylation blockade (focal) -UREASE Note: UREASE encodes urea carboxylase in the GEM. Both canonical cores contain exactly four deletions. | Phenotypic verification Both cores reduce urea-supported growth from to zero while preserving nitrate-supported growth of and all remaining targets. Reintroducing the focal deletion restores the failed urea-growth phenotype, confirming its functional necessity. Search behavior diagnosis Vanilla produces successful raw proposals containing this intracellular branch ( ), but completely eliminates it during canonical reduction ( ). In contrast, QuotientPO successfully retains both functional intervention sites ( ). |
| Case 2 (C2): Streptomyces inusitatus Target: Alternative proline supply Focal pipeline ( ): : vs : Yield ( ): : : | |
| Shared deletion (1): -UREAt Pathway 1: Direct conversion +ORNCD Pathway 2: Semialdehyde/P5C route (focal) +G5SADs Note: GSA denotes glutamate-5-semialdehyde; P5C denotes -pyrroline-5-carboxylate. Both canonical cores contain two edits. | Phenotypic verification Both two-edit cores restore proline supply and recover citrate-supported growth from zero to while preserving the inactive urea target. Flux witnesses independently confirm active throughput along both pathways. Search behavior diagnosis Vanilla already contains the focal G5SADs route in successful proposals but loses it after canonicalization ( ); QuotientPO retains it ( ). Thus, this is a retention rather than discovery gap. ORNCD and ORNCD 2 are treated as one functional pathway family. |
| Repair alternatives & encoded chemistry | Functional verification & policy diagnosis |
| Case 3 (C3): Lactobacillus farraginis Target: Alternative phosphate donors for UTP synthesis Focal pipeline ( ): : vs : Yield ( ): : : | |
| Shared edits (3): +MOHMT , +NTRIR2y , -NTRIR3pp Pathway 1: ATP donor +NDPK2 Pathway 2: PEP donor (focal) +PYK2 Both complete cores contain four edits. | Phenotypic verification The alternatives provide UDP-to-UTP phosphorylation using different phosphate donors. Both QuotientPO cores satisfy all six targets, with glucose growth and , respectively, and nitrate growth . Flux witnesses carry approximately through NDPK2 or through PYK2 . Search behavior diagnosis Vanilla proposes PYK2 four times, but these proposals lack the coordinated repair required for success. Adding NTRIR2y restores the nitrate-positive target. |
| Case 4 (C4): Streptomyces ferralitis Target: Alternative nitrogen donors for asparagine Focal pipeline ( ): : vs : Yield ( ): : : | |
| Shared edits (3): +ARGDC , -FFSD , -SUCRtex Pathway 1: Glutamine donor +ASNS1 Pathway 2: Ammonium donor (focal) +ASNS2 | Phenotypic verification The two reactions complete asparagine synthesis using different nitrogen donors. Both QuotientPO cores satisfy the glucose-positive and sucrose-negative targets, with glucose growth and . Search behavior diagnosis Vanilla frequently proposes the focal ASNS2 branch, but the corresponding repair sets remain unsuccessful. One such proposal also introduces extracellular sucrose hydrolysis through SUCRe ; removing this additional bypass restores the sucrose-negative target. Thus, successful use of the focal biosynthetic branch depends on its surrounding repair combination. |
| Repair alternatives & encoded chemistry | Functional verification & policy diagnosis |
| Case 5 (C5): Streptomyces olivochromogenes Target: Alternative control of citrate-supported growth Focal pipeline ( ): : vs : Yield ( ): : : | |
| Shared edit (1): -UREAt Alternative 1: Entry restriction -CITt14 , -FE3t CITt14 directly supports citrate uptake, while FE3t provides a reversible iron route coupled to the remaining FEDCabc citrate-entry pathway. Alternative 2: Intracellular node (focal) -SUCOAS | Phenotypic verification Both cores eliminate the initial urea and citrate false positives while preserving glucose- and sucrose-supported growth. The entry-restriction core retains glucose/sucrose growth at , whereas the SUCOAS core retains . Restoring either CITt14 or FE3t breaks only the citrate-negative target in the first core. Restoring SUCOAS does the same in the focal core. Restoring UREAt affects only the urea-negative target. Search behavior diagnosis Vanilla never proposes -SUCOAS in this rollout panel. QuotientPO discovers it twice. Both samples are successful and canonicalize to the same focal core. Thus, this comparison exposes an intracellular repair alternative absent from Vanilla raw proposals, rather than two distinct mechanisms. |
| Case 6 (C6): Streptomyces durbertensis Target: Alternative quinone-precursor supply Focal pipeline ( ): : vs : Yield ( ): : : | |
| Shared edits (5): +HSERTA , +HSST , +THZPSN2_SC , -NO3abc , -NO3tex Pathway 1: Aromatic branch +PHPYROX Pathway 2: Three-ADD branch (focal) +SHCHCS , +SUCBZS , +SUCBZL_1 Note: SSA denotes succinic semialdehyde and OSB denotes -succinylbenzoate. | Phenotypic verification Both branches converge on OSB-CoA and downstream menaquinone-precursor supply. The resulting cores satisfy all four targets and yield glucose growth of approximately . Minimum-total-flux witnesses use the corresponding precursor branch at approximately , matching the quinone demand induced by biomass production. Search behavior diagnosis The complete three-ADD module does not occur in Vanilla raw proposals, whereas QuotientPO proposes it five times and retains one successful canonical core. Its five total distinct cores therefore correspond to multiple edit sets, rather than five distinct quinone-supply mechanisms. |