PubMed Health⌕ Search

Biomedical subjects

M A Chaplain

Publications and source records attributed to M A Chaplain.

At least 19 recordsLinked to original sources
PubMed2001

Spatio-temporal pattern formation on spherical surfaces: numerical simulation and application to solid tumour growth.

In this paper we examine spatio-temporal pattern formation in reaction-diffusion systems on the surface of the unit sphere in 3D. We first generalise the usual linear stability analysis for a two-chemical system to this geometrical context. Noting the limitations of this approach (in terms of rigorous prediction of spatially heterogeneous steady-states) leads us to develop, as an alternative, a novel numerical method which can be applied to systems of any dimension with any reaction kinetics. This numerical method is based on the method of lines with spherical harmonics and uses fast Fourier transforms to expedite the computation of the reaction kinetics. Numerical experiments show that this method efficiently computes the evolution of spatial patterns and yields numerical results which coincide with those predicted by linear stability analysis when the latter is known. Using these tools, we then investigate the rĵle that pre-pattern (Turing) theory may play in the growth and development of solid tumours. The theoretical steady-state distributions of two chemicals (one a growth activating factor, the other a growth inhibitory factor) are compared with the experimentally and clinically observed spatial heterogeneity of cancer cells in small, solid spherical tumours such as multicell spheroids and carcinomas. Moreover, we suggest a number of chemicals which are known to be produced by tumour cells (autocrine growth factors), and are also known to interact with one another, as possible growth promoting and growth inhibiting factors respectively. In order to connect more concretely the numerical method to this application, we compute spatially heterogeneous patterns on the surface of a growing spherical tumour, modelled as a moving-boundary problem. The numerical results strongly support the theoretical expectations in this case. Finally in an appendix we give a brief analysis of the numerical method.

Animals↗
PubMed2001

A new mathematical model for avascular tumour growth.

The early development of solid tumours has been extensively studied, both experimentally via the multicellular spheroid assay, and theoretically using mathematical modelling. The vast majority of previous models apply specifically to multicell spheroids, which have a characteristic structure of a proliferating rim and a necrotic core, separated by a band of quiescent cells. Many previous models represent these as discrete layers, separated by moving boundaries. Here, the authors develop a new model, formulated in terms of continuum densities of proliferating, quiescent and necrotic cells, together with a generic nutrient/growth factor. The model is oriented towards an in vivo rather than in vitro setting, and crucially allows for nutrient supply from underlying tissue, which will arise in the two-dimensional setting of a tumour growing within an epithelium. In addition, the model involves a new representation of cell movement, which reflects contact inhibition of migration. Model solutions are able to reproduce the classic three layer structure familiar from multicellular spheroids, but also show that new behaviour can occur as a result of the nutrient supply from underlying tissue. The authors analyse these different solution types by approximate solution of the travelling wave equations, enabling a detailed classification of wave front solutions.

Cell Division↗
PubMed2000

Analytical solutions of a minimal model of species migration in a bounded domain.

A minimal model of species migration is presented which takes the form of a parabolic equation with boundary conditions and initial data. Solutions to the differential problem are obtained that can be used to describe the small- and large-time evolution of a species distribution within a bounded domain. These expressions are compared with the results of numerical simulations and are found to be satisfactory within appropriate temporal regimes. The solutions presented can be used to describe existing observations of nematode distributions, can be used as the basis for further work on nematode migration, and may also be interpreted more generally.

Animals↗
PubMed1998

Continuous and discrete mathematical models of tumor-induced angiogenesis.

