Elsevier

Journal of Theoretical Biology

Volume 300, 7 May 2012, Pages 309-316
Journal of Theoretical Biology

Wound healing angiogenesis: The clinical implications of a simple mathematical model

https://doi.org/10.1016/j.jtbi.2012.01.043Get rights and content

Abstract

Nonhealing wounds are a major burden for health care systems worldwide. In addition, a patient who suffers from this type of wound usually has a reduced quality of life. While the wound healing process is undoubtedly complex, in this paper we develop a deterministic mathematical model, formulated as a system of partial differential equations, that focusses on an important aspect of successful healing: oxygen supply to the wound bed by a combination of diffusion from the surrounding unwounded tissue and delivery from newly formed blood vessels. While the model equations can be solved numerically, the emphasis here is on the use of asymptotic methods to establish conditions under which new blood vessel growth can be initiated and wound-bed angiogenesis can progress. These conditions are given in terms of key model parameters including the rate of oxygen supply and its rate of consumption in the wound. We use our model to discuss the clinical use of treatments such as hyperbaric oxygen therapy, wound bed debridement, and revascularisation therapy that have the potential to initiate healing in chronic, stalled wounds.

Highlights

► In this paper, a deterministic mathematical model of wound healing is developed. ► We consider oxygen supply, blood vessels and newly formed capillary buds. ► Conditions are established for successful wound healing. ► We discuss the use of treatments that can initiate healing in chronic, stalled wounds.

Introduction

Lazarus et al. (1994) define a chronic wound to be one which either fails to proceed through an orderly and timely process to produce anatomic and functional integrity, or proceeds through the repair process without establishing a sustained anatomic and functional result. A recent estimate suggests that the US health care system spends in excess of US$25 billion annually treating patients with nonhealing wounds (Sen et al., 2009). Fig. 1 shows typical data illustrating how wound closure varies between healing and nonhealing wounds (Roy et al., 2009). It should be noted that these data are from porcine wounds, such animal wounds being thought accurately to reflect equivalent human wound data (Sullivan et al., 2001).

Wound healing is a highly regulated and complex process, consisting of four stages (haemostasis, inflammation, proliferation and remodelling using the terminology of Grey and Harding, 2006) and requiring the coordination of the activities of many chemical and cellular species. Haemostasis typically lasts for a few hours and involves the control of blood loss in the damaged region. The inflammation stage lasts several days and coincides with inflammatory cell migration into the wound space and the release of chemical factors such as vascular endothelial growth factor (VEGF). These chemicals provide the stimulus that ultimately leads to the formation of new blood vessels (angiogenesis), an important step in the proliferative stage of healing. During this healing phase, there is a surge in the proliferation rate of fibroblasts, endothelial and epithelial cells and the rate at which collagen is deposited by fibroblasts (Jeffcoate et al., 2004). The final stage of healing sees the wound increase in tensile strength via remodelling of the extracellular matrix. The healing process is tightly regulated by many factors including oxygen supply and new capillary development. Grey and Harding (2006) provide a recent review of human wound healing and the factors that modulate it.

Most mathematical models of wound healing can be categorised as either population-based (or continuum) or cell-based (or discrete). An advantage of adopting a discrete approach is that it is possible to incorporate details that cannot easily be included within a continuum framework (see for example, Anderson and Chaplain, 1998, Cumming et al., 2010). These features include cell–cell interactions, individual cell cycles, positioning of daughter cells after proliferation and discrete collagen fibres with individual orientations. However, continuum models are often more amenable to analysis than discrete ones. Geris et al. (2010) review in silico treatment strategies for wound healing together with some of the mathematical models that have been used to describe the healing process. Mathematical models that closely relate to experimental and/or clinical data have been the focus of recent modelling research applied to angiogenesis and vasculature development (Machado et al., 2011, Aubert et al., 2011).

