FitMultiCell: simulating and parameterizing computational models of multi-scale and multi-cellular processes.

MOTIVATION: Biological tissues are dynamic and highly organized. Multi-scale models are helpful tools to analyse and understand the processes determining tissue dynamics. These models usually depend on parameters that need to be inferred from experimental data to achieve a quantitative understanding, to predict the response to perturbations, and to evaluate competing hypotheses. However, even advanced inference approaches such as approximate Bayesian computation (ABC) are difficult to apply due to the computational complexity of the simulation of multi-scale models. Thus, there is a need for a scalable pipeline for modeling, simulating, and parameterizing multi-scale models of multi-cellular processes. RESULTS: Here, we present FitMultiCell, a computationally efficient and user-friendly open-source pipeline that can handle the full workflow of modeling, simulating, and parameterizing for multi-scale models of multi-cellular processes. The pipeline is modular and integrates the modeling and simulation tool Morpheus and the statistical inference tool pyABC. The easy integration of high-performance infrastructure allows to scale to computationally expensive problems. The introduction of a novel standard for the formulation of parameter inference problems for multi-scale models additionally ensures reproducibility and reusability. By applying the pipeline to multiple biological problems, we demonstrate its broad applicability, which will benefit in particular image-based systems biology. AVAILABILITY AND IMPLEMENTATION: FitMultiCell is available open-source at https://gitlab.com/fitmulticell/fit.

Influence of cell type specific infectivity and tissue composition on SARS-CoV-2 infection dynamics within human airway epithelium.

Human airway epithelium (HAE) represents the primary site of viral infection for SARS-CoV-2. Comprising different cell populations, a lot of research has been aimed at deciphering the major cell types and infection dynamics that determine disease progression and severity. However, the cell type-specific replication kinetics, as well as the contribution of cellular composition of the respiratory epithelium to infection and pathology are still not fully understood. Although experimental advances, including Air-liquid interface (ALI) cultures of reconstituted pseudostratified HAE, as well as lung organoid systems, allow the observation of infection dynamics under physiological conditions in unprecedented level of detail, disentangling and quantifying the contribution of individual processes and cells to these dynamics remains challenging. Here, we present how a combination of experimental data and mathematical modelling can be used to infer and address the influence of cell type specific infectivity and tissue composition on SARS-CoV-2 infection dynamics. Using a stepwise approach that integrates various experimental data on HAE culture systems with regard to tissue differentiation and infection dynamics, we develop an individual cell-based model that enables investigation of infection and regeneration dynamics within pseudostratified HAE. In addition, we present a novel method to quantify tissue integrity based on image data related to the standard measures of transepithelial electrical resistance measurements. Our analysis provides a first aim of quantitatively assessing cell type specific infection kinetics and shows how tissue composition and changes in regeneration capacity, as e.g. in smokers, can influence disease progression and pathology. Furthermore, we identified key measurements that still need to be assessed in order to improve inference of cell type specific infection kinetics and disease progression. Our approach provides a method that, in combination with additional experimental data, can be used to disentangle the complex dynamics of viral infection and immunity within human airway epithelial culture systems.

Guidance by followers ensures long-range coordination of cell migration through α-catenin mechanoperception.

Morphogenesis, wound healing, and some cancer metastases depend upon the migration of cell collectives that need to be guided to their destination as well as coordinated with other cell movements. During zebrafish gastrulation, the extension of the embryonic axis is led by the mesendodermal polster that migrates toward the animal pole, followed by the axial mesoderm that undergoes convergence and extension. Here, we investigate how polster cells are guided toward the animal pole. Using a combination of precise laser ablations, advanced transplants, and functional as well as in silico approaches, we establish that each polster cell is oriented by its immediate follower cells. Each cell perceives the migration of followers, through E-cadherin/α-catenin mechanotransduction, and aligns with them. Therefore, directional information propagates from cell to cell over the whole tissue. Such guidance of migrating cells by followers ensures long-range coordination of movements and developmental robustness.

Cell-cell adhesion and 3D matrix confinement determine jamming transitions in breast cancer invasion.

Plasticity of cancer invasion and metastasis depends on the ability of cancer cells to switch between collective and single-cell dissemination, controlled by cadherin-mediated cell-cell junctions. In clinical samples, E-cadherin-expressing and -deficient tumours both invade collectively and metastasize equally, implicating additional mechanisms controlling cell-cell cooperation and individualization. Here, using spatially defined organotypic culture, intravital microscopy of mammary tumours in mice and in silico modelling, we identify cell density regulation by three-dimensional tissue boundaries to physically control collective movement irrespective of the composition and stability of cell-cell junctions. Deregulation of adherens junctions by downregulation of E-cadherin and p120-catenin resulted in a transition from coordinated to uncoordinated collective movement along extracellular boundaries, whereas single-cell escape depended on locally free tissue space. These results indicate that cadherins and extracellular matrix confinement cooperate to determine unjamming transitions and stepwise epithelial fluidization towards, ultimately, cell individualization.

