In the previous paper , the authors made a bilayer 3D artificial tumor tissue (BxPC3/NHDF3D) using human pancreatic cancer cell line BxPC3 and normal human dermal-derived fibroblast NHDF3D, and pharmacokinetic study for 10BPA-BNCT treatment  -  by neutron irradiation from reactor Went. As a result, optical observation of the irradiated tissue gave a T/N ratio of 3.19 at a neutron dose of 2.52 × 1012 n/cm2, indicating the effectiveness of the BNCT pharmacokinetics test using 3D artificial tissue  - .
On the other hand, what is important in BNCT treatment of tumor patients is to secure the accumulation in the tumor affected area after the boron drug administration and to optimize the irradiation dose during the treatment.
Therefore, in this paper, we examined the applicability of the 3D artificial tumor tissue used in the previous paper  as a testing technique for optimizing the radiation dose at the time of BNCT treatment while verifying the sensitivity performance to radiation dose.
2. Materials and Methods
2.1. Cells, Reagents and Instruments
Normal human dermal-derived fibroblast (NHDFs) and red fluorescent protein (RFP) -labeled human pancreatic cancer cell line BxPC3 used in the experiment were purchased from LONZA (Walkersville, MD) and Anti-Cancer Japan (Ibaraki, Japan), respectively. Dulbecco’s modified Eagle’s medium (DMEM) (Wako, Osaka, Japan) containing 10% fetal bovine serum (FBS) (Nichirei, Tokyo, Japan) was used to proliferate cells prior to construction of the tumor tissue model. The cells were cultivated at 37˚C, 5% carbon dioxide. Bovine plasma-derived fibronectin (FN) and porcine skin gelatin (G) were purchased from Sigma-Aldrich (St. Louis, MO) and Wako Pure Chemical Industries, Ltd. (Osaka, Japan), respectively. Transwell inserts with porous polyester bottom (pore size: 0.4 µm) for 12-well culture plate (12 mm diameter, 112 mm2 area, cat. No. 3401) were purchased from CORNING Inc. (New York, NY). The boron drug 10BPA (boronophenylalanine; C9H12BNO4, molecular weight 209.01, fructose complex) was kindly gifted from Interpharma Praha, a.s. (Komorany, Czechia). A solid-state nuclear track detector CR-39, an optical plastic material with composition C12H18O7  , was purchased from Cokin (Tokyo, Japan) .
2.2. Preparation of in Vitro Human Three-Dimensional Tumor Tissue Model
Human three-dimensional tumor tissue was prepared by cell accumulation method    . First, connective tissue-like structures were fabricated by three-dimensional lamination of NHDFs. As shown in Figure 1, NHDFs are cultured under the conditions mentioned above, then ECM-nano film (about 10 nm thick) was formed on each cell surface by coating cells with fibronectin and gelatin dissolved in Tris-HCl buffer according to previously published methods
Figure 1. Procedures for detection of BNCT reaction by using in vitro three-dimensional artificial human tumor tissue model. (a) Fabrication process of cancer cells-loaded human artificial tissue model. (b)-(e) Operation of BNCT for the tissue models, and detection of cancer cell and nuclear track distribution. In present paper, we adopted combination of an artificial tumor tissue model, comprised of normal human dermal-derived fibroblast (NHDF) and human pancreatic cancer cell line BxPC3.
   . The cells are seeded on the transwell inserts at a density of 27.2 × 105 cells/insert (8 layers) and cultured under the conditions of 5% carbon dioxide at 37˚C for 12 to 24 hours. We regarded this connective tissue-like structure as an artificial human normal tissue model, termed as NHDF3D .
