Statistical approach for parameter identification by Turing patterns

Abstract

Prevailing theories in biological pattern formation, such as in morphogenesis or multicellular structures development, have been based on purely chemical processes, with the Turing models as the prime example. Recent studies have challenged the approach, by underlining the role of mechanical forces. A quantitative discrimination of competing theories is difficult, however, due to the elusive character of the processes: different mechanisms may result in similar patterns, while patterns obtained with a fixed model and fixed parameter values, but with small random perturbations of initial values, will significantly differ in shape, while being of the “same” type. In this sense each model parameter value corresponds to a family of patterns, rather than a fixed solution. For this situation we create a likelihood that allows a statistically sound way to distinguish the model parameters that correspond to given patterns. The method allows us to identify model parameters of reaction-diffusion systems by using Turing patterns only, i.e., the steady-state solutions of the respective equations without the use of transient data or initial values. The method is tested with three classical models of pattern formation: the FitzHugh-Nagumo model, Gierer-Meinhardt system and Brusselator reaction-diffusion system. We quantify the accuracy achieved by different amounts of training data by Bayesian sampling methods. We demonstrate how a large enough ensemble of patterns leads to detection of very small but systematic structural changes, practically impossible to distinguish with the naked eye.