1
|
Lorenzo G, Ahmed SR, Hormuth DA, Vaughn B, Kalpathy-Cramer J, Solorio L, Yankeelov TE, Gomez H. Patient-Specific, Mechanistic Models of Tumor Growth Incorporating Artificial Intelligence and Big Data. Annu Rev Biomed Eng 2024; 26:529-560. [PMID: 38594947 DOI: 10.1146/annurev-bioeng-081623-025834] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/11/2024]
Abstract
Despite the remarkable advances in cancer diagnosis, treatment, and management over the past decade, malignant tumors remain a major public health problem. Further progress in combating cancer may be enabled by personalizing the delivery of therapies according to the predicted response for each individual patient. The design of personalized therapies requires the integration of patient-specific information with an appropriate mathematical model of tumor response. A fundamental barrier to realizing this paradigm is the current lack of a rigorous yet practical mathematical theory of tumor initiation, development, invasion, and response to therapy. We begin this review with an overview of different approaches to modeling tumor growth and treatment, including mechanistic as well as data-driven models based on big data and artificial intelligence. We then present illustrative examples of mathematical models manifesting their utility and discuss the limitations of stand-alone mechanistic and data-driven models. We then discuss the potential of mechanistic models for not only predicting but also optimizing response to therapy on a patient-specific basis. We describe current efforts and future possibilities to integrate mechanistic and data-driven models. We conclude by proposing five fundamental challenges that must be addressed to fully realize personalized care for cancer patients driven by computational models.
Collapse
Affiliation(s)
- Guillermo Lorenzo
- Oden Institute for Computational Engineering and Sciences, University of Texas, Austin, Texas, USA
- Department of Civil Engineering and Architecture, University of Pavia, Pavia, Italy
| | - Syed Rakin Ahmed
- Geisel School of Medicine, Dartmouth College, Hanover, New Hampshire, USA
- Health Sciences and Technology, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA
- Harvard Graduate Program in Biophysics, Harvard Medical School, Harvard University, Cambridge, Massachusetts, USA
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Boston, Massachusetts, USA
| | - David A Hormuth
- Livestrong Cancer Institutes, University of Texas, Austin, Texas, USA
- Oden Institute for Computational Engineering and Sciences, University of Texas, Austin, Texas, USA
| | - Brenna Vaughn
- Weldon School of Biomedical Engineering, Purdue University, West Lafayette, Indiana, USA;
| | | | - Luis Solorio
- Weldon School of Biomedical Engineering, Purdue University, West Lafayette, Indiana, USA;
| | - Thomas E Yankeelov
- Department of Imaging Physics, MD Anderson Cancer Center, Houston, Texas, USA
- Department of Biomedical Engineering, Department of Oncology, and Department of Diagnostic Medicine, University of Texas, Austin, Texas, USA
- Livestrong Cancer Institutes, University of Texas, Austin, Texas, USA
- Oden Institute for Computational Engineering and Sciences, University of Texas, Austin, Texas, USA
| | - Hector Gomez
- School of Mechanical Engineering and Purdue Center for Cancer Research, Purdue University, West Lafayette, Indiana, USA
- Weldon School of Biomedical Engineering, Purdue University, West Lafayette, Indiana, USA;
| |
Collapse
|
2
|
Pajic-Lijakovic I, Eftimie R, Milivojevic M, Bordas SPA. Segregation of co-cultured multicellular systems: review and modeling consideration. Q Rev Biophys 2024; 57:e5. [PMID: 38351868 DOI: 10.1017/s0033583524000015] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/04/2024]
Abstract
Cell segregation caused by collective cell migration (CCM) is crucial for morphogenesis, functional development of tissue parts, and is an important aspect in other diseases such as cancer and its metastasis process. Efficiency of the cell segregation depends on the interplay between: (1) biochemical processes such as cell signaling and gene expression and (2) physical interactions between cells. Despite extensive research devoted to study the segregation of various co-cultured systems, we still do not understand the role of physical interactions in cell segregation. Cumulative effects of these physical interactions appear in the form of physical parameters such as: (1) tissue surface tension, (2) viscoelasticity caused by CCM, and (3) solid stress accumulated in multicellular systems. These parameters primarily depend on the interplay between the state of cell-cell adhesion contacts and cell contractility. The role of these physical parameters on the segregation efficiency is discussed on model systems such as co-cultured breast cell spheroids consisting of two subpopulations that are in contact. This review study aims to: (1) summarize biological aspects related to cell segregation, mechanical properties of cell collectives, effects along the biointerface between cell subpopulations and (2) describe from a biophysical/mathematical perspective the same biological aspects summarized before. So that overall it can illustrate the complexity of the biological systems that translate into very complex biophysical/mathematical equations. Moreover, by presenting in parallel these two seemingly different parts (biology vs. equations), this review aims to emphasize the need for experiments to estimate the variety of parameters entering the resulting complex biophysical/mathematical models.
Collapse
Affiliation(s)
- Ivana Pajic-Lijakovic
- Faculty of Technology and Metallurgy, Department of Chemical Engineering, University of Belgrade, Beograd, Serbia
| | - Raluca Eftimie
- Laboratoire Mathematiques de Besançon, UMR-CNRS 6623, Université de Bourgogne Franche-Comte, Besançon, France
| | - Milan Milivojevic
- Faculty of Technology and Metallurgy, Department of Chemical Engineering, University of Belgrade, Beograd, Serbia
| | - Stéphane P A Bordas
- Faculty of Science, Technology and Communication, University of Luxembourg, Institute for Computational Engineering, Esch-sur-Alzette, Luxembourg
- Department of Medical Research, China Medical University Hospital, Taichung, Taiwan
| |
Collapse
|
3
|
Eftimie R, Rolin G, Adebayo OE, Urcun S, Chouly F, Bordas SPA. Modelling Keloids Dynamics: A Brief Review and New Mathematical Perspectives. Bull Math Biol 2023; 85:117. [PMID: 37855947 DOI: 10.1007/s11538-023-01222-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2023] [Accepted: 10/02/2023] [Indexed: 10/20/2023]
Abstract
Keloids are fibroproliferative disorders described by excessive growth of fibrotic tissue, which also invades adjacent areas (beyond the original wound borders). Since these disorders are specific to humans (no other animal species naturally develop keloid-like tissue), experimental in vivo/in vitro research has not led to significant advances in this field. One possible approach could be to combine in vitro human models with calibrated in silico mathematical approaches (i.e., models and simulations) to generate new testable biological hypotheses related to biological mechanisms and improved treatments. Because these combined approaches do not really exist for keloid disorders, in this brief review we start by summarising the biology of these disorders, then present various types of mathematical and computational approaches used for related disorders (i.e., wound healing and solid tumours), followed by a discussion of the very few mathematical and computational models published so far to study various inflammatory and mechanical aspects of keloids. We conclude this review by discussing some open problems and mathematical opportunities offered in the context of keloid disorders by such combined in vitro/in silico approaches, and the need for multi-disciplinary research to enable clinical progress.
Collapse
Affiliation(s)
- R Eftimie
- Laboratoire de Mathématiques de Besançon, Université de Franche-Comté, 25000, Besançon, France.
| | - G Rolin
- INSERM CIC-1431, CHU Besançon, F-25000, Besançon, France
- EFS, INSERM, UMR 1098 RIGHT, Université de Franche-Comté, F-25000, Besançon, France
| | - O E Adebayo
- Laboratoire de Mathématiques de Besançon, Université de Franche-Comté, 25000, Besançon, France
| | - S Urcun
- Institute for Computational Engineering, Faculty of Science, Technology and Communication, University of Luxembourg, Esch-sur-Alzette, Luxembourg
| | - F Chouly
- Institut de Mathématiques de Bourgogne, Université de Franche-Comté, 21078, Dijon, France
- Center for Mathematical Modelling and Department of Mathematical Engineering, University of Chile and IRL 2807 - CNRS, Santiago, Chile
- Departamento de Ingeniería Matemática, CI2MA, Universidad de Concepción, Casilla 160-C, Concepción, Chile
| | - S P A Bordas
- Institute for Computational Engineering, Faculty of Science, Technology and Communication, University of Luxembourg, Esch-sur-Alzette, Luxembourg
| |
Collapse
|
4
|
Lavigne T, Urcun S, Rohan PY, Sciumè G, Baroli D, Bordas SPA. Single and bi-compartment poro-elastic model of perfused biological soft tissues: FEniCSx implementation and tutorial. J Mech Behav Biomed Mater 2023; 143:105902. [PMID: 37209595 DOI: 10.1016/j.jmbbm.2023.105902] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/13/2023] [Revised: 04/26/2023] [Accepted: 05/07/2023] [Indexed: 05/22/2023]
Abstract
Soft biological tissues demonstrate strong time-dependent and strain-rate mechanical behavior, arising from their intrinsic visco-elasticity and fluid-solid interactions. The time-dependent mechanical properties of soft tissues influence their physiological functions and are related to several pathological processes. Poro-elastic modeling represents a promising approach because it allows the integration of multiscale/multiphysics data to probe biologically relevant phenomena at a smaller scale and embeds the relevant mechanisms at the larger scale. The implementation of multiphase flow poro-elastic models however is a complex undertaking, requiring extensive knowledge. The open-source software FEniCSx Project provides a novel tool for the automated solution of partial differential equations by the finite element method. This paper aims to provide the required tools to model the mixed formulation of poro-elasticity, from the theory to the implementation, within FEniCSx. Several benchmark cases are studied. A column under confined compression conditions is compared to the Terzaghi analytical solution, using the L2-norm. An implementation of poro-hyper-elasticity is proposed. A bi-compartment column is compared to previously published results (Cast3m implementation). For all cases, accurate results are obtained in terms of a normalized Root Mean Square Error (RMSE). Furthermore, the FEniCSx computation is found three times faster than the legacy FEniCS one. The benefits of parallel computation are also highlighted.
Collapse
Affiliation(s)
- Thomas Lavigne
- Institute of Computational Engineering, Department of Engineering, University of Luxembourg, 6, avenue de la Fonte, Esch-sur-Alzette, L-4364, Luxembourg; Arts et Metiers Institute of Technology, IBHGC, 151 bd de l'hopital, Paris, 75013, France; Arts et Metiers Institute of Technology, Univ. of Bordeaux, CNRS, Bordeaux INP, INRAE, I2M Bordeaux, Avenue d'Aquitaine, Pessac, 33607, France
| | - Stéphane Urcun
- Institute of Computational Engineering, Department of Engineering, University of Luxembourg, 6, avenue de la Fonte, Esch-sur-Alzette, L-4364, Luxembourg
| | - Pierre-Yves Rohan
- Arts et Metiers Institute of Technology, IBHGC, 151 bd de l'hopital, Paris, 75013, France
| | - Giuseppe Sciumè
- Arts et Metiers Institute of Technology, Univ. of Bordeaux, CNRS, Bordeaux INP, INRAE, I2M Bordeaux, Avenue d'Aquitaine, Pessac, 33607, France
| | - Davide Baroli
- Università della Svizzera Italiana, Euler Institute, Lugano, Switzerland.
| | - Stéphane P A Bordas
- Institute of Computational Engineering, Department of Engineering, University of Luxembourg, 6, avenue de la Fonte, Esch-sur-Alzette, L-4364, Luxembourg; Clyde Visiting Fellow, Department of Mechanical Engineering, The University of Utah, Salt Lake City, UT, United States; Visiting Department of Medical Research, China Medical University Hospital, China Medical University, Taichung, Taiwan
| |
Collapse
|
5
|
Hervas-Raluy S, Wirthl B, Guerrero PE, Robalo Rei G, Nitzler J, Coronado E, Font de Mora Sainz J, Schrefler BA, Gomez-Benito MJ, Garcia-Aznar JM, Wall WA. Tumour growth: An approach to calibrate parameters of a multiphase porous media model based on in vitro observations of Neuroblastoma spheroid growth in a hydrogel microenvironment. Comput Biol Med 2023; 159:106895. [PMID: 37060771 DOI: 10.1016/j.compbiomed.2023.106895] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2023] [Revised: 03/09/2023] [Accepted: 04/09/2023] [Indexed: 04/17/2023]
Abstract
To unravel processes that lead to the growth of solid tumours, it is necessary to link knowledge of cancer biology with the physical properties of the tumour and its interaction with the surrounding microenvironment. Our understanding of the underlying mechanisms is however still imprecise. We therefore developed computational physics-based models, which incorporate the interaction of the tumour with its surroundings based on the theory of porous media. However, the experimental validation of such models represents a challenge to its clinical use as a prognostic tool. This study combines a physics-based model with in vitro experiments based on microfluidic devices used to mimic a three-dimensional tumour microenvironment. By conducting a global sensitivity analysis, we identify the most influential input parameters and infer their posterior distribution based on Bayesian calibration. The resulting probability density is in agreement with the scattering of the experimental data and thus validates the proposed workflow. This study demonstrates the huge challenges associated with determining precise parameters with usually only limited data for such complex processes and models, but also demonstrates in general how to indirectly characterise the mechanical properties of neuroblastoma spheroids that cannot feasibly be measured experimentally.
Collapse
Affiliation(s)
- Silvia Hervas-Raluy
- Multiscale in Mechanical and Biological Engineering, Department of Mechanical Engineering, University of Zaragoza, Aragon Institute for Engineering Research (I3A), Maria de Luna 3, Zaragoza, 50018, Spain.
| | - Barbara Wirthl
- Institute for Computational Mechanics, Technical University of Munich, TUM School of Engineering and Design, Department of Engineering Physics & Computation, Boltzmannstraße 15, Garching b. Munich, 85748, Germany
| | - Pedro E Guerrero
- Multiscale in Mechanical and Biological Engineering, Department of Mechanical Engineering, University of Zaragoza, Aragon Institute for Engineering Research (I3A), Maria de Luna 3, Zaragoza, 50018, Spain
| | - Gil Robalo Rei
- Institute for Computational Mechanics, Technical University of Munich, TUM School of Engineering and Design, Department of Engineering Physics & Computation, Boltzmannstraße 15, Garching b. Munich, 85748, Germany
| | - Jonas Nitzler
- Institute for Computational Mechanics, Technical University of Munich, TUM School of Engineering and Design, Department of Engineering Physics & Computation, Boltzmannstraße 15, Garching b. Munich, 85748, Germany; Professorship for Data-Driven Materials Modeling, Technical University of Munich, TUM School of Engineering and Design, Department of Engineering Physics & Computation, Boltzmannstraße 15, Garching b. Munich, 85748, Germany
| | - Esther Coronado
- Clinical and Translational Oncology Research Group, Instituto de Investigación La Fe,, Fernando Abril Martorell 106, Valencia, 46026, Spain
| | - Jaime Font de Mora Sainz
- Clinical and Translational Oncology Research Group, Instituto de Investigación La Fe,, Fernando Abril Martorell 106, Valencia, 46026, Spain
| | - Bernhard A Schrefler
- Department of Civil, Environmental and Architectural Engineering, University of Padua, Marzolo 9, Padua, 35131, Italy; Institute for Advanced Study, Technical University of Munich, Boltzmannstraße 15, Garching b. Munich, 85748, Germany
| | - Maria Jose Gomez-Benito
- Multiscale in Mechanical and Biological Engineering, Department of Mechanical Engineering, University of Zaragoza, Aragon Institute for Engineering Research (I3A), Maria de Luna 3, Zaragoza, 50018, Spain
| | - Jose Manuel Garcia-Aznar
- Multiscale in Mechanical and Biological Engineering, Department of Mechanical Engineering, University of Zaragoza, Aragon Institute for Engineering Research (I3A), Maria de Luna 3, Zaragoza, 50018, Spain
| | - Wolfgang A Wall
- Institute for Computational Mechanics, Technical University of Munich, TUM School of Engineering and Design, Department of Engineering Physics & Computation, Boltzmannstraße 15, Garching b. Munich, 85748, Germany
| |
Collapse
|
6
|
Urcun S, Baroli D, Rohan PY, Skalli W, Lubrano V, Bordas SP, Sciumè G. Non-operable glioblastoma: Proposition of patient-specific forecasting by image-informed poromechanical model. BRAIN MULTIPHYSICS 2023. [DOI: 10.1016/j.brain.2023.100067] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/03/2023] Open
|
7
|
Lavigne T, Sciumè G, Laporte S, Pillet H, Urcun S, Wheatley B, Rohan PY. Société de Biomécanique Young Investigator Award 2021: Numerical investigation of the time-dependent stress-strain mechanical behaviour of skeletal muscle tissue in the context of pressure ulcer prevention. Clin Biomech (Bristol, Avon) 2022; 93:105592. [PMID: 35151107 DOI: 10.1016/j.clinbiomech.2022.105592] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 11/09/2021] [Revised: 01/31/2022] [Accepted: 02/01/2022] [Indexed: 02/07/2023]
Abstract
BACKGROUND Pressure-induced tissue strain is one major pathway for Pressure Ulcer development and, especially, Deep Tissue Injury. Biomechanical investigation of the time-dependent stress-strain mechanical behaviour of skeletal muscle tissue is therefore essential. In the literature, a viscoelastic formulation is generally assumed for the experimental characterization of skeletal muscles, with the limitation that the underlying physical mechanisms that give rise to the time dependent stress-strain behaviour are not known. The objective of this study is to explore the capability of poroelasticity to reproduce the apparent viscoelastic behaviour of passive muscle tissue under confined compression. METHODS Experimental stress-relaxation response of 31 cylindrical porcine samples tested under fast and slow confined compression by Vaidya and collaborators were used. An axisymmetric Finite Element model was developed in ABAQUS and, for each sample a one-to-one inverse analysis was performed to calibrate the specimen-specific constitutive parameters, namely, the drained Young's modulus, the void ratio, hydraulic permeability, the Poisson's ratio, the solid grain's and fluid's bulk moduli. FINDINGS The peak stress and consolidation were recovered for most of the samples (N=25) by the poroelastic model (normalised root-mean-square error ≤0.03 for fast and slow confined compression conditions). INTERPRETATION The strength of the proposed model is its fewer number of variables (N=6 for the proposed poroelastic model versus N=18 for the viscohyperelastic model proposed by Vaidya and collaborators). The incorporation of poroelasticity to clinical models of Pessure Ulcer formation could lead to more precise and mechanistic explorations of soft tissue injury risk factors.
Collapse
Affiliation(s)
- T Lavigne
- Arts et Metiers Institute of Technology, IBHGC, 151 bd de l'hopital, Paris, 75013, France; Arts et Metiers Institute of Technology, Univ. of Bordeaux, CNRS, Bordeaux INP, INRAE, I2M Bordeaux, Avenue d'Aquitaine, Pessac, 33607, France.
| | - G Sciumè
- Arts et Metiers Institute of Technology, Univ. of Bordeaux, CNRS, Bordeaux INP, INRAE, I2M Bordeaux, Avenue d'Aquitaine, Pessac, 33607, France
| | - S Laporte
- Arts et Metiers Institute of Technology, IBHGC, 151 bd de l'hopital, Paris, 75013, France
| | - H Pillet
- Arts et Metiers Institute of Technology, IBHGC, 151 bd de l'hopital, Paris, 75013, France
| | - S Urcun
- Arts et Metiers Institute of Technology, IBHGC, 151 bd de l'hopital, Paris, 75013, France; Arts et Metiers Institute of Technology, Univ. of Bordeaux, CNRS, Bordeaux INP, INRAE, I2M Bordeaux, Avenue d'Aquitaine, Pessac, 33607, France; Institute for Computational Engineering Sciences, Department of Engineering Sciences, Faculte des Sciences, de la Technologie et de Medecine, Universite du Luxembourg, Campus Kirchberg, 6, rue Coudenhove-Kalergi, Luxembourg, L-1359, Luxembourg
| | - B Wheatley
- Department of Mechanical Engineering, Bucknell University, 1 Dent Drive, Lewisburg 17837, PA, USA
| | - P-Y Rohan
- Arts et Metiers Institute of Technology, IBHGC, 151 bd de l'hopital, Paris, 75013, France
| |
Collapse
|