JBO Letters

Mathematical model that describes the transition from thermal to photochemical damage in retinal pigment epithelial cell culture

[+] Author Affiliations
Clifton D. Clark, Michael L. Denton

TASC Inc., Biomedical Sciences and Technologies Department, Brooks City-Base, Texas 78235

Robert J. Thomas

Air Force Research Laboratory, 711 HPW/RHDO, Brooks City-Base, Texas 78235

J. Biomed. Opt. 16(2), 020504 (February 17, 2011). doi:10.1117/1.3544504
History: Received October 29, 2010; Accepted December 30, 2010; Revised December 28, 2010; Published February 17, 2011; Online February 17, 2011
Text Size: A A A

Open Access Open Access

* Address all correspondence to: Robert J. Thomas, Air Force Research Laboratory, AFRL/RHDO, 2624 Louis Bauer Drive, Brooks City-Base, Texas 78235-5128. Tel: 200 536 6558; Fax: 210 536 6558; E-mail: Robert.Thomas@brooks.af.mil.

We propose a rate process model for describing photochemical damage to retinal cells by short wavelength laser exposures. The rate equation for photochemical damage contains a positive rate that is temperature independent, and a negative (quenching) rate that is temperature dependent. Using the traditional Arrhenius integral to describe thermal damage, we derive damage threshold doses for both thermal and photochemical mechanisms, and show that the model accounts for the sharp transition from thermal to photochemical damage thresholds that have recently been observed in an in-vitro retinal model.

Figures in this Article

For long pulse durations (

>
1 s), laser induced tissue damage can be divided into two broad categories: thermal and nonthermal (photochemical). Thermal damage has been modeled by the Arrhenius damage integral, first proposed by Henriques1 in 1947, with great success.25 Other models have been developed to address the shortcomings of the Arrhenius integral,67 but it continues to be a popular method for predicting thermal damage thresholds due to its relative simplicity.

No such damage model exists for predicting photochemical damage thresholds. Perhaps the demand for such a model has been low because photochemical damage is defined by the principle of reciprocity, where the total radiant exposure remains constant with exposure duration, allowing damage to be predicted by simply considering the total energy delivered.8 However, in his in-vitro retinal model, Denton et al. showed that at 413 nm, damage thresholds follow a thermal trend line for exposure durations up to 60 s and can be a factor of 2 higher than the reciprocity radiant exposure.9 From this work, the damage threshold mechanism seems to transition sharply from thermal to photochemical somewhere between 60 and 100 s, and a similar trend has been described in an animal model as well.10

In this paper, we propose a rate equation to describe photochemical damage that includes a temperature dependence. The model assumes the creation of some activated molecule (B*) through a photon absorption that is independent of temperature and a temperature dependent quenching rate that removes B* (Fig. 1). We treat thermal and photochemical damage mechanisms as independent processes, the damage threshold is defined as the lowest threshold of the two. The temperature independence of B* creation and the temperature dependence of quenching leads to a sharp cutoff for photochemical damage. Above some critical temperature, the photochemical damage mechanism is completely negated, allowing the damage threshold to be determined by thermal damage at energies much higher than the reciprocity energy.

Graphic Jump LocationF1 :

Photochemical damage mechanism. A single photon (γ) is absorbed by B to produce B* and some other molecule ξ reacts with B* to regenerate B, or a non-reactive product of B. Any accumulation of B* can be available to trigger a chain reaction resulting in cell death (dashed line)

Thermal Damage Rate

We assume thermal damage to be described by the usual Arrhenius rate equation, where some molecule A (typically assumed to be a protein) is converted to a denatured state A*. The incident light provides a heat source which raises the temperature of the cells, increasing the rate of this reaction. The damage rate is given as Display Formula

1d[A*]dt=Z1eEa1/RTt.
Here the brackets, [ · ], denote molar concentrations, and R and T(t) are the gas constant and system temperature (in Kelvin), respectively. The two constants, Ea1 and Z1, are termed activation energy and collision frequency. These rate coefficients must be experimentally determined, and are strongly dependent on the system and damage endpoint assessed. Several sets of coefficients exist in the literature.11

Photochemical Damage Rate

We propose a two-process rate equation to describe photochemical damage. Some molecule B absorbs a single photon to create an activated molecule B*. This activated molecule may then be deactivated through an interaction with some other molecule, ξ, provided by the cell. The activated molecule, B*, may also initiate a chain of events that eventually leads to cell death (Fig. 1).

The total rate of creation for B* can then be written as a positive rate dependent only on the incident photon flux plus a negative rate that accounts for the chemical reaction-based removal of B*, Display Formula

