Free and Interfacial Boundaries in Individual-Based Models of Multicellular Biological systems.

Coordination of cell behaviour is key to a myriad of biological processes including tissue morphogenesis, wound healing, and tumour growth. As such, individual-based computational models, which explicitly describe inter-cellular interactions, are commonly used to model collective cell dynamics. However, when using individual-based models, it is unclear how descriptions of cell boundaries affect overall population dynamics. In order to investigate this we define three cell boundary descriptions of varying complexities for each of three widely used off-lattice individual-based models: overlapping spheres, Voronoi tessellation, and vertex models. We apply our models to multiple biological scenarios to investigate how cell boundary description can influence tissue-scale behaviour. We find that the Voronoi tessellation model is most sensitive to changes in the cell boundary description with basic models being inappropriate in many cases. The timescale of tissue evolution when using an overlapping spheres model is coupled to the boundary description. The vertex model is demonstrated to be the most stable to changes in boundary description, though still exhibits timescale sensitivity. When using individual-based computational models one should carefully consider how cell boundaries are defined. To inform future work, we provide an exploration of common individual-based models and cell boundary descriptions in frequently studied biological scenarios and discuss their benefits and disadvantages.

Mathematical models of developmental vascular remodelling: A review.

Over the past 40 years, there has been a strong focus on the development of mathematical models of angiogenesis, while developmental remodelling has received little such attention from the mathematical community. Sprouting angiogenesis can be seen as a very crude way of laying out a primitive vessel network (the raw material), while remodelling (understood as pruning of redundant vessels, diameter control, and the establishment of vessel identity and hierarchy) is the key to turning that primitive network into a functional network. This multiscale problem is of prime importance in the development of a functional vasculature. In addition, defective remodelling (either during developmental remodelling or due to a reactivation of the remodelling programme caused by an injury) is associated with a significant number of diseases. In this review, we discuss existing mathematical models of developmental remodelling and explore the important contributions that these models have made to the field of vascular development. These mathematical models are effectively used to investigate and predict vascular development and are able to reproduce experimentally observable results. Moreover, these models provide a useful means of hypothesis generation and can explain the underlying mechanisms driving the observed structural and functional network development. However, developmental vascular remodelling is still a relatively new area in mathematical biology, and many biological questions remain unanswered. In this review, we present the existing modelling paradigms and define the key challenges for the field.

Choice of spatial discretisation influences the progression of viral infection within multicellular tissues.

There has been an increasing recognition of the utility of models of the spatial dynamics of viral spread within tissues. Multicellular models, where cells are represented as discrete regions of space coupled to a virus density surface, are a popular approach to capture these dynamics. Conventionally, such models are simulated by discretising the viral surface and depending on the rate of viral diffusion and other considerations, a finer or coarser discretisation may be used. The impact that this choice may have on the behaviour of the system has not been studied. Here we demonstrate that under realistic parameter regimes - where viral diffusion is small enough to support the formation of familiar ring-shaped infection plaques - the choice of spatial discretisation of the viral surface can qualitatively change key model outcomes including the time scale of infection. Importantly, we show that the choice between implementing viral spread as a cell-scale process, or as a high-resolution converged PDE can generate distinct model outcomes, which raises important conceptual questions about the strength of assumptions underpinning the spatial structure of the model. We investigate the mechanisms driving these discretisation artefacts, the impacts they may have on model predictions, and provide guidance on the design and implementation of spatial and especially multicellular models of viral dynamics. We obtain our results using the simplest TIV construct for the viral dynamics, and therefore anticipate that the important effects we describe will also influence model predictions in more complex models of virus-cell-immune system interactions. This analysis will aid in the construction of models for robust and biologically realistic modelling and inference.

Multiscale modelling of desquamation in the interfollicular epidermis.

