Abstract

The kinetics of photosensitizer 2-devinyl-2-(1-methoxylethyl) chlorin f (CPD4) uptake in MCF-7 human breast cancer cells is described by a diffusion kinetics model and experimentally investigated using laser scanning confocal microscopy (LSCM). CPD4 permeated into MCF-7 cells with increasing incubation time, which was followed by its binding to cell organelles. Subcellular distribution study revealed that CPD4 was primarily localized on the mitochondria and membranes, supporting that the mode of transmembrane transport was diffusion. A kinetics model describing CPD4 passing through the plasma membrane of MCF-7 cells was proposed based on Fick's first law of diffusion. The kinetics of cellular uptake of CPD4 was studied by three-dimensional LSCM. By fitting the experimental data using the above model, important cellular uptake and distribution parameters were obtained, which are of clinical significance in photodynamic therapy.

1. Introduction

Photodynamic therapy (PDT) is a promising therapeutic technique for the treatment of cancerous or other unwanted tissues based on the combined use of light irradiation and photosensitizer (PS) [19]. After photoactivation of PS, the interaction between the excited triplet PS and ground-state triplet oxygen (3O2) generates singlet oxygen (1O2), which is the main cytotoxic species that can lead to photodynamic damage of tumor cells and tissues [1016]. In comparison with surgical resection followed by radio- and/or chemotherapy, PDT offers the advantage of being able to site-selectively destroy diseased tissues with minimal damage to the surrounding healthy tissues [1721]. These outstanding features pave the way for clinical application of PDT in many developed or developing countries.

Photodynamic damage can be only induced in the close vicinity of 1O2, because 1O2 in cells has a short lifetime (<40 ns) and a limited diffusion in space (<20 nm) [22]. This indicates that the efficacy of PDT is highly dependent on the amount and distribution of intracellular PS [23]. To evaluate the photodynamic activity of a certain PS, two properties need to be addressed. The first is the intracellular accumulation level of 1O2, which relies on both the concentration of PS and the quantum yield of 1O2 production. The second is the subcellular distribution of PS, which is a crucial parameter influencing the extent of cell damage, as well as cell response and the mechanism of the resultant cell death. From a kinetic point of view, the binding between PS and organelles could also lead to loading more PS into the cells by means of maintaining a concentration gradient beneficial to transport of PS from the outside to the inside of cells and meanwhile suppressing the reverse process.

Previous investigations have shown that the partition of PS among subcellular structures is most likely dependent on the dynamics of its transport through membranes [2428]. For example, photofrin II permeated into cells by diffusion and then localized on the mitochondria, while mono-𝐿-aspartyl chlorin 𝑒6 was internalized into cells via endocytosis and then accumulated on the lysosomes [25]. Thus, revealing the subcellular distribution of PS is useful for determining the transmembrane transport modes, which provides the basis for the development of physical models toward quantitative understanding of the transport process. The use of physical models to obtain novel information about biological processes has been reported in the literature. For example, to estimate the threshold dose of 1O2 for tissue necrosis, Stratonnikov et al. proposed a model to describe the oxygen consumption and photobleaching in whole blood samples incubated with PS [29]. Quantitative analyses of molecular motions in cells were also successfully carried out based on the use of diffusion models [30, 31]. In this contribution, we aim at establishing a model to quantitatively describe the transmembrane transport process of PS.

Most of the PSs currently used in research and clinics are derived from the tetrapyrrole aromatic ring. First-generation PSs include hematoporphyrin derivative (HpD), while second-generation PSs are composed of metalloporphyrins, porphycenes, pheophorbides, purpurins, chlorins, protoporphyrin, and phthalocyanines [3]. Chlorins belong to the second-generation PDT sensitizers and are characterized by the presence of one reduced peripheral double bond of the tetrapyrrole ring, leading to a strong absorption of red light (𝜆max=650–680 nm), where tissue transparency is optimal [32]. Prominent features of chlorins also include efficient triplet formation at high energy to allow 1O2 production, absence of dark toxicity, and selective retention in the selected target tissue combined with rapid clearance from the body [33]. Chlorophyll 𝑎, 𝑏, and 𝑐 are the most abundant among the naturally occurring tetrapyrrolic pigments and serve as the starting materials for the preparation of various chlorins [34]. The PS used in this work is 2-devinyl-2-(1-methoxylethyl) chlorin f (hereafter abbreviated as CPD4, Figure 1), which is a chlorophyll derivative synthesized from the extract of silkworm excrement [35]. Attractive properties of CPD4 have been disclosed by spectral measurements, 1O2 quantum yield determination, and photoinduced damage studies [3638]. In this paper, the process of CPD4 uptake in MCF-7 human breast cancer cells is described by a diffusion kinetics model and experimentally investigated using laser scanning confocal microscopy (LSCM).

