Toward Predictive Multiscale Modeling of Vascular Tumor Growth

New directions in medical and biomedical sciences have gradually emerged over recent years that will change the way diseases are diagnosed and treated and are leading to the redirection of medicine toward patient-specific treatments. We refer to these new approaches for studying biomedical systems as predictive medicine, a new version of medical science that involves the use of advanced computer models of biomedical phenomena, high-performance computing, new experimental methods for model data calibration, modern imaging technologies, cutting-edge numerical algorithms for treating large stochastic systems, modern methods for model selection, calibration, validation, verification, and uncertainty quantification, and new approaches for drug design and delivery, all based on predictive models. The methodologies are designed to study events at multiple scales, from genetic data, to sub-cellular signaling mechanisms, to cell interactions, to tissue physics and chemistry, to organs in living human subjects. The present document surveys work on the development and implementation of predictive models of vascular tumor growth, covering aspects of what might be called modeling-and-experimentally based computational oncology. The work described is that of a multi-institutional team, centered at ICES with strong participation by members at M. D. Anderson Cancer Center and University of Texas at San Antonio. This exposition covers topics on signaling models, cell and cell-interaction models, tissue models based on multi-species mixture theories, models of angiogenesis, and beginning work of drug effects. A number of new parallel computer codes for implementing finite-element methods, multi-level Markov Chain Monte Carlo sampling methods, data classification methods, stochastic PDE solvers, statistical inverse algorithms for model calibration and validation, models of events at different spatial and temporal scales is presented. Importantly, new methods for model selection in the presence of uncertainties fundamental to predictive medical science, are described which are based on the notion of Bayesian model plausibilities. Also, as part of this general approach, new codes for determining the sensitivity of model outputs to variations in model parameters are described that provide a basis for assessing the importance of model parameters and controlling and reducing the number of relevant model parameters. Model specific data is to be accessible through careful and model-specific platforms in the Tumor Engineering Laboratory. We describe parallel computer platforms on which large-scale calculations are run as well as specific time-marching algorithms needed to treat stiff systems encountered in some phase-field mixture models. We also cover new non-invasive imaging and data classification methods that provide in vivo data for model validation. The study concludes with a brief discussion of future work and open challenges.

[1]  Brian B. Avants,et al.  Symmetric diffeomorphic image registration with cross-correlation: Evaluating automated labeling of elderly and neurodegenerative brain , 2008, Medical Image Anal..

[2]  C. Fortuin,et al.  Study of the sensitivity of coupled reaction systems to uncertainties in rate coefficients. I Theory , 1973 .

[3]  M. Swartz,et al.  Lymphatic and interstitial flow in the tumour microenvironment: linking mechanobiology with immunity , 2012, Nature Reviews Cancer.

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

[5]  V. Sakkalis,et al.  A complete mathematical study of a 3D model of heterogeneous and anisotropic glioma evolution , 2009, 2009 Annual International Conference of the IEEE Engineering in Medicine and Biology Society.

[6]  Min Wu,et al.  The effect of interstitial pressure on tumor growth: coupling with the blood and lymphatic vascular systems. , 2013, Journal of theoretical biology.

[7]  Dirk Drasdo,et al.  Individual-based and continuum models of growing cell populations: a comparison , 2009, Journal of mathematical biology.

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

[9]  Rui-Sheng Wang,et al.  Boolean modeling in systems biology: an overview of methodology and applications , 2012, Physical biology.

[10]  Luigi Preziosi,et al.  An agent-based model for elasto-plastic mechanical interactions between cells, basement membrane and extracellular matrix. , 2012, Mathematical biosciences and engineering : MBE.

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

[12]  Alain Trouvé,et al.  Computing Large Deformation Metric Mappings via Geodesic Flows of Diffeomorphisms , 2005, International Journal of Computer Vision.

[13]  Alfons G. Hoekstra,et al.  Cellular automata models of tumour natural shrinkage , 2011 .

[14]  Kevin Burrage,et al.  Stochastic simulation in systems biology , 2014, Computational and structural biotechnology journal.

