Skip to main content

Theoretical analysis of the dose dependence of the oxygen enhancement ratio and its relevance for clinical applications

Abstract

Background

The increased resistance of hypoxic cells to ionizing radiation is usually believed to be the primary reason for treatment failure in tumors with oxygen-deficient areas. This oxygen effect can be expressed quantitatively by the oxygen enhancement ratio (OER). Here we investigate theoretically the dependence of the OER on the applied local dose for different types of ionizing irradiation and discuss its importance for clinical applications in radiotherapy for two scenarios: small dose variations during hypoxia-based dose painting and larger dose changes introduced by altered fractionation schemes.

Methods

Using the widespread Alper-Howard-Flanders and standard linear-quadratic (LQ) models, OER calculations are performed for T1 human kidney and V79 Chinese hamster cells for various dose levels and various hypoxic oxygen partial pressures (pO2) between 0.01 and 20 mmHg as present in clinical situations in vivo. Our work comprises the analysis for both low linear energy transfer (LET) treatment with photons or protons and high-LET treatment with heavy ions. A detailed analysis of experimental data from the literature with respect to the dose dependence of the oxygen effect is performed, revealing controversial opinions whether the OER increases, decreases or stays constant with dose.

Results

The behavior of the OER with dose per fraction depends primarily on the ratios of the LQ parameters alpha and beta under hypoxic and aerobic conditions, which themselves depend on LET, pO2 and the cell or tissue type. According to our calculations, the OER variations with dose in vivo for low-LET treatments are moderate, with changes in the OER up to 11% for dose painting (1 or 3 Gy per fraction compared to 2 Gy) and up to 22% in hyper-/hypofractionation (0.5 or 20 Gy per fraction compared to 2 Gy) for oxygen tensions between 0.2 and 20 mmHg typically measured clinically in hypoxic tumors. For extremely hypoxic cells (0.01 mmHg), the dose dependence of the OER becomes more pronounced (up to 36%). For high LET, OER variations up to 4% for the whole range of oxygen tensions between 0.01 and 20 mmHg were found, which were much smaller than for low LET.

Conclusions

The formalism presented in this paper can be used for various tissue and radiation types to estimate OER variations with dose and help to decide in clinical practice whether some dose changes in dose painting or in fractionation can bring more benefit in terms of the OER in the treatment of a specific hypoxic tumor.

Background

The poor treatment prognosis for tumors with oxygen-deficient areas is usually attributed to the decreased radiosensitivity of hypoxic cells. Hypoxia has been identified in many tumor types and the disadvantageous impact of hypoxia on local tumor control has been well recognized [1–4]. Due to the rapid development of noninvasive imaging methods to estimate the spatial distribution of the oxygen partial pressure within the tumor using different hypoxia markers [5–7], patient-specific treatment planning using dose painting or dose escalation in multifraction regimes based on functional hypoxia imaging is now considered a promising approach to overcome hypoxia in radiotherapy [8, 9].

The radioprotective effect of hypoxia can be expressed quantitatively by the oxygen enhancement ratio (OER), which is defined as the ratio of doses D(p 1 ) and D(p 2 ) given under two different oxygenation conditions with oxygen partial pressures p 1 ≤p 2 to produce the same biological effect. Depending on the interest of the researcher, the comparison can be done between hypoxic-aerobic, anoxic-hypoxic, aerobic-hypobaric etc. conditions. Historically the OER has been defined as the radiation dose under anoxic conditions (p 1 = 0 mmHg = const) divided by the dose under conditions of some partial pressure of oxygen p 2 = p. In this way, the OER typically increases from unity at anoxic conditions to approximately 3 for normoxic conditions in vitro. Alternatively, OER is often stated as the ratio of the dose to hypoxic cells (at different levels of p 1 = p) to the dose to aerobic cells (usually in air, p 2 = 160 mmHg = const). In this case the OER decreases with increasing oxygen partial pressure in the cell environment, and the term hypoxia reduction factor (HRF) is sometimes used instead of OER here. In this paper, we employ the most general definition of OER with variable p 1 and p 2 at the same time.

It has been observed that the OER depends on many factors such as the oxygen partial pressures (pO2) for hypoxic and aerobic conditions, the tissue type, the linear energy transfer (LET) of the radiation and the chosen cell survival level (or alternatively the applied local dose). In this work we focus on the dose dependence of the OER and investigate if and how this dependence is important for clinical applications in radiotherapy. To open the discussion we review the experimental data of the dose dependence of the OER published from 1975 up to 2010. In contrary to the relative biological effectiveness (RBE), which was observed to decrease with increasing dose per fraction independent of radiation type and cell line, no clear tendency for the dose dependence of the OER (decreasing, increasing or constant) was reported in the literature.

The purpose of this paper is to estimate the amount of potential OER variations with dose in low- and high-LET radiotherapy and to assess the potential impact for clinical applications. Based on the previously developed OER model [10] and parameters for several tissue types, we will evaluate two clinically relevant scenarios: the first scenario deals with dose painting strategies, where an inhomogeneous dose within the target volume according to the oxygenation status is prescribed. Typically this involves relatively small local changes of the dose in the order of ± 50% of the mean dose per fraction. The second scenario considers much larger changes in the dose per fraction which occur if the overall fractionation scheme is modified (hyper- or hypofractionation).

