Information geometry of chemical reaction networks: Cramer-Rao bound and absolute sensitivity revisited

Information geometry is based on classical Legendre duality but allows to incorporate additional structure such as algebraic constraints and Bregman divergence functions. It is naturally suited, and has been successfully used, to describe the thermodynamics of chemical reaction networks (CRNs) based on the Legendre duality between concentration and potential spaces, where algebraic constraints are enforced by the stoichiometry. In this article, the Riemannian geometrical aspects of the theory are explored. It is shown that duality between concentration and potential spaces and the natural parametrizations of equilibrium subspace are isometries, which leads to a multivariate Cramer-Rao bound through the comparison of two Riemannian metric tensors. In the subsequent part, the theory is applied to the recently introduced concept of absolute sensitivity. Using the Riemannian geometric tools, it is proven that the absolute sensitivity is a projection operator onto the tangent bundle of the equilibrium manifold. A linear algebraic characterization and explicit results on first order corrections to the thermodynamics of ideal solutions are provided. Finally, the theory is applied to the IDHKP-IDH glyoxylate bypass regulation system. The novelty of the theory is that it is applicable to CRNs with non-ideal thermodynamical behavior, which are prevalent in highly crowded cellular environments due to various interactions between the chemicals. Indeed, the analyzed example shows remarkable behavior ranging from hypersensitivity to negative-self regulations. These are effects which usually require strongly nonlinear reaction kinetics. However, here, they are obtained by tuning thermodynamical interactions providing a complementary, and physically well-founded, viewpoint on such phenomena.