Dataset Information

Practical parameter identifiability for spatio-temporal models of cell invasion.

ABSTRACT: We examine the practical identifiability of parameters in a spatio-temporal reaction-diffusion model of a scratch assay. Experimental data involve fluorescent cell cycle labels, providing spatial information about cell position and temporal information about the cell cycle phase. Cell cycle labelling is incorporated into the reaction-diffusion model by treating the total population as two interacting subpopulations. Practical identifiability is examined using a Bayesian Markov chain Monte Carlo (MCMC) framework, confirming that the parameters are identifiable when we assume the diffusivities of the subpopulations are identical, but that the parameters are practically non-identifiable when we allow the diffusivities to be distinct. We also assess practical identifiability using a profile likelihood approach, providing similar results to MCMC with the advantage of being an order of magnitude faster to compute. Therefore, we suggest that the profile likelihood ought to be adopted as a screening tool to assess practical identifiability before MCMC computations are performed.

SUBMITTER: Simpson MJ

PROVIDER: S-EPMC7115235 | biostudies-literature | 2020 Mar

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Practical parameter identifiability for spatio-temporal models of cell invasion.

Simpson Matthew J MJ Baker Ruth E RE Vittadello Sean T ST Maclaren Oliver J OJ

Journal of the Royal Society, Interface 20200304 164

We examine the practical identifiability of parameters in a spatio-temporal reaction-diffusion model of a scratch assay. Experimental data involve fluorescent cell cycle labels, providing spatial information about cell position and temporal information about the cell cycle phase. Cell cycle labelling is incorporated into the reaction-diffusion model by treating the total population as two interacting subpopulations. Practical identifiability is examined using a Bayesian Markov chain Monte Carlo ...[more]

PMID: 32126193

Similar Datasets

Project description:BackgroundHeidenreich et al. (Risk Anal 1997 17 391-399) considered parameter identifiability in the context of the two-mutation cancer model and demonstrated that combinations of all but two of the model parameters are identifiable. We consider the problem of identifiability in the recently developed carcinogenesis models of Little and Wright (Math Biosci 2003 183 111-134) and Little et al. (J Theoret Biol 2008 254 229-238). These models, which incorporate genomic instability, generalize a large number of other quasi-biological cancer models, in particular those of Armitage and Doll (Br J Cancer 1954 8 1-12), the two-mutation model (Moolgavkar et al. Math Biosci 1979 47 55-77), the generalized multistage model of Little (Biometrics 1995 51 1278-1291), and a recently developed cancer model of Nowak et al. (PNAS 2002 99 16226-16231).Methodology/principal findingsWe show that in the simpler model proposed by Little and Wright (Math Biosci 2003 183 111-134) the number of identifiable combinations of parameters is at most two less than the number of biological parameters, thereby generalizing previous results of Heidenreich et al. (Risk Anal 1997 17 391-399) for the two-mutation model. For the more general model of Little et al. (J Theoret Biol 2008 254 229-238) the number of identifiable combinations of parameters is at most less than the number of biological parameters, where is the number of destabilization types, thereby also generalizing all these results. Numerical evaluations suggest that these bounds are sharp. We also identify particular combinations of identifiable parameters.Conclusions/significanceWe have shown that the previous results on parameter identifiability can be generalized to much larger classes of quasi-biological carcinogenesis model, and also identify particular combinations of identifiable parameters. These results are of theoretical interest, but also of practical significance to anyone attempting to estimate parameters for this large class of cancer models.

