Fri, 06 Dec 2024 in Brazilian Journal of Radiation Sciences
Effect of radiobiological parameters on the TCP for breast cancer radiotherapy
Abstract
Breast cancer remains the most prevalent malignancy affecting women globally. Among the various treatment modalities, radiotherapy stands out as a cornerstone for tumor eradication. This research explores the impact of radiosensitivity parameters on Tumor Control Probability (TCP) in breast cancer, with an emphasis on distinct radiotherapeutic techniques such as conventional, hypofractionated, and FAST, as well as the role of tumor repopulation. Based on the literature review, we obtained data on α and β radiosensitivity parameters, cell repopulation rates and standard breast cancer treatment protocols. These parameters informed the calculation of the fraction of cells surviving irradiation via the linear-quadratic model, facilitating an assessment of treatment efficacy through the Poissonian TCP model. Our findings underscore the critical influence of radiosensitivity parameters α and β on treatment outcomes, with β emerging as the predominant factor due to its quadratic contribution to the survival fraction. Moreover, our analysis indicates that tumor growth is negligible relative to the substantial cell mortality induced by radiation in the case of breast cancer. Techniques such as FAST and hypofractionated radiotherapy were identified as particularly effective, offering expedited tumor control, especially with elevated α and β values. The quadratic term β significantly enhances treatment success, while tumor repopulation exerts minimal influence on TCP, corroborating previous model comparisons. Notably, higher doses per fraction, rather than increased cumulative doses, were associated with improved TCP, providing a critical insight for optimizing radiotherapy protocols. Currently, radiobiology is not systematically integrated into clinical practice, and its analysis through PCT optimizes radiotherapy treatments, improving patient quality of life and healthcare delivery.
Main Text
1. INTRODUCTION
Breast cancer is the most common cancer in the world, regardless of the country's level of development [1]. Furthermore, excluding non-melanoma skin tumors, breast cancer ranks first in terms of cancer incidence and mortality in women, with 66,280 new cases estimated for 2022 in Brazil [2]. Early diagnosis of the disease is essential for a positive prognosis, and mammography is considered the most effective method of detection [3].
Cancer is a disease of uncontrolled proliferation by transformed cells subject to evolution by natural selection [4]. Treatments include surgery, radiotherapy, and systemic chemotherapy, whether or not these modalities are used together [5, 6]. There are approximately 66,000 new cases for the three-year period 2020-2023 [7].
The goal of radiotherapy is to use ionizing radiation to destroy tumor cells in a planned manner, using protocols that minimize damage to adjacent healthy tissues. Radiotherapy techniques can be divided into three: conventional, which has a lower daily dose, ranging from 1.8 to 2.0 Gy, with an average total dose of 45 to 60 Gy, and the technique with the highest number of fractions, divided into 25 to 30 [8]. Hypofractionated is the technique with doses above 2.0 Gy, resulting in a lower accumulated dose of 40 Gy, as it is fractionated into 15 and 16 fractions [9]. Finally, the FAST technique uses the highest daily dose of the three, 5.2 Gy, with a total dose of 26 Gy delivered in 5 uninterrupted fractions [10].
Mathematical modeling plays an important role in understanding and optimizing treatment, whether in comprehending tumor dynamics or the quality of diagnostic and therapeutic systems [11, 12], as it provides insight into tumor dynamics and the treatment itself without experimental costs [13, 14], making it very important for patient care. The linear quadratic model (LQM) is the most widely used representation for quantifying the fraction of cells that survive a given dose of radiation [14]. Although it was discovered empirically by experimental fitting, over time its parameters have gained a mechanistic interpretation related to the damage caused to DNA strands [15-17]. The LQM is expressed as a function of the alpha (α) and beta (β) radiosensitivity parameters of the tissue according to Equation 1:
F(D)=e-αD-βD2
where α refers to the one hit event on the DNA and is designated as a linear parameter (αD), and β refers to the multiple hit on the DNA and is designated as a quadratic parameter (βD²) [18].
Some mathematical probability models make it possible to estimate the treatment success of a given protocol, widely known in the literature as the Tumor Control Probability or TCP. The TCP makes it possible to compare radiotherapy protocols and techniques to optimize clinical decisions [20, 21]. However, to obtain more realistic responses, the rate of cell repopulation between irradiations must be taken into account, since cells are constantly proliferating. This study investigates the effects of radiosensitivity on the Poissonian Tumor Control Probability (TCP), in terms of α and β parameters, considering conventional, hypofractionated and FAST radiotherapy techniques, as well as tumor repopulation. Therefore, this study is concerned with the effect of radiobiological parameters on the response to TCP.
2. MATERIALS AND METHODS
2.1. Gathering data from the literature
To perform the calculations and visualize the results graphically, this work was based on the literature review conducted by van Leeuwen et al.[21], which in turn was based on the work of Qi et al.[22]. From which the values of the radiosensitivity parameters α and β derived from the linear-quadratic model were obtained, as shown in Table 1. All parameters were taken from the work of Leeuwen et al. [21], as they are based on a large review of randomized clinical trials for breast irradiation, while the literature is vague regarding more recent studies.
In addition, values were collected in the literature [24] for the clinical dose deliberation parameters used in practice for the conventional, hypofractionated, and FAST techniques in terms of total dose, fractions, and dose per fraction, as shown in Table 2. A escolha desse documento [24] se deu por se tratar de um documento guia para prática de oncologia clínica, baseada nas revisões da mais recente literatura.
2.2. TCP calculation
The probability calculations were implemented in Python on a conventional computer using the free Google Colab platform. From the survival fraction described in equation 1, one can derive the Poisson TCP model in equation 2,
TCP=e-N0F(D)
in which Gong and colleagues [26] showed that for low proliferation tumors, such as breast and prostate, more sophisticated TCP approaches lead to basically the same results as simple approaches, such as the Poisson approach. Therefore, this paper focuses on the use of the Poisson TCP.
TCPs were calculated by varying the radiosensitivity, growth rate, and dose deliberation parameters according to the protocols described in Tables 1 and 2, respectively. The values were substituted into Equation 1, which corresponds to the function F(D) in Equation 2, multiplied by the initial number of cells (N0 = 109) [20].
3. RESULTS AND DISCUSSION
The plots in Figure 1 A, B, C, D, E and F evaluate the effect of the α parameter on TCP for three different breast radiotherapy protocols: conventional, hypofractionated, and FAST. Thus, α = [0.04; 0.16] Gy-1, while the β parameter is kept constant and corresponds to 0.028 Gy-2. These values for α e β were used according to Leeuwen et al. [21] and Qi et al.[22]. The solid curves represent the TCPs without repopulation, while the dashed curves include cell repopulation.
The plots in Figure 2 A, B, C, D, E and F evaluate the effect of the β parameter on TCP for three different breast radiotherapy protocols: conventional, hypofractionated, and FAST. The β parameter varies in the range β=[0.01; 0.06] Gy-2, while α is kept constant and corresponds to 0.08 Gy-1. Again, the solid curves exclude tumor repopulation during treatment, while the dashed curves include tumor repopulation.
When tumor repopulation is taken into account, TCP values decrease in the absence of treatment, leaving a sloping profile as tumor growth becomes dominant in the face of radiation death. As an example, we have the blue and purple dotted curves in Figure 1D. The same behavior is repeated for all dashed curves in Figures 1 and 2.
For the curves without repopulation, represented by the solid lines, there is only a less accentuated horizontal plateau or staircase profile: the TCP remains constant because there is no death due to radiation or tumor growth. This behavior is more evident in the purple curves in Figures 1D and 1E, although it is repeated in the others, except for the cyan curve in Figure 1F, which quickly reaches TCP=1 due to the high values of α and β.
Repopulation has a direct impact on treatment time, as the tumor may grow again over the weekend or between radiation sessions during the week. Eradication is delayed due to repopulation. However, as shown in a previous study [25], tumor repopulation was observed to be subtle mainly in tumors with high α and β values, such as glioblastoma, head and neck cancer.
Figures 1A, 1B and 1C show the conventional, hypofractionated and FAST techniques, respectively. For the same protocol, the curves follow the same behavior in terms of time to TCP for different values of the parameter. One can observe that the pink curve with the highest α value (0.16 Gy-1) is always the first to reach treatment success (TCP=100%) compared to the cyan curve (= 0.04 Gy-1). The horizontal time difference for TCP=1 between these curves is Δt ~ 4 days. This parameter weights the effect of radiation in the low-dose regime such that the higher its value, the faster TCP=1 is reached. In this regime, according to the mechanistic interpretation, one DNA hits predominate.
In Figures 2A, 2B, and 2C, one can see that in each protocol, the effect of the β parameter on TCP is more drastic compared to the effect of the parameter (from Figure 1). This is most evident when looking at the horizontal distance between the cyan and pink curves in Figure 2A: when β varies from 0.01 to 0.06 Gy-2, the time to tumor control varies from Δt ~ 20 days. In its equivalent, Figure 1A, the horizontal difference is only Δt ~ 4 days, showing the drastic effect of the β parameter on TCP. This is typically related to the multiple-hit events on the DNA [31]. In Figures 2B and 2C, such difference is not evident for the dose accumulates very quickly due to the hypofractionation and FAST protocols.
In terms of protocols, comparing Figures 1A, 1B and 1C with their equivalents 2A, 2B and 2C, one can see that TCPs achieve success more quickly when the protocol provides a higher dose per fraction and more fractions. In particular, the effect of a higher dose per fraction outweighs the effect of a higher cumulative dose, contradicting the intuition that the total dose would dominate treatment success. Thus, the FAST technique guarantees faster tumor control, followed by the hypofractionated technique and finally the conventional technique. However, it should be noted that successful treatment must include the optimization of maximum doses in tumors and minimum doses in healthy tissue.
4. CONCLUSIONS
This analysis highlights the significant influence of the radiosensitivity parameters α and β in the calculation of Tumor Control Probability (TCP), where subtle variations result in significant changes in tumor control times, with Δt ~ 4 days and Δt ~ 20 days for variations in α and β, respectively, in conventional protocols.
In addition, the β parameter emerges as the most important factor in TCP formulation due to its quadratic weight in the survival fraction, because it decreases the function faster and thus increases tumor control faster.
The tumor repopulation rate has minimal influence on tumor control, as it is superimposed on the cell death rates induced by the radiation doses according to Equation 1. This observation is consistent with studies comparing TCP models in tumors of varying aggressiveness.
In protocol comparisons, higher doses per fraction accelerate treatment, contrary to the intuition that a higher total dose would lead to faster control, which is consistent with findings from previous studies. TCP analysis can aid in the selection of treatment protocols by highlighting techniques or protocols that are more likely to lead to faster treatment success.
Abstract
Main Text
1. INTRODUCTION
2. MATERIALS AND METHODS
2.1. Gathering data from the literature
2.2. TCP calculation
3. RESULTS AND DISCUSSION
4. CONCLUSIONS