2d[B*]dt=εΦ[B]Z2eEa2/RTt[ξ][B*].
Here, Φ is the incident photon flux density (number of photons per area per time), ε denotes an efficiency factor, which is the probability that an incident photon is absorbed and creates the activated molecule, B*, and Ea2 and Z2 are activation energy and collision frequency associated with ξ and B* reacting to give B. We will assume that the molecules B and ξ are abundant, so their concentrations can be considered constant. We then use the simplified rate equation for B* Display Formula
3d[B*]dt=εΦZ2eEa2/RTt[B*].
We assume photochemical damage to be caused by a buildup of B*, and so the damage threshold will correlate to a concentration of B*. The relationships between positive and negative reaction rates, the temperature independence of damage creation, and the temperature dependence of damage quenching, are unique to the photochemical mechanism.

We seek to find an expression for the threshold radiant exposure for causing damage to cells. In our analysis here, the thermal and photochemical damage mechanisms compete. Thermal damage is caused when A* reaches some threshold concentration, [A*]th, and photochemical damage is caused when B* reaches a threshold concentration, [B*]th. Both thermal and photochemical mechanisms will have a threshold radiant exposure to cause damage. The two rate equations, 1,3, are integrated to give [A*] and [B*] as a function of time and then solved for the radiant exposure (H) that gives [A*]th or [B*]th for the exposure duration τ. The lower of these two thresholds will be the observed damage threshold for the exposure.

In the experiments conducted by Denton et al. ,9 retinal pigment epithelial (RPE) cells were exposed to 413 nm laser radiation for durations of 0.1, 1.0, 20, 40, 60, 100, and 200 s. For exposures greater than about 20 s, the system reached steady-state early in the exposure, and the temperature can be assumed constant throughout. The transition from thermal to photochemical damage occurred somewhere between 60 and 100 s, which we seek to describe, so we can assume the exposure temperature to be constant.

Due to the linearity of the heat equation, the peak steady-state temperature can be written as a linear function of the peak incident power, E0, Display Formula

4T=mE0+T0,
where T0 is the initial (ambient or body) temperature. The slope of this line, m, will depend on the specific laser parameters and system configuration.

Thermal Damage Threshold

Integrating Eq. 1 with the temperature held constant gives Display Formula

5[A*]=Z1eEa1/RTτ,
where τ is the exposure duration and T is the exposure temperature. The threshold temperature, Tth, is the temperature that leads to the threshold damage accumulation, [A*]th. However, the exposure temperature is related to the exposure irradiance through Eq. 4. Inserting T = mEth + T0 into Eq. 5 and solving for Eth gives the damage threshold irradiance. Most rate coefficients (Z1 and Ea1 ) for thermal damage are measured assuming that [A*]th = 1, so the damage threshold irradiance is Display Formula
6H th =τE th =τmEa1RlnZ1τT0.

Photochemical Damage Threshold

To integrate Eq. 3, we define

K2TZ2eEa2/RT
and multiply by the integration factor eK2t which gives, Display Formula
7[B*]=εΦK21eK2Tτ.
We must define [B*]th to obtain a damage threshold exposure. For this, we turn to the principle of reciprocity. Here, the radiant exposure remains constant as the exposure duration varies. It follows that as the exposure duration increases the incident irradiance is reduced, which in turn implies that the exposure temperature is lowered. Reciprocity suggests that the quenching rate, K2(T), is small for low temperature rises.

We define the reciprocity radiant exposure as Hr. The photon flux density can be written in terms of the radiant exposure as Φ = Hhν, where hν is the photon energy and h is Planck's constant. At reciprocity then, Eq. 7 can be written as, Display Formula

8[B*] th =εK2THrτhν1eK2Tτ.
Since K2(T) is small here, we may expand the exponential in a power series and ignore all second order and higher terms ((K2τ)2 ≪ 1) and Eq. 8 simplifies toDisplay Formula
9[B*] th =εHrhν.
Using this definition, we can solve Eq. 7 for the threshold radiant exposure, Display Formula
10H th =HrK2Tτ1eK2Tτ1.
This gives the damage threshold radiant exposure for some fixed temperature. However, the exposure temperature is also dependent on H through Eq. 4 so Eq. 10 becomes Display Formula
11H th =HrK2mH th τ+T0τ1eK2[mH th /τ+T0]τ1.
We use a numerical search algorithm to determine the radiant exposure, Hth, that satisfies this equation.

Measurements for Z1 and Ea1 exist in the literature, but no such measurements exist for Z2 and Ea2. However, we can use the critical temperature (the temperature for which the Arrhenius rate is equal to one11), Tc = Ea2/Rln Z2, to define the relationship between Z2 and Ea2. We previously showed9 that the expected temperature for the 100 s threshold exposure was around 38°C, while the expected temperature for the 60 s threshold exposure was about 48°C. This is in-line with the rule of thumb that thermal damage requires about a 10°C temperature rise. Here, we assume that the critical temperature associated with the photochemical rate process is 47°C (320.15 K).

