ChemChaste: Simulating spatially inhomogenous biochemical reaction-diffusion systems for modelling cell-environment feedbacks

Abstract

MotivationSpatial organisation plays an important role in the function of many biological systems, from cell fate specification in animal development to multi-step metabolic conversions in microbial communities. The study of such systems benefits from the use of spatially explicit computational models that combine a discrete description of cells with a continuum description of one or more chemicals diffusing within a surrounding bulk medium. These models allow the in silico testing and refinement of mechanistic hypotheses. However, most existing models of this type do not account for concurrent bulk and intracellular biochemical reactions and their possible coupling.ResultsHere, we describe ChemChaste, an extension for the open-source C++ computational biology library Chaste. ChemChaste enables the spatial simulation of both multicellular and bulk biochemistry by expanding on Chaste’s existing capabilities. In particular, ChemChaste enables: (i) simulation of an arbitrary number of spatially diffusing chemicals; (ii) spatially heterogeneous chemical diffusion coefficients; and (iii) inclusion of both bulk and intracellular biochemical reactions and their coupling. ChemChaste also introduces a file-based interface that allows users to define the parameters relating to these functional features without the need to interact directly with Chaste’s core C++ code. We describe ChemChaste and demonstrate its functionality using a selection of chemical and biochemical exemplars, with a focus on demonstrating increased ability in modelling bulk chemical reactions and their coupling with intracellular reactions.Availability and implementationChemChaste is a free, open-source C++ library, available via GitHub at https://github.com/OSS-Lab/ChemChaste under the BSD license.ContactO.Soyer@warwick.ac.uk or a.g.fletcher@sheffield.ac.uk

ChemChaste: Simulating spatially inhomogeneous biochemical reaction–diffusion systems for modeling cell–environment feedbacks

“…Snapshots of our code and other data further supporting this work are openly available in the GigaScience repository GigaDB [ 51 ].…”