[15]  T. Seyfried,et al.  Metabolic therapy: a new paradigm for managing malignant brain cancer. , 2015, Cancer letters.

[16]  J. Berger Statistical Decision Theory and Bayesian Analysis , 1988 .

[17]  D. Xiu Numerical Methods for Stochastic Computations: A Spectral Method Approach , 2010 .

[18]  Paul Macklin,et al.  Improved patient-specific calibration for agent-based cancer modeling. , 2013, Journal of theoretical biology.

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

[20]  Marissa Nichole Rylander,et al.  Review of collagen I hydrogels for bioengineered tissue microenvironments: characterization of mechanics, structure, and transport. , 2014, Tissue engineering. Part B, Reviews.

[21]  Trachette L. Jackson,et al.  Modeling Tumor Vasculature: Molecular, Cellular, and Tissue Level Aspects and Implications , 2012 .

[22]  I. Shmulevich,et al.  Analysis of angiogenesis using in vitro experiments and stochastic growth models. , 2005, Physical review. E, Statistical, nonlinear, and soft matter physics.

[23]  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.

[24]  D. Wallace A Mitochondrial Paradigm of Metabolic and Degenerative Diseases, Aging, and Cancer: A Dawn for Evolutionary Medicine , 2005, Annual review of genetics.

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

[26]  Alissa M. Weaver,et al.  Mathematical modeling of cancer: the future of prognosis and treatment. , 2005, Clinica chimica acta; international journal of clinical chemistry.

[27]  Zhihui Wang,et al.  Cross-scale sensitivity analysis of a non-small cell lung cancer model: Linking molecular signaling properties to cellular behavior , 2007, Biosyst..

[28]  T. Ricken,et al.  A triphasic model of transversely isotropic biological tissue with applications to stress and biologically induced growth , 2007 .

[29]  Denis Wirtz,et al.  The physics of cancer: the role of physical interactions and mechanical forces in metastasis , 2011, Nature Reviews Cancer.

[30]  Ronald E. Miller,et al.  Atomistic/continuum coupling in computational materials science , 2003 .

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

[32]  J. Tinsley Oden,et al.  Solution verification, goal-oriented adaptive methods for stochastic advection-diffusion problems , 2010 .

[33]  W. Y. Tan,et al.  Stochastic modeling of carcinogenesis: Some new insights , 1998 .

[34]  G. Tierra,et al.  Numerical Methods for Solving the Cahn–Hilliard Equation and Its Applicability to Related Energy-Based Models , 2015 .

[35]  William A. Curtin,et al.  Multiscale plasticity modeling: coupled atomistics and discrete dislocation mechanics , 2004 .

[36]  Ina Koch,et al.  Genome Scale Modeling in Systems Biology: Algorithms and Resources , 2014, Current genomics.

[37]  M. Iruela-Arispe,et al.  Extracellular matrix, inflammation, and the angiogenic response. , 2010, Cardiovascular research.

[38]  M. Chaplain Avascular growth, angiogenesis and vascular growth in solid tumours: The mathematical modelling of the stages of tumour development , 1996 .

[39]  L. M. Franks,et al.  Introduction to the cellular and molecular biology of cancer , 1986 .

[40]  Kumbakonam R. Rajagopal,et al.  On the development and generalizations of Cahn–Hilliard equations within a thermodynamic framework , 2012 .

[41]  L E Shilkrot,et al.  Coupled atomistic and discrete dislocation plasticity. , 2002, Physical review letters.

[42]  Triantafyllos Stylianopoulos,et al.  Role of Constitutive Behavior and Tumor-Host Mechanical Interactions in the State of Stress and Growth of Solid Tumors , 2014, PloS one.

[43]  J. Miyake,et al.  Measurement of Biomolecular Diffusion in Extracellular Matrix Condensed by Fibroblasts Using Fluorescence Correlation Spectroscopy , 2013, PloS one.

[44]  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.

[45]  Michael Bergdorf,et al.  A hybrid model for three-dimensional simulations of sprouting angiogenesis. , 2008, Biophysical journal.