Angiogenesis, the formation of blood vessels from a pre-existing vasculature, is a process whereby capillary sprouts are formed in response to externally supplied chemical stimuli. The sprouts then grow and develop, driven initially by endothelial-cell migration, and organize themselves into a dendritic structure. Subsequent cell proliferation near the sprout tip permits further extension of the capillary and ultimately completes the process. Angiogenesis occurs during embryogenesis, wound healing, arthritis and during the growth of solid tumors. In this paper we present both continuous and discrete mathematical models which describe the formation of the capillary sprout network in response to chemical stimuli (tumor angiogenic factors, TAF) supplied by a solid tumor. The models also take into account essential endothelial cell-extracellular matrix interactions via the inclusion of the matrix macromolecule fibronectin. The continuous model consists of a system of nonlinear partial differential equations describing the initial migratory response of endothelial cells to the TAF and the fibronectin. Numerical simulations of the system, using parameter values based on experimental data, are presented and compared qualitatively with in vivo experiments. We then use a discretized form of the partial differential equations to develop a biased random-walk model which enables us to track individual endothelial cells at the sprout tips and incorporate anastomosis, mitosis and branching explicitly into the model. The theoretical capillary networks generated by computer simulations of the discrete model are compared with the morphology of capillary networks observed in in vivo experiments.

Angiogenesis Inducing Agents↗
PubMed1997

Two-dimensional models of tumour angiogenesis and anti-angiogenesis strategies.

There is a very strong link between the vascularization of a tumour and the spread of the disease, both locally and to distant sites (Gimbrone et al., 1974, J. Natl. Cancer Inst. 52, 413-27; Muthukkaruppan et al., 1982, J. Natl. Cancer Inst. 69, 699-704; Ellis & Fiddler, 1995, Lancet 346, 388-9). A tumour becomes vascularized by a process known as angiogenesis. Tumour angiogenesis is initiated by the release of diffusible substances by the tumour, whereby neighbouring capillary vessels are stimulated to grow and eventually penetrate the tumour. Anti-angiogenesis has been proposed as a potential strategy for the treatment of cancer (Folkman, 1995, Nature Med. 1, 21-31; Harris et al., 1996, Breast Cancer Res. Treat. 38, 97-108). In this paper, a mathematical model of the development of the tumour vasculature is presented. By suitable manipulation of the model parameters, we simulate various anti-angiogenesis strategies and we examine the roles that haptotaxis and chemotaxis may play during the growth of the neovasculature. The model is simulated in two space dimensions (on a square domain) so that it is, in theory, experimentally reproducible and any predictions of the model can therefore be tested.

Angiogenesis Inducing Agents↗
PubMed1996

On the rôle of angiogenesis in wound healing.

Angiogenesis, the formation of blood vessels, may be described as a process whereby capillary sprouts are formed in response to externally supplied chemical stimuli. The sprouts then develop and organize themselves into a dendritic structure. Angiogenesis occurs during embryogenesis, wound healing, arthritis and during the growth of solid tumours. In this paper we present a mathematical model which describes the rŏle of angiogenesis as observed during (soft-tissue) wound healing. We focus attention on certain principal players involved in this complex process, namely capillary tips, capillary sprouts, fibroblasts, macrophage-derived chemical attractants, oxygen and extracellular matrix. The model consists of a system of nonlinear partial differential equations describing the interactions in space and time of the above substances. Numerical simulations are presented which are in very good qualitative agreement with experimental observations.

Animals↗
PubMed1996

A mathematical model of the first steps of tumour-related angiogenesis: capillary sprout formation and secondary branching.

The growth of a solid tumour is dependent on an adequate supply of nutrients. A tumour can establish a blood supply by inducing neighbouring blood vessels to sprout and grow towards it, a process known as angiogenesis. The tumour cells may secrete a number of diffusible chemicals which stimulate endothelial cell to migrate, to rearrange themselves into capillary tubes or sprouts, and to proliferate. In this paper we focus firstly upon the early stage of angiogenesis wherein the endothelial cells group together in the parent vessel to form the initial capillary-sprout buds. A mathematical model for the formation of the capillary buds is presented which focuses on the potential role that haptotaxis may play. In Section 2 we turn attention to the endothelial cells within the growing and developing capillary sprouts as they migrate towards the tumour cells. Once again the potential role of haptotaxis is focused upon.

Animals↗
PubMed1996

The role of N-methyl-D-aspartate (NMDA) receptors in wind-up: a mathematical model.

