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.     

Multiscale brain modelling

A central difficulty of brain modelling is to span the range of spatio-temporal scales from synapses to the whole brain. This paper overviews results from a recent model of the generation of brain electrical activity that incorporates both basic microscopic neurophysiology and large-scale brain anatomy to predict brain electrical activity at scales from a few tenths of a millimetre to the whole brain. This model incorporates synaptic and dendritic dynamics, nonlinearity of the firing response, axonal propagation and corticocortical and corticothalamic pathways. Its relatively few parameters measure quantities such as synaptic strengths, corticothalamic delays, synaptic and dendritic time constants, and axonal ranges, and are all constrained by independent physiological measurements. It reproduces quantitative forms of electroencephalograms seen in various states of arousal, evoked response potentials, coherence functions, seizure dynamics and other phenomena. Fitting model predictions to experimental data enables underlying physiological parameters to be inferred, giving a new non-invasive window into brain function that complements slower, but finer-resolution, techniques such as fMRI. Because the parameters measure physiological quantities relating to multiple scales, and probe deep structures such as the thalamus, this will permit the testing of a range of hypotheses about vigilance, cognition, drug action and brain function. In addition, referencing to a standardized database of subjects adds strength and specificity to characterizations obtained.     

A multiphase model describing vascular tumour growth

In this paper we present a new model framework for studying vascular tumour growth, in which the blood vessel density is explicitly considered. Our continuum model comprises conservation of mass and momentum equations for the volume fractions of tumour cells, extracellular material and blood vessels. We include the physical mechanisms that we believe to be dominant, namely birth and death of tumour cells, supply and removal of extracellular fluid via the blood and lymph drainage vessels, angiogenesis and blood vessel occlusion. We suppose that the tumour cells move in order to relieve the increase in mechanical stress caused by their proliferation. We show how to reduce the model to a system of coupled partial differential equations for the volume fraction of tumour cells and blood vessels and the phase averaged velocity of the mixture. We consider possible parameter regimes of the resulting model. We solve the equations numerically in these cases, and discuss the resulting behaviour. The model is able to reproduce tumour structure that is found in vivo in certain cases. Our framework can be easily modified to incorporate the effect of other phases, or to include the effect of drugs.     

Multiscale simulation of transmembrane proteins

Multiscale simulation is employed to examine changes in atomistic-level protein structure due to long wavelength membrane undulations and plane stress fields. An ensemble of atomistic-level simulations of a model of a transmembrane influenza A virus M2 proton channel in a dimyristoylphosphatidylcholine (DMPC) bilayer is coupled to a corresponding mesoscopic model of a DMPC bilayer in an explicit mesoscopic solvent. Structural variations in the key proton gating His37 residues of the M2 channel are examined. Small, but distinct variations in the structure of the His37 residues are observed in both the open and closed states of the channel as a result of the coupling to mesoscopic-level membrane motions.     

一个良性肿瘤细胞生长的计算机仿真模型

本文建立了一个良性肿瘤在正常细胞组织中生长的计算机仿真模型。初始细胞模型采用二维Voronoi结构,使用繁特卡罗法将其离散成x^*y个小细胞,对肿瘤与正常细胞分别施加蒙特卡罗工关过程使细胞生长。根据细胞组织内营养水平可以决定肿瘤细胞（分别为休眠细胞和有生长繁殖能力的细胞）的生长状态。     

Multiscale simulation in polymer science

The paper gives a short overview on recent approaches to link several time and length scales in soft matter simulations. Special attention is given to the fact that in contrast to low molecular weight compounds, intramolecular entropy is as important as the energetic contribution. First applications to industrially relevant problems are mentioned.     

Recent progress in modelling and simulation of three-dimensional tumor growth and treatment

This paper illustrates how to apply methods of systems analysis, control theory and simulation to the field of biology and medicine. For this purpose normal and abnormal cell growth has been modelled at different levels. It was possible to simulate three-dimensional tumor growth and different kinds of treatment. The paper shows how tumor treatment may be optimized in the long run using computer simulation experiments as a powerful new tool prior to clinical therapy.     

Multiscale simulation of actin filaments and actin-associated proteins

Actin is an important cytoskeletal protein that serves as a building block to form filament networks that span across the cell. These networks are orchestrated by a myriad of other cytoskeletal entities including the unbranched filament–forming protein formin and branched network–forming protein complex Arp2/3. Computational models have been able to provide insights into many important structural transitions that are involved in forming these networks, and into the nature of interactions essential for actin filament formation and for regulating the behavior of actin-associated proteins. In this review, we summarize a subset of such models that focus on the atomistic features and those that can integrate atomistic features into a larger picture in a multiscale fashion.     

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.     