2. Experimental Section

2.1. Photosensitizer and Cell Culture

CPD4 was synthesized from the extract of silkworm feces with the purity greater than 98%. 10−4 M of CPD4 stock solution was prepared by dissolving CPD4 powder in 0.1 M NaOH solution. MCF-7 human breast cancer cells were offered by the Cancer Hospital of Tianjin Medical University, Tianjin, China. MCF-7 cells were grown in RPMI 1640 supplemented with 10% fetal bovine serum (FBS) and 1% penicillin/streptomycin (Gibco BRL, USA) and were cultured under the atmosphere of 5% CO2 and 95% air at 37°C.

2.2. Subcellular Distribution of CPD4 in MCF-7 Cells

4×104 MCF-7 cells were plated into two 35 mm cell culture dishes (Corning, USA). CPD4 was added to each dish to a final concentration of 4.4×106 M and incubated with cells for 2 hours. Rhodamine 123 (Sigma, USA) was added to a dish to a final concentration of 25 nM and incubated for 15 min; DiI-7 (Sigma, USA) was added to a second dish to a final concentration of 20 nM and incubated for 30 min. Using LSCM (Olympus FV300, Olympus Co., Tokyo, Japan), the samples were excited by 488 nm Ar+ laser (for CPD4 and Rhodamine 123) and 543 nm He-Ne laser (for CPD4 and DiI-7), respectively. Combinations of dichroic mirrors and narrow bandpass filters were used for the detection of fluorescence. The 40x objective was used to visualize features smaller than 30 μm. LSCM images were scanned at 1.65 s/frame and each image contained 512×512 pixels.

2.3. Kinetics of CPD4 Uptake in MCF-7 Cells

4×104 MCF-7 cells were plated into a 35 mm cell culture dish. CPD4 was added to the dish to a final concentration of 4.4×106 M. Fluorescence images of CPD4 at 37°C were taken with LSCM under 488 nm Ar+ laser excitation. The 60x objective was used to get the best image quality at high magnification. The laser beam was scanned in the lateral (𝑥-𝑦) plane generating 2D optically sliced images comprised of 512×512 pixels with an excellent depth resolution. The objective was moved up the z-axis after each lateral plane scanning to obtain a series of 2D images with an increment Δ𝑧=0.5μm (−10 μm ≤ 𝑧 ≤ 20 μm). LSCM images were scanned at 1.65 s/frame. The time interval between two adjacent images was 15 min.

3. Results and Discussion

3.1. Subcellular Distribution of CPD4 in MCF-7 Cells

It is well known that the photodynamic efficacy and the cell death mechanism are closely associated with the localization sites of PS in cells, while the subcellular distribution of PS is determined by its transmembrane transport mode [27]. Therefore, identifying the localization sites of CPD4 in MCF-7 cells will help us to elucidate the transport mode. Since CPD4 is fluorescent, confocal fluorescence microscopy techniques coupled with organelle-specific markers permit the determination of the subcellular distribution with high spatial resolution. In this work, the subcellular distribution of CPD4 in MCF-7 cells is studied by LSCM using two fluorescent probes—the mitochondrial probe Rhodamine 123 and the membrane probe DiI-7.

