Ultrasensitive dual phosphorylation dephosphorylation cycle kinetics exhibits canonical competition behavior

We establish a mathematical model for a cellular biochemical signaling module in terms of a planar differential equation system. The signaling process is carried out by two phosphorylation-dephosphorylation reaction steps that share common kinase and phosphatase with saturated enzyme kinetics. The pair of equations is particularly simple in the present mathematical formulation, but they are singular. A complete mathematical analysis is developed based on an elementary perturbation theory. The dynamics exhibits the canonical competition behavior in addition to bistability. Although widely understood in ecological context, we are not aware of a full range of biochemical competition in a simple signaling network. The competition dynamics has broad implications to cellular processes such as cell differentiation and cancer immunoediting. The concepts of homogeneous and heterogeneous multisite phosphorylation are introduced and their corresponding dynamics are compared: there is no bistability in a heterogeneous dual phosphorylation system. A stochastic interpretation is also provided that further gives intuitive understanding of the bistable behavior inside the cells.