Selection, calibration, and validation of models of tumor growth.

This paper presents general approaches for addressing some of the most important issues in predictive computational oncology concerned with developing classes of predictive models of tumor growth. First, the process of developing mathematical models of vascular tumors evolving in the complex, heterogeneous, macroenvironment of living tissue; second, the selection of the most plausible models among these classes, given relevant observational data; third, the statistical calibration and validation of models in these classes, and finally, the prediction of key Quantities of Interest (QOIs) relevant to patient survival and the effect of various therapies. The most challenging aspects of this endeavor is that all of these issues often involve confounding uncertainties: in observational data, in model parameters, in model selection, and in the features targeted in the prediction. Our approach can be referred to as "model agnostic" in that no single model is advocated; rather, a general approach that explores powerful mixture-theory representations of tissue behavior while accounting for a range of relevant biological factors is presented, which leads to many potentially predictive models. Then representative classes are identified which provide a starting point for the implementation of OPAL, the Occam Plausibility Algorithm (OPAL) which enables the modeler to select the most plausible models (for given data) and to determine if the model is a valid tool for predicting tumor growth and morphology (in vivo). All of these approaches account for uncertainties in the model, the observational data, the model parameters, and the target QOI. We demonstrate these processes by comparing a list of models for tumor growth, including reaction-diffusion models, phase-fields models, and models with and without mechanical deformation effects, for glioma growth measured in murine experiments. Examples are provided that exhibit quite acceptable predictions of tumor growth in laboratory animals while demonstrating successful implementations of OPAL.

[1]  J. Tinsley Oden,et al.  Bayesian calibration, validation, and uncertainty quantification of diffuse interface models of tumor growth , 2012, Journal of Mathematical Biology.

[2]  G. Box Robustness in the Strategy of Scientific Model Building. , 1979 .

[3]  D. Hanahan,et al.  Hallmarks of Cancer: The Next Generation , 2011, Cell.

[4]  Thomas E Yankeelov,et al.  Predicting the Response of Breast Cancer to Neoadjuvant Therapy Using a Mechanically Coupled Reaction-Diffusion Model. , 2015, Cancer research.

[5]  Ivo Babuška,et al.  Predictive Computational Science: Computer Predictions in the Presence of Uncertainty , 2017 .

[6]  Yusheng Feng,et al.  Toward Predictive Multiscale Modeling of Vascular Tumor Growth , 2015, Archives of Computational Methods in Engineering.

[7]  J. Tinsley Oden,et al.  Analysis and numerical solution of stochastic phase‐field models of tumor growth , 2015 .

[8]  G. Chow A comparison of the information and posterior probability criteria for model selection , 1981 .

[9]  Benjamin S. Kirk,et al.  Library for Parallel Adaptive Mesh Refinement / Coarsening Simulations , 2006 .

[10]  J. Tinsley Oden,et al.  GENERAL DIFFUSE-INTERFACE THEORIES AND AN APPROACH TO PREDICTIVE TUMOR GROWTH MODELING , 2010 .

[11]  Omar M. Knio,et al.  Spectral Methods for Uncertainty Quantification , 2010 .

[12]  J. Tinsley Oden,et al.  Selection, calibration, and validation of coarse-grained models of atomistic systems , 2015 .

[13]  Mario Bunge,et al.  Classical Field Theories , 1967 .

[14]  Sai Hung Cheung,et al.  PARALLEL ADAPTIVE MULTILEVEL SAMPLING ALGORITHMS FOR THE BAYESIAN ANALYSIS OF MATHEMATICAL MODELS , 2012 .

[15]  Nicola Bellomo,et al.  Toward a mathematical theory of Keller–Segel models of pattern formation in biological tissues , 2015 .

[16]  D. Hanahan,et al.  The Hallmarks of Cancer , 2000, Cell.

[17]  Vittorio Cristini,et al.  Multiscale Modeling of Cancer: An Integrated Experimental and Mathematical Modeling Approach , 2010 .

[18]  G. Kitagawa,et al.  Information Criteria and Statistical Modeling , 2007 .

[19]  Thomas E Yankeelov,et al.  Clinically Relevant Modeling of Tumor Growth and Treatment Response , 2013, Science Translational Medicine.

[20]  R. Jain,et al.  Micro-Environmental Mechanical Stress Controls Tumor Spheroid Size and Morphology by Suppressing Proliferation and Inducing Apoptosis in Cancer Cells , 2009, PloS one.

[21]  Karl W. Schulz,et al.  The Parallel C++ Statistical Library 'QUESO': Quantification of Uncertainty for Estimation, Simulation and Optimization , 2011, Euro-Par Workshops.

[22]  Clifford Ambrose Truesdell,et al.  Mechanical Basis of Diffusion , 1962 .

[23]  Paolo A. Netti,et al.  Solid stress inhibits the growth of multicellular tumor spheroids , 1997, Nature Biotechnology.

[24]  Vittorio Cristini,et al.  Multiscale cancer modeling. , 2010, Annual review of biomedical engineering.

[25]  Nicola Bellomo,et al.  On the foundations of cancer modelling: Selected topics, speculations, and perspectives , 2008 .

[26]  J. Tinsley Oden,et al.  SELECTION AND ASSESSMENT OF PHENOMENOLOGICAL MODELS OF TUMOR GROWTH , 2013 .

[27]  A. Concha Theory of Mixtures , 2014 .

[28]  Thomas E Yankeelov,et al.  Toward a science of tumor forecasting for clinical oncology. , 2015, Cancer research.

[29]  James D. Murray,et al.  A generalized diffusion model for growth and dispersal in a population , 1981 .

[30]  Triantafyllos Stylianopoulos,et al.  Causes, consequences, and remedies for growth-induced solid stress in murine and human tumors , 2012, Proceedings of the National Academy of Sciences.

[31]  David A. Hormuth II,et al.  Predicting in vivo glioma growth with the reaction diffusion equation constrained by quantitative magnetic resonance imaging data , 2015, Physical biology.

[32]  J. Beck Bayesian system identification based on probability logic , 2010 .

[33]  Ernesto A. B. F. Lima,et al.  A hybrid ten-species phase-field model of tumor growth , 2014 .

[34]  L. M. M.-T. Theory of Probability , 1929, Nature.

[35]  Masahiko Morinaga,et al.  Phase-Field Simulation of the Effect of Elastic Inhomogeneity on Microstructure Evolution in Ni-Based Superalloys , 2009 .

[36]  John D. Ingram,et al.  A Continuum theory of chemically reacting media—I , 1965 .

[37]  Thomas E Yankeelov,et al.  A mechanically coupled reaction–diffusion model for predicting the response of breast tumors to neoadjuvant chemotherapy , 2013, Physics in medicine and biology.

[38]  H. Frieboes,et al.  Three-dimensional multispecies nonlinear tumor growth--I Model and numerical method. , 2008, Journal of theoretical biology.

[39]  R. M. Bowen Part I – Theory of Mixtures , 1976 .

[40]  Saltelli Andrea,et al.  Global Sensitivity Analysis: The Primer , 2008 .