We use the thermal damage rate coefficients reported by Welch and Polhamus for thermal damage:12Ea1 = 6.28 × 105J mol−1, Z1 = 3.1 × 1099s−1. The expected temperature rises reported by Denton et al. 9 give m = 0.26, and the reciprocity energy was 940 J cm−2. We have computed damage threshold predictions for three different values of Z2: 1080, 10100, and 10120. With Ea2 = 320.15Rln Z2, this gives 4.9 × 105, 6.13 × 105, and 7.35 × 105 for Ea2. The model predicts the same sharp transition from thermal to photochemical damage as was observed (Fig. 2), but the transition region depends on the value of Z2. The smallest collision frequency (Z2 = 1080) gives photochemical damage at exposure duration less than 60 s, while the largest (Z2 = 10120) leads to no photochemical damage at 100 s. The best agreement was with Z2 = 10100.

Graphic Jump LocationF2 :

Predicted temporal action profile (line) compared to the measured damage thresholds by Denton et al. 9 (points). The inset shows photochemical damage predictions for three different collision frequencies: 1080 (top), 10100 (middle), and 10120 (bottom).

We have presented a theoretical model that assumes a temperature dependent photochemical damage rate and provides an explanation for the sharp transition between thermal and photochemical damage recently observed in in-vitro and in-vivo data for short wavelength damage thresholds. As the exposure duration is decreased, the irradiance must be increased to deliver the same number of photons. But the increase in irradiance leads to a larger temperature rise, which eventually disables the photochemical damage. We have not attempted to identify the underlying pathways responsible for causing photochemical damage; the model only assumes that the damage rate can be described by a rate equation. A discussion of possible photochemical mechanisms at the molecular level can be found in Ref. 13.

The damage threshold predictions reproduce the sharp transition between thermal and photochemical damage observed in the in-vitro data reported by Denton et al. 9 The two rate coefficients for photochemical damage quenching have been estimated, but these would need to be determined experimentally to make quantitative predictions. In principle, these coefficients would need to be determined for each tissue and endpoint to be described, but could then be used to predict whether or not photochemical damage is caused under various exposure conditions, similar to how the thermal damage Arrhenius coefficients are used today.

Henriques  F. C., “ Studies of thermal injury v the predictability and the significance of thermally induced rate process leading to irreversible epidermal injury. ,” Arch. Pathol.. 43, , 489–502  ((1947)).
Mainster  M. A., , While  T. J., , Tips  J. H., , and Wilson  P. W., “ Transient thermal behavior in biological systems. ,” Bull. Math. Biol.. 32, , 302–315  ((1970)).
Clarke  A. M., , Geeraets  W. J., , and Ham  W. T.  Jr., “ An equilibrium thermal model for retinal injury from optical sources. ,” Appl. Opt.. 8, (5 ), 1051–1054  ((1969)).
Takata  A. N., , Zaneveld  L., , and Richter  W., “ Laser-induced thermal damage of skin. ,” Technical report, IIT Research Institute ((1977)).
Chen  B., “ Histological and modeling study of skin thermal injury to 2.0 micron laser irradiation. ,” Lasers Surg. Med.. 40, , 358–370  ((2007)).
Feng  Y., , Oden  J. T., , and Rylande  M. N., “ A two-state cell damage model under hyperthermic conditions: theory and in vitro experiments. ,” J. Biomech. Eng.. 130, , 041016  ((2008)).
Till  S. J., , Till  J., , Milsom  P. K., and Rowlands  G., “ A new model for laser-induced thermal damage in the retina. ,” Bull. Math. Biol.. 65, , 731–746  ((2003)).
Wu  J., , Seregard  S., , and Algvere  P. V., “ Photochemical damage of the retina. ,” Surv. Ophthalmol.. 51, , 461–481  ((2006)).
Denton  M. L., , Clark  C. D.  III, , Foltz  M. S., , Schuster  K. J., , Noojin  G. D., , Estlack  L. E., , and Thomas  R. J., “ In vitro retinal model reveals a sharp transition between laser damage mechanisms. ,” J. Biomed. Opt.. 15, , 030532  ((2010)).
Lund  D. J., , Stuck  B., , and Edsall  P., “ Retinal injury threshold for blue wavelength lasers. ,” Health Phys.. 90, , 477–484  ((2006)).
Welch  A., and van Germert  M.,  Optical-Thermal Response of Laser-Irradiated Tissue. ,  Plenum Press ,  New York  ((1995)).
Welch  A., and Polhamus  G., “ Measurement and prediction of thermal injury in the retina of the rhesus monkey. ,” IEEE Trans. Biomed. Eng.. 31, , 633–644  ((1984)).
Gerstman  B. S., and Glickman  R. D., “ Activated rate processes and a specific biochemical mechanism for explaining delayed laser induced thermal damage to the retina. ,” J. Biomed. Opt.. 4, , 345–351  ((1999)).
© 2011 Society of Photo-Optical Instrumentation Engineers (SPIE)