Morpheus: a user-friendly modeling environment for multiscale and multicellular systems biology.

Morpheus is a modeling environment for the simulation and integration of cell-based models with ordinary differential equations and reaction-diffusion systems. It allows rapid development of multiscale models in biological terms and mathematical expressions rather than programming code. Its graphical user interface supports the entire workflow from model construction and simulation to visualization, archiving and batch processing.

Collective motion and nonequilibrium cluster formation in colonies of gliding bacteria.

We characterize cell motion in experiments and show that the transition to collective motion in colonies of gliding bacterial cells confined to a monolayer appears through the organization of cells into larger moving clusters. Collective motion by nonequilibrium cluster formation is detected for a critical cell packing fraction around 17%. This transition is characterized by a scale-free power-law cluster-size distribution, with an exponent 0.88±0.07, and the appearance of giant number fluctuations. Our findings are in quantitative agreement with simulations of self-propelled rods. This suggests that the interplay of self-propulsion and the rod shape of bacteria is sufficient to induce collective motion.

Pattern-formation mechanisms in motility mutants of Myxococcus xanthus.

Formation of spatial patterns of cells is a recurring theme in biology and often depends on regulated cell motility. Motility of the rod-shaped cells of the bacterium Myxococcus xanthus depends on two motility machineries, type IV pili (giving rise to S-motility) and the gliding motility apparatus (giving rise to A-motility). Cell motility is regulated by occasional reversals. Moving M. xanthus cells can organize into spreading colonies or spore-filled fruiting bodies, depending on their nutritional status. To ultimately understand these two pattern-formation processes and the contributions by the two motility machineries, as well as the cell reversal machinery, we analyse spatial self-organization in three M. xanthus strains: (i) a mutant that moves unidirectionally without reversing by the A-motility system only, (ii) a unidirectional mutant that is also equipped with the S-motility system, and (iii) the wild-type that, in addition to the two motility systems, occasionally reverses its direction of movement. The mutant moving by means of the A-engine illustrates that collective motion in the form of large moving clusters can arise in gliding bacteria owing to steric interactions of the rod-shaped cells, without the need of invoking any biochemical signal regulation. The two-engine strain mutant reveals that the same phenomenon emerges when both motility systems are present, and as long as cells exhibit unidirectional motion only. From the study of these two strains, we conclude that unidirectional cell motion induces the formation of large moving clusters at low and intermediate densities, while it results in vortex formation at very high densities. These findings are consistent with what is known from self-propelled rod models, which strongly suggests that the combined effect of self-propulsion and volume exclusion interactions is the pattern-formation mechanism leading to the observed phenomena. On the other hand, we learn that when cells occasionally reverse their moving direction, as observed in the wild-type, cells form small but strongly elongated clusters and self-organize into a mesh-like structure at high enough densities. These results have been obtained from a careful analysis of the cluster statistics of ensembles of cells, and analysed in the light of a coagulation Smoluchowski equation with fragmentation.

Early embryonic vascular patterning by matrix-mediated paracrine signalling: a mathematical model study.

During embryonic vasculogenesis, endothelial precursor cells of mesodermal origin known as angioblasts assemble into a characteristic network pattern. Although a considerable amount of markers and signals involved in this process have been identified, the mechanisms underlying the coalescence of angioblasts into this reticular pattern remain unclear. Various recent studies hypothesize that autocrine regulation of the chemoattractant vascular endothelial growth factor (VEGF) is responsible for the formation of vascular networks in vitro. However, the autocrine regulation hypothesis does not fit well with reported data on in vivo early vascular development. In this study, we propose a mathematical model based on the alternative assumption that endodermal VEGF signalling activity, having a paracrine effect on adjacent angioblasts, is mediated by its binding to the extracellular matrix (ECM). Detailed morphometric analysis of simulated networks and images obtained from in vivo quail embryos reveals the model mimics the vascular patterns with high accuracy. These results show that paracrine signalling can result in the formation of fine-grained cellular networks when mediated by angioblast-produced ECM. This lends additional support to the theory that patterning during early vascular development in the vertebrate embryo is regulated by paracrine signalling.

Modeling of pattern development during fibronectin nanofibril formation.

Formation of fibrillar fibronectin networks is a major process during embryogenesis and tissue formation, but the molecular details of fibril assembly remain poorly understood. Based on current ideas of fibronectin fibrillogenesis, a stochastic model was developed to enlighten the mechanism of the formation of paired fibronectin nanofibrils by adherent endothelial cells, which has been observed recently. The development of fibronectin clusters and fibrils was investigated by means of Monte Carlo simulations, including diffusion-controlled aggregation and myosin-driven transport of fibronectin-integrin complexes in the vicinity of a focal adhesion. Different evolving growth patterns were summarized in a morphological diagram as a function of the fibronectin substrate and fibronectin-fibronectin interaction energies. The formation of paired nanofibrils was found to occur in a certain region of interaction parameters. Beyond this region branched fibronectin clusters as well as tear-off of fibronectin fibrils were observed.