Next, RFP-labeled BxPC3 cells which have been cultured and proliferated were collected by trypsin treatment, washed, and uniformly seeded on the upper surface of NHDF3D at a density of 300 cells/mm2, then further cultured for 24 hours under the above culture conditions. BxPC3 cells proliferated on the surface of NHDF3D forming the groups with a flat shape. We regarded this cancer-loaded NHDF3D as an artificial human tumor tissue model, termed as NHDF3D/BxPC3. In the present study, one case of NHDF3D and three cases of NHDF3D/BxPC3 were prepared by the methods mentioned above and were used for the experiments.
2.3. BPA Immersion Treatment and Fixation
The boronophenylalanine (BPA) solution (3% w/v) was diluted to a concentration of 40 ppm with DMEM containing 10% FBS. This is referred to as a BPA treatment solution. After removing the culture solution of NHDF3D or NHDF3D/BxPC3, 750 μl of BPA treatment solution was added, and incubated for 2 hours (BPA exposure) under conditions of 5% carbon dioxide at 37˚C. After that, the BPA treatment solution was removed and the tissues were washed 3 times with 0.01 M phosphate buffered saline (PBS, pH 7.3). Subsequently, the tissues were fixed by 4% paraformaldehyde/0.1 M phosphate buffer (pH 7.3) for 30 minutes at room temperature shading the light. After the fixation, the cellular nucleus was stained by 4’,6-diamidino-2-phenylindole (DAPI) .
2.4. Observation of BxPC3 cell Distribution on the Artificial Tissues
After the fixation, three holes were provided on the tissue using an 18 G injection needle in order to provide alignment marks. Then the top surface of NHDF3D/BxPC3 was observed by fluorescence microscope BZ-X700 (Keyence, Osaka, Japan). The distribution of BxPC3 cells was visualized as a fluorescence image by RFP excitation. A low magnification image including the entire tissue and high magnification images of various parts of each tissue were respectively obtained. A fluorescence image of NHDF3D was also obtained as a control .
2.5. Neutron Irradiation Experiment
The above-mentioned NHDF3D or NHDF3D/BxPC3 in the transwell inserts were cut out with a knife together with the polyester base, mounted on the solid track detector CR-39   with close contact, and used as a sample for track image acquisition (Figure 1(d)). Irradiation experiments using these samples were conducted at the Heavy Water Neutron Irradiation Facility of Kyoto University Reactor (KUR), and irradiation was performed for 60 minutes under an irradiation flux of 1.4 × 109 n/cm2/s (total flux = 4.5 × 1012 n/cm2). After the neutron irradiation, the above sample was etched (6N NaOH, 70˚C × 2 hours) to visualize the α-ray/recoiled Li particle tracks generated on the CR-39 surface .
2.6. Analysis of α-Ray/Recoiled Li Particle Track Images
The α-ray/recoiled Li particle track image of etched CR-39 was taken using the bright field function of the fluorescence microscope BZ-X700. A low magnification image including the entire tissue mount and 10 random high magnification images were obtained respectively, and the following analysis was performed.
1) The whole tissue images of fluorescent BxPC3 cell distribution and the α-ray/ recoiled Li particle track distribution were compared referring to the position of three-hole markers, and the relationship of these distributions was observed.
2) Alpha-ray/recoiled Li particle tracks in the high magnification images were regarded as the particles, and quantitatively analyzed by using software FIJI (https://fiji.sc). Briefly, after binarizing the image, the tracks with more than 4 μm diameter were detected, and their number per unit area (0.01 mm2) and the track size were measured. A significant difference between the data from the artificial normal tissue model (NHDF3D) and artificial tumor tissue model (NHDF3D/BxCP3) was confirmed using student t test .
3. Results and Discussions
3.1. Distribution of the BxPC3 Cells and Alpha-Ray/Recoiled Li Particle Tracks on the Tissue Models
Low-magnification images of NHDF3D and NHDF3D/BxPC3 (three cases) detecting RFP fluorescence were shown in Figure 2(a) (NHDF3D) and Figures 2(b)-(d) (NHDF3D/BxPC3) before irradiation. Although the BxPC3 cells (Figures 2(b)-(d)) were seeded over the artificial tissue (Figure 2(a)) in regular cell concentration, slight heterogeneity of cell distribution was observed after the cultivation for 24 hours.
The α-ray/recoiled Li particle tracks in corresponded parts of Figures 2(a)-(d) were shown in Figures 3(a)-(d). The marked #1 - #12 in the figure are measurement points that have been observed by high-magnification optical observation in Figure 3.
Figure 4 group (a) (3DNHDF) and group (b) (NHDF3D/BxPC3), and Figure 4 group (c) (d) (NHDF3D/BxPC3) showed the distribution of α-ray/recoiled Li particle tracks in 3DNHDF and NHDF3D/BxPC3, respectively.
In these figures, the α-ray/recoiled Li particle tracks were observed as small dots at #1 - #12, reflecting the intensity of the emission. These tracks distribution and density at different observation points in each tissue are approximated for each group. This means that the α-ray/recoiled Li particle damage occurs almost uniformly in each tissue.
Since the contrast between NHDF3D (Figure 4(a)) and NHDF3D/BxPC3
Figure 2. Comparision of Fluorescent staining of NHDF3D and BxPC3 cells. (a)-(d) Microscopic images of NHDF3D and BxPC3 cells on culture dish. (a) NHDF3D. (b)-(d) red fluorescence image of RFP-labeled BxPC3 cells.
Figure 3. Comparison of α-ray/recoiled Li particle track distribution in low magnification. (a) NHDF3D as a control without BxPC3 cells. (b)-(d) three cases of cancer-loaded tissue; NHDF3D/BxPC3 (1 - 3). The marked #1 to #12 in the figure are measurement points that has been observed by high-magnification optical observation in Figure 4.
3.2. Radiation Sensitivity of NHDF3D and NHDF3D/BxCP3
Here, we verified the suitability as a tool to determine the dose optimization during BNCT treatment by using NHDF3D/BxPC3.
To investigate detailed α-ray/recoiled Li particle track density in NHDF3D and NHDF3D/BxCP3 samples, high magnification images of α-ray/recoiled Li particle tracks on CR-39 [Figure 5(A) Figure 5(B) and Figure 5(a’) Figure 5(b’), photo] were quantitatively analyzed. Here, Figure 5(a’) Figure 5(b’) correspond to Figure 4(a) and Figure 4(c), respectively. Figure 5(A) and Figure 5(B) show α-ray/recoiled Li particle tracks observed in NHDF3D and NHDF3D/BxPC3, when irradiated with Case 1 (2.52 × 1012)  and Figure 5(a’) Figure 5(b’) with Case 2 (4.5 × 1012 n/cm2 neutron dose), respectively. Table 1 shows the results of measuring the density of α-ray/recoiled Li particle tracks in the unit area (0.01 mm2) in Case 1 and Case 2.
There is a clear difference in track density between the two, and it can be seen that the higher the neutron irradiation amount, the higher the track density. Similarly, when the track density of NHDF3D/BxPC3 was compared between Figure 5(B) and Figure 5(b’), it was found that the track density was high at the high dose NHDF3D/BxPC3.
Figure 4. Comparison of NHDF3D and BxPC3 cell distribution and α-ray/recoiled Li particle track distribution in higher magnification. Groupe (a) and (b): High-magnification optical observation α-ray/recoiled Li particle track distribution at observation points #1 - 6. Groupe (a): NHDF3D. Groupe (b): NHDF3D/BxPC3 (1); Groupe (c) (d): High-magnification optical observation α-ray/recoiled Li particle track distribution at observation points #7 - 12. Groupe (c) and Groupe (d): NHDF3D/BxPC3 (2 - 3).
Table 1. α-ray/recoiled Li particle tracks density (counts/0.01 mm2), and T/N ratio of NHDF3D/BxPC3.
Figure 5. Comparison of α-ray/recoiled Li particle tracks distribution in NHDF3D and NHDF3D/BxCP3 samples. High magnification photo images of α-ray/recoiled Li particle track distribution were processed by using software Fiji for quantitative analysis of α-ray/recoiled Li particle tracks. (A): NHDF3D of Case 1 and (a’): NHDF3D of Case 2. (B): NHDF3D/BxPC3 of Case 1 and (b’) NHDF3D/BxPC3 (2) of Case 2. The binary data were prepared from the photo images and particle analysis was performed to count the number of the tracks and measure the of their size and the number of the particle tracks in 0.01 mm2 (track density). By further expanding the tissue change sites in the normal tissue and tumor tissue in the artificial tumor tissue after BNCT treatment, it was possible to grasp in detail the α track shape and distribution that occurred after BNCT treatment in Case 1 and Case 2.
The track density increases as the irradiation dose increases in this way, meaning that 10B absorbed in the sample remains in the tissue in an unreacted state even after neutron irradiation for a short time of about 30 minutes .
Figure 6 shows the relationship between neutron irradiation dose (n/cm2) and track density (counts/0.01 mm2). According to this, the increase rate of NHDF3D and BxPC3 is different with the increase of the irradiation amount, and the increase rate of the track density of NHDF3D is remarkable. On the other hand, it can be seen that the increase rate of BxPC3 is almost saturated, however, in Case 2, the rate of increase is similar to that of normal cells.
This means that when neutrons dose reached up to 4.25 × 1012 n/cm2 or more, the effect of radiation damage on normal cells must be considered in addition to the original BNCT treatment effect.
Figure 7 shows the relationship between the neutron irradiation dose and the T/N ratio, where the definition of T/N ratio here is defined as the track density ratio of NHDF3D and BxPC3 in Table 1. According to this, although the T/N ratio increases as the neutron irradiation dose increases, it turns to decrease when the irradiation dose in Case 1 is exceeded. In other words, in Case 2 from the viewpoint of BNCT treatment, radiation damage to the normal tissue becomes significant and there is a possibility of over-irradiation.
Figure 6. The relationship between irradiation dose and track density. The increase rate of NHDF3D and BxPC3 is different with the increase of the irradiation amount, and the increase rate of the track density of NHDF3D is remarkable. On the other hand, it can be seen that the increase rate of BxPC3 is almost saturated, however, in Case 2, the rate of increase is similar to that of normal cells.
Figure 7. The relationship between irradiation dose and T/N ratio. The T/N ratio increases as the neutron irradiation dose increases, it turns to decrease when the irradiation dose in Case 1 is exceeded. In other words, in Case 2 from the viewpoint of BNCT treatment, radiation damage to the normal tissue becomes significant and there is a possibility of over-irradiation.
In conclusion, our in vitro model of tumor tissue for BNCT demonstrated the pharmacokinetics of BPA and the efficacy of neutron irradiation by direct observation of α-ray/recoiled Li particle tracks that are corresponding to the distribution of BxPC3 cells . Furthermore, in this paper, following the previous paper, heavy irradiation was performed using the same artificial tumor tissue. As a result, it became clear that selection of the optimal dose, which is one of the most important treatment conditions for BNCT treatment, is possible to some extent by an evaluation method using an artificial tumor tissue.
 Ishiyama, Asano, Y., Suzuki, M., et al. (2019) In Vitro Evaluation System of Pharmacokinetics and Irradiation Effect in Boron Neutron Capture Therapy (BNCT) Using Three-Dimentional Artificial Human Tumor Tissue Model. Journal of Cancer Therapy, 10, 835-845.
 Masunaga, S., Sakurai, Y., Tanaka, H., et al. (2014) The Dependency of Compound Biological Effectiveness Factors on the Type and Concentration of Administered Neutron Capture Agents in Boron Neutron Capture Therapy. Springerplus, 3, 128.
 Ishiyama, S., Baba, Y., Fujii, R., Nakamura, M. and Imahori, Y. (2012) Synthesis of Lithium Nitride for Neutron Production Target of BNCT by In-Situ Lithium Deposition and Ion Implantation. Nuclear Instruments and Methods in Physics Research, 293, 42-47.
 Ishiyama, S., Baba, Y., Fujii, R., Nakamura, M. and Imahori, Y. (2013) Thermal Stability of BNC Neutron Production Target Lithium Synthesized by In-Situ Lithium Deposition and Ion Implantation. Materials Transactions, 54, 1760-1764.
 Ishiyama, S. and Imahori, Y. (2014) Deterministic Parsing Model of CBE Factor for Intra-Organ 10B Distribution in Boron Neutron Capture Therapy. Journal of Cancer Therapy, 5, 1388-1398.
 Ishiyama, S., Imahori, Y., Itami, J. and Hanna, V. (2015) Determination of the Compound Biological Effectiveness (CBE) Factors Based on the Ishiyama-Imahori Deterministic Parsing Model with the Dynamic PET Technique. Journal of Cancer Therapy, 6, 759-766.
 Suzuki, M., et al. (2007) First Attempt of Boron Neutron Capture Therapy (BNCT) for Hepatocellular Carcinoma. Japanese Journal of Clinical Oncology, 37, 376-381.
 Ichihashi, M., et al. (1982) Specific Killing Effect of 10B1-Para-boronophenylalanine in Thermal Neutron Capture Therapy of Malignant Melanoma, in Vitro Radiobiological Evaluation. Journal of Investigative Dermatology, 78, 215-218.
 Hiratsuka, J., et al. (1982) RBEs of Thermal Neutron Capture Therapy and 10B(n, alpha)7Li Reaction on Melanoma Bearing Hamsters. Pigment Cell Research, 2, 352-355.
 Matsusaki, M. (2012) Development of Three-Dimensional Tissue Models Based on Hierarchical Cell Manipulation Using Nanofilms. Bulletin of the Chemical Society of Japan, 85, 401-414.
 Matsusaki, M., Case, C.P. and Akashi, M. (2014) Three-Dimensional Cell Culture Technique and Pathophysiology. Advanced Drug Delivery Reviews, 74, 95-103.
 Asano, Y., Odagiri, T., Oikiri, H., Matsusaki, M., Alashi, M. and Shimoda, H. (2017) Construction of Artificial Human Peritoneal Tissue by Cell-Accumulation Technique and Its Application for Visualizing Morphological Dynamics of Cancer Peritoneal Metastasis. Biochemical and Biophysical Research Communications, 494, 213-219.
 Nishiguchi, A., Matsusaki, M., Kano, M.R., Nishihara, H., Okano, D., Asano, Y., Shimoda, H., Kishimoto, S., Iwai, S. and Akashi, M. (2018) In Vitro 3D Blood/ Lymph-Vascularized Human Stromal Tissues for Preclinical Assays of Cancer Metastasis. Biomaterials, 179, 144-155.