1
|
Mohammadi M, Soltani M, Aghanajafi C, Kohandel M. Investigation of the evolution of tumor-induced microvascular network under the inhibitory effect of anti-angiogenic factor, angiostatin: A mathematical study. MATHEMATICAL BIOSCIENCES AND ENGINEERING : MBE 2023; 20:5448-5480. [PMID: 36896553 DOI: 10.3934/mbe.2023252] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/18/2023]
Abstract
Anti-angiogenesis as a treatment strategy for normalizing the microvascular network of tumors is of great interest among researchers, especially in combination with chemotherapy or radiotherapy. According to the vital role that angiogenesis plays in tumor growth and in exposing the tumor to therapeutic agents, this work develops a mathematical framework to study the influence of angiostatin, a plasminogen fragment that shows the anti-angiogenic function, in the evolutionary behavior of tumor-induced angiogenesis. Angiostatin-induced microvascular network reformation is investigated in a two-dimensional space by considering two parent vessels around a circular tumor by a modified discrete angiogenesis model in different tumor sizes. The effects of imposing modifications on the existing model, i.e., the matrix-degrading enzyme effect, proliferation and death of endothelial cells, matrix density function, and a more realistic chemotactic function, are investigated in this study. Results show a decrease in microvascular density in response to the angiostatin. A functional relationship exists between angiostatin's ability to normalize the capillary network and tumor size or progression stage, such that capillary density decreases by 55%, 41%, 24%, and 13% in tumors with a non-dimensional radius of 0.4, 0.3, 0.2, and 0.1, respectively, after angiostatin administration.
Collapse
Affiliation(s)
- Mahya Mohammadi
- Department of Mechanical Engineering, K. N. Toosi University of Technology, Tehran 19919-43344, Iran
- Department of Applied Mathematics, University of Waterloo, Waterloo, ON N2L 3G1, Canada
| | - M Soltani
- Department of Mechanical Engineering, K. N. Toosi University of Technology, Tehran 19919-43344, Iran
- Department of Electrical and Computer Engineering, University of Waterloo, Waterloo, ON N2L 3G1, Canada
- School of Optometry and Vision Science, Faculty of Science, University of Waterloo, Waterloo, ON N2L 3G1, Canada
- Advanced Bioengineering Initiative Center, Multidisciplinary International Complex, K. N. Toosi University of Technology, Tehran 19697-64499, Iran
- Centre for Biotechnology and Bioengineering (CBB), University of Waterloo, Waterloo, ON N2L 3G1, Canada
| | - Cyrus Aghanajafi
- Department of Mechanical Engineering, K. N. Toosi University of Technology, Tehran 19919-43344, Iran
| | - Mohammad Kohandel
- Department of Applied Mathematics, University of Waterloo, Waterloo, ON N2L 3G1, Canada
| |
Collapse
|
2
|
Abstract
Cancer cells require higher oxygen levels and nutrition than normal cells. Cancer cells induce angiogenesis (the development of new blood vessels) from preexisting vessels. This biological process depends on the special, chemical, and physical properties of the microenvironment surrounding tumor tissues. The complexity of these properties hinders an understanding of their mechanisms. Various mathematical models have been developed to describe quantitative relationships related to angiogenesis. We developed a three-dimensional mathematical model that incorporates angiogenesis and tumor growth. We examined angiopoietin, which regulates the spouting and branching events in angiogenesis. The simulation successfully reproduced the transient decrease in new vessels during vascular network formation. This chapter describes the protocol used to perform the simulations.
Collapse
Affiliation(s)
- Masahiro Sugimoto
- Institute of Medical Science, Tokyo Medical University, Shinjuku, Tokyo, Japan.
- Institute for Advanced Biosciences, Keio University, Tsuruoka, Yamagata, Japan.
| |
Collapse
|
3
|
Numerical Investigation on the Anti-Angiogenic Therapy-Induced Normalization in Solid Tumors. Pharmaceutics 2022; 14:pharmaceutics14020363. [PMID: 35214095 PMCID: PMC8877966 DOI: 10.3390/pharmaceutics14020363] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/12/2021] [Revised: 01/29/2022] [Accepted: 01/31/2022] [Indexed: 01/27/2023] Open
Abstract
This study numerically analyzes the fluid flow and solute transport in a solid tumor to comprehensively examine the consequence of normalization induced by anti-angiogenic therapy on drug delivery. The current study leads to a more accurate model in comparison to previous research, as it incorporates a non-homogeneous real-human solid tumor including necrotic, semi-necrotic, and well-vascularized regions. Additionally, the model considers the effects of concurrently chemotherapeutic agents (three macromolecules of IgG, F(ab′)2, and F(ab′)) and different normalization intensities in various tumor sizes. Examining the long-term influence of normalization on the quality of drug uptake by necrotic area is another contribution of the present study. Results show that normalization decreases the interstitial fluid pressure (IFP) and spreads the pressure gradient and non-zero interstitial fluid velocity (IFV) into inner areas. Subsequently, wash-out of the drug from the tumor periphery is decreased. It is also demonstrated that normalization can improve the distribution of solute concentration in the interstitium. The efficiency of normalization is introduced as a function of the time course of perfusion, which depends on the tumor size, drug type, as well as normalization intensity, and consequently on the dominant mechanism of drug delivery. It is suggested to accompany anti-angiogenic therapy by F(ab′) in large tumor size (Req=2.79 cm) to improve reservoir behavior benefit from normalization. However, IgG is proposed as the better option in the small tumor (Req=0.46 cm), in which normalization finds the opportunity of enhancing uniformity of IgG average exposure by 22%. This study could provide a perspective for preclinical and clinical trials on how to take advantage of normalization, as an adjuvant treatment, in improving drug delivery into a non-homogeneous solid tumor.
Collapse
|
4
|
Bayesian calibration of a stochastic, multiscale agent-based model for predicting in vitro tumor growth. PLoS Comput Biol 2021; 17:e1008845. [PMID: 34843457 PMCID: PMC8659698 DOI: 10.1371/journal.pcbi.1008845] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/24/2021] [Revised: 12/09/2021] [Accepted: 11/02/2021] [Indexed: 12/31/2022] Open
Abstract
Hybrid multiscale agent-based models (ABMs) are unique in their ability to simulate individual cell interactions and microenvironmental dynamics. Unfortunately, the high computational cost of modeling individual cells, the inherent stochasticity of cell dynamics, and numerous model parameters are fundamental limitations of applying such models to predict tumor dynamics. To overcome these challenges, we have developed a coarse-grained two-scale ABM (cgABM) with a reduced parameter space that allows for an accurate and efficient calibration using a set of time-resolved microscopy measurements of cancer cells grown with different initial conditions. The multiscale model consists of a reaction-diffusion type model capturing the spatio-temporal evolution of glucose and growth factors in the tumor microenvironment (at tissue scale), coupled with a lattice-free ABM to simulate individual cell dynamics (at cellular scale). The experimental data consists of BT474 human breast carcinoma cells initialized with different glucose concentrations and tumor cell confluences. The confluence of live and dead cells was measured every three hours over four days. Given this model, we perform a time-dependent global sensitivity analysis to identify the relative importance of the model parameters. The subsequent cgABM is calibrated within a Bayesian framework to the experimental data to estimate model parameters, which are then used to predict the temporal evolution of the living and dead cell populations. To this end, a moment-based Bayesian inference is proposed to account for the stochasticity of the cgABM while quantifying uncertainties due to limited temporal observational data. The cgABM reduces the computational time of ABM simulations by 93% to 97% while staying within a 3% difference in prediction compared to ABM. Additionally, the cgABM can reliably predict the temporal evolution of breast cancer cells observed by the microscopy data with an average error and standard deviation for live and dead cells being 7.61±2.01 and 5.78±1.13, respectively. The calibration of agent-based models of tumor cell growth to experimental data remains a challenge in computational oncology. Besides the computational cost of modeling thousands of agents, the model’s intrinsic stochasticity demands numerous realizations of the simulations to accurately represent the statistical features of the model predictions. We developed a hybrid, multiscale, coarse-grain, agent-based model that captures the growth and decline of human breast carcinoma cells under different initial conditions. We determined the effects of coarse-graining the ABM on the multiscale model output and the number of repetitions necessary to capture the stochastic transitions present in the model. We identified the most influential parameters on the model prediction through a sensitivity analysis and selected which parameters can be fixed and which ones should be calibrated. Using Bayesian calibration, we show that the model can accurately represent the experimental data. The validation step indicates that our model can reliably predict the in vitro temporal data, depending on the choice of the training (calibration data) sets.
Collapse
|
5
|
Hormuth DA, Phillips CM, Wu C, Lima EABF, Lorenzo G, Jha PK, Jarrett AM, Oden JT, Yankeelov TE. Biologically-Based Mathematical Modeling of Tumor Vasculature and Angiogenesis via Time-Resolved Imaging Data. Cancers (Basel) 2021; 13:3008. [PMID: 34208448 PMCID: PMC8234316 DOI: 10.3390/cancers13123008] [Citation(s) in RCA: 26] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2021] [Revised: 06/07/2021] [Accepted: 06/13/2021] [Indexed: 01/03/2023] Open
Abstract
Tumor-associated vasculature is responsible for the delivery of nutrients, removal of waste, and allowing growth beyond 2-3 mm3. Additionally, the vascular network, which is changing in both space and time, fundamentally influences tumor response to both systemic and radiation therapy. Thus, a robust understanding of vascular dynamics is necessary to accurately predict tumor growth, as well as establish optimal treatment protocols to achieve optimal tumor control. Such a goal requires the intimate integration of both theory and experiment. Quantitative and time-resolved imaging methods have emerged as technologies able to visualize and characterize tumor vascular properties before and during therapy at the tissue and cell scale. Parallel to, but separate from those developments, mathematical modeling techniques have been developed to enable in silico investigations into theoretical tumor and vascular dynamics. In particular, recent efforts have sought to integrate both theory and experiment to enable data-driven mathematical modeling. Such mathematical models are calibrated by data obtained from individual tumor-vascular systems to predict future vascular growth, delivery of systemic agents, and response to radiotherapy. In this review, we discuss experimental techniques for visualizing and quantifying vascular dynamics including magnetic resonance imaging, microfluidic devices, and confocal microscopy. We then focus on the integration of these experimental measures with biologically based mathematical models to generate testable predictions.
Collapse
Affiliation(s)
- David A. Hormuth
- Oden Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, TX 78712, USA; (C.M.P.); (C.W.); (E.A.B.F.L.); (G.L.); (P.K.J.); (J.T.O.); (T.E.Y.)
- Livestrong Cancer Institutes, Dell Medical School, The University of Texas at Austin, Austin, TX 78712, USA
| | - Caleb M. Phillips
- Oden Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, TX 78712, USA; (C.M.P.); (C.W.); (E.A.B.F.L.); (G.L.); (P.K.J.); (J.T.O.); (T.E.Y.)
| | - Chengyue Wu
- Oden Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, TX 78712, USA; (C.M.P.); (C.W.); (E.A.B.F.L.); (G.L.); (P.K.J.); (J.T.O.); (T.E.Y.)
| | - Ernesto A. B. F. Lima
- Oden Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, TX 78712, USA; (C.M.P.); (C.W.); (E.A.B.F.L.); (G.L.); (P.K.J.); (J.T.O.); (T.E.Y.)
- Texas Advanced Computing Center, The University of Texas at Austin, Austin, TX 78758, USA
| | - Guillermo Lorenzo
- Oden Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, TX 78712, USA; (C.M.P.); (C.W.); (E.A.B.F.L.); (G.L.); (P.K.J.); (J.T.O.); (T.E.Y.)
- Department of Civil Engineering and Architecture, University of Pavia, Via Ferrata 3, 27100 Pavia, Italy
| | - Prashant K. Jha
- Oden Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, TX 78712, USA; (C.M.P.); (C.W.); (E.A.B.F.L.); (G.L.); (P.K.J.); (J.T.O.); (T.E.Y.)
| | - Angela M. Jarrett
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, TX 78712, USA;
| | - J. Tinsley Oden
- Oden Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, TX 78712, USA; (C.M.P.); (C.W.); (E.A.B.F.L.); (G.L.); (P.K.J.); (J.T.O.); (T.E.Y.)
- Department of Aerospace Engineering and Engineering Mechanics, The University of Texas at Austin, Austin, TX 78712, USA
- Department of Mathematics, The University of Texas at Austin, Austin, TX 78712, USA
- Department of Computer Science, The University of Texas at Austin, Austin, TX 78712, USA
| | - Thomas E. Yankeelov
- Oden Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, TX 78712, USA; (C.M.P.); (C.W.); (E.A.B.F.L.); (G.L.); (P.K.J.); (J.T.O.); (T.E.Y.)
- Livestrong Cancer Institutes, Dell Medical School, The University of Texas at Austin, Austin, TX 78712, USA
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, TX 78712, USA;
- Department of Diagnostic Medicine, The University of Texas at Austin, Austin, TX 78712, USA
- Department of Oncology, The University of Texas at Austin, Austin, TX 78712, USA
- Department of Imaging Physics, The University of Texas MD Anderson Cancer Center, Houston, TX 77030, USA
| |
Collapse
|
6
|
Mathematical simulation of tumour angiogenesis: angiopoietin balance is a key factor in vessel growth and regression. Sci Rep 2021; 11:419. [PMID: 33432093 PMCID: PMC7801613 DOI: 10.1038/s41598-020-79824-8] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/25/2020] [Accepted: 12/10/2020] [Indexed: 12/12/2022] Open
Abstract
Excessive tumour growth results in a hypoxic environment around cancer cells, thus inducing tumour angiogenesis, which refers to the generation of new blood vessels from pre-existing vessels. This mechanism is biologically and physically complex, with various mathematical simulation models proposing to reproduce its formation. However, although temporary vessel regression is clinically known, few models succeed in reproducing this phenomenon. Here, we developed a three-dimensional simulation model encompassing both angiogenesis and tumour growth, specifically including angiopoietin. Angiopoietin regulates both adhesion and migration between vascular endothelial cells and wall cells, thus inhibiting the cell-to-cell adhesion required for angiogenesis initiation. Simulation results showed a regression, i.e. transient decrease, in the overall length of new vessels during vascular network formation. Using our model, we also evaluated the efficacy of administering the drug bevacizumab. The results highlighted differences in treatment efficacy: (1) earlier administration showed higher efficacy in inhibiting tumour growth, and (2) efficacy depended on the treatment interval even with the administration of the same dose. After thorough validation in the future, these results will contribute to the design of angiogenesis treatment protocols.
Collapse
|
7
|
A hybrid model of tumor growth and angiogenesis: In silico experiments. PLoS One 2020; 15:e0231137. [PMID: 32275674 PMCID: PMC7147760 DOI: 10.1371/journal.pone.0231137] [Citation(s) in RCA: 35] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2019] [Accepted: 03/16/2020] [Indexed: 12/18/2022] Open
Abstract
Tumor associated angiogenesis is the development of new blood vessels in response to proteins secreted by tumor cells. These new blood vessels allow tumors to continue to grow beyond what the pre-existing vasculature could support. Here, we construct a mathematical model to simulate tumor angiogenesis by considering each endothelial cell as an agent, and allowing the vascular endothelial growth factor (VEGF) and nutrient fields to impact the dynamics and phenotypic transitions of each tumor and endothelial cell. The phenotypes of the endothelial cells (i.e., tip, stalk, and phalanx cells) are selected by the local VEGF field, and govern the migration and growth of vessel sprouts at the cellular level. Over time, these vessels grow and migrate to the tumor, forming anastomotic loops to supply nutrients, while interacting with the tumor through mechanical forces and the consumption of VEGF. The model is able to capture collapsing and breaking of vessels caused by tumor-endothelial cell interactions. This is accomplished through modeling the physical interaction between the vasculature and the tumor, resulting in vessel occlusion and tumor heterogeneity over time due to the stages of response in angiogenesis. Key parameters are identified through a sensitivity analysis based on the Sobol method, establishing which parameters should be the focus of subsequent experimental efforts. During the avascular phase (i.e., before angiogenesis is triggered), the nutrient consumption rate, followed by the rate of nutrient diffusion, yield the greatest influence on the number and distribution of tumor cells. Similarly, the consumption and diffusion of VEGF yield the greatest influence on the endothelial and tumor cell numbers during angiogenesis. In summary, we present a hybrid mathematical approach that characterizes vascular changes via an agent-based model, while treating nutrient and VEGF changes through a continuum model. The model describes the physical interaction between a tumor and the surrounding blood vessels, explicitly allowing the forces of the growing tumor to influence the nutrient delivery of the vasculature.
Collapse
|
8
|
Datta P, Barui A, Wu Y, Ozbolat V, Moncal KK, Ozbolat IT. Essential steps in bioprinting: From pre- to post-bioprinting. Biotechnol Adv 2018; 36:1481-1504. [PMID: 29909085 DOI: 10.1016/j.biotechadv.2018.06.003] [Citation(s) in RCA: 68] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/07/2018] [Revised: 05/15/2018] [Accepted: 06/10/2018] [Indexed: 12/17/2022]
Abstract
An increasing demand for directed assembly of biomaterials has inspired the development of bioprinting, which facilitates the assembling of both cellular and acellular inks into well-arranged three-dimensional (3D) structures for tissue fabrication. Although great advances have been achieved in the recent decade, there still exist issues to be addressed. Herein, a review has been systematically performed to discuss the considerations in the entire procedure of bioprinting. Though bioprinting is advancing at a rapid pace, it is seen that the whole process of obtaining tissue constructs from this technique involves multiple-stages, cutting across various technology domains. These stages can be divided into three broad categories: pre-bioprinting, bioprinting and post-bioprinting. Each stage can influence others and has a bearing on the performance of fabricated constructs. For example, in pre-bioprinting, tissue biopsy and cell expansion techniques are essential to ensure a large number of cells are available for mass organ production. Similarly, medical imaging is needed to provide high resolution designs, which can be faithfully bioprinted. In the bioprinting stage, compatibility of biomaterials is needed to be matched with solidification kinetics to ensure constructs with high cell viability and fidelity are obtained. On the other hand, there is a need to develop bioprinters, which have high degrees of freedom of movement, perform without failure concerns for several hours and are compact, and affordable. Finally, maturation of bioprinted cells are governed by conditions provided during the post-bioprinting process. This review, for the first time, puts all the bioprinting stages in perspective of the whole process of bioprinting, and analyzes their current state-of-the art. It is concluded that bioprinting community will recognize the relative importance and optimize the parameter of each stage to obtain the desired outcomes.
Collapse
Affiliation(s)
- Pallab Datta
- Centre for Healthcare Science and Technology, Indian Institute of Engineering Science and Technology Shibpur, Howrah 711103, West Bengal, India
| | - Ananya Barui
- Centre for Healthcare Science and Technology, Indian Institute of Engineering Science and Technology Shibpur, Howrah 711103, West Bengal, India
| | - Yang Wu
- Engineering Science and Mechanics Department, Penn State University, University Park, PA 16802, USA; The Huck Institutes of the Life Sciences, Penn State University, University Park, PA 16802, USA
| | - Veli Ozbolat
- Engineering Science and Mechanics Department, Penn State University, University Park, PA 16802, USA; The Huck Institutes of the Life Sciences, Penn State University, University Park, PA 16802, USA; Ceyhan Engineering Faculty, Cukurova University, Adana 01950, Turkey
| | - Kazim K Moncal
- Engineering Science and Mechanics Department, Penn State University, University Park, PA 16802, USA; The Huck Institutes of the Life Sciences, Penn State University, University Park, PA 16802, USA
| | - Ibrahim T Ozbolat
- Engineering Science and Mechanics Department, Penn State University, University Park, PA 16802, USA; The Huck Institutes of the Life Sciences, Penn State University, University Park, PA 16802, USA; Biomedical Engineering Department, Penn State University, University Park, PA 16802, USA; Materials Research Institute, Penn State University, University Park, PA 16802, USA.
| |
Collapse
|