列決定逆問題におけるパラメタ推定法と薬物動態モデルへの適用
薬 の 効果や 副作用 は 体内 の 薬物濃 度 に 影響 され るの で ,薬物 が 体 内 に お い て どの よ うに 吸 収 ,拡 散そ して 排 出 され る か と い う薬物動 態 の 理 解 は薬 品 開発 に不 可 欠で あ る.体 内 の 臓 器 を微 分方程 式で 表 し ,人体 モ デ ル を構 成す る Whole Body Physiologically Based Pharmaco-Kinetics model (PBPK model )は 製 薬 の 現 場 で は 広 く認 知 され た シ ミ ュ レ ーシ ョ ン技 法 で あ る が、実験から得られるデータに比べてモデルとなる連立方程式のパラメタ数のほうが多く、解空間を求めることが困難な典型的な列決定逆問題であった。この逆問題を効率よく解くアルゴリズムを小長谷研究室と2011年当時にUppsala大学から国立情報学研究所速水謙研究室にMOUインターシップ生として来日していた青木康憲氏(現AstraZeneca)と開発し、Cluster Newton Method (CNM)と名付けた。そして、CNMをSlatter ら が 行 っ た 抗がん剤lrinotecanおよびその代謝物の肝排泄量と腎排泄量の比較データ( mass −balance−study)を解析するために作成したPBPKモデル(有熊 威、小長谷研究室)に適用した。これにより、胆管腫瘍患者の肝排泄量と腎排泄量の比が他の腫瘍患者と比べて大きく異なる理由をPBPKモデルを用いて推察することに成功した。CNMはSIAM Journal on Scientific Computingに採録されると共に、日本応用数理学会より2015年度ベストオーサー賞を受賞した。
CNMは、理研GSCおよび東京工業大学計算工学専攻小長谷研究室で推進していた薬物相互作用オントロジーの研究(吉川澄美, 理研GSC, 2004)および粒子シミュレーションによる分子間相互作用の研究(我妻竜三, 理研GSC,2006) が土台となっている。理研GSCおよび東工大小長谷研究室(計算工学科)における薬物間相互作用の解析を目指したオントロジーと数理モデルを融合する研究の中から、薬物相互作用オントロジーから連立微分方程式モデルを自動生成する研究(有熊威, 2007年度修論)が生まれた。そこに数理モデリングにおけるパラメタ-パラメタ相互作用の研究(我妻竜三,理研GSC, 2007)が加わり、薬物動態モデルの解空間を逆問題として解くことの重要性が認識された。そして、逆問題を解くためのアルゴリズムとして、機械学習による仮想患者集団収束法の研究(渡辺崇, 2008年度修論)ならびにKernel PLSの零空間サンプリングによる仮想患者集団収束法の研究(尾崎 賢人, 2008年度卒論)が始まった。CNMはこの仮想集団収束法と零空間サンプリングを兼ね備えたアルゴリズムであり、仮想患者集団をニュートン法で集団ごと零空間にマッピングすることで、ロバストな線形近似と効率的な零空間サンプリングを実現した。吉田 健太(東大杉山研究室から小長谷研究室に研究指導を委託された博士課程学生(当時))は青木らが開発したCNMを有熊が作成した数理モデルに適用し、10種程度の少ない観測実験データを50種以上の理論パラメタで補完した薬物動態モデルで胆管腫瘍患者の代謝パラメタが他の患者と大きく異なることを示した(2013, BMC Bioinformatics)。また、2012年のthe 25th JSAAE Meeting においてBest Presentation Awardを受賞し、2014年にCNMを薬物間相互作用の解析に応用し、「クラスターニュートン法を用いた薬物間相互作用の生理学的薬物速度論モデル解析に関する研究」で学位を取得した。
近年、CNMは、薬物動態のみならず、システム生物学の分野でも逆問題を解くために活用され始めている。また、後継であるCluster-Gauss Newton Methodが開発され、薬物動態研究で活用されはじめている。
薬物代謝モデルにおけるパラメタ-パラメタ相互作用の研究
オントロジーとシミュレーションを用いたハイブリッド型薬物間相互作用予測システムの研究
薬物間相互作用オントロジーからの数理モデルの自動生成の研究
仮想集団収束法の研究
Kernel PLSを用いた零空間サンプリングの研究
Cluster Newton Methodの研究
CNMを用いた胆管腫瘍患者PBPKモデルの解空間推定の研究
CNMを用いた薬物間相互作用解析の研究
Azuma R., Umetsu R., Ohki S., Konishi F., Yoshikawa S., Konagaya A., and Matsumura K.
Discovering dynamic characteristics of biochemical pathways using geometric patterns among parameter-parameter dependencies in differential equations.
New Generation Computing Vol.25, No.4, pp.425-441 (2007)
https://doi.org/10.1007/s00354-007-0027-3
小長谷明彦、有熊威、吉川澄美、我妻竜三、梅津亮、小西史一, オントロジーとシミュレーションを用いたハイブリッド型薬物間相互作用予測システムの設計思想, シミュレーション, 日本シミュレーション学会, Vol.26, No.4, pp.31-41 (2007).(2008年論文賞受賞)
Takeshi Arikuma, Sumi Yoshikawa, Ryuzo Azuma, Kentaro Watanabe, Kazumi Matsumura, Akihiko Konagaya: Drug interaction prediction using ontology-driven hypothetical assertion framework for pathway generation followed by numerical simulation.
BMC Bioinformatics 9(S-6) (2008).
https://link.springer.com/article/10.1186/1471-2105-9-S6-S11
サポートベクトルマシンを用いた仮想集団収束方法に関する研究
修士論文 渡辺 崇 (2008年度)
・赤点:生成仮想患者集団
・青枠:胆管腫瘍患者の肝腎排泄量実測値の許容範囲
kernel PLS上の零空間サンプリングを用いた仮想集団の収束方法に関する研究
卒業論文 尾崎 賢人 (2008年度)
・kernel PLSによる微分方程式の入力パラメタ空間と出力パラメタ空間の線形近似
・線形近似モデルの零空間からの解候補のサンプリング
・左図は、SN-38(イリノテカン)、SN-38G(グルクロン酸抱合体)の肝排泄量について、文献値、有熊推定値、渡辺推定値と比較したグラフ
Aoki Yasunori, Hayami Ken, SterckHans De, Konagaya Akihiko.
Cluster Newton Method for Sampling Multiple Solutions of an IUnderdetermined nverse Problem: Parameter Identification for Pharmacokinetics.
NII Report
https://www.nii.ac.jp/TechReports/public_html/11-002E.pdf
SIAM Journal on Scientific Computing, 36 (1) pp.B14-B44 (2014).
https://doi.org/10.1137/120885462
応用数理, 24(4) pp.151-159 (2014) に発刊した論文「劣決定逆問題に対するCluster Newton法とその薬物動態モデルへの応用」が日本応用数理学会よりベストオーサー賞受賞(2015年)
Kenta Yoshida, Kazuya Maeda, Hiroyuki Kusuhara, Akihiko Konagaya:
Estimation of feasible solution space using Cluster Newton Method: application to pharmacokinetic analysis of irinotecan with physiologically-based pharmacokinetic models. BMC Systems Biology 7(S-3): S3 (2013)
https://doi.org/10.1186/1752-0509-7-S3-S3
Kenta Yoshida, Kazuya Maeda, Akihiko Konagaya, Hiroyuki Kusuhara: “Accurate Estimation of In Vivo Inhibition Constants of Inhibitors and Fraction Metabolized of Substrates with Physiologically Based Pharmacokinetic Drug–Drug Interaction Models Incorporating Parent Drugs and Metabolites of Substrates with Cluster Newton Method“, DRUG METABOLISM AND DISPOSITION, 46, pp.1805-1816 (2018)
https://dmd.aspetjournals.org/article/S0090-9556(24)05590-9/abstract