mass-action-kinetics-formulation
Summary
Computes Reaction Propensity Scores (RPS) from intracellular metabolomics abundance data using mass action law formulation, quantifying the expected relative metabolic flux changes across cell lines based on substrate availability alone. This provides a metabolic-regulation baseline independent of gene expression effects.
When to use
Apply this skill when you have intracellular metabolomics concentration measurements across multiple cell lines or conditions, a stoichiometric metabolic network model with reaction-metabolite associations, and you need to predict how differences in substrate availability (independent of enzyme expression) translate into differences in metabolic flux propensity. Use it specifically to isolate the metabolic regulation layer (substrate-driven) from the transcriptional regulation layer (enzyme abundance-driven) for concordance analysis.
When NOT to use
- When extracellular flux data or direct metabolic flux measurements (e.g., from labeled substrate tracing) are available and preferred; RPS is a proxy, not direct flux measurement.
- When substrate concentrations are not measured intracellularly or are only available in extracellular medium; model assumes intracellular substrate availability drives flux, not external concentration gradients.
- When the metabolic network lacks stoichiometric coefficients or reaction-metabolite associations; the mass action law formulation requires complete stoichiometric data.
Inputs
- Intracellular metabolomics abundance dataset (metabolite concentrations per cell line, replicated measurements)
- Stoichiometric metabolic network model with reaction-metabolite associations and stoichiometric coefficients (e.g., ENGRO2 in SBML or similar format)
- Conversion mapping between metabolomics metabolite identifiers and model metabolite identifiers
Outputs
- Normalized RPS dataset (table with reactions as rows, cell lines as columns, normalized RPS values ∈ [0,1] as entries)
- Pre-aggregation RPS scores (per-reaction, per-cell-line, per-replicate, before normalization)
- Metadata log of omitted reactions (those with unmeasured substrates, per cell line)
How to apply
Load intracellular metabolomics abundance data (metabolite concentrations per cell line) and the stoichiometric metabolic network model (e.g., ENGRO2) with reaction-metabolite associations and stoichiometric coefficients. For each reaction and each cell line, retrieve all substrate metabolites from the model stoichiometry and their measured intracellular concentrations; omit any reaction missing one or more substrate measurements. Compute the RPS for each eligible reaction in each cell line as the product of substrate concentrations each raised to their stoichiometric coefficient power: RPS_r^c = ∏(q=1 to N) [X_q]^(s_r,q), where [X_q] is substrate q concentration and s_r,q is its stoichiometric coefficient. Aggregate RPS values across biological replicates within each cell line using median or mean as appropriate to produce a cell-line-level RPS score. Normalize RPS scores within each reaction across all cell lines by dividing by the maximum RPS value observed for that reaction, yielding dimensionless propensity scores bounded in [0,1] suitable for cross-reaction and cross-cell-line comparison.
Related tools
- ENGRO2 (Stoichiometric metabolic network model providing reaction-metabolite associations and stoichiometric coefficients for RPS computation)
- createMetabolicDataset.py (Preprocessing and preparation of metabolomics data for RPS calculation; quality filtering and mean computation across replicates) — https://github.com/qLSLab/integrate
- concordanceAnalysis.py (Post-processing integration and comparison of RPS output against RAS (transcription-derived scores) to identify metabolic vs. transcriptional regulation) — https://github.com/qLSLab/integrate
- constraint-based stoichiometric metabolic models (Formal scaffold providing stoichiometry, reaction definitions, and metabolite associations necessary for mass action law application)
Evaluation signals
- All reactions in the output RPS table have at least one substrate measurement in at least one cell line; check that no reactions are missing that should have valid data.
- Normalized RPS scores fall within [0, 1] for all reactions and cell lines; the maximum RPS per reaction equals 1.0 (verification of normalization).
- RPS values correlate positively with intracellular substrate concentrations (mass action relationship: higher substrate → higher RPS); audit a sample of 5–10 reactions visually.
- Missing substrates cause consistent reaction omission across the pipeline; verify that reactions with unmeasured substrates are reported in the omission log and absent from the final RPS table.
- RPS concordance with RAS (Reaction Activity Scores from transcriptomics) is significantly less than perfect (Cohen κ and Pearson r < 0.8), confirming that substrate availability captures a distinct regulatory layer from gene expression.
Limitations
- If any substrate is unmeasured, the entire reaction is omitted from the RPS dataset; this may remove important metabolic reactions and limits coverage to well-measured metabolites.
- RPS is a propensity score based on mass action law, not a direct measurement of metabolic flux; it does not account for enzyme kinetics (Km, Vmax), allosteric regulation, feedback inhibition, or post-translational modifications.
- Substrate concentrations are assumed to be at steady state within the 48-hour measurement window; dynamic changes or compartmentalization effects are not captured.
- RPS assumes all reactions are mass-action limited and does not distinguish reactions controlled by regulatory proteins, cofactor availability, or other non-concentration-dependent mechanisms.
- Cross-cell-line normalization by maximum RPS may obscure absolute differences in metabolic propensity; normalized scores are relative and suitable only for concordance analysis, not absolute flux prediction.
Evidence
- [other] For each reaction r and cell line c, the Reaction Propensity Score (RPS) is computed as the product of substrate concentrations each raised to their stoichiometric coefficients: RPS_r^c = ∏(q=1 to N) [X_q]^(s_r,q): "For each reaction r and cell line c, the Reaction Propensity Score (RPS) is computed as the product of substrate concentrations each raised to their stoichiometric coefficients: RPS_r^c = ∏(q=1 to N)"
- [other] If any substrate is unmeasured, the reaction is omitted from the RPS dataset.: "If any substrate is unmeasured, the reaction is omitted from the RPS dataset."
- [other] For each metabolic reaction and each cell line, identify all substrate metabolites from the model's stoichiometry and retrieve their measured intracellular concentrations; omit reactions missing one or more substrate abundance measurements.: "For each metabolic reaction and each cell line, identify all substrate metabolites from the model's stoichiometry and retrieve their measured intracellular concentrations; omit reactions missing one"
- [other] Aggregate RPS values across biological replicates within each cell line (using median or mean as appropriate) to produce a cell-line-level RPS score for each reaction.: "Aggregate RPS values across biological replicates within each cell line (using median or mean as appropriate) to produce a cell-line-level RPS score for each reaction."
- [other] Normalize RPS scores within each reaction across all cell lines by dividing by the maximum RPS value observed for that reaction.: "Normalize RPS scores within each reaction across all cell lines by dividing by the maximum RPS value observed for that reaction."
- [intro] INTEGRATE uses intracellular metabolomics datasets and the mass action law formulation to predict how differences in substrate availability translate into differences in metabolic fluxes: "INTEGRATE uses intracellular metabolomics datasets and the mass action law formulation to predict how differences in substrate availability translate into differences in metabolic fluxes"
- [full_text] If one single reaction substrate is missing from the metabolomics measurements, the reaction is omitted from the dataset: "If one single reaction substrate is missing from the metabolomics measurements, the reaction is omitted from the dataset"
1---2name: mass-action-kinetics-formulation3description: Use when you have intracellular metabolomics concentration measurements across multiple cell lines or conditions, a stoichiometric metabolic network model with reaction-metabolite associations, and you need to predict how differences in substrate availability (independent of enzyme expression).4license: CC-BY-4.05---67# mass-action-kinetics-formulation89## Summary1011Computes Reaction Propensity Scores (RPS) from intracellular metabolomics abundance data using mass action law formulation, quantifying the expected relative metabolic flux changes across cell lines based on substrate availability alone. This provides a metabolic-regulation baseline independent of gene expression effects.1213## When to use1415Apply this skill when you have intracellular metabolomics concentration measurements across multiple cell lines or conditions, a stoichiometric metabolic network model with reaction-metabolite associations, and you need to predict how differences in substrate availability (independent of enzyme expression) translate into differences in metabolic flux propensity. Use it specifically to isolate the metabolic regulation layer (substrate-driven) from the transcriptional regulation layer (enzyme abundance-driven) for concordance analysis.1617## When NOT to use1819- When extracellular flux data or direct metabolic flux measurements (e.g., from labeled substrate tracing) are available and preferred; RPS is a proxy, not direct flux measurement.20- When substrate concentrations are not measured intracellularly or are only available in extracellular medium; model assumes intracellular substrate availability drives flux, not external concentration gradients.21- When the metabolic network lacks stoichiometric coefficients or reaction-metabolite associations; the mass action law formulation requires complete stoichiometric data.2223## Inputs2425- Intracellular metabolomics abundance dataset (metabolite concentrations per cell line, replicated measurements)26- Stoichiometric metabolic network model with reaction-metabolite associations and stoichiometric coefficients (e.g., ENGRO2 in SBML or similar format)27- Conversion mapping between metabolomics metabolite identifiers and model metabolite identifiers2829## Outputs3031- Normalized RPS dataset (table with reactions as rows, cell lines as columns, normalized RPS values ∈ [0,1] as entries)32- Pre-aggregation RPS scores (per-reaction, per-cell-line, per-replicate, before normalization)33- Metadata log of omitted reactions (those with unmeasured substrates, per cell line)3435## How to apply3637Load intracellular metabolomics abundance data (metabolite concentrations per cell line) and the stoichiometric metabolic network model (e.g., ENGRO2) with reaction-metabolite associations and stoichiometric coefficients. For each reaction and each cell line, retrieve all substrate metabolites from the model stoichiometry and their measured intracellular concentrations; omit any reaction missing one or more substrate measurements. Compute the RPS for each eligible reaction in each cell line as the product of substrate concentrations each raised to their stoichiometric coefficient power: RPS_r^c = ∏(q=1 to N) [X_q]^(s_r,q), where [X_q] is substrate q concentration and s_r,q is its stoichiometric coefficient. Aggregate RPS values across biological replicates within each cell line using median or mean as appropriate to produce a cell-line-level RPS score. Normalize RPS scores within each reaction across all cell lines by dividing by the maximum RPS value observed for that reaction, yielding dimensionless propensity scores bounded in [0,1] suitable for cross-reaction and cross-cell-line comparison.3839## Related tools4041- **ENGRO2** (Stoichiometric metabolic network model providing reaction-metabolite associations and stoichiometric coefficients for RPS computation)42- **createMetabolicDataset.py** (Preprocessing and preparation of metabolomics data for RPS calculation; quality filtering and mean computation across replicates) — https://github.com/qLSLab/integrate43- **concordanceAnalysis.py** (Post-processing integration and comparison of RPS output against RAS (transcription-derived scores) to identify metabolic vs. transcriptional regulation) — https://github.com/qLSLab/integrate44- **constraint-based stoichiometric metabolic models** (Formal scaffold providing stoichiometry, reaction definitions, and metabolite associations necessary for mass action law application)4546## Evaluation signals4748- All reactions in the output RPS table have at least one substrate measurement in at least one cell line; check that no reactions are missing that should have valid data.49- Normalized RPS scores fall within [0, 1] for all reactions and cell lines; the maximum RPS per reaction equals 1.0 (verification of normalization).50- RPS values correlate positively with intracellular substrate concentrations (mass action relationship: higher substrate → higher RPS); audit a sample of 5–10 reactions visually.51- Missing substrates cause consistent reaction omission across the pipeline; verify that reactions with unmeasured substrates are reported in the omission log and absent from the final RPS table.52- RPS concordance with RAS (Reaction Activity Scores from transcriptomics) is significantly less than perfect (Cohen κ and Pearson r < 0.8), confirming that substrate availability captures a distinct regulatory layer from gene expression.5354## Limitations5556- If any substrate is unmeasured, the entire reaction is omitted from the RPS dataset; this may remove important metabolic reactions and limits coverage to well-measured metabolites.57- RPS is a propensity score based on mass action law, not a direct measurement of metabolic flux; it does not account for enzyme kinetics (Km, Vmax), allosteric regulation, feedback inhibition, or post-translational modifications.58- Substrate concentrations are assumed to be at steady state within the 48-hour measurement window; dynamic changes or compartmentalization effects are not captured.59- RPS assumes all reactions are mass-action limited and does not distinguish reactions controlled by regulatory proteins, cofactor availability, or other non-concentration-dependent mechanisms.60- Cross-cell-line normalization by maximum RPS may obscure absolute differences in metabolic propensity; normalized scores are relative and suitable only for concordance analysis, not absolute flux prediction.6162## Evidence6364- [other] For each reaction r and cell line c, the Reaction Propensity Score (RPS) is computed as the product of substrate concentrations each raised to their stoichiometric coefficients: RPS_r^c = ∏(q=1 to N) [X_q]^(s_r,q): "For each reaction r and cell line c, the Reaction Propensity Score (RPS) is computed as the product of substrate concentrations each raised to their stoichiometric coefficients: RPS_r^c = ∏(q=1 to N)"65- [other] If any substrate is unmeasured, the reaction is omitted from the RPS dataset.: "If any substrate is unmeasured, the reaction is omitted from the RPS dataset."66- [other] For each metabolic reaction and each cell line, identify all substrate metabolites from the model's stoichiometry and retrieve their measured intracellular concentrations; omit reactions missing one or more substrate abundance measurements.: "For each metabolic reaction and each cell line, identify all substrate metabolites from the model's stoichiometry and retrieve their measured intracellular concentrations; omit reactions missing one"67- [other] Aggregate RPS values across biological replicates within each cell line (using median or mean as appropriate) to produce a cell-line-level RPS score for each reaction.: "Aggregate RPS values across biological replicates within each cell line (using median or mean as appropriate) to produce a cell-line-level RPS score for each reaction."68- [other] Normalize RPS scores within each reaction across all cell lines by dividing by the maximum RPS value observed for that reaction.: "Normalize RPS scores within each reaction across all cell lines by dividing by the maximum RPS value observed for that reaction."69- [intro] INTEGRATE uses intracellular metabolomics datasets and the mass action law formulation to predict how differences in substrate availability translate into differences in metabolic fluxes: "INTEGRATE uses intracellular metabolomics datasets and the mass action law formulation to predict how differences in substrate availability translate into differences in metabolic fluxes"70- [full_text] If one single reaction substrate is missing from the metabolomics measurements, the reaction is omitted from the dataset: "If one single reaction substrate is missing from the metabolomics measurements, the reaction is omitted from the dataset"
HolobiomicsLab 