Calibration of Multi-Parameter Models of Avascular Tumor Growth Using Time Resolved Microscopy Data

Two of the central challenges of using mathematical models for predicting the spatiotemporal development of tumors is the lack of appropriate data to calibrate the parameters of the model, and quantitative characterization of the uncertainties in both the experimental data and the modeling process itself. We present a sequence of experiments, with increasing complexity, designed to systematically calibrate the rates of apoptosis, proliferation, and necrosis, as well as mobility, within a phase-field tumor growth model. The in vitro experiments characterize the proliferation and death of human liver carcinoma cells under different initial cell concentrations, nutrient availabilities, and treatment conditions. A Bayesian framework is employed to quantify the uncertainties in model parameters. The average difference between the calibration and the data, across all time points is between 11.54% and 14.04% for the apoptosis experiments, 7.33% and 23.30% for the proliferation experiments, and 8.12% and 31.55% for the necrosis experiments. The results indicate the proposed experiment-computational approach is generalizable and appropriate for step-by-step calibration of multi-parameter models, yielding accurate estimations of model parameters related to rates of proliferation, apoptosis, and necrosis.

[1]  Kristoffer G. van der Zee,et al.  Numerical simulation of a thermodynamically consistent four‐species tumor growth model , 2012, International journal for numerical methods in biomedical engineering.

[2]  Thomas E. Yankeelov,et al.  Multi-scale Modeling in Clinical Oncology: Opportunities and Barriers to Success , 2016, Annals of Biomedical Engineering.

[3]  Burton Ac,et al.  Rate of growth of solid tumours as a problem of diffusion. , 1966, Growth.

[4]  Elena Faggiano,et al.  A personalized mathematical tool for neuro-oncology: A clinical case study , 2018, International Journal of Non-Linear Mechanics.

[5]  Ingo Roeder,et al.  Metabolism is the tie: The Bertalanffy-type cancer growth model as common denominator of various modelling approaches , 2018, Biosyst..

[6]  A C Burton,et al.  Rate of growth of solid tumours as a problem of diffusion. , 1966, Growth.

[7]  Harald Garcke,et al.  A Cahn–Hilliard–Darcy model for tumour growth with chemotaxis and active transport , 2015, 1508.00437.

[8]  J. Oden,et al.  Selection and Validation of Predictive Models of Radiation Effects on Tumor Growth Based on Noninvasive Imaging Data. , 2017, Computer methods in applied mechanics and engineering.

[9]  Wei Li,et al.  Effects of CD73 on human colorectal cancer cell growth in vivo and in vitro. , 2016, Oncology reports.

[10]  Elizabeth E. Antoine,et al.  Tunable Collagen I Hydrogels for Engineered Physiological Tissue Micro-Environments , 2015, PloS one.

[11]  H. Lee,et al.  Mitomycin C induces apoptosis in cultured corneal fibroblasts derived from type II granular corneal dystrophy corneas , 2008, Molecular vision.

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

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

[14]  A. Brú,et al.  The universal dynamics of tumor growth. , 2003, Biophysical journal.

[15]  T E Yankeelov,et al.  Selection, calibration, and validation of models of tumor growth. , 2016, Mathematical models & methods in applied sciences : M3AS.

[16]  Thomas E Yankeelov,et al.  A Predictive Mathematical Modeling Approach for the Study of Doxorubicin Treatment in Triple Negative Breast Cancer , 2017, Scientific Reports.

[17]  Matthew R. DeWitt,et al.  Vascularized microfluidic platforms to mimic the tumor microenvironment , 2018, Biotechnology and bioengineering.

[18]  Amir Shahmoradi,et al.  Multilevel Bayesian Parameter Estimation in the Presence of Model Inadequacy and Data Uncertainty , 2017, 1711.10599.

[19]  Thomas J. R. Hughes,et al.  Hierarchically refined and coarsened splines for moving interface problems, with particular application to phase-field models of prostate tumor growth , 2017 .

[20]  Kwou-Yeung Wu,et al.  Mechanism of mitomycin-induced apoptosis in cultured corneal endothelial cells , 2008, Molecular vision.

[21]  J. Crowston,et al.  T lymphocyte mediated lysis of mitomycin C treated Tenon’s capsule fibroblasts , 2004, British Journal of Ophthalmology.

