Studying the effects of adding spatiality to a process algebra model

We use NetLogo to create simulations of two models of disease transmission originally expressed in WSCCS. This allows us to introduce spatiality into the models and explore the consequences of having different contact structures among the agents. In previous work, mean field equations were derived from the WSCCS models, giving a description of the aggregate behaviour of the overall population of agents. These results turned out to differ from results obtained by another team using cellular automata models, which differ from process algebra by being inherently spatial. By using NetLogo we are able to explore whether spatiality, and resulting differences in the contact structures in the two kinds of models, are the reason for this different results. Our tentative conclusions, based at this point on informal observations of simulation results, are that space does indeed make a big difference. If space is ignored and individuals are allowed to mix randomly, then the simulations yield results that closely match the mean field equations, and consequently also match the associated global transmission terms (explained below). At the opposite extreme, if individuals can only contact their immediate neighbours, the simulation results are very different from the mean field equations (and also do not match the global transmission terms). These results are not surprising, and are consistent with other cellular automata-based approaches. We found that it was easy and convenient to implement and simulate the WSCCS models within NetLogo, and we recommend this approach to anyone wishing to explore the effects of introducing spatiality into a process algebra model

Multi-scale modelling of biological systems in process algebra

There is a growing interest in combining different levels of detail of biological phenomena into unique multi-scale models that represent both biochemical details and higher order structures such as cells, tissues or organs. The state of the art of multi-scale models presents a variety of approaches often tailored around specific problems and composed of a combination of mathematical techniques. As a result, these models are difficult to build, compose, compare and analyse. In this thesis we identify process algebra as an ideal formalism to multi-scale modelling of biological systems. Building on an investigation of existing process algebras, we define process algebra with hooks (PAH), designed to be a middle-out approach to multi-scale modelling. The distinctive features of PAH are: the presence of two synchronisation operators, distinguishing interactions within and between scales, and composed actions, representing events that occur at multiple scales. A stochastic semantics is provided, based on functional rates derived from kinetic laws. A parametric version of the algebra ensures that a model description is compact. This new formalism allows for: unambiguous definition of scales as processes and interactions within and between scales as actions, compositionality between scales using a novel vertical cooperation operator and compositionality within scales using a traditional cooperation operator, and relating models and their behaviour using equivalence relations that can focus on specified scales. Finally, we apply PAH to define, compose and relate models of pattern formation and tissue growth, highlighting the benefits of the approach

Relating PDEs in Cylindrical Coordinates and CTMCs with Levels of Concentration

We present the derivation of a CTMC with levels model of diusion in cylindrical coordinates from the partial dierential equation for Fick's law. The resulting model abstracts both molar concentration, by discrete levels, and spatial location, by discrete compartments. We apply the results to the diusion of nitric oxide in human vessels and illustrate with simulations in the PRISM tool