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.

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.

MOTIVATION: Live-cell microscopy has become an essential tool for analyzing dynamic processes in various biological applications. Thereby, high-throughput and automated tracking analyses allow the simultaneous evaluation of large numbers of objects. However, to critically assess the influence of individual objects on calculated summary statistics, and to detect heterogeneous dynamics or possible artifacts, such as misclassified or -tracked objects, a direct mapping of gained statistical information onto the actual image data would be necessary. RESULTS: We present VisuStatR as a platform independent software package that allows the direct visualization of time-resolved summary statistics of morphological characteristics or motility dynamics onto raw images. The software contains several display modes to compare user-defined summary statistics and the underlying image data in various levels of detail. AVAILABILITY AND IMPLEMENTATION: VisuStatR is a free and open-source R-package, containing a user-friendly graphical-user interface and is available via GitHub at https://github.com/grrchrr/VisuStatR/ under the MIT+ license. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Mathematical models in epidemiology are an indispensable tool to determine the dynamics and important characteristics of infectious diseases. Apart from their scientific merit, these models are often used to inform political decisions and interventional measures during an ongoing outbreak. However, reliably inferring the epidemical dynamics by connecting complex models to real data is still hard and requires either laborious manual parameter fitting or expensive optimization methods which have to be repeated from scratch for every application of a given model. In this work, we address this problem with a novel combination of epidemiological modeling with specialized neural networks. Our approach entails two computational phases: In an initial training phase, a mathematical model describing the epidemic is used as a coach for a neural network, which acquires global knowledge about the full range of possible disease dynamics. In the subsequent inference phase, the trained neural network processes the observed data of an actual outbreak and infers the parameters of the model in order to realistically reproduce the observed dynamics and reliably predict future progression. With its flexible framework, our simulation-based approach is applicable to a variety of epidemiological models. Moreover, since our method is fully Bayesian, it is designed to incorporate all available prior knowledge about plausible parameter values and returns complete joint posterior distributions over these parameters. Application of our method to the early Covid-19 outbreak phase in Germany demonstrates that we are able to obtain reliable probabilistic estimates for important disease characteristics, such as generation time, fraction of undetected infections, likelihood of transmission before symptom onset, and reporting delays using a very moderate amount of real-world observations.

Most biological systems are difficult to analyse due to a multitude of interacting components and the concomitant lack of information about the essential dynamics. Finding appropriate models that provide a systematic description of such biological systems and that help to identify their relevant factors and processes can be challenging given the sheer number of possibilities. Model selection algorithms that evaluate the performance of a multitude of different models against experimental data provide a useful tool to identify appropriate model structures. However, many algorithms addressing the analysis of complex dynamical systems, as they are often used in biology, compare a preselected number of models or rely on exhaustive searches of the total model space which might be unfeasible dependent on the number of possibilities. Therefore, we developed an algorithm that is able to perform model selection on complex systems and searches large model spaces in a dynamical way. Our algorithm includes local and newly developed non-local search methods that can prevent the algorithm from ending up in local minima of the model space by accounting for structurally similar processes. We tested and validated the algorithm based on simulated data and showed its flexibility for handling different model structures. We also used the algorithm to analyse experimental data on the cell proliferation dynamics of CD4+ and CD8+ T cells that were cultured under different conditions. Our analyses indicated dynamical changes within the proliferation potential of cells that was reduced within tissue-like 3D ex vivo cultures compared to suspension. Due to the flexibility in handling various model structures, the algorithm is applicable to a large variety of different biological problems and represents a useful tool for the data-oriented evaluation of complex model spaces.

