BERIL Knowledge Wiki

Method in Metabolic Models and Pathway-Level Evidence

Flux Balance Analysis

Well corroborated. AI-generated wiki from AI-conducted research reports. Partially reviewed by scientists. See About This Wiki.

Sources cited
5 project reports and 1 cross-project digest, listed in the margin
Relations
supports 10refines 7contradicts 1
Linked from
12 pages5 topics

What this entity is

Canonical name: Flux Balance Analysis.
Known alias: FBA. 1

Flux Balance Analysis is a metabolic-modeling method represented in the ADP1 database through predicted metabolic fluxes and growth phenotypes. 1 No stable external identifier for this method is specified in the source document. 1

Key facts and concordance evidence

FBA flux data covered 15% of the 5,852 genes in the central genome_features table. 1

Among 866 genes with both FBA flux predictions and TnSeq essentiality calls, 639 genes, or 73.8%, had concordant results, while 227 genes were discordant. 1 The discordant genes were identified as candidates for metabolic-model refinement or as possible examples of regulatory effects not represented by FBA. 1

FBA flux classes changed between rich and minimal media for 177 of the 866 genes, or 20%, indicating condition-dependent metabolic rewiring among the genes with flux data. 1

The newer analysis refines this earlier concordance claim by separating prediction tasks: among 478 TnSeq-dispensable genes, FBA class did not predict quantitative growth defects (chi-squared = 0.93, p = 0.63, 2 df; Kruskal-Wallis H = 1.67, p = 0.43). 2 In contrast, FBA showed moderate concordance with complete-knockout lethality: rich media yielded recall = 60.8%, precision = 64.0%, specificity = 79.7%, F1 = 0.624, and Cohen’s kappa = 0.486 across 724 genes; minimal media yielded recall = 65.6%, precision = 69.2%, specificity = 78.9%, F1 = 0.673, and Cohen’s kappa = 0.493 across 833 genes. 2

Thus, the earlier 73.8% concordance and the newer null result are not directly contradictory: they concern different comparisons. The newer analysis supports FBA as a moderate classifier at the lethal-versus-dispensable boundary but contradicts using FBA class to explain growth variation among TnSeq-dispensable genes. 1,2

Condition-specific FBA flux correlations with measured growth were weak and mixed, including glucose ρ = -0.021 (p = 0.677; n = 387), acetate ρ = -0.153 (p = 0.004; n = 352), asparagine ρ = -0.257 (p < 0.001; n = 286), and glucarate ρ = +0.246 (p = 0.005; n = 127); the opposite-than-expected glucarate direction suggests condition-specific model inaccuracies. 2 Aromatic degradation was enriched among FBA-discordant genes (9 of 11; OR = 9.70; FDR-adjusted q = 0.012), with directional enrichment for FBA under-prediction (OR = 12.0; q = 0.004), supporting environmental-assumption gaps as a refinement target. 2

The aromatic-catabolism analysis supports and refines this environmental-gap interpretation: FBA captured 1.76× higher Complex I flux on aromatic substrates than on the comparison condition, with fluxes of 0.55 versus 0.31, but predicted 0% essentiality for Complex I. 3 This indicates that growth-optimizing flux redistribution can represent increased respiratory demand without identifying disruption of a multi-subunit respiratory complex as a bottleneck. 3 Of 51 quinate-specific genes, 30 had no FBA reaction mappings, including cofactor-supply, iron-acquisition, regulatory, and newly identified Complex I-associated functions, showing that the model represents core metabolism more completely than some supporting infrastructure. 3 Independent quinate-specific growth defects in 10/13 Complex I operon subunits further contradict interpreting 0% model essentiality as evidence that the corresponding functions are biologically dispensable. 3

The respiratory-chain wiring study supports and refines this limitation: FBA predicted zero flux through NDH-2 and ACIAD3522 on all standard media because growth optimization preferentially routed NADH through the higher-ATP-yielding Complex I. 4 Consequently, the model can miss alternative respiratory pathways and capacity constraints that force use of suboptimal routes, even when gene-phenotype data indicate condition-specific requirements. 4 This interpretation is based on model behavior and theoretical pathway stoichiometry rather than measured flux distributions. 4

The discoveries synthesis supports the condition-specific interpretation: in ADP1, each of 8 carbon sources required a distinct respiratory configuration; quinate required only Complex I, acetate required Complex I plus cytochrome bo3 and ACIAD3522 among other components, and glucose had no specifically required component. 5 Quinate produced 0.57 NADH per carbon versus 1.50 for glucose, yet Complex I was more essential because simultaneous β-ketoadipate products created a concentrated NADH burst. 5 This refines interpretation of flux magnitude alone: a lower total reductant yield can still impose a sharper respiratory bottleneck. 5

The same synthesis reports that FBA lacked reaction mappings for 30/51 quinate-specific genes (59%) and predicted 0% Complex I essentiality despite predicting 1.76× higher aromatic flux. 5 These results support treating missing respiratory, cofactor-supply, iron-acquisition, and regulatory reactions as model gaps rather than treating the essentiality prediction as evidence against the measured phenotype. 5