Methods

The survival of cells after exposure to a radiation dose is often described by the linear-quadratic (LQ) model [11]. This model is now in widespread use in both experimental and clinical radiobiology and generally works well in reproducing experimental results both in vitro and in vivo. We employed the standard LQ model with its two radiosensitivity parameters α and β in its most simplified form (without reoxygenation between fractions). Dose rate effects [12, 13] were also not taken into account.

In this simplified case the fraction S of cells that survive an applied dose D may be written as:

S = exp ( - α D - β D 2 )
(1)

The applied dose has to be changed to achieve the same biological effect under different irradiation conditions (e.g. high vs. low LET, aerobic vs. hypoxic environment etc.). This can be expressed by an enhancement factor (EF) defined as the ratio of doses given to a biological system under two different conditions, where the dose applied under the first condition 1 must be modified compared to the dose applied under the second condition 2 to obtain the same cell survival level:

E F = D 1 D 2 S = c o n s t
(2)

Due to the nonlinear form of the initial shoulder of the survival curves, the EF can vary depending on the choice of the specific survival level used as endpoint. This means the EF depends on the applied local dose (Figure 1A). When comparing low- and high-LET radiation (e.g. x-rays and carbon ion beams), this enhancement factor is called RBE and Eq. 1 can be written as:

Figure 1
figure 1

Dose dependence of Survival and Enhancement Factor. Panel A: Schematic illustration of the dependence of the enhancement factor (EF) on survival level and on dose. Panel B: Representative examples of the dose dependence of the EF: RBE (dashed-dotted line) and OER (solid lines) as a function of dose for high-LET radiation with argon ions (OERAr) and low-LET radiation with x-rays (OERX) calculated using experimental data for V79 Chinese hamster cells in vitro. The OER curves were obtained using Eq. 6 and the experimental data for LQ parameters α a , α h , β a and β h from [26] (OERAr, OERX1), [27] (OERX2) and [28] (OERX3). The RBE was calculated according to [29] using parameters α a and β a for argon ions and x-rays from [26].

R B E = D x - r a y s D c a r b o n S = c o n s t
(3)

If the same radiation type is compared under different oxygenation conditions (p 1 and p 2 , with p 1 ≤p 2 ), this factor defines the OER:

O E R = D ( p 1 ) D ( p 2 ) S = c o n s t
(4)

It was observed in the vast majority of experiments using various cell lines/tissue types and various radiation types that the RBE decreases with increasing dose per fraction (e.g. [14, 15]), see Figure 1B for a typical example. The reason is that the x-ray survival curve usually has a relatively large initial shoulder, whereas the shoulder for high-LET radiation is usually smaller and the initial slope steeper (larger α value). For the OER, the situation is more complex. Depending on the cell line and the radiation type, some investigators found an increasing, some a decreasing and others a constant OER with increasing dose (see Figure 1B, and a list of relevant experimental studies in Tables 1 and 2), and the changes in the shape of the survival curves are quite variable. Table 2 adds experimental OER studies for high LET, which was here limited to the range 80-120 keV/μm.

Table 1 Summary of references to published mammalian cell survival data measured in experiments at low LET
Table 2 Summary of references to published mammalian cell survival data measured in experiments at high LET

By equating the LQ-predicted surviving fractions (Eq. 1) for cells irradiated under aerobic (a) and hypoxic (h) conditions

S ( D a ) = S ( D h ) exp ( - α a D a - β a D a 2 ) = exp ( - α h D h - β h D h 2 )
(5)

and using Eq. 4 one obtains a simple formula for the OER depending on dose D h given to hypoxic cells (which needs to be the dose per fraction in a multifraction regime) and on tissue specific parameters α a , α h , β a and β h :

O E R ( D h , α a , β a , α h , β h ) = 2 D h β a / ( α a 2 + 4 β a ( α h D h + β h D h 2 ) − α a )
(6)

In the limit of very small doses per fraction (D h →0), the OER is given by the ratio α a /α h , whereas for very large doses per fraction (D h →∞), OER is determined by the β parameters only:

O E R ( D h → 0 ) = α a α h O E R ( D h → ∞ ) = β a β h
(7)

Since the first derivative of OER with respect to D never vanishes for D > 0 (unless α a /α h = (β a /β h )1/2), the OER increases with dose per fraction if the ratio α a /α h is smaller than (β a /β h )1/2 for a specific cell line, and decreases with dose for α a /α h > (β a /β h )1/2. In the case of α a /α h = (β a /β h )1/2, the OER is independent of dose. In contrast to the RBE, where the ratio of the α values is typically much larger than the square root of the ratio of the β values (leading to a decreasing RBE with dose), these ratios show larger variability for the OER.