Cytotoxic T lymphocytes (CTLs) are important agents in the control of intracellular pathogens, which specifically recognize and kill infected cells. Recently developed experimental methods allow the estimation of the CTL's efficacy in detecting and clearing infected host cells. One method, the in vivo killing assay, utilizes the adoptive transfer of antigen displaying target cells into the bloodstream of mice. Surprisingly, killing efficacies measured by this method are often much higher than estimates obtained by other methods based on, for instance, the dynamics of escape mutations. In this study, we investigated what fraction of this variation can be explained by differences in peptide loads employed in in vivo killing assays. We addressed this question in mice immunized with lymphocytic choriomeningitis virus (LCMV). We conducted in vivo killing assays varying the loads of the immunodominant epitope GP33 on target cells. Using a mathematical model, we determined the efficacy of effector and memory CTL, as well as CTL in chronically infected mice. We found that the killing efficacy is substantially reduced at lower peptide loads. For physiological peptide loads, our analysis predicts more than a factor 10 lower CTL efficacies than at maximum peptide loads. Assuming that the efficacy scales linearly with the frequency of CTL, a clear hierarchy emerges among the groups across all peptide antigen concentrations. The group of mice with chronic LCMV infections shows a consistently higher killing efficacy per CTL than the acutely infected mouse group, which in turn has a consistently larger efficacy than the memory mouse group. We conclude that CTL killing efficacy dependence on surface epitope frequencies can only partially explain the variation in in vivo killing efficacy estimates across experimental methods and viral systems, which vary about four orders of magnitude. In contrast, peptide load differences can explain at most two orders of magnitude.

Chronic liver infection by hepatitis C virus (HCV) is a major public health concern. Despite partly successful treatment options, several aspects of intrahepatic HCV infection dynamics are still poorly understood, including the preferred mode of viral propagation, as well as the proportion of infected hepatocytes. Answers to these questions have important implications for the development of therapeutic interventions. In this study, we present methods to analyze the spatial distribution of infected hepatocytes obtained by single cell laser capture microdissection from liver biopsy samples of patients chronically infected with HCV. By characterizing the internal structure of clusters of infected cells, we are able to evaluate hypotheses about intrahepatic infection dynamics. We found that individual clusters on biopsy samples range in size from 4-50 infected cells. In addition, the HCV RNA content in a cluster declines from the cell that presumably founded the cluster to cells at the maximal cluster extension. These observations support the idea that HCV infection in the liver is seeded randomly (e.g. from the blood) and then spreads locally. Assuming that the amount of intracellular HCV RNA is a proxy for how long a cell has been infected, we estimate based on models of intracellular HCV RNA replication and accumulation that cells in clusters have been infected on average for less than a week. Further, we do not find a relationship between the cluster size and the estimated cluster expansion time. Our method represents a novel approach to make inferences about infection dynamics in solid tissues from static spatial data.

Secondary lymphoid organs (SLO), such as lymph nodes and the spleen, display a complex micro-architecture. In the T cell zone the micro-architecture is provided by a network of fibroblastic reticular cells (FRC) and their filaments. The FRC network is thought to enhance the interaction between immune cells and their cognate antigen. However, the effect of the FRC network on cell interaction cannot be quantified to date because of limitations in immunological methodology. We use computational models to study the influence of different densities of FRC networks on the probability that two cells meet. We developed a 3D cellular automaton model to simulate cell movements and interactions along the FRC network inside lymphatic tissue. We show that the FRC network density has only a small effect on the probability of a cell to come into contact with a static or motile target. However, damage caused by a disruption of the FRC network is greatest at FRC densities corresponding to densities observed in the spleen of naïve mice. Our analysis suggests that the FRC network as a guiding structure for moving T cells has only a minor effect on the probability to find a corresponding dendritic cell. We propose alternative hypotheses by which the FRC network might influence the functionality of immune responses in a more significant way.

Cytotoxic T lymphocytes (CTLs) are important immune effectors against intra-cellular pathogens. These cells search for infected cells and kill them. Recently developed experimental methods in combination with mathematical models allow for the quantification of the efficacy of CTL killing in vivo and, hence, for the estimation of parameters that characterize the effect of CTL killing on the target cell populations. It is not known how these population-level parameters relate to single-cell properties. To address this question, we developed a three-dimensional cellular automaton model of the region of the spleen where CTL killing takes place. The cellular automaton model describes the movement of different cell populations and their interactions. Cell movement patterns in our cellular automaton model agree with observations from two-photon microscopy. We find that, despite the strong spatial nature of the kinetics in our cellular automaton model, the killing of target cells by CTLs can be described by a term which is linear in the target cell frequency and saturates with respect to the CTL levels. Further, we find that the parameters describing CTL killing on the population level are most strongly impacted by the time a CTL needs to kill a target cell. This suggests that the killing of target cells, rather than their localization, is the limiting step in CTL killing dynamics given reasonable frequencies of CTL. Our analysis identifies additional experimental directions which are of particular importance to interpret estimates of killing rates and could advance our quantitative understanding of CTL killing.