[46]  R Foroni,et al.  An exponential-Gompertzian description of LoVo cell tumor growth from in vivo and in vitro data. , 1989, Cancer research.

[47]  L. D. de Pillis,et al.  A cellular automata model of tumor-immune system interactions. , 2006, Journal of theoretical biology.

[48]  Daniela Calvetti,et al.  Introduction to Bayesian Scientific Computing: Ten Lectures on Subjective Computing , 2007 .

[49]  V. Giorno,et al.  A stochastic model in tumor growth. , 2006, Journal of theoretical biology.

[50]  M J Clugston,et al.  A Dictionary Of Science , 1950 .

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

[52]  Michael S. Eldred,et al.  DAKOTA : a multilevel parallel object-oriented framework for design optimization, parameter estimation, uncertainty quantification, and sensitivity analysis. Version 5.0, user's reference manual. , 2010 .

[53]  D. Lauffenburger,et al.  Physicochemical modelling of cell signalling pathways , 2006, Nature Cell Biology.

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

[55]  Charles M. Elliott,et al.  On the Cahn-Hilliard equation , 1986 .

[56]  Noam Bernstein,et al.  Spanning the continuum to quantum length scales in a dynamic simulation of brittle fracture , 1998 .

[57]  P. Selby,et al.  Introduction to the cellular and molecular biology of cancer , 2005 .

[58]  Philippe R.B. Devloo,et al.  Adaptive finite element methods for the analysis of inviscid compressible flow: Part I. Fast refinement/unrefinement and moving mesh methods for unstructured meshes , 1986 .

[59]  R. Jain,et al.  Solid stress generated by spheroid growth estimated using a linear poroelasticity model. , 2003, Microvascular research.

[60]  Andrea Hawkins-Daarud,et al.  Toward a predictive model of tumor growth , 2011 .

[61]  T. Deisboeck,et al.  Simulating non-small cell lung cancer with a multiscale agent-based model , 2007, Theoretical Biology and Medical Modelling.

[62]  Gyan Bhanot,et al.  A 2D mechanistic model of breast ductal carcinoma in situ (DCIS) morphology and progression. , 2010, Journal of theoretical biology.

[63]  Robert L. Winkler,et al.  An Introduction to Bayesian Inference and Decision , 1972 .

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

[65]  Marissa Nichole Rylander,et al.  Microfluidic culture models to study the hydrodynamics of tumor progression and therapeutic response , 2013, Biotechnology and bioengineering.

[66]  C. Lo,et al.  Stochastic Gompertz model of tumour cell growth. , 2007, Journal of theoretical biology.

[67]  S Torquato,et al.  Simulated brain tumor growth dynamics using a three-dimensional cellular automaton. , 2000, Journal of theoretical biology.

[68]  V. A. Cimmelli,et al.  A generalized Coleman–Noll procedure for the exploitation of the entropy principle , 2010, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[69]  M. Ortiz,et al.  Quasicontinuum analysis of defects in solids , 1996 .

[70]  Robert E. Rudd,et al.  COARSE-GRAINED MOLECULAR DYNAMICS AND THE ATOMIC LIMIT OF FINITE ELEMENTS , 1998 .

[71]  P. Maini,et al.  Multiscale modeling in biology , 2007 .

[72]  Michael I. Miller,et al.  Deformable templates using large deformation kinematics , 1996, IEEE Trans. Image Process..

[73]  Shannon M. Mumenthaler,et al.  Modeling Multiscale Necrotic and Calcified Tissue Biomechanics in Cancer Patients: Application to Ductal Carcinoma In Situ (DCIS) , 2013 .

[74]  Steven M. Wise,et al.  An adaptive multigrid algorithm for simulating solid tumor growth using mixture models , 2011, Math. Comput. Model..

[75]  Paul T. Bauman,et al.  A computational framework for dynamic data‐driven material damage control, based on Bayesian inference and model selection , 2015 .

[76]  Hana El-Samad,et al.  A master equation and moment approach for biochemical systems with creation-time-dependent bimolecular rate functions. , 2014, The Journal of chemical physics.

