Rules, hypergraphs, and probabilities: the three-level analysis of chemical reaction systems and other stochastic stoichiometric population processes

Abstract

We consider problems in the functional analysis and evolution of combinatorial chemical reaction networks as rule-based, or three-level systems. The first level consists of rules, realized here as graph-grammar representations of reaction mechanisms. The second level consists of stoichiometric networks of molecules and reactions, modeled as hypergraphs. At the third level is the stochastic population process on molecule counts, solved for dynamics of population trajectories or probability distributions. Earlier levels in the hierarchy generate later levels combinatorially, and as a result constraints imposed in earlier and smaller layers can propagate to impose order in the architecture or dynamics in later and larger layers. We develop general methods to study rule algebras, emphasizing system consequences of symmetry; decomposition methods of flows on hypergraphs including the stoichiometric counterpart to Kirchhoff’s current decomposition and work/dissipation relations studied in [1]; and the large-deviation theory for currents in a stoichiometric stochastic population process, deriving additive decompositions of the large-deviation function that relate a certain Kirchhoff flow decomposition to the extended Pythagorean theorem from information geometry. The latter result allows us to assign a natural probabilistic cost to topological changes in a reaction network of the kind produced by selection for catalyst-substrate specificity. We develop as an example a model of biological sugar-phosphate chemistry from a rule system published in [2]. It is one of the most potentially combinatorial reaction systems used by biochemistry, yet one in which two ancient, widespread and nearly unique pathways have evolved in the Calvin-Benson cycle and the Pentose Phosphate pathway, which are additionally nearly reverses of one another. We propose a probabilistic accounting in which physiological costs can be traded off against the fitness advantages that select them, and which suggests criteria under which these pathways may be optimal.