2145 - Measurements and Modeling of Cancer Cell Population Dynamics at Diagnosis, Treatment and Follow-Up
Presenter(s)
S. Li1, C. T. Miyamoto1, J. Price1, B. Wang1, T. G. Giaddui1, N. Burbure1, S. E. Weiss2, and E. M. Horwitz2; 1Department of Radiation Oncology, Fox Chase Cancer Center at Temple University Hospital, Philadelphia, PA, 2Department of Radiation Oncology, Fox Chase Cancer Center, Philadelphia, PA
Purpose/Objective(s): Non-genetic phenotypic heterogeneity had lately challenged the classic view of genetic determinacy in cell growing and dying processes. Here, we simplified in-vivo cell dynamics with a logistic cancer growing and stable normal-tissue cell model pretreatment, a unified multi-activation model for the cell-killing during treatment and a growing and arresting subpopulation model at the follow-up.
Materials/Methods: Coupled ODEs of cancer cells and sugar glucose co-existing with stable normal-tissue cells results a logistic growth of cell number: N (t) = NoK/[No-(K-No)exp(-rt)] with the initial number No, carry capacity of K and growth coefficient of r measurable from three times’ pretreatment images which can predict cancer cell number NX at treatment time. During the relative short course of treatment, cell subpopulations satisfying a set of ODEs of dN/dD = -cNT and dT/dD = dcNT with T being the integral toxic concentration (e.g. free radicals) per each of total m-fractions, resulting survived cell number at EOT: NE =NX[n/[1+(n-1)exp(D/Do)]m with n and Do estimated from measured cell survival curves. After treatment, a surviving cancer cell would either actively grow as Ng ?Ng + Ng or become an arrested cell as Ng ? Na so that coupled ODEs of dNg /dt = (g– a)Ng and dNa /dt = a*Ng would predict the total number of cancer cell at the follow-up NF = Na+Ng = NE*{a/(g-a) exp[g-a)tF] - a/(g-a)} with the arrested rate, a, and growing rate, g, measured with two or more times’ follow-up data
Results: Some model values to various tumors with diagnostic PET images are listed in the table. We had found (1) fast growing GTV with high PET SUV would had NX approaching to a stable value of K and the GTV quickly shrunk on the late treatment or follow up images; (2) Late-stage CTV had non-uniform cancer cell density, although GTVs of the H&N and rectal masses were removed prior to irradiation, the residual cancel cells grew slowly with dominant normal tissue or arrested cancer cells and response slowly to radiation; and (3) follow-up images had already shown recurrence after conventional fractionated RT with insufficient dose for NE >1 and SBRT of the small metastatic nodule possibly due to early arrested cancer cells that are persistent to radiation.
Conclusion: Dynamic models of cell populations can reasonably explain and predict the clinical outcomes which would be further verified with AI data mining.
Abstract 2145 - Table 1| Tumor | Dose/Fx (Gy) | total Fx | No | K | 1/r (day) | n | Do (Gy) | NX* | NE |
| H&N mass | 2 | 35 | 3.2E+09 | 3.5E+09 | 200 | 4 | 3 | 2.8E+09 | 19.3 |
| Lung mass | 3 | 15 | 4.6E+08 | 5.2E+08 | 150 | 3.2 | 3 | 5.2E+08 | 4317.3 |
| Lung nod. | 18 | 3 | 2.0E+06 | 3.0E+06 | 14 | 3.2 | 3 | 3.0E+06 | 0.1 |
| Prostate | 7.4 | 5 | 2.4E+06 | 2.3E+06 | -623 | 2.8 | 2.2 | 3.3E+06 | 0.9 |
| Mets nod. | 5 | 10 | 1.4E+06 | 1.6E+06 | 66 | 2.2 | 2.8 | 1.6E+06 | 3.3 |
| Rectum | 2 | 27 | 6.5E+08 | 6.4E+08 | 312 | 15 | 3.2 | 3.3E+08 | 35.7 |