2At the crossroad between biology and mathematical modelling, computational systems biology 3 can contribute to a mechanistic understanding of high-level biological phenomenon. But as our 4 knowledge accumulates, the size and complexity of mathematical models increase, calling for the 5 development of efficient dynamical analysis methods. 6 Here, we take advantage of generic computational techniques to enable the dynamical study of 7 complex cellular network models. 8 A first approach, called "model verification" and inspired by unitary testing in software 9 development, enables the formalisation and automated verification of validation criteria for 10 whole models or selected sub-parts, thereby greatly facilitating model development. 11 A second approach, called "value propagation", enables efficient analytical computation of the 12 impact of specific environmental or genetic conditions on model dynamics. 13 We apply these two approaches to the delineation and the analysis of a comprehensive model 14 for T cell activation, taking into account CTLA4 and PD-1 checkpoint inhibitory pathways. While 15 the use of model verification greatly eased the delineation of logical rules complying with a set 16 of dynamical specifications, the use of value propagation provided interesting insights into the 17 different potential of CTLA4 and PD-1 immunotherapies. 18 Both methods are implemented and made available in the all-inclusive CoLoMoTo Docker image, 19 while the different steps of the model analysis are fully reported in two companion interactive 20 jupyter notebooks, thereby ensuring the reproduction of our results.
21Recent technical developments have allowed scientists to study immunology and health-related issues 23 from a variety of angles. For many diseases, especially for cancer, the current trend consists in aggregating 24 data coming from different sources to gain a global view of cell, tissue, or organ dysfunction. Over the 25 last decades, diverse mathematical frameworks have been proposed to seize a multiplicity of biological 26 questions (Eftimie et al., 2016; Chakraborty, 2017). However, the increasing complexity of biological 27 questions implies the development of more sophisticated models, which in turn bring serious computational 28 challenges.
29Among the mathematical approaches proposed for the modelling of cellular networks, the logical 30 modelling framework appeared particularly adapted (see e.g. (Thomas, 1991)). In particular, it has been 31 successfully applied to immunology and cancer, leading to the creation of models encompassing dozens 32 of components, some including many inputs components (Grieco et al., 2013; Abou-Jaoudé et al., 2014; 33 Flobak et al., 2015; Oyeyemi et al., 2015). However, the large size of these models hinders the complete 34 exploration of their dynamical behaviour through simulation, especially in non-deterministic settings.
35To address these difficulties, we hereby propose the application of two generic computational methods.
36First, we d...