シグナルパスウェイネットワークの特性を推定する方法およびコンピュータプログラム
【課題】パスウェイネットワークの特性を効率よく推定できるようにする。
【解決手段】この方法は、パスウェイネットワークを構成する反応物と、反応物間の相互作用とをそれぞれノードおよびエッジとするグラフ表現を準備するステップ１６２と、シグナル伝達を、反応物間の相互作用に所定のパラメータ値を割当ててシミュレートした結果を検証して、所定の条件を充足する相互作用に対応するエッジのみを選択的に残すようにグラフを更新するステップ１６４、１６６と、ステップ１６４、１６６を、各パラメータに複数通りの値を設定して繰返し実行し、各パラメータの各値に対して、当該パラメータに対応するエッジであってかつ更新するステップの結果、グラフ中に残ったエッジの数の分布を推定するステップ１６８、１７２、１７４と、推定するステップで推定されたエッジの数の分布に基づいて、パラメータの値の範囲を推定するステップ１７６とを含む。
【解決手段】この方法は、パスウェイネットワークを構成する反応物と、反応物間の相互作用とをそれぞれノードおよびエッジとするグラフ表現を準備するステップ１６２と、シグナル伝達を、反応物間の相互作用に所定のパラメータ値を割当ててシミュレートした結果を検証して、所定の条件を充足する相互作用に対応するエッジのみを選択的に残すようにグラフを更新するステップ１６４、１６６と、ステップ１６４、１６６を、各パラメータに複数通りの値を設定して繰返し実行し、各パラメータの各値に対して、当該パラメータに対応するエッジであってかつ更新するステップの結果、グラフ中に残ったエッジの数の分布を推定するステップ１６８、１７２、１７４と、推定するステップで推定されたエッジの数の分布に基づいて、パラメータの値の範囲を推定するステップ１７６とを含む。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は細胞内におけるシグナル伝達のメカニズムを効率的に推定するための方法に関し、特に、細胞内におけるシグナル伝達のパスウェイネットワークの各経路において、クロストークが存在する場合のシグナル伝達の特性を推定するための方法に関する。
【背景技術】
【０００２】
生物細胞には、外部から与えられた信号を細胞核またはその他の細胞内要素に伝達するためのシグナル伝達の仕組みが備えられている。シグナル伝達は、酵素によるたんぱく質のリン化および脱リン化による活性化・不活性化に伴う連鎖的な生化学変化によって行なわれる。図１に、その中で低分子量ＧＴＰａｓｅ（ＧＴＰ結合たんぱく質）であるＲｈｏの活性化によるＭＬＣ（ミオシン軽鎖）のリン酸化を伴うシグナル伝達の例を示す。
【０００３】
図１を参照して、外部シグナルである特定の活性化因子により活性化されたＲｈｏ２２は、酵素であるＲｈｏキナーゼ２０と結合して、これを活性化させる（Ｒｈｏキナーゼ２４で示す。）。活性化されたＲｈｏキナーゼ２４はＭＬＣ２６をリン酸化し、リン酸化されたＭＬＣ３０が形成される。このようにしてＭＬＣのリン酸化レベルが上昇すると、細胞内にストレスファイバが形成されたりする結果になると考えられている。
【０００４】
このように、シグナル伝達にはたんぱく質のリン酸化および脱リン酸化が大きな役割を果たしている。ところで、あるたんぱく質と別のたんぱく質との間でのこのシグナル伝達の経路（パスウェイ）において、これらたんぱく質の間にどのような生化学変化上の関係があるかを考える。そのために、図２に、二種のたんぱく質（ｉ，ｊ）のリン酸化（ＰｉおよびＰｊ）、および脱リン酸化（ｄＰｉおよびｄＰｊ）の間での相互関係を示す。
【０００５】
図２を参照して、たんぱく質ｉのリン酸化４０およびたんぱく質ｊのリン酸化４２、たんぱく質ｉのリン酸化４０およびたんぱく質ｊの脱リン酸化４６、たんぱく質ｉの脱リン酸化４４およびたんぱく質ｊのリン酸化４２、ならびにたんぱく質ｉの脱リン酸化４４およびたんぱく質ｊの脱リン酸化４６の間には、それぞれ生化学変化により及ぼされる相互関係５０、５２、５４、および５６が存在する。これら４個の相互関係５０、５２、５４および５６はいずれも二方向の関係である。したがって、たんぱく質ｉとたんぱく質ｊとの間の一つのパスウェイは、４×２＝８＝２3個の関係によって定まる。
【０００６】
実際の細胞内では、シグナル伝達に関係するたんぱく質の種類は２種類ではなくはるかに多い。その場合、パスウェイを形成するたんぱく質相互の関係は非常に複雑になる。例えば、上記したＲｈｏに関連するＲｈｏファミリと、ＭＬＣおよびＭＢＳ（ミオシン結合サブユニット）との３種のたんぱく質の間での関係を図示すると図３のとおりとなる。
【０００７】
図３を参照して、Ｒｈｏ６０およびＭＬＣ６２、ＭＬＣ６２およびＭＢＳ６４、ならびにＭＢＳ６４およびＲｈｏ６０の間には、それぞれパスウェイ７０、７４および７２が考えられる。これらパスウェイ７０、７４および７２の各々について、図２に示すとおり、生化学変化と考慮すべき８＝２3個の相互関係がある。したがって、図３に示す３種のたんぱく質が関係するパスウェイネットワークでは、２3×２3×２3＝２9個、すなわち５１２個の相互関係が存在する。
【発明の開示】
【発明が解決しようとする課題】
【０００８】
上述したとおり、一つのシグナルパスウェイについても考慮すべき相互関係が多いのに加え、パスウェイには他の反応からの影響（クロストーク）も存在する。図３に示すＲｈｏ６０、ＭＬＣ６２、およびＭＢＳ６４の各々に対して２3通りのクロストークを考えると、図３のパスウェイネットワークで考慮すべき関係は２9×２9＝２18個、すなわち２６２，１４４となり、非常に複雑なものとなる。
【０００９】
実際のパスウェイネットワークでは、さらに多数のシグナルパスウェイが関係し、考慮しなければならない要素がさらに多く、さらに複雑になる。図４を参照して、図２に示したようなシグナルパスウェイを一つのモチーフで表すものとする。すなわち、あるたんぱく質１２０のリン酸化１２４が酵素などの要因１２２により生ずるものとし、脱リン酸化１２８が別の要因１２６により生じるものとする。すると、一つのモチーフにつき前述したとおり、考慮すべき関係が２9通り存在する。実際のパスウェイネットワーク９０が、複数個（ｎ個とする。）のモチーフ１００、１０２、１０４，…，１０８を含むものとすれば、考慮すべき関係は２3nとなる。このパスウェイネットワークにおいて、各モチーフにおける生化学変化の間の関係を推定することはＮＰ完全問題であって、解くことができるようなアルゴリズムが得られたとしても、その実行には莫大な時間がかかり事実上不可能である。一方で、例えば生体においてこうしたシグナル伝達の仕組みを用いて、ある種のたんぱく質の生成を防止したりすることができる薬剤などを開発しようとする場合、パスウェイネットワークの特性がどのようなものであるかを短時間に解析することが必要となる。
【００１０】
しかし、従来はそのようにパスウェイネットワークの特性を効率よく推定する方法はなく、したがってシグナル伝達の仕組みを利用した薬剤を効率よく得ることはできなかった。
【００１１】
それゆえに本発明の目的は、パスウェイネットワークの特性を効率よく推定することができる方法を提供することである。
【課題を解決するための手段】
【００１２】
本発明の第１の局面に係る方法は、シグナルパスウェイネットワークの特性を推定する方法であって、パスウェイネットワークを構成する反応物と、反応物間の相互作用とをそれぞれノードおよびエッジとするグラフ表現を準備するステップと、パスウェイネットワークにおけるシグナル伝達を、反応物間の相互作用に所定のパラメータ値を割当ててシミュレートした結果を検証して、所定の条件を充足する相互作用に対応するエッジのみを選択的に残すようにグラフを更新するステップと、グラフを更新するステップを、各パラメータに複数通りの値を設定して繰返し実行し、各パラメータの各値に対して、当該パラメータに対応するエッジであってかつ更新するステップの結果グラフ中に残ったエッジの数の分布を推定するステップと、推定するステップで推定されたエッジの数の分布に基づいて、パラメータの値の範囲を推定するステップとを含む。
【００１３】
シグナルパスウェイネットワークの特性を推定するにあたり、最初に当該シグナルパスウェイネットワークのグラフ表現を準備する。各エッジに、反応物間の相互作用に対するパラメータを割当て、各パラメータにある値を設定してシグナル伝達をシミュレートし、所定の条件を充足する相互作用に対応するエッジのみを残すようにグラフを残す。この操作を、パラメータの値を変えて複数回行なう。その結果、各パラメータごとに、各値に対しグラフ中に対応のエッジが残っている場合の数が得られ、その分布が推定される。分布が推定されるとその分布に対し、そのパラメータの値が、その分布にしたがって適切に選ばれた範囲に設定される。したがって、グラフを用いた有限回の操作により、シグナルパスウェイネットワークの特性を信頼性を持って推定できる。
【００１４】
好ましくは、反応物は、Ｒｈｏファミリのシグナルパスウェイを構成する反応物である。
【００１５】
さらに好ましくは、グラフを更新するステップは、反応物間の相互作用に所定のパラメータ値を割当ててパスウェイネットワークにおけるシグナル伝達をシミュレートするステップと、シミュレートするステップで得られた結果を検証し、クロストークのあるパスウェイの制約にしたがうパスウェイのみを選択的に残すようにグラフを更新するステップとを含む。
【００１６】
本発明の第２の局面に係るコンピュータプログラムは、コンピュータにより実行されると、上記したいずれかのシグナルパスウェイネットワークの特性を推定する方法の全てのステップを実行するようにコンピュータを制御するものである。
【発明を実施するための最良の形態】
【００１７】
以下、本発明の一実施の形態に係るパスウェイネットワークの特性の推定方法について図面を参照しながら説明する。以下の図において、同じ部品または処理工程（ステップ）については同じ参照番号を付し、それらの詳細な説明の繰返しは省略する。
【００１８】
一般に、図５に示したように、パスウェイを構成する反応物Ｘと反応物Ｙとの関係１４０、および反応物Ｘ’と反応物Ｙ’との関係１４２が存在する場合、式１４４のように各反応物の濃度ｘ、ｙ、ｘ’およびｙ’の時間的な変化の間には、所定の相互関係が存在する。例えば反応物Ｘの濃度ｘが減少するのに応じて反応物Ｙの濃度ｙが増加するような場合であれば、図５に示す式１４４のうち上の式の定数αは負となる。同様に、反応物Ｘ’の濃度ｘ’が増加すれば反応物Ｙ’の濃度ｙ’も増加するのであれば、式１４４の下式の定数βは正となる。こうした変化は一般に反応物のリン酸化および脱リン酸化という生化学反応を介して生じる。
【００１９】
本実施の形態では、パスウェイを構成するたんぱく質などの反応物の間に、上記したようなリン酸化および脱リン酸化と、それに伴う濃度の変化という相互関係が生じるものと想定する。そして、対象となるパスウェイネットワークを構成するたんぱく質相互の間に式１４４で示されるような、パラメータを含む関係を最初に想定し、それをグラフ表現する。すなわち、反応物をノード、反応物間の関係をエッジとするグラフ表現を考える。各エッジには、その相互関係の特性を表すパラメータが変数として割当てられる。
【００２０】
次に、各パラメータについてある値の組合せで、上記したグラフにしたがってシグナル伝達をシミュレートする。シミュレートした結果を、実際の実験結果を用いて検証し、上記した反応物間の相互関係の中で、クロストークによる制約を充足する相互関係のみ選択し、それ以外のあり得ない相互関係を取除く。そして、グラフ中ではその取り除かれた相互関係に対応するエッジを削除し、クロストークによる制約を充足するエッジを残してグラフを再構成する。その結果に基づいて、そのパラメータの値の組合せに対するグラフの各エッジの存在／不存在を記録する。
【００２１】
こうした処理を、パラメータの値の組合せを様々に変えて繰返し行ない、それらの各々に対しグラフの各エッジの存在／不存在を記録する。その結果、各パラメータに対し、その各値について、対応するエッジがグラフ中で存在していた場合の数が得られる。
【００２２】
パラメータの値の組合せについて、予め定めていた全ての値を用いて上記した処理を行なった後、その結果から各パラメータについて、その値に対し残っていたエッジの数の分布を推定し、その分布に対し統計的に適切と思われるパラメータの値の範囲を決定し、最初のパスウェイネットワーク内の各パスウェイ（グラフ内のエッジ）に特性として付与する。
【００２３】
以上が本実施の形態に係るパスウェイネットワーク特性の推定方法の概略である。この方法は、実際上、コンピュータを用いて実現することになる。図６に、そのためのコンピュータプログラムのメイン処理のフローチャートを示す。図６を参照して、まずステップ１６０で、パスウェイネットワークのパラメータの初期値を定める。次に、ステップ１６２において、予め準備されていた、パスウェイネットワークをグラフ表現する場合の論理形式を読込む。ステップ１６４で、このパスウェイネットワークで、設定されたパラメータを用いた各反応物の濃度変化をシミュレートする。ステップ１６６で、このシミュレーションの結果を、予め準備している実際の実験で部分的に得られた結果を用いて検証する。そして、実験結果から見てクロストークの制約を充足せず、存在しないと考えられる経路を特定し、その経路に対応するエッジをグラフから削除し、グラフを更新する。
【００２４】
更新したグラフの形状に基づき、パラメータ値と、グラフのパス形状とを記録する（ステップ１６８）。より具体的には、図７に示すように、予めパス形状とパラメータ値との記憶領域２００を設けておく。シミュレーションごとに、そのときのパラメータの値と、各エッジ（パス）の存在／不存在とを記録する。例えばパラメータαについての記憶領域２１０では、α＝α１，α２，…，αｍのときのパスＰＡＴＨ１〜ＰＡＴＨｎの存在／不存在が１／０の形式で記録される。α１〜αｍについては、互いに同じ値であることもあり得る。同様に他のパラメータについても記憶領域２１２，…，２２０にパラメータの値と、そのときのパスＰＡＴＨ１〜ＰＡＴＨｎの存在／不存在が記録される。
【００２５】
再び図６を参照して、ステップ１７０において、パラメータの値の全ての組合せについてステップ１６２〜１６８の処理が終了したか否かを判定する。もしも終了していれば制御はステップ１７４に進む。それ以外の場合には制御はステップ１７２に進む。
【００２６】
ステップ１７２では、パラメータの値の組合せを変えるために、特定のパラメータを所定の変化量で変化させる。その後ステップ１６２に戻る。
【００２７】
ステップ１７０で全てのパラメータの組合せに対して処理が終了したと判定されると、ステップ１７４で各パラメータの値に対して存在するパスの分布を推定する。具体的には、図８に示すように、例えばあるパラメータαの値ごとに、シミュレーションとその検証との結果、検証後のグラフにこのパラメータαに対応するエッジ（パス）が残っていたか否かを調べる。そして、パラメータαの各値に対して対応するエッジがどれだけのケースにおいて残っていたかをヒストグラム化する。そして、このヒストグラムからその分布曲線２４０がどのような分布を表しているかを推定する。例えば分布としては正規分布、ポアソン分布、ｔ分布、Ｆ分布などを想定し、分布曲線２４０がそれらのうちのどの分布を表しているかを決定する。
【００２８】
ステップ１７６において、推定された分布の型に応じ、それぞれパラメータの値の範囲として適切な範囲を分布曲線から決定する。例えば、図９に示すように、分布曲線２４０が正規分布を表し、その平均が一点鎖線２５０で示される値ｍであり、標準偏差がσと計算できたものとする。この場合には、αの値を、ｍ±３σの範囲２５２に決定する。このようにして各パラメータごとに決定した値の範囲をステップ１７８で出力する。すなわち、最初に想定されたパスウェイネットワークの各エッジに、それらエッジに対応するパラメータの範囲を割当てる。これにより、パスウェイネットワークの特性が推定されたことになる。
【００２９】
このようにして、本実施の形態に依れば、パスウェイネットワークの各シグナルパスについて、反応の特性を示すパラメータを決定できる。この決定のアルゴリズムのための手順は有限回で済み、かつその回数はパラメータの値の変化の選び方により制御できる。したがって、パスウェイネットワークのように複雑な系の振る舞いを特定する値を、比較的短い時間の間に推定することができる。このパラメータの値にしたがってシグナル伝達をシミュレートすることで、薬剤により細胞内に生じる変化を効率よく推定できる。
【００３０】
［コンピュータによる実現］
この実施の形態に係るパスウェイネットワーク特性の推定方法は、コンピュータハードウェアと、そのコンピュータハードウェアにより実行されるプログラムと、コンピュータハードウェアに格納されるデータとにより実現することができる。図１０はこのコンピュータシステム３３０の外観を示し、図１１はコンピュータシステム３３０の内部構成を示す。
【００３１】
図１０を参照して、このコンピュータシステム３３０は、ＦＤ（フレキシブルディスク）ドライブ３５２およびＣＤ−ＲＯＭ（コンパクトディスク読出専用メモリ）ドライブ３５０を有するコンピュータ３４０と、キーボード３４６と、マウス３４８と、モニタ３４２とを含む。
【００３２】
図１１を参照して、コンピュータ３４０は、ＦＤドライブ３５２およびＣＤ−ＲＯＭドライブ３５０に加えて、ＣＰＵ（中央処理装置）３５６と、ＣＰＵ３５６、ＦＤドライブ３５２およびＣＤ−ＲＯＭドライブ３５０に接続されたバス３６６と、ブートアッププログラム等を記憶する読出専用メモリ（ＲＯＭ）３５８と、バス３６６に接続され、プログラム命令、システムプログラム、および作業データ等を記憶するランダムアクセスメモリ（ＲＡＭ）３６０とを含む。コンピュータシステム３３０はさらに、プリンタ３４４を含んでいる。
【００３３】
ここでは示さないが、コンピュータ３４０はさらにローカルエリアネットワーク（ＬＡＮ）への接続を提供するネットワークアダプタボードを含んでもよい。
【００３４】
コンピュータシステム３３０に上記したパスウェイネットワーク特性の推定方法を実行させるためのコンピュータプログラムは、ＣＤ−ＲＯＭドライブ３５０またはＦＤドライブ３５２に挿入されるＣＤ−ＲＯＭ３６２またはＦＤ３６４に記憶され、さらにハードディスク３５４に転送される。または、プログラムは図示しないネットワークを通じてコンピュータ３４０に送信されハードディスク３５４に記憶されてもよい。プログラムは実行の際にＲＡＭ３６０にロードされる。ＣＤ−ＲＯＭ３６２から、ＦＤ３６４から、またはネットワークを介して、直接にＲＡＭ３６０にプログラムをロードしてもよい。
【００３５】
このプログラムは、コンピュータ３４０にこの実施の形態に係る方法を実行させる複数の命令を含む。ただし、この方法を行なわせるのに必要な基本的機能のいくつかはコンピュータ３４０上で動作するオペレーティングシステム（ＯＳ）またはサードパーティのプログラム、もしくはコンピュータ３４０にインストールされる各種ツールキットのモジュールにより提供される。従って、このプログラムはこの実施の形態に係る方法を実現するのに必要な機能全てを必ずしも含まなくてよい。このプログラムは、命令のうち、所望の結果が得られるように制御されたやり方で適切な機能または「ツール」を呼出すことにより、上記した方法を実行する命令のみを含んでいればよい。コンピュータシステム３３０の動作は周知であるので、ここでは繰り返さない。
【００３６】
今回開示された実施の形態は単に例示であって、本発明が上記した実施の形態のみに制限されるわけではない。本発明の範囲は、発明の詳細な説明の記載を参酌した上で、特許請求の範囲の各請求項によって示され、そこに記載された文言と均等の意味および範囲内でのすべての変更を含む。
【図面の簡単な説明】
【００３７】
【図１】シグナル伝達の一例を示す図である。
【図２】パスウェイにおいて考慮すべき反応物間の関係を模式的に示す図である。
【図３】Ｒｈｏファミリにおけるシグナルパスウェイのネットワークの例を示す図である。
【図４】パスウェイネットワーク特性の推定の複雑さを説明するための図である。
【図５】パスウェイネットワーク特性を推定するために想定される、反応物の濃度の時間的変化の間に想定される式およびパラメータの例を示す図である。
【図６】パスウェイネットワーク特性を推定するためのコンピュータプログラムのメインルーチンを示すフローチャートである。
【図７】パスの存在／不存在とパラメータの値との記憶領域を示す図である。
【図８】パラメータの値に対するパスの分布曲線の推定を示すグラフである。
【図９】推定された分布曲線からのパラメータの値の範囲の決定方法を説明するための図である。
【図１０】本発明の一実施の形態に係る方法およびシステムを実現するコンピュータシステムの外観図である。
【図１１】図１０に示すコンピュータのブロック図である。
【符号の説明】
【００３８】
２０，２４ Ｒｈｏキナーゼ、２２，６０ Ｒｈｏ、２６，６２ ＭＬＣ、３０ リン酸化されたＭＬＣ、６４ ＭＢＳ
【技術分野】
【０００１】
この発明は細胞内におけるシグナル伝達のメカニズムを効率的に推定するための方法に関し、特に、細胞内におけるシグナル伝達のパスウェイネットワークの各経路において、クロストークが存在する場合のシグナル伝達の特性を推定するための方法に関する。
【背景技術】
【０００２】
生物細胞には、外部から与えられた信号を細胞核またはその他の細胞内要素に伝達するためのシグナル伝達の仕組みが備えられている。シグナル伝達は、酵素によるたんぱく質のリン化および脱リン化による活性化・不活性化に伴う連鎖的な生化学変化によって行なわれる。図１に、その中で低分子量ＧＴＰａｓｅ（ＧＴＰ結合たんぱく質）であるＲｈｏの活性化によるＭＬＣ（ミオシン軽鎖）のリン酸化を伴うシグナル伝達の例を示す。
【０００３】
図１を参照して、外部シグナルである特定の活性化因子により活性化されたＲｈｏ２２は、酵素であるＲｈｏキナーゼ２０と結合して、これを活性化させる（Ｒｈｏキナーゼ２４で示す。）。活性化されたＲｈｏキナーゼ２４はＭＬＣ２６をリン酸化し、リン酸化されたＭＬＣ３０が形成される。このようにしてＭＬＣのリン酸化レベルが上昇すると、細胞内にストレスファイバが形成されたりする結果になると考えられている。
【０００４】
このように、シグナル伝達にはたんぱく質のリン酸化および脱リン酸化が大きな役割を果たしている。ところで、あるたんぱく質と別のたんぱく質との間でのこのシグナル伝達の経路（パスウェイ）において、これらたんぱく質の間にどのような生化学変化上の関係があるかを考える。そのために、図２に、二種のたんぱく質（ｉ，ｊ）のリン酸化（ＰｉおよびＰｊ）、および脱リン酸化（ｄＰｉおよびｄＰｊ）の間での相互関係を示す。
【０００５】
図２を参照して、たんぱく質ｉのリン酸化４０およびたんぱく質ｊのリン酸化４２、たんぱく質ｉのリン酸化４０およびたんぱく質ｊの脱リン酸化４６、たんぱく質ｉの脱リン酸化４４およびたんぱく質ｊのリン酸化４２、ならびにたんぱく質ｉの脱リン酸化４４およびたんぱく質ｊの脱リン酸化４６の間には、それぞれ生化学変化により及ぼされる相互関係５０、５２、５４、および５６が存在する。これら４個の相互関係５０、５２、５４および５６はいずれも二方向の関係である。したがって、たんぱく質ｉとたんぱく質ｊとの間の一つのパスウェイは、４×２＝８＝２3個の関係によって定まる。
【０００６】
実際の細胞内では、シグナル伝達に関係するたんぱく質の種類は２種類ではなくはるかに多い。その場合、パスウェイを形成するたんぱく質相互の関係は非常に複雑になる。例えば、上記したＲｈｏに関連するＲｈｏファミリと、ＭＬＣおよびＭＢＳ（ミオシン結合サブユニット）との３種のたんぱく質の間での関係を図示すると図３のとおりとなる。
【０００７】
図３を参照して、Ｒｈｏ６０およびＭＬＣ６２、ＭＬＣ６２およびＭＢＳ６４、ならびにＭＢＳ６４およびＲｈｏ６０の間には、それぞれパスウェイ７０、７４および７２が考えられる。これらパスウェイ７０、７４および７２の各々について、図２に示すとおり、生化学変化と考慮すべき８＝２3個の相互関係がある。したがって、図３に示す３種のたんぱく質が関係するパスウェイネットワークでは、２3×２3×２3＝２9個、すなわち５１２個の相互関係が存在する。
【発明の開示】
【発明が解決しようとする課題】
【０００８】
上述したとおり、一つのシグナルパスウェイについても考慮すべき相互関係が多いのに加え、パスウェイには他の反応からの影響（クロストーク）も存在する。図３に示すＲｈｏ６０、ＭＬＣ６２、およびＭＢＳ６４の各々に対して２3通りのクロストークを考えると、図３のパスウェイネットワークで考慮すべき関係は２9×２9＝２18個、すなわち２６２，１４４となり、非常に複雑なものとなる。
【０００９】
実際のパスウェイネットワークでは、さらに多数のシグナルパスウェイが関係し、考慮しなければならない要素がさらに多く、さらに複雑になる。図４を参照して、図２に示したようなシグナルパスウェイを一つのモチーフで表すものとする。すなわち、あるたんぱく質１２０のリン酸化１２４が酵素などの要因１２２により生ずるものとし、脱リン酸化１２８が別の要因１２６により生じるものとする。すると、一つのモチーフにつき前述したとおり、考慮すべき関係が２9通り存在する。実際のパスウェイネットワーク９０が、複数個（ｎ個とする。）のモチーフ１００、１０２、１０４，…，１０８を含むものとすれば、考慮すべき関係は２3nとなる。このパスウェイネットワークにおいて、各モチーフにおける生化学変化の間の関係を推定することはＮＰ完全問題であって、解くことができるようなアルゴリズムが得られたとしても、その実行には莫大な時間がかかり事実上不可能である。一方で、例えば生体においてこうしたシグナル伝達の仕組みを用いて、ある種のたんぱく質の生成を防止したりすることができる薬剤などを開発しようとする場合、パスウェイネットワークの特性がどのようなものであるかを短時間に解析することが必要となる。
【００１０】
しかし、従来はそのようにパスウェイネットワークの特性を効率よく推定する方法はなく、したがってシグナル伝達の仕組みを利用した薬剤を効率よく得ることはできなかった。
【００１１】
それゆえに本発明の目的は、パスウェイネットワークの特性を効率よく推定することができる方法を提供することである。
【課題を解決するための手段】
【００１２】
本発明の第１の局面に係る方法は、シグナルパスウェイネットワークの特性を推定する方法であって、パスウェイネットワークを構成する反応物と、反応物間の相互作用とをそれぞれノードおよびエッジとするグラフ表現を準備するステップと、パスウェイネットワークにおけるシグナル伝達を、反応物間の相互作用に所定のパラメータ値を割当ててシミュレートした結果を検証して、所定の条件を充足する相互作用に対応するエッジのみを選択的に残すようにグラフを更新するステップと、グラフを更新するステップを、各パラメータに複数通りの値を設定して繰返し実行し、各パラメータの各値に対して、当該パラメータに対応するエッジであってかつ更新するステップの結果グラフ中に残ったエッジの数の分布を推定するステップと、推定するステップで推定されたエッジの数の分布に基づいて、パラメータの値の範囲を推定するステップとを含む。
【００１３】
シグナルパスウェイネットワークの特性を推定するにあたり、最初に当該シグナルパスウェイネットワークのグラフ表現を準備する。各エッジに、反応物間の相互作用に対するパラメータを割当て、各パラメータにある値を設定してシグナル伝達をシミュレートし、所定の条件を充足する相互作用に対応するエッジのみを残すようにグラフを残す。この操作を、パラメータの値を変えて複数回行なう。その結果、各パラメータごとに、各値に対しグラフ中に対応のエッジが残っている場合の数が得られ、その分布が推定される。分布が推定されるとその分布に対し、そのパラメータの値が、その分布にしたがって適切に選ばれた範囲に設定される。したがって、グラフを用いた有限回の操作により、シグナルパスウェイネットワークの特性を信頼性を持って推定できる。
【００１４】
好ましくは、反応物は、Ｒｈｏファミリのシグナルパスウェイを構成する反応物である。
【００１５】
さらに好ましくは、グラフを更新するステップは、反応物間の相互作用に所定のパラメータ値を割当ててパスウェイネットワークにおけるシグナル伝達をシミュレートするステップと、シミュレートするステップで得られた結果を検証し、クロストークのあるパスウェイの制約にしたがうパスウェイのみを選択的に残すようにグラフを更新するステップとを含む。
【００１６】
本発明の第２の局面に係るコンピュータプログラムは、コンピュータにより実行されると、上記したいずれかのシグナルパスウェイネットワークの特性を推定する方法の全てのステップを実行するようにコンピュータを制御するものである。
【発明を実施するための最良の形態】
【００１７】
以下、本発明の一実施の形態に係るパスウェイネットワークの特性の推定方法について図面を参照しながら説明する。以下の図において、同じ部品または処理工程（ステップ）については同じ参照番号を付し、それらの詳細な説明の繰返しは省略する。
【００１８】
一般に、図５に示したように、パスウェイを構成する反応物Ｘと反応物Ｙとの関係１４０、および反応物Ｘ’と反応物Ｙ’との関係１４２が存在する場合、式１４４のように各反応物の濃度ｘ、ｙ、ｘ’およびｙ’の時間的な変化の間には、所定の相互関係が存在する。例えば反応物Ｘの濃度ｘが減少するのに応じて反応物Ｙの濃度ｙが増加するような場合であれば、図５に示す式１４４のうち上の式の定数αは負となる。同様に、反応物Ｘ’の濃度ｘ’が増加すれば反応物Ｙ’の濃度ｙ’も増加するのであれば、式１４４の下式の定数βは正となる。こうした変化は一般に反応物のリン酸化および脱リン酸化という生化学反応を介して生じる。
【００１９】
本実施の形態では、パスウェイを構成するたんぱく質などの反応物の間に、上記したようなリン酸化および脱リン酸化と、それに伴う濃度の変化という相互関係が生じるものと想定する。そして、対象となるパスウェイネットワークを構成するたんぱく質相互の間に式１４４で示されるような、パラメータを含む関係を最初に想定し、それをグラフ表現する。すなわち、反応物をノード、反応物間の関係をエッジとするグラフ表現を考える。各エッジには、その相互関係の特性を表すパラメータが変数として割当てられる。
【００２０】
次に、各パラメータについてある値の組合せで、上記したグラフにしたがってシグナル伝達をシミュレートする。シミュレートした結果を、実際の実験結果を用いて検証し、上記した反応物間の相互関係の中で、クロストークによる制約を充足する相互関係のみ選択し、それ以外のあり得ない相互関係を取除く。そして、グラフ中ではその取り除かれた相互関係に対応するエッジを削除し、クロストークによる制約を充足するエッジを残してグラフを再構成する。その結果に基づいて、そのパラメータの値の組合せに対するグラフの各エッジの存在／不存在を記録する。
【００２１】
こうした処理を、パラメータの値の組合せを様々に変えて繰返し行ない、それらの各々に対しグラフの各エッジの存在／不存在を記録する。その結果、各パラメータに対し、その各値について、対応するエッジがグラフ中で存在していた場合の数が得られる。
【００２２】
パラメータの値の組合せについて、予め定めていた全ての値を用いて上記した処理を行なった後、その結果から各パラメータについて、その値に対し残っていたエッジの数の分布を推定し、その分布に対し統計的に適切と思われるパラメータの値の範囲を決定し、最初のパスウェイネットワーク内の各パスウェイ（グラフ内のエッジ）に特性として付与する。
【００２３】
以上が本実施の形態に係るパスウェイネットワーク特性の推定方法の概略である。この方法は、実際上、コンピュータを用いて実現することになる。図６に、そのためのコンピュータプログラムのメイン処理のフローチャートを示す。図６を参照して、まずステップ１６０で、パスウェイネットワークのパラメータの初期値を定める。次に、ステップ１６２において、予め準備されていた、パスウェイネットワークをグラフ表現する場合の論理形式を読込む。ステップ１６４で、このパスウェイネットワークで、設定されたパラメータを用いた各反応物の濃度変化をシミュレートする。ステップ１６６で、このシミュレーションの結果を、予め準備している実際の実験で部分的に得られた結果を用いて検証する。そして、実験結果から見てクロストークの制約を充足せず、存在しないと考えられる経路を特定し、その経路に対応するエッジをグラフから削除し、グラフを更新する。
【００２４】
更新したグラフの形状に基づき、パラメータ値と、グラフのパス形状とを記録する（ステップ１６８）。より具体的には、図７に示すように、予めパス形状とパラメータ値との記憶領域２００を設けておく。シミュレーションごとに、そのときのパラメータの値と、各エッジ（パス）の存在／不存在とを記録する。例えばパラメータαについての記憶領域２１０では、α＝α１，α２，…，αｍのときのパスＰＡＴＨ１〜ＰＡＴＨｎの存在／不存在が１／０の形式で記録される。α１〜αｍについては、互いに同じ値であることもあり得る。同様に他のパラメータについても記憶領域２１２，…，２２０にパラメータの値と、そのときのパスＰＡＴＨ１〜ＰＡＴＨｎの存在／不存在が記録される。
【００２５】
再び図６を参照して、ステップ１７０において、パラメータの値の全ての組合せについてステップ１６２〜１６８の処理が終了したか否かを判定する。もしも終了していれば制御はステップ１７４に進む。それ以外の場合には制御はステップ１７２に進む。
【００２６】
ステップ１７２では、パラメータの値の組合せを変えるために、特定のパラメータを所定の変化量で変化させる。その後ステップ１６２に戻る。
【００２７】
ステップ１７０で全てのパラメータの組合せに対して処理が終了したと判定されると、ステップ１７４で各パラメータの値に対して存在するパスの分布を推定する。具体的には、図８に示すように、例えばあるパラメータαの値ごとに、シミュレーションとその検証との結果、検証後のグラフにこのパラメータαに対応するエッジ（パス）が残っていたか否かを調べる。そして、パラメータαの各値に対して対応するエッジがどれだけのケースにおいて残っていたかをヒストグラム化する。そして、このヒストグラムからその分布曲線２４０がどのような分布を表しているかを推定する。例えば分布としては正規分布、ポアソン分布、ｔ分布、Ｆ分布などを想定し、分布曲線２４０がそれらのうちのどの分布を表しているかを決定する。
【００２８】
ステップ１７６において、推定された分布の型に応じ、それぞれパラメータの値の範囲として適切な範囲を分布曲線から決定する。例えば、図９に示すように、分布曲線２４０が正規分布を表し、その平均が一点鎖線２５０で示される値ｍであり、標準偏差がσと計算できたものとする。この場合には、αの値を、ｍ±３σの範囲２５２に決定する。このようにして各パラメータごとに決定した値の範囲をステップ１７８で出力する。すなわち、最初に想定されたパスウェイネットワークの各エッジに、それらエッジに対応するパラメータの範囲を割当てる。これにより、パスウェイネットワークの特性が推定されたことになる。
【００２９】
このようにして、本実施の形態に依れば、パスウェイネットワークの各シグナルパスについて、反応の特性を示すパラメータを決定できる。この決定のアルゴリズムのための手順は有限回で済み、かつその回数はパラメータの値の変化の選び方により制御できる。したがって、パスウェイネットワークのように複雑な系の振る舞いを特定する値を、比較的短い時間の間に推定することができる。このパラメータの値にしたがってシグナル伝達をシミュレートすることで、薬剤により細胞内に生じる変化を効率よく推定できる。
【００３０】
［コンピュータによる実現］
この実施の形態に係るパスウェイネットワーク特性の推定方法は、コンピュータハードウェアと、そのコンピュータハードウェアにより実行されるプログラムと、コンピュータハードウェアに格納されるデータとにより実現することができる。図１０はこのコンピュータシステム３３０の外観を示し、図１１はコンピュータシステム３３０の内部構成を示す。
【００３１】
図１０を参照して、このコンピュータシステム３３０は、ＦＤ（フレキシブルディスク）ドライブ３５２およびＣＤ−ＲＯＭ（コンパクトディスク読出専用メモリ）ドライブ３５０を有するコンピュータ３４０と、キーボード３４６と、マウス３４８と、モニタ３４２とを含む。
【００３２】
図１１を参照して、コンピュータ３４０は、ＦＤドライブ３５２およびＣＤ−ＲＯＭドライブ３５０に加えて、ＣＰＵ（中央処理装置）３５６と、ＣＰＵ３５６、ＦＤドライブ３５２およびＣＤ−ＲＯＭドライブ３５０に接続されたバス３６６と、ブートアッププログラム等を記憶する読出専用メモリ（ＲＯＭ）３５８と、バス３６６に接続され、プログラム命令、システムプログラム、および作業データ等を記憶するランダムアクセスメモリ（ＲＡＭ）３６０とを含む。コンピュータシステム３３０はさらに、プリンタ３４４を含んでいる。
【００３３】
ここでは示さないが、コンピュータ３４０はさらにローカルエリアネットワーク（ＬＡＮ）への接続を提供するネットワークアダプタボードを含んでもよい。
【００３４】
コンピュータシステム３３０に上記したパスウェイネットワーク特性の推定方法を実行させるためのコンピュータプログラムは、ＣＤ−ＲＯＭドライブ３５０またはＦＤドライブ３５２に挿入されるＣＤ−ＲＯＭ３６２またはＦＤ３６４に記憶され、さらにハードディスク３５４に転送される。または、プログラムは図示しないネットワークを通じてコンピュータ３４０に送信されハードディスク３５４に記憶されてもよい。プログラムは実行の際にＲＡＭ３６０にロードされる。ＣＤ−ＲＯＭ３６２から、ＦＤ３６４から、またはネットワークを介して、直接にＲＡＭ３６０にプログラムをロードしてもよい。
【００３５】
このプログラムは、コンピュータ３４０にこの実施の形態に係る方法を実行させる複数の命令を含む。ただし、この方法を行なわせるのに必要な基本的機能のいくつかはコンピュータ３４０上で動作するオペレーティングシステム（ＯＳ）またはサードパーティのプログラム、もしくはコンピュータ３４０にインストールされる各種ツールキットのモジュールにより提供される。従って、このプログラムはこの実施の形態に係る方法を実現するのに必要な機能全てを必ずしも含まなくてよい。このプログラムは、命令のうち、所望の結果が得られるように制御されたやり方で適切な機能または「ツール」を呼出すことにより、上記した方法を実行する命令のみを含んでいればよい。コンピュータシステム３３０の動作は周知であるので、ここでは繰り返さない。
【００３６】
今回開示された実施の形態は単に例示であって、本発明が上記した実施の形態のみに制限されるわけではない。本発明の範囲は、発明の詳細な説明の記載を参酌した上で、特許請求の範囲の各請求項によって示され、そこに記載された文言と均等の意味および範囲内でのすべての変更を含む。
【図面の簡単な説明】
【００３７】
【図１】シグナル伝達の一例を示す図である。
【図２】パスウェイにおいて考慮すべき反応物間の関係を模式的に示す図である。
【図３】Ｒｈｏファミリにおけるシグナルパスウェイのネットワークの例を示す図である。
【図４】パスウェイネットワーク特性の推定の複雑さを説明するための図である。
【図５】パスウェイネットワーク特性を推定するために想定される、反応物の濃度の時間的変化の間に想定される式およびパラメータの例を示す図である。
【図６】パスウェイネットワーク特性を推定するためのコンピュータプログラムのメインルーチンを示すフローチャートである。
【図７】パスの存在／不存在とパラメータの値との記憶領域を示す図である。
【図８】パラメータの値に対するパスの分布曲線の推定を示すグラフである。
【図９】推定された分布曲線からのパラメータの値の範囲の決定方法を説明するための図である。
【図１０】本発明の一実施の形態に係る方法およびシステムを実現するコンピュータシステムの外観図である。
【図１１】図１０に示すコンピュータのブロック図である。
【符号の説明】
【００３８】
２０，２４ Ｒｈｏキナーゼ、２２，６０ Ｒｈｏ、２６，６２ ＭＬＣ、３０ リン酸化されたＭＬＣ、６４ ＭＢＳ
【特許請求の範囲】
【請求項１】
シグナルパスウェイネットワークの特性を推定する方法であって、
前記パスウェイネットワークを構成する反応物と、反応物間の相互作用とをそれぞれノードおよびエッジとするグラフ表現を準備するステップと、
前記パスウェイネットワークにおけるシグナル伝達を、反応物間の相互作用に所定のパラメータ値を割当ててシミュレートした結果を検証して、所定の条件を充足する相互作用に対応するエッジのみを選択的に残すように前記グラフを更新するステップと、
前記グラフを更新するステップを、各パラメータに複数通りの値を設定して繰返し実行し、各パラメータの各値に対して、当該パラメータに対応するエッジであってかつ更新するステップの結果前記グラフ中に残ったエッジの数の分布を推定するステップと、
前記推定するステップで推定されたエッジの数の分布に基づいて、パラメータの値の範囲を推定するステップとを含む、シグナルパスウェイネットワークの特性を推定する方法。
【請求項２】
前記反応物は、Ｒｈｏファミリのシグナルパスウェイを構成する反応物である、請求項１に記載のシグナルパスウェイネットワークの特性を推定する方法。
【請求項３】
前記グラフを更新するステップは、
反応物間の相互作用に所定のパラメータ値を割当てて前記パスウェイネットワークにおけるシグナル伝達をシミュレートするステップと、
前記シミュレートするステップで得られた結果を検証し、クロストークのあるパスウェイの制約にしたがうパスウェイのみを選択的に残すように前記グラフを更新するステップとを含む、請求項１または請求項２に記載のシグナルパスウェイネットワークの特性を推定する方法。
【請求項４】
コンピュータにより実行されると、請求項１〜請求項３のいずれかに記載のシグナルパスウェイネットワークの特性を推定する方法の全てのステップを実行するように前記コンピュータを制御する、コンピュータプログラム。
【請求項１】
シグナルパスウェイネットワークの特性を推定する方法であって、
前記パスウェイネットワークを構成する反応物と、反応物間の相互作用とをそれぞれノードおよびエッジとするグラフ表現を準備するステップと、
前記パスウェイネットワークにおけるシグナル伝達を、反応物間の相互作用に所定のパラメータ値を割当ててシミュレートした結果を検証して、所定の条件を充足する相互作用に対応するエッジのみを選択的に残すように前記グラフを更新するステップと、
前記グラフを更新するステップを、各パラメータに複数通りの値を設定して繰返し実行し、各パラメータの各値に対して、当該パラメータに対応するエッジであってかつ更新するステップの結果前記グラフ中に残ったエッジの数の分布を推定するステップと、
前記推定するステップで推定されたエッジの数の分布に基づいて、パラメータの値の範囲を推定するステップとを含む、シグナルパスウェイネットワークの特性を推定する方法。
【請求項２】
前記反応物は、Ｒｈｏファミリのシグナルパスウェイを構成する反応物である、請求項１に記載のシグナルパスウェイネットワークの特性を推定する方法。
【請求項３】
前記グラフを更新するステップは、
反応物間の相互作用に所定のパラメータ値を割当てて前記パスウェイネットワークにおけるシグナル伝達をシミュレートするステップと、
前記シミュレートするステップで得られた結果を検証し、クロストークのあるパスウェイの制約にしたがうパスウェイのみを選択的に残すように前記グラフを更新するステップとを含む、請求項１または請求項２に記載のシグナルパスウェイネットワークの特性を推定する方法。
【請求項４】
コンピュータにより実行されると、請求項１〜請求項３のいずれかに記載のシグナルパスウェイネットワークの特性を推定する方法の全てのステップを実行するように前記コンピュータを制御する、コンピュータプログラム。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【公開番号】特開２００６−５９２８５（Ｐ２００６−５９２８５Ａ）
【公開日】平成１８年３月２日（２００６．３．２）
【国際特許分類】
【出願番号】特願２００４−２４３１２９（Ｐ２００４−２４３１２９）
【出願日】平成１６年８月２４日（２００４．８．２４）
【国等の委託研究の成果に係る記載事項】（出願人による申告）平成１６年度独立行政法人情報通信研究機構、研究テーマ「人間情報コミュニケーションの研究開発」に関する委託研究、産業活力再生特別措置法第３０条の適用を受ける特許出願
【出願人】（３９３０３１５８６）株式会社国際電気通信基礎技術研究所 (905)
【公開日】平成１８年３月２日（２００６．３．２）
【国際特許分類】
【出願日】平成１６年８月２４日（２００４．８．２４）
【国等の委託研究の成果に係る記載事項】（出願人による申告）平成１６年度独立行政法人情報通信研究機構、研究テーマ「人間情報コミュニケーションの研究開発」に関する委託研究、産業活力再生特別措置法第３０条の適用を受ける特許出願
【出願人】（３９３０３１５８６）株式会社国際電気通信基礎技術研究所 (905)
[ Back to top ]