The simplest way to account for the influence of diffusion on the kinetics of multisite phosphorylation is to modify the rate constants in the conventional rate equations of chemical kinetics. We have previously shown that this is not enough and new transitions between the reactants must also be introduced. Here we extend our results by considering enzymes that are inactive after modifying the substrate and need time to become active again. This generalization leads to a surprising result. The introduction of enzyme reactivation results in a diffusion-modified kinetic scheme with a new transition that has a negative rate constant. The reason for this is that mapping non-Markovian rate equations onto Markovian ones with time-independent rate constants is not a good approximation at short times. We then developed a non-Markovian theory that involves memory kernels instead of rate constants. This theory is now valid at short times, but is more challenging to use. As an example, the diffusion-modified kinetic scheme with new connections was used to calculate kinetics of double phosphorylation and steady-state response in a phosphorylation-dephosphorylation cycle. We have reproduced the loss of bistability in the phosphorylation-dephosphorylation cycle when the enzyme reactivation time decreases, which was obtained by particle-based computer simulations.
Read full abstract