[22]  Martin A. Nowak,et al.  A spatial model predicts that dispersal and cell turnover limit intratumour heterogeneity , 2015, Nature.

[23]  Yusheng Feng,et al.  A fully coupled space-time multiscale modeling framework for predicting tumor growth. , 2017, Computer methods in applied mechanics and engineering.

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

[25]  V. Quaranta,et al.  Fractional Proliferation: A method to deconvolve cell population dynamics from single-cell data , 2012, Nature Methods.

[26]  Thomas E. Yankeelov,et al.  Variable Cell Line Pharmacokinetics Contribute to Non-Linear Treatment Response in Heterogeneous Cell Populations , 2018, Annals of Biomedical Engineering.

[27]  Rong Li,et al.  Tumor necrosis factor-α attenuates starvation-induced apoptosis through upregulation of ferritin heavy chain in hepatocellular carcinoma cells , 2013, BMC Cancer.

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

[29]  H. Byrne.  Dissecting cancer through mathematics: from the cell to the animal model , 2010, Nature Reviews Cancer.

[30]  R. Kamm,et al.  Three-dimensional microfluidic model for tumor cell intravasation and endothelial barrier function , 2012, Proceedings of the National Academy of Sciences.

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

[32]  Yolanda Vasquez,et al.  A three-dimensional in vitro model to demonstrate the haptotactic effect of monocyte chemoattractant protein-1 on atherosclerosis-associated monocyte migration. , 2017, International journal of biological macromolecules.

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

[34]  Joseph W Freeman,et al.  3D in vitro bioengineered tumors based on collagen I hydrogels. , 2011, Biomaterials.

[35]  L. H. Gray,et al.  The Histological Structure of Some Human Lung Cancers and the Possible Implications for Radiotherapy , 1955, British Journal of Cancer.

[36]  T. Yankeelov,et al.  In vitro vascularized liver and tumor tissue microenvironments on a chip for dynamic determination of nanoparticle transport and toxicity , 2018, Biotechnology and bioengineering.

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

[38]  A. C. M. Resende,et al.  A HYBRID THREE-SCALE MODEL OF TUMOR GROWTH. , 2018, Mathematical models & methods in applied sciences : M3AS.

[39]  J. Tinsley Oden,et al.  Adaptive multiscale predictive modelling , 2018, Acta Numerica.

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

[41]  Coryn A. L. Bailer-Jones,et al.  Practical Bayesian Inference: A Primer for Physical Scientists , 2017 .

[42]  J. Rathmell,et al.  Cell metabolism: an essential link between cell growth and apoptosis. , 2011, Biochimica et biophysica acta.

[43]  Konstantinos Konstantopoulos,et al.  Cancer cell motility: lessons from migration in confined spaces , 2016, Nature Reviews Cancer.

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

[45]  G. Lorusso,et al.  The matricellular protein CYR61 promotes breast cancer lung metastasis by facilitating tumor cell extravasation and suppressing anoikis , 2016, Oncotarget.

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

[47]  K S Iyer,et al.  A stochastic model for the growth of cells in cancer. , 1970, Biometrics.

[48]  Harald Garcke,et al.  A coupled surface-Cahn--Hilliard bulk-diffusion system modeling lipid raft formation in cell membranes , 2015, 1509.03655.

[49]  D. J. Eyre.  Unconditionally Gradient Stable Time Marching the Cahn-Hilliard Equation , 1998 .

[50]  Xiangrong Li,et al.  Nonlinear simulations of solid tumor growth using a mixture model: invasion and branching , 2009, Journal of mathematical biology.

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

[52]  P. Vlachos,et al.  Three-dimensional microfluidic collagen hydrogels for investigating flow-mediated tumor-endothelial signaling and vascular organization. , 2014, Tissue engineering. Part C, Methods.

[53]  De-guang Sun,et al.  Activated δ‑opioid receptors inhibit hydrogen peroxide‑induced apoptosis in liver cancer cells through the PKC/ERK signaling pathway. , 2014, Molecular medicine reports.

[54]  David A Hormuth,et al.  Incorporating drug delivery into an imaging-driven, mechanics-coupled reaction diffusion model for predicting the response of breast cancer to neoadjuvant chemotherapy: theory and preliminary clinical results , 2018, Physics in medicine and biology.