Maintenance of epidermal thickness is critical to the barrier function of the skin. Decreased tissue thickness, specifically in the stratum corneum (the outermost layer of the tissue), causes discomfort and inflammation, and is related to several severe diseases of the tissue. In order to maintain both stratum corneum thickness and overall tissue thickness it is necessary for the system to balance cell proliferation and cell loss. Cell proliferation in the epidermis occurs in the basal layer and causes constant upwards movement in the tissue. Cell loss occurs when dead cells at the top of the tissue are lost to the environment through a process called desquamation. Desquamation is thought to occur through a gradual reduction in adhesion between cells, due to the cleaving of adhesion proteins by enzymes, in the stratum corneum. In this paper we will investigate combining a (mass action) subcellular model of desquamation with a three dimensional (cell centre based) multicellular model of the interfollicular epidermis to better understand maintenance of epidermal thickness. Specifically, our aim is to determine if a hypothesised biological model for the degradation of cell-cell adhesion, from the literature, is sufficient to maintain a steady state tissue thickness. These investigations show the model is able to provide a consistent rate of cell loss in the multicellular model. This loss balances proliferation, and hence maintains a homeostatic tissue thickness. Moreover, we find that multiple proliferative cell populations in the basal layer can be represented by a single proliferative cell population, simplifying investigations with this model. The model is used to investigate a disorder (Netherton Syndrome) which disrupts desquamation. The model shows how biochemical changes can cause disruptions to the tissue, resulting in a reduced tissue thickness and consequently diminishing the protective role of the tissue. A hypothetical treatment result is also investigated: we compare the cases of a partially effective homogeneous treatment (where all cells partially recover) and a totally effective heterogeneous treatment (in which a proportion of the cells totally recover) with the aim to determine the difference in the response of the tissue to these different scenarios. Results show an increased benefit to corneum thickness from the heterogeneous treatment over the homogeneous treatment.

Modelling realistic 3D deformations of simple epithelia in dynamic homeostasis.

The maintenance of tissue and organ structures during dynamic homeostasis is often not well understood. In order for a system to be stable, cell renewal, cell migration and cell death must be finely balanced. Moreover, a tissue's shape must remain relatively unchanged. Simple epithelial tissues occur in various structures throughout the body, such as the endothelium, mesothelium, linings of the lungs, saliva and thyroid glands, and gastrointestinal tract. Despite the prevalence of simple epithelial tissues, there are few models which accurately describe how these tissues maintain a stable structure. Here, we present a novel, 3D, deformable, multilayer, cell-centre model of a simple epithelium. Cell movement is governed by the minimisation of a bending potential across the epithelium, cell-cell adhesion, and viscous effects. We show that the model is capable of maintaining a consistent tissue structure while undergoing self renewal. We also demonstrate the model's robustness under tissue renewal, cell migration and cell removal. The model presented here is a valuable advancement towards the modelling of tissues and organs with complex and generalised structures.

Traffic Patterns of the Migrating Endothelium: How Force Transmission Regulates Vascular Malformation and Functional Shunting During Angiogenic Remodelling.

Angiogenesis occurs in distinct phases: initial spouting is followed by remodelling in which endothelial cells (ECs) composing blood vessels rearrange by migrating against the direction of flow. Abnormal remodelling can result in vascular malformation. Such is the case in mutation of the Alk1 receptor within the mouse retina which disrupts flow-migration coupling, creating mixed populations of ECs polarised with/against flow which aggregate into arteriovenous malformations (AVMs). The lack of live imaging options in vivo means that the collective EC dynamics that drive AVM and the consequences of mixed populations of polarity remain a mystery. Therefore, our goal is to present a novel agent-based model to provide theoretical insight into EC force transmission and collective dynamics during angiogenic remodelling. Force transmission between neighbouring agents consists of extrusive forces which maintain spacing and cohesive forces which maintain the collective. We performed migration simulations within uniformly polarised populations (against flow) and mixed polarity (with/against flow). Within uniformly polarised populations, extrusive forces stabilised the plexus by facilitating EC intercalation which ensures that cells remained evenly distributed. Excess cohesion disrupts intercalation, resulting in aggregations of cells and functional shunting. Excess cohesion between ECs prevents them from resolving diameter balances within the plexus, leading to prolonged flow reversals which exert a critical behaviour change within the system as they switch the direction of cell migration and traffic patterns at bifurcations. Introducing mixtures of cell polarity dramatically changed the role of extrusive forces within the system. At low extrusion, opposing ECs were able to move past each other; however, at high extrusion the pushing between cells resulted in migration speeds close to zero, forming traffic jams and disrupting migration. In our study, we produced vascular malformations and functional shunting with either excess cohesion between ECs or mixtures of cell polarity. At the centre of both these mechanisms are cell-cell adherens junctions, which are involved in flow sensing/polarity and must remodelling dynamically to allow rearrangements of cells during vascular patterning. Thus, our findings implicate junctional dysfunction as a new target in the treatment and prevention of vascular disease and AVMs.

How quickly does a wound heal? Bayesian calibration of a mathematical model of venous leg ulcer healing.

