Multi-scale model of lumen formation via inverse membrane blebbing mechanism during sprouting angiogenesis process.

Cancer is one of the leading causes of mortality and morbidity among people worldwide. Cancer appears as solid tumors in many cases. Angiogenesis is the growth of blood vessels from the existing vasculature and is one of the imperative processes in tumor growth. Another vital phenomenon for formation and functionality of this vasculature network is lumen formation. The results of recent studies indicate the importance of blood pressure in this mechanism. Computational modeling can study these processes in different scales. Hence, wide varieties of these models have been proposed during recent years. In this research, a multi-scale model is developed for the angiogenesis process. In the extracellular scale, the growth factor concentration is calculated via the reaction diffusion equation. At the cellular scale, growth, migration, and the adhesion of endothelial cells are modeled by the Potts cellular model. At the intra-cellular scale by considering biochemical signals, a Boolean network model describes migration, division, or apoptosis of endothelial cells. A stochastic model developed for lumen formation via inverse membrane blebbing mechanism. A CFD simulation was also used to investigate the role of pulsated blood pressure in the inverse membrane blebbing mechanism. The lumen formation model shows stochastic behavior in blebs expansion and lumen expansion. Comparing the stochastic model's results with the CFD simulation also shows the vital role of pressure pulse and the topology of the blebs in bleb retraction.

Correction: The Vital Role of Blood Flow-Induced Proliferation and Migration in Capillary Network Formation in a Multiscale Model of Angiogenesis.

[This corrects the article DOI: 10.1371/journal.pone.0128878.].

Spatiotemporal distribution modeling of PET tracer uptake in solid tumors.

OBJECTIVE: Distribution of PET tracer uptake is elaborately modeled via a general equation used for solute transport modeling. This model can be used to incorporate various transport parameters of a solid tumor such as hydraulic conductivity of the microvessel wall, transvascular permeability as well as interstitial space parameters. This is especially significant because tracer delivery and drug delivery to solid tumors are determined by similar underlying tumor transport phenomena, and quantifying the former can enable enhanced prediction of the latter. METHODS: We focused on the commonly utilized FDG PET tracer. First, based on a mathematical model of angiogenesis, the capillary network of a solid tumor and normal tissues around it were generated. The coupling mathematical method, which simultaneously solves for blood flow in the capillary network as well as fluid flow in the interstitium, is used to calculate pressure and velocity distributions. Subsequently, a comprehensive spatiotemporal distribution model (SDM) is applied to accurately model distribution of PET tracer uptake, specifically FDG in this work, within solid tumors. RESULTS: The different transport mechanisms, namely convention and diffusion from vessel to tissue and in tissue, are elaborately calculated across the domain of interest and effect of each parameter on tracer distribution is investigated. The results show the convection terms to have negligible effect on tracer transport and the SDM can be solved after eliminating these terms. CONCLUSION: The proposed framework of spatiotemporal modeling for PET tracers can be utilized to comprehensively assess the impact of various parameters on the spatiotemporal distribution of PET tracers.

Blood flow and endothelial cell phenotype regulation during sprouting angiogenesis.

The role of the endothelial cell environment and shear stress induced by blood flow in phenotype determination and lumen formation has been clearly illustrated in recent studies. In the present work, a model is developed to map environmental and flow induced signals in sprouting angiogenesis to endothelial cell phenotype and lumen formation. To follow the endothelial cell lumen formation, its signaling pathway is incorporated in the present work within the phenotype determination pathway that has been recently utilized to model endothelial cell migration, proliferation, and apoptosis. Moreover, a signaling cascade for shear stress activation of endothelial cells is proposed and used for phenotype determination with activation of blood flow. A Boolean network model is employed to build a hybrid map for the relation between the endothelial cell environmental signals and the endothelial cell fate in sprouting angiogenesis with and without blood flow. This map is very useful in the development of models for sprouting angiogenesis. Moreover, this study shows that inhibition of intracellular signaling molecules, solely or in pairs, blocks angiogenic-signaling pathways and can be used to inhibit angiogenesis.

Activation of Apoptotic Signal in Endothelial Cells through Intracellular Signaling Molecules Blockade in Tumor-Induced Angiogenesis.

