Effective sampling trajectory optimisation for sensitivity analysis of biological systems

Sensitivity analysis has been widely applied to study the biological systems, including metabolic networks, signalling pathways, and genetic circuits. The Morris method is a kind of screening sensitivity analysis approach, which can fast identify a few key factors from numerous biological parameters and inputs. The parameter or input space is randomly sampled to produce a very limited number of trajectories for the calculation of elementary effects. It is clear that the sampled trajectories are not enough to cover the whole uncertain space, which eventually causes unstable sensitivity measures. This paper presents a novel trajectory optimisation algorithm for the Morris-based sensitivity calculation to ensure a good scan throughout the whole uncertain space. The paper demonstrates that this presented method gets more consistent sensitivity results through a benchmark example. The application to a previously published ordinary differential equation model of a cellular signalling network is presented. In detail, the parameter sensitivity analysis verifies the good agreement with data of the literatures.