Mutual Linearity of Complexes in Chemical Reaction Networks

Abstract

Exact statements about how a system far from equilibrium responds to changes in its rate constants are scarce beyond linear response, yet they are central to the control of biochemical networks. A notable exception is mutual linearity: in a jump process, controlling the rates of a single transition makes any two stationary probabilities obey an exact affine relation, arbitrarily far from equilibrium. With nonlinear mass-action kinetics, and conservation laws in place of a normalization, chemical reaction networks seem to defy mutual linearity since species concentrations obey no such relation. We show that it was there all along, in the complexes' activities (mass-action monomials) rather than in the concentrations of species. For zero-deficiency networks with a single linkage class, controlling one reaction makes any three stationary activities obey a linear relation whose coefficients are independent of the controlled rate constants, the conserved quantities (hence the initial state), and even of the stoichiometry of complexes. This predicts hidden activities and bounds the activity ratios the control can reach, while a violation signals nonzero deficiency or multiple linkage classes.