The ultimate goal of cancer treatment utilizing thermotherapy is to eradicate tumors and minimize damage to surrounding host tissues. To achieve this goal, it is important to develop an accurate cell damage model to characterize the population of cell death under various thermal conditions. The traditional Arrhenius model is often used to characterize the damaged cell population under the assumption that the rate of cell damage is proportional to $exp(−Ea∕RT)$, where $Ea$ is the activation energy, $R$ is the universal gas constant, and $T$ is the absolute temperature. However, this model is unable to capture transition phenomena over the entire hyperthermia and ablation temperature range, particularly during the initial stage of heating. Inspired by classical statistical thermodynamic principles, we propose a general two-state model to characterize the entire cell population with two distinct and measurable subpopulations of cells, in which each cell is in one of the two microstates, viable (live) and damaged (dead), respectively. The resulting cell viability can be expressed as $C(τ,T)=exp(−Φ(τ,T)∕kT)∕(1+exp(−Φ(τ,T)∕kT))$, where $k$ is a constant. The in vitro cell viability experiments revealed that the function $Φ(τ,T)$ can be defined as a function that is linear in exposure time $τ$ when the temperature $T$ is fixed, and linear as well in terms of the reciprocal of temperature $T$ when the variable $τ$ is held as constant. To determine parameters in the function $Φ(τ,T)$, we use in vitro cell viability data from the experiments conducted with human prostate cancerous (PC3) and normal (RWPE-1) cells exposed to thermotherapeutic protocols to correlate with the proposed cell damage model. Very good agreement between experimental data and the derived damage model is obtained. In addition, the new two-state model has the advantage that is less sensitive and more robust due to its well behaved model parameters.