We present a mathematical model for the phenomenon of wind-up (Mendell, 1966, Exper. Neur. 16,316-22) which occurs in many neurons. We concentrate on its occurrence in the substantia gelatinosa of the dorsal horns of the spinal cord, where it is connected with certain pathological and nonpathological pain states. The model is a development of the model by Britton & Skevington (1989, J. Theor. Biol. 137, 91-105) for Melzack & Wall's gate control theory of pain (1965, Science, New York, 150, 971-9; 1982, The Challenge of Pain, Penguin: Harmondsworth), modified to take account of more recent information. Its variables are the electric potentials of various cells in the midbrain and the spinal cord. Britton & Skevington's original model simulated many of the phenomena observed in acute pain in humans, but not the wind-up mechanism. This is not surprising, since this model did not include the N-methyl-D-aspartate (NMDA) receptors that are now recognized as being crucial to the phenomenon. Here we rectify this omission, and obtain good agreement between the model and experimental data on wind-up. The positive feedback that NMDA receptors exhibit is shown to be the essential feature in producing wind-up. As an independent test of the model we simulate a completely different experimental set-up, and obtain good qualitative agreement with data there. Finally, we present a prediction of the model that has yet to be tested experimentally.

Animals↗
PubMed1996

Mathematical modelling, simulation and prediction of tumour-induced angiogenesis.

Angiogenesis, the formation of blood vessels from a pre-existing vasculature, is a process whereby capillary sprouts are formed in response to externally supplied chemical stimuli. The sprouts then grow and develop, driven by endothelial cell migration and proliferation, and organise themselves into a dendritic structure. Angiogenesis occurs during embryogenesis, wound healing, arthritis and during the growth of solid tumours. In this paper we present a novel mathematical model which describes the formation of the capillary sprout network in response to chemical stimuli (tumour angiogenesis factors, TAF) supplied by a solid tumour. The model also takes into account endothelial cell-extracellular matrix interactions via the inclusion of fibronectin in the model. The model consists of a system of nonlinear partial differential equations describing the response in space and time of endothelial cells to the TAF and the fibronectin (migration, proliferation, anastomosis, branching). Using the discretized system of partial differential equations, we use a deterministic cellular automata (DCA) model, which enables us to track individual endothelial cells and incorporate branching explicity into the model. Numerical simulations are presented which are in very good qualitative agreement with experimental observations. Certain experiments are suggested which could be used to test the hypotheses of the model and various extensions and developments of the model with particular applications to anti-angiogenesis strategies are discussed.

Animals↗
PubMed1995

A mathematical analysis of a model for tumour angiogenesis.

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. Neighbouring endothelial cells 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. A mathematical model is presented which takes into account two of the most important events associated with the endothelial cells as they form capillary sprouts and make their way towards the tumour i.e. cell migration and proliferation. The numerical simulations of the model compare very well with the actual experimental observations. We subsequently investigate the model analytically by making some relevant biological simplifications. The mathematical analysis helps to clarify the particular contributions to the model of the two independent processes of endothelial cell migration and proliferation.

Angiogenesis Inducing Agents↗
PubMed1995

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.

Angiogenesis Inducing Agents↗
PubMed1995

Mathematical models for tumour angiogenesis: numerical simulations and nonlinear wave solutions.

To ensure its sustained growth, a tumour may secrete chemical compounds which cause neighbouring capillaries to form sprouts which then migrate towards it, furnishing the tumour with an increased supply of nutrients. In this paper a mathematical model is presented which describes the migration of capillary sprouts in response to a chemoattractant field set up by a tumour-released angiogenic factor, sometimes termed a tumour angiogenesis factor (TAF). The resulting model admits travelling wave solutions which correspond either to successful neovascularization of the tumour or failure of the tumour to secure a vascular network, and which exhibit many of the characteristic features of angiogenesis. For example, the increasing speed of the vascular front, and the evolution of an increasingly developed vascular network behind the leading capillary tip front (the brush-border effect) are both discernible from the numerical simulations. Through the development and analysis of a simplified caricature model, valuable insight is gained into how the balance between chemotaxis, tip proliferation and tip death affects the tumour's ability to induce a vascular response from neighbouring blood vessels. In particular, it is possible to define the success of angiogenesis in terms of known parameters, thereby providing a potential framework for assessing the viability of tumour neovascularization in terms of measurable quantities.