Chronic wounds, such as venous leg ulcers, are difficult to treat and can reduce the quality of life for patients. Clinical trials have been conducted to identify the most effective venous leg ulcer treatments and the clinical factors that may indicate whether a wound will successfully heal. More recently, mathematical modelling has been used to gain insight into biological factors that may affect treatment success but are difficult to measure clinically, such as the rate of oxygen flow into wounded tissue. In this work, we calibrate an existing mathematical model using a Bayesian approach with clinical data for individual patients to explore which clinical factors may impact the rate of wound healing for individuals. Although the model describes group-level behaviour well, it is not able to capture individual-level responses in all cases. From the individual-level analysis, we propose distributions for coefficients of clinical factors in a linear regression model, but ultimately find that it is difficult to draw conclusions about which factors lead to faster wound healing based on the existing model and data. This work highlights the challenges of using Bayesian methods to calibrate partial differential equation models to individual patient clinical data. However, the methods used in this work may be modified and extended to calibrate spatiotemporal mathematical models to multiple data sets, such as clinical trials with several patients, to extract additional information from the model and answer outstanding biological questions.

Seven challenges in the multiscale modeling of multicellular tissues.

The growth and dynamics of multicellular tissues involve tightly regulated and coordinated morphogenetic cell behaviors, such as shape changes, movement, and division, which are governed by subcellular machinery and involve coupling through short- and long-range signals. A key challenge in the fields of developmental biology, tissue engineering and regenerative medicine is to understand how relationships between scales produce emergent tissue-scale behaviors. Recent advances in molecular biology, live-imaging and ex vivo techniques have revolutionized our ability to study these processes experimentally. To fully leverage these techniques and obtain a more comprehensive understanding of the causal relationships underlying tissue dynamics, computational modeling approaches are increasingly spanning multiple spatial and temporal scales, and are coupling cell shape, growth, mechanics, and signaling. Yet such models remain challenging: modeling at each scale requires different areas of technical skills, while integration across scales necessitates the solution to novel mathematical and computational problems. This review aims to summarize recent progress in multiscale modeling of multicellular tissues and to highlight ongoing challenges associated with the construction, implementation, interrogation, and validation of such models. This article is categorized under: Reproductive System Diseases > Computational Models Metabolic Diseases > Computational Models Cancer > Computational Models.

Modelling Cellular Interactions and Dynamics During Kidney Morphogenesis.

Kidney disease and renal disorders account for a significant proportion of health complications in mid-late adulthood worldwide. Many renal deficiencies are due to improper formation of the kidneys before birth, which are caused by disorders in the developmental process that arise from genetic and/or environmental factors. Mathematical modelling can help build on experimental knowledge to increase our understanding of the complexities of kidney organogenesis. In this paper, we present a discrete cell-based model of kidney development. Specifically, we model the tip of the developing ureteric tree to investigate the behaviours of cap mesenchyme cells which are required to sustain ureteric tip growth. We find that spatial regulation of the differentiation of cap mesenchyme cells through cellular signalling is sufficient to ensure robust ureteric tip development. Additionally, we find that increased adhesion interactions between cap mesenchyme cells and the ureteric tip surface can lead to a more stable tip-cap unit. Our analysis of the various processes on this scale highlights essential components for healthy kidney growth and provides insight into mechanisms to be studied further in order to replicate the process in vitro.

Maintaining the proliferative cell niche in multicellular models of epithelia.

The maintenance of the proliferative cell niche is critical to epithelial tissue morphology and function. In this paper we investigate how current modelling methods can result in the erroneous loss of proliferative cells from the proliferative cell niche. Using an established model of the inter-follicular epidermis we find there is a limit to the proliferative cell densities that can be maintained in the basal layer (the niche) if we do not include additional mechanisms to stop the loss of proliferative cells from the niche. We suggest a new methodology that enables maintenance of a desired homeostatic population of proliferative cells in the niche: a rotational force is applied to the two daughter cells during the mitotic phase of division to enforce a particular division direction. We demonstrate that this new methodology achieves this goal. This methodology reflects the regulation of the orientation of cell division.

A mathematical model of cell fate selection on a dynamic tissue.

Multicellular tissues are the building blocks of many biological systems and organs. These tissues are not static, but dynamically change over time. Even if the overall structure remains the same there is a turnover of cells within the tissue. This dynamic homeostasis is maintaned by numerous governing mechanisms which are finely tuned in such a way that the tissue remains in a homeostatic state, even across large timescales. Some of these governing mechanisms include cell motion, and cell fate selection through inter cellular signalling. However, it is not yet clear how to link these two processes, or how they may affect one another across the tissue. In this paper, we present a multicellular, multiscale model, which brings together the two phenomena of cell motility, and inter cellular signalling, to describe cell fate selection on a dynamic tissue. We find that the affinity for cellular signalling to occur greatly influences a cells ability to differentiate. We also find that our results support claims that cell differentiation is a finely tuned process within dynamic tissues at homeostasis, with excessive cell turnover rates leading to unhealthy (undifferentiated and unpatterned) tissues.

