(57)【要約】【課題】比較的シンプルな手順により、細胞内の反応について信頼性のあるシミュレーション結果を得ることが可能な生物科学シミュレータを提供する。【解決手段】シミュレーション方法は、細胞をシミュレーション空間であるフィールドとして定義するステップと、少なくともｍＲＮＡおよびタンパクを含む、細胞内の要素をオブジェクトとして定義するステップと、定義されたオブジェクトのうち、選択されたもの同士が近接した場合の反応を定義するステップと、定義されたオブジェクトの少なくとも一部を単位時間ごとに移動させるステップと、オブジェクトのうちの少なくとも複数のオブジェクトが所定の位置にある場合に、当該複数のオブジェクト間の干渉により、定義された反応のうち、所定の反応を生じさせるステップと、生じた反応にしたがって所定の処理を実行するステップとを備えている。【選択図】  図７
2005229805.tif 000002
【特許請求の範囲】
【請求項１】
細胞を、シミュレーション空間であるフィールドとして定義するステップと、
少なくともｍＲＮＡおよびタンパクを含む、細胞内の要素をオブジェクトとして定義するステップと、
定義されたオブジェクトのうち、選択されたもの同士が近接した場合の反応を定義するステップと、
定義されたオブジェクトの少なくとも一部を単位時間ごとに移動させるステップと、
前記オブジェクトのうちの少なくとも複数のオブジェクトが所定の位置にある場合に、当該複数のオブジェクト間の干渉により、前記定義された反応のうち、所定の反応を生じさせるステップと、
前記生じた反応にしたがって、所定の処理を実行するステップとを備えたことを特徴とする生命科学のシミュレーション方法。
【請求項２】
前記所定の反応が、新たなオブジェクトの発生、オブジェクトの消滅、および、オブジェクトの生成レートの制御の何れかであることを特徴とする請求項１に記載の生命科学のシミュレーション方法。
【請求項３】
前記所定の反応が、前記オブジェクトの、安定性、結合特性、酵素活性、および、細胞内局在の変化のうち、何れか一つであることを特徴とする請求項１または２に記載の生命科学のシミュレーション方法。
【請求項４】
前記干渉による所定の反応を生じさせるステップが、
フィールドを分割して得られた所定の範囲の各々において、複数のオブジェクトが存在するか否かを判定するステップと、
前記所定の範囲に複数のオブジェクトが存在する場合に、当該複数のオブジェクトによる反応が、定義された反応中に存在するか否かを判定するステップとを備えたことを特徴とする請求項１ないし３の何れか一項に記載のシミュレーション方法。
【請求項５】
前記所定の範囲が、前記フィールドを分割した矩形により画定されることを特徴とする請求項４に記載のシミュレーション方法。
【請求項６】
さらに、前記オブジェクトのうち、移動しないものをサイトとして、その位置を定義するステップを備えたことを特徴とする請求項１ないし５の何れか一項に記載のシミュレーション方法。
【請求項７】
細胞をシミュレーション空間であるフィールドとして定義し、当該細胞内の要素をオブジェクトとして定義し、かつ、定義されたオブジェクトのうち選択されたもの同士が近接した場合の反応を定義するシミュレーション条件設定手段と、
前記オブジェクトの少なくとも位置およびその寿命を含む情報を管理するオブジェクト管理テーブルを有するオブジェクト管理手段と、
前記定義された反応に関する種々の情報を管理するテーブルを有するモデル管理手段と、
前記オブジェクト管理部およびモデル管理部にて管理された情報に基づき、単位時間ごとに選択されたオブジェクトを移動して、オブジェクト間の反応をシミュレートするシミュレーション実行手段とを備えたことを特徴とする生命科学シミュレータ。
【請求項８】
前記定義された反応が、新たなオブジェクトの発生、オブジェクトの消滅、および、オブジェクトの生成レートの制御の何れかであることを特徴とする請求項７に記載の生命科学シミュレータ。
【請求項９】
前記定義された反応が、前記オブジェクトの、安定性、結合特性、酵素活性、および、細胞内局在の変化であることを特徴とする請求項７または８に記載の生命科学シミュレータ。
【請求項１０】
前記オブジェクト管理テーブルが、前記オブジェクトの各々の属性、座標、および、残存する寿命を含み、
シミュレーションにおいて単位時間が経過するごとに、前記残存する寿命が減じられるように構成されたことを特徴とする請求項７ないし９の何れか一項に記載の生命科学シミュレータ。
【発明の詳細な説明】
【０００１】
【産業上の技術分野】
本発明は、生命科学シミュレータに関し、より詳細には、分子生物的な現象を適切にシミュレートすることが可能な生命科学シミュレータに関する。
【０００２】
【従来の技術】
現在、シグナル伝達などの生命現象を、ミカエリス−メンテン式、素反応式或いはＨｉｌｌの方程式を利用し、微分方程式を、従来の解法（オイラー法、ルンゲ−クッタ法など）にて近似積分して、タンパクやｍＲＮＡの濃度を求める手法が知られている。
たとえば、ミカエリス−メンテン式（ｖ＝Ｖｍａｘ［Ｓ］／（［Ｓ］＋Ｋｍ）ｖ：反応速度、［Ｓ］：基質濃度、Ｖｍａｘ：最大反応速度、Ｋｍ：ミカエリス定数）に基づき、酵素反応のキネティクスを表すことができる。また、Ｈｉｌｌの式（ｖ＝Ｖｍａｘ［Ｓ］ｎ／Ｋｍ＋［Ｓ］ｎ）によりアロステリック酵素など親和性が変化する酵素反応を記述することができる。
【０００３】
さらに、素反応式を用いることにより、酵素反応以外の細胞内で起こる生化学反応を微分方程式の形で記述できる。
たとえば、
【式１】
2005229805.tif 000003
という平衡反応であれば、
ｄ［Ｃ］／ｄｔ＝ｋ１［Ａ］［Ｂ］−ｋ−１［Ｃ］
という微分方程式で表すことができる。
このようにして作られた微分方程式を連立して解を求めることにより、細胞内の要素の各々の濃度を計算することができる。
【０００４】
【発明が解決しようとする課題】
ミカエリス−メンテンの式などにおいては、そのパラメータであるＶｍａｘ、Ｋｍは、実測可能である。しかしながら、その式自体が多くの前提に基づいており、細胞内において上記式の前提となる仮定が成立しているか否かが確かでないという問題点があった。
【０００５】
また、素反応式におけるパラメータである反応速度定数は、平衡反応であれば求めることができない。非可逆反応では、その反応の自由エネルギーの値がわかれば、アレニウス式より求めることができる。しかしながら、細胞内でその値を求めることは不可能である。つまり、素反応式においては、入力すべきパラメータが実測不可能である。
【０００６】
すなわち、従来のシミュレータにおいては、細胞内にて式が成立するための条件を満たしているか否かを検証できず、したがって、式を解いてもその解により、細胞内の反応が正しく記述できているかどうかが明らかでなかったり、或いは、パラメータの実測ができず、正確な解自体が求められないなど、信頼性のあるシミュレーション結果を得ることが不可能であった。
さらに、従来のシミュレータは、細胞内の物質濃度に関する連立微分方程式を解くものであるため、細胞内におけるタンパク質などの局在化をうまく表現することが不可能であった。
【０００７】
本発明は、比較的シンプルな手順により、細胞内の反応について信頼性のあるシミュレーション結果を得ることが可能な生物科学シミュレータを提供することを目的とする。
【０００８】
【課題を解決するための手段】
本発明の目的は、細胞を、シミュレーション空間であるフィールドとして定義するステップと、少なくともｍＲＮＡおよびタンパクを含む、細胞内の要素をオブジェクトとして定義するステップと、定義されたオブジェクトのうち、選択されたもの同士が近接した場合の反応を定義するステップと、定義されたオブジェクトの少なくとも一部を単位時間ごとに移動させるステップと、前記オブジェクトのうちの少なくとも複数のオブジェクトが所定の位置にある場合に、当該複数のオブジェクト間の干渉により、前記定義された反応のうち、所定の反応を生じさせるステップと、前記生じた反応にしたがって、所定の処理を実行するステップとを備えたことを特徴とする生命科学のシミュレーション方法により達成される。
【０００９】
本発明によれば、細胞内のｍＲＮＡやタンパクをオブジェクトとして定義し、これらオブジェクトのうち所定のものを移動させて、オブジェクト相互間の干渉である所定の反応を生じさせている。すなわち、現実の細胞のモデルに非常に近い形態でシミュレーションを実行する。したがって、実測不可能なパラメータを使用したり、細胞内において数式が成立する条件を満たしているか否かを考えたりする必要なく、所望の細胞内の反応をシミュレートすることが可能となり、かつ、正確なシミュレート結果を得ることが可能となる。その際に、実際に実験的に得られたｍＲＮＡのコピー数などをパラメータとして入力し、これを利用することができる。
【００１０】
本発明の好ましい実施態様においては、所定の反応が、新たなオブジェクトの発生、オブジェクトの消滅、および、オブジェクトの生成レートの制御の何れかである。また、所定の反応として、オブジェクトの、安定性、結合特性、酵素活性、および、細胞内局在の変化のうち、何れか一つを定義しても良い。
【００１１】
本発明のさらに好ましい実施態様においては、干渉による所定の反応を生じさせるステップが、フィールドを分割して得られた所定の範囲の各々において、複数のオブジェクトが存在するか否かを判定するステップと、所定の範囲に複数のオブジェクトが存在する場合に、当該複数のオブジェクトによる反応が、定義された反応中に存在するか否かを判定するステップとを備えている。これにより、オブジェクトが相互に近接することによる反応を適切にシミュレートすることが可能となる。また、上記所定の範囲は、フィールドを分割した矩形により画定されているのがより望ましい。
本発明のさらに他の実施態様においては、オブジェクトのうち、移動しないものをサイトとして、その位置を定義している。
【００１２】
本発明の別の実施態様において、生物科学シミュレータは、細胞をシミュレーション空間であるフィールドとして定義し、当該細胞内の要素をオブジェクトとして定義し、かつ、定義されたオブジェクトのうち選択されたもの同士が近接した場合の反応を定義するシミュレーション条件設定手段と、前記オブジェクトの少なくとも位置およびその寿命を含む情報を管理するオブジェクト管理テーブルを有するオブジェクト管理手段と、前記定義された反応に関する種々の情報を管理するテーブルを有するモデル管理手段と、前記オブジェクト管理部およびモデル管理部にて管理された情報に基づき、単位時間ごとに選択されたオブジェクトを移動して、オブジェクト間の反応をシミュレートするシミュレーション実行手段とを備えている。
上記実施態様において、オブジェクト管理テーブルが、前記オブジェクトの各々の属性、座標、および、残存する寿命を含み、シミュレーションにおいて単位時間が経過するごとに、前記残存する寿命が減じられるように構成されているのが望ましい。
【００１３】
【発明の実施の形態】
以下、添付図面を参照して、本発明の実施の形態につき説明を加える。図１は、本発明の実施の形態にかかるシミュレータのハードウェア構成を示す図である。図１に示すように、シミュレータ１０は、ＣＰＵ１２と、メモリ１４と、外部記憶装置１６と、インタフェース（Ｉ／Ｆ）１８と、マウス、キーボードなどの入力装置２０と、ＣＲＴディスプレイやＬＣＤからなる表示装置２２と、プリンタ２４とから構成されている。ＣＰＵ１２、メモリ１４およびＩ／Ｆ１８は、データバス２６を介して接続されている。また、外部記憶装置１６、入力装置２０、表示装置２２およびプリンタ２４は、Ｉ／Ｆ１８を介してデータバス２６と接続でき、これにより、ＣＰＵ１２との間で、データの授受が可能となっている。
【００１４】
メモリ１４は、シミュレーション中のデータや、作動中のシミュレーションプログラムを一時的に記憶する主記憶を有している。また、外部記憶装置１６には、本実施の形態にかかるシミュレータ１０を作動させるのに必要なプログラムが記憶されている。上記プログラムは、ＣＤ−ＲＯＭやフロッピーディスクに格納され、ＣＤ−ＲＯＭ駆動装置やフロッピーディスク駆動装置（図示せず）により読み出され、メモリ１４の主記憶に記憶され、これにより、必要な処理が実行されるように構成しても良い。
【００１５】
図２は、本実施の形態にかかるシミュレータの機能を示すブロックダイヤグラムである。このシミュレータは、機能的にみると、シミュレーション空間内を移動するオブジェクトを管理するオブジェクト管理部２０２と、オブジェクトの生成およびオブジェクト間の干渉による反応を制御するモデル制御部２０４と、モデルの初期状態やシミュレーション時間を含むシミュレーション条件を設定するモデル設定部２０６と、各種条件にしたがってシミュレーションの計算を実行するシミュレーション実行部２０８と、モデル設定部２０６にて設定されたシミュレーション条件を記憶する条件記憶部２１０と、シミュレーション結果を記憶するシミュレーション記憶部２１２とを有している。
このように構成されたシミュレータにて実行されるシミュレーション例の原理につき以下に説明する。第1の実施の形態においては、生命科学シミュレーションの一例として、サーカディアンリズムのシミュレーションを実行する。このサーカディアンリズムは、以下のようなものであることが知られている。
【００１６】
例として、Ｄｒｏｓｏｐｈｉｌａにみられるサーカディアンリズム（日慨周期）について説明する。定常的な暗やみの中で飼われたＤｒｏｓｏｐｈｉａの細胞内には、２４．７時間という周期で自立的な周期を描く濃度変化をもつタンパク質、ｍＲＮＡが存在する。これらのタンパク質、ｍＲＮＡが自立的な周期を描くこと（サーカディアンリズム）については、詳しい研究が行なわれ、そのメカニズムにはフィードバックループという機構が存在していることが分かっている。このフィードバックループという機構には、実際に、多くの細胞内要素が関わっている。そこで、このサーカディアンリズムのシミュレーションにおいては、そのフィードバックループの重要な要素のみ考えたモデルを構築している。
具体的に、このモデル化されたフィードバックループの機構を説明する。このモデルには、２種類のｍＲＮＡと、そのから翻訳される２種類のタンパク質、そして、これをコードする２つのＤＮＡの転写サイト、および、ｍＲＮＡが翻訳されるリボソームにより構成される。ｍＲＮＡは、核内の転写サイトにおいて生成されるが、この生成されたｍＲＮＡは、細胞質に移動し、リボソーム上で翻訳されタンパク質を生成する。生成されたタンパク質は再び、核内に戻り、ｍＲＮＡの生成の場であるＤＮＡの転写サイトを制御する。その制御法には２種類あり、アクチベータタンパク質が転写サイトに入ると、一定時間、ｍＲＮＡの生成レートを増加させる。その一方、リプレッサタンパク質が転写サイトに入ると、ｍＲＮＡの生成を一定時間完全に停止させる。このような制御が、交互に転写サイトにかかることにより、それぞれのタンパク質、ｍＲＮＡの量に周期的変動が生じる。
【００１７】
図３は、本実施の形態にて用いられたモデルを示す図である。図３に示すように、細胞に対応するフィールド３００において、核３０１の内部には、転写サイト３０２、３０３（たとえば、符号３０２がアクチベータのｍＲＮＡの転写サイト、符号３０３がリプレッサのｍＲＮＡの転写サイト）が設けられ、その一方、核外にはリボソーム３０４が設けられる。アクチベータ３１０が、転写サイト３０２、３０３に入ることにより、当該転写サイトにおけるｍＲＮＡの生成レート（プロダクションレート）が増大する。ｍＲＮＡは、転写サイトにおける生成レートにしたがって、転写サイトにて作り出される。その一方、リプレッサ３１１は、何れかの転写サイト３０２、３０３に入ることによって、そのサイトにおける転写をある一定の時間だけ停止させる。これにより、そのサイトにて転写されるｍＲＮＡの数は一定時間だけ「０」となる。
【００１８】
フィールド３００の核外には、リボソーム３０４が設けられる。リボソーム３０４は、タンパク質生成の場であり、リボソーム３０４にｍＲＮＡが位置すると、当該ｍＲＮＡが示す情報にしたがったタンパクが生成される。たとえば、アクチベータのｍＲＮＡがリボソームに位置すると、所定量のアクチベータが生成される。リプレッサに関しても同様な現象が起こる。また、ｍＲＮＡは分解（degradation）により消滅するまで、タンパク生成の役目を果たす。
リボソーム３０４にて生成されたアクチベータやリプレッサが再度核内に移動すると、転写サイト３０２、３０３において上述したようなｍＲＮＡの生成レートの増加や生成の制限（停止）など種々の作用を果たす。
【００１９】
本実施の形態においては、このようなモデルを以下のように捉えている。なお、本明細書において、場合によっては、移動可能なオブジェクトのみを「オブジェクト」と称し、位置の固定されたオブジェクト（リボソームや転写サイト）を「サイト」と称する。まず、図４に示すように、シミュレーション空間であるフィールドを、所定の大きさ（２００×１００）の平面として捉えている。すなわち、（０，０）から（２００，１００）までの座標平面４００を設定して、この座標平面上を、アクチベータ、リプレッサおよびｍＲＮＡが移動すると考えた。また、上記座標平面４００上の所定の位置に、アクチベータのための転写サイト（第１の転写サイト）４０１、レセプタのための転写サイト（第２の転写サイト）４０２およびリボソーム４０３を配置した。本実施の形態においては、これら転写サイトを点（たとえば、第１の転写サイトは、（９０，５０））と考えた。無論、転写サイトをある平面領域と考えてよいことは言うまでもない。
【００２０】
また、アクチベータ、レセプタ、転写サイトなどの要素（物質）間の関係を規定し、シミュレーションを進めるようにしている。図５は、要素間の基本的な関係の一例を説明するための図である。図５（ａ）に示すように、物質Ａ（符号５０１参照）が物質Ｂ（符号５０２参照）の基質であり、その一方、物質Ｂが、物質Ａを物質Ｃに変える酵素であるとする。これらが、図５（ｂ）に示すように、空間的に（つまり座標上で）反応可能な位置にある（符号５０３参照）と、図５（ｃ）に示すように、物質Ａの情報が物質Ｃの情報に書き換えられ（符号５０４参照）、その一方、酵素である物質Ｂにおいては、その情報が維持される（符号５０５参照）。上記関係以外にも、種々の関係を規定することができる。
【００２１】
上記関係として、おおむね、以下の４つを考えることができる。すなわち、
（１）安定性
（２）結合特性
（３）酵素活性
（４）細胞内局在の変化
である。
安定性には、たとえば、ＲＮＡに、その分解酵素であるＲＮＡａｓｅが近接することにより、ＲＮＡの寿命を減少させるという反応がある。
結合特性には、転写因子ＣＲＥＢが、リン酸化酵素であるＡキナーゼと近接することにより、ＣＲＥＢがリン酸化され、ＣＲＥＢの特性が変化するという反応が知られている。
酵素活性として、インアクティブなＡキナーゼが、ｃＡＭＰと近接することにより、Ａキナーゼが活性化（アクティブ）される反応が知られている。
また、細胞内局在には、Ａキナーゼが、調節サブユニットと触媒サブユニットとから構成され、上記ｃＡＭＰと反応することにより、二つの調節サブユニットと２つの触媒サブユニットとに分けられ、それぞれが、核外と核内に局在するようになるような反応が知られている。
【００２２】
これらの反応を、後述するように定義しておくことにより、シミュレーションにおいて、オブジェクト間が近接したときに所望の反応を生じさせることが可能となる。
また、本実施の形態において、アクチベータ、レセプタおよびｍＲＮＡは、二次元のランダムウォークにて移動させている。ここで、単位時間あたりで移動できる距離は、「１」としている。したがって、ある時間「ｔ」において、アクチベータ等の要素（物質）が、座標（ｘｔ,ｙｔ）に位置していた場合に、次の時間「ｔ＋１」における要素の座標（ｘｔ＋１,ｙｔ＋１）は、以下のように表される。
ｘｔ＋１＝ｘｔ＋cosθ
ｙｔ＋１＝ｙｔ＋sinθ
（ここで、θは、０≦θ＜２πの乱数）
【００２３】
無論、上記ランダムウォークに限定されるものではなく、必要な場合には、異方性をもったランダムウォークにてオブジェクトを移動させても良いことは明らかである。
このような原理に基づくシミュレータ１０の作動について以下に説明を加える。図６は、本実施の形態にかかるシミュレータ１０のモデルやシミュレーション条件の設定に関する処理手順を示すフローチャートである。図６に示すように、まず、オペレータは、入力装置２０を操作して、シミュレーション空間のサイズや、オブジェクトの位置をシミュレータ１０に与える。シミュレータ１０のモデル設定部２０６は、これに応答して、シミュレーション空間サイズや所定のオブジェクト（サイト）の位置等を含むシミュレーションモデルのデータを生成し、これを条件記憶部２１０の所定の領域に記憶する（ステップ６０１、６０２）。ここでは、たとえば、移動しないオブジェクト（サイト：本実施の形態では、転写サイトやリボソーム）の位置が決定される。また、オブジェクト（サイト）を点ではなく所定の平面領域とする場合には、当該平面領域のサイズやその位置が決定され、これらデータが条件記憶部２１０に記憶される。
【００２４】
次いで、シミュレーション条件の入力に応答して、シミュレーション条件に関する種々のデータが生成され、これらデータも条件記憶部２１０の所定の領域に記憶される（ステップ６０３、６０４）。本実施の形態においては、アクチベータ、リプレッサ、アクチベータｍＲＮＡおよびリプレッサｍＲＮＡの４つのオブジェクトに関して、それぞれ、以下に述べるデータがシミュレーション条件として設定される。
（ａ）初期値：シミュレーション開始時の当該オブジェクトの数
（ｂ）寿命：オブジェクトが発生してから消失するまでのステップ数
（ｃ）発生サイト：当該オブジェクトを発生させるオブジェクト（サイト）、或いは、当該オブジェクトを発生させるオブジェクト（サイト）の位置
（ｄ）エフェクトサイト：当該オブジェクトが入ると、反応のイベントが生じるようなオブジェクト（サイト）、或いは、当該オブジェクト（サイト）の位置（ｅ）デフォルト生成レート：デフォルトでの当該オブジェクトの生成レート
（ｆ）反応：各オブジェクト間の直接的な制御
（ｇ）反応継続時間：上記反応が継続するステップ数
【００２５】
たとえば、本実施の形態において、アクチベータｍＲＮＡの発生サイトは「アクチベータ転写サイト（第１の転写サイト）」であり、エフェクトサイトは「リボソーム」であり、かつ、その反応は、「所定の数のアクチベータを発生させること」、より詳細には、「アクチベータの生成レートをある値にすること」となる。リプレッサｍＲＮＡの発生サイト、エフェクトサイト、および、反応は、それぞれ、「リプレッサ転写サイト（第２の転写サイト）」、「リボソーム」および「所定数のリプレッサを発生させること」となる。また、アクチベータおよびリプレッサに関しても、上記（ａ）〜（ｇ）に示すシミュレーション条件が設定される。表１は、アクチベータｍＲＮＡ、リプレッサｍＲＮＡ、アクチベータおよびリプレッサの、それぞれのシミュレーション条件の一例を示している。
【表１】
2005229805.tif 000004
また、オブジェクトの動き（ランダムウォーク）に異方性を持たせる場合には、その条件も与えておけばよい。
【００２６】
図７は、本実施の形態にかかるシミュレーションの処理の概略を示すフローチャートである。上述したシミュレーション条件等の設定が終了した後に、オペレータがシミュレーションを実行すべき旨の指示を与えると（ステップ７０１）、オブジェクト管理部２０２およびモデル制御部２０４が、条件記憶部２１０に記憶されたシミュレーション条件を読み出し、シミュレーション条件にしたがって、シミュレーション実行部２０８にて、シミュレーションが実行される（ステップ７０２以下参照）。ここで、オブジェクト管理部２０２は、条件記憶部２１０に記憶されたシミュレーション条件、特にオブジェクトに関するものにしたがって、図８に示すように、初期的或いはシミュレーション作動中に生成されるオブジェクトの各々の、座標、ライフタイム、属性、および、後述するオブジェクトの粒子の存在の有無を示す存在フラグを把握した管理リストを有し、管理リスト中の値をシミュレーション実行部２０８に伝達している。また、モデル制御部２０４は、後述するように、オブジェクト間の反応（相互作用）に関するデータなどを、シミュレーション実行部２０８に与える。
【００２７】
シミュレーション実行部７０２は、まず、オブジェクト管理部２０２から与えられたオブジェクトに関する情報にしたがって、当該オブジェクトを移動させる処理（オブジェクト処理）を実行する。図９は、オブジェクト処理をより詳細に示すフローチャートである。図９に示すように、この処理においては、まず、処理対象となるオブジェクトのライフタイムが判定される（ステップ９０１）。このライフタイムは、分解による消滅時間に対応する。したがって、このライフタイムが０であり、そのため、存在フラグが「１」となっているようなオブジェクト（ステップ９０１にてイエス(Yes)）は、もはや存在しないものであるため、座標を算出する必要がない。ライフタイムが０でない場合（ステップ９０１にてノー(No)）には、上述したランダムウォークを用いて、オブジェクトの新たな座標値が算出される（ステップ９０２）。次いで、当該オブジェクトに関するライフタイムがデクリメントされる（ステップ９０３）。このような処理の後、オブジェクト管理部２０２は、算出された座標値およびライフタイムを用いて、管理リストを更新する（ステップ９０４）。なお、ライフタイムのデクリメントにより、値が「０」になったオブジェクトに関しては、ステップ９０４において、その存在フラグが「１」となり、以下に述べるマトリクス処理が実行されないことに留意されたい。このような処理が、移動する全てのオブジェクトに関して実行される（ステップ９０５、９０６参照）。このように管理リストを用いた処理を行なうのは、メモリを有効に利用するためである。この処理のほかにも、自己参照構造体を用いたデータ形式の処理なども有効である。
【００２８】
上記オブジェクト処理が終了した後に、オブジェクト間の干渉である反応などを含むマトリクス処理が実行される。本実施の形態において、サイト（転写サイトやリボソーム）をフィールド上の点として表しているが、マトリクス処理において、当該サイトを、ある大きさを有する矩形（マトリクス）として把握し、移動後のオブジェクトが、マトリクス内に入ったときに、反応が生じると考えている。たとえば、図１０に示すように、本実施の形態において（ｘｊ，ｙｊ）（たとえば、ｘｉ、ｙｉは整数）に位置するサイト（符号１００１参照）に関して、当該サイトの位置を中心とする大きさ「１」のマトリクス（符号１０００参照）、移動後のオブジェクト（符号１００２参照）のｘ座標ｘｔが、ｘｊ≦ｘｔ＜ｘｊ＋１．０であり、かつ、そのｙ座標ｙｔが、ｙｊ≦ｙｔ＜ｙｊ＋１．０であれば、当該オブジェクトが、（ｘｊ、ｙｊ）に位置するサイトに近接し、反応が生じるとしている。
【００２９】
図１１は、モデル制御部２０４にて把握されているマトリクスに関する情報を説明するための図である。図１１に示すように、各マトリクスに関して、サイトが存在する場合には、その属性、反応による生成物、生成レート、デフォルト生成レート、および、生成レートが継続する時間が対応付けられている。したがって、モデル管理部２０４は、これらテーブルのデータを、必要に応じて、シミュレーション実行部２０８に供給する。
【００３０】
図１２は、マトリクス処理をより詳細に示すフローチャートである。シミュレーション実行部２１０は、最初のマトリクスから順次、以下の処理を実行する。まず、処理対象となるマトリクスの属性を調べ、何らかの属性（この実施の形態ではサイト）を有しているか否かを判断する（ステップ１２０１）。ステップ１２０１でイエス(Yes)と判断された場合には、当該マトリクス中に、オブジェクト処理により移動したオブジェクトが存在しているか否かを判断する（ステップ１２０３）。オブジェクトが存在する場合には、マトリクスの属性（つまり、どのサイトが位置しているか）と、オブジェクトの属性とに基づき、必要に応じて値等を変更する。
【００３１】
本実施の形態においては、前述したように、アクチベータｍＲＮＡがリボソームに近接した場合には、アクチベータの生成レートが所定の値だけ増大し、その一方、リプレッサｍＲＮＡがリボソームに近接した場合には、リプレッサの生成レートが所定の値だけ増大する。また、アクチベータが、何れかの転写サイトに近接した場合には、当該転写サイトにて生成されるｍＲＮＡ（アクチベータＭＲＮＡ或いはリプレッサｍＲＮＡ）の生成レートが所定の値だけ増大する。その一方、リプレッサが、何れかの転写サイトに近接した場合には、当該転写サイトにて生成されるｍＲＮＡ（アクチベータｍＲＮＡ或いはリプレッサｍＲＮＡ）の生成レートが、「０」となる。前述したように、これら作用が生じる時間も予め規定されており、モデル制御部２０４に把握されている（図１１の継続時間参照）。
【００３２】
したがって、モデル制御部２０４は、シミュレーション処理部２０８による処理に応答して、処理対象たるマトリクスに関するデータを更新する（ステップ１２０４）。なお、ステップ１２０３にてノー(No)と判断された場合には、データの更新は実行されず、直接、以下に説明するステップ１２０５に進む。
【００３３】
シミュレーション処理部２０８は、生成レート等にしたがって、当該マトリクスにて生成すべきオブジェクトを生成させる（ステップ１２０５）。たとえば、図１１に示すリプレッサ転写サイトにおいては、生成レートにしたがって、１０個のリプレッサｍＲＮＡが生成されることになる。オブジェクトの生成など反応に関する処理が終了すると、モデル管理部２０４は、当該マトリクスに関する継続時間をデクリメントする（ステップ１２０６）。これにより、継続時間が「０」、すなわち、反応が終了したと判断される場合には、生成レートがデフォルトレートに書き換えられる。さらに、生成されたオブジェクトに関するデータが、オブジェクト管理部２０２に伝達され、これらデータが管理テーブルに配置される（ステップ１２０７）。
【００３４】
これらステップ１２０２ないしステップ１２０７の処理が、すべてのオブジェクトに関して実行される（ステップ１２０８およびステップ１２０９参照）。
このようにして、オブジェクト処理（図７のステップ７０２）およびマトリクス処理（ステップ７０３）が終了すると、シミュレーション実行部２０８は、等がシミュレーション結果をシミュレーション結果記憶部２１２に記憶するとともに、シミュレーション結果を示す画像データを生成し、当該画像データを表示装置２２に出力する（ステップ７０４）。シミュレーションの終了（たとえば、予め設定されたシミュレーション時間が経過した場合や、シミュレーション停止を示すデータの受理）があった場合(ステップ７０５でイエス(Yes))には、処理が終了し、その一方、終了すべき場合でなければ（ステップ７０５でノー(No)）、シミュレーション時間をインクリメントして、次の時刻におけるオブジェクト処理が実行される。
【００３５】
図１３ないし図１８は、第１の実施の形態にかかるシミュレータを利用してシミュレーションを実行した場合の結果を示す図である。図１３は、シミュレーション中におけるフィールドの部分拡大図である。実際には、アクチベータ、リプレッサ、アクチベータｍＲＮＡおよびリプレッサｍＲＮＡなどのオブジェクトは、空間中の点（つまりある位置）を占めているが、図１３においては、オブジェクトが、当該空間中の点を含む単位矩形（すなわち１辺の長さが「１」である矩形）中に位置する場合に、当該矩形にオブジェクトが存在することを示す印（たとえば、「網掛け」や「ハッチング」）を付している。また、図中、黒塗りされた矩形は、リボソーム、アクチベータの転写サイトおよびリプレッサの転写サイトが、位置していることを示している。
【００３６】
図１４は、表１に示すシミュレーション条件を与えた場合のアクチベータおよびリプレッサの数を示す図である。図１４によれば、これらの数がそれぞれ周期性をもって増減していることが理解できる。その一方、表２に示すような初期条件を与えると、図１５に示すように、アクチベータおよびリプレッサの数が定常状態となる。
【表２】
2005229805.tif 000005
図１６ないし図１８は、フィールド中のオブジェクトの状態、並びに、アクチベータおよびリプレッサの数を含む画像の例である。図１６の例においては、表１と同様のシミュレーション条件にてシミュレーションを実行している。その一方、図１７の例においては、アクチベータｍＲＮＡの寿命（ライフタイム）を「２０００」に設定し、また、図１８の例においては、表２と同様に、アクチベータｍＲＮＡの寿命を「８０００」に設定した。
【００３７】
このように、本実施の形態によれば、細胞内の移動可能な要素（たとえば、アクチベータ、リプレッサ、ｍＲＮＡなど）をオブジェクトとして捉えるとともに、移動しない要素（たとえば、リボソームなど）を、上記オブジェクトと別種のオブジェクト（サイト）として捉え、オブジェクト同士の近接、オブジェクトとサイトとの近接により、所定の反応を生じさせている。したがって、細胞内の任意の要素に関する反応を定義すれば、所望のシミュレーションを実現することが可能となる。
【００３８】
次に、本発明の第２の実施の形態につき説明を加える。第１の実施の形態においては、移動するオブジェクトが、定められた位置のサイトに近接することで、所定の反応が生じるようなシミュレーションを実行したが、第２の実施の形態においては、移動するオブジェクト同士の反応が生じるようなシミュレーションを実現している。この実施の形態においても、シミュレータのハードウェア構成およびシミュレータの機能ブロックダイヤグラムは、図１および図２に示すものと同様である。
【００３９】
第２の実施の形態においては、移動するオブジェクト同士の近接による反応を実現するために、第２の実施の形態においては、以下のような手法をとっている。
シミュレーションモデルを設定する際（図６参照）に、フィールドのサイズ、オブジェクトの種類と数、存在する場合にはサイトの種類と数、それぞれのサイトの位置に加えて、反応判定関数を定義する。この反応判定関数とは、どのオブジェクトとどのオブジェクト（或いはサイト）とが反応するのか、また、反応によりどのような作用が生じるのかを示すものである。
【００４０】
たとえば、あるタンパク質「Ａ」とある酵素「Ｂ」とが近接することにより、タンパク質「Ｃ」が生じる場合、下記の式に示すものが反応判定関数となる。
タンパク質「Ａ」＋酵素「Ｂ」→たんぱく質「Ｃ」＋酵素「Ｂ」
反応は、オブジェクト（要素）の生成に限定されず、他のオブジェクトの生成レートの変更やオブジェクトやサイトの不活性化を示すものであっても良い。たとえば、第1の実施の形態において説明したサーカディアンリズムについて考えると、
アクチベータ「Ｄ」＋第１の転写サイト「Ｅ」
→アクチベータｍＲＮＡ「Ｆ」の生成レートの書き換え
リプレッサ「Ｇ」＋第１の転写サイト「Ｅ」
→第１の転写サイトの不活性化
というように捉えることもできる。
【００４１】
また、各オブジェクトにその種別を示すオブジェクトＩＤを付して、反応判定関数を、ＩＤを用いてあらわすのが望ましい。たとえば、上述したタンパク質「Ａ」と酵素「Ｂ」との反応を、（１）式に示すように定義しておく。
Ａ＋Ｂ→Ｃ＋Ｂ                      ・・・・（１）
このように、反応判定関数を定義しておくことにより、シミュレーション中のマトリクス処理にて、単位矩形中に位置するオブジェクトのオブジェクトＩＤを調べ、これらＩＤを有するオブジェクトからなる反応判定関数が存在するか否かを調べれば、オブジェクト間（或いはオブジェクトとサイトとの間）における反応を決定することが可能となる。
【００４２】
第２の実施の形態においても、シミュレーション処理の概略は図７に示すものと同様であり、オブジェクト処理も図９に示すものと同様である。図１９は、第２の実施の形態にかかるマトリクス処理を示すフローチャートである。図１９に示すように、シミュレーション実行部２１０は、処理対象となるマトリクス中に位置するオブジェクトのオブジェクトＩＤを検出する（ステップ１９０２）。次いで、当該マトリクス中に複数のマトリクスが存在すれば（ステップ１９０３でイエス(Yes)）、反応判定関数を参照して、存在した複数のオブジェクトのオブジェクトＩＤにて定義されているものが存在するか否かを判断する（ステップ１９０４）。たとえば、反応判定関数として、上述した（１）式が含まれる場合には、マトリクス中の複数のオブジェクトＩＤが、オブジェクトＩＤ「Ａ」および「Ｂ」を含むものであれば、このステップにて「関数あり」と判断される。
【００４３】
マトリクス中の複数のオブジェクトに関連する反応判定関数が存在する場合（ステップ１９０４にて「あり」）には、反応判定関数にて定義された反応を生じさせる処理を実行する（ステップ１９０５）。
その一方、ステップ１９０３にてノー(No)と判断され、或いは、ステップ１９０４にて「なし」と判断された場合であっても、当該マトリクスにて継続されているべき反応があれば、この反応に対応する処理を実行する（ステップ１９０６）。たとえば、一方のオブジェクトがサイトであれば、移動するオブジェクトとサイトが近接した後であっても、一定の期間だけ当該サイトにて生じる反応が継続されている場合がある。
【００４４】
このような処理の後に、シミュレーション処理部２１０は、反応継続時間をデクリメントする（ステップ１９０７）。処理によって生じたオブジェクトなどに関するデータは、オブジェクト管理部２０７に伝達され、これらデータが管理テーブルに配置される（ステップ１９０８）。上記ステップ１９０２〜ステップ１９０８の処理は、全てのマトリクスに関して実行される（ステップ１９０９およびステップ１９１０参照）。
【００４５】
このように本実施の形態によれば、各オブジェクトに、その種別を示すオブジェクトＩＤを付与し、かつ、マトリクスにて生じる反応を、オブジェクトＩＤを用いた反応判定関数にて定義したため、移動するオブジェクト同士が近接するような反応であっても、マトリクスに位置するオブジェクトを検出し、これらオブジェクトを含む反応判定関数を調べることにより、適切な反応を生じさせることが可能となる。
【００４６】
本発明は、以上の実施の形態に限定されることなく、特許請求の範囲に記載された発明の範囲内で、種々の変更が可能であり、それらも本発明の範囲内に包含されるものであることは言うまでもない。
たとえば、第１の実施の形態においては、サーカディアンリズムを引き起こすタンパク質等のうち、アクチベータ、リプレッサ、ｍＲＮＡを移動するオブジェクトと定義し、転写サイトおよびリボソームを移動しないオブジェクト（サイト）と定義してシミュレーションを実行しているが、このようなモデルに限定されない。上記モデルに、所定の酵素反応を加えたものを定義したモデルを作成し、これに基づくシミュレーションを実行できることはいうまでもない。
【００４７】
また、前記実施の形態において、シミュレーション空間であるフィールドを二次元としたが、これに限定されるものではなく、三次元としても良い。
さらに、前記実施の形態においては、オブジェクトの位置（座標）を移動させるとともに、オブジェクト間の近接或いはオブジェクトとサイトとの近接については、同一のマトリクス中に位置することにより判断している。しかしながら、このような手法に限定されるものではなく、種々の手法を利用しても良い。たとえば、オブジェクト間の距離が所定の閾値より小さくなった場合に、これらが近接したと判断しても良い。
【００４８】
また、本発明は、サーカディアンリズムに限定されず、細胞内の種々の要素（タンパク、酵素など）の反応をシミュレートするのに利用することができる。たとえば、上述した、安定性、結合特性、酵素活性および／または細胞内局在の変化を利用して、所望の反応を定義して、シミュレーションを実行することができる。
なお、本明細書において、手段とは必ずしも物理的手段を意味するものではなく、各手段の機能が、ソフトウェアによって実現される場合も無論包含する。また、一つの手段の機能が、二つ以上の物理的手段により実現されても、若しくは、二つ以上の手段の機能が、一つの物理的手段により実現されてもよい。
【００４９】
【発明の効果】
本発明によれば、比較的シンプルな手順により、細胞内の反応について信頼性のあるシミュレーション結果を得ることが可能な生物科学シミュレータを提供することが可能である。
【図面の簡単な説明】
【図１】  図１は、本発明の実施の形態にかかるシミュレータのハードウェア構成を示す図である。
【図２】  図２は、本実施の形態にかかるシミュレータの機能を示すブロックダイヤグラムである。
【図３】  図３は、本実施の形態にて用いられたモデルを示す図である。
【図４】  図４は、本実施の形態にかかるシミュレーション空間を説明するための図である。
【図５】  図５は、要素間の基本的な関係の一例を説明するための図である。
【図６】  図６は、本実施の形態にかかるシミュレータのモデルやシミュレーション条件の設定に関する処理手順を示すフローチャートである。
【図７】  図７は、本実施の形態にかかるシミュレーションの処理の概略を示すフローチャートである。
【図８】  図８は、第１の実施の形態にかかるオブジェクトの属性を示す図である。
【図９】  図９は、第1の実施の形態にかかるオブジェクト処理をより詳細に示すフローチャートである。
【図１０】  図１０は、第１の実施の形態にかかるマトリクス処理を説明するための図である。
【図１１】  図１１は、モデル制御部にて把握されているマトリクスに関する情報を説明するための図である。
【図１２】  図１２は、第１の実施の形態にかかるマトリクス処理をより詳細に示すフローチャートである。
【図１３】  図１３は、第１の実施の形態にかかるシミュレータを利用してシミュレーションを実行した場合の結果を示す図である。
【図１４】  図１４は、第１の実施の形態にかかるシミュレータを利用してシミュレーションを実行した場合の結果を示す図である。
【図１５】  図１５は、第１の実施の形態にかかるシミュレータを利用してシミュレーションを実行した場合の結果を示す図である。
【図１６】  図１６は、第１の実施の形態にかかるシミュレータを利用してシミュレーションを実行した場合の結果を示す図である。
【図１７】  図１７は、第１の実施の形態にかかるシミュレータを利用してシミュレーションを実行した場合の結果を示す図である。
【図１８】  図１８は、第１の実施の形態にかかるシミュレータを利用してシミュレーションを実行した場合の結果を示す図である。
【図１９】  図１９は、第２の実施の形態にかかるマトリクス処理を示すフローチャートである。
【符号の説明】
１０      シミュレータ
１２      ＣＰＵ
１４      メモリ
１６      外部記憶装置
１８      Ｉ／Ｆ
２０      入力装置
２２      表示装置
２４      プリンタ
２０２    オブジェクト管理部
２０４    モデル制御部
２０６    モデル設定部
【図１】
2005229805.tif 000006
【図２】
2005229805.tif 000007
【図３】
2005229805.tif 000008
【図４】
2005229805.tif 000009
【図５】
2005229805.tif 000010
【図６】
2005229805.tif 000011
【図７】
2005229805.tif 000012
【図８】
2005229805.tif 000013
【図９】
2005229805.tif 000014
【図１０】
2005229805.tif 000015
【図１１】
2005229805.tif 000016
【図１２】
2005229805.tif 000017
【図１３】
2005229805.tif 000018
【図１４】
2005229805.tif 000019
【図１５】
2005229805.tif 000020
【図１６】
2005229805.tif 000021
【図１７】
2005229805.tif 000022
【図１８】
2005229805.tif 000023
【図１９】
2005229805.tif 000024