— Models of biochemical reaction networks can be decomposed into a stoichiometric part and a kinetic part. The stoichiometric part describes the structural mass flows while the kinetic part describes how the flow rates vary with substrate concentrations and regulatory interactions. Herein a method for analyzing the robustness of biochemical networks with respect to perturbations of the kinetic part is proposed. In particular, we consider a class of perturbations that modify the local kinetic slopes while leaving the reaction flow rates in steady state unchanged. A method for computing the associated robustness radii for perturbations of single or multiple kinetic slopes is devised. The corresponding non-robust perturbations can be implemented in the original nonlinear model through specific parameter variations described by the perturbation class. The proposed method is illustrated through application to the Huang–Ferrell model of MAPK signaling cascades. In particular, we comp...
Steffen Waldherr, Frank Allgöwer, Elling W. J