Obviously, the OER depends not only on dose and tissue type but also on LET and pO2. Based on Eq. 6 we assume that these dependences are determined by the LQ parameters α(LET, pO 2 ) and β(LET, pO 2 ). The following calculations were performed using an OER model that we presented previously [10]. Briefly, our model is based on the experimental data from the literature. For an analysis of these data in order to obtain the dependence of the LQ parameters α and β on pO2 we used (like many other investigators [6, 16, 17]) the concept of Relative Radiosensitivity (RR) according to the Alper and Howard-Flanders model [18]. The RR describes the response of a biological system to radiation dependent on oxygen tension p in the cell environment. RR is maximized at high oxygen concentrations and approaches unity for low oxygen levels:

R R ( p ) = m â‹… p + 1 â‹… K p + K
(8)

Here m is the maximum radiosensitivity and K is the oxygen concentration at which the RR equals half of its maximum. Based on this concept we assume an Alper-Howard-Flanders dependence of α(p) and β1/2(p) on pO2. These functions have a similar shape as the relative radiosensitivity, although the maxima and minima are different. Furthermore, we suggest in the clinically relevant LET region (where the RBE increases with LET) a simple linear dependence of α on LET and suppose β to be independent of LET (the dependence of β on LET was also discussed in [10]). Taking the pO2 and LET dependence together, one obtains:

α ( L , p ) = ( a 1 + a 2 ⋅ L ) ⋅ p + ( a 3 + a 4 ⋅ L ) ⋅ K p + K β ( L , p ) = β ( p ) = b 1 ⋅ p + b 2 ⋅ K p + K
(9)

where L is LET and a 1 , a 2 , a 3 , a 4 , b 1 , and b 2 are constant coefficients that were estimated by fitting the experimental data in vitro from the literature [10]. Finally, equation (6) can be written as:

O E R ( D h , L , p a , p h ) = 2 D h β ( L , p a ) α 2 ( L , p a ) + 4 β ( L , p a ) ( α ( L , p h ) D h + β ( L , p h ) D h 2 ) - α ( L , p a )
(10)

with oxygen partial pressures p a and p h under hypoxic and aerobic conditions (p h ≤p a ). As mentioned in Background these oxygen partial pressures are two independent values of the oxygen tension and the OER depends on the choice of both of them. Equation (10) can be used to describe the OER for different radiation types (low-LET and high-LET) and for various oxygen levels relevant for cell experiments in vitro and clinical situations in vivo. This OER model is a simple tool to quantify the oxygen effect in a practical way. The results of our model for the dependence of OER on LET and pO2 as discussed in our previous paper are in good agreement with preclinical and clinical studies [10].

In this paper we deal with the dependence of OER on dose per fraction for different irradiation types and degrees of hypoxia. The OER calculations were done for two cell types (V79 Chinese hamster cells and T1 human kidney cells) because only for these cell lines there was sufficient experimental data both in the low- and high-LET area (Tables 1 and 2). For low LET (0.2-2 keV/μm), the mean values of α and β found in the literature were taken (V79: α a = 0.135 Gy-1, β a = 0.032 Gy-2, α h = 0.06 Gy-1, β h = 0.003 Gy-2; T1: α a = 0.10 Gy-1, β a = 0.047 Gy-2, α h = 0.02 Gy-1, β h = 0.009 Gy-2). The values for V79 are similar to the data used by Carlson et al. [19] in an OER modeling study for prostate cancer and might therefore be relevant for clinical applications as well. The tissue parameters for high LET were estimated by an LET-dependent fitting of the aerobic and hypoxic experimental data in vitro for α and β in the full high-LET region between 10 and 260 keV/μm [10] to determine the parameters at 100 keV/μm (V79: α a = 0.75 Gy-1, β a = 0.061 Gy-2, α h = 0.41 Gy-1, β h = 0.014 Gy-2; T1: α a = 0.62 Gy-1, β a = 0.067 Gy-2, α h = 0.35 Gy-1, β h = 0.019 Gy-2).

Results

OER values as a function of dose for V79 and T1 cells are shown in Figure 2 for low LET = 1 keV/μm (the middle of the range between 0.2 and 2 keV/μm) and high LET = 100 keV/μm (the middle of the range between 80 and 120 keV/μm). At low LET, OER is given for doses up to 20 Gy per fraction, whereas for high LET the physical absorbed dose is typically lower (due to the RBE), and the relevant dose is set up to 10 Gy per fraction. Of course, the predictions of the LQ model for high doses should be taken with great care. This model is well validated, both experimentally and theoretically, up to doses per fraction of about 10 Gy, and may be applicable to higher doses as well [20].

Figure 2
figure 2

Dose dependence of OER for V79 and T1 cells. Dependence of OER on dose per fraction given to hypoxic cells for V79 (Panels A and B) and T1 cell lines (Panels C and D) at low-LET (Panels A and C) and high-LET (B and D). Solid lines show the model calculations for oxygen partial pressures p h in vivo between 0.01 and 20 mmHg (p a = 30 mmHg). Dashed lines correspond to cell experiments in vitro under extreme hypoxia (p h = 0.01 mmHg, p a = 160 mmHg).