The fluorescence images of CPD4 (Figures 2(a), 2(d)), Rhodamine 123 (Figure 2(b)), and DiI-7 (Figure 2(e)) in MCF-7 cells have been assigned red, green, and orange false colors, respectively. The overlaid image of CPD4 and Rhodamine 123 shows a yellowish-green color (Figure 2(c)), while that of CPD4 and DiI-7 shows a reddish-orange color (Figure 2(f)). It can be clearly seen that there are very nice fluorescence emission overlaps between CPD4 and Rhodamine 123 and also between CPD4 and DiI-7. To further examine the degrees of these overlaps, a quantitative analysis was carried out by the calculation of correlation coefficient 𝑟 [39], which is a measure of the linear relationship between two variables: 𝑟=0 indicates no linear relationship between the two variables, while 𝑟=1 indicates a perfect linear relationship. The 𝑟 values for CPD4 with Rhodamine 123 and with DiI-7 are 0.897 and 0.627, respectively, which support that the major localization sites of CPD4 in MCF-7 cells are on the mitochondria and membranes. It is noteworthy that similar experiments using Lyso Tracker dyes as the fluorescent probes were also done. Results confirmed that the lysosomes are not the target for CPD4 binding. The details of these observations will be reported elsewhere, together with other results. It has been pointed out that the cellular uptake of PS via diffusion transport generally leads to its accumulation on mitochondria and membranes [25, 26]. Our present study supports the passive diffusion transport mode of CPD4. This result also implies that the binding between CPD4 and organelles helps to maintain a concentration gradient necessary for the continuous loading of PS into MCF-7 cells, which becomes a basic consideration when trying to establish a model to quantitatively describe the transmembrane transport process.

3.2. Diffusion Model of Cellular Uptake of CPD4

The above subcellular distribution study reveals that CPD4 permeates into the cells mostly by diffusion along the concentration gradient, the rates of which are determined by the diffusion coefficient. A practical way to obtain the diffusion coefficient and other parameters is to establish a diffusion model and then fit experimental data using this model. Here we propose a model based on Fick’s first law of diffusion.

To facilitate understanding of the model, Figure 3 presents a cartoon illustrating the diffusion of CPD4 from culture medium to the inside of a cell. In this figure, CPD4 diffuses into the cells along the x-axis, where the diffusion length 𝑙 is equal to the thickness of the plasma membrane. Because the amount of CPD4 permeated into the cells is considerably smaller than that of CPD4 in the culture medium, we assume that the concentration of CPD4 in the culture medium, defined here as 𝐶𝐸(=4.4×106 M), will not change during the whole diffusion process. Therefore, in our diffusion model, the total concentration of CPD4 inside a cell can be written as 𝐶total(𝑥𝑙,𝑡)=𝐶free(𝑡)+𝐶bound(𝑡),(1) where 𝐶free(𝑡) and 𝐶bound(𝑡) are the concentrations of intracellular free and organelle-bound CPD4, respectively. Although practically the distribution of bound CPD4 inside a cell is not even, 𝐶bound(𝑡) refers to the mean concentration here for simplicity. In this equation, 𝑥𝑙 refers to CPD4 inside the cells as depicted in Figure 3. When CPD4 permeates into the cells from the culture medium, we suppose that it immediately binds organelles, thus becoming bound CPD4. During this process, 𝐶free(𝑡)=0 and 𝐶total(𝑡)=𝐶bound(𝑡). As more and more CPD4 permeates into the cells, we assume, at time 𝑡𝑆, the binding between CPD4 and organelles reaches a saturated state. We define the CPD4 concentration at this saturated binding state as 𝐶𝑆, and thus we have 𝐶free(𝑡𝑡𝑆)=0 and 𝐶total(𝑡𝑆)=𝐶𝑆.

In the case of 𝑡𝑡𝑆, according to Fick’s first law of diffusion [40], the diffusion equation of the system can be written as𝐽(𝑥,𝑡)=𝐷𝜕𝐶(𝑥,𝑡)𝜕𝑥,(2) where 𝐽 is the flux, which represents the number of CPD4 passing through unit area of the plasma membrane within unit time, 𝐷 is the diffusion coefficient, and 𝜕𝐶(𝑥,𝑡)/𝜕𝑥 is the concentration gradient along the 𝑥 direction. At 𝑡=0, 𝐶total(𝑥𝑙,𝑡=0)=0 and the concentration gradient of CPD4 across the membrane has a maximal value: 𝜕𝐶(𝑥,𝑡)𝐶𝜕𝑥=𝐸𝐶free(𝑡)𝑙𝐶=𝐸𝑙.(3)

Substituting (3) into (2), we get𝐶𝐽(𝑥=𝑙,𝑡)=𝐷𝐸𝑙.(4)

