Multiscale modelling and nonlinear simulation of vascular tumour growth

In this article, we present a new multiscale mathematical model for solid tumour growth which couples an improved model of tumour invasion with a model of tumour-induced angiogenesis. We perform nonlinear simulations of the multi-scale model that demonstrate the importance of the coupling between the development and remodeling of the vascular network, the blood flow through the network and the tumour progression. Consistent with clinical observations, the hydrostatic stress generated by tumour cell proliferation shuts down large portions of the vascular network dramatically affecting the flow, the subsequent network remodeling, the delivery of nutrients to the tumour and the subsequent tumour progression. In addition, extracellular matrix degradation by tumour cells is seen to have a dramatic affect on both the development of the vascular network and the growth response of the tumour. In particular, the newly developing vessels tend to encapsulate, rather than penetrate, the tumour and are thus less effective in delivering nutrients.     

A history of the study of solid tumour growth: The contribution of mathematical modelling

A miscellany of new strategies, experimental techniques and theoretical approaches are emerging in the ongoing battle against cancer. Nevertheless, as new, ground-breaking discoveries relating to many and diverse areas of cancer research are made, scientists often have recourse to mathematical modelling in order to elucidate and interpret these experimental findings. Indeed, experimentalists and clinicians alike are becoming increasingly aware of the possibilities afforded by mathematical modelling, recognising that current medical techniques and experimental approaches are often unable to distinguish between various possible mechanisms underlying important aspects of tumour development. This short treatise presents a concise history of the study of solid tumour growth, illustrating the development of mathematical approaches from the early decades of the twentieth century to the present time. Most importantly these mathematical investigations are interwoven with the associated experimental work, showing the crucial relationship between experimental and theoretical approaches, which together have moulded our understanding of tumour growth and contributed to current anti-cancer treatments. Thus, a selection of mathematical publications, including the influential theoretical studies by Burton, Greenspan, Liotta et al., McElwain and co-workers, Adam and Maggelakis, and Byrne and co-workers are juxtaposed with the seminal experimental findings of Gray et al. on oxygenation and radio-sensitivity, Folkman on angiogenesis, Dorie et al. on cell migration and a wide variety of other crucial discoveries. In this way the development of this field of research through the interactions of these different approaches is illuminated, demonstrating the origins of our current understanding of the disease.     

Multiscale modelling of vascular tumour growth in 3D: the roles of domain size and boundary conditions

We investigate a three-dimensional multiscale model of vascular tumour growth, which couples blood flow, angiogenesis, vascular remodelling, nutrient/growth factor transport, movement of, and interactions between, normal and tumour cells, and nutrient-dependent cell cycle dynamics within each cell. In particular, we determine how the domain size, aspect ratio and initial vascular network influence the tumour's growth dynamics and its long-time composition. We establish whether it is possible to extrapolate simulation results obtained for small domains to larger ones, by constructing a large simulation domain from a number of identical subdomains, each subsystem initially comprising two parallel parent vessels, with associated cells and diffusible substances. We find that the subsystem is not representative of the full domain and conclude that, for this initial vessel geometry, interactions between adjacent subsystems contribute to the overall growth dynamics. We then show that extrapolation of results from a small subdomain to a larger domain can only be made if the subdomain is sufficiently large and is initialised with a sufficiently complex vascular network. Motivated by these results, we perform simulations to investigate the tumour's response to therapy and show that the probability of tumour elimination in a larger domain can be extrapolated from simulation results on a smaller domain. Finally, we demonstrate how our model may be combined with experimental data, to predict the spatio-temporal evolution of a vascular tumour.     

Towards whole-organ modelling of tumour growth

Multiscale approaches to modelling biological phenomena are growing rapidly. We present here some recent results on the formulation of a theoretical framework which can be developed into a fully integrative model for cancer growth. The model takes account of vascular adaptation and cell-cycle dynamics. We explore the effects of spatial inhomogeneity induced by the blood flow through the vascular network and of the possible effects of p27 on the cell cycle. We show how the model may be used to investigate the efficiency of drug-delivery protocols.     

Coupled modelling of tumour angiogenesis, tumour growth and blood perfusion

We propose a mathematical modelling system to investigate the dynamic process of tumour cell proliferation, death and tumour angiogenesis by fully coupling the vessel growth, tumour growth and blood perfusion. Tumour growth and angiogenesis are coupled by the chemical microenvironment and the cell-matrix interaction. The haemodynamic calculation is carried out on the updated vasculature. The domains of intravascular, transcapillary and interstitial fluid flow were coupled in the model to provide a comprehensive solution of blood perfusion variables. An estimation of vessel collapse is made according to the wall shear stress criterion to provide feedback on vasculature remodelling. The simulation can show the process of tumour angiogenesis and the spatial distribution of tumour cells for periods of up to 24 days. It can show the major features of tumour and tumour microvasculature during the period such as the formation of a large necrotic core in the tumour centre with few functional vessels passing through, and a well circulated tumour periphery regions in which the microvascular density is high and associated with more aggressive proliferating cells of the growing tumour which are all consistent with physiological observations. The study also demonstrated that the simulation results are not dependent on the initial tumour and networks, which further confirms the application of the coupled model feedback mechanisms. The model enables us to examine the interactions between angiogenesis and tumour growth, and to study the dynamic response of a solid tumour to the changes in the microenvironment. This simulation framework can be a foundation for further applications such as drug delivery and anti-angiogenic therapies.     

