Reaction systems are a formal framework for modeling processes driven by biochemical reactions. They are based on the mechanisms of facilitation and inhibition. A main assumption is that if a resource is available, then it is present in sufficient amounts and as such, several reactions using the same resource will not compete concurrently against each other; this makes reaction systems very different as a modeling framework than traditional frameworks such as ODEs or continuous time Markov chains. We demonstrate in this paper that reaction systems are rich enough to capture the essential characteristics of ODE-based models. We construct a reaction system model for the heat shock response in such a way that its qualitative behavior correlates well with the quantitative behavior of the corresponding ODE model. We construct our reaction system model based on a novel concept of dominance graph that captures the competition on resources in the ODE model. We conclude with a discussion on the expressivity of reaction systems as compared to that of ODEbased models.

Reaction systems are a new mathematical formalism inspired by the living cell and driven by only two basic mechanisms: facilitation and inhibition. As a modeling framework, they differ from the traditional approaches based on ODEs and CTMCs in two fundamental aspects: their qualitative character and the nonpermanency of resources. In this article we introduce to reaction systems several notions of central interest in biomodeling: mass conservation, invariants, steady states, stationary processes, elementary fluxes, and periodicity. We prove that the decision problems related to these properties span a number of complexity classes from P to NP-and coNP-complete to PSPACE-complete.

No abstract

Available online xxxx Communicated by M.P. JimenezReaction systems is a new mathematical formalism inspired by the biological cell, which focuses on an abstract set-based representation of chemical reactions via facilitation and inhibition. In this article we focus on the property of mass conservation for reaction systems. We show that conservation of sets gives rise to a relation between the species, which we capture in the concept of the conservation dependency graph. We then describe an application of this relation to the problem of listing all conserved sets. We further give a sufficient negative polynomial criterion which can be used for proving that a set is not conserved. Finally, we present a simulator of reaction systems, which also includes an implementation of the algorithm for listing the conserved sets of a given reaction system.

No abstract

Quantitative models may exhibit sophisticated behaviour that includes having multiple steady states, bistability, limit cycles, and period-doubling bifurcation. Such behaviour is typically driven by the numerical dynamics of the model, where the values of various numerical parameters play the crucial role. We introduce in this paper natural correspondents of these concepts to reaction systems modelling, a framework based on elementary set theoretical, forbidding/enforcing-based mechanisms. We construct several reaction systems models exhibiting these properties.

No abstract

We present the EVONANO platform for the evolution of nanomedicines with application to anti-cancer treatments. Our work aims to decrease both the time and cost required to develop nanoparticle designs. EVONANO includes a simulator to grow tumours, extract representative scenarios, and simulate nanoparticle transport through these scenarios in order to predict nanoparticle distribution. The nanoparticle designs are optimised using machine learning to efficiently find the most effective anti-cancer treatments. We demonstrate EVONANO with two examples optimising the properties of nanoparticles and treatment to selectively kill cancer cells over a range of tumour environments. Our platform shows how in silico models that capture both tumour and tissue-scale dynamics can be combined with machine learning to optimise nanomedicine.