Since in our experiments the cancer cells adhere to the surface of the culture dish, we assume that, for simplicity, only the upper half of the plasma membrane can accept the permeation of CPD4, while the close contact between the lower half of the plasma membrane and the culture dish blocks CPD4 from entering into the cells. Therefore, based on (4), the total number of CPD4 permeated into a cell at the time 𝑡 can be written as𝐴𝑁(𝑡)=𝐽(𝑥=𝑙,𝑡)2𝐶𝑡=𝐷𝐸𝑙𝐴2𝑡,(5) where 𝐴 is the total area of the cell membrane. Assuming Δ𝑥=𝑉/𝐴, where 𝑉 is the average volume of the cell, from the above deductions we get𝐶total(𝑡)=𝐶bound(𝑡)=𝑁(𝑡)𝑉𝐶=𝐷𝐸2𝑙Δ𝑥𝑡,𝑡𝑡𝑆.(6)

In the case of 𝑡𝑡𝑆, all binding sites in the organelles have been occupied by CPD4, and the concentration of free CPD4 inside the cells begins to increase from zero. Thus, we have 𝐶bound(𝑡)=𝐶𝑆 and 𝐶free(𝑡)=𝐶total(𝑡)𝐶𝑆. If the time interval (Δ𝑡) is considerably short, we can neglect the change of the concentration gradient of CPD4 across the plasma membrane during Δ𝑡. Therefore, we have𝜕𝐶(𝑥,𝑡+Δ𝑡)𝜕𝑥𝜕𝐶(𝑥,𝑡)𝐶𝜕𝑥=𝐸𝐶free(𝑡)𝑙𝐶=𝐸𝐶total(𝑡)𝐶𝑆𝑙.(7)

Substituting (7) into (2), we can express Δ𝐶total(𝑡,Δ𝑡), the concentration increment of intracellular CPD4 during Δ𝑡, asΔ𝐶total(𝑡,Δ𝑡)=Δ𝑁(𝑡,Δ𝑡)𝑉=𝐶𝐷𝐸𝐶total(𝑡)𝐶𝑆/𝑙(𝐴/2)Δ𝑡=𝐷𝐶𝐴Δ𝑥𝐸𝐶total(𝑡)𝐶𝑆Δ𝑡.2𝑙Δ𝑥(8)

𝐶total(𝑡+Δ𝑡) can be calculated by using recursion as shown in the following:𝐶total(𝑡+Δ𝑡)=𝐶total(𝑡)+Δ𝐶total=(𝑡,Δ𝑡)𝐷Δ𝑡𝐶2𝑙Δ𝑥𝐸+𝐶𝑆+1𝐷Δ𝑡𝐶2𝑙Δ𝑥total(𝑡),𝑡𝑡𝑆.(9)

Combining (6) and (9), we now get the equations for the diffusion kinetics model, which can be used to calculate the diffusion coefficient and other parameters:𝐶total(𝐷𝐶𝑡)=𝐸2𝑙Δ𝑥𝑡,𝑡𝑡𝑆𝐷Δ𝑡𝐶2𝑙Δ𝑥𝐸+𝐶𝑆+1𝐷Δ𝑡𝐶2𝑙Δ𝑥total(𝑡Δ𝑡),𝑡Δ𝑡𝑡𝑆.(10)

In the next section, we will carry out the experimental study of the cellular uptake kinetics of CPD4 and fit the data using (10).

3.3. Kinetics of CPD4 Uptake in MCF-7 Cells

Imaging the distributions of CPD4 in MCF-7 cells at different incubation times is done using three-dimensional LSCM. Figure 4 shows the observed images, which are recorded on the same focal plane of 𝑧=0. The variation of the amount of CPD4 is proportional to the change of its relative fluorescence intensities [41]. Figure 4(a) clearly shows that CPD4 accumulates on the plasma membrane of MCF-7 cells 5 min after adding CPD4. As the incubation time increases from 15 to 45 min, CPD4 keeps on permeating into the cells and migrates from the plasma membrane to various sensitive stores within the cells such as the mitochondria and the nucleus membranes (Figures 4(b)4(d)). From 60 to 90 min, the amount of CPD4 localized on these cell organelles increases obviously (Figures 4(e)4(g)). After 90 minutes, the fluorescence intensities of CPD4 in these images remain unchanged indicating that a saturated binding state has been reached (Figures 4(h)-4(i)).