For V79 cells, the OER increases with dose at low and high LET, whereas it decreases with dose for T1 cells at low LET and stays almost constant for high LET. The dashed lines in Figure 2 show the OER for a typical cell experiment in vitro under extreme hypoxia (p h = 0.01 mmHg, p a = 160 mmHg). The ratios of the employed LQ parameters in vitro under hypoxic and aerobic conditions, which determine how and how strong the OER changes with increasing dose per fraction, were α a /α h = 2.25, (β a /β h )1/2 = 3.27 (V79 cells) and α a /α h = 5.0, (β a /β h )1/2 = 2.29 (T1 cells) at low LET and α a /α h = 1.83, (β a /β h )1/2 = 2.09 (V79) and α a /α h = 1.77, (β a /β h )1/2 = 1.88 (T1) at high LET. They can be compared with the ratios listed in Tables 1 and 2. If for example the value α a /α h for some experiment from Table 1 or 2 is lower than the value given above for V79 and the value (β a /β h )1/2 higher as above for the same cell line and LET range, the curves of the dose dependence of OER will be steeper than in Figure 2 and the dose dependence is more pronounced.

To make the situation clinically more realistic, OER is also calculated for a set of hypoxic oxygen partial pressures p h between 0.2 and 20 mmHg, whereas p a = 30 mmHg is referred to as aerobic (well oxygenated cells in vivo) [15, 21]. Since typical threshold values for the tumor oxygenation status in clinical and preclinical practice to estimate the treatment outcome of patients with a hypoxic tumor are < 2.5 mmHg, < 5 mmHg, < 10 mmHg and < 20 mmHg (see also Discussion in [10]), we calculated the dose dependence of the OER for these pO2 levels (solid lines in Figure 2). Because of the properties of clinically used Eppendorf histographs it is not possible to determine pO2 values for extreme hypoxia [22]. The lowest median pO2 values measured experimentally with Eppendorf histographs are typically 0.2 mmHg [6, 23]. Figure 2 shows our calculation also for this oxygen level. Of course, cells at much lower oxygen concentrations can exist within a tumor and contribute to treatment failure. As a "worst case approximation", we also calculated the dependence of OER on dose for such extreme hypoxic cells in vivo (p h = 0.01 mmHg, p a = 30 mmHg). The OER values relevant for two scenarios of clinical dose variations are detailed in Table 3. For scenario 1 (dose painting), a mean dose per fraction of 2 Gy (low LET) or 1 Gy (high LET) is assumed as the baseline, and the OERs for ± 50% of this dose are given (which can be considered as extreme values for the dose variation within the target volume). For scenario 2 (hyper-/hypofractionation), the same baseline is used and OERs at one fourth (extreme hyperfractionation) and ten times the baseline (extreme hypofractionation) are evaluated. The values in Table 3 are calculated for a low, but still clinically measurable value of p h = 0.2 mmHg. The higher the oxygen partial pressure in a tumor, the less pronounced will be the variations of the OER with changing dose per fraction. For extremely hypoxic cells within a tumor, these variations could be more distinct (see Figure 2). For p h = 0.01 mmHg and p a = 30 mmHg, the changes in OER amount up to 6% (19%) for V79 und 26% (36%) for T1 cells if the dose per fraction is varied from 2 to 0.5 Gy (20 Gy) for low-LET radiation. For high LET, the changes in OER for extreme hypoxia are comparable with the changes at p h = 0.2 mmHg (Table 3).

Table 3 OER variations with dose per fraction

Discussion

In previous modelling studies regarding the OER [6, 16, 17, 24, 25], both dose-independent and dose-dependent implementations of the OER were used, mostly based on the Alper-Howard-Flanders formula of relative radiosensitivity (Eq. 8) or some modifications of this equation. With respect to experimental data published in the literature, decreasing, increasing and constant OERs with increasing dose for various cell lines and various radiation types were reported (Tables 1 and 2). Corresponding to our calculation in the framework of the LQ model, the OER depends on dose and its behavior is determined by the ratios of the LQ parameters α a /α h and (β a /β h )1/2 under aerobic and hypoxic conditions (Eq. 7). Since these ratios can vary considerably with tissue/cell type, LET and pO2 (at least experimentally), this could explain the controversial findings from the publications in Tables 1 and 2. Only if these ratios are equal, the dose dependence disappears. This can be implemented on the modelling side if the pO2 dependence of α and β is given by α(p)= α a /f(p) and β(p) = β a /f(p)2 (with the same function f(p) for both α and β) as used by Malinen et al. [6] or Carlson et al. [25]. One can argue whether or not this special case is actually realized in all cell lines or tissue types [25] (at least within experimental uncertainties), and whether a mechanistic interpretation of the LQ parameters and the underlying microscopic processes of radiation damage (see e.g. [12]) supports this situation. Given the caveats of mechanistic interpretations in radiation biology and the relatively large experimental error bars, a dose dependence of the OER can certainly not be excluded, which motivated us to study the potential clinical impact of a dose dependent OER.

The calculations performed in this work can be used to estimate - for a certain cell line or tissue type, irradiation type and oxygenation condition - whether, how and with which magnitude the OER varies with dose. Although the underlying OER model [10] used here was primarily based on experimental data from the literature in vitro, due to the implemented Alper-Howard-Flanders concept it can also provide reasonable predictions for preclinical and clinical situations in vivo (as discussed in detail in [10]). We therefore conclude that our analysis of the dependence of OER on dose per fraction has also some value for the assessment of realistic clinical situations, at least qualitatively.

