Enforcing Thermodynamic Consistency in Parameter Optimization and Bayesian Inference of Biochemical Reaction Networks
Abstract
Subcellular signaling pathways regulating, e.g., cell division or synaptic plasticity, are often modeled as dynamical networks of reactions satisfying the law of mass action. To remain physically grounded, these models must satisfy different physical conservation laws, such as conservation of energy. This gives rise to thermodynamic constraints formalized by the Wegscheider conditions, which ensure that the system is in detailed balance at thermodynamic equilibrium. Explicitly imposing Wegscheider conditions requires identifying cycles in the reaction network. This can become computationally cumbersome, especially in large, complex systems. Furthermore, even when these constraints are explicitly available, enforcing them requires choosing a reparameterization of the network, which may introduce asymmetries among the parameters and affect numerical conditioning. Moreover, this approach often causes inverse uncertainty quantification methods like Markov chain Monte Carlo (MCMC) to struggle. Instead of imposing Wegscheider conditions directly in the parameter space, we propose an approach for promoting physical consistency by penalizing parameters that violate detailed balance. We first illustrate the approach on a simple reaction network containing a single cycle. We show that parameter values inferred without thermodynamic constraints can reproduce the experimental data seemingly accurately, while substantially violating the Wegscheider condition. In contrast, the proposed approach preserves the agreement with the observed dynamics while yielding posterior samples that approximately satisfy the thermodynamic constraint. We further apply the approach to a model of PKA-calcineurin signaling in synaptic long-term depression, demonstrating its applicability to a more complex biochemical reaction network.
Related
- Multi-parameter optimization of recombinant protein secretion in Bacillus subtilis: a review.
- Parameter optimization of IPT systems based on genetic algorithms
- Design and parameter optimization of a bilateral staggered wheel-spoon precision seed-metering device for millet plot breeding
- NR2F6 functions as a key transcriptional checkpoint enforcing dysfunction of human CD3+ T cells and limiting CAR-T cell potency
- How Institutional Distrust Breeds One-Sided Leadership: A Consistency-Based Account
BibTeX
@misc{milinanni2026enforcing,
title = {{Enforcing Thermodynamic Consistency in Parameter Optimization and Bayesian Inference of Biochemical Reaction Networks}},
author = {Federica Milinanni and Andrei Kramer-Miehe and Jan Glaubitz and Olivia Eriksson},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.17155},
eprint = {2610.17155},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.17155},
}
RIS
TY - UNPB
TI - Enforcing Thermodynamic Consistency in Parameter Optimization and Bayesian Inference of Biochemical Reaction Networks
AU - Federica Milinanni
AU - Andrei Kramer-Miehe
AU - Jan Glaubitz
AU - Olivia Eriksson
PY - 2026
DA - 2026/10/12
PB - arXiv
DO - 10.48550/arxiv.2610.17155
UR - https://arxiv.org/abs/2610.17155
AB - Subcellular signaling pathways regulating, e.g., cell division or synaptic plasticity, are often modeled as dynamical networks of reactions satisfying the law of mass action. To remain physically grounded, these models must satisfy different physical conservation laws, such as conservation of energy. This gives rise to thermodynamic constraints formalized by the Wegscheider conditions, which ensure that the system is in detailed balance at thermodynamic equilibrium. Explicitly imposing Wegscheider conditions requires identifying cycles in the reaction network. This can become computationally cumbersome, especially in large, complex systems. Furthermore, even when these constraints are explicitly available, enforcing them requires choosing a reparameterization of the network, which may introduce asymmetries among the parameters and affect numerical conditioning. Moreover, this approach often causes inverse uncertainty quantification methods like Markov chain Monte Carlo (MCMC) to struggle. Instead of imposing Wegscheider conditions directly in the parameter space, we propose an approach for promoting physical consistency by penalizing parameters that violate detailed balance. We first illustrate the approach on a simple reaction network containing a single cycle. We show that parameter values inferred without thermodynamic constraints can reproduce the experimental data seemingly accurately, while substantially violating the Wegscheider condition. In contrast, the proposed approach preserves the agreement with the observed dynamics while yielding posterior samples that approximately satisfy the thermodynamic constraint. We further apply the approach to a model of PKA-calcineurin signaling in synaptic long-term depression, demonstrating its applicability to a more complex biochemical reaction network.
ER -