Abnormal morphology biases hematocrit distribution in tumor vasculature and contributes to heterogeneity in tissue oxygenation.

Oxygen heterogeneity in solid tumors is recognized as a limiting factor for therapeutic efficacy. This heterogeneity arises from the abnormal vascular structure of the tumor, but the precise mechanisms linking abnormal structure and compromised oxygen transport are only partially understood. In this paper, we investigate the role that red blood cell (RBC) transport plays in establishing oxygen heterogeneity in tumor tissue. We focus on heterogeneity driven by network effects, which are challenging to observe experimentally due to the reduced fields of view typically considered. Motivated by our findings of abnormal vascular patterns linked to deviations from current RBC transport theory, we calculated average vessel lengths [Formula: see text] and diameters [Formula: see text] from tumor allografts of three cancer cell lines and observed a substantial reduction in the ratio [Formula: see text] compared to physiological conditions. Mathematical modeling reveals that small values of the ratio λ (i.e., [Formula: see text]) can bias hematocrit distribution in tumor vascular networks and drive heterogeneous oxygenation of tumor tissue. Finally, we show an increase in the value of λ in tumor vascular networks following treatment with the antiangiogenic cancer agent DC101. Based on our findings, we propose λ as an effective way of monitoring the efficacy of antiangiogenic agents and as a proxy measure of perfusion and oxygenation in tumor tissue undergoing antiangiogenic treatment.

Chaste: Cancer, Heart and Soft Tissue Environment.

Chaste (Cancer, Heart And Soft Tissue Environment) is an open source simulation package for the numerical solution of mathematical models arising in physiology and biology. To date, Chaste development has been driven primarily by applications that include continuum modelling of cardiac electrophysiology ('Cardiac Chaste'), discrete cell-based modelling of soft tissues ('Cell-based Chaste'), and modelling of ventilation in lungs ('Lung Chaste'). Cardiac Chaste addresses the need for a high-performance, generic, and verified simulation framework for cardiac electrophysiology that is freely available to the scientific community. Cardiac chaste provides a software package capable of realistic heart simulations that is efficient, rigorously tested, and runs on HPC platforms. Cell-based Chaste addresses the need for efficient and verified implementations of cell-based modelling frameworks, providing a set of extensible tools for simulating biological tissues. Computational modelling, along with live imaging techniques, plays an important role in understanding the processes of tissue growth and repair. A wide range of cell-based modelling frameworks have been developed that have each been successfully applied in a range of biological applications. Cell-based Chaste includes implementations of the cellular automaton model, the cellular Potts model, cell-centre models with cell representations as overlapping spheres or Voronoi tessellations, and the vertex model. Lung Chaste addresses the need for a novel, generic and efficient lung modelling software package that is both tested and verified. It aims to couple biophysically-detailed models of airway mechanics with organ-scale ventilation models in a package that is freely available to the scientific community. Chaste is designed to be modular and extensible, providing libraries for common scientific computing infrastructure such as linear algebra operations, finite element meshes, and ordinary and partial differential equation solvers. This infrastructure is used by libraries for specific applications, such as continuum mechanics, cardiac models, and cell-based models. The software engineering techniques used to develop Chaste are intended to ensure code quality, re-usability and reliability. Primary applications of the software include cardiac and respiratory physiology, cancer and developmental biology.

Comparison of computational modelling techniques for braided stent analysis.

Braided stents are associated with a number of complications in vivo. Accurate computational modelling of these devices is essential for the design and development of the next generation of these stents. In this study, two commonly utilised methods of computationally modelling filament interaction in braided stents are investigated: the join method and the weave method. Three different braided stent designs are experimentally tested and computationally modelled in both radial and v-block configurations. The results of the study indicate that while both methods are capable of capturing braided stent performance to some degree, the weave method is much more robust.

Computational modelling of magnesium stent mechanical performance in a remodelling artery: Effects of multiple remodelling stimuli.

Significant research has been conducted in the area of coronary stents/scaffolds made from resorbable metallic and polymeric biomaterials. These next-generation bioabsorbable stents have the potential to completely revolutionise the treatment of coronary artery disease. The primary advantage of resorbable devices over permanent stents is their temporary presence which, from a theoretical point of view, means only a healed coronary artery will be left behind following degradation of the stent potentially eliminating long-term clinical problems associated with permanent stents. The healing of the artery following coronary stent/scaffold implantation is crucial for the long-term safety of these devices. Computational modelling can be used to evaluate the performance of complex stent devices in silico and assist in the design and development and understanding of the next-generation resorbable stents. What is lacking in computational modelling literature is the representation of the active response of the arterial tissue in the weeks and months following stent implantation, ie, neointimal remodelling, in particular for the case of biodegradable stents. In this paper, a computational modelling framework is developed, which accounts for two major physiological stimuli responsible for neointimal remodelling and combined with a magnesium corrosion model that is capable of simulating localised pitting (realistic) stent corrosion. The framework is used to simulate different neointimal growth patterns and to explore the effects the neointimal remodelling has on the mechanical performance (scaffolding support) of the bioabsorbable magnesium stent.

