- Research
- Open access
- Published:
A local effect model-based interpolation framework for experimental nanoparticle radiosensitisation data
Cancer Nanotechnology volume 8, Article number: 1 (2017)
Abstract
A local effect model (LEM)-based framework capable of interpolating nanoparticle-enhanced photon-irradiated clonogenic cell survival fraction measurements as a function of nanoparticle concentration was developed and experimentally benchmarked for gold nanoparticle (AuNP)-doped bovine aortic endothelial cells (BAECs) under superficial kilovoltage X-ray irradiation. For three different superficial kilovoltage X-ray spectra, the BAEC survival fraction response was predicted for two different AuNP concentrations and compared to experimental data. The ability of the developed framework to predict the cell survival fraction trends is analysed and discussed. This developed framework is intended to fill in the existing gaps of individual cell line response as a function of NP concentration under photon irradiation and assist the scientific community in planning future pre-clinical trials of high Z nanoparticle-enhanced photon radiotherapy.
Background
Photon radiotherapy has undergone significant evolution with the development of new technologies and increased understanding of radiobiology (Mayles et al. 2007; Joiner and van der Kogel 2009). Over the last 15 years, one of the most promising refinements of this cancer treatment modality has been the development and functionalisation of high Z nanoparticles to target cancerous small animals/humans cell lines (Hainfeld et al. 2004, 2008; Jain et al. 2011). This class of novel nanomedicines, of which gold nanoparticles (AuNP) are the most popular (Jain et al. 2012), is thought to increase the local energy deposition and, in-turn, water radiolysis free-radical yield with a few 10–100 nms surrounding each NP (Jones et al. 2010; McMahon et al. 2011; Lechtman et al. 2013; Lin et al. 2014; Sicard-Roselli et al. 2014; Tran et al. 2016). Whilst this technology is still in development and its exact biological action pathway is under intensive investigation, it has already been shown that NP radiosensitising agents utilised in conjunction with radiotherapy are able to provide increased tumour control and life expectancy in small animal models (Hainfeld et al. 2004, 2013; Joh et al. 2013; Xing et al. 2013).
Development and experimental testing of functionalisation high Z NP radiosensitisers for a given cell line is a complex process which can take significant time and resources. Over the last decade, the scientific community has shifted towards exploring the potential of a developed high Z NP radiosensitiser for photon radiotherapy through mechanistic characterisation utilising radiation transport codes such as EGSnrc (Kawrakow 2000), Geant4/Geant4-DNA (Agostinelli et al. 2003; Allison et al. 2006, 2016; Incerti et al. 2010; Bernal et al. 2015), MCNPX (Pelowitz 2005) and PENELOPE (Baro et al. 1995; Salvat et al. 2006). Originally, the scientific community tried to predict the increased effect of high Z NPs through the use of a variety of dose enhancement figures of merit (DEFM) known via a number of different names. All of these DEFMs were based on the assumption that expected biological outcome of cells/tumours could be described via the ratio of dose deposition with and without high Z NP doping under uniform photon irradiation (Cho 2005; Roeske et al. 2007; Ngwa et al. 2010). This underlying assumption neglects two of the key physical factors which determine the action of high Z NP within cells under photon irradiation: (1) the increased localised energy deposition within the first few 10–100 nms of the NP surface (Jones et al. 2010; McMahon et al. 2011; Lechtman et al. 2013; Lin et al. 2014; Sicard-Roselli et al. 2014; Tran et al. 2016), and (2) NP distribution within the irradiated cells (Lechtman et al. 2013; Brun et al. 2009; Coulter et al. 2012; Cui et al. 2014; McQuaid et al. 2016). An alternative to these DEFMs, the local effect model (LEM) (Scholz and Kraft 1996, 2004) was first applied 5 years ago to photon radiotherapy in an attempt to account for one of these two key physical factors: the increased dose localisation within the first few 10–100 nm of the NP surface (McMahon et al. 2011). Two years later, Lechtman et al. (2013) proposed an extension specifically for AuNPs, the AuNP radiosensitisation predictive (ARP) model, in an attempt to account for both of these two physical factors neglected via DEFMs (Lechtman et al. 2013). Both these models were shown to be able to predict specific cell survival fraction behaviour under photon irradiation observed through clonogenic assay (McMahon et al. 2011; Lechtman et al. 2013).
The following work builds on the success of the LEM and presents a new experimentally benchmarked framework capable of interpolating NP-enhanced photon-irradiated clonogenic cell survival fraction measurements as a function of NP concentration. This LEM-based framework was developed to fill in the existing gaps of individual cell line response as a function of NP concentration under photon irradiation to assist the scientific community in planning future pre-clinical trials of high Z nanoparticle-enhanced photon radiotherapy.
Local effect model-based interpolation framework
The developed LEM-based interpolation framework is intended to be used in conjunction with the existing wealth of available experimental survival fraction data for high Z NP-undoped and NP-doped specific cell line studies (Jain et al. 2012). At a minimum each of these studies possesses a set of in vitro clonogenic assays of a cell line undoped and doped with high Z NPs that have been irradiated by a gamma-/X-ray source with a known energy spectra. The following derivation outlines how these data can be interpolated as a function of NP concentration, up to a maximum concentration corresponding to the NP-doped cell line survival data, within the LEM formalism for a given cell line/incident photon energy spectra combination.
The LEM can be constructed utilising three main assumptions. First, the survival fraction of a cellular colony/system under photon irradiation (SF) can be described via a linear-quadratic response:
where \(\alpha\) and \(\beta\) are characteristics of the target cell line, and D is the mean dose delivered to the entire volume of the cellular colony/system (McMahon et al. 2011; Douglas and Fowler 1976). Second, that cell “inactivation”, e.g. cell death, can be attributed to the creation of a number of lethal lesions within a sensitive small sub-cellular volume such as the cell nucleus (Scholz and Kraft 1996, 2004). Here, a lethal lesion is defined as the local modification of DNA generated from the direct and indirect action of ionisation radiation (i.e. a double-strand break). And finally, any contribution of sub-lethal damage at distances larger than the order of a few microns is ignored as it is assumed that there is no interaction between distant sites (Scholz and Kraft 1996, 2004).
Using these assumptions, it is possible to describe the survival fraction for a cell under photon irradiation in terms of the mean number of lethal lesions (\(\langle N (D) \rangle\)):
and inversely:
Within each cell under photon irradiation, lethal lesions are generated inhomogeneously and the probability of their creation is a direct function of local dose deposition. These properties mean that total lesion number in a cell’s sensitive region can be given via integration over its whole volume:
where \({{d}}(x,y,z)\) is the local dose deposited for a given position within the sensitive region of the cell and \(V_{{\text{sens}}}\) is the total volume of the sensitive region of interest.
For a cellular colony/system doped with a concentration of high Z NPs (C), the LEM framework allows for the total local dose deposition within the sensitive region of the cell to be separated into two parts:
where \({{d}}_{\rm U}(x,y,z)\) and \({{d}}_{{\text{NP}}}(C,x,y,z)\) are the dose distributions generated within the sensitive region from the direct interaction of radiation with the bulk cell and high Z NPs, respectively. With this separation, Eq. 4 can be expressed as:
In addition, over the range of validity of dose in the linear-quadratic model, 1–6 Gy (Joiner and van der Kogel 2009), the probability of two energy deposits within \({{d}}_{\rm U}(x,y,z)\) and \({{d}}_{{\text{NP}}}(C,x,y,z)\) at the same location can be assumed to be negligible. Therefore, their product term in Eq. 6 can be set to zero such that:
where \(\langle N_{U} (D) \rangle\) is the mean number of lethal lesion generated via photon interaction within an undoped cellular region, and \(\langle N_{{\text{NP}}} (C,D) \rangle\) is the mean number of lethal lesion generated via high Z NP action within the doped cellular region. Here, \(\langle N_{{\text{NP}}} (C,D) \rangle\) encompasses the lethal lesion generated from direct photon interaction with NPs, secondary electron generated from photon–cellular medium interaction collisions with NPs, and secondary electron/photons generated from photon–NP interactions collision with other NPs. If the spatial distribution of NP uptake within the cell line remains approximately constant with concentration, then from a mechanistic perspective the mean number of lethal lesions generated from these effects can be scaled with average NP density up to a critical saturation threshold (McKinnon et al. 2016). Under these assumptions, Eq. 7 can be manipulated to yield:
where \(\langle N_{{\text{total}}} (C_{0},D) \rangle\) is the mean number of lethal lesions for a given dose D at a known reference concentration \(C_{0}\). With this, Eq. 7 can be expressed as:
where \(\Delta \alpha = \alpha _{{\text{total}}}(C_{0}) - \alpha _{\rm{U}}\) and \(\Delta \beta = \beta _{{\text{total}}}(C_{0}) - \beta _{\rm{U}}\). The final form of the interpolation framework is then given via the substitution of Eq. 9 into Eq. 2:
Multiple concentration and incident photon spectra experimental benchmarking
Experimental benchmarking of the develop framework was undertaken using the only published multiple concentration and incident photon spectra experimental NP radiosensitisation study; the Ph.D. thesis of Rahman, RMIT University (Australia) (Rahman 2010). Within this thesis the radiosensitisation of 1.9 nm AuNP (Nanoprobes Inc., Yaphank, NY 11980, USA) in Bovine Aortic Endothelial Cells (BAECs) under superficial kilovoltage X-ray was studied as a surrogate model for human endothelial cells. The radiosensitivity of four different AuNP concentrations (0, 0.25, 0.5 and 1.0 mMol/L) was explored in triplicate trials for three different kilovoltage X-ray spectra (80, 100 and 150 kVp) delivered via a superficial X-ray therapy (SXRT) machine (Therapax 3 Series, Pantak Inc., Branford, CT, USA) at the William Buckland Radiotherapy Centre (The Alfred Hospital, Australiaρ) (Rahman 2010). Each of these 12 different cell survival curves were composed of a control and five different dose irradiations that were assessed via a CellTiter 96 AQueous One Solution Cell Proliferation Assay (Promega Corp., Madison, Wisconsin). The mean survival fraction, uncertainty (± cell survival standard deviation) and fitted linear-quadratic response of the control (0 mMol/L) and highest concentration (1 mMol/L) data for all three different incident photon spectra are presented in Fig. 1. Each data set’s linear-quadratic response was fitted using least-squares regression in Python, restricting \(\alpha\) and \(\beta\) to positive values, and their corresponding parameters can be found in Table 1. Further information regarding experimental procedure, AuNP cellular localisation, AuNP cytotoxicity, cell viability, and cell mobility can be found in Rahman’s thesis (Rahman 2010).
The developed interpolation framework was applied to the control and AuNP-doped fitted linear-quadratic parameters contained in Table 1 to predict the BAEC survival fraction response as a function of dose for AuNP concentrations of 0.25 and 0.5 mMol/L for all three different incident photon spectra. Figure 2 presents these predicted data sets in conjunction with the 0.25 and 0.5 mMol/L experimental data from Rahman (2010). Comparison of the predicted response and experimental data sets shows that the developed interpolation framework is able to accurately predict the BAEC survival fraction response to within experimental uncertainties for all dose points in the 100 and 150 kVp data sets. For the 80 kVp data, the predicted survival fraction response is within experimental uncertainty for three data points out of six in both the tested 0.25 and 0.5 mMol/L cases. This poor performance of the developed interpolation framework at 80 kVp can be attributed to the high level of statistical fluctuation within the base 80 kVp experimental data seen in Fig. 1.
Figure 3 presents the percentage difference between the control and highest concentration experimental data sets with respect to their fitted linear-quadratic responses shown in Fig. 1. In this figure, it can be seen that the level of difference in the 80 kVp data exceeds both the 100 and 150 kVp data sets. However, the magnitude of the observed difference in Fig. 2 cannot be explained via Fig. 3 alone. Figure 4 presents the percentage difference of the 0.25 and 0.5 mMol/L experimental data in Fig. 2 with respect to their fitted linear-quadratic responses obtained utilising the same protocols as Table 1. The level of difference in the 80 kVp data again exceeds the 100 and 150 kVp data sets, and their combined respective magnitudes with those seen in Fig. 3 correlate with the observation deviation between the experimental and predicted 80 kVp data seen in Fig. 2. These observations indicate that the performance of the developed interpolation framework is directly dependent on the quality of input data, a characteristic common to many interpolative frameworks.
Discussion
A LEM-based framework capable of interpolating NP-enhanced photon irradiated clonogenic cell survival fraction measurements as a function of NP concentration was developed and experimentally benchmarked for 1.9 nm AuNP-doped BAECs under superficial kilovoltage X-ray irradiation. It was illustrated that the performance of the developed framework is directly dependent on the quality of input experimental data. However, further inspection of the percentage differences between experimental data and their respective fitted linear-quadratic responses shown in Figs. 3 and 4 also illustrates that there are limits to which statistical fluctuation can be suppressed via a linear-quadratic fitting approach. Another observation with respect to linear-quadratic response fitting and the present work is that the resultant \(\alpha\) and \(\beta\) values must be restricted to being positive. Without these restrictions, the predicted survival fraction response would be incorrectly estimated. For example, if either value of \(\alpha _{\text{total}}(C_{0})\) or \(\beta _{\text{total}}(C_{0})\) was negative, it would result in an underestimation of the predicted survival fraction response. Whereas if either value of \(\alpha _{\rm{U}}\) or \(\beta _{\rm{U}}\) was negative, it would result in an overestimation of the predicted survival fraction response. Either of these outcomes in the context of high Z NP-enhanced photon radiotherapy treatment planning is unacceptable as it would pose a significant risk to the patient.
The LEM-based interpolation framework presented in this work was developed to fill in the existing gaps within individual cell line response data as a function of NP concentration under photon irradiation. These interpolated data sets will be used in conjunction with another predictive framework that has been developed at Queen’s University Belfast which expresses the enhanced biological response of NP-doped cells/systems in terms of standard photon radiotherapy dose. These two predictive frameworks form the basis of a novel methodology which is intended to assist the scientific community in planning future pre-clinical trials of high Z NP-enhanced photon radiotherapy. Further work is presently underway to illustrate the potential of these two frameworks in the context of AuNP-enhanced breast cancer MV photon radiotherapy as a medical exemplar.
Conclusion
A LEM-based framework capable of interpolating NP-enhanced photon irradiated clonogenic cell survival fraction measurements as a function of NP concentration was developed and experimentally benchmarked for 1.9 nm AuNP-doped BAECs under superficial kilovoltage X-ray irradiation. For three different superficial kilovoltage X-ray spectra (80, 100 and 150 kVp), the BAEC survival fraction response was predicted for two different AuNP concentrations (0.25 and 0.5 mMol/L). Two of the three predicted spectra data sets (100 and 150 kVp) were within experimental uncertainties for all data points, whereas the other data set (80 kVp) was within experimental uncertainties half of the time. The observed poor performance for the 80 kVp data set was found to be due to a high level of statistical fluctuation within the base data and this illustrated that the performance of developed interpolation framework is directly dependent on the quality of the input experimental data. It is anticipated that this interpolation framework will serve as an important tool for planning future pre-clinical and clinical trials of high Z NP-enhanced photon radiotherapy.
Abbreviations
- ARP Model:
-
gold nanoparticle radiosensitisation predictive model
- AuNP:
-
gold nanoparticle
- BAEC:
-
bovine aortic endothelial cell
- DEFM:
-
dose enhancement figures of merit
- DNA:
-
deoxyribonucleic acid
- LEM:
-
local effect model
- NP:
-
nanoparticle
References
Agostinelli S, et al. GEANT4—a simulation toolkit. Nucl Instrum Methods Phys Res Sect A. 2003;506(3):250–303.
Allison J, et al. Geant4 developments and applications. IEEE Trans Nucl Sci. 2006;53(1):270–8.
Allison J, et al. Recent developments in Geant4. Nucl Instrum Methods Phys Res Sect A. 2016;835:186–225.
Baro J, Sempau J, Fernández-Varea JM, Salvat F. PENELOPE: an algorithm for Monte Carlo simulation of the penetration and energy loss of electrons and positrons in matter. Nucl Instrum Methods Phys Res Sect B. 1995;100(1):31–46.
Bernal MA, et al. Track structure modeling in liquid water: a review of the Geant4-DNA very low energy extension of the Geant4 Monte Carlo simulation toolkit. Phys Med. 2015;31(8):861–74.
Brun E, Sanche L, Sicard-Roselli C. Parameters governing gold nanoparticle X-ray radiosensitization of DNA in solution. Colloids Surf B Biointerfaces. 2009;72(1):128–34.
Cho SH. Estimation of tumour dose enhancement due to gold nanoparticles during typical radiation treatments: a preliminary Monte Carlo study. Phys Med Biol. 2005;50(15):N163.
Coulter JA, Jain S, Butterworth KT, Taggart LE, Dickson GR, McMahon SJ, Hyland WB, Muir MF, Trainor C, Hounsell AR, O’Sullivan JM. Cell type-dependent uptake, localization, and cytotoxicity of 1.9 nm gold nanoparticles. Int J Nanomedicine. 2012;7(1):2673–85.
Cui L, Tse K, Zahedi P, Harding SM, Zafarana G, Jaffray DA, Bristow RG, Allen C. Hypoxia and cellular localization influence the radiosensitizing effect of gold nanoparticles (AuNPs) in breast cancer cells. Radiat Res. 2014;182(5):475–88.
Douglas BG, Fowler JF. The effect of multiple small doses of X-rays on skin reactions in the mouse and a basic interpretation. Radiat Res. 1976;66(2):401–26.
Hainfeld JF, Dilmanian FA, Slatkin DN, Smilowitz HM. Radiotherapy enhancement with gold nanoparticles. J Pharm Pharmacol. 2008;60(8):977–85.
Hainfeld JF, Slatkin DN, Smilowitz HM. The use of gold nanoparticles to enhance radiotherapy in mice. Phys Med Biol. 2004;49(18):N309–15.
Hainfeld JF, Smilowitz HM, O’Connor MJ, Avraham’Dilmanian F, Slatkin DN. Gold nanoparticle imaging and radiotherapy of brain tumors in mice. Nanomedicine. 2013;8(10):1601–9.
Incerti S, et al. Comparison of GEANT4 very low energy cross section models with experimental data in water. Med Phys. 2010;37(9):4692–708.
Jain S, Coulter JA, Hounsell AR, Butterworth KT, McMahon SJ, Hyland WB, Muir MF, Dickson GR, Prise KM, Currell FJ, O’Sullivan JM, Hirst DG. Cell-specific radiosensitization by gold nanoparticles at megavoltage radiation energies. Int J Radiat Oncol Biol Phys. 2011;79(2):531–9.
Jain S, Hirst DG, O’Sullivan JM. Gold nanoparticles as novel agents for cancer therapy. Br J Radiol. 2012;85:101–13.
Joh DY, Sun L, Stangl M, Zaki, Murty S, Santoiemma PP, Davis JJ, Baumann BC, Alonso-Basanta M, Bhang D, Kao GD, Tsourkas A, Dorsey JF. Selective targeting of brain tumors with gold nanoparticle-induced radiosensitization. PLoS ONE. 2013;8(4):e62425.
Joiner M, van der Kogel A. Basic clinical radiobiological. 4th ed. Boca Raton: CRC Press, Taylor & Francis Group; 2009.
Jones BL, Krishnan S, Cho SH. Estimation of microscopic dose enhancement factor around gold nanoparticles by Monte Carlo calculations. Med Phys. 2010;37(7):3809–16.
Kawrakow I. Accurate condensed history Monte Carlo simulation of electron transport. I. EGSnrc, the new EGS4 version. Med Phys. 2000;27(3):485–98.
Lechtman E, Mashouf S, Chattopadhyay N, Keller BM, Lai P, Cai Z, Reilly RM, Pignol JP. A Monte Carlo-based model of gold nanoparticle radiosensitization accounting for increased radiobiological effectiveness. Phys Med Biol. 2013;58(10):3075–87.
Lin Y, McMahon SJ, Scarpelli M, Paganetti H, Schuemann J. Comparing gold nano-particle enhanced radiotherapy with protons, megavoltage photons and kilovoltage photons: a Monte Carlo simulation. Phys Med Biol. 2014;59(24):7675–89.
Mayles P, Nahum A, Rosenwald J-C. Handbook of radiotherapy physics: theory and practice. New York: CRC press, Taylor & Francis Group; 2007.
McKinnon S, Engels E, Tehei M, Konstantinov K, Corde S, Oktaria S, Incerti S, Lerch M, Rosenfeld A, Guatelli S. Study of the effect of ceramic Ta2O5 nanoparticle distribution on cellular dose enhancement in a kilovoltage photon field. Phys Med. 2016;32(10):1216–24.
McMahon SJ, Hyland WB, Muir MF, Coulter JA, Jain S, Butterworth KT, Schettino G, Dickson GR, Hounsell AR, O’Sullivan JM, Prise KM, Hirst DG, Currell FJ. Biological consequences of nanoscale. Sci Rep. 2011;1:18.
McQuaid HN, Muir MF, Taggart LE, McMahon SJ, Coulter JA, Hyland WB, Jain S, Butterworth KT, Schettino G, Prise KM, Hirst DG. Imaging and radiation effects of gold nanoparticles in tumour cells. Sci Rep. 2016;6:19442.
Ngwa W, Makrigiorgos GM, Berbeco RI. Applying gold nanoparticles as tumor-vascular disrupting agents during brachytherapy: estimation of endothelial dose enhancement. Phys Med Biol. 2010;55(21):6533–48.
Pelowitz DB. MCNPX user’s manual version 2.5, vol. 76. Los Alamos: Los Alamos National Laboratory; 2005.
Rahman WNWA. Gold nanoparticles: novel radiobiological dose enhancement studies for radiation therapy, synchrotron based microbeam and stereotactic radiotherapy. Doctoral Dissertation, RMIT University (Australia). 2010.
Roeske JC, Nuñez L, Hoggarth M, Labay E, Weichselbaum RR. Characterization of the theoretical radiation dose enhancement from nanoparticles. Technol Cancer Res Treat. 2007;6(5):395–401.
Salvat F, Fernández-Varea JM, Sempau J. PENELOPE-2006: a code system for Monte Carlo simulation of electron and photon transport. In: Workshop proceedings OECD/NEA, vol. 7. 2006.
Scholz M, Kraft G. Track structure and the calculation of biological effects of heavy charged particles. Adv Space Res. 1996;18(1):5–14.
Scholz M, Kraft G. The physical and radiobiological basis of the local effect model: a response to the commentary by R. Katz. Radiat Res. 2004;161(5):612–20.
Sicard-Roselli C, Brun E, Gilles M, Baldacchino G, Kelsey C, McQuaid H, Polin C, Wardlow N, Currell F. A new mechanism for hydroxyl radical production in irradiated nanoparticle solutions. Small. 2014;10(16):3338–46.
Tran HN, Karamitros M, Ivanchenko VN, Guatelli S, McKinnon S, Murakami K, Sasaki T, Okada S, Bordage MC, Francis Z, El Bitar Z, Bernal MA, Shin JI, Lee SB, Barberet Ph, Tran TT, Brown JMC, Nhan Hao TV, Incerti S. Geant4 Monte Carlo simulation of absorbed dose and radiolysis yields enhancement from a gold nanoparticle under MeV proton irradiation. Nucl Instrum Methods Phys Res Sect B. 2016;373:126–39.
Xing H, Zheng X, Ren Q, Bu W, Ge W, Xiao Q, Zhang S, Wei C, Qu H, Wang Z, Hua Y, Zhou L, Peng W. Computed tomography imaging-guided radiotherapy by targeting upconversion nanocubes with significant imaging and radiosensitization enhancements. Sci Rep. 2013;3:1751.
Authors' contributions
JMCB developed the theory, designed and undertook the analysis, and wrote/edited the manuscript. FJC commissioned the project, mentored JMCB, and edited the manuscript. Both authors read and approved the final manuscript.
Acknowledgements
JMCB acknowledges N. Lampe of LPC Clermont-Ferrand, France for his helpful comments and suggestions.
Competing interests
The authors declare that they have no competing interests.
Data and materials
Supporting data sets are accessible via the Queen’s University Belfast Research Portal (http://pure.qub.ac.uk/portal).
Funding
This work was supported by EPSRC Grant CEP/K039342/1.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
About this article
Cite this article
Brown, J.M.C., Currell, F.J. A local effect model-based interpolation framework for experimental nanoparticle radiosensitisation data. Cancer Nano 8, 1 (2017). https://doi.org/10.1186/s12645-016-0025-6
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12645-016-0025-6