研究課題/領域番号 |
20K12064
|
研究種目 |
基盤研究(C)
|
配分区分 | 基金 |
応募区分 | 一般 |
審査区分 |
小区分62010:生命、健康および医療情報学関連
|
研究機関 | 九州大学 |
研究代表者 |
嶋吉 隆夫 九州大学, 情報基盤研究開発センター, 学術研究員 (60373510)
|
研究分担者 |
竹田 有加里 福井大学, 学術研究院医学系部門, 助教 (20582159)
|
研究期間 (年度) |
2020-04-01 – 2024-03-31
|
研究課題ステータス |
交付 (2022年度)
|
配分額 *注記 |
4,290千円 (直接経費: 3,300千円、間接経費: 990千円)
2022年度: 910千円 (直接経費: 700千円、間接経費: 210千円)
2021年度: 780千円 (直接経費: 600千円、間接経費: 180千円)
2020年度: 2,600千円 (直接経費: 2,000千円、間接経費: 600千円)
|
キーワード | 計算生理学 / データ科学 / イオンチャネル数理モデル / ベイズ推論 / MCMC法 / レプリカ交換法 / 生物物理学的細胞数理モデル / 個体差分布 / パラメータ解析 |
研究開始時の研究の概要 |
細胞を構成するタンパク質などに対する生理学実験データには個体差によるばらつきが生じる.生理学実験結果に基づき細胞の機能を数理的にモデル化してシミュレーションにより解析する研究では従来,個体差に起因するモデル内定数のばらつきを考慮できていなかった.本研究課題は,データサイエンス分野の方法論を応用することで,個体差を含んだ生理学実験データからモデル内定数の値の個体差による分布を推定する手法を開発し,モデル内定数値分布を考慮してモデル開発とシミュレーションを行い,細胞ごとの表現型の差異を解析するための手法を開発することを目的とする.
|
研究実績の概要 |
細胞を構成するタンパク質などに対する生理学実験データには個体差によるばらつきが生じる.分子レベルの生物物理学的現象についての生理学実験データから同定した要素モデルを集約して構築した包括的細胞モデルでは,従来,実験データに含まれる個体差によるばらつきを考慮せずにパラメータ値を点推定しており,細胞レベルのマクロな振る舞いにおける個体差による差異は再現できなかった. 前年度では,あるモデルパラメータが与えられたときにヒストグラムデータが得られる確率を算出する方法を導出したが,その算出には線形代数計算が必要であり,一般的なMCMC法の実装では確率モデルは陽的に記述されることを前提としていることから,線形代数計算の実装方法について検討・試行を進めていた. 令和3年度は,まず,一般的なMCMC法の実装であるStanに対して,MCMC計算の内部で線形代数計算を実行するC++言語による実装を完成させた.次いで,その実装を用いて,既知パラメータから生成した模擬データを対象に推定を試みたが,対象とするパラメータ空間は評価値(対数確率)の凹凸が非常に多く,一般的なMCMC法ではパラメータが推定できないことが分かった.そこで,Stanを用いてレプリカ交換法を実行する実装を行い,さらにレプリカ交換法を一部変形してテストデータに適用することで,真値の推定に成功した.しかし,実際の実験データに対して適用したところ,推定結果が安定しないことが分かった.さらなる調査検討を行っているが,対象としている実験データが既に縮約されたヒストグラムデータであることから,パラメータ推定に対しては不十分である可能性も考えられる.
|
現在までの達成度 (区分) |
現在までの達成度 (区分)
3: やや遅れている
理由
今年度は精力的に研究を進め,実装面では大きな前進があったが,当初予想より対象データの推定難易度が高いことなどから,令和3年度に新型コロナ対策への対応に起因して生じた遅れを十分には挽回できていない.
|
今後の研究の推進方策 |
まず,対象とする実験データについて,パラメータ空間の性質を可視化し,今後の対応を検討する.また,パラメータの値に分布を持たせ,分布に従ってランダムに選択したパラメータセットから実験データを模擬するテストデータを生成し,そのデータに対してパラメータ推定を試みる. また,本年度が最終年度であることから,研究成果のとりまとめを行う.
|