Our investigation of the oxygen effect at clinically measured median oxygen tensions p h between 0.2 and 20 mmHg shows a relatively moderate dependence of the OER on dose per fraction (Table 3). The higher the oxygen partial pressure in a tumor, the less pronounced is the variation of the OER with dose (Figure 2). For hypoxia-based dose painting in the target volume (scenario 1, with local dose variations in the target between 1 and 3 Gy for low LET), the changes of the local OER in vivo were below ± 11% relative to the baseline of 2 Gy per fraction (Table 3). For high LET (dose variations between 0.5 and 1.5 Gy), these variations were much smaller (up to 1%). If the dose per fraction is varied over a larger range (by changing the fractionation scheme, scenario 2), OER variations in the order of ± 21% (0.5 Gy vs. 2 Gy per fraction) and ± 22% (20 Gy vs. 2 Gy per fraction) can be seen for low LET. Again, this was much smaller for high LET (up to 1% difference in OER for 0.25 Gy vs. 1 Gy and up to 4% for 10 Gy vs. 1 Gy per fraction). For a potential portion of extremely hypoxic cells within a tumor the dose dependence of the OER becomes more pronounced in low-LET treatment (see Results), but stays still moderate for high LET. One has to note that the values above depend strongly on the oxygenation conditions and tissue types, and the direction and magnitude of the OER variation could differ from the OERs for V79 and T1 cell lines. Nevertheless, the formalism presented in this paper can be used to estimate OER variations with dose and help to decide in clinical practice whether some changes in fractionation (hyper- or hypofractionation) can bring more benefit in the treatment of patients with a specific hypoxic tumor. For example, for human salivary gland (HSG) tumor cells the OER decreases with increasing dose under high-LET irradiation (Table 2). This means that a hypofractionated treatment in high-LET radiotherapy with heavy ions could be more advantageous with respect to the oxygen effect for patients with this kind of tumor. However, this can only be confirmed by clinical studies using high LET and various fractionation schemes, which (to the best of our knowledge) are currently not available for this tumor type. Alternatively, small animal experiments with implanted human tumors (which often exhibit large hypoxic fractions) could be a means to evaluate the impact of OER for various fractions schemes in a pre-clinical setting, and to validate the findings of our modeling study.

Another example can be based on the data published by Malinen et al. [6]. The authors investigated a spontaneous sarcoma of a dog with different hypoxic compounds within the tumor. In the most hypoxic area the measured average oxygen partial pressure p h was 0.2 mmHg. Using their data we calculated for this pO2 value α a /α h = 1.75 and (β a /β h )1/2 = 3.25 (Eq. 4 and Table 1 in [6]). Compared with our analysis for V79 cell line under low-LET irradiation (α a /α h = 2.25 and (β a /β h h )1/2 = 3.27 and Figure 2A) the dependence of OER on dose per fraction for sarcomas from the study by Malinen et al. is also relatively moderate. The higher the oxygen partial pressure in a tumor, the less significant is the variation in the OER with the choice of dose per fraction. However, if hypoxia based treatment planning shall be performed, a larger number of fractions with smaller doses (hyperfractionation) could - due to the increasing OER with dose - bring more benefit with respect to the oxygen effect for the treatment of the sarcoma presented by Malinen et al. [6].

Although the OER seems to be only moderately dependent on the choice of dose per fraction, tumor hypoxia itself has a large negative effect on cell killing and there is potential for large errors in calculation of alternative dose fractionation schemas using models that do not account for tumor hypoxia at all [19], even if the dose dependence is not considered explicitly. In the case where measurements of oxygen partial pressures in a tumor are possible (e.g. with an Eppendorf histograph or noninvasive methods), the additional dose to hypoxic areas of the tumor required to achieve a constant biological effect in the whole target can be calculated in the framework of our model (Eq. 10) for various radiation types, oxygen tensions and tissue types. In the long run, this may help to overcome the adverse effects of low oxygen concentrations in many tumors.

Conclusions

Since tumors with hypoxic areas exist and the treatment outcome of patients with hypoxic tumors is relatively poor, new predictive models are required to individualize and improve the treatment strategies for radiotherapy. In this work we investigated theoretically the importance of the dose dependence of the OER for clinical applications in radiotherapy. The analysis was performed for two scenarios: small dose variations within the target during dose painting and larger changes of the dose per fraction for different fractionation schemes. The calculations were performed for both low-LET treatment with photons or protons and high-LET treatment with heavy ions.

According to our analysis, the OER in clinical practice is moderately sensitive to the choice of dose per fraction. The OER can increase, decrease or remain a constant with increasing dose per fraction, and this behavior is determined by the ratios of the LQ parameters under hypoxic and aerobic conditions. These effects should be taken into account in hypoxia based treatment plan optimization. The formalism presented in this paper can be used to estimate OER variations with dose and to help with clinical decisions about any changes in dose prescription or treatment planning with respect to the oxygen effect. If simplified models without explicit consideration of the dose dependence are used for optimization in dose painting or changes of the fractionation scheme, our methods can be used to estimate the potential error in OER due to dose variations.

