Complexity reduction preserving dynamical behavior of biochemical networks.

ABSTRACT

The complexity of biochemical systems, stemming from both the large number of components and the intricate interactions between these components, may hinder us in understanding the behavior of these systems. Therefore, effective methods are required to capture their key components and interactions. Here, we present a novel and efficient reduction method to simplify mathematical models of biochemical systems. Our method is based on the exploration of the so-called admissible region, that is the set of parameters for which the mathematical model yields some required output. From the shape of the admissible region, parameters that are really required in generating the output of the system can be identified and hence retained in the model, whereas the rest is removed. To describe the idea, first the admissible region of a very small artificial network with only three nodes and three parameters is determined. Despite its simplicity, this network reveals all the basic ingredients of our reduction method. The method is then applied to an epidermal growth factor receptor (EGFR) network model. It turns out that only about 34% of the network components are required to yield the correct response to the epidermal growth factor (EGF) that was measured in the experiments, whereas the rest could be considered as redundant for this purpose. Furthermore, it is shown that parameter sensitivity on its own is not a reliable tool for model reduction, because highly sensitive parameters are not always retained, whereas slightly sensitive parameters are not always removable.