Mathematical modelling of the Warburg effect in tumour cords

The model proposed here links together two approaches to describe tumours: a continuous medium to describe the movement and the mechanical properties of the tissue, and a population dynamics approach to represent internal genetic inhomogeneity and instability of the tumour. In this way one can build models which cover several stages of tumour progression. In this paper we focus on describing transition from aerobic to purely glycolytic metabolism (the Warburg effect) in tumour cords. From the mathematical point of view this model leads to a free boundary problem where domains in contact are characterized by different sets of equations. Accurate stitching of the solution was possible with a modified ghost fluid method. Growth and death of the cells and uptake of the nutrients are related through ATP production and energy costs of the cellular processes. In the framework of the bi-population model this allowed to keep the number of model parameters relatively small.     

Forecasting the growth of multicell tumour spheroids: implications for the dynamic growth of solid tumours

The growth dynamics of multicell tumour spheroids (MTS) were analysed by means of mathematical techniques derived from signal processing theory. Volume vs. time trajectories of individual spheroids were fitted with the Gompertz growth equation and the residuals (i.e. experimental volume determinations minus calculated values by fitting) were analysed by fast fourier transform and power spectrum. Residuals were not randomly distributed around calculated growth trajectories demonstrating that the Gompertz model partially approximates the growth kinetics of three-dimensional tumour cell aggregates. Power spectra decreased with increasing frequency following a 1/f(delta) power-law. Our findings suggest the existence of a source of 'internal' variability driving the time-evolution of MTS growth. Based on these observations, a new stochastic Gompertzian-like mathematical model was developed which allowed us to forecast the growth of MTS. In this model, white noise is additively superimposed to the trend described by the Gompertz growth equation and integrated to mimic the observed intrinsic variability of MTS growth. A correlation was found between the intensity of the added noise and the particular upper limit of volume size reached by each spheroid within two MTS populations obtained with two different cell lines. The dynamic forces generating the growth variability of three-dimensional tumour cell aggregates also determine the fate of spheroid growth with a strong predictive significance. These findings suggest a new approach to measure tumour growth potential.     

Deterministic material-based averaging theory model of collagen gel micromechanics

Mechanics of collagen gels, like that of many tissues, is governed by events occurring on a length scale much smaller than the functional scale of the material. To deal with the challenge of incorporating deterministic micromechanics into a continuous macroscopic model, we have developed an averaging-theory-based modeling framework for collagen gels. The averaging volume, which is constructed around each integration point in a macroscopic finite-element model, is assumed to experience boundary deformations homogeneous with the macroscopic deformation field, and a micromechanical problem is solved to determine the average stress at the integration point. A two-dimensional version was implemented with the microstructure modeled as a network of nonlinear springs, and 500 segments were found to be sufficient to achieve statistical homogeneity. The method was then used to simulate the experiments of Tower et al. (Ann. Biomed. Eng., 30, pp. 1221-1233) who performed uniaxial extension of prealigned collagen gels. The simulation captured many qualitative features of the experiments, including a toe region and the realignment of the fibril network during extension. Finally, the method was applied to an idealized wound model based on the characterization measurements of Bowes et al. (Wound Repair Regen., 7, pp. 179-186). The model consisted of a strongly aligned "wound" region surrounded by a less strongly aligned "healthy" region. The alignment of the fibrils in the wound region led to reduced axial strains, and the alignment of the fibrils in the healthy region, combined with the greater effective stiffness of the wound region, caused rotation of the wound region during uniaxial stretch. Although the microscopic model in this study was relatively crude, the multiscale framework is general and could be employed in conjunction with any microstructural model.     

Multiscale stochastic modelling of gene expression

Stochastic phenomena in gene regulatory networks can be modelled by the chemical master equation for gene products such as mRNA and proteins. If some of these elements are present in significantly higher amounts than the rest, or if some of the reactions between these elements are substantially faster than others, it is often possible to reduce the master equation to a simpler problem using asymptotic methods. We present examples of such a procedure and analyse the relationship between the reduced models and the original.     

Spatial-temporal modelling of bacterial colony growth on solid media