[77]  Mark A. J. Chaplain,et al.  Mathematical Modeling of Cancer Invasion: The Role of Membrane-Bound Matrix Metalloproteinases , 2013, Front. Oncol..

[78]  Matija Snuderl,et al.  Coevolution of solid stress and interstitial fluid pressure in tumors during progression: implications for vascular collapse. , 2013, Cancer research.

[79]  S. Chib,et al.  Understanding the Metropolis-Hastings Algorithm , 1995 .

[80]  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.

[81]  Antonio Criminisi,et al.  Decision Forests for Computer Vision and Medical Image Analysis , 2013, Advances in Computer Vision and Pattern Recognition.

[82]  R. Thomas,et al.  Boolean formalization of genetic control circuits. , 1973, Journal of theoretical biology.

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

[84]  G M Saidel,et al.  Stochastic model of metastases formation. , 1976, Biometrics.

[85]  Albert Tarantola,et al.  Inverse problem theory - and methods for model parameter estimation , 2004 .

[86]  C. Szot,et al.  In vitro angiogenesis induced by tumor-endothelial cell co-culture in bilayered, collagen I hydrogel bioengineered tumors. , 2013, Tissue engineering. Part C, Methods.

[87]  E. Jaynes Probability theory : the logic of science , 2003 .

[88]  H. Othmer,et al.  Mathematical modeling of tumor-induced angiogenesis , 2004, Journal of mathematical biology.

[89]  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.

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

[91]  Paul Steinmann,et al.  On higher gradients in continuum-atomistic modelling , 2003 .

[92]  H. Frieboes,et al.  Computer simulation of glioma growth and morphology , 2007, NeuroImage.

[93]  Alexander R. A. Anderson,et al.  Mathematical modelling of cancer cell invasion of tissue , 2008, Math. Comput. Model..

[94]  van der Erik Giessen,et al.  Discrete dislocation plasticity: a simple planar model , 1995 .

[95]  Xi-Qiao Feng,et al.  Biomechanical Modeling of Surface Wrinkling of Soft Tissues with Growth-Dependent Mechanical Properties , 2012 .

[96]  J. Hinchion,et al.  Releasing pressure in tumors: what do we know so far and where do we go from here? A review. , 2014, Cancer research.

[97]  Steffen Klamt,et al.  Bridging the layers: towards integration of signal transduction, regulation and metabolism into mathematical models. , 2013, Molecular bioSystems.

[98]  M. Piotrowska,et al.  A quantitative cellular automaton model of in vitro multicellular spheroid tumour growth. , 2009, Journal of theoretical biology.

[99]  Xiangpeng Zheng,et al.  Setup and characterization of a human head and neck squamous cell carcinoma xenograft model in nude rats , 2006, Otolaryngology--head and neck surgery : official journal of American Academy of Otolaryngology-Head and Neck Surgery.

[100]  R. Hesketh Introduction to Cancer Biology , 2013 .

[101]  Paul C. Fife,et al.  Models for phase separation and their mathematics. , 2000 .

[102]  Yusheng Feng,et al.  Biological Systems: Multiscale Modeling Based on Mixture Theory , 2015 .

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

[104]  Eugenia Corvera Poiré,et al.  Tumor Angiogenesis and Vascular Patterning: A Mathematical Model , 2011, PloS one.

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

[106]  D A Lauffenburger,et al.  Analysis of the roles of microvessel endothelial cell random motility and chemotaxis in angiogenesis. , 1991, Journal of theoretical biology.

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

[108]  James O. Berger,et al.  A Bayesian analysis of the thermal challenge problem , 2008 .

[109]  T. Belytschko,et al.  A bridging domain method for coupling continua with molecular dynamics , 2004 .

[110]  J. Beck,et al.  Model Selection using Response Measurements: Bayesian Probabilistic Approach , 2004 .

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

[112]  Marissa Nichole Rylander,et al.  Flow shear stress regulates endothelial barrier function and expression of angiogenic factors in a 3D microfluidic tumor vascular model , 2014, Cell adhesion & migration.