Mathematical modelling of microenvironment and growth in EMT6/Ro multicellular tumour spheroids

In order to determine the role of micromilieu in tumour spheroid growth, a mathematical model was developed to predict EMT6/Ro spheroid growth and microenvironment based upon numerical solution of the diffusion/reaction equation for oxygen, glucose, lactate ion, carbon dioxide, bicarbonate ion, chlorine ion and hydrogen ion along with the equation of electroneutrality. This model takes into account the effects of oxygen concentration, glucose concentration and extracellular pH on cell growth and metabolism. Since independent measurements of EMT6/Ro single cell growth and metabolic rates, spheroid diffusion constants, and spinner flask mass transfer coefficients are available, model predictions using these parameters were compared with published data on EMT6/Ro spheroid growth and micro-environment. The model predictions of reduced spheroid growth due to reduced cell growth rates and cell shedding fit experimental spheroid growth data below 700 microns, but overestimated the spheroid growth rate at larger diameters. Predicted viable rim thicknesses based on predicted near zero glucose concentrations fit published viable rim thickness data for 1000 microns spheroids grown at medium glucose concentrations of 5.5 mM or less. However, the model did not accurately predict the onset of necrosis. Moreover, the model could not predict the observed decreases in oxygen and glucose metabolism seen in spheroids with time, nor could it predict the observed growth plateau. This suggests that other unknown factors, such as inhibitors or cell-cell contact effects, must also be important in affecting spheroid growth and cellular metabolism.     

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     

Multiscale simulation of protein mediated membrane remodeling

Proteins interacting with membranes can result in substantial membrane deformations and curvatures. This effect is known in its broadest terms as membrane remodeling. This review article will survey current multiscale simulation methodologies that have been employed to examine protein mediated membrane remodeling.     

Cellular automaton simulation of tumour growth -- equivocal relationships between simulation parameters and morphologic pattern features.

OBJECTIVE: To develop an interpretation procedure which estimates simulation parameters (tumour cell motility, tumour cell adhesion, autocrine and paracrine growth control, stroma destruction) of simulated patterns solely based on morphometric features of the morphologic pattern. METHODS: A cellular automaton computer simulation program was developed which produces morphologic patterns by growth of a seed of tumour cells. At the beginning of each simulation run certain simulation parameters are assigned to the tumour cells. After the run has been completed, the resulting pattern is evaluated by a set of morphometric features. Simulation parameters and resulting morphometric features of 27,800 simulations were stored in a database and were used for the evaluation of potential relationships. RESULTS: Correlation analysis showed highly significant correlations between morphometric features on the one hand and the preset simulation parameters (tumour cell motility, tumour cell adhesion, autocrine and paracrine growth control, stroma destruction) on the other. Correlation coefficients, however, varied from 0.72 to 0.99. When only one simulation parameter varied while all others were kept constant, morphometric features yielded a highly reliable estimate of the particular simulation parameter. When variability was extended to 4 simulation parameters, morphometric features were less effective in estimating the setting of the parameters. Though in all patterns tested several possible simulation parameter constellations could be ruled out, morphometric features were usually compatible with more than one set of simulation parameters thus preventing a straightforward interpretation. CONCLUSIONS: Though simulation parameters significantly and reproducibly influence the resulting morphologic pattern as characterized by morphometric features, estimates of the simulation parameters based on morphometric features yield equivocal results.     

The mathematical modelling of tumour angiogenesis and invasion

In order to accomplish the transition from avascular to vascular growth, solid tumours secrete a diffusible substance known as tumour angiogenesis factor (TAF) into the surrounding tissue. Endothelial cells which form the lining of neighbouring blood vessels respond to this chemotactic stimulus in a well-ordered sequence of events comprising, at minimum, of a degradation of their basement membrane, migration and proliferation. Capillary sprouts are formed which migrate towards the tumour eventually penetrating it and permitting vascular growth to take place. It is during this stage of growth that the insidious process of invasion of surrounding tissues can and does take place. A model mechanism for angiogenesis is presented which includes the diffusion of the TAF into the surrounding host tissue and the response of the endothelial cells to the chemotactic stimulus. Numerical simulations of the model are shown to compare very well with experimental observations. The subsequent vascular growth of the tumour is discussed with regard to a classical reaction-diffusion pre-pattern model.