FluxSprout: Single-cell metabolic modeling of vascular sprouting in zebrafish

Genome-scale metabolic modeling offers a powerful framework for identifying cell-state-specific metabolic vulnerabilities, yet its application to heterogeneous biological systems remains limited by uncertainty in network annotations, thermodynamic feasibility, transcriptome-to-flux inference, and target specificity. Here, we introduce FluxSprout, a single-cell transcriptomics-informed genome-scale metabolic modeling framework to resolve metabolic heterogeneity during zebrafish sprouting angiogenesis. By integrating single-cell RNA sequencing of 10,090 quality-controlled cells with Zebrafish-GEM, we reconstructed network ensembles for 14 endothelial and hemogenic states. We established four safeguards: stoichiometric auditing, loop auditing, cell-cycle-controlled testing, and ensemble specificity screening. First, we identified four inorganic pyrophosphatase reactions misannotated as oxidative phosphorylation, showing that apparent pathway essentiality arose only with unbounded exchange fluxes. Second, 19.5% of reactions could carry thermodynamically infeasible loop flux, including every candidate of a conventional screen. Third, tip-like and stalk-like cells differed mainly in cell-cycle phase, and elevated glycolysis separated sprouting from established vessels only in gene expression. In contrast, sprouting vessels had higher predicted growth capacity at both constraint strengths, with serine and nucleotide synthesis following in fixed proportion. Fourth, screening 825 genes across four network-extraction settings and 13 non-endothelial lineages left no candidate called by more than one setting. Notably, SLC25A20, that screen's lead candidate, had no predicted effect in any FASTCORE network. Although the data derive from a single pooled clutch and require experimental validation, FluxSprout provides a transferable, uncertainty-aware framework for distinguishing genuine metabolic dependencies from artifacts of network reconstruction, model assumptions, and cellular heterogeneity, and identifies biosynthetic capacity as a testable feature of sprouting vessels.