|
最初に、前年度までにNMR法を用いて構築した、b2アドレナリン受容体の不活性型・部分活性型・活性型の構造平衡のモデルに対して、先行論文で報告されているGタンパク質の活性化および不活性化の数理モデル、GRKによるリン酸化とアレスチンの内在化の数理モデル、cAMPの生成およびPKAとPDEによるフィードバック阻害の数理モデルを組み合わせて、GPCRシグナルの数理モデルを構築した。
次に、前年度までに開発した交換モンテカルロ法のプログラムを使って、b2AR発現細胞に作動薬を添加した時のcAMP濃度の経時変化をシグナルの様々な箇所に阻害がかかった条件で観測した文献値が再現できるように、構築した数理モデルに含まれるパラメータを網羅的に探索して、ベイズ推定により最適値と誤差を算出することを試みた。その結果、構築した数理モデルに含まれるパラメータの数が49個と多いため、ベイズ推定を行うのに十分な量の計算を行うことが困難であった。そこで、交換モンテカルロ法のプログラムの改良を試みたところ、並列計算の手法を導入することにより、計算速度を5倍以上上昇させることに成功した。高速化したプログラムを使って数理モデルに含まれるパラメータを網羅的に探索した結果、各パラメータの最適値および誤差を得ることに成功した。得られたパラメータを使って、各種阻害がかかった条件におけるcAMP濃度の経時変化を計算した結果、計算値が実験値と良く対応することが示された。
|