1.
Roemer
,
R. B.
, 1999, “
Engineering Aspects of Hyperthermia Therapy
,”
Annu. Rev. Biomed. Eng.
1523-9829,
1
(
1
), pp.
347
376
.
2.
Shinohara
,
K.
, 2004, “
Thermal Ablation of Prostate Diseases: Advantages and Limitations
,”
Int. J. Hyperthermia
0265-6736,
20
(
7
), pp.
679
697
.
3.
Arrhenius
,
S.
, 1889, “
On the Reaction Velocity of the Inversion of Cane Sugar by Acids
,”
Z. Phys. Chem.
0942-9352,
4
, pp.
226
248
.
4.
Laidler
,
K. J.
, 1987,
Chemical Kinetics
, 3rd ed.,
Harper & Row
,
New York
.
5.
Henriques
,
F. C.
, and
Moritz
,
A. R.
, 1947, “
Studies of Thermal Injury. I. The Conduction of Heat to and Through Skin and the Temperatures Attained Therein. A Theoretical and an Experimental Investigation
,”
Am. J. Pathol.
0002-9440,
23
(
53
), pp.
1
549
.
6.
Moritz
,
A. R.
, and
Henriques
,
F. C.
, 1947, “
Studies of Thermal Injury. II. The Relative Importance of Time and Surface Temperature in the Causation of Cutaneous Burns
,”
Am. J. Pathol.
0002-9440,
23
, pp.
695
720
.
7.
He
,
X.
, and
Bischof
,
J. C.
, 2003, “
Quantification of Temperature and Injury Response in Thermal Therapy and Cryosurgery
,”
Crit. Rev. Biomed. Eng.
0278-940X,
31
, pp.
355
422
.
8.
Roti Roti
,
J. L.
, and
Henle
,
K. J.
, 1980, “
Comparison of Two Mathematical Models for Describing Heat-Induced Cell Killing
,”
0033-7587,
81
(
3
), pp.
374
383
.
9.
Lepock
,
J. R.
, 1987, “
How Do Cells Respond to Their Thermal Environment?
,”
Int. J. Hyperthermia
0265-6736,
21
(
8
), pp.
681
687
.
10.
Pearce
,
J.
, and
Thomsen
,
S.
, 1995, “
Rate Process Analysis of Thermal Damage
,”
,
A. J.
Welch
and
M. J. C.
van Gemert
, eds.,
Plenum Press
,
New York
, pp.
561
606
.
11.
Philibert
,
J.
, 2006, “
Some Thoughts and/or Questions About Activation Energy and Pre-Exponential Factor
,”
Defect Diffus. Forum
1012-0386,
249
, pp.
61
72
.
12.
Dewey
,
W. C.
,
Hopwood
,
L. E.
,
Sapareto
,
S. A.
, and
Gerweck
,
L. E.
, 1977, “
Cellular Responses to Combinations of Hyperthermia and Radiation
,”
0033-8419,
123
(
2
), pp.
463
474
.
13.
Rylander
,
M. N.
,
Diller
,
K. R.
,
Wang
,
S.
, and
Aggarwal
,
S. J.
, 2005, “
Correlation of HSP70 Expression and Cell Viability Following Thermal Stimulation of Bovine Aortic Endothelial Cells
,”
ASME J. Biomech. Eng.
0148-0731,
127
, pp.
751
757
.
14.
Rylander
,
M. N.
,
Feng
,
Y.
,
Bass
,
J.
, and
Diller
,
K. R.
, 2005, “
Thermally Induced Injury and Heat-Shock Protein Expression in Cells and Tissues
,”
0077-8923,
1066
, pp.
222
242
.
15.
Tsang
,
W.
,
Bedanov
,
V.
, and
Zachariah
,
M. R.
, 1996, “
Master Equation Analysis of Thermal Activation Reactions: Energy-Transfer Constraints on Falloff Behavior in the Decomposition of Reactive Intermediates With Low Thresholds
,”
J. Phys. Chem.
0022-3654,
100
, pp.
4011
4018
.
16.
Chan
,
H. S.
, and
Dill
,
K. A.
, 2006, “
Transition States and Folding Dynamics of Proteins and Heteropolymers
,”
J. Chem. Phys.
0021-9606,
100
(
12
), pp.
9238
9257
.
17.
Dill
,
K. A.
, 1985, “
Theory for the Folding and Stability of Globular Proteins
,”
Biochemistry
0006-2960,
24
(
6
), pp.
1501
1509
.
18.
He
,
X.
, and
Bischof
,
J. C.
, 2005, “
The Kinetics of Thermal Injury in Human Renal Carcinoma Cells
,”
Ann. Biomed. Eng.
0090-6964,
33
(
4
), pp.
502
510
.
19.
Merlo
,
C.
,
Dill
,
K. A.
, and
Weikl
,
T. R.
, 2005, “
Φ Values in Protein-Folding Kinetics Have Energetic and Structural Components
,”
0027-8424,
102
(
29
), pp.
10171
10175
.
20.
Weikl
,
T. R.
,
Palassini
,
M.
, and
Dill
,
K. A.
, 2004, “
Cooperativity in Two-State Protein Folding Kinetics
,”
Protein Sci.
0961-8368,
13
(
3
), pp.
822
829
.
21.
Moussa
,
N. A.
,
Tell
,
E. N.
, and
Cravalho
,
E. G.
, 1979, “
Time Progression of Hemolysis of Erythrocyte Populations Exposed to Supraphysiological Temperatures
,”
ASME J. Biomech. Eng.
0148-0731,
101
, pp.
213
217
.
22.
Jung
,
H. A.
, 1986, “
A Generalized Concept for Cell Killing by Heat
,”
0033-7587,
106
, pp.
56
72
.
23.
Jung
,
H. A.
, 1991, “
A Generalized Concept for Cell Killing by Heat: Effect of Chronically Induced Thermotolerance
,”
0033-7587,
127
, pp.
235
242
.
24.
Schrödinger
,
E.
, 1989,
Statistical Thermodynamics
,
Dover
,
New York
.
25.
McQuarrie
,
D. A.
, 2000,
Statistical Mechanics
,
University Science Books
,
Sausalito, CA
.
26.
Wright
,
N. T.
, 2003, “
On a Relationship Between the Arrhenius Parameters From Thermal Damage Studies
,”
ASME J. Biomech. Eng.
0148-0731,
125
, pp.
300
304
.