Citation

Clifton D. Clark, III ; Michael L. Denton and Robert J. Thomas
"Mathematical model that describes the transition from thermal to photochemical damage in retinal pigment epithelial cell culture", J. Biomed. Opt. 16(2), 020504 (February 17, 2011). ; http://dx.doi.org/10.1117/1.3544504


Figures

Graphic Jump LocationF1 :

Photochemical damage mechanism. A single photon (γ) is absorbed by B to produce B* and some other molecule ξ reacts with B* to regenerate B, or a non-reactive product of B. Any accumulation of B* can be available to trigger a chain reaction resulting in cell death (dashed line)

Graphic Jump LocationF2 :

Predicted temporal action profile (line) compared to the measured damage thresholds by Denton et al. 9 (points). The inset shows photochemical damage predictions for three different collision frequencies: 1080 (top), 10100 (middle), and 10120 (bottom).

Tables

References

Henriques  F. C., “ Studies of thermal injury v the predictability and the significance of thermally induced rate process leading to irreversible epidermal injury. ,” Arch. Pathol.. 43, , 489–502  ((1947)).
Mainster  M. A., , While  T. J., , Tips  J. H., , and Wilson  P. W., “ Transient thermal behavior in biological systems. ,” Bull. Math. Biol.. 32, , 302–315  ((1970)).
Clarke  A. M., , Geeraets  W. J., , and Ham  W. T.  Jr., “ An equilibrium thermal model for retinal injury from optical sources. ,” Appl. Opt.. 8, (5 ), 1051–1054  ((1969)).
Takata  A. N., , Zaneveld  L., , and Richter  W., “ Laser-induced thermal damage of skin. ,” Technical report, IIT Research Institute ((1977)).
Chen  B., “ Histological and modeling study of skin thermal injury to 2.0 micron laser irradiation. ,” Lasers Surg. Med.. 40, , 358–370  ((2007)).
Feng  Y., , Oden  J. T., , and Rylande  M. N., “ A two-state cell damage model under hyperthermic conditions: theory and in vitro experiments. ,” J. Biomech. Eng.. 130, , 041016  ((2008)).
Till  S. J., , Till  J., , Milsom  P. K., and Rowlands  G., “ A new model for laser-induced thermal damage in the retina. ,” Bull. Math. Biol.. 65, , 731–746  ((2003)).
Wu  J., , Seregard  S., , and Algvere  P. V., “ Photochemical damage of the retina. ,” Surv. Ophthalmol.. 51, , 461–481  ((2006)).
Denton  M. L., , Clark  C. D.  III, , Foltz  M. S., , Schuster  K. J., , Noojin  G. D., , Estlack  L. E., , and Thomas  R. J., “ In vitro retinal model reveals a sharp transition between laser damage mechanisms. ,” J. Biomed. Opt.. 15, , 030532  ((2010)).
Lund  D. J., , Stuck  B., , and Edsall  P., “ Retinal injury threshold for blue wavelength lasers. ,” Health Phys.. 90, , 477–484  ((2006)).
Welch  A., and van Germert  M.,  Optical-Thermal Response of Laser-Irradiated Tissue. ,  Plenum Press ,  New York  ((1995)).
Welch  A., and Polhamus  G., “ Measurement and prediction of thermal injury in the retina of the rhesus monkey. ,” IEEE Trans. Biomed. Eng.. 31, , 633–644  ((1984)).
Gerstman  B. S., and Glickman  R. D., “ Activated rate processes and a specific biochemical mechanism for explaining delayed laser induced thermal damage to the retina. ,” J. Biomed. Opt.. 4, , 345–351  ((1999)).

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging & repositioning the boxes below.

Related Book Chapters

Topic Collections

PubMed Articles
Advertisement


 

  • Don't have an account?
  • Subscribe to the SPIE Digital Library
  • Create a FREE account to sign up for Digital Library content alerts and gain access to institutional subscriptions remotely.
Access This Article
Sign in or Create a personal account to Buy this article ($20 for members, $25 for non-members).
Access This Proceeding
Sign in or Create a personal account to Buy this article ($15 for members, $18 for non-members).
Access This Chapter

Access to SPIE eBooks is limited to subscribing institutions and is not available as part of a personal subscription. Print or electronic versions of individual SPIE books may be purchased via SPIE.org.