The total concentration of CPD4 at time t can be determined by fluorescence measurement:𝐶total𝐼(𝑡)=total(𝑡)𝐼𝐸𝐶𝐸,(11) where 𝐼total(𝑡) and 𝐼𝐸 represent the average of relative fluorescence intensities of CPD4 inside the cells and in the culture medium, respectively, and can be obtained by summing the fluorescence intensities of each corresponding pixel in all scanned images. 𝐶𝐸(=4.4×106 M) is known; thus 𝐶total(𝑡) at different incubation time 𝑡 can be calculated using (11) and the results are shown in Figure 5.

Using (10), cellular uptake and distribution parameters of CPD4 can be acquired by fitting the data given in Figure 5. In the fitting, 𝑙=7.5 nm [42], Δ𝑥=3.33μm, and Δ𝑡=1 s are used. Fitting procedures are repeated iteratively until the degree of fitting exceeds 0.99. The fitting curve is also shown in Figure 5. From the fitting, we obtain that 𝐷 is 8.11×1016 m2/s, 𝐶𝑆 is 4.018×104 M, and the CPD4 concentration at the diffusion equilibrium state inside the cells is 4.062×104 M. These results indicate that the concentration of intracellular CPD4 is much higher than that of CPD4 in the culture medium (4.4×106 M). The number of binding sites in an MCF-7 cell can be calculated: 𝑁=𝐶𝑆𝑉=1.01×109. It is obvious that the rates of cellular uptake of CPD4 before and after reaching the saturated binding state are quite different: the former is equal to the slope of the fitting curve, 4.3×106 M/min (=𝐶𝑆/𝑡𝑆, from which we can also calculate the time of reaching the saturated binding state: 𝑡𝑆=93.7 min); the latter is 5.5×107 M/min, which is attributed to the concentration difference of free CPD4 on both sides of the plasma membrane. It is noteworthy that the fitting curve shows a smooth curvature (the inset of Figure 5), which indicates the choice of Δ𝑡 is appropriate. Taken together, these results suggest that (1) the diffusion kinetics model established here can be used to quantitatively describe the transport process of CPD4 through the plasma membrane of MCF-7 cells and (2) a large number of CPD4 can bind inside MCF-7 cells, which is an important factor in achieving optimal effect of PDT.

4. Conclusions

In summary, the kinetics of photosensitizer CPD4 uptake in MCF-7 human breast cancer cells is described by a diffusion kinetics model and experimentally investigated using LSCM. CPD4 permeated into MCF-7 cells with increasing incubation time, which was followed by the binding between CPD4 and organelles. Subcellular distribution study revealed that CPD4 was primarily localized on the mitochondria and cell membranes, supporting that the mode of transmembrane transport was diffusion. A kinetics model describing CPD4 passing through the plasma membrane of MCF-7 cells was proposed based on Fick’s first law of diffusion. The kinetics of cellular uptake of CPD4 was studied by three-dimensional LSCM. By fitting fluorescence intensity data using the above model, important cellular uptake and distribution parameters were obtained: the diffusion coefficient of CPD4 across the plasma membrane was 8.11×1016 m2/s; the concentration of organelle-bound CPD4 at the saturated binding state was 4.018×104 M; the concentration of CPD4 at the diffusion equilibrium state inside the cells was 4.062×104 M; the average number of CPD4 binding sites in an MCF-7 cell is 1.01×109; the rates of cellular uptake of CPD4 before and after reaching the saturated binding state were 4.3×106 and 5.5×107 M/min, respectively. These results indicate that the diffusion kinetics model established in this work can be used to obtain quantitative information about the cellular transport process of PS, which will contribute to the better understanding of PDT mechanisms and the development of new PS with optimized interactions with target cells and tissues. Future works that are being planned include investigating the physiochemical properties of CPD4 and performing concentration-dependent uptake experiments using CPD4 as well as other compounds with specific known properties.

Authors’ Contribution

P. Chen and J.-Z. Yao contributed equally to this work.

Acknowledgments

The authors gratefully acknowledge Professor Tong and Mr. Qiao (Cancer Hospital of Tianjin Medical University, Tianjin, China) for offering MCF-7 human breast cancer cells. This research was supported by National Natural Science Foundation of China (nos. 60508004, 60778043, 30371737, and 81172950), the National High Technology Research and Development Program (“863” Program) of China (no. 2011AA030205), and Tianjin Municipal Science and Technology Commission (no. 08ZCDFGX09400).