efficient estimation of the robustness region of biological models with oscillatory behavior有效估计的鲁棒性地区的生物模型与振荡行为.pdfVIP

Efficient Estimation of the Robustness Region of Biological Models with Oscillatory Behavior 1,2 1,2 1 1 Mochamad Apri *, Jaap Molenaar , Maarten de Gee , George van Voorn 1 Biometris, Department of Mathematical and Statistical Methods, Wageningen University, Wageningen, The Netherlands, 2 Netherlands Consortium for Systems Biology, Amsterdam, The Netherlands Abstract Robustness is an essential feature of biological systems, and any mathematical model that describes such a system should reflect this feature. Especially, persistence of oscillatory behavior is an important issue. A benchmark model for this phenomenon is the Laub-Loomis model, a nonlinear model for cAMP oscillations in Dictyostelium discoideum. This model captures the most important features of biomolecular networks oscillating at constant frequencies. Nevertheless, the robustness of its oscillatory behavior is not yet fully understood. Given a system that exhibits oscillating behavior for some set of parameters, the central question of robustness is how far the parameters may be changed, such that the qualitative behavior does not change. The determination of such a ‘‘robustness region’’ in parameter space is an intricate task. If the number of parameters is high, it may be also time consuming. In the literature, several methods are proposed that partially tackle this problem. For example, some methods only detect particular bifurcations, or only find a relatively small box- shaped estimate for an irregularly shaped robustness region. Here, we present an approach that is much more general, and is especially designed to be efficient for systems with a large number of parameters. As an illustration, we apply the method first to a well