One of the earliest continuum models of angiogenesis in wound healing is a six-species, partial differential equation (PDE) model (chemoattractant, tips, oxygen blood vessels, fibroblasts and extracellular matrix) developed by Pettet et al. (1996b) to simulate the ingrowth of blood vessels. By performing numerical simulations they identified parameter sets for which healing stalled (by looking at mechanisms such as new capillary tip sprouting and chemoattractant production). In other work, Pettet et al. (1996a) used perturbation methods to derive approximate expressions for the wave speed of a soft-tissue healing and Byrne et al. (2000) used the same model to compare to experimental data. More recently, Schugart et al. (2008) developed a seven-species model of angiogenesis (the six species from the model by Pettet et al., 1996b and additionally macrophages) in order to investigate the role of oxygen tension in cutaneous wound healing. Schugart et al. used their model to generate several predictions. For example, they claim that wounds will not heal in extremely hypoxic environments and that the use of hyperbaric oxygen therapy may stimulate angiogenesis. These predictions are consistent with earlier experimental work (Hopf et al., 2005, Kim et al., 2007). Xue et al. (2009) extended the model of Schugart et al. (2008) by incorporating mechanical effects in that they treated the ECM as a viscoelastic material. To the best of our knowledge, Xue et al. (2009) is the first mechanochemical model of wound healing angiogenesis. Vermolen (2009) has developed a system of nonlinear reaction-diffusion equations for oxygen, growth factors, epidermal cells and capillaries and simulated healing in two spatial dimensions. While these models give considerable insight into wound healing, they do not lend themselves to the derivation of simple expressions relating the success (or failure) of wound healing to the system parameters.

There have been several discrete models of the related phenomenon of tumour-induced angiogenesis, including important work by Chaplain and Anderson (1999) that is based on a finite difference approximation of PDEs and Bauer et al. (2007) who use the cellular Potts model framework. It is worth noting that wound healing angiogenesis is a regulated process whereas tumour-induced angiogenesis is uncontrolled (Byrne and Chaplain, 1996). For a detailed review of mathematical models of tumour-induced angiogenesis, see Mantzaris et al. (2004).

In this paper, we develop a mathematical model based on the assumption that revascularisation of the wound bed is the rate-limiting step in successful healing. While other processes, such as ECM deposition and remodelling, are undoubtedly important, restoration of a good oxygen supply is vital for the repair of damaged tissue (Hunt and Gimbel, 2002). Our aim is to derive simple criteria, in terms of the model parameters, for which successful healing will be initiated. We then use these criteria to assess common treatment therapies such as debridement, revascularistion, and hyperbaric oxygen therapy (Thackham et al., 2008, Flegg et al., 2009).

In the next section we develop our mathematical model and present the governing partial differential equations. In Section 3, we present typical numerical simulations while in Section 4 we analyse the model and identify regions of parameter space in which we predict whether healing will succeed or fail. In Section 5, we discuss the implications of our results and make suggestions for further work.

Section snippets

Description of the mathematical model

Our mathematical model comprises three partial differential equations: one for the oxygen concentration, w, one for the capillary tip density, n, and one for the blood vessel density, b. The equations are based on the principle of mass conservation and are stated below in dimensional form. For simplicity, we consider a one-dimensional wound whose edge is located at x=0 and whose centre lies at x=L, with symmetry about x=L. It is worth noting that in our one-dimensional model, the capillary tip

Numerical results

Eqs. (3a), (3b), (3c), (4a), (4b) are solved numerically using a finite volume method, with a Roe-flux limiting approach employed to discretise the chemotaxis term in Eq. (3b) (Thackham et al., 2009). Unless otherwise stated, the dimensionless parameters are fixed at the following values: Dw=10,χ=1,k2=150,k4=150,k5=100,k6=100,wH=0.5,wL=0.3,ε=0.1,L=2.

We stress that this choice of parameter values is employed for illustrative purposes only, since our main goal is to derive conditions, in terms of

Establishing necessary conditions for initiating healing

