Stochastic modelling of the eukaryotic heat shock response

Andrzej Mizera (1, 2) , Barbara Gambin (1)

1 PAN - Polska Akademia Nauk = Polish Academy of Sciences = Académie polonaise des sciences
2 Turku Centre for Computer Science

The heat shock response (HSR) is a highly evolutionarily conserved defence mechanism allowing the cell to promptly react to elevated temperature conditions and other forms of stress. It has been subject to intense research for at least two main reasons. First, it is considered a promising candidate for deciphering the engineering principles underlying regulatory networks. Second, heat shock proteins (main actors of the HSR) play crucial role in many fundamental cellular processes. Therefore, profound understanding of the heat shock response would have far-reaching ramifications for the cell biology. Recently, a new deterministic model of the eukaryotic heat shock response has been proposed in the literature. It is very attractive since it consists of only the minimum number of components required by any functional regulatory network, while yet being capable of biological validation. However, it admits small molecule populations of some of the considered metabolites. In this paper a stochastic model corresponding to the deterministic one is constructed and the outcomes of these two models are confronted. The aim with this comparison is two show that, in the case of the heat shock response, the approximation of a discrete system with a continuous model is a reasonable approach. This is not always the truth, especially when the numbers of molecules of the considered species are small. By making the effort of performing and analysing 1000 stochastic simulations, we investigate the range of behaviour the stochastic model is likely to exhibit. We demonstrate that the obtained results agree well with the dynamics displayed by the continuous model, which strengthens the trust in the deterministic description. A proof of the existence and uniqueness of the stationary distribution of the Markov chain underlying the stochastic model is given. Moreover, the obtained view of the stochastic dynamics and the performed comparison to the outcome of the continuous formulation provide more insight into the dynamics of the heat shock response mechanism.

Mots clés: Mathematical modelling; Stochastic model; Computer simulations; Markov chain; Gillespie algorithm; Stationary distribution; Steady-state; Clustering

Andrzej Mizera, Barbara Gambin. Stochastic modelling of the eukaryotic heat shock response. Journal of Theoretical Biology, 2010, 265 (3), pp.455. ⟨10.1016/j.jtbi.2010.04.029⟩. ⟨hal-00608418⟩