A rigid body framework for multicellular modeling.

Off-lattice models are a well-established approach in multicellular modeling, where cells are represented as points that are free to move in space. The representation of cells as point objects is useful in a wide range of settings, particularly when large populations are involved; however, a purely point-based representation is not naturally equipped to deal with objects that have length, such as cell boundaries or external membranes. Here we introduce an off-lattice modeling framework that exploits rigid body mechanics to represent objects using a collection of conjoined one-dimensional edges in a viscosity-dominated system. This framework can be used to represent cells as free moving polygons, to allow epithelial layers to smoothly interact with themselves, to model rod-shaped cells such as bacteria and to robustly represent membranes. We demonstrate that this approach offers solutions to the problems that limit the scope of current off-lattice multicellular models.

A modular framework for multiscale, multicellular, spatiotemporal modeling of acute primary viral infection and immune response in epithelial tissues and its application to drug therapy timing and effectiveness.

Simulations of tissue-specific effects of primary acute viral infections like COVID-19 are essential for understanding disease outcomes and optimizing therapies. Such simulations need to support continuous updating in response to rapid advances in understanding of infection mechanisms, and parallel development of components by multiple groups. We present an open-source platform for multiscale spatiotemporal simulation of an epithelial tissue, viral infection, cellular immune response and tissue damage, specifically designed to be modular and extensible to support continuous updating and parallel development. The base simulation of a simplified patch of epithelial tissue and immune response exhibits distinct patterns of infection dynamics from widespread infection, to recurrence, to clearance. Slower viral internalization and faster immune-cell recruitment slow infection and promote containment. Because antiviral drugs can have side effects and show reduced clinical effectiveness when given later during infection, we studied the effects on progression of treatment potency and time-of-first treatment after infection. In simulations, even a low potency therapy with a drug which reduces the replication rate of viral RNA greatly decreases the total tissue damage and virus burden when given near the beginning of infection. Many combinations of dosage and treatment time lead to stochastic outcomes, with some simulation replicas showing clearance or control (treatment success), while others show rapid infection of all epithelial cells (treatment failure). Thus, while a high potency therapy usually is less effective when given later, treatments at late times are occasionally effective. We illustrate how to extend the platform to model specific virus types (e.g., hepatitis C) and add additional cellular mechanisms (tissue recovery and variable cell susceptibility to infection), using our software modules and publicly-available software repository.

A modular framework for multiscale, multicellular, spatiotemporal modeling of acute primary viral infection and immune response in epithelial tissues and its application to drug therapy timing and effectiveness: A multiscale model of viral infection in epithelial tissues.

Simulations of tissue-specific effects of primary acute viral infections like COVID-19 are essential for understanding disease outcomes and optimizing therapies. Such simulations need to support continuous updating in response to rapid advances in understanding of infection mechanisms, and parallel development of components by multiple groups. We present an open-source platform for multiscale spatiotemporal simulation of an epithelial tissue, viral infection, cellular immune response and tissue damage, specifically designed to be modular and extensible to support continuous updating and parallel development. The base simulation of a simplified patch of epithelial tissue and immune response exhibits distinct patterns of infection dynamics from widespread infection, to recurrence, to clearance. Slower viral internalization and faster immune-cell recruitment slow infection and promote containment. Because antiviral drugs can have side effects and show reduced clinical effectiveness when given later during infection, we studied the effects on progression of treatment potency and time-of-first treatment after infection. In simulations, even a low potency therapy with a drug which reduces the replication rate of viral RNA greatly decreases the total tissue damage and virus burden when given near the beginning of infection. Many combinations of dosage and treatment time lead to stochastic outcomes, with some simulation replicas showing clearance or control (treatment success), while others show rapid infection of all epithelial cells (treatment failure). Thus, while a high potency therapy usually is less effective when given later, treatments at late times are occasionally effective. We illustrate how to extend the platform to model specific virus types (e.g., hepatitis C) and add additional cellular mechanisms (tissue recovery and variable cell susceptibility to infection), using our software modules and publicly-available software repository.

Modelling the effect of subcellular mutations on the migration of cells in the colorectal crypt.

