- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents
Computational and Mathematical Methods in Medicine
Volume 2013 (2013), Article ID 141087, 9 pages
Modelling Tumour Oxygenation, Reoxygenation and Implications on Treatment Outcome
1Medical Radiation Physics, Stockholm University and Karolinska Institutet, 171 76 Stockholm, Sweden
2Department of Radiation Physics UHL, County Council of Östergötland, Radiation Physics and Department of Medical and Health Sciences, Faculty of Health Sciences, Linköping University, 581 83 Linköping, Sweden
Received 11 December 2012; Accepted 26 December 2012
Academic Editor: Eva Bezak
Copyright © 2013 Iuliana Toma-Dasu and Alexandru Dasu. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Oxygenation is an important component of the tumour microenvironment, having a significant impact on the progression and management of cancer. Theoretical determination of tissue oxygenation through simulations of the oxygen transport process is a powerful tool to characterise the spatial distribution of oxygen on the microscopic scale and its dynamics and to study its impact on the response to radiation. Accurate modelling of tumour oxygenation must take into account important aspects that are specific to tumours, making the quantitative characterisation of oxygenation rather difficult. This paper aims to discuss the important aspects of modelling tumour oxygenation, reoxygenation, and implications for treatment.
Oxygenation is an important component of the tumour microenvironment and has a significant impact on the progression and management of cancer. Hypoxia, or reduced oxygenation, has been associated with restrained proliferation, apoptosis, necrosis, angiogenesis, the development of an aggressive phenotype, and resistance to treatment . These aspects make the quantification of tumour oxygenation and the study of its impact on treatment approaches highly relevant topics in oncology research. Consequently, several methods have been proposed for measuring tissue oxygenation, employing polarographic, fluorescent, radionuclide, or magnetic resonance techniques . The results of such measurements have generally been used for qualitative correlations with treatment outcome [3–9]. However, many of these methods are invasive or suffer from sampling or resolution limitations and these aspects may reduce the significance of the results. From this perspective, theoretical determination of tissue oxygenation through simulations of the oxygen transport process might be a practical alternative to characterise the spatial distribution of oxygen on the microscopic scale and its dynamics and to study its impact on the response to radiation. This paper aims to discuss the important aspects of modelling tumour oxygenation, reoxygenation, and implications for treatment.
2. Origins of Tumour Hypoxia
One of the most striking features of solid tumours in comparison to normal tissues is their architecture of vascular networks. In normal tissues, the vasculature has developed to provide adequate supply of oxygen to all of the cells, while in tumours, the vascular network is inadequate, almost chaotic and the quality of the blood flowing in it can be quite poor [10, 11]. This is the result of the generally high proliferation rates of tumour cells in comparison to normal cells that allow the tumours to quickly outgrow the existing blood and nutrient supply. As sustained angiogenesis is one of the hallmarks of cancer , oxygen supply is restored in tumours through the formation of new blood vessels. This is a relatively slow process  and newly formed vessels are also outgrown by the proliferating cells. Furthermore, newly formed vessels in tumours usually originate on the venous side of the vasculature . Newly formed vessels in most solid tumours are often dilated, tortuous and have incomplete or missing endothelial lining resulting in an increased vascular permeability . They are also unable to respond to local or systemic signals for vasodilatation or vasoconstriction . In addition, tumour vasculature is not accompanied by a lymphatic system hampering the drainage of interstitial fluid and by-products of cell metabolism which can further influence vascular flow. Indeed, changes in interstitial pressure are thought to be responsible for the temporary cessation or resumption of flow through tumour microvessels. Further, the preference of tumour cells to glucose metabolism leading to the formation of lactic acid as a by-product leads to an acidification of the extracellular fluid in tumours . All of these effects lead to a unique microenvironment in tumours characterised by poor supply and steep gradients of oxygen and other nutrients.
Depending on the underlying mechanisms and their duration, two main types of hypoxia have been identified in tumours, chronic and acute. Limitations in oxygen diffusion from the blood vessels into the tissue lead to the development of diffusion-limited or chronic hypoxia, first described by Thomlinson and Gray . In contrast, local disturbances in vessel perfusion lead to the appearance of perfusion-limited or acute hypoxia, described by Brown  and verified experimentally by Chaplin et al. . This has also been referred to as cyclic or fluctuating hypoxia due to the temporary character of the disturbances and the restoration of the perfusion through the affected vessels. Further subtypes of hypoxia have been identified depending on the mechanisms behind the development of hypoxic regions in tumours . An important feature of tumour oxygenation is that it is not static and that several mechanisms determine its dynamics on timescales ranging from a few minutes for acute hypoxia to several days for chronic hypoxia. Accurate modelling of tumour oxygenation would have to take into account all of these aspects. However, their underlying complexity makes the study of oxygen transport to tumours quite difficult, especially when the aim is to accurately quantify the oxygen distribution.
3. Intravascular Oxygen Transport
In vertebrates, oxygen transport from the lungs to tissues takes place through the circulatory system. Most oxygen is transported bound to the heme groups of the haemoglobin (Hb) molecule in red blood cells (RBCs). The binding of the oxygen to haemoglobin takes places in 4 steps corresponding to the 4 heme groups of the metalloprotein, since binding one oxygen molecule determines a conformational change of the Hb molecule that influences its ability to further bind other molecules. From this perspective, the kinetics of the oxygen-haemoglobin reaction is described by the Adair equation  giving the equilibrium relationship between the saturation of the Hb molecule and the oxygen tension. However, the Adair equation is seldom used for the mathematical modelling of oxygen transport since it cannot be inverted analytically to give the oxygen tension as a function of the saturation of the haemoglobin. Instead, the simpler Hill equation  is used: where is the fractional haemoglobin saturation, p is the local partial pressure of oxygen (or pO2), , is the pO2 value at which the haemoglobin is 50% saturated, and is the Hill parameter.
The sigmoid relationship described by the Adair and Hill equations is usually known as the oxyhaemoglobin dissociation curve or the oxygen-haemoglobin equilibrium curve. This curve is affected by a number of factors including temperature, pH, and CO2 concentration . Given the unique microenvironment of tumours with a low pH, the oxygen-haemoglobin equilibrium may be achieved at different oxygen tensions in tumours compared to normal tissues and there may also be differences between various normal tissues. In light of these expected differences, caution is advised at the generalisation of parameters determined for a particular tissue or during the extrapolation of the results achieved using these parameters.
Besides the bound component, modelled oxygen transport would also have to account for the free and dissolved oxygen in the red blood cells and blood plasma. Unbound oxygen represents only 1.5% of the total oxygen content of the blood in normal tissues, provided Hb is completely saturated . However, this fraction might be higher for tumours, where vessels often originate on the venous side of the vasculature containing desaturated blood. The association between anaemia and cancer would further decrease the relative importance of the bound component.
Modelling of intravascular transport has to account for both the free and the Hb-bound oxygen, as well as for the transport of oxygen through the membrane of the red blood cells. Popel  provided a comprehensive review of the earlier theoretical work on the particulate nature of RBC-plasma transport of oxygen. More recent approaches, however, have incorporated RBC-plasma transport into intravascular transport models accounting for intracapillary O2 gradients and giving the average flux of oxygen through the capillary wall . These models have mainly been used for studying the oxygen delivery to muscles and have found that blood-tissue transport encounters an effective intravascular resistance (IVR) that could be quantified as a mass transfer coefficient. This is influenced by the assumptions for intravascular transport, the shape of the RBC or vessel permeability. This suggests that care has to be employed when extrapolating normal tissue models to tumours, as the poor endothelial lining of tumour microvessels increases their permeability, while the acidic tumour microenvironment could influence the shape of the RBC as well as the Hb dissociation curve.
4. Tissue Transport Models
The first attempts to theoretically describe the oxygen transport into tissue date from the beginning of the 20th century [23, 24] and were based on systems with simple geometries for which an analytical solution could be obtained. One of these early models, the Krogh tissue cylinder model, has been particularly important as it has served as the foundation for many subsequent studies. Its main assumptions were the following.(i)pO2 distribution in tissue has a cylindrical symmetry and axial diffusion is not significant;(ii)tissue pO2 at capillary wall equals capillary pO2 (no intravascular resistance);(iii)oxygen transport in tissue takes place through passive diffusion;(iv)tissue diffusivity is independent of spatial position;(v)transport phenomena are steady state.
Many models derived from the Krogh model have focused on the study of oxygenation in normal tissues as reviewed by Popel  and Goldman . However, the assumptions of the Krogh tissue cylinder model have also been used by Thomlinson and Gray  and by Tannock  to describe theoretically the spatial distribution of tumour hypoxia. Later developments saw the introduction of models with more realistic vascular geometries to reflect the complexity of tumour vasculature, using Green functions or numerical methods to describe the tissue oxygenation [2, 26–31].
A general expression for the reaction-diffusion equation describing the oxygen transport in homogeneous tissues is shown in where is the local pO2, is the diffusion coefficient, is the flow velocity, q represents the local consumption, and is a source term.
In (2), the term describes the movement of oxygen through diffusion and the term describes the movement through convection. The convection term is usually neglected as the convective currents in the interstitial compartment of the normal tissues are estimated to represent only about 0.5–1% of the plasma flow. This may, however, not be the case of tumours, where the convective compartment could reach up to 15% . The impact of neglecting the convective compartment in modelling tumour oxygenation is yet to be studied.
The rate of oxygen consumption in tissues, q, is thought to follow the Michaelis-Menten kinetics (3) as originally proposed by Tang , although variations of the equation have also been considered : where is the maximum consumption rate, p is the local pO2, and is a parameter describing the pO2 at which the consumption rate halves. Equation (3) predicts an almost constant consumption rate at high oxygen tensions that is in agreement with the zero-order kinetics assumed by some of the earlier models [23, 25, 28, 34]. The Michaelis-Menten kinetics has, however, the advantage of more accurately describing the consumption at low pO2 and has been used in recent models [2, 27, 29, 31, 35–38].
The source term, , in (2) has been used in different ways in various models. Some models assume that the source term could be set to zero since cells do not produce oxygen [29, 37–40], while others implement a different approach where the oxygen supply from the capillaries is described as a distributed source throughout the tissue with localised spikes at the vessel positions [30, 40].
Most of the studies have focused on a steady state of the system where . This is in agreement with the observations by Hill , that the nonstationary evolution of the system is very short. Changes of oxygenation patterns due to perfusion limitations have time scales of minutes , while the formation of new blood vessels through angiogenesis takes place over several days  and therefore the steady state solution would apply in most situations. However, there might be some cases where the time-dependent oxygen transport would have to be considered, such as the onset of acute hypoxia after the cessation of blood flow, its disappearance following the reperfusion of capillaries, or cases when the oxygen supply through tumour capillaries has variations for longer time periods.
Solving (2) to determine tissue oxygenation requires setting boundary conditions for the limits of the modelled domain and even in this case various approaches have been used. For the interface between capillaries and tissue, one could neglect the intravascular resistance and use a Dirichlet boundary condition assuming that the pO2 value at the capillary wall equals the capillary pO2 [23, 29]. Alternatively, one could use a Robin boundary condition, assuming that the oxygen flux across the interface, that is, the vessel wall, must be continuous [37, 39, 40]. Assuming cylindrical symmetry for a capillary with radius and thickness w, the diffusive flux on the tissue side of the vessel wall could be written as where is the diffusivity in the wall and is the capillary pO2. Assuming that the wall thickness is much smaller than the external radius R, (4) becomes where is the vessel permeability.
The expression in (5) could be integrated over the surface of the capillary to obtain the net rate of oxygen entering the tissue. Thus, the source term in (2) could be written for each point representing the capillary as where is the capillary volume and is its area. Explicitly writing the volume and area of a cylinder, (6) becomes
As expected, different boundary conditions would lead to different solutions of the reaction-diffusion equation, especially in the regions close to the vessels. However, Skeldon et al.  showed that for vessels with very high permeability, the Robin boundary condition leads to similar results as the Dirichlet boundary condition. If this is the case of tumour capillaries with increased permeability, it could provide an advantage when solving (2), since Dirichlet boundary conditions are easier to program.
Similarly, different approaches may be used for the outer boundaries of the simulated domain. These include Neuman boundary conditions setting the flux through the outer boundaries to zero, imposing periodic boundary conditions or other approaches [23, 29, 37]. The choice of the border conditions for the outer boundaries could also influence the amount of calculated hypoxia in the tissue as shown by Secomb et al.  who discussed the advantages and disadvantages of these approaches.
5. Modelling Parameters
An important aspect for the theoretical modelling of tissue oxygenation is the choice of parameters in (1)–(7). Experimental determination of these parameters is rather challenging and in many cases it involves the assumption of a model relating a measured quantity to the quantity to be determined. This highlights the importance of using parameters in the limit of the models for which they have been determined.
The least debated parameter in the literature is the average diffusivity of oxygen into tissue, with most authors agreeing on D = 2 × 10−5 cm2 s−1 [25, 28]. This could be explained from the point of view of the passive transport of oxygen through tissue and the relative homogeneity of soft tissues in the body.
In contrast, a relatively large range of values have been proposed for the maximum consumption rate, , in (3) [10, 41–43], predominantly explained by the variability of the metabolic characteristics of tissues. As early as the beginning of the 20th century, Hill  observed significant differences in the consumption rates in active versus inactive muscle, while more recently Kallinowski et al.  showed that oxygen consumption in tumour cells decreases as cell quiescence develops. These are key observations because oxygen consumption is an important determinant of tissue oxygenation as shown by Secomb et al. , who found that a relatively modest reduction of tissue consumption may abolish tissue hypoxia, while a considerable increase of capillary oxygen content would be required for the same purpose. These results illustrate the difficulties in selecting the relevant consumption parameter for modelling. Indeed, significant differences may appear between tumour types or individuals depending on their metabolic activity, and intratumour heterogeneity could also be expected. Further confounding factors when determining the tissue consumption rate may reside in the method used for the determination and the assumptions made with respect to the oxygen supply rate or the invasiveness of the method used. Nevertheless, a relevant value for the maximum oxygen consumption rate in tumours that has been used in many studies is 15 mmHg s−1 [15, 25, 29].
Another parameter that is needed for the Robin boundary conditions (5) or for the distributed source models (7) is the permeability of tumour capillaries. Most of the permeability values available are relevant for normal tissue capillaries and these values have been used in studies investigating the oxygenation of tumours [30, 39, 40]. However, it has to be borne in mind that tumour vasculature has poor endothelial lining resulting in increased vascular permeability  and that normal tissue values may lead to an overestimation of the intravascular resistance of tumour capillaries.
Simulations of tumour oxygenation also require vascular oxygenation for the boundary conditions. Given the venous origin of newly formed vessels in tumours , a common assumption is that vascular pO2 in tumours is rather low [15, 25, 29, 37] and this is in agreement with experimental determinations of vascular oxygenations [45, 46]. Some studies have assumed a uniform oxygenation of the vessels [29, 30, 37, 39], while others have considered distributions of values [2, 28, 31]. Comparisons between simulations with the two assumptions have shown microscale differences between the resulting distributions, but no significant differences in the global characterisation of tissue oxygenation .
The vascular geometry is another important factor that has to be taken into account when simulating tumour oxygenation. Given the complexity of tumour vasculature, three-dimensional (3D) data would be required for the full characterisation of tumour oxygenation, but such images are difficult to obtain with existing in vivo imaging modalities due to their limited spatial resolution. Skin flap window chambers could be a solution to provide 3D geometries , but the thickness of the tissue in this model is much smaller than the chamber diameter. Alternatively, two-dimensional (2D) sections of ex vivo tissue samples could be used for the characterisation of vasculature and oxygenation. This approach was originally used by Thomlinson and Gray  for comparing the appearance of necrotic regions in lung tumours with calculations of tumour oxygenation and also in more recent studies employing tissue sections as input for simulations [28, 39]. An alternative approach is to use generated 2D maps of capillaries, either reflecting measured distributions of intervascular distances [2, 29, 31] or randomly generated distributions [30, 38, 40]. The former approach should be preferred as it has been shown that a full description of tissue vasculature is required for an accurate characterisation of tissue oxygenation. Average intervascular distance or the related quantity mean vascular density giving the number of vessels per unit area is one determinant of tissue oxygenation, the other being the shape of the distribution of intravascular distances .
Oxygen transport has also been coupled with the transport of other substances. Kirkpatrick et al.  used a mathematical model of oxygen and glucose mass transport in tumours to study the influence of kinetic and physical factors on metabolism. Other models are concerned with the transport and metabolism of tracers used in advanced imaging methods capable of visualising tumour hypoxia [30, 39, 40].
6. Modelling the Dynamics of Hypoxia
The dynamics of tumour hypoxia is an important aspect that has to be taken into consideration for many modelling studies. Indeed, angiogenesis and fluctuations in the perfusion of the blood vessels supplying the tumours have the potential to alter the oxygenation pattern, even in the absence of a therapeutic action upon the tumour cell population. These changes could in principle be modelled by assuming different spatial distributions of blood vessels over time, as the impaired perfusion of some vessels considered responsible for acute hypoxia would be equivalent to a decrease of the effective vascular density, while the appearance of new vessels would increase it. In this context it is important to note that even if the resulting distributions from repeated simulations are similar, in the case of acute hypoxia, different cells will be in different hypoxic subcompartments at each modelled time point and this dynamic process needs to be reflected in the modelled radiosensitivity of the cells.
Besides fluctuations in perfusion that are responsible for the cycling variation of the acutely hypoxic compartment, there are other mechanisms that may lead to changes in oxygenation during the course of a treatment. Radiation and other cytotoxic agents may induce growth inhibition in the irradiated cells that would subsequently lead to a decreased oxygen consumption in the surviving tumour cells  and thus improved oxygenation for previously hypoxic cells [35, 48]. On a longer timescale, the removal of the doomed cells near the blood vessels damaged by the cytotoxic agents might also improve the oxygenation of cells situated further away.
Another important aspect related to the timescale of hypoxia is the radiosensitivity of the affected cells. It is rather well known that decreased oxygen availability could lead to increased radioresistance of the hypoxic cells. Less known, however, is the effect of nutrient starvation that often accompanies chronic hypoxia, but not acute hypoxia. Indeed, it has been suggested [49, 50] that the two types of hypoxia might have different radiosensitivities due to changes in the repair capacity of the chronically hypoxic cells starved of oxygen and nutrients [51–55]. Consequently, differentiating between the two types of hypoxia as well as the time-dependent variation of the oxygen and nutrient availability might be needed for studies aiming to quantify the impact of certain treatment approaches on tumour containing both types of hypoxic cells. In this case, simulations would have to account for the perfusion status of the blood vessels in the tumour at various time points during the treatment and the impact on cellular radiosensitivity .
7. Tumour Oxygenation and Treatment Modelling
Theoretical simulations of oxygen transport to tumours have played a key role in understanding many aspects of tumour microenvironment that could influence treatment outcome. For example, theoretical modelling has allowed the characterisation of the gradients that may appear around the blood vessels, the study of the impact of the temporary closure of some capillaries, the radiobiological distinction between the two hypoxic compartments, or the effectiveness of various measurements methods. Comparisons between the maximum diffusion distance of oxygen into tissue and size of the viable rims of tumour cells around stroma or individual blood vessels led to the first descriptions of the appearance of chronic hypoxia in tumours [15, 25]. Simulations have also shown that several factors influence the relationship between tumour vasculature and oxygenation, indicating that there is an equivocal relationship between tissue oxygenation and vascular oxygenation, mean intravascular distance, or the number of unperfused vessels. Thus, mean vascular density is an important determinant of tumour oxygenation, but the relationship is modulated by the shape of the distribution of intervascular distances, the vascular oxygen content, and the number of closed capillaries [2, 29, 31]. This explains the failure to find in experimental studies a direct correlation between tumour oxygenation and vascular density .
Calculating the expected response from cell populations with different oxygenation requires a cell survival model with parameters that could be modified according to the radiation sensitivity of each compartment. According to the linear quadratic LQ model [58, 59], cell survival in a fully oxic population after a single radiation dose is given by where and are model parameters for oxic conditions. The modification of radiosensitivity in hypoxia could be accounted for through the use of modifying factors (OMFs) that are dependent on the local oxygen tension and the duration of hypoxia as in
A noteworthy expression for the oxygen modification factors has been proposed by Alper and Howard-Flanders : where is the maximum protection achieved in the absence of oxygen and is a reaction constant around 2.5–3 mmHg [61, 62].
Cell survival for populations of cells characterised by a distribution of oxygenations could be obtained as the weighted sum of cell survival for each compartment: where is the relative weight of each oxygenation compartment.
For fractionated treatments delivering fractions of size , different equations may be employed depending on the interfraction dynamics of oxygenation in the assumed population. Thus, for populations with static oxygenations, the cell survival is given by while for populations with full dynamics of the oxygenation, cell survival is given by
Theoretical simulations using (9)–(13) have shown significant differences in the predicted response depending on the assumed oxygenation of the tissue and these results indicate that full distributions of pO2 values, not only the hypoxic fraction, are needed for modelling [25, 56, 63]. This is particularly important in light of the equivocal relationship between tumour oxygenation and the above-mentioned parameters.
Equally interesting have been the results of simulations taking into account the possible radiobiological differences between acute and chronic hypoxic cells. Experimental studies have shown that chronically hypoxic cells may have low or depleted energy reserves which would impair their repair mechanisms [52–55], in contrast with acutely hypoxic cells that are repair competent. Correlating these findings with calculated oxygen distributions in tumours has suggested that the added biochemical sensitisation of chronically hypoxic cells could provide an explanation for the success of radiotherapy, as the chemical radioresistance conferred by hypoxia in general would otherwise require unrealistically high doses to control tumours containing hypoxic cells [49, 50, 56]. This comes as a confirmation of studies showing an improvement of control in experimental systems in the presence of chronically hypoxic cells . In contrast, the presence of acute hypoxia in tumours leads to a worsening of the response to radiation therapy [56, 65, 66]. However, while the fluctuating character of the perfusion disturbances in tumour capillaries does not lead to a full “washout” of the radioresistance of acutely hypoxic cells, the fluctuations are essential for the success of fractionated therapy. Indeed, treatment schedules employing too few fractions (5 or fewer) might not allow enough interfraction reoxygenation opportunities to change the sensitivity of the hypoxic cells in tumours which would in turn lead to poor outcome . This could have significant implications for stereotactic radiotherapy where extremely hypofractionated schedules should be pursued with caution .
Theoretical simulations have also been used to study the effectiveness of polarographic electrodes used to characterise tissue oxygenation [65, 69–72] or imaging with hypoxia-specific tracers [30, 39, 73–75]. It has thus been shown that averaging characteristic both to polarographic methods and to imaging methods with finite spatial resolution as well as differences in sensitivity of the measurement probes could lead to measured oxygenations that would be quite different from the tissue oxygenation [70, 73, 76], that in turn would translate into different predictions regarding tissue responsiveness to radiation therapy. This highlights the importance of carefully extrapolating measurements of tissue oxygenation to quantitative assessments of the expected response from radiation therapy.
In vivo imaging of hypoxia with dedicated tracers has received considerable attention in recent years as it could be used to quantify the spatial distribution and the severity of hypoxic areas in tumours for individualising or adapting the treatment according to biological factors [67, 77, 78]. In this context it has been shown that the dynamic character of hypoxia observed in clinical and experimental studies [79–84] has important implications for radiation therapy and it needs to be taken into account for treatment simulations. Indeed, neglecting it for prescriptions of highly heterogeneous dose distributions could lead to mismatches between planned dose hotspots and radioresistant hypoxic domains in tumours that would lead to poor outcome . Furthermore, the limited spatial resolution of imaging methods could lead to poor rendering of the small hypoxic regions  that would translate into an underestimation of their radiosensitivity. However, this might be less of a problem for rapid interfraction reoxygenation that might ensure the success of dose painting techniques that do not take into account the microscale heterogeneity of hypoxia .
Besides these direct applications, theoretical modelling has proven a valuable tool in supplementing experimental data and for understanding the effects of factors that could not be studied individually or experimentally, in spite of uncertainties regarding tissue geometry, vascular O2 distributions, and in vivo parameters such as oxygen consumption rate [21, 22, 40]. It could also be used for more complex models studying the interplay between tissue oxygenation and cell growth and the macroscopic dynamics of tumours [86–88].
Far from being exhaustive, this review illustrates the importance of simulating the many complex facets tumour oxygenation. Considerable amount of work remains to be performed, both in updating the models of oxygen transport into tissue according to new findings and validations regarding the involved mechanisms and parameters and in terms of new opportunities for simulating treatments. Therefore, the increased interest in light of ion therapy  with its included potential of reducing the hypoxic radioresistance opens up new directions for research in terms of studying the impact of hypoxia on radiation treatment. New models for taking into account the response of hypoxic cells to high LET radiation are being developed  and will be tested with the ultimate aim to counteract hypoxic radioresistance and achieve maximum clinical benefit for patients.
Conflict of Interests
The authors declare that they have no conflict of interests.
- P. Vaupel and A. Mayer, “Hypoxia in cancer: significance and impact on clinical outcome,” Cancer and Metastasis Reviews, vol. 26, no. 2, pp. 225–239, 2007.
- A. Daşu and I. Toma-Daşu, “Vascular oxygen content and the tissue oxygenation—a theoretical analysis,” Medical Physics, vol. 35, no. 2, pp. 539–545, 2008.
- M. Hockel, C. Knoop, K. Schlenger et al., “Intratumoral pO2 predicts survival in advanced cancer of the uterine cervix,” Radiotherapy and Oncology, vol. 26, no. 1, pp. 45–50, 1993.
- D. M. Brizel, R. K. Dodge, R. W. Clough, and M. W. Dewhirst, “Oxygenation of head and neck cancer: changes during radiotherapy and impact on treatment outcome,” Radiotherapy and Oncology, vol. 53, no. 2, pp. 113–117, 1999.
- J. Bussink, J. H. A. M. Kaanders, A. M. Strik, and A. J. Van Der Kogel, “Effects of nicotinamide and carbogen on oxygenation in human tumor xenografts measured with luminescense based fiber-optic probes,” Radiotherapy and Oncology, vol. 57, no. 1, pp. 21–30, 2000.
- N. A. Mayr, W. T. Yuh, J. C. Arnholt, et al., “Pixel analysis of MR perfusion imaging in predicting radiation therapy outcome in cervical cancer,” Journal of Magnetic Resonance Imaging, vol. 12, no. 6, pp. 1027–1033, 2000.
- J. H. A. M. Kaanders, K. I. E. M. Wijffels, H. A. M. Marres et al., “Pimonidazole binding and tumor vascularity predict for treatment outcome in head and neck cancer,” Cancer Research, vol. 62, no. 23, pp. 7066–7074, 2002.
- M. Kikuchi, T. Yamane, S. Shinohara, et al., “18F-fluoromisonidazole positron emission tomography before treatment is a predictor of radiotherapy outcome and survival prognosis in patients with head and neck squamous cell carcinoma,” Annals of Nuclear Medicine, vol. 25, no. 9, pp. 625–633, 2011.
- A. Sun, S. Johansson, I. Turesson, A. Daşu, and J. Sörensen, “Imaging tumor perfusion and oxidative metabolism in patients with head-and-neck cancer using 1- [11C]-acetate PET during radiotherapy: preliminary results,” International Journal of Radiation Oncology, Biology, Physics, vol. 82, no. 2, pp. 554–560, 2012.
- P. Vaupel, F. Kallinowski, and P. Okunieff, “Blood flow, oxygen and nutrient supply, and metabolic microenvironment of human tumors: a review,” Cancer Research, vol. 49, no. 23, pp. 6449–6465, 1989.
- J. Denekamp, A. Daşu, and A. Waites, “Vasculature and microenvironmental gradients: the missing links in novel approaches to cancer therapy?” Advances in Enzyme Regulation, vol. 38, no. 1, pp. 281–299, 1998.
- D. Hanahan and R. A. Weinberg, “The hallmarks of cancer,” Cell, vol. 100, no. 1, pp. 57–70, 2000.
- J. Folkman, “The vascularization of tumors,” Scientific American, vol. 234, no. 5, pp. 58–73, 1976.
- P. Vaupel, “Tumor microenvironmental physiology and its implications for radiation oncology,” Seminars in Radiation Oncology, vol. 14, no. 3, pp. 198–206, 2004.
- R. H. Thomlinson and L. H. Gray, “The histological structure of some human lung cancers and the possible implications for radiotherapy,” British Journal of Cancer, vol. 9, no. 4, pp. 539–549, 1955.
- J. M. Brown, “Evidence for acutely hypoxic cells in mouse tumours, and a possible mechanism of reoxygenation,” British Journal of Radiology, vol. 52, no. 620, pp. 650–656, 1979.
- D. J. Chaplin, R. E. Durand, and P. L. Olive, “Acute hypoxia in tumors: implications for modifiers of radiation effects,” International Journal of Radiation Oncology Biology Physics, vol. 12, no. 8, pp. 1279–1282, 1986.
- C. Bayer, K. Shi, S. T. Astner, C. A. Maftei, and P. Vaupel, “Acute versus chronic hypoxia: why a simplified classification is simply not enough,” International Journal of Radiation Oncology Biology Physics, vol. 80, no. 4, pp. 965–968, 2011.
- G. S. Adair, “The hemoglobin system. VI. The oxygen dissociation curve of hemoglobin,” Journal of Biological Chemistry, vol. 63, pp. 529–545, 1925.
- A. V. Hill, “The possible effects of the aggregation of the molecules of haemoglobin on its disassociation curve,” Journal of Physiology, vol. 41, pp. 4–7, 1910.
- A. S. Popel, “Theory of oxygen transport to tissue,” Critical Reviews in Biomedical Engineering, vol. 17, no. 3, pp. 257–321, 1989.
- D. Goldman, “Theoretical models of microvascular oxygen transport to tissue,” Microcirculation, vol. 15, no. 8, pp. 795–811, 2008.
- A. Krogh, “The number and distribution of capillaries in muscles with calculations of the oxygen pressure head necessary for supplying the tissue,” Journal of Physiology, vol. 52, no. 6, pp. 409–415, 1919.
- A. V. Hill, “The diffusion of oxygen and lactic acid through tissues,” Proceedings of the Royal Society B, vol. 104, pp. 39–96, 1928.
- I. F. Tannock, “Oxygen diffusion and the distribution of cellular radiosensitivity in tumours,” British Journal of Radiology, vol. 45, no. 535, pp. 515–524, 1972.
- T. W. Secomb and R. Hsu, “Analysis of oxygen delivery to tissue by microvascular networks,” Advances in Experimental Medicine and Biology, vol. 222, pp. 95–103, 1988.
- T. W. Secomb, R. Hsu, M. W. Dewhirst, B. Klitzman, and J. F. Gross, “Analysis of oxygen transport to tumor tissue by microvascular networks,” International Journal of Radiation Oncology Biology Physics, vol. 25, no. 3, pp. 481–489, 1993.
- B. W. Pogue, J. A. O'Hara, C. M. Wilmot, K. D. Paulsen, and H. M. Swartz, “Estimation of oxygen distribution in RIF-1 tumors by diffusion model-based interpretation of pimonidazole hypoxia and eppendorf measurements,” Radiation Research, vol. 155, no. 1, pp. 15–25, 2001.
- A. Daşu, I. Toma-Daşu, and M. Karlsson, “Theoretical simulation of tumour oxygenation and results from acute and chronic hypoxia,” Physics in Medicine and Biology, vol. 48, no. 17, pp. 2829–2842, 2003.
- C. J. Kelly and M. Brady, “A model to simulate tumour oxygenation and dynamic [18F]-Fmiso PET data,” Physics in Medicine and Biology, vol. 51, no. 22, pp. 5859–5873, 2006.
- A. Daşu and I. Toma-Daşu, “The relationship between vascular oxygen distribution and tissue oxygenation,” Advances in Experimental Medicine and Biology, vol. 645, pp. 255–260, 2009.
- P. S. Tang, “On the rate of oxygen consumption by tissues and lower organisms as a function of oxygen tension,” Quarterly Review of Biology, vol. 8, pp. 260–274, 1933.
- P. K. Nair, D. G. Buerk, and W. J. Whalen, “Cat carotid body oxygen metabolism and chemoreception described by a two-cytochrome model,” American Journal of Physiology, vol. 250, no. 2, pp. H202–H207, 1986.
- R. Hsu and T. W. Secomb, “A green's function method for analysis of oxygen delivery to tissue by microvascular networks,” Mathematical Biosciences, vol. 96, no. 1, pp. 61–78, 1989.
- T. W. Secomb, R. Hsu, E. T. Ong, J. F. Gross, and M. W. Dewhirst, “Analysis of the effects of oxygen supply and demand on hypoxic fraction in tumors,” Acta Oncologica, vol. 34, no. 3, pp. 313–316, 1995.
- A. P. Bassom, A. Ilchmann, and H. Voss, “Oxygen diffusion in tissue preparations with Michaelis-Menten kinetics,” Journal of Theoretical Biology, vol. 185, no. 1, pp. 119–127, 1997.
- T. W. Secomb, R. Hsu, E. Y. H. Park, and M. W. Dewhirst, “Green's function methods for analysis of oxygen delivery to tissue by microvascular networks,” Annals of Biomedical Engineering, vol. 32, no. 11, pp. 1519–1529, 2004.
- S. F. Petit, A. L. A. J. Dekker, R. Seigneuric et al., “Intra-voxel heterogeneity influences the dose prescription for dose-painting with radiotherapy: a modelling study,” Physics in Medicine and Biology, vol. 54, no. 7, pp. 2179–2196, 2009.
- D. Mönnich, E. G. C. Troost, J. H. A. M. Kaanders, W. J. G. Oyen, M. Alber, and D. Thorwarth, “Modelling and simulation of [18F]fluoromisonidazole dynamics based on histology-derived microvessel maps,” Physics in Medicine and Biology, vol. 56, no. 7, pp. 2045–2057, 2011.
- A. C. Skeldon, G. Chaffey, D. J. Lloyd, V. Mohan, D. A. Bradley, and A. Nisbet, “Modelling and detecting tumour oxygenation levels,” PLoS ONE, vol. 7, no. 6, Article ID e38597, 2012.
- P. M. Gullino, F. H. Grantham, and A. H. Courtney, “Utilization of oxygen by transplanted tumors in vivo,” Cancer Research, vol. 27, no. 6, pp. 1020–1030, 1967.
- P. L. Olive, C. Vikse, and M. J. Trotter, “Measurement of oxygen diffusion distance in tumor cubes using a fluorescent hypoxia probe,” International Journal of Radiation Oncology Biology Physics, vol. 22, no. 3, pp. 397–402, 1992.
- M. W. Dewhirst, T. W. Secomb, E. T. Ong, R. Hsu, and J. F. Gross, “Determination of local oxygen consumption rates in tumors,” Cancer Research, vol. 54, no. 13, pp. 3333–3336, 1994.
- F. Kallinowski, G. Tyler, W. Mueller-Klieser, and P. Vaupel, “Growth-related changes of oxygen consumption rates of tumor cells grown in vitro and in vivo,” Journal of Cellular Physiology, vol. 138, no. 1, pp. 183–191, 1989.
- D. F. Wilson, S. A. Vinogradov, B. W. Dugan, D. Biruski, L. Waldron, and S. A. Evans, “Measurement of tumor oxygenation using new frequency domain phosphorometers,” Comparative Biochemistry and Physiology, vol. 132, no. 1, pp. 153–159, 2002.
- L. S. Ziemer, W. M. F. Lee, S. A. Vinogradov, C. Sehgal, and D. F. Wilson, “Oxygen distribution in murine tumors: characterization using oxygen-dependent quenching of phosphorescence,” Journal of Applied Physiology, vol. 98, no. 4, pp. 1503–1510, 2005.
- J. P. Kirkpatrick, D. M. Brizel, and M. W. Dewhirst, “A mathematical model of tumor oxygen and glucose mass transport and metabolism with complex reaction kinetics,” Radiation Research, vol. 159, no. 3, pp. 336–344, 2003.
- B. M. Fenton, E. M. Lord, and S. F. Paoni, “Effects of radiation on tumor intravascular oxygenation, vascular configuration, development of hypoxia, and clonogenic survival,” Radiation Research, vol. 155, no. 2, pp. 360–368, 2001.
- J. Denekamp and A. Dau, “Inducible repair and the two forms of tumour hypoxia—time for a paradigm shift,” Acta Oncologica, vol. 38, no. 7, pp. 903–918, 1999.
- A. Daşu and J. Denekamp, “The impact of tissue microenvironment on treatment simulation,” Advances in Experimental Medicine and Biology, vol. 510, pp. 63–67, 2003.
- E. J. Hall, J. S. Bedford, and R. Oliver, “Extreme hypoxia; its effect on the survival of mammalian cells irradiated at high and low dose-rates,” British Journal of Radiology, vol. 39, no. 460, pp. 302–307, 1966.
- I. J. Spiro, K. A. Kennedy, R. Stickler, and C. C. Ling, “Cellular and molecular repair of X-ray-induced damage: dependence on oxygen tension and nutritional status,” Radiation Research, vol. 101, no. 1, pp. 144–155, 1985.
- C. C. Ling, E. Robinson, and D. C. Shrieve, “Repair of radiation induced damage—dependence on oxygen and energy status,” International Journal of Radiation Oncology Biology Physics, vol. 15, no. 5, pp. 1179–1186, 1988.
- F. Zölzer and C. Streffer, “Increased radiosensitivity with chronic hypoxia in four human tumor cell lines,” International Journal of Radiation Oncology Biology Physics, vol. 54, no. 3, pp. 910–920, 2002.
- N. Chan, M. Koritzinsky, H. Zhao et al., “Chronic hypoxia decreases synthesis of homologous recombination proteins to offset chemoresistance and radioresistance,” Cancer Research, vol. 68, no. 2, pp. 605–614, 2008.
- A. Daşu and I. Toma-Daşu, “Treatment modelling: the influence of micro-environmental conditions,” Acta Oncologica, vol. 47, no. 5, pp. 896–905, 2008.
- K. I. E. M. Wÿffels, J. H. A. M. Kaanders, P. F. J. W. Rijken et al., “Vascular architecture and hypoxic profiles in human head and neck squamous cell carcinomas,” British Journal of Cancer, vol. 83, no. 5, pp. 674–683, 2000.
- G. W. Barendsen, “Dose fraction, dose rate and iso-effect relationships for normal tissue responses,” International Journal of Radiation Oncology Biology Physics, vol. 8, no. 11, pp. 1981–1997, 1982.
- J. F. Fowler, “The linear-quadratic formula and progress in fractionated radiotherapy,” British Journal of Radiology, vol. 62, no. 740, pp. 679–694, 1989.
- T. Alper and P. Howard-Flanders, “Role of oxygen in modifying the radiosensitivity of E. coli B,” Nature, vol. 178, no. 4540, pp. 978–979, 1956.
- T. Alper, Cellular Radiobiology, Cambridge University Press, Cambridge, UK, 1979.
- E. J. Hall and A. J. Giaccia, Radiobiology For the Radiologist, Lippincott Williams & Wilkins, Philadelphia, Pa, USA, 6th edition, 2006.
- B. G. Wouters and J. M. Brown, “Cells at intermediate oxygen levels can be more important than the 'hypoxic fraction' in determining tumor response to fractionated radiotherapy,” Radiation Research, vol. 147, no. 5, pp. 541–550, 1997.
- F. M. Buffa, C. West, K. Byrne, J. V. Moore, and A. E. Nahum, “Radiation response and cure rate of human colon adenocarcinoma spheroids of different size: the significance of hypoxia on tumor control modelling,” International Journal of Radiation Oncology Biology Physics, vol. 49, no. 4, pp. 1109–1118, 2001.
- I. Toma-Daşu, A. Daşu, and M. Karlsoon, “The relationship between temporal variation of hypoxia, polarographic measurements and predictions of tumour response to radiation,” Physics in Medicine and Biology, vol. 49, no. 19, pp. 4463–4475, 2004.
- I. Toma-Dasu, A. Dasu, and A. Brahme, “Influence of acute hypoxia on radiation therapy outcome,” in Physical, Chemical and Biological Targeting in Radiatio Oncology, M. Mehta, B. R. Paliwal, and S. Bentzen, Eds., vol. 14 of AAPM Symposium Proceedings, pp. 111–117, Medical Physics Publishing, Madison, Wis, USA, 2005.
- I. Toma-Daşu, A. Daşu, and A. Brahme, “Dose prescription and optimisation based on tumour hypoxia,” Acta Oncologica, vol. 48, no. 8, pp. 1181–1192, 2009.
- E. Lindblom, The impact of hypoxia on tumour control probability in the high-dose range used in stereotactic body radiation therapy [M.S. thesis], Stockholm University, 2012.
- I. Toma-Daşu, A. Waites, A. Daşu, and J. Denekamp, “Theoretical simulation of oxygen tension measurement in tissues using a microelectrode: I. The response function of the electrode,” Physiological Measurement, vol. 22, no. 4, pp. 713–725, 2001.
- I. Toma-Daşu, A. Daşu, A. Waites, J. Denekamp, and J. F. Fowler, “Theoretical simulation of oxygen tension measurement in the tissue using a microelectrode: II. Simulated measurements in tissues,” Radiotherapy and Oncology, vol. 64, no. 1, pp. 109–118, 2002.
- I. Toma-Daşu, A. Daşu, A. Waites, and J. Denekamp, “Computer simulation of oxygen microelectrode measurements in tissues,” Advances in Experimental Medicine and Biology, vol. 510, pp. 157–161, 2003.
- G. Powathil, M. Kohandel, M. Milosevic, and S. Sivaloganathan, “Modeling the spatial distribution of chronic tumor hypoxia: implications for experimental and clinical studies,” Computational and Mathematical Methods in Medicine, vol. 2012, Article ID 410602, 11 pages, 2012.
- I. Toma-Daşu, A. Daşu, and A. Brahme, “Quantifying tumour hypoxia by PET imaging—a theoretical analysis,” Advances in Experimental Medicine and Biology, vol. 645, pp. 267–272, 2009.
- I. Toma-Dasu, J. Uhrdin, A. Daşu, and A. Brahme, “Therapy optimization based on non-linear uptake of PET tracers versus “linear dose painting”,” in World Congress on Medical Physics and Biomedical Engineering: Radiation Oncology, vol. 25/1 of IFMBE Proceedings, pp. 221–224, September 2009.
- I. Toma-Dasu and A. Dasu, “Biologically-optimised IMRT based on molecular imaging of tumour hypoxia—the impact of the tracer used,” in World Congress on Medical Physics and Biomedical Engineering, vol. 39 of IFMBE Proceedings, pp. 1742–1745.
- M. Busk, M. R. Horsman, and J. Overgaard, “Resolution in PET hypoxia imaging: voxel size matters,” Acta Oncologica, vol. 47, no. 7, pp. 1201–1210, 2008.
- C. C. Ling, J. Humm, S. Larson et al., “Towards multidimensional radiotherapy (MD-CRT): biological imaging and biological conformality,” International Journal of Radiation Oncology Biology Physics, vol. 47, no. 3, pp. 551–560, 2000.
- A. Brahme, “Biologically optimized 3-dimensional in vivo predictive assay-based radiation therapy using positron emission tomography-computerized tomography imaging,” Acta Oncologica, vol. 42, no. 2, pp. 123–136, 2003.
- D. J. Chaplin, P. L. Olive, and R. E. Durand, “Intermittent blood flow in a murine tumor: radiobiological effects,” Cancer Research, vol. 47, no. 2, pp. 597–601, 1987.
- J. Bussink, J. H. A. M. Kaanders, P. F. J. W. Rijken et al., “Vascular architecture and microenvironmental parameters in human squamous cell carcinoma xenografts: effects of carbogen and nicotinamide,” Radiotherapy and Oncology, vol. 50, no. 2, pp. 173–184, 1999.
- P. F. J. W. Rijken, H. J. J. A. Bernsen, J. P. W. Peters, R. J. Hodgkiss, J. A. Raleigh, and A. J. Van Der Kogel, “Spatial relationship between hypoxia and the (perfused) vascular network in a human glioma xenograft: a quantitative multi-parameter analysis,” International Journal of Radiation Oncology Biology Physics, vol. 48, no. 2, pp. 571–582, 2000.
- W. J. Koh, K. S. Bergman, J. S. Rasey et al., “Evaluation of oxygenation status during fractionated radiotherapy in human nonsmall cell lung cancers using [F-18]fluoromisonidazole positron emission tomography,” International Journal of Radiation Oncology Biology Physics, vol. 33, no. 2, pp. 391–398, 1995.
- S. M. Eschmann, F. Paulsen, C. Bedeshem et al., “Hypoxia-imaging with 18F-Misonidazole and PET: changes of kinetics during radiotherapy of head-and-neck cancer,” Radiotherapy and Oncology, vol. 83, no. 3, pp. 406–410, 2007.
- S. Roels, P. Slagmolen, J. Nuyts et al., “Biological image-guided radiotherapy in rectal cancer: is there a role for FMISO or FLT, next to FDG?” Acta Oncologica, vol. 47, no. 7, pp. 1237–1248, 2008.
- A. Dasu and I. Toma-Dasu, “Dose painting by numbers—do the practical limitations of the technique decrease or increase the probability of controlling tumours?” in World Congress on Medical Physics and Biomedical Engineering, vol. 39 of IFMBE Proceedings, pp. 1731–1734, 2012.
- M. Kohandel, M. Kardar, M. Milosevic, and S. Sivaloganathan, “Dynamics of tumor growth and combination of anti-angiogenic and cytotoxic therapies,” Physics in Medicine and Biology, vol. 52, no. 13, pp. 3665–3677, 2007.
- C. Harting, P. Peschke, K. Borkenstein, and C. P. Karger, “Single-cell-based computer simulation of the oxygen-dependent tumour response to irradiation,” Physics in Medicine and Biology, vol. 52, no. 16, pp. 4775–4789, 2007.
- W. Tuckwell, E. Bezak, E. Yeoh, and L. Marcu, “Efficient Monte Carlo modelling of individual tumour cell propagation for hypoxic head and neck cancer,” Physics in Medicine and Biology, vol. 53, no. 17, pp. 4489–4507, 2008.
- A. Brahme, “Recent advances in light ion radiation therapy,” International Journal of Radiation Oncology Biology Physics, vol. 58, no. 2, pp. 603–616, 2004.
- L. Antonovic, A. Brahme, Y. Furusawa, and I. Toma-Dasu, “Radiobiological description of the LET dependence of the cell survival of oxic and anoxic cells irradiated by carbon ions,” Journal of Radiation Research. In press.