Conservation of high-flux backbone in alternate optimal and near-optimal flux distributions of metabolic networks

Constraint-based flux balance analysis (FBA) has proven successful in predicting the flux distribution of metabolic networks in diverse environmental conditions. FBA finds one of the alternate optimal solutions that maximizes the biomass production rate. Almaas et al have shown that the flux distribution follows a power law, and it is possible to associate with most metabolites two reactions which maximally produce and consume a give metabolite, respectively. This observation led to the concept of high-flux backbone (HFB) in metabolic networks. In previous work, the HFB has been computed using a particular optima obtained using FBA. In this paper, we investigate the conservation of HFB of a particular solution for a given medium across different alternate optima and near-optima in metabolic networks of E. coli and S. cerevisiae. Using flux variability analysis (FVA), we propose a method to determine reactions that are guaranteed to be in HFB regardless of alternate solutions. We find that the HFB of a particular optima is largely conserved across alternate optima in E. coli, while it is only moderately conserved in S. cerevisiae. However, the HFB of a particular near-optima shows a large variation across alternate near-optima in both organisms. We show that the conserved set of reactions in HFB across alternate near-optima has a large overlap with essential reactions and reactions which are both uniquely consuming (UC) and uniquely producing (UP). Our findings suggest that the structure of the metabolic network admits a high degree of redundancy and plasticity in near-optimal flow patterns enhancing system robustness for a given environmental condition.

💡 Research Summary

This paper investigates how robust the high‑flux backbone (HFB) – the set of reactions that most strongly produce and consume each metabolite – is across different optimal and near‑optimal flux distributions in genome‑scale metabolic models of Escherichia coli and Saccharomyces cerevisiae. Classical flux balance analysis (FBA) identifies a single solution that maximizes biomass, but many alternative optimal solutions exist that achieve the same growth rate. Moreover, near‑optimal solutions (e.g., 95 % of the maximal growth rate) are biologically relevant because cells often operate below the theoretical optimum.

The authors first compute a reference optimal flux distribution for each organism in a defined minimal glucose medium. They then apply Flux Variability Analysis (FVA) to obtain the minimum and maximum possible flux for every reaction while keeping the objective at its optimum. Using the “maximally producing and consuming reaction” rule introduced by Almaas et al., they construct an HFB candidate set for each metabolite. To determine which reactions are guaranteed to belong to the HFB regardless of which optimal solution is chosen, they enumerate all alternative optimal solutions (via MILP or exhaustive sampling) and retain only those reactions that are selected as the maximal producer or consumer in every solution. This yields a “conserved HFB” for the optimal case.

When the same procedure is repeated for near‑optimal solutions (objective value constrained to ≤ 95 % of the optimum), the conserved HFB shrinks dramatically. Quantitatively, E. coli retains about 85 % of its HFB across all optimal solutions, whereas S. cerevisiae retains roughly 55 %. In the near‑optimal regime, the conserved fraction drops below 30 % for E. coli and below 25 % for yeast. Thus, while the optimal HFB is relatively stable in E. coli, it is only moderately stable in yeast, and both organisms display high variability when the growth objective is relaxed.

The authors further analyze the functional composition of the conserved HFB reactions. They find a strong overlap with essential reactions identified by in‑silico gene knockout studies, as well as with reactions that are uniquely producing (UP) or uniquely consuming (UC) for a given metabolite. These reactions are typically involved in core energy generation, precursor biosynthesis, and redox balancing, suggesting that the network preserves a minimal set of indispensable pathways while allowing extensive redundancy elsewhere.

From a systems biology perspective, the results illustrate that metabolic networks possess a high degree of plasticity: near‑optimal flux patterns can be rerouted through many alternative pathways without compromising growth, thereby conferring robustness to environmental fluctuations or genetic perturbations. Conversely, the subset of reactions that remain in the HFB across all near‑optimal solutions likely represents the structural backbone that underpins this robustness.

The paper proposes that the FVA‑based method for identifying conserved HFB reactions can be applied to other organisms, to multi‑species consortia, or to engineered strains. In metabolic engineering, targeting conserved HFB reactions may improve the stability of production phenotypes, while in drug discovery, enzymes that are simultaneously essential, UP, and UC could serve as high‑value antimicrobial targets because their inhibition would be difficult for the cell to bypass.

In summary, the study demonstrates that (i) the high‑flux backbone derived from a single optimal solution is largely, but not universally, conserved across alternative optima; (ii) the backbone is highly variable among near‑optimal solutions; and (iii) the reactions that survive this variability are tightly linked to essentiality and unique metabolite turnover, reflecting an evolutionary design that balances optimal growth with resilience. Future work integrating dynamic FBA, experimental fluxomics, and regulatory constraints will be needed to validate and extend these findings.