JJSLSM Vol.4 No.4 (22) 3 原著 ORIGINAL ARTICLE 熱損傷シミュレーションに基づく皮膚良性色素疾患治療用ナノ秒パルスレーザー装置の計算機臨床試験方法 西村隆宏, 下条裕, 間久直, 粟津邦男,2,3 大阪大学大学院工学研究科 2 大阪大学大学院生命機能研究科 3 大阪大学国際医工情報センター (29 年 月 22 日受付,29 年 7 月 9 日掲載決定,29 年 月 26 日 J-STAGE 早期公開 ) A Method of Computational Clinical Trial of a Nanosecond Pulsed Laser Skin Treatment Device by Numerical Simulation of Photothermal Damage Takahiro Nishimura, Yu Shimojo, Hisanao Hazama, Kunio Awazu,2,3 Graduate School of Engineering, Osaka University 2 Graduate School of Frontier Biosciences, Osaka University 3 Global Center for Medical Engineering and Informatics, Osaka University (Received January 22, 29, Accepted July 9, 29, Advance Publication released online October 26, 29) 要旨本研究では, 新規レーザー治療装置に対する迅速かつ低コストな評価手法として, レーザー治療における計算機臨床試験を提案する. レーザー光照射により生じる生体組織内の光熱伝搬とそれに伴う熱損傷過程をモデル化し, レーザー治療における安全性を数値シミュレーションにより評価する. 皮膚良性色素疾患治療用ナノ秒パルスレーザー装置の承認申請機器を対象とし, パルス幅の異なる既承認機器と比較して熱損傷の観点から安全性が同等であることが, 計算機臨床試験により評価可能であることを示した. キーワード : 計算機臨床試験, レギュラトリーサイエンス, 医療機器審査, ナノ秒レーザー, 熱損傷 Abstract This paper proposes a computational clinical trial of laser treatment for rapid and low cost evaluation of novel medical laser devices. The optical propagation and photothermal damage process after pulse laser irradiation are modeled and the safety of the laser treatment with the evaluation target devices is evaluated by numerical simulation. The computational clinical trial result shows the safety of a novel nanosecond laser device for skin treatment is equivalent to the device that has been approved for clinical use. Key words: computational clinical trial, regulatory science, medical device review, nanosecond laser, photothermal damage 565-87 大阪府吹田市山田丘 2- TEL: 6-6879-852, FAX: 6-6879-852 (2-, Yamadaoka, Suita, Osaka 565-87, Japan) Corresponding author: nishimura-t@see.eng.osaka-u.ac.jp( 西村隆宏 )
32 日レ医誌 (JJSLSM) 第 4 巻第 4 号 (22). 背景超短パルスレーザーをはじめとしたレーザー技術の進展は, レーザー治療の高度化や適用範囲の拡大につながる,2). ミリ秒からピコ秒オーダーのパルス幅を持つレーザー装置を用いた治療は, レーザー光が組織や病変に吸収され,2 組織や病変の温度が上昇し,3 組織や病変が凝固, 溶融, 膨張, 蒸散することにより, 治療効果が得られる 3). 各作用の物理モデルやそのパラメータは, レーザー光の波長, 照射時間幅, 照射フルエンス ( あるいはピーク強度 ) により変わる. そのため, 特性の異なる新規レーザー装置の積極的な臨床利用により, 新たな治療法や, 従来と比較してより安全でより効果的なレーザー治療の展開が期待できる. レーザー治療では外部から光エネルギーを注入して治療効果を得るため, リスクが全くないことはありえない. そのため, レーザー装置を臨床利用するためには, その装置のリスク評価が必要となる. 医療機器におけるリスクの評価には, 装置自体の品質と生体への安全性の評価が求められる. レーザー治療装置としての品質に関しては,JIS 規格や IEC 規格など標準化が進められており, 安全規格や基準が定められている 4). 一方, 生体に対する安全性評価には, 主に, 臨床試験や動物実験のデータをエビデンスとした統計的評価手法がとられている. 臨床試験にはコストと時間が必要になるため, レーザー装置の安全性を十分に評価するための障壁となる. これは, 日進月歩のレーザー装置開発と臨床応用とに時間差を生じさせるだけでなく, レーザー医療の進歩を停滞させる要因となり得る. また, 未承認レーザー装置を用いた施術は, 医療過誤の要因の一つともいえる. 新規レーザー装置に対して, 低コストで迅速に十分な評価を行うためのレギュラトリーサイエンスが重要となる. レーザー装置を医療現場へ迅速に展開する取り組みは, 日本の医療機器承認審査を担う独立行政法人医薬品医療機器総合機構においても進められている. 例えば平成 28 年に, レーザー医療機器の承認申請では, 使用目的と作用原理の既承認機器が存在し, 性能及び使用方法に既承認機器との同等性が認められる場合, 臨床試験の成績の提出が不要と定められた 5). 臨床試験の成績提出が不要とされるパスも用意されたが, 性能及び使用方法に既承認機器との同等性の証明が必要である. 現状では, 追加の評価試験として動物やファントムを用いた ex vivo 実験が必要なままである. 低コスト 迅速な評価手法の確立は, 医薬品や医療機 器開発における共通の課題である. その中で, 前臨床試験や臨床試験に加えて第三の試験として, モデリングとシミュレーションに基づく Computational clinical trials( 計算機臨床試験 ) の重要性が高まっている 6). これまで物理空間で実施されていた前臨床試験 臨床試験を計算機上に再現することにより, 迅速 低コストな医療機器評価が可能となる. また, これまでの評価試験や実験では, 例えば, 臨床試験データにおける患者の個人差や手技者によるばらつきや, 動物実験データのヒトへの外挿性の科学的根拠が課題と考えられるが, シミュレーションに基づく計算機臨床試験ではそれらを排除して, 明確な前提条件のもと生体への作用を定量的に扱うことができる. レーザー治療においては, 種々の物理パラメーターにおける光 生体組織相互作用のシミュレーションがこれまでに多く示されている 7-). レーザー照射による蒸散深さシミュレーションでは, 非臨床試験結果との同等性が確認されており, シミュレーションに基づく評価の有効性がすでに実証されている ). これは, 光 生体組織相互作用における熱的影響のモデル化により, 高精度に数値シミュレーションできることを意味する. そのため, レーザー医療機器の承認申請に必要となる既承認機器との同等性の証明においても,ex vivo 実験に替わる, モデリングとシミュレーションに基づく計算機臨床試験の適用が期待できる. レーザー医療機器に対する計算機臨床試験により, 迅速かつ低コストに定量的安全評価が実現でき, 新規光学機器の医療 医用応用の促進につながる. 2. 目的本研究の目的は, 新規レーザー医用装置の計算機臨床試験の一例を示し, 新たなレーザー治療装置おける評価のための方法論を示すことである. 一例として, 本研究では, 皮膚良性色素疾患治療用ナノ秒パルスレーザー装置を用いた治療に対する計算機臨床試験を行う. ナノ秒パルスレーザー照射における熱損傷過程をモデル化し, 熱損傷体積を数値シミュレーションにより求める. 正常皮膚組織へパルス光照射した際の温度上昇と熱損傷体積を安全性の指標として, 既承認機器に対する承認申請機器の安全性を評価する. 3. 対象と方法 3. 対象本研究では, 既承認レーザー装置 ( 装置 A, 装置 B) の 2 機種を用いた. それぞれの仕様を Table に示す. い Table Specifications of Q-switched alexandrite laser devices that are numerically evaluated in this computational clinical trial Device A Device B Wavelength [nm] 755 ± 2 755 ± 5 Repetition rate [Hz], 2, 3, 4, 5, 2, 5 Spot diameter [mm] 2, 3, 4 2, 3, 4 Fluence [J/cm 2 ] 8, 5, 3. 3.5 8, 5, 3 5 Pulse width [ns] 54 5 Output stability (3σ) [%] ±5.5 ±5.5 Output accuracy [%] ±2 ±5.5
JJSLSM Vol.4 No.4 (22) 33 ずれもフラッシュランプ励起式 Q スイッチアレキサンドライトレーザーであり, 皮膚良性色素疾患治療用として承認されている. アレキサンドライトレーザーから出力される波長 755 nm は, 紫外域, 他の可視域および中赤外域の波長の光と比較して, ヒト皮膚組織への光侵達深さが比較的大きい. 装置 A が先に承認された経緯から, 本研究では装置 A を既承認機器, 装置 B を承認申請機器と想定して, 装置 A に対する装置 B の安全性を計算機臨床試験にて評価した. 両機器にてレーザー光を皮膚組織へ照射した際の, 組織内の温度分布の時間変化とそれに伴う熱損傷体積を比較項目とした. 3.2 方法ナノ秒パルスレーザーを照射した際に生じる皮膚組織への熱損傷作用のシミュレーションモデルを構築した. Optical properties (μ a, μ s, g) Thermal properties (k, ρ, c) Skin structure model Light propagation model S Thermal diffusion model T Thermal damage model Tissue damage (τ) Boundary condition Damage parameters (A, E a ) Fig. Outline of calculating photothermal damage by a laser pulse irradiation. Light propagation in skin structure model is calculated with optical property parameters. The absorbed energy distribution is used for the initial value of a heat source to calculate thermal diffusion. From the calculated spatiotemporal temperature change, spatial distribution of thermal damage ratio is obtained numerically. 熱損傷計算の概要を Fig. に示す. まず, 皮膚組織モデルを構築し, モンテカルロ法による光伝搬計算によって, 皮膚組織に吸収される光エネルギー分布を求めた ). 皮膚組織中の光伝搬速度と比較すると, 熱拡散速度は無視できるほど低い. この近似の下, 全パルスエネルギーがパルス幅の間に均一に入射すると仮定し, 組織内の温度上昇とその後の熱伝搬による温度変化を計算した. その後, 温度分布を元に, アレニウスモデルにより, 凝固や細胞死領域を算出した. 以下ではその詳細を示す. 3.2. 皮膚組織モデル本計算機臨床試験では,Fig.2 に示す三次元皮膚組織モデルを仮定した 2). 組織の面積を 5 5 mm 2 とし, 厚さを 4 mm とした. 皮膚組織は平面かつ層状の構造とし, 表皮, 真皮, 皮下組織の 3 層モデルを仮定した. より詳細な多層皮膚モデル 3) などの多層モデルとシミュレーション結果を比較したところ差異がなかったため, 本研究では 3 層モデルを採用した.Fig.2 に示すように, 組織の深さ方向を z 軸とし, 表皮上面を z = とした.z 軸に垂直に xy 面をとり, 三次元皮膚組織モデルの xy 面中心を (x, y) = (, ) とした. 各層は文献 2) を参考に, 表皮を. mm, 真皮を. mm, 皮下脂肪を 2.8 mm と設定した. 真皮には毛細血管, 真皮上部血管, 真皮下部血管が含まれる血管叢を仮定した. 各血管の構造は文献 4) より, Table 2 に示す値を用いた. 血液の体積分率と血管の体積分率が等しいとし, 血管の体積分率から組織中に含まれる血管本数を決定した.Tables 3,4 に示す光学特性値, 熱伝導パラメータを, 皮膚組織モデルの各組織の対応する位置へ割り当てた三次元モデルを構築した 5-23). 皮膚組織モデルの作成には, 数値計算ソフトウェアの Matlab (MathWorks Inc.) を使用した. 3.2.2 光伝搬シミュレーション皮膚組織内の光伝搬計算を, モンテカルロ法にて計算した ). レーザー光の波長は, フラッシュランプ励起式 Q スイッチアレキサンドライトレーザーの 755 nm とした. レーザースポットサイズをメーカー公称値である 2, 3,4 mm とし, それぞれの場合において, 皮膚組織表面に対して均一な強度で垂直に入射した際のエネルギー密 2. 4. -2.5 Epidermis Capillary vessel Upper dermis vessel Dermis Deep dermis vessel Subcutaneous fat y [mm] 2.5 2.5-2.5 Fig.2 Three dimensional skin model for the computational clinical trial.
34 日レ医誌 (JJSLSM) 第 4 巻第 4 号 (22) Table 2 Parameters for vessels 4) Center depth [μm] Diameter [μm] Thickness [μm] Volume fraction [%] Capillary vessel 75 5 4 Upper blood net plexus 29 5 8 3 Deep blood net plexus 6 8 8 Table 3 Optical properties of skin model. μ a, μ s and g are absorption coefficient, scattering coefficient, and anisotropy factor 5-8) μ a [cm ] μ s [cm ] g [-] Epidermis 22.6 394.9 Dermis.25 238.9 Subcutaneous.9 8.9 Blood 3. 32.9 Table 4 Thermal properties of skin tissue 9-23) Density [g/cm 3 ] Specific heat [J/(g C)] Thermal conductivity [W/(cm C)] Epidermis.9 3.6 2.9 3 Dermis.2 3.22 3. 3 Subcutaneous fat.97 2.3 2.5 3 Blood. 3.84 4.92 3 度を計算した. 空間グリッドサイズを μm とし, 入射光子数は.2 7 とした.Oregon Medical Laser Center が公開している, 組織内の光伝搬シミュレーションコード ) を使用した. 3.2.3 温度分布シミュレーション光伝搬シミュレーションにより得た算出した吸収エネルギーの空間分布を熱源エネルギーの空間分布として, 各位置における温度上昇, 熱拡散過程を計算した. 三次元熱伝導方程式は, 皮膚組織の温度 T, 熱源のエネルギー密度 S を用いて, ρc p T t = λ 2 T + S, () と表せる ). ここで,ρ,C p,λ は皮膚組織の密度, 比熱, 熱伝導率を表す. 皮膚組織, 空気層の初期温度は 32 C,24 C とした. 皮膚モデル外への熱伝達はないものと仮定した. 空間グリッドサイズを μm とし, 有限要素法により算出した. 時間ステップ間隔は, パルス光照射中はパルス幅を 分割した時間とし, パルス照射後は ns とした. 温度分布シミュレーションコードは C 言語を用いて実装した. 3.2.4 熱損傷シミュレーション熱損傷はアレニウスモデルに従ってモデル化し 24), 温度分布シミュレーション結果を元に熱損傷を計算した. 組織の熱損傷指標 はアレニウスの式から, = ln C C t t = A exp ΔE dt (2) RT と表すことができる.C はある時間における熱変性を受けていないタンパク分子濃度である.A,ΔE,R,T はそれぞれ, 頻度因子 (/s), 活性化エネルギー (J/mol), 気体定数 (J/(mol k)) 絶対温度 (K) である. が大きいほど熱変性タンパク質の割合が高くなり, 生体組織の熱損傷割合が高くなる. となる場合, すなわち熱変性を受けていないタンパク質の割合 C(t)/C() が.37 より低くなると, 不可逆的損傷となる 25). 本研究においても, この不可逆的損傷を有意な損傷として, となる領域を熱損傷領域として計算により求めた. このとき, 皮膚組織の A,ΔE は実験的に求められた, A = 8.82 94 T 53 C.3 3 T > 53 C ΔE = 6.3 5 T 53 C 2.4 5 T > 53 C を用いた 26). (3) (4) 4. 結果装置 A, 装置 B ともに設定可能なスポットサイズが同じであったため, 各スポットサイズにおいて両者の皮膚
JJSLSM Vol.4 No.4 (22) 35 (a) (i) (ii) (iii) 2.5... 2 2. 2. 2..5 3. 3. 3.. (b) 4. 4. -2. -.. 2. -2. -. 4.. 2. (i) (ii) (iii) -2. -.. 2..5 S 2.5... 2 2. 2. 2..5 3. 3. 3.. 4. -2. -.. 2. 4. -2. -.. 2. 4. -2. -.. 2..5 Fig.3 Spatial distributions of (a) light fluence ϕ and (b) absorbed energy S on zx plane (y = ) after a laser pulse irradiation with device A. The spot sizes are (i) 2, (ii) 3, and (iii) 4 mm. 損傷を比較した. パルス光エネルギーは, 各スポットサイズにおいて公称値の最大となるように設定した. この設定のもと, 装置の出力しうる照射パワー密度となる条件として, 出力安定性を考慮して公称値から 3σ 加算した照射フルエンスかつ装置の出力しうる最小のパルス幅となる条件を過酷条件と定義した. 一方, 公称値の照射フルエンスかつ装置の出力しうる最大のパルス幅の条件を穏和条件と定義した. レーザー光は, 組織モデルの組織表面の中心に垂直方向に照射した. 装置 A により, 各スポットサイズにて照射した際の組織内光エネルギー分布と吸収エネルギー分布を Fig.3 に示す. 装置 B においても, 波長, スポットサイズは同じであり, パルスのエネルギーに大きな違いがないため, 同様の組織内光エネルギー分布と吸収エネルギー分布が得られ, 光侵達深さはほぼ等しい結果が得られた. この場合の光侵逹深さは,.95 mm( スポットサイズ :2 mm),. mm( スポットサイズ :3 mm),.6 mm( スポットサイズ :4 mm) と求まった. 一方, 装置 B においても同様に組織内光エネルギー分布と吸収エネルギー分布を計算した結果, 光侵逹深さは.95 mm( スポットサイズ :2 mm),. mm ( スポットサイズ :3 mm),.6 mm( スポットサイズ : 4 mm) であった. 装置 A, 装置 B においても, 波長, スポットサイズが同じであるため, 光侵逹深さは同じ結果となった. スポットサイズ間の吸収エネルギー分布を比較すると, スポットサイズが 4 mm の場合, ほぼ表皮層で光吸収しており深部への熱影響は小さいことがわかる. 装置 A, 装置 B いずれにおいても, 最も照射エネルギー密度が高くなるスポットサイズ 2 mm の場合の結果を挙げる.Fig.4 に, 装置 A と装置 B のレーザーパルスを照射した際の, 直後 ( s),5,,5,2 ms 後の組織中心における深さ方向の温度分布を示す. 装置 A, 装置 B ともにレーザーパルス照射直後 s では, 波長 (a) Temperature [ C] 2 5 5.5..5 2. s 5 ms ms 5 ms 2 ms Temperature [ C].5..5 2. 775 nm の吸収係数が大きい表皮, 血管領域に対する z =.,.8,.29,.6 mm 付近で温度上昇のピークが観察できる. 表皮組織における最高温度は, 装置 A と装置 B ともに, レーザー光照射直後の表皮層で計測された 26 C であった. その後, 周辺組織への熱拡散により, 表皮, 血管領域の温度は低下した.2 ms 後には, 表皮層の最高温度は 9 C まで低下した. 装置 A と装置 B ともに温度上昇に差はなかった. 装置 A と装置 B に各スポットサイズにおける, 最高温度を Fig.5 にまとめる. いずれも照射直後の表皮組織で最大となり, 最高温度の差は, スポットサイズ 4 mm における,6.9 C であった. レーザーパルス照射 2 ms 後の熱損傷指標の y = となる,zx 平面での空間分布を Fig.6 に示す. 装置 A の穏和条件と装置 B の過酷条件では, いずれのスポットサイズにおいても, 熱損傷指標の分布は深部までわずかに広がっているものの, 両者に大きな差はみられない. 装置 A, 装置 B の, 熱損傷領域の計算結果を Fig.7 にまとめ (b) 2 5 5 s 5 ms ms 5 ms 2 ms Fig.4 Temperature profiles of skin tissue on z axis (x =, y = ) (immediately after irradiation), 5,, 5, and 2 ms after a single pulse irradiation with (a) device A and (b) device B.
36 日レ医誌 (JJSLSM) 第 4 巻第 4 号 (22) る. 装置 A の穏和条件における損傷体積は, スポットサイズ 2,3,4 mm では,.88,.38, mm 3 であった. 装 Maximum temperature ( C) 25 2 5 5 Device A Nominal maximum value Nominal maximum value + 3σ Device B 2 3 4 2 3 4 Spot size (mm) Fig.5 Maximum temperature in skin tissues after one single pulse irradiation. (a) (b).5. -2.5.5 (i) (ii) -2. -.5 -. -.5.5..5 2. 2.5. -2.5-2. -.5 -. -.5.5..5 2. 2.5 (iii).5. -2.5-2. -.5 -. -.5.5..5 2. 2.5.5. -2.5.5 (i) (ii) -2. -.5 -. -.5.5..5 2. 2.5. -2.5-2. -.5 -. -.5.5..5 2. 2.5 (iii).5. -2.5-2. -.5 -. -.5.5..5 2. 2.5 Fig.6 Spatial distributions of thermal damage index on zx plane (y = ) by irradiation with (i) 2, (ii) 3, and (iii) 4 mm sized laser spot of (a) device A and (b) device B. 置 A の過酷条件における損傷体積は, スポットサイズ 2, 3,4 mm では,.97,.53,.99 mm 3 であった. 装置 B の穏和条件における損傷体積は, スポットサイズ 2,3, 4 mm では,.88,.38, mm 3 であった. 装置 B の過酷条件における損傷体積は, スポットサイズ 2,3,4 mm では,.97,.53, mm 3 であった. 各条件においても, 熱損傷部位体積は, 装置 A を用いてスポットサイズ 3 mm とした場合に最大となった. 穏和条件, 過酷条件ともに, 既承認機器である装置 A を用いた場合に生じ得る熱損傷体積に対して, 装置 B を用いた場合に生じる熱損傷体積は同等または小さくなった. 計算機臨床試験から, 正常組織に生じる熱損傷の観点から装置 B は装置 A と比較して同等もしくは高い安全性を有するといえる. 5. 考察 Fig.5 と Fig.7 に示したように, 計算機臨床試験の結果は, パルス幅による温度上昇, 熱損傷に大きな差はみられなかった. 今回比較対象とした装置 A, 装置 B では波長と光照射フルエンスが同じであり, 比較したパルス幅が異なった. 各層の熱拡散係数 α = λ/ρc p は, 表皮, 真皮, 皮下組織中で,4.88 4,8.32 4,9.9 4 [cm 2 /s] となる. 皮膚組織モデルにおける熱閉じ込め条件を満たすパルス幅 τ は,τ = d 2 /4α, と表せる.d は各層の厚さであり,τ は, 表皮, 真皮, 皮下組織で, それぞれ,.5,3.6,2.3[s] と求まる. 装置 A, 装置 B の両条件とも熱閉じ込め条件を満たしていたことがあげられる. 今回のシミュレーションでは, 熱閉じ込め条件近似を用いず熱伝搬を計算していたが, 温度上昇, 熱損傷が両条件ほぼ等しい結果となったのはこのためである. ナノ秒パルス領域においては組織の体積膨張による熱弾性波も生じ得るが, 今回の計算機臨床試験では熱損傷が最大となる条件として, 組織が吸収した光エネルギーは全て熱変換され作用するものとした. この場合, 組織の光吸収後のモデリングに必要なパラメータは温度のみであり, 本研究でも利用したように三次元熱伝導方程式により精度よく計算機上で再現できる. また, 熱損傷の Damaged volume (mm ) 3 2.5.5 Device A Nominal maximum value Nominal maximum value + 3σ Device B 2 3 4 2 3 4 Spot size (mm) Fig.7 Thermal damages after one single pulse irradiation.
JJSLSM Vol.4 No.4 (22) 37 生じるしきい値などのパラメータなど, これまでに十分な知見が得られており, 熱損傷評価には数値シミュレーションを取りれることが標準となっている 24). 装置 A, 装置 B の両条件とも正常生体組織への作用として熱損傷を計算することにより, 安全性が高精度に評価できる. Fig.4 に示すように, 両条件においても 2 ms 後の温度は初期値である 32 C まで戻らず,9 C と温度上昇したままであった. 装置 A, 装置 B ともに最大繰り返し周波数は 5 Hz であるため, 同一箇所への連続照射した場合, 熱の蓄積による温度上昇が想定される. そのため, 同一箇所へのパルス光照射間隔は十分にとる, 施術中の皮膚組織を冷却するなどの対応が求められる. ただし, 本皮膚組織モデルでは光吸収の大きい血管の血流が考慮されていないため, 温度上昇に対しては過剰評価となっている. 血管系による冷却効果が大きな温度変化を与えることは多くの報告例があるが 24), 数値シミュレーションにより現象を再現するための物理パラメータは十分に蓄積されていない. 今後, 実験的な解析により, 適切な冷却機能等の設計が可能となる. 外科的な処置を使用目的とするレーザー機器の性能評価指針 5) において, 実際の手技を反映した使用方法により, 機器の設定可能なモード, 出力全体で行うことが求められている. これまでは, 動物やファントムを用いた試験が必要であったのに対して, 計算機臨床試験においては, 計算機上でパラメータを設定するのみの簡便な操作で試験を完結することができる. 本研究で利用したモンテカルロ法による光伝搬計算, 熱伝導方程式による温度計算はすでに広く利用される確立された手法である. また, アレニウスモデルによる熱損傷計算では, 生体組織における熱的作用を高精度にシミュレーション可能であることが示されている ). そのため, ヒト皮膚組織においても高精度にシミュレーションできていると考えられるが, 今後, 実測実験との比較などによる精度保証が課題となる. 本論文では, 安全性評価として正常組織に対してパルス光を照射した際の熱損傷評価にとどまっているが, 病変部位を含む皮膚組織モデルを導入することにより, 今後, 計算機臨床試験による有効性評価への適用が期待できる. 6. 結論本研究では, 計算機臨床試験により, 承認申請機器と想定したナノ秒パルスレーザーの安全性を評価した. 既承認機器における穏和条件と承認申請機器における過酷条件の熱損傷を計算機シミュレーションにより比較し, 両機器の安全性が同等であると結論づけた. 今回対象とした機器は実際に承認されており, 本計算機臨床試験の示す結果と同じである. 今後, 臨床研究データや物理パラメータの高精度な取得により, 計算機臨床試験の確度や信頼性を向上させることにより, 新たな医療機器評価手法として展開できる. 利益相反の開示開示すべき利益相反なし. 引用文献 ) 葛西健一郎. ピコ秒レーザーによる刺青と皮膚良性色素性病変の治療. 日レ医誌.27; 37: 44-446. 2) 大城貴史, 大城俊夫, 佐々木克己. ピコ秒レーザーによる刺青および良性色素性病変に対する治療.27; 37: 427-434. 3) 石原美弥. レーザー生体相互作用とその治療と診断への応用. 光学.22; 4: 548-555. 4) 橋新裕一. 光 レーザーの安全基準とその動向. 日レ医誌.22; 33: 27-33. 5) 厚生労働省医薬 生活衛生局医療機器審査管理課長. レーザ医療機器の承認申請の取扱いについて. 薬生機審発 629 第 4 号.26. 6) Morrison TM, Pathmanathan P, Adwan M, Margerrison E. Advancing Regulatory Science With Computational Modeling for Medical Devices at the FDA s Office of Science and Engineering Laboratories. Front Med. 28; 5: 24. 7) Dua R, Chakraborty S. A novel modeling and simulation technique of photo thermal interactions between lasers and living biological tissues undergoing multiple changes in phase. J Biophoton. 25; 35: 447-462. 8) Cain CP, Polhamus GD, Roach WP, Stolarski DJ, Schuster KJ, Stockton K, Rockwell BA, Chen B, Welch AJ. Porcine skin visible lesion thresholds for near-infrared lasers including modeling at two pulse durations and spot sizes. J Biomed Opt. 26 Jul-Aug; : 49. 9) Bogni S, Stumpp O, Reinert M, Frenz M. Thermal model for optimization of vascular laser tissue soldering. J Biophoton. 2; 3: 284-295. ) 髙田隼也, 本多典広, 間久直, 庵谷尚正, 粟津邦男. 前立腺肥大症に対するレーザー蒸散術における熱影響深さの光 熱伝搬シミュレーションを用いた解析. 日レ医誌.25; 36: 44-449. ) Jacques SL, Li T. Monte Carlo simulations of light transport in 3D heterogenous tissues (mcxyz.c).( 参照 28 年 2 月 25 日 ),http://omlc.org/software/mc/mcxyz/index.html. 2) Lee Y, Hwang K. Skin thickness of Korean adults. Surg Radiol Anat. 22; 24: 83-89. 3) Maeda T, Arakawa N, Takahashi M, Aizu Y. Monte Carlo simulation of spectral reflectance using a multilayered skin tissue model. Opt Rev. 2; 7: 223-229. 4) Braverman IM. The Cutaneous Microcirculation. J Investig Dermatol Symp Proc. 2; 5: 3-9. 5) Jacques SL. Optical properties of biological tissues: a review. Phys Med Biol. 23; 58: R37-6. 6) Simpson CR, Kohl M, Essenpreisk M, Cope M. Near-infrared optical properties of ex vivo human skin and subcutaneous tissues measured using the Monte Carlo inversion technique. Phys Med Biol. 998; 43: 2465-2478. 7) Meinke M, Müller G, Helfmann J, Friebel M. Empirical model functions to calculate hematocrit-dependent optical properties of human blood. Appl Opt. 27; 46: 742-753. 8) Jacques SL. Oregon Medical Laser Center News. Skin Optics. ( 参照 28 年 2 月 25 日 ),https://omlc.org/news/jan98/ skinoptics.html. 9) Duck FA. Physical Properties of Tissues. Academic press; 99. p. 43-7. 2) Cohen ML. Measurement of the thermal properties of human skin. A review. J Invest Dermatol. 977; 69: 333-338. 2) Balasubramaniam TA, Bowman HF. Thermal Conductivity
38 日レ医誌 (JJSLSM) 第 4 巻第 4 号 (22) and Thermal Diffusivity of Biomaterials: A Simultaneous Measurement Technique. J Biomech Eng. 997; 99: 48-54. 22) Henriques FC, Moritz AR. 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. 947; 23: 53-549. 23) Atzler E, Richter F. Über Wärmekapazität des arteriellen und venösen Blutes. Biochem Zeit. 92; 2: 3-32. 24) 石島歩, 小関義彦, 東隆, 山本恵介, 立石敬介, 鎮西清行, 佐久間一郎, 小林英津子. 手術機器が与える熱損傷の非臨床的評価手法に関する調査. 日本コンピュータ外科学会誌.27; 9: 63-73. 25) Marqa MF, Colin P, Nevoux P, Mordon SR, Betrouni N. Focal laser ablation of prostate cancer: numerical simulation of temperature and damage distribution. Biomed Eng Online. 2; : 45. 26) Pearce JA. Relationship between Arrhenius models of thermal damage and the CEM 43 thermal dose. Proc SPIE. 29: 78; 784. 著者紹介西村隆宏 (Takahiro Nishimura) 23 年 3 月, 大阪大学大学院情報科学研究科情報数理学専攻博士後期課程修了 ( 博士 ( 情報科学 )).23 年 4 月, 日本学術振興会特別研究員 ( 大阪大学大学院医学系研究科 ).24 年 7 月, 大阪大学大学院情報科学研究科特任助教.28 年 6 月, 大阪大学大学院工学研究科環境 エネルギー工学専攻助教. 専門は情報光学, 生体光学. 日本レーザー医学会, 日本光学会, 応用物理学会, 各会員. 粟津邦男 (Kunio Awazu) 984 年 3 月, 神戸大学大学院工学研究科修士課程修了.984 年 4 月, 住友電気工業株式会社研究開発本部研究員. 993 年 6 月, テキサス州立大学 MD アンダーソン癌研究所客員研究員.995 年 9 月, 自由電子レーザ研究所 ( 通産省系プロジェクト研究所 ) 主任研究員,999 年 月, 超分子工学研究部部長.996 年 4 月, 博士 ( 工学 ) 神戸大学,997 年 3 月, 博士 ( 医学 ) 順天堂大学.2 年 5 月, 大阪大学大学院工学研究科教授となり現在に至る. 27 年 4 月, 大阪大学国際医工情報センター教授兼任. 2 年 月, 大阪大学大学院生命機能研究科教授兼任. 九州大学大学院客員教授, 光産業創成大学院大学客員教授, 千歳科学技術大学客員教授. 専門は, 医用工学, 生体組織光学, レーザー医療機器のレギュラトリーサイエンス. 日本レーザー医学会 ( 副理事長 編集委員長 医用レーザー機器開発推進委員長 ), 日本レーザー歯学会 ( 理事 ), 日本光線力学学会 ( 幹事 ), 日本レーザー治療学会 ( 理事 ), 日本レーザー獣医学研究会 ( 副会長 ) レーザー学会 ( 諮問員 ), 日本生体医工学会,SPIE 等. 下条裕 (Yu Shimojo) 25 年 4 月大阪大学工学部環境 エネルギー工学専攻入学. 日本レーザー医学会, レーザー学会, 日本光学会, 各会員. 間久直 (Hisanao Hazama) 2 年 3 月, 東京理科大学大学院理工学研究科電気工学専攻博士後期課程修了 ( 博士 ( 工学 )).2 年 4 月, 川崎重工業株式会社技術研究所.26 年 6 月, 大阪大学大学院工学研究科環境 エネルギー工学専攻, 特任研究員. 29 年 6 月, 同研究科附属高度人材育成センター助教. 24 年 4 月, 同研究科環境 エネルギー工学専攻講師. 25 年 4 月, 同専攻准教授. 日本レーザー医学会 ( 評議員, 安全教育委員会委員 ), 日本レーザー歯学会 ( 編集委員会査読委員 ), 日本レーザー治療学会 ( 評議員 ), 日本光線力学学会, 応用物理学会, 電気学会, レーザー学会, 日本質量分析学会,The American Society for Mass Spectrometry, 各会員.