Piecewise linear models of chemical reaction networks

We show that certain non-linear dynamical systems with non-linearities in the form of Hill functions, can be approximated by piecewise linear dynamical systems. The resulting piecewise systems have closed form solutions that can be used to understand the behavior of the fully nonlinear system. We justify the reduction using geometric singular perturbation theory, and illustrate the results in networks modeling a genetic switch and a genetic oscillator.