The enormous diversity among bacterial colonies of different species, and among colonies of the same species under different environmental conditions, has long interested microbiologists. Yet it is only comparatively recently that quantitative, rather than merely conceptual, models have been developed to explain the dynamics of bacterial colony formation and growth. Understanding the fundamental processes that drive these dynamics is still at a rudimentary level, though a number of advances have been made. This review traces the history of bacterial colony growth modelling, from the pioneering work of Pirt in the late 1960s, through experimental investigations by Wimpenny and his colleagues in the 1970s, and further models extending from that work to understand complex bacterial colony formations. It concludes with recent results which find that both diameter and height of the colony follow simple power-law behaviour over the entire active growth period (from a few hours to several days old), and that the results of Pirt, Wimpenny and their contemporaries can be re-interpreted.     

The influence of tumour vasculature on fluid flow in solid tumours:a mathematical modelling study

Tumour vasculature is known to be aberrant,tortuous and erratic which can have significant implications for fluid flow.Fluid dynamics in tumour tissue plays an ...     

An evolutionary hybrid cellular automaton model of solid tumour growth

We propose a cellular automaton model of solid tumour growth, in which each cell is equipped with a micro-environment response network. This network is modelled using a feed-forward artificial neural network, that takes environmental variables as an input and from these determines the cellular behaviour as the output. The response of the network is determined by connection weights and thresholds in the network, which are subject to mutations when the cells divide. As both available space and nutrients are limited resources for the tumour, this gives rise to clonal evolution where only the fittest cells survive. Using this approach we have investigated the impact of the tissue oxygen concentration on the growth and evolutionary dynamics of the tumour. The results show that the oxygen concentration affects the selection pressure, cell population diversity and morphology of the tumour. A low oxygen concentration in the tissue gives rise to a tumour with a fingered morphology that contains aggressive phenotypes with a small apoptotic potential, while a high oxygen concentration in the tissue gives rise to a tumour with a round morphology containing less evolved phenotypes. The tissue oxygen concentration thus affects the tumour at both the morphological level and on the phenotype level.     

The role of acidity in solid tumour growth and invasion

Acidic pH is a common characteristic of human tumours. It has a significant impact on tumour progression and response to therapies. In this paper, we develop a simple model of three-dimensional tumour growth to examine the role of acidosis in the interaction between normal and tumour cell populations. Both vascular and avascular tumour dynamics are investigated, and a number of different behaviours are observed. Whilst an avascular tumour always proceeds to a benign steady state, a vascular tumour may display either benign or invasive dynamics, depending on the value of a critical parameter. Analysis of the model allows us to assess novel therapies directed towards changing the level of acidity within the tumour.     

Multiscale model of fatigue of collagen gels

Biomechanics and Modeling in Mechanobiology - Fatigue as a mode of failure becomes increasingly relevant with age in tissues that experience repeated fluctuations in loading. While there has been a...     

Multiscale modelling of material failure using particles

In this paper, we explore a prototypical model of material fragmentation – the penetration of a plate by a fast moving ball in two dimensions using both smooth particle applied mechanics (SPAM) and molecular dynamics (MD). Using a family of short-ranged pair potentials, we investigate the fracture as a function of impact velocity and plate strength. The MD provides data from which we make a direct comparison with SPAM simulations of the same system. The parameters for the SPAM simulations: equation of state, bulk and shear moduli, yield strength and tensile strength are obtained from the pair potentials and via non-equilibrium MD (tension test). By careful parameterisation of the SPAM model, we are able to observe qualitatively similar behaviour to that shown in the MD simulations of a single crystal (brittle) specimen; bending at low penetrator speed, cracking at higher speeds and eventually, complete penetration. The behaviour depends on the plate strength, which is directly related to details of the force law. For plates characterised as high strength, the agreement between MD and SPAM is evidently poor. We discuss the possible causes for this discrepancy. Nevertheless, SPAM, when used in conjunction with accurate atomistic simulation, is shown to be a highly promising tool for studying failure in materials.     

Mathematical modelling of tumour acidity

Acid-mediated tumour invasion is receiving increasing experimental and clinical attention. Previous models proposed to describe this phenomenon failed to capture key properties of the system, such as the existence of the benign steady state, or predicted incorrectly the size of the inter-tissue gap. Here we show that taking proper account of quiescence ameliorates these drawbacks as well as revealing novel behaviour. The simplicity of the model allows us to fully identify the key parameters controlling different aspects of behaviour.     

A weakly nonlinear analysis of a model of avascular solid tumour growth

 In this paper we study a mathematical model that describes the growth of an avascular solid tumour. Our analysis concentrates on the stability of steady, radially-symmetric model solutions with respect to perturbations taken from the class of spherical harmonics. Using weakly nonlinear analysis, previous results are extended to show how the amplitudes of the asymmetric modes interact. Attention focuses on a special case for which the model equations simplify. Analysis of the simplified model equations leads to the identification of a two-parameter family of asymmetric steady solutions, the dimensions of whose stable and unstable manifolds depend on the system parameters. The asymmetric steady solutions limit the basin of attraction of the radially-symmetric steady state when it is linearly stable. On the basis of these numerical and analytical results we postulate the existence of fully nonlinear steady solutions which are stable with respect to time-dependent perturbations. Received: 25 October 1998 / Revised version: 20 June 1998