Reaction conservation and gapfilling

Across 1,330 unique metabolic reactions, 1,248, or 94%, were shared across all 14 analyzed genomes and classified as core, 62 were variable, and 20 were genome-unique. 1

Gapfilling accounted for 7.7% of reactions on average, and 243 missing functions were cataloged. 1 Of 121,519 growth phenotype predictions across 14 genomes, 105,376, or 87%, required at least one gapfilled reaction, making growth-prediction interpretation tightly dependent on gapfilling quality in this analysis. 1

The annotation-gap study supports the importance of this dependency: across 574 organism-carbon source combinations, baseline FBA achieved 42.5% overall accuracy, with recall of 86.5% (244 of 282 growth-positive conditions correctly predicted) and precision of 42.5% (244 of 574 growth predictions correct), producing 330 false positives. 6 Conditional gapfilling for 38 false-negative cases added 219 reactions: 201 enzymatic, 14 transport, and 12 exchange, averaging 5.8 reactions per case. 6

The same study refines gapfilling interpretation by showing that evidence triangulation assigned candidate genes to 96 of 201 gapfilled enzymatic reaction-organism pairs (47.8%), including 44 high-confidence, 19 medium-confidence, and 33 low-confidence pairs; 105 pairs (52.2%) remained unresolved. 6 BLAST alone resolved 70 pairs (34.8%), whereas the full pipeline resolved 96 (47.8%), so no single evidence stream exceeded 35% resolution. 6 Of 50 EC-less “dark reactions,” only 8 (16%) were resolved, compared with 88 of 151 (58.3%) reactions with known EC numbers, identifying a specific annotation limitation rather than evidence that all gapfills are equally uncertain. 6

FBA validation in the annotation-gap study was inconclusive: 23 gene-protein-reaction (GPR) rules—formal links between genes and reactions—were inserted into SBML models, but knockout simulations produced zero wildtype growth because the models required the gapfilled reactions themselves to grow on the tested minimal media. 6 This refines the earlier use of FBA for essentiality assessment by showing that validation can become circular when the tested reaction is itself required for model growth. 6

Across 7 Fitness Browser organisms and 23 GapMind pathways, 35.4% of pathway-organism pairs were Active Dependencies, 41.0% Latent Capabilities, 14.9% Incomplete but Important, and 8.7% Missing; all Latent Capabilities became fitness-important under condition-specific analyses. 5 This supports using measured condition-specific fitness to test whether apparently complete FBA pathways are latent rather than biologically irrelevant capabilities.

The aromatic analysis supports adding respiratory-chain capacity, PQQ biosynthesis, and iron-homeostasis constraints to ADP1 FBA, because these supporting functions were not fully represented despite their measured association with quinate growth. 3

The FBA results support comparison of metabolic-model predictions with tnseq essentiality measurements and feed into gene-essentiality. 1

The reaction-conservation and gapfilling results support analysis in metabolic-model-gapfilling. 1 The aromatic-discordance result refines that analysis by identifying missing aromatic substrates or mismatched environmental assumptions as testable model gaps. 2 The respiratory-chain study further supports this connection by identifying growth optimization and omitted capacity constraints as specific causes of false-negative pathway-use predictions. 4

The annotation-gap workflow further connects FBA to condition-specific-fitness, pangenome-integration, and multi-omics-integration through fitness, pangenome, pathway, annotation, and homology evidence. 6 The discoveries synthesis supports these connections by showing that pathway completeness and condition-specific fitness can diverge, particularly for latent capabilities. 5 The respiratory-chain analysis also supports multi-omics-integration by combining FBA with gene-phenotype measurements, theoretical stoichiometry, cross-species fitness, and proteomics. 4

The full source-level accounts are Acinetobacter baylyi ADP1 Data Explorer, Triple Essentiality Concordance Analysis for *Acinetobacter baylyi* ADP1, Annotation-Gap Discovery via Phenotype-Fitness-Pangenome-Gapfilling Integration, Aromatic Catabolism Support Network in ADP1, discoveries, and Condition-Specific Respiratory Chain Wiring in ADP1. 1,2,6,3,5,4

Sources cited

  1. Acinetobacter baylyi ADP1 Data Exploreracinetobacter_adp1_explorer · cited 14× · primary source · supports · contradicts · refines
  2. Triple Essentiality Concordance Analysis for *Acinetobacter baylyi* ADP1adp1_triple_essentiality · cited 7× · refines · supports · contradicts
  3. Aromatic Catabolism Support Network in ADP1aromatic_catabolism_network · cited 6× · supports · refines
  4. Condition-Specific Respiratory Chain Wiring in ADP1respiratory_chain_wiring · cited 6× · supports · refines
  5. Discoveries Logdiscoveries · cited 8× · cross-project digest · supports · refines
  6. Annotation-Gap Discovery via Phenotype-Fitness-Pangenome-Gapfilling Integrationannotation_gap_discovery · cited 9× · supports · refines