The importance of geometry in the corneal micropocket angiogenesis assay.

The corneal micropocket angiogenesis assay is an experimental protocol for studying vessel network formation, or neovascularization, in vivo. The assay is attractive due to the ease with which the developing vessel network can be observed in the same animal over time. Measurements from the assay have been used in combination with mathematical modeling to gain insights into the mechanisms of angiogenesis. While previous modeling studies have adopted planar domains to represent the assay, the hemispherical shape of the cornea and asymmetric positioning of the angiogenic source can be seen to affect vascular patterning in experimental images. As such, we aim to better understand: i) how the geometry of the assay influences vessel network formation and ii) how to relate observations from planar domains to those in the hemispherical cornea. To do so, we develop a three-dimensional, off-lattice mathematical model of neovascularization in the cornea, using a spatially resolved representation of the assay for the first time. Relative to the detailed model, we predict that the adoption of planar geometries has a noticeable impact on vascular patterning, leading to increased vessel 'merging', or anastomosis, in particular when circular geometries are adopted. Significant differences in the dynamics of diffusible aniogenesis simulators are also predicted between different domains. In terms of comparing predictions across domains, the 'distance of the vascular front to the limbus' metric is found to have low sensitivity to domain choice, while metrics such as densities of tip cells and vessels and 'vascularized fraction' are sensitive to domain choice. Given the widespread adoption and attractive simplicity of planar tissue domains, both in silico and in vitro, the differences identified in the present study should prove useful in relating the results of previous and future theoretical studies of neovascularization to in vivo observations in the cornea.

Microvessel Chaste: An Open Library for Spatial Modeling of Vascularized Tissues.

Spatial models of vascularized tissues are widely used in computational physiology. We introduce a software library for composing multiscale, multiphysics models for applications including tumor growth, angiogenesis, osteogenesis, coronary perfusion, and oxygen delivery. Composition of such models is time consuming, with many researchers writing custom software. Recent advances in imaging have produced detailed three-dimensional (3D) datasets of vascularized tissues at the scale of individual cells. To fully exploit such data there is an increasing need for software that allows user-friendly composition of efficient, 3D models of vascularized tissues, and comparison of predictions with in vivo or in vitro experiments and alternative computational formulations. Microvessel Chaste can be used to build simulations of vessel growth and adaptation in response to mechanical and chemical stimuli; intra- and extravascular transport of nutrients, growth factors and drugs; and cell proliferation in complex 3D geometries. In addition, it can be used to develop custom software for integrating modeling with experimental data processing workflows, facilitated by a comprehensive Python interface to solvers implemented in C++. This article links to two reproducible example problems, showing how the library can be used to build simulations of tumor growth and angiogenesis with realistic vessel networks.

Predicting the Influence of Microvascular Structure On Tumor Response to Radiotherapy.

OBJECTIVE: The purpose of this study is to investigate how theoretical predictions of tumor response to radiotherapy (RT) depend on the morphology and spatial representation of the microvascular network. METHODS: A hybrid multiscale model, which couples a cellular automaton model of tumor growth with a model for oxygen transport from blood vessels, is used to predict the viable fraction of cells following one week of simulated RT. Both artificial and biologically derived three-dimensional (3-D) vessel networks of well vascularized tumors are considered and predictions compared with 2-D descriptions. RESULTS: For literature-derived values of the cellular oxygen consumption rate there is little difference in predicted viable fraction when 3-D network representations of biological or artificial vessel networks are employed. Different 2-D representations are shown to either over- or under-estimate viable fractions relative to the 3-D cases, with predictions based on point-wise descriptions shown to have greater sensitivity to vessel network morphology. CONCLUSION: The predicted RT response is relatively insensitive to the morphology of the microvessel network when 3-D representations are adopted, however, sensitivity is greater in certain 2-D representations. SIGNIFICANCE: By using realistic 3-D vessel network geometries this study shows that real and artificial network descriptions and assumptions of spatially uniform oxygen distributions lead to similar RT response predictions in relatively small tissue volumes. This suggests that either a more detailed description of oxygen transport in the microvasculature is required or that the oxygen enhancement ratio used in the well known linear-quadratic RT response model is relatively insensitive to microvascular structure.

Optimizing the design of a bioabsorbable metal stent using computer simulation methods.

Computer simulation is used extensively in the design of permanent stents. In order to address new challenges that arise in the design of absorbable metal stents (AMSs), such as corrosion and the limited mechanical properties of bioabsorbable alloys, new simulation and design techniques are needed. In this study a new method for simulating AMS corrosion is developed to study the effects of corrosion on the mechanical performance of a range of stent designs. The corrosion model is combined with an optimization strategy to identify AMS features that give optimal corrosion performance in the body. It is found that strut width is the predominant geometrical factor in determining long-term AMS scaffolding performance. An AMS with superior scaffolding performance to a commercial design is identified, based on deployment and corrosion simulations in stenosed vessels. These simulation and design techniques give new insights into in-vivo AMS performance and the role of device geometry in determining long-term scaffolding performance.