BACKGROUND: Many cancers arise from mutations in cells within epithelial tissues. Mutations manifesting at the subcellular level influence the structure and function of the tissue resulting in cancer. Previous work has proposed how cell level properties can lead to mutant cell invasion, but has not incorporated detailed subcellular modelling RESULTS: We present a framework that allows the straightforward integration and simulation of SBML representations of subcellular dynamics within multiscale models of epithelial tissues. This allows us to investigate the effect of mutations in subcellular pathways on the migration of cells within the colorectal crypt. Using multiple models we find that mutations in APC, a key component in the Wnt signalling pathway, can bias neutral drift and can also cause downward invasion of mutant cells in the crypt. CONCLUSIONS: Our framework allows us to investigate how subcellular mutations, i.e. knockouts and knockdowns, affect cell-level properties and the resultant migration of cells within epithelial tissues. In the context of the colorectal crypt, we see that mutations in APC can lead directly to mutant cell invasion.

Chaste: Cancer, Heart and Soft Tissue Environment.

Chaste (Cancer, Heart And Soft Tissue Environment) is an open source simulation package for the numerical solution of mathematical models arising in physiology and biology. To date, Chaste development has been driven primarily by applications that include continuum modelling of cardiac electrophysiology ('Cardiac Chaste'), discrete cell-based modelling of soft tissues ('Cell-based Chaste'), and modelling of ventilation in lungs ('Lung Chaste'). Cardiac Chaste addresses the need for a high-performance, generic, and verified simulation framework for cardiac electrophysiology that is freely available to the scientific community. Cardiac chaste provides a software package capable of realistic heart simulations that is efficient, rigorously tested, and runs on HPC platforms. Cell-based Chaste addresses the need for efficient and verified implementations of cell-based modelling frameworks, providing a set of extensible tools for simulating biological tissues. Computational modelling, along with live imaging techniques, plays an important role in understanding the processes of tissue growth and repair. A wide range of cell-based modelling frameworks have been developed that have each been successfully applied in a range of biological applications. Cell-based Chaste includes implementations of the cellular automaton model, the cellular Potts model, cell-centre models with cell representations as overlapping spheres or Voronoi tessellations, and the vertex model. Lung Chaste addresses the need for a novel, generic and efficient lung modelling software package that is both tested and verified. It aims to couple biophysically-detailed models of airway mechanics with organ-scale ventilation models in a package that is freely available to the scientific community. Chaste is designed to be modular and extensible, providing libraries for common scientific computing infrastructure such as linear algebra operations, finite element meshes, and ordinary and partial differential equation solvers. This infrastructure is used by libraries for specific applications, such as continuum mechanics, cardiac models, and cell-based models. The software engineering techniques used to develop Chaste are intended to ensure code quality, re-usability and reliability. Primary applications of the software include cardiac and respiratory physiology, cancer and developmental biology.

Addressing Interdisciplinary Difficulties in Developmental Biology/Mathematical Collaborations: A Neural Crest Example.

Mathematical modeling can allow insight into the biological processes that can be difficult to access by conventional biological means alone. Such projects are becoming increasingly attractive with the appearance of faster and more powerful quantitative techniques in both biological data acquisition and data storage, manipulation, and presentation. However, as is frequent in interdisciplinary research, the main hurdles are not within the mindset and techniques of each discipline but are usually encountered in attempting to meld the different disciplines together. Based upon our experience in applying mathematical methods to investigate how neural crest cells interact to form the enteric nervous system, we present our views on how to pursue biomathematical modeling projects, what difficulties to expect, and how to overcome, or at least survive, these hurdles. The main advice being: persevere.

Cooperation, competition and antibiotic resistance in bacterial colonies.

Bacteria commonly live in dense and genetically diverse communities associated with surfaces. In these communities, competition for resources and space is intense, and yet we understand little of how this affects the spread of antibiotic-resistant strains. Here, we study interactions between antibiotic-resistant and susceptible strains using in vitro competition experiments in the opportunistic pathogen Pseudomonas aeruginosa and in silico simulations. Selection for intracellular resistance to streptomycin is very strong in colonies, such that resistance is favoured at very low antibiotic doses. In contrast, selection for extracellular resistance to carbenicillin is weak in colonies, and high doses of antibiotic are required to select for resistance. Manipulating the density and spatial structure of colonies reveals that this difference is partly explained by the fact that the local degradation of carbenicillin by ß-lactamase-secreting cells protects neighbouring sensitive cells from carbenicillin. In addition, we discover a second unexpected effect: the inducible elongation of cells in response to carbenicillin allows sensitive cells to better compete for the rapidly growing colony edge. These combined effects mean that antibiotic treatment can select against antibiotic-resistant strains, raising the possibility of treatment regimes that suppress sensitive strains while limiting the rise of antibiotic resistance. We argue that the detailed study of bacterial interactions will be fundamental to understanding and overcoming antibiotic resistance.