Abbreviations

EF:

Enhancement factor

HRF:

Hypoxia reduction factor

LET:

Linear energy transfer

LQ model:

Linear-quadratic model

OER:

Oxygen enhancement ratio

pO2:

Oxygen partial pressure

RBE:

Relative biological effectiveness

RR:

Relative Radiosensitivity.

References

  1. Höckel M, Schlenger K, Mitze M, Schäffer U, Vaupel P: Hypoxia and radiation response in human tumors. Semin Radiat Oncol. 1996, 6: 3-9. 10.1016/S1053-4296(96)80031-2.

    Article  PubMed  Google Scholar 

  2. Rofstad EK, Sundfør K, Lyng H, Tropé CG: Hypoxia-induced treatment failure in advanced squamous cell carcinoma of the uterine cervix is primarily due to hypoxia-induced radiation resistance rather than hypoxia-induced metastasis. Br J Cancer. 2000, 83: 354-359. 10.1054/bjoc.2000.1266.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  3. Vaupel P, Mayer A: Hypoxia in cancer: significance and impact on clinical outcome. Cancer Metastasis Rev. 2007, 26: 225-239. 10.1007/s10555-007-9055-1.

    Article  CAS  PubMed  Google Scholar 

  4. Wang JZ, Li XA, Mayr NA: Dose escalation to combat hypoxia in prostate cancer: a radiobiological study on clinical data. Br J Radiol. 2006, 79: 905-911. 10.1259/bjr/18700614.

    Article  CAS  PubMed  Google Scholar 

  5. Cooper RA, Carrington BM, Loncaster JA, Todd SM, Davidson SE, Logue JP, Luthra AD, Jones AP, Stratford I, Hunter RD, West CM: Tumour oxygenation levels correlate with dynamic contrast-enhanced magnetic resonance imaging parameters in carcinoma of the cervix. Radiaother Oncol. 2000, 57: 53-59. 10.1016/S0167-8140(00)00259-0.

    Article  CAS  Google Scholar 

  6. Malinen E, Søvik Å, Hristov D, Bruland ØS, Olsen DR: Adapting radiotherapy to hypoxic tumours. Phys Med Biol. 2006, 51: 4903-4921. 10.1088/0031-9155/51/19/012.

    Article  CAS  PubMed  Google Scholar 

  7. Mönnich D, Troost EG, Kaanders JH, Oyen WJ, Alber M, Thorwarth D: Modelling and simulation of [18F]fluoromisonidazole dynamics based on histology-derived microvessel maps. Phys Med Biol. 2011, 56: 2045-2057. 10.1088/0031-9155/56/7/009.

    Article  PubMed  Google Scholar 

  8. Søvik Å, Malinen E, Olsen DR: Strategies for biologic image-guided dose escalation: a review. Int J Radiat Oncol Biol Phys. 2009, 73: 650-658. 10.1016/j.ijrobp.2008.11.001.

    Article  PubMed  Google Scholar 

  9. Thorwarth D, Alber M: Implementation of hypoxia imaging into treatment planning and delivery. Radiother Oncol. 2010, 97: 172-175. 10.1016/j.radonc.2010.05.012.

    Article  PubMed  Google Scholar 

  10. Wenzl T, Wilkens JJ: Modelling of the oxygen enhancement ratio for ion beam radiation therapy. Phys Med Biol. 2011, 56: 3251-3268. 10.1088/0031-9155/56/11/006.

    Article  PubMed  Google Scholar 

  11. Kellerer AM, Rossi HH: A generalized formulation of dual radiation action. Radiat Res. 1978, 75: 471-488. 10.2307/3574835.

    Article  Google Scholar 

  12. Sachs RK, Hahnfeld P, Brenner DJ: The link between low-LET dose-response relations and the underlying kinetics of damage production/repair/misrepair. Int J Radiat Biol. 1997, 72: 351-374. 10.1080/095530097143149.

    Article  CAS  PubMed  Google Scholar 

  13. Carlson DJ, Stewart RD, Li XA, Jennings K, Wang JZ, Guerrero M: Comparison of in vitro and in vivo alpha/beta ratios for prostate cancer. Phys Med Biol. 2004, 49: 4477-4491. 10.1088/0031-9155/49/19/003.

    Article  PubMed  Google Scholar 

  14. Paganetti H: Relative biological effectiveness (RBE) values for proton beam therapy. Int J Radiat Oncol Biol Phys. 2002, 53: 407-421. 10.1016/S0360-3016(02)02754-2.

    Article  PubMed  Google Scholar 

  15. Hall EJ, Giaccia AJ: Radiobiology for the radiologist. 2006, Philadelphia: Lippincott Williams and Wilkins, 6

    Google Scholar 

  16. Wouters BG, Brown JM: Cells at intermediate oxygen levels can be more important than the "hypoxic fraction" in determining tumor response to fractionated radiotherapy. Radiat Res. 1997, 147: 541-550. 10.2307/3579620.

    Article  CAS  PubMed  Google Scholar 

  17. Petit SF, Dekker AL, Seigneuric R, Murrer L, van Riel NA, Nordsmark M, Overgaard J, Lambin P, Wouters BG: Intra-voxel heterogeneity influences the dose prescription for dose-painting with radiotherapy: a modelling study. Phys Med Biol. 2009, 54: 2179-2196. 10.1088/0031-9155/54/7/022.

    Article  PubMed  Google Scholar 

  18. Alper T, Howard-Flanders P: Role of oxygen in modifying the radiosensitivity of E. coli B. Nature. 1956, 178: 978-979. 10.1038/178978a0.

    Article  CAS  PubMed  Google Scholar 

  19. Carlson DJ, Keall PJ, Loo BW, Chen ZJ, Brown JM: Hypofractionation results in reduced tumor cell kill compared to conventional fractionation for tumors with regions of hypoxia. Int J Radiat Oncol Biol Phys. 2011, 79: 1188-1195. 10.1016/j.ijrobp.2010.10.007.

    Article  PubMed Central  PubMed  Google Scholar 

  20. Brenner DJ: The linear-quadratic model is an appropriate methodology for determining isoeffective doses at large doses per fraction. Semin Radiat Oncol. 2008, 18: 234-239. 10.1016/j.semradonc.2008.04.004.

    Article  PubMed Central  PubMed  Google Scholar 

  21. Höckel M, Vaupel P: Tumor hypoxia: definitions and current clinical, biologic, and molecular aspects. J Natl Cancer Inst. 2001, 93: 266-276. 10.1093/jnci/93.4.266.

    Article  PubMed  Google Scholar 

  22. Stone HB, Brown JM, Phillips TL, Sutherland RM: Oxygen in human tumors: correlations between methods of measurement and response to therapy. Radiat Res. 1997, 136: 422-434.

    Article  Google Scholar 

  23. Brizel DM, Sibley GS, Prosnitz LR, Scher RL, Dewhirst MW: Tumor hypoxia adversely affects the prognosis of carcinoma of the head and neck. Int J Radiat Oncol Biol Phys. 1997, 38: 285-289. 10.1016/S0360-3016(97)00101-6.

    Article  CAS  PubMed  Google Scholar 

  24. Antipas VP, Stamatakos GS, Uzunoglu NK, Dionysiou DD, Dale RG: A spatio-temporal simulation model of the response of solid tumours to radiotherapy in vivo: parametric validation concerning oxygen enhancement ratio and cell cycle duration. Phys Med Biol. 2004, 49: 1485-1504. 10.1088/0031-9155/49/8/008.

    Article  PubMed  Google Scholar 

  25. Carlson DJ, Stewart RD, Semenenko VA: Effects of oxygen on intrinsic radiation sensitivity: A test of the relationship between aerobic and hypoxic linear-quadratic (LQ) model parameters. Med Phys. 2006, 33: 3105-3115. 10.1118/1.2229427.

    Article  CAS  PubMed  Google Scholar 

  26. Raju MR, Amols HI, Bain E, Carpenter SG, Cox RA, Robertson JB: A heavy particle comparative study. Part III. OER and RBE. Br J Radiol. 1978, 51: 712-719. 10.1259/0007-1285-51-609-712.

    Article  CAS  PubMed  Google Scholar 

  27. Prise KM, Folkard M, Davieset S, Michael BD: The irradiation of V79 mammalian cells by protons with energies below 2 MeV. II. Measurement of oxygen enhancement ratios and DNA damage. Int J Radiat Biol. 1990, 58: 261-277. 10.1080/09553009014551611.

    Article  CAS  PubMed  Google Scholar 

  28. Skarsgard LD, Harrison I: Dose dependence of the oxygen enhancement ratio (OER) in radiation inactivation of Chinese hamster V79-171 cells. Radiat Res. 1991, 127: 243-247. 10.2307/3577937.

    Article  CAS  PubMed  Google Scholar 

  29. Wilkens JJ, Oelfke U: A phenomenological model for the relative biological effectiveness in therapeutic proton beams. Phys Med Biol. 2004, 49: 2811-2825. 10.1088/0031-9155/49/13/004.

    Article  CAS  PubMed  Google Scholar 

  30. Chapman JD, Gillespie CJ, Reuvers AP, Dugle DL: The inactivation of Chinese hamster cells by × rays: the effects of chemical modifiers on single- and double-events. Radiat Res. 1975, 64: 365-375. 10.2307/3574272.

    Article  CAS  PubMed  Google Scholar 

  31. Watts ME, Hodgkissn RJ, Jones R, Fowler JF: Radiosensitization of Chinese hamster cells by oxygen and misonidazole at low X-ray doses. Int J Radiat Biol. 1986, 50: 1009-1021. 10.1080/09553008614551421.

    Article  CAS  Google Scholar 

  32. Sapora O, Barone F, Belli M, Maggi A, Quintiliani M, Tabocchini MA: Relationship between cell killing, mutation induction and DNA damage in X-irradiated V79 cells; the influence of oxygen und DMSO. Int J Radiat Biol. 1991, 60: 467-482. 10.1080/09553009114552321.

    Article  CAS  PubMed  Google Scholar 

  33. Sprong D, Janssen HL, Vens C, Begg AC: Resistance of hypoxic cells to ionizing radiation is influenced by homologous recombination status. Int J Radiat Oncol Biol Phys. 2006, 64: 562-572. 10.1016/j.ijrobp.2005.09.031.

    Article  CAS  PubMed  Google Scholar 

  34. Raju MR, Amols HI, Bain E, Carpenter SG, Cox RA, Robertson JB: OER and RBE for negative pion beams of different peak widths. Br J Radiol. 1979, 52: 494-498. 10.1259/0007-1285-52-618-494.

    Article  CAS  PubMed  Google Scholar 

  35. Freyer JP, Jarrett K, Carpenter S, Raju MR: Oxygen enhancement ratio as a function of dose and cell cycle phase for radiation-resistant and sensitive CHO cells. Radiat Res. 1991, 127: 297-307. 10.2307/3577945.

    Article  CAS  PubMed  Google Scholar 

  36. Ma N, Durante M, Lai S, Tinganelli W, Weyrather WK: Cellular response of CHO-K1 cells to X-ray irradiation under different states of oxygenation. GSI Scientific Report 2010. Edited by: Große K. 2010, Darmstadt: GSI Helmholtzzentrum für Schwerionenforschung GmbH, 460-

    Google Scholar 

  37. Hall EJ: The effect of Hypoxia on the repair of sublethal radiation damage in cultured mammalian cells. Radiat Res. 1972, 49: 405-415. 10.2307/3573277.

    Article  CAS  PubMed  Google Scholar 

  38. Curtis SB: The OER of mixed high- and low-LET radiation. Radiat Res. 1976, 65: 566-572. 10.2307/3574388.

    Article  CAS  PubMed  Google Scholar 

  39. Tenforde TS, Curtis SB, Crabtree KE, Tenforde SD, Schilling WA, Howard J, Lyman JT: In vivo cell survival and volume response characteristics of rat rhabdomyosarcoma tumors irradiated in the extended peak region of carbon- and neon-ion beams. Radiat Res. 1980, 83: 42-56. 10.2307/3575257.

    Article  CAS  PubMed  Google Scholar 

  40. Urano M, Goitein M, Verhey L, Mendiondo O, Suit HD, Koehler A: Relative biological effectiveness of a high energy modulated proton beam using a spontaneous murine tumor in vivo. Int J Radiat Oncol Biol Phys. 1980, 6: 1187-1193. 10.1016/0360-3016(80)90172-8.

    Article  CAS  PubMed  Google Scholar 

  41. Blakely EA, Tobias CA, Yang TC, Smith KC, Lyman JT: Inactivation of human kidney cells by high-energy monoenergetic heavy-ion beams. Radiat Res. 1979, 80: 122-160. 10.2307/3575121.

    Article  CAS  PubMed  Google Scholar 

  42. Guichard M, Lachet B, Malaise EP: Measurement of RBE, OER and recovery of potentially lethal damage of a 645 MeV helium ion beam using EMT6 cells. Radiat Res. 1977, 71: 413-29. 10.2307/3574684.

    Article  CAS  PubMed  Google Scholar 

  43. Born R: Die Abhängigkeit des Sauerstoffeffekts von der Dosisleistung. Strahlentherapie. 1978, 154: 650-652.

    CAS  PubMed  Google Scholar 

  44. Shridhar R, Estabrook W, Yudelev M, Rakowski J, Burmeister J, Wilson GD, Joiner MC: Characteristic 8 keV × rays possess radiobiological properties of higher-LET radiation. Rad Res. 2010, 173: 290-297. 10.1667/RR1782.1.

    Article  CAS  Google Scholar 

  45. Furusawa Y, Fukutsu K, Aoki M, Itsukaichi H, Eguchi-Kasai K, Ohara H, Yatagai F, Kanai T, Ando K: Inactivation of aerobic and hypoxic cells from three different cell lines by accelerated 3He-, 12C- and 20Ne-ion beams. Radiat Res. 2000, 154: 485-496. 10.1667/0033-7587(2000)154[0485:IOAAHC]2.0.CO;2.

    Article  CAS  PubMed  Google Scholar 

  46. Curtis SB, Schilling WA, Tenforde TS, Crabtree KE, Tenforde SD, Howard J, Lyman JT: Survival of oxygenated and hypoxic tumor cells in the extended-peak regions of heavy charged-particle beams. Radiat Res. 1982, 90: 292-309. 10.2307/3575707.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

Supported by DFG Cluster of Excellence: Munich-Centre for Advanced Photonics and by DFG grant WI 3745/1-1.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Tatiana Wenzl.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

All authors have made substantive intellectual contributions to the published research. Both TW and JJW have been involved in analyzing and interpreting the data from the literature, developing the model and preparation of the manuscript. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Wenzl, T., Wilkens, J.J. Theoretical analysis of the dose dependence of the oxygen enhancement ratio and its relevance for clinical applications. Radiat Oncol 6, 171 (2011). https://doi.org/10.1186/1748-717X-6-171

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1748-717X-6-171

Keywords