- Open Access
Modelling the correlation between EGFr expression and tumour cell radiosensitivity, and combined treatments of radiation and monoclonal antibody EGFr inhibitors
© Pedicini et al.; licensee BioMed Central Ltd. 2012
- Received: 4 March 2012
- Accepted: 26 May 2012
- Published: 19 June 2012
The Erratum to this article has been published in Theoretical Biology and Medical Modelling 2012 9:37
To estimate the effects of heterogeneity on tumour cell sensitivity to radiotherapy combined with radiosensitizing agents attributable to differences in expression levels of Epidermal Growth Factor Receptor (EGFr).
Materials and methods
Differences in radiosensitivity are not limited to cells of different cancer histotypes but also occur within the same cancer, or appear during radiotherapy if radiosensitizing drugs are combined with ionizing radiation. A modified biologically effective dose (MBED), has been introduced to account for changes in radiosensitivity parameters (α and α/β) rather than changes in dose/fraction or total dose as normally done with standard biologically effective dose (BED). The MBED approach was applied to cases of EGFr over-expression and cases where EGFr inhibitors were combined with radiation. Representative examples in clinical practice were considered.
Assuming membrane EGFr over-expression corresponds to reduced radiosensitivity (α H = 0.15 Gy -1 and α H /β H = 7.5 Gy) relative to normal radiosensitivity (α = 0.2 Gy -1 and α/β = 10 Gy), an increased dose per fraction of 2.42 Gy was obtained through the application of MBED, which is equivalent to the effect of a reference schedule with 30 fractions of 2 Gy. An equivalent hypo-fractionated regime with a dose per fraction of 2.80 Gy is obtained if 25 fractions are set. Dose fractionations modulated according to drug pharmacokinetics are estimated for combined treatments with biological drugs. Soft and strong modulated equivalent hypo-fractionations result from subtraction of 5 or 10 fractions, respectively.
During this computational study, a new radiobiological tool has been introduced. The MBED allows the required dose per fraction to be estimated when tumour radiosensitivity is reduced because EGFr is over-expressed. If radiotherapy treatment is combined with EGFr inhibitors, MBED suggests new treatment strategies, with schedules modulated according to drug pharmacokinetics.
- Epidermal Growth Factor Receptor
- Epidermal Growth Factor Receptor Expression
- Epidermal Growth Factor Receptor Inhibitor
- Dose Fractionation
- Biological Drug
Recently, radiobiology has been transformed thanks to new knowledge concerning cellular activation processes in response to an external stimulus. This knowledge has led to the identification of promising new drug therapies called "targeted therapy”.
Epidermal growth factor receptor (EGFr) has emerged as a central molecular target for modulation during cancer therapy. The correlation between over-expression of EGFr and clinically aggressive malignant disease suggested that EGFr was a promising target for several epithelial tumours, which represent approximately two thirds of all human cancers. Furthermore, the favourable interaction profile for EGFr blocking agents combined with radiation has stimulated clinical trials in diverse anatomical sites including head and neck, colorectal region, pancreas and lung, where molecular inhibition of EGFr signalling in combination with radiation represents a highly promising area[3, 4].
where α and β represent intrinsic and repair cell radiosensitivity, respectively, d represents the dose per fraction and D is the total dose delivered during the radiation treatment. The BED is considered a “biological dose” delivered by a particular combination of dose per fraction and total dose to a given tissue, characterized by a given α/β ratio, and is commonly used to equate or compare various fractionation schedules.
However, eq. (1) demonstrates that the same number of cells killed – the equivalent effect – could be obtained equating the BED not only for schedules with different numbers of fractions and various doses per fraction, but also for schedules where the dose per fraction is increased if a reduction in radiosensitivity results (i.e. α or β is reduced).
This could be applicable for subsets of cells that over-express EGFr, representing a source of heterogeneity closely connected with the repopulation rate and intrinsic radiosensitivity. However, the heterogeneous population of EGFr expression cannot be represented by a single equation of tumour control probability (TCP), as it is intrinsically linked to a group of tumours with identical characteristics.
Furthermore, equations considering the radiation response that take into account different compartments of sensitivity within tumours or a Gaussian distribution of individual radio sensitivities[13, 14] cannot be used because various levels of radiosensitivity coexist in the tumors or in the statistical sample.
Therefore, during this computational study, a new mathematical interpretation of radiosensitivity parameters that are normally used in standard radiobiological models (i.e. as functions of EGFr expression) is proposed using simple examples.
The final aim of the current study is to provide an additional mathematical tool that can be used to carry out radiobiological analysis, taking into account the radioresistance effects due to EGFr over-expression and/or radiosensitization effects due to EGFr inhibitors when they are combined with radiation.
These examples are not intended to simulate a particular type of radiotherapy treatment, but are designed to demonstrate a general effect.
Modified BED: Effects due to a change in EGFr expression levels EGFr expression has been assessed through intensity of staining (i.e., absent, minimal, moderate, or intense staining) in clinical practice. During the present analysis, normal and high expression levels of EGFr (i.e. below and above 50% staining) were distinguished. The subscript H was added to indicate high EGFr expression.
In eq. (3) the LHS provides a measure of treatment effect under standard conditions of radiosensitivity, while the RHS represents the same effect achieved under non standard conditions of radiosensitivity.
Therefore, the MBED distinguishes between changes in biological effect due to irreparable and/or reparable damage variations, rather than changes due to dose/fraction or total dose variations. A reduction in radiosensitivity due to increased membrane EGFr expression[11, 18] implies equivalence between treatments by increasing the dose per fraction with an equal number of fractions.
Modified BED: Effects due to biological drugs
Combined treatment comprising radiation and radiosensitizing EGFr inhibitor drugs requires the daily dose that achieves the same effect without drugs to be calculated. This will result in a calculation of the daily radiosensitivity conditions induced by the drug compared with standard radiosensitivity.
On the basis of a preclinical assessment, we propose a method to estimate the daily radiosensitivity when radiotherapy treatment is combined with biological drugs. Subsequently, the MBED method is applied to assess the changes required in terms of dose fractionation when such daily radiosensitivity is considered.
This solution leads to a modulated hypo-fractionation if the number of weeks is less than the standard fractionation (vice versa for the hyper-fractionation).
Eq. (7) and eq. (9) represent dose values that have the same effect. However, as in the drug is also absorbed by normal tissue cells, these cells will show increased radiosensitivity. Therefore, modulated dose fractionation with a reduced dose of radiation corresponding to higher radiosensitivity could lead to a reduction in harmful effects.
This proposal could be verified through clinical trials.
This section discusses results from representative examples occurring in clinical practice. Schedules with the equivalent effect of 30 fractions of 2 Gy/fraction (assumed as a reference standard regime) were calculated. To analyze an increase in radiosensitivity, a change in α or β, and consequently a change in α/β, has been assumed to simplify the calculations without losing generality.
For examples 3, 4 and 5, substantial changes in β alone has been adopted, assuming that data were obtained from the polynomial regressions of curves depicted in Figure3.
Of note, the unchanged α, β (without polynomial regression) and the fractionation schemes assumed in these examples are plausible but should not be considered as recommendations for real clinical situations.
Numerical results for Examples 1 and 2
d ex1 (Gy)
d ex2 (Gy)
Numerical results for Examples 3, 4 and 5
d ex3 (Gy)
d ex4 (Gy)
d ex5 (Gy)
In this example a selection of patients that should be treated with the reference schedule (consisting of 30 fractions of 2 Gy/fraction on PTV) was assumed. Patients in the first subset (G1) were considered to have normal EGFr expression on clonogenic tumour cells, with radiosensitivity corresponding to α = 0.2 Gy -1 , β = 0.02 Gy -2 (α/β = 10 Gy). In addition, we considered a second subset of patients (G2) as presenting with EGFr cell membrane over-expression, resulting in a reduction of radiosensitivity with α H = 0.15 Gy -1 , β H = 0.02 Gy -2 (α/β H = 7.5 Gy).
with a noticeable reduction in the effect of overall treatment.
The new schedule will be not equivalent in terms of toxicity to organs at risk (OAR). Therefore, the plan will require re evaluation of the harmful effects for OARs. In the opposing situation, that is for an increase of radiosensitivity in the clonogens of G2 compared with G1 (owing to a radiosensitizing drug), one can adopt the same procedure. In such cases, the equivalent effect on the PTV, with the same number of fractions, will be reached by reducing the fraction dose.
For the same subsets of patients used in Example 1, we analyzed a hypo-fractionated schedule that lasted for one week less for patients in G2, with the same effect as the standard schedule for patients in G1. In the hypo-fractionation schedule, the number of fractions was m = 5·(n w -1) = 5·5 = 25 fractions.
In this example we refer to group G2 having substantial membrane EGFr over-expression, with α H = 0.2 Gy -1 and α H /β H = 16 Gy (similar estimated α/β values are reported in the literature). We compare the reference treatment with a combined treatment comprising radiation and biological drugs that produce an increase in radiosensitivity.
In addition, we assume a weekly drug dosage with a pharmacokinetics curve showing maximum absorption during the first day of treatment. The weekly radiosensitivity is assumed to be that described by the set of values reported in Table2.
Examples 4 and 5
The equivalent global effect of the reference schedules could be obtained by subtracting one or two weeks of treatment from eq. (9), with a modulated soft hypo-fractionation (5 weeks) and with a modulated strong hypo-fractionation (4 weeks), respectively. The results are presented in Table2 and Figure5.
During practical applications of radiobiological models, the main difficulty is to decide which parameter values should be included in individual calculations. It is important to clarify that population based estimates of the α/β value represent averages, and that values are likely to vary between and within tumour types. It is clear that the assumption of a single value for α or α/β is a simplification and this could have a considerable impact on the predictive use of BED when deciding on dose fractionation.
However, recent knowledge concerning molecular mechanisms allows new developments to be explored and provides important information relating to the intrinsic radiosensitivity and fractionation sensitivity. Cell studies in vitro demonstrate that differences in radiosensitivity occur among cell lines derived from different types of tumours or from the same type of tumour, and during irradiation when combined treatments using radiation and radiosensitizing drugs are utilised[16, 23–25].
These considerations may lead the way for new studies concerning evaluation of α and β, in which cellular radiosensitivity is modified using known concentrations of radiosensitizing drugs, as described in Figure4 and Figure5.
Therefore, the historical inability to distinguish among effects resulting in differences in radiosensitivity could be overcome through new knowledge concerning heterogeneity[26, 27]. These effects are well known from preclinical studies, and could be used to reduce uncertainties and investigated through clinical trials. The ideal situation could be to use assay methods to allocate patients to various treatment schedules on the basis of individual measurements of tumour cell radiosensitivity (for example, due to varied expression of EGFr) or absorption of drugs. This approach is expected to be applied in the foreseeable future.
On the basis of these considerations, a new method to interpret BED expression, named MBED, was introduced during this computational study to take account of intrinsic differences in radiosensitivity.
The requirement to introduce MBED arises because radiosensitivity is usually considered to be fixed for a cell type and constant during any radiation treatment. For this reason, α and β are considered fixed values with considerable uncertainty. Therefore, in the standard use of the BED, the hypothesis that one fractionation is equivalent to another underlies the assumption that the values of α and β are the same: to have the same effect – resulting in the same number of cells being killed – changing the dose per fraction, one must alter the number of fractions.
Herein, it is argued that for various values of radiosensitivity, the same number of cells can be killed with the same number of fractions by varying the dose per fraction. This requires identification of prognostic parameters such as the over-expression of EGFr, which allows the radiosensitivity of the individual patient to be classified and the most appropriate radiation dose fractionation to be identified.
The results of this study demonstrate that for a subset of patients presenting with EGFr cell membrane over-expression, resulting in reduced radiosensitivity with respect to a subset of patients with normal EGFr expression of clonogenic tumour cells, the dose per fraction should be increased to produce the same therapeutic effect with the same number of fractions taken in the reference treatment.
When radiation is combined with a biological drug that produces an increase in radiosensitivity, depending on the drug dosage, the equivalent treatment with the same number of fractions is obtained by a dose of radiation modulated according to drug pharmacokinetics.
The dose needs to be increased if the number of fractions is reduced.
In the examples reported herein, the absorption of EGFr inhibitors was considered for cancer cells alone. In general, cells of normal tissues also absorb the drug. In particular, EGFr is over-expressed in skin cells. Therefore, the effect of increased radiosensitivity will affect these cells, and modulated fractionations with a lower dose of radiation corresponding to higher radiosensitivity could lead to a reduction of harmful effects.
With MBED, this study was not intended to implement a finely tuned model based on accurate data obtained from preclinical analysis. The aim was to demonstrate the potential of the model and its malleability in terms of including further information that selective preclinical studies may provide.
In addition, previous analyses have depended on the validity of the LQ model, which has limitations. In particular, the LQ model used during this study does not include the time factor. In the generalized LQ model[5, 10] the temporal factor is affected by differences in EGFr expression due to its influence on potential doubling time, T D [29–32]. This temporal factor can be particularly important when the MBED model is used to compare treatment schedules that differ in terms of overall treatment times, tumour control or acute effects (where time dependent repopulation may be important). The difference of doubling time between the High EGFr group and the Low EGFr group identified during the current study will be investigated further in new studies. This difference in terms of T D can be transformed into an equivalent dose that would be required to offset the modified proliferation occurring in one day. The value of this equivalent dose can be taken into account during the previous analysis.
Overall, in practical applications of the MBED concept, there should be careful consideration of the relevant physical dose variations, the possible range of biological parameters and pertinent clinical factors. The prudent clinical oncologist should use MBED as a guide during clinical decisions rather than as an absolute indicator. The advice of acknowledged experts in radiobiological modelling should be sought in more complicated clinical situations.
Despite these limitations, the MBED model provides a valid means of accounting for modulated intrinsic radiosensitivity effects, which is preferable to neglecting them by using a biologically uncorrected physical dose.
Furthermore, the method is not intrinsically associated with the disease, and can be applied to any case by integrating traditional treatment plans and improving the overall radiotherapy performances in combined treatments comprising radiosensitizing drugs.
During this computational study, the MBED method was introduced. The MBED provides a new tool to estimate the effects of heterogeneity on tumour radiosensitivity and to assess the dose per fraction required for increased tumour radiosensitivity due to EGFr over-expression. Where radiotherapy treatment is combined with radiosensitizing drugs, MBED suggests that the fraction sizes modulated according to drug pharmacokinetics will allow new schedules of dose fractionation to be more effective.
In conclusion, the MBED method could improve overall radiotherapy performances and be utilised to perform more appropriate radiobiological analysis, particularly when combined treatment comprising radiation and biological drugs is employed.
This work resulted from collaborative research between the CROB of Rionero in Vulture and IEO of Milan. CROB funded this collaboration.
We thank Cossu Rocca Maria – Senior Medical oncologist IEO Milan - for expert assistance concerning pharmacokinetics.
- Lammering G: Molecular predictor and promising target: will EGFr now become a star in radiotherapy?. Radioth Oncol. 2005, 74: 89-91. 10.1016/j.radonc.2005.01.004.View ArticleGoogle Scholar
- Harari P, Huang SM: Modulation of molecular targets to enhance radiation. Clin Cancer Res. 2000, 6: 323-335.PubMedGoogle Scholar
- Baumann M, Krause M, Dikomey E, Dittmann K, Dörr W, Kasten Pisula U, Rodemann HP: EGFr targeted anticancer drugs in radiotherapy: preclinical evaluation of mechanisms. Radioth Oncol. 2007, 83: 238-248. 10.1016/j.radonc.2007.04.006.View ArticleGoogle Scholar
- Morgan MA, Parsels LA, Kollar LE, Normolle DP, Maybaum J, Lawrence TS: The combination of epidermal growth factor receptor inhibitors with gemcitabine and radiation in pancreatic cancer. Clin Cancer Res. 2008, 14: 5142-5149. 10.1158/1078-0432.CCR-07-4072.PubMed CentralView ArticlePubMedGoogle Scholar
- Fowler JF: The linear quadratic formula and progress in fractionated radiotherapy. Brit J Rad. 1989, 62: 679-694. 10.1259/0007-1285-62-740-679.View ArticleGoogle Scholar
- Niemierko A, Goitein M: Implementation of model for estimating tumour control probability for an inhomogeneously irradiated tumour. Radioth Oncol. 1993, 29: 140-147. 10.1016/0167-8140(93)90239-5.View ArticleGoogle Scholar
- Zaider M, Hanin L: Biologically equivalent dose and long term survival time in radiation treatments. Phys Med Biol. 2007, 52: 6355-6362. 10.1088/0031-9155/52/20/017.View ArticlePubMedGoogle Scholar
- Dubray BM, Thames HD: The clinical significance of ratios of radiobiological parameters. Int J Radiat Oncol Biol Phys. 1996, 35: 1099-1111. 10.1016/0360-3016(96)00255-6.View ArticlePubMedGoogle Scholar
- Barendsen GW: Dose fractionation, dose rate and isoeffect relationships for normal tissue responses. Int J Radiat Oncol Biol Phys. 1982, 8: 779-790.Google Scholar
- Hendry JH, Bentzen SM, Dale RG, Fowler JF, Wheldon TE, Jones B, Munro AJ, Slevin NJ, Robertson AG: A modelled comparison of the effects of using different ways to compensate for missed treatment days in radiotherapy. Clin Oncol. 1996, 8: 297-307. 10.1016/S0936-6555(05)80715-0.View ArticleGoogle Scholar
- Zagars GK, Schultheiss TE, Peters LJ: Inter tumour heterogeneity and radiation dose control curves. Radioth Oncol. 1987, 8: 353-362. 10.1016/S0167-8140(87)80186-X.View ArticleGoogle Scholar
- Eriksen JG, Steiniche T, Overgaard J, Danish Head and Neck Cancer Study group: The role of epidermal growth factor receptor and E-cadherin for the outcome of reduction in the overall treatment time of radiotherapy of supraglottic larynx squamous cell carcinoma. Acta Oncol. 2005, 44: 50-58. 10.1080/02841860510007396.View ArticlePubMedGoogle Scholar
- Bengt KL, Brahme A: The radiation response of heterogeneous tumours. Phys Med. 2007, 23: 91-99. 10.1016/j.ejmp.2007.08.001.View ArticleGoogle Scholar
- Webb S, Nahum AE: A model for calculating tumour control probability in radiotherapy including the effects of inhomogeneous distribution of dose and clonogenic cell density. Phys Med Biol. 1993, 38: 653-666. 10.1088/0031-9155/38/6/001.View ArticlePubMedGoogle Scholar
- Ang KK, Berkey BA, Tu X, Zhang HZ, Katz R, Hammond EH, Fu KK, Milas L: Impact of Epidermal Growth Factor Receptor expression on survival and pattern of relapse in patients with advanced head and neck carcinoma. Cancer Res. 2002, 62: 7350-7356.PubMedGoogle Scholar
- Dittmann K, Mayer C, Rodemann HP: Inhibition of radiation induced EGFr nuclear import by C225 (Cetuximab) suppresses DNA PK activity. Radioth Oncol. 2005, 76: 157-161. 10.1016/j.radonc.2005.06.022.View ArticleGoogle Scholar
- Akashi Y, Okamoto I, Iwasa T, Yoshida T, Hatashita Y, Yamada Y, Satoh T, Fukuoka M, Ono K, Nakagawa K: Enhancement of the antitumour activity of ionising radiation by nimotuzumab, a humanised monoclonal antibody to the epidermal growth factor receptor, in non-small cell lung cancer cell lines of differing epidermal growth factor receptor status. Br J Cancer. 2008, 98: 749-755. 10.1038/sj.bjc.6604222.PubMed CentralView ArticlePubMedGoogle Scholar
- Bonner JA, Maihle NJ, Folven BR, Christianson TJ, Spain K: The interaction of epidermal growth factor and radiation in human head and neck squamous cell carcinoma cell lines with vastly different radiosensitivities. Int J Radiat Oncol Biol Phys. 1994, 29: 243-247. 10.1016/0360-3016(94)90269-0.View ArticlePubMedGoogle Scholar
- Barendsen GW: Differences in radiosensitivity among cells in culture and in experimental tumours: significance for the effectiveness of human cancer therapy. Radioth Oncol. 1987, 8: 285-300. 10.1016/S0167-8140(87)80178-0.View ArticleGoogle Scholar
- Fracasso PM, Burris H, Arquette MA, Govindan R, Gao F, Wright LP, Goodner SA, Greco FA, Ones SF, Willcut N, Chodkiewcz C, Pathak A, Springett GM, Simon GR, Sullivan DM, Marcelpoil R, Mayfield SD, Mauro D, Garrett CR: A Phase 1 Escalating Single Dose and Weekly Fixed Dose Study of Cetuximab: Pharmacokinetic and Pharmacodynamic Rationale for Dosing. Clin Cancer Res. 2007, 13: 986-993. 10.1158/1078-0432.CCR-06-1542.View ArticlePubMedGoogle Scholar
- Joiner MC, van der Kogel AJ: The linear-quadratic approach to fractionation and calculation of isoeffect relationships. Basic clinical radiobiology. Edited by: Steel GG. 1997, London: Arnold, 107-122.Google Scholar
- Barendsen GW: RBE as a function of dose for effects on tissues and tumours assessed by the linear quadratic model. Int J Radiat Oncol Biol Phys. 2000, 46: 684-685.View ArticlePubMedGoogle Scholar
- Carlson DJ, Stewart RD, Li XA, Jennings K, Wang JZ, Guerrero M: Comparison of in vitro and in vivo α/β ratios for prostate cancer. Phys Med Biol. 2004, 49: 4477-4491. 10.1088/0031-9155/49/19/003.View ArticlePubMedGoogle Scholar
- Begg AC, Vens C: Genetic manipulation of radiosensitivity. Int J Radiat Oncol Biol Phys. 2001, 49: 367-371. 10.1016/S0360-3016(00)01497-8.View ArticlePubMedGoogle Scholar
- Brown JM: Therapeutic targets in radiotherapy. Int J Radiat Oncol Biol Phys. 2001, 49: 319-326. 10.1016/S0360-3016(00)01482-6.View ArticlePubMedGoogle Scholar
- Strigari L, D'Andrea M, Abate A, Benassi M: A heterogeneous dose distribution in simultaneous integrated boost: the role of the clonogenic cell density on the tumour control probability. Phys Med Biol. 2008, 53: 5257-5273. 10.1088/0031-9155/53/19/001.View ArticlePubMedGoogle Scholar
- Roberts SA, Hendry JH: A realistic closed-form radiobiological model of clinical tumour-control data incorporating intertumour heterogeneity. Int J Radiat Oncol Biol Phys. 1998, 41: 689-699. 10.1016/S0360-3016(98)00100-X.View ArticlePubMedGoogle Scholar
- Thariat J, Milas L, Ang KK: Integrating radiotherapy with epidermal growth factor receptor antagonists and other molecular therapeutics for the treatment of the head and neck cancer. Int J Radiat Oncol Biol Phys. 2007, 69: 974-984. 10.1016/j.ijrobp.2007.06.076.PubMed CentralView ArticlePubMedGoogle Scholar
- Suwinski R, Jaworska M, Nikiel B, Grzegorz W, Bankowska-Wwozniak M, Wojciech M, Krysztof S, Dariusz L: Predicting the effect of accelerated fractionation in postoperative radiotherapy for head and neck cancer based on molecular marker profiles: data from a randomized clinical trial. Int J Radiat Oncol Biol Phys. 2010, 77: 438-446. 10.1016/j.ijrobp.2009.05.021.View ArticlePubMedGoogle Scholar
- Chung CH, Zhang Q, Hammond EM, Trotti AM, Wang H, Spencer S, Zhang HZ, Cooper J, Jordan R, Rotman MH, Ang KK: Integrating epidermal growth factor receptor assay with clinical parameters improves risk classification for relapse and survival in head-and-neck squamous cell carcinoma. Int J Radiat Oncol Biol Phys. 2010, 81: 331-338.PubMed CentralView ArticlePubMedGoogle Scholar
- Bentzen SM, Atasoy BM, Daley FM, Dische S, Richman PI, Saunders MI, Trott KR, Wilson GD: Epidermal growth factor receptor expression in pretreatment biopsies from head and neck squamous cell carcinoma as a predictive factor for a benefit from accelerated radiation therapy in a randomized controlled trial. J Clin Oncol. 2005, 23: 5560-5567. 10.1200/JCO.2005.06.411.View ArticlePubMedGoogle Scholar
- Smid EJ, Stoter TR, Bloemena E, Lafleur MV, Leemans CR, van der Waal I, Slotman BJ, Langendijk JA: The importance of immunohistochemical expression of EGFr in squamous cell carcinoma of the oral cavity treated with surgery and postoperative radiotherapy. Int J Radiat Oncol Biol Phys. 2006, 65: 1323-1329. 10.1016/j.ijrobp.2006.03.011.View ArticlePubMedGoogle Scholar
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.