The early behaviour is known to be an accurate predictor of the ultimate success or failure of a wound to heal (Margolis et al., 2004). In this section we investigate whether the success or failure of healing can be predicted from the early time dynamics of our model. We can then use our results to predict those types of nonhealing wounds that would benefit from treatment. Guided by our numerical simulations we consider cases for which the parameters Dw, k2, k4, k5 and k6 are large. More

Clinical implications and discussion

Regardless of their aetiology, many chronic wounds are hypoxic and this can compromise healing (Mathieu, 2002). We can use our model to compare the efficacy of different treatments that are routinely used to treat nonhealing wounds. We note that if a wound has stalled due to a lack of oxygen (see point U1 in Fig. 4), it should be possible to initiate healing by either sufficiently increasing the oxygen supply rate, as characterised by the parameter k2 (the point S2) or sufficiently decreasing

Acknowledgements

This work was supported by the award of a doctoral scholarship to J.A.F. from the Institute of Health and Biomedical Innovation at Queensland University of Technology and was funded by Australian Research Council's Discovery Projects funding scheme (Project no. DP0878011). This research was carried out while H.M.B. was visiting Queensland University of Technology, funded by the Institute of Health and Biomedical Innovation and the Discipline of Mathematical Sciences. Computational resources and

References (37)

  • H. Byrne et al.

    Mathematical modelling of angiogenesis in wound healing: comparison of theory and experiment

    J. Theor. Med.

    (2000)
  • M. Chaplain et al.

    Modelling the Growth and Form of Capillary Networks

    (1999)
  • B. Cumming et al.

    A mathematical model of wound healing and subsequent scarring

    J. R. Soc. Interface

    (2010)
  • L. Edelstein et al.

    A model for fungal colony growth applied to Sclerotium reolfsii

    J. Gen. Microbiol.

    (1983)
  • J.A. Flegg et al.

    Mathematical model of hyperbaric oxygen therapy applied to chronic diabetic wounds

    Bull. Math. Biol.

    (2010)
  • J.A. Flegg et al.

    A three species model to simulate application of hyperbaric oxygen therapy to chronic wounds

    PLoS Comput. Biol.

    (2009)
  • E. Gaffney et al.

    Investigating a simple model of cutaneous wound healing angiogenesis

    J. Math. Biol.

    (2002)
  • L. Geris et al.

    In silico design of treatment strategies in wound healing and bone fracture healing

    Philos. Trans. R. Soc. A: Math. Phys. Eng. Sci.

    (2010)
  • Cited by (35)

    • A mathematical model of the use of supplemental oxygen to combat surgical site infection

      2019, Journal of Theoretical Biology
      Citation Excerpt :

      Meanwhile, others have investigated the potential therapeutic effects of oxygen. By modelling the impact of hyperbaric oxygen on chronic wounds and performing a parameter sensitivity analysis, they identifed parameter ranges for specific patient groups who may benefit from such treatment (Flegg et al., 2012; 2009). Conversely, there exist mathematical models of infection which do not consider the effect of oxygen levels.

    • A mathematical model for lymphangiogenesis in normal and diabetic wounds

      2015, Journal of Theoretical Biology
      Citation Excerpt :

      Some authors have used moving boundary methods to study the movement of the wound edge during healing (Javierre et al., 2009), and attempts have been made to understand more specific aspects of the healing process, for example macrophage dynamics in diabetic wounds (Waugh and Sherratt, 2006) and the resolution of the inflammatory phase (Dunster et al., 2014). One aspect of wound healing that has received considerable attention from mathematicians is (blood) angiogenesis (see, for instance, Byrne et al. (2000); Chaplain (2000); Flegg et al. (2012); Gaffney et al. (2002); Machado et al. (2011); Pettet et al. (1996); Schugart et al. (2008); Zheng et al. (2013)). A comprehensive review of mathematical models for vascular network formation can be found in Scianna et al.׳s (2013) review.

    View all citing articles on Scopus
    View full text