[113]  Brian B. Avants,et al.  The Multimodal Brain Tumor Image Segmentation Benchmark (BRATS) , 2015, IEEE Transactions on Medical Imaging.

[114]  R. Jain Normalizing tumor microenvironment to treat cancer: bench to bedside to biomarkers. , 2013, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

[115]  C. Patrick,et al.  Nonlinear behaviors of capillary formation in a deterministic angiogenesis model , 2005 .

[116]  R. T. Cox Probability, frequency and reasonable expectation , 1990 .

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

[118]  K. Miller,et al.  Modeling Three-Dimensional Avascular Tumor Growth Using Lattice Gas Cellular Automata , 2014 .

[119]  J. Geschwind,et al.  Tumor glycolysis as a target for cancer therapy: progress and prospects , 2013, Molecular Cancer.

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

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

[122]  D. Schoenfeld,et al.  Comparison of postoperative radiotherapy and combined postoperative radiotherapy and chemotherapy in the multidisciplinary management of malignant gliomas . A joint radiation therapy oncology group and eastern cooperative oncology group study , 1983, Cancer.

[123]  Charles M. Elliott,et al.  The Cahn-Hilliard Model for the Kinetics of Phase Separation , 1989 .

[124]  Harold S. Park,et al.  Nano Mechanics and Materials: Theory, Multiscale Methods and Applications , 2006 .

[125]  Luigi Preziosi,et al.  Cancer Modelling and Simulation , 2003 .

[126]  Shuyu Sun,et al.  A deterministic model of growth factor-induced angiogenesis , 2005, Bulletin of mathematical biology.

[127]  Vittorio Cristini,et al.  Patient-calibrated agent-based modelling of ductal carcinoma in situ (DCIS): from microscopic measurements to macroscopic predictions of clinical progression. , 2012, Journal of theoretical biology.

[128]  Zhihui Wang,et al.  Cross-scale, cross-pathway evaluation using an agent-based non-small cell lung cancer model , 2009, Bioinform..

[129]  Vittorio Cristini,et al.  Three-dimensional multispecies nonlinear tumor growth-II: Tumor invasion and angiogenesis. , 2010, Journal of theoretical biology.

[130]  Ilya M. Sobol,et al.  Sensitivity Estimates for Nonlinear Mathematical Models , 1993 .

[131]  H. Frieboes,et al.  Nonlinear modelling of cancer: bridging the gap between cells and tumours , 2010, Nonlinearity.

[132]  Lev Naumov,et al.  The influence of mitoses rate on growth dynamics of a cellular automata model of tumour growth , 2010, ICCS.

[133]  Christian Jacob,et al.  Hierarchical self-organized learning in agent-based modeling of the MAPK signaling pathway , 2011, 2011 IEEE Congress of Evolutionary Computation (CEC).

[134]  M. Chaplain,et al.  Continuous and discrete mathematical models of tumor-induced angiogenesis , 1998, Bulletin of mathematical biology.

[135]  Paul Dupuis,et al.  Variational problems on ows of di eomorphisms for image matching , 1998 .

[136]  Kevin P. Murphy,et al.  Machine learning - a probabilistic perspective , 2012, Adaptive computation and machine learning series.

[137]  S Torquato,et al.  Cellular automaton of idealized brain tumor growth dynamics. , 2000, Bio Systems.

[138]  Xiaobo Zhou,et al.  Computational Modeling of 3D Tumor Growth and Angiogenesis for Chemotherapy Evaluation , 2014, PloS one.

[139]  Walter Noll,et al.  The thermodynamics of elastic materials with heat conduction and viscosity , 1963 .

[140]  E. Somersalo,et al.  Statistical and computational inverse problems , 2004 .

[141]  Virginia Giorno,et al.  Inference on a stochastic two-compartment model in tumor growth , 2012, Comput. Stat. Data Anal..

[142]  Ted Belytschko,et al.  A bridging domain and strain computation method for coupled atomistic–continuum modelling of solids , 2007 .