Animals↗
PubMed1995

Growth of nonnecrotic tumors in the presence and absence of inhibitors.

In this article a model for the evolution of a spherically symmetric, nonnecrotic tumor is presented. The effects of nutrients and inhibitors on the existence and stability of time-independent solutions are studied. With a single nutrient and no inhibitors present, the trivial solution, which corresponds to a state in which no tumor is present, persists for all parameter values, whereas the nontrivial solution, which corresponds to a tumor of finite size, exists for only a prescribed range of parameters, which corresponds to a balance between cell proliferation and cell death. Stability analysis, based on a two-timing method, suggests that, where it exists, the nontrivial solution is stable and the trivial solution unstable. Otherwise, the trivial solution is stable. Modification to these characteristic states brought about by the presence of different types of inhibitors are also investigated and shown to have significant effect. Implications of the model for the treatment of cancer are also discussed.

Cell Division↗
PubMed1994

Nonlinear diffusion of a growth inhibitory factor in multicell spheroids.

A mathematical model is presented for the production of a growth inhibitory factor (GIF) within a multicell spheroid. The model is based on the assumption that the GIF diffuses within the spheroid in a nonlinear spatially dependent manner. This is in contrast with previous models, in which the nonlinearity was assumed in the production term. The results of the new model are compared with those of previous models and with experimental data.

Animals↗
PubMed1993

Modelling the growth of solid tumours and incorporating a method for their classification using nonlinear elasticity theory.

Medically, tumours are classified into two important classes--benign and malignant. Generally speaking, the two classes display different behaviour with regard to their rate and manner of growth and subsequent possible spread. In this paper, we formulate a new approach to tumour growth using results and techniques from nonlinear elasticity theory. A mathematical model is given for the growth of a solid tumour using membrane and thick-shell theory. A central feature of the model is the characterisation of the material composition of the model through the use of a strain-energy function, thus permitting a mathematical description of the degree of differentiation of the tumour explicitly in the model. Conditions are given in terms of the strain-energy function for the processes of invasion and metastasis occurring in a tumour, being interpreted as the bifurcation modes of the spherical shell which the tumour is essentially modelled as. Our results are compared with actual experimental results and with the general behaviour shown by benign and malignant tumours. Finally, we use these results in conjunction with aspects of surface morphogenesis of tumours (in particular, the Gaussian and mean curvatures of the surface of a solid tumour) in an attempt to produce a mathematical formulation and description of the important medical processes of staging and grading cancers. We hope that this approach may form the basis of a practical application.

Cell Division↗
PubMed1993

A qualitative analysis of some models of tissue growth.

Using maximum principles for parabolic and elliptic operators, we examine, in a general way, some models of tissue growth. These typically consist of a model mechanism for the diffusion of a mitotic inhibitor (growth inhibitory factor, GIF) throughout the tissue. Central to the modeling is the inclusion of a source function that models the production of GIF throughout the tissue. We examine the effect this term has on the resulting distribution of GIF in the tissue and comment on the appropriateness of different source functions, in particular a uniform production rate or a nonuniform production rate of inhibitor. Given that it is more appropriate to infer from the patterns of mitosis that are observed experimentally in various tissues the GIF concentration profile rather than the source function profile, it may be more appropriate to use these types of models to determine the qualitative form of the source term rather than proposing this function a priori.

Animals↗
PubMed1993

On the concentration profile of a growth inhibitory factor in multicell spheroids.

A mathematical model is presented for the production of a growth inhibitory factor (GIF) within a multicell spheroid. The main assumption of the model is that the GIF is produced by cells within the spheroid in some prescribed nonlinear, spatially dependent manner. Given that the diffusion of the GIF is known to take place over a much shorter time scale than that of spheroid growth, the steady-state profile of the GIF in various spheroids of differing radii is examined and theoretical results are compared with actual experimental data.

Animals↗