Deterministic function computation with chemical reaction networks

Abstract

Chemical reaction networks (CRNs) formally model chemistry in a well-mixed solution. CRNs are widely used to describe information processing occurring in natural cellular regulatory networks, and with upcoming advances in synthetic biology, CRNs are a promising language for the design of artificial molecular control circuitry. Nonetheless, despite the widespread use of CRNs in the natural sciences, the range of computational behaviors exhibited by CRNs is not well understood. CRNs have been shown to be efficiently Turing-universal (i.e., able to simulate arbitrary algorithms) when allowing for a small probability of error. CRNs that are guaranteed to converge on a correct answer, on the other hand, have been shown to decide only the semilinear predicates (a multi-dimensional generalization of “eventually periodic” sets). We introduce the notion of function, rather than predicate, computation by representing the output of a function (Formula presented.) by a count of some molecular species, i.e., if the CRN starts with (Formula presented.) molecules of some “input” species (Formula presented.) the CRN is guaranteed to converge to having (Formula presented.) molecules of the “output” species (Formula presented.). We show that a function (Formula presented.) is deterministically computed by a CRN if and only if its graph (Formula presented.) is a semilinear set. Finally, we show that each semilinear function f (a function whose graph is a semilinear set) can be computed by a CRN on input x in expected time (Formula presented.).