Project description:BackgroundKinetic models of biochemical systems usually consist of ordinary differential equations that have many unknown parameters. Some of these parameters are often practically unidentifiable, that is, their values cannot be uniquely determined from the available data. Possible causes are lack of influence on the measured outputs, interdependence among parameters, and poor data quality. Uncorrelated parameters can be seen as the key tuning knobs of a predictive model. Therefore, before attempting to perform parameter estimation (model calibration) it is important to characterize the subset(s) of identifiable parameters and their interplay. Once this is achieved, it is still necessary to perform parameter estimation, which poses additional challenges.MethodsWe present a methodology that (i) detects high-order relationships among parameters, and (ii) visualizes the results to facilitate further analysis. We use a collinearity index to quantify the correlation between parameters in a group in a computationally efficient way. Then we apply integer optimization to find the largest groups of uncorrelated parameters. We also use the collinearity index to identify small groups of highly correlated parameters. The results files can be visualized using Cytoscape, showing the identifiable and non-identifiable groups of parameters together with the model structure in the same graph.ResultsOur contributions alleviate the difficulties that appear at different stages of the identifiability analysis and parameter estimation process. We show how to combine global optimization and regularization techniques for calibrating medium and large scale biological models with moderate computation times. Then we evaluate the practical identifiability of the estimated parameters using the proposed methodology. The identifiability analysis techniques are implemented as a MATLAB toolbox called VisId, which is freely available as open source from GitHub ( https://github.com/gabora/visid ).ConclusionsOur approach is geared towards scalability. It enables the practical identifiability analysis of dynamic models of large size, and accelerates their calibration. The visualization tool allows modellers to detect parts that are problematic and need refinement or reformulation, and provides experimentalists with information that can be helpful in the design of new experiments.

Project description:BackgroundDifferentiating infiltrating myeloid cells from resident microglia in neuroinflammatory disease is challenging, because bone marrow-derived inflammatory monocytes infiltrating the inflamed brain adopt a 'microglia-like' phenotype. This precludes the accurate identification of either cell type without genetic manipulation, which is important to understand their temporal contribution to disease and inform effective intervention in its pathogenesis. During West Nile virus (WNV) encephalitis, widespread neuronal infection drives substantial CNS infiltration of inflammatory monocytes, causing severe immunopathology and/or death, but the role of microglia in this remains unclear.MethodsUsing high-parameter cytometry and dimensionality-reduction, we devised a simple, novel gating strategy to identify microglia and infiltrating myeloid cells during WNV-infection. Validating our strategy, we (1) blocked the entry of infiltrating myeloid populations from peripheral blood using monoclonal blocking antibodies, (2) adoptively transferred BM-derived monocytes and tracked their phenotypic changes after infiltration and (3) labelled peripheral leukocytes that infiltrate into the brain with an intravenous dye. We demonstrated that myeloid immigrants populated only the identified macrophage gates, while PLX5622 depletion reduced all 4 subsets defined by the microglial gates.ResultsUsing this gating approach, we identified four consistent microglia subsets in the homeostatic and WNV-infected brain. These were P2RY12hi CD86-, P2RY12hi CD86+ and P2RY12lo CD86- P2RY12lo CD86+. During infection, 2 further populations were identified as 'inflammatory' and 'microglia-like' macrophages, recruited from the bone marrow. Detailed kinetic analysis showed significant increases in the proportions of both P2RY12lo microglia subsets in all anatomical areas, largely at the expense of the P2RY12hi CD86- subset, with the latter undergoing compensatory proliferation, suggesting replenishment of, and differentiation from this subset in response to infection. Microglia altered their morphology early in infection, with all cells adopting temporal and regional disease-specific phenotypes. Late in disease, microglia produced IL-12, downregulated CX3CR1, F4/80 and TMEM119 and underwent apoptosis. Infiltrating macrophages expressed both TMEM119 and P2RY12 de novo, with the microglia-like subset notably exhibiting the highest proportional myeloid population death.ConclusionsOur approach enables detailed kinetic analysis of resident vs infiltrating myeloid cells in a wide range of neuroinflammatory models without non-physiological manipulation. This will more clearly inform potential therapeutic approaches that specifically modulate these cells.

Dataset Information

Practical parameter identifiability for spatio-temporal models of cell invasion.

Publications

Practical parameter identifiability for spatio-temporal models of cell invasion.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets