Mechanistic models have been used for centuries to describe complex interconnected processes, including biological ones. As the scope of these models has widened, so have their computational demands. This complexity can limit their suitability when running many simulations or when real-time results are required. Surrogate machine learning (ML) models can be used to approximate the behaviour of complex mechanistic models, and once built, their computational demands are several orders of magnitude lower. This paper provides an overview of the relevant literature, both from an applicability and a theoretical perspective. For the latter, the paper focuses on the design and training of the underlying ML models. Application-wise, we show how ML surrogates have been used to approximate different mechanistic models. We present a perspective on how these approaches can be applied to models representing biological processes with potential industrial applications (e.g., metabolism and whole-cell modelling) and show why surrogate ML models may hold the key to making the simulation of complex biological systems possible using a typical desktop computer.

Agent-based models (ABMs) provide a number of advantages relative to traditional continuum modeling approaches, permitting incorporation of great detail and realism into simulations, allowing in silico tracking of single-cell behaviors and correlation of these with emergent effects at the macroscopic level. In this study we present BSim 2.0, a radically new version of BSim, a computational ABM framework for modeling dynamics of bacteria in typical experimental environments including microfluidic chemostats. This is facilitated through the implementation of new methods including cells with capsular geometry that are able to physically and chemically interact with one another, a realistic model of cellular growth, a delay differential equation solver, and realistic environmental geometries.

Organoids are engineered three-dimensional tissue cultures derived from stem cells and capable of self-renewal and self-organization into a variety of progenitors and differentiated cell types. An organoid resembles the cellular structure of an organ and retains some of its functionality, while still being amenable to in vitro experimental study. Compared with two-dimensional cultures, the three-dimensional structure of organoids provides a more realistic environment and structural organization of in vivo organs. Similarly, organoids are better suited to reproduce signaling pathway dynamics in vitro, due to a more realistic physiological environment. As such, organoids are a valuable tool to explore the dynamics of organogenesis and offer routes to personalized preclinical trials of cancer progression, invasion, and drug response. Complementary to experiments, mathematical and computational models are valuable instruments in the description of spatiotemporal dynamics of organoids. Simulations of mathematical models allow the study of multiscale dynamics of organoids, at both the intracellular and intercellular levels. Mathematical models also enable us to understand the underlying mechanisms responsible for phenotypic variation and the response to external stimulation in a cost- and time-effective manner. Many recent studies have developed laboratory protocols to grow organoids resembling different organs such as the intestine, brain, liver, pancreas, and mammary glands. However, the development of mathematical models specific to organoids remains comparatively underdeveloped. Here, we review the mathematical and computational approaches proposed so far to describe and predict organoid dynamics, reporting the simulation frameworks used and the models’ strengths and limitations.

In many biotechnological applications, it is useful for gene expression to be regulated by multiple signals, as this allows the programming of complex behavior. Here we implement, in Escherichia coli, a system that compares the concentration of two signal molecules, and tunes GFP expression proportionally to their relative abundance. The computation is performed via molecular titration between an orthogonal σ factor and its cognate anti-σ factor. We use mathematical modeling and experiments to show that the computation system is predictable and able to adapt GFP expression dynamically to a wide range of combinations of the two signals, and our model qualitatively captures most of these behaviors. We also demonstrate in silico the practical applicability of the system as a reference-comparator, which compares an intrinsic signal (reflecting the state of the system) with an extrinsic signal (reflecting the desired reference state) in a multicellular feedback control strategy.