Tumor-induced angiogenesis is the bridge between avascular and vascular tumor growth phases. In tumor-induced angiogenesis, endothelial cells start to migrate and proliferate toward the tumor and build new capillaries toward the tumor. There are two stages for sprout extension during angiogenesis. The first stage is prior to anastomosis, when single sprouts extend. The second stage is after anastomosis when closed flow pathways or loops are formed and blood flows in the closed loops. Prior to anastomosis, biochemical and biomechanical signals from extracellular matrix regulate endothelial cell phenotype; however, after anastomosis, blood flow is the main regulator of endothelial cell phenotype. In this study, the critical signaling pathways of each stage are introduced. A Boolean network model is used to map environmental and flow induced signals to endothelial cell phenotype (proliferation, migration, apoptosis, and lumen formation). Using the Boolean network model, blockade of intracellular signaling molecules of endothelial cell is investigated prior to and after anastomosis and the cell fate is obtained in each case. Activation of apoptotic signal in endothelial cell can prevent the extension of new vessels and may inhibit angiogenesis. It is shown that blockade of a few signaling molecules in endothelial cell activates apoptotic signal that are proposed as antiangiogenic strategies.

The Vital Role of Blood Flow-Induced Proliferation and Migration in Capillary Network Formation in a Multiscale Model of Angiogenesis.

Sprouting angiogenesis and capillary network formation are tissue scale phenomena. There are also sub-scale phenomena involved in angiogenesis including at the cellular and intracellular (molecular) scales. In this work, a multiscale model of angiogenesis spanning intracellular, cellular, and tissue scales is developed in detail. The key events that are considered at the tissue scale are formation of closed flow path (that is called loop in this article) and blood flow initiation in the loop. At the cellular scale, growth, migration, and anastomosis of endothelial cells (ECs) are important. At the intracellular scale, cell phenotype determination as well as alteration due to blood flow is included, having pivotal roles in the model. The main feature of the model is to obtain the physical behavior of a closed loop at the tissue scale, relying on the events at the cellular and intracellular scales, and not by imposing physical behavior upon it. Results show that, when blood flow is considered in the loop, the anastomosed sprouts stabilize and elongate. By contrast, when the loop is modeled without consideration of blood flow, the loop collapses. The results obtained in this work show that proper determination of EC phenotype is the key for its survival.

Effect of fluid friction on interstitial fluid flow coupled with blood flow through solid tumor microvascular network.

A solid tumor is investigated as porous media for fluid flow simulation. Most of the studies use Darcy model for porous media. In Darcy model, the fluid friction is neglected and a few simplified assumptions are implemented. In this study, the effect of these assumptions is studied by considering Brinkman model. A multiscale mathematical method which calculates fluid flow to a solid tumor is used in this study to investigate how neglecting fluid friction affects the solid tumor simulation. The mathematical method involves processes such as blood flow through vessels and solute and fluid diffusion, convective transport in extracellular matrix, and extravasation from blood vessels. The sprouting angiogenesis model is used for generating capillary network and then fluid flow governing equations are implemented to calculate blood flow through the tumor-induced capillary network. Finally, the two models of porous media are used for modeling fluid flow in normal and tumor tissues in three different shapes of tumors. Simulations of interstitial fluid transport in a solid tumor demonstrate that the simplifications used in Darcy model affect the interstitial velocity and Brinkman model predicts a lower value for interstitial velocity than the values that Darcy model predicts.

Effect of tumor shape, size, and tissue transport properties on drug delivery to solid tumors.

BACKGROUND: The computational methods provide condition for investigation related to the process of drug delivery, such as convection and diffusion of drug in extracellular matrices, drug extravasation from microvessels or to lymphatic vessels. The information of this process clarifies the mechanisms of drug delivery from the injection site to absorption by a solid tumor. In this study, an advanced numerical method is used to solve fluid flow and solute transport equations simultaneously to investigate the effect of tumor shape and size on drug delivery to solid tumor. METHODS: The advanced mathematical model used in our previous work is further developed by adding solute transport equation to the governing equations. After applying appropriate boundary and initial conditions on tumor and surrounding tissue geometry, the element-based finite volume method is used for solving governing equations of drug delivery in solid tumor. Also, the effects of size and shape of tumor and some of tissue transport parameters such as effective pressure and hydraulic conductivity on interstitial fluid flow and drug delivery are investigated. RESULTS: Sensitivity analysis shows that drug delivery in prolate shape is significantly better than other tumor shapes. Considering size effect, increasing tumor size decreases drug concentration in interstitial fluid. This study shows that dependency of drug concentration in interstitial fluid to osmotic and intravascular pressure is negligible. CONCLUSIONS: This study shows that among diffusion and convection mechanisms of drug transport, diffusion is dominant in most different tumor shapes and sizes. In tumors in which the convection has considerable effect, the drug concentration is larger than that of other tumors at the same time post injection.