(57)【要約】      （修正有）【課題】  タンパク質の確からしい立体構造を効率的に探索することを可能とするタンパク質の立体構造探索技術を提供する。【解決手段】  タンパク質の連結構造の内、立体構造の探索に使用するための立体構造が既知な部分を決定し、このタンパク質の連結構造を、上記決定された立体構造既知部分を含むフラグメントに分割し、上記決定された立体構造既知部分についてはその立体構造を初期状態として採用し、他の部分については任意の立体構造を初期状態として採用し、動径分布の実験値を真の動径分布とし、逆モンテカルロ法を用いて、与えられた各フラグメントの立体構造から計算によって求められる動径分布が真の動径分布に近づくように、各フラグメントの立体構造を最適化する。【選択図】      図４
2006171982.tif 000002
【特許請求の範囲】
【請求項１】
タンパク質の連結構造の内、立体構造の探索に使用するための立体構造が既知な部分を決定し、
当該タンパク質の連結構造を、当該決定された立体構造既知部分を含むフラグメントに分割し、
当該決定された立体構造既知部分についてはその立体構造を初期状態として採用し、他の部分については任意の立体構造を初期状態として採用し、
動径分布の実験値を真の動径分布とし、逆モンテカルロ法を用いて、与えられた各フラグメントの立体構造から計算によって求められる動径分布が当該真の動径分布に近づくように、各フラグメントの立体構造を最適化する
タンパク質の立体構造探索方法。
【請求項２】
前記フラグメントの少なくともいずれか二つが、隣接するフラグメントと重複する部分を有するように、前記フラグメントへの分割を行う、請求項１に記載のタンパク質の立体構造探索方法。
【請求項３】
フラグメントへの分割をランダムに決定する、請求項１または２に記載のタンパク質の立体構造探索方法。
【請求項４】
コンピュータに、
入力されたタンパク質の連結構造の内、立体構造の探索に使用するための立体構造が既知な部分を決定する手順と、
当該タンパク質の連結構造を、当該決定された立体構造既知部分を含むフラグメントに分割する手順と、
当該決定された立体構造既知部分についてはその立体構造を初期状態として採用し、他の部分については任意の立体構造を初期状態として採用する手順と、
入力された動径分布の実験値を真の動径分布とし、逆モンテカルロ法を用いて、与えられた各フラグメントの立体構造から計算によって求められる動径分布が当該真の動径分布に近づくように、各フラグメントの立体構造を最適化する手順と
を含む手順を実行させるための、タンパク質の立体構造探索用コンピュータプログラム。
【請求項５】
タンパク質の連結構造、立体構造が既知な部分および動径分布の実験値を含む探索条件を入力する入力部と、
立体構造の探索に使用するための立体構造既知部分を決定し、当該タンパク質の連結構造を、当該決定された立体構造既知部分を含むフラグメントに分割する入力処理部と、
当該決定された立体構造既知部分についてはその立体構造を初期状態として採用し、他の部分については任意の立体構造を初期状態として採用して、当該動径分布の実験値を真の動径分布とし、逆モンテカルロ法を用いて、与えられた各フラグメントの立体構造から計算によって求められる動径分布が当該真の動径分布に近づくように、各フラグメントの立体構造を最適化する計算部と
を含む、タンパク質の立体構造探索装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、連結構造データの分かっているタンパク質の立体構造を探索する方法に関する。
【背景技術】
【０００２】
連結構造データの分かっているタンパク質の立体構造を探索する方法としては、そのエネルギー的側面に注目して、配列の指定されたペプチド鎖の自由エネルギー最小の立体構造を求める方法がある。
【０００３】
このうち、経験的な方法としてはホモロジーモデリングや３Ｄ−１Ｄ法などがあげられる。これらの方法は、予測すべき立体構造に類似の立体構造がすでに実験的に解かれていない場合には適用することができない。
【０００４】
このような制限のない非経験的な方法としては、シンプルなモデルを用いて明快かつ厳密な理論的解析を目指す格子モデル法、本来はタンパク質の動力学的性質を調べるためのモレキュラーダイナミクス法（ＭＤ法）、統計力学的性質を調べるためのモンテカルロ法（ＭＣ法）等を折り畳みシミュレーションに応用したもの等があげられる。
【０００５】
格子モデル法はタンパク質の立体構造をより単純な格子モデルに置き換え、様々なエネルギー計算やＭＣ法と組み合わせてシミュレーションを行う方法であって、理論的扱いが他の方法と比較して容易である。
【０００６】
ＭＤ法は、単純化されたモデルではない実際の分子立体構造に忠実なモデルを用いるシミュレーション方法であり、ニュートン力学に基づいたシミュレーションにより、タンパク質を含む系がとりうるエネルギーの低い立体構造を求める。
【０００７】
ＭＣ法（たとえば特許文献１および非特許文献１参照）も、単純化されたモデルではない実際の分子立体構造に忠実なモデルを用いるシミュレーション方法であり、計算上で分子を人為的に変形させそのエネルギー変化に応じて立体構造を変えていく。この分子変形は 現実とは必ずしも対応しないという意味では欠点となっているが、この変形によって自由エネルギーの低い立体構造の探索を素早く、広範囲に行えるという意味では利点となっている。
【特許文献１】特開２００３−２０６２４６号公報（特許請求の範囲）
【非特許文献１】ハンスマン（Ｈａｎｓｍａｎｎ）・Ｕ・Ｈ・Ｅ，オカモト・Ｙ，「タンパク質折り畳み構造に関する新モンテカルロ法アルゴリズム。構造生物学の最近の見解（Ｎｅｗ  Ｍｏｎｔｅ  Ｃａｒｌｏ  ａｌｇｏｒｉｔｈｍｓ  ｆｏｒ  ｐｒｏｔｅｉｎ  ｆｏｌｄｉｎｇ．  Ｃｕｒｒｅｎｔ  Ｏｐｉｎｉｏｎ  ｉｎ  Ｓｔｒｕｃｔｕｒａｌ  Ｂｉｏｌｏｇｙ）」，１９９年，第９巻，ｐ．１７７−１８３
【発明の開示】
【発明が解決しようとする課題】
【０００８】
これらの非経験的な方法は、いずれもペプチド鎖の自由エネルギー最小の立体構造を求めるものであるが、この自由エネルギーは理論計算に基づくものであり、実際のエネルギー極小点と一致する保証はない。しかも、実際の立体構造にはエネルギー極小点が多数存在するものと考えられ、これらの方法で立体構造を決定しようとする場合、理論計算では非常に多くの候補が挙げられてしまうことが多い。
【０００９】
本発明は、この問題を解決し、立体構造をより推定しやすい非経験的方法を提供することを目的としている。本発明のさらに他の目的および利点は、以下の説明から明らかになるであろう。
【課題を解決するための手段】
【００１０】
本発明の一態様によれば、タンパク質の連結構造の内、立体構造の探索に使用するための立体構造が既知な部分を決定し、このタンパク質の連結構造を、上記決定された立体構造既知部分を含むフラグメントに分割し、上記決定された立体構造既知部分についてはその立体構造を初期状態として採用し、他の部分については任意の立体構造を初期状態として採用し、動径分布の実験値を真の動径分布とし、逆モンテカルロ法を用いて、与えられた各フラグメントの立体構造から計算によって求められる動径分布が真の動径分布に近づくように、各フラグメントの立体構造を最適化するタンパク質の立体構造探索方法が提供される。
【００１１】
本発明の他の一態様によれば、コンピュータに、入力されたタンパク質の連結構造の内、立体構造の探索に使用するための立体構造が既知な部分を決定する手順と、このタンパク質の連結構造を、上記決定された立体構造既知部分を含むフラグメントに分割する手順と、上記決定された立体構造既知部分についてはその立体構造を初期状態として採用し、他の部分については任意の立体構造を初期状態として採用する手順と、入力された動径分布の実験値を真の動径分布とし、逆モンテカルロ法を用いて、与えられた各フラグメントの立体構造から計算によって求められる動径分布が真の動径分布に近づくように、各フラグメントの立体構造を最適化する手順とを含む手順を実行させるための、タンパク質の立体構造探索用コンピュータプログラムが提供される。
【００１２】
本発明の更に他の一態様によれば、タンパク質の連結構造、立体構造が既知な部分および動径分布の実験値を含む探索条件を入力する入力部と、立体構造の探索に使用するための立体構造既知部分を決定し、上記タンパク質の連結構造を、上記決定された立体構造既知部分を含むフラグメントに分割する入力処理部と、上記決定された立体構造既知部分についてはその立体構造を初期状態として採用し、他の部分については任意の立体構造を初期状態として採用して、動径分布の実験値を真の動径分布とし、逆モンテカルロ法を用いて、与えられた各フラグメントの立体構造から計算によって求められる動径分布が真の動径分布に近づくように、各フラグメントの立体構造を最適化する計算部とを含む、タンパク質の立体構造探索装置が提供される。
【００１３】
本発明の諸態様により、タンパク質の確からしい立体構造を効率的に探索することが可能となる。
【００１４】
上記いずれの態様においても、フラグメントの少なくともいずれか二つが、隣接するフラグメントと重複する部分を有するように、フラグメントへの分割を行うこと、フラグメントへの分割をランダムに決定すること、フラグメントへの分割を複数回行うこと、逆モンテカルロ法において、ペプチド構造の二面角ψおよび二面角φを変数として使用すること、が好ましい。
【発明の効果】
【００１５】
本発明により、タンパク質の確からしい立体構造を効率的に探索することが可能となる。
【発明を実施するための最良の形態】
【００１６】
以下に、本発明の実施の形態を図、実施例等を使用して説明する。なお、これらの図、実施例等および説明は本発明を例示するものであり、本発明の範囲を制限するものではない。本発明の趣旨に合致する限り他の実施の形態も本発明の範疇に属し得ることは言うまでもない。
【００１７】
本発明に係るタンパク質の立体構造探索方法では、タンパク質の連結構造の内、立体構造の探索に使用するための立体構造が既知な部分を決定し、このタンパク質の連結構造を、上記決定された立体構造既知部分を含むフラグメントに分割し、上記決定された立体構造既知部分についてはその立体構造を初期状態として採用し、他の部分については任意の立体構造を初期状態として採用し、動径分布の実験値を真の動径分布とし、逆モンテカルロ法を用いて、与えられた各フラグメントの立体構造から計算によって求められる動径分布が真の動径分布に近づくように、各フラグメントの立体構造を最適化する。
【００１８】
このようにして、実験データおよびタンパク質の既知構造情報をもとに、確からしい立体構造を効率的に精度よく探索することが可能となる。
【００１９】
一般に、乱数発生によって初期立体構造モデルを生成して逆モンテカルロ法で実験データを再現するような立体構造モデルを探索していくことは、タンパク質のような複雑な形状を持った数千個原子のものでは不可能であった。これに対し、本発明では、上記に述べたように実験データを真として、この真の値とモデルの値との差を小さくするように立体構造モデルを探索することで効率的に最適化された立体構造が得られる。
【００２０】
なお、本発明において「タンパク質の連結構造」とはタンパク質分子を構成する各原子の繋がる順序を意味する。結合手の長さおよび結合角が判明している場合が多いので、通常、これらも含めたものをタンパク質の連結構造とするが、不明の場合には、これらの因子をタンパク質の連結構造には含めず、逆モンテカルロ法を実施する場合の変数として扱うこともできる。
【００２１】
また、上記において、立体構造既知部分とは、その構造部分を含む、任意の構造体中でその構造部分の立体構造が既知であることを意味するものであり、検討対象であるタンパク質中の構造として既知であることを条件とするものではない。このことを図１を使用して説明する。図１の（ａ）では、タンパク質の連結構造を一本の線で表している。この内、Ｘの部分については、たとえば、図１の（ｂ）のように両端に−ＮＨ2と−ＣＯＯＨとを持つ構造について、その立体構造が知られていたとするとき、Ｘを立体構造既知部分と見なし、図１の（ｂ）のＸの部分の立体構造を、検討対象のタンパク質のＸ部分の初期条件として採用するのである。
【００２２】
本発明においては、このような立体構造を全て採用する必要はなく、収集された立体構造既知部分についての情報の中から立体構造の探索に実際に使用するための立体構造既知部分を選択決定して使用する。その決定方法については特に制限はなく、任意の方法を採用することができる。タンパク質によっては、更に別の分子または構造部分の一部としてＸ部分の立体構造が既知な場合もあり得る。このような場合は、そのいずれを採用するかも、任意に定めることができる。一般的に言えば、検討対象の連結構造との類似性の高い分子または構造部分のデータを採用することが好ましい。同一の個所について複数の立体構造を採用することもできる。
【００２３】
動径分布の実験値を真の動径分布として採用することにより、非経験的な方法であっても、実際に即した立体構造を迅速に探索することが可能となる。また、動径分布は一義的に定まり、エネルギー極小点のように複数の値を取らないため、この点でも立体構造を迅速に探索することが可能となる。動径分布の計算の基礎となる「与えられた各フラグメントの立体構造」は、最初は、初期状態の立体構造であり、その後は逆モンテカルロ法によって与えられる。
【００２４】
図２に本発明の原理図を示す。まず、図２の（ａ）のＡ，Ｂ，Ｃ，・・・に示すように、立体構造の探索に使用するための立体構造既知部分を決め、求めたいタンパク質の連結構造情報に対して、タンパク質の一分子をこれらの立体構造既知部分を含むフラグメントに分割する。この際、どのようにフラグメントに分割するかは任意に定めることができる。何らかの規則に従って自動的に分割してもよい。またランダムに行ってもよい。ランダムに行う方が、逆モンテカルロ法による探索がよりスムーズに行く場合が多い。一つのフラグメント中に含まれる立体構造既知部分は１個でも複数個でもよい。また立体構造既知部分を含まないフラグメントが存在していてもよい。
【００２５】
フラグメントへの分割は、一分子の全体について、一度におこなっても、フラグメントの立体構造の最適化を図りつつ、次のフラグメントの大きさを決めても、その他のどのような方法を採用してもよい。
【００２６】
また隣り合うフラグメントが重複する部分を有するようになっていてもよい。このようにするとセグメント同士の接続をより実状に近い状態として推定することができる。
【００２７】
なお、一般的には、検討対象であるタンパク質の一分子の全部をフラグメント化するが、部分的にフラグメント化しない部分を残してもよい。たとえば、検討対象であるタンパク質の一部として、立体構造が既知な部分あるいは立体構造がほぼ確実に解析されていると思われる部分については、フラグメント化せず、逆モンテカルロ法を実行する際に立体構造を固定しておいてもよい。
【００２８】
フラグメントへの分割は一種類に固定しておく必要はなく、逆モンテカルロ法を行う際に、新たな分割条件に変更することも可能である。分割条件の変更は、本発明に係るコンピュータプログラムにより自動的に行うこともできる。複数回分割を行うことで、最終的に得られる立体構造モデルの信頼性が向上する場合が多い。
【００２９】
フラグメントへの分割を複数回行うやり方も任意に定めることができる。一分子の全体のフラグメントを一時に変更しても、一部ずつ変更しても、その他のどのような方法を採用してもよい。
【００３０】
次に、立体構造の探索に使用するための立体構造既知部分についてはその立体構造を初期状態として採用し、他の部分については任意の立体構造を初期状態として採用して、初期の立体構造を構築する。図２の（ｂ）は、一分子の初期の立体構造を模式的に表した部分図である。この例では、立体構造既知部分以外は一分子の主鎖が一平面上にある状態を模式的に表している。
【００３１】
一方、中性子散乱等の実験結果を解析して得られた動径分布関数を真の動径分布関数として採用する。本発明では、図２の（ｃ）に示すように、逆モンテカルロ法を用いて、与えられた立体構造から計算によって求められる動径分布関数が真の動径分布に近づくように、各フラグメントの立体構造を最適化する。すなわち、計算された動径分布関数と真の動径分布関数とを比較する。そして、この一連の操作を十分な収束が得られるまで繰り返して、図２の（ｄ）に示すように確からしい立体構造の構造モデルを生成する。
【００３２】
所与の立体構造から動径分布関数を求める方法については、公知のどのような方法を採用してもよい。たとえば次のようにして計算できる。
【００３３】
【数１】
2006171982.tif 000003
逆モンテカルロ法で採用する条件はどのようなものでもよい。最も一般的な例として、図３に示すペプチド構造の主軸の二面角φとψとを変更することにより、各フラグメントの立体構造を変化させるやり方が挙げられる。なお、通常のペプチド結合においては、二面角ωに関わるＣ−Ｎ結合は二重結合性をおびており自由に回転できず、固定されることになるので、既知データとして入力することになる。二面角φ，ψ，ω、原子間の結合距離（結合手の長さ）および結合角が決まればタンパク質の立体構造が定まるが、本例ではω、原子間の結合距離（結合手の長さ）および結合角が既知であるので、逆モンテカルロ法により二面角φとψとを変更して、立体構造を最適化することになる。
【００３４】
このような立体構造探索方法は、コンピュータに、入力されたタンパク質の連結構造の内、立体構造の探索に使用するための立体構造が既知な部分を決定する手順と、このタンパク質の連結構造を、上記決定された立体構造既知部分を含むフラグメントに分割する手順と、上記決定された立体構造既知部分についてはその立体構造を初期状態として採用し、他の部分については任意の立体構造を初期状態として採用する手順と、入力された動径分布の実験値を真の動径分布とし、逆モンテカルロ法を用いて、与えられた各フラグメントの立体構造から計算によって求められる動径分布が真の動径分布に近づくように、各フラグメントの立体構造を最適化する手順とを含む手順を実行させるための、タンパク質の立体構造探索用コンピュータプログラムを使用して、コンピュータに行わせることができる。
【００３５】
また、上記のようにしてタンパク質の立体構造の探索を実行するためには、タンパク質の立体構造探索装置を使用することができる。この装置には、タンパク質の連結構造、立体構造が既知な部分および動径分布の実験値を含む探索条件を入力する入力部と、立体構造の探索に使用するための立体構造既知部分を決定し、上記タンパク質の連結構造を、上記決定された立体構造既知部分を含むフラグメントに分割する入力処理部と、上記決定された立体構造既知部分についてはその立体構造を初期状態として採用し、他の部分については任意の立体構造を初期状態として採用して、動径分布の実験値を真の動径分布とし、逆モンテカルロ法を用いて、与えられた各フラグメントの立体構造から計算によって求められる動径分布が真の動径分布に近づくように、各フラグメントの立体構造を最適化する計算部とが含まれる。
【００３６】
この様子を図５にブロック図として示してある。計算部には、図５に示すように、フラグメントを記憶するためのフラグメント記憶部、初期、中間および最終の立体構造を記憶しておくための構造データ記憶部、入力情報、中間情報、その他の諸情報を記憶するための各種記憶部や、タンパク質分子をフラグメントに分割するためのフラグメント分割処理部、どのフラグメントについて、逆モンテカルロ法によるシミュレーションを行うかを決めるフラグメント選択部、動径分布を計算するための動径分布計算部、計算された動径分布を真の動径分布と比較し、その差が充分小さくなったかどうかを判定するための比較判定部、逆モンテカルロ法によるシミュレーションを何回行ったかを判定するステップ数計数部、スィープ数計数部等を設けることができる。
【００３７】
これらの各部は独立に存在している必要はなく、一つの部で複数の機能を有していてもよい。なお、本発明では、一つのフラグメントについて、二面角ψとφとをいくつのグループに分けてシミュレーションを行う場合のそのグループ数と、一分子全体についてシミュレーションを行う場合のフラグメント数を掛け合わせたものを１モンテカルロステップと呼ぶ。
【００３８】
また、探索結果や中間および最終の立体構造を出力するための出力処理部および出力部が付設されていてもよい。
【００３９】
上記の好ましい方法を実行するためには、フラグメントの少なくともいずれか二つが、隣接するフラグメントと重複する部分を有するようにすること、フラグメントへの分割をランダムに決定すること、およびフラグメントへの分割を複数回行うことが可能であることが好ましいが、いずれもコンピュータプログラムにそのようなステップを設けることにより、あるいは、計算部にそのような機能を与えることにより、容易に実現することができる。
【００４０】
この様子を図４，５を使用して説明する。なお、この説明はあくまで例示的なものであり、他の方法を採用してもよいことは言うまでもない。図４は、本発明に係る方法を行うためのフローチャートであり、図５は、本発明に係るタンパク質の立体構造探索装置のブロック図である。
【００４１】
まず、図４のステップＳ１に従って、図５の入力部に探索条件を入力する。たとえば、タンパク質の連結構造、立体構造既知部分の情報および動径分布の実験値を含む探索条件を入力する。
【００４２】
例えば、タンパク質の連結構造の情報として、原子間の結合距離および結合角を入力する。また二面角φ，ψ，ωの内のωを既知のデータとして入力する。また、立体構造既知部分の情報としては、たとえば、原子間の結合距離、結合角、二面角φ，ψ，ωを入力する。
【００４３】
このほか、モンテカルロ法における条件を入力することも必要である。この例では、結合角と主鎖の二面角ωを固定し、二面角φとψとを変化させる。たとえば１０°ずつ変化させることが多いがその他の条件であってもよい。さらに、フラグメントへの分割を行う条件やどの立体構造既知部分の立体構造を採用するかの情報、立体構造既知部分以外の部分の初期の立体構造条件等を含めることも有用である。
【００４４】
ついで、図４のステップＳ２に従って、図５の入力処理部で、立体構造の探索に使用するための立体構造既知部分を決定し、タンパク質の連結構造を、上記決定された立体構造既知部分を含むフラグメントに分割する。
【００４５】
ついで、あるフラグメントについて、図５の計算部で、図４のステップＳ３に従って、決定された立体構造既知部分についてはその立体構造を初期状態として採用し、その他の部分については任意の立体構造を初期状態として採用する。立体構造は、具体的には、例えば、原子間の結合距離、結合角、二面角φ，ψ，ωで定義される。
【００４６】
次に、図４のステップＳ４に従って、各フラグメントの立体構造から計算によって動径分布を求め、図４のステップＳ５に従って、上記動径分布の実験値と比較する。そして、二面角φとψとをランダムに変化させて、各フラグメントの立体構造を変え、このようにして与えられた各フラグメントの立体構造から計算によって得られた動径分布が動径分布の実験値に対し一定の条件を満たすようになるまで探索を続ける。たとえば、計算によって得られた動径分布と動径分布の実験値との差が一定値以下になるまで、あるいは計算によって得られた動径分布が動径分布の実験値に向けて収束するまで探索を続ける。
【００４７】
具体的には、たとえば、図４のステップＳ６に従って、計算によって得られた動径分布が動径分布の実験値に対し収束しているかどうかを判定し、ＮＯであれば、図４のステップＳ７に従って、フラグメント内に１以上存在する二面角φとψとをランダムに変化させて、各フラグメントの立体構造を変え、図４のステップＳ４に戻って、各フラグメントの立体構造から計算によって動径分布を求める。
【００４８】
ステップＳ６でＹＥＳであれば、そのフラグメントの立体構造について探索が完了したものとして、次のフラグメントについて同様の探索を行う。
【００４９】
このようにして、実験データおよびタンパク質の既知構造情報をもとに、確からしい立体構造を効率的に探索することが可能となる。なお、全フラグメントについて、計算によって得られた動径分布が動径分布の実験値に向けて収束した後、あるいはその途中で、フラグメントの分割条件を変更してもよい。この変更は外部から入力することによって行ってもよく、タンパク質の立体構造探索装置が自動的に行うようにプログラムされていてもよい。
【００５０】
なお、上記に開示した内容から、下記の付記に示した発明が導き出せる。
【００５１】
（付記１）
タンパク質の連結構造の内、立体構造の探索に使用するための立体構造が既知な部分を決定し、
当該タンパク質の連結構造を、当該決定された立体構造既知部分を含むフラグメントに分割し、
当該決定された立体構造既知部分についてはその立体構造を初期状態として採用し、他の部分については任意の立体構造を初期状態として採用し、
動径分布の実験値を真の動径分布とし、逆モンテカルロ法を用いて、与えられた各フラグメントの立体構造から計算によって求められる動径分布が当該真の動径分布に近づくように、各フラグメントの立体構造を最適化する
タンパク質の立体構造探索方法。
【００５２】
（付記２）
前記フラグメントの少なくともいずれか二つが、隣接するフラグメントと重複する部分を有するように、前記フラグメントへの分割を行う、付記１に記載のタンパク質の立体構造探索方法。
【００５３】
（付記３）
フラグメントへの分割をランダムに決定する、付記１または２に記載のタンパク質の立体構造探索方法。
【００５４】
（付記４）
前記フラグメントへの分割を複数回行う、付記１〜３のいずれかに記載のタンパク質の立体構造探索方法。
【００５５】
（付記５）
前記逆モンテカルロ法において、ペプチド構造の二面角ψおよび二面角φを変数として使用する、付記１〜４のいずれかに記載のタンパク質の立体構造探索方法。
【００５６】
（付記６）
コンピュータに、
入力されたタンパク質の連結構造の内、立体構造の探索に使用するための立体構造が既知な部分を決定する手順と、
当該タンパク質の連結構造を、当該決定された立体構造既知部分を含むフラグメントに分割する手順と、
当該決定された立体構造既知部分についてはその立体構造を初期状態として採用し、他の部分については任意の立体構造を初期状態として採用する手順と、
入力された動径分布の実験値を真の動径分布とし、逆モンテカルロ法を用いて、与えられた各フラグメントの立体構造から計算によって求められる動径分布が当該真の動径分布に近づくように、各フラグメントの立体構造を最適化する手順と
を含む手順を実行させるための、タンパク質の立体構造探索用コンピュータプログラム。
【００５７】
（付記７）
前記フラグメントに分割する手順において、前記フラグメントの少なくともいずれか二つが、隣接するフラグメントと重複する部分を有するように分割が行われる、付記６に記載のコンピュータプログラム。
【００５８】
（付記８）
前記フラグメントに分割する手順において、フラグメントへの分割がランダムに決定される、付記６または７に記載のコンピュータプログラム。
【００５９】
（付記９）
前記フラグメントに分割する手順において、フラグメントへの分割が複数回行われる、付記６〜８のいずれかに記載のコンピュータプログラム。
【００６０】
（付記１０）
前記逆モンテカルロ法において、ペプチド構造の二面角ψおよび二面角φが変数として使用される、付記６〜９のいずれかに記載のコンピュータプログラム。
【００６１】
（付記１１）
タンパク質の連結構造、立体構造が既知な部分および動径分布の実験値を含む探索条件を入力する入力部と、
立体構造の探索に使用するための立体構造既知部分を決定し、当該タンパク質の連結構造を、当該決定された立体構造既知部分を含むフラグメントに分割する入力処理部と、
当該決定された立体構造既知部分についてはその立体構造を初期状態として採用し、他の部分については任意の立体構造を初期状態として採用して、当該動径分布の実験値を真の動径分布とし、逆モンテカルロ法を用いて、与えられた各フラグメントの立体構造から計算によって求められる動径分布が当該真の動径分布に近づくように、各フラグメントの立体構造を最適化する計算部と
を含む、タンパク質の立体構造探索装置。
【００６２】
（付記１２）
前記入力処理部が、前記フラグメントの少なくともいずれか二つが隣接するフラグメントと重複する部分を有するように、前記フラグメントへの分割を行う、付記１１に記載のタンパク質の立体構造探索装置。
【００６３】
（付記１３）
前記入力処理部が、フラグメントへの分割をランダムに決定する、付記１１または１２に記載のタンパク質の立体構造探索装置。
【００６４】
（付記１４）
前記入力処理部が、前記フラグメントへの分割を複数回行う、付記１１〜１３のいずれかに記載のタンパク質の立体構造探索装置。
【００６５】
（付記１５）
前記計算部が、前記逆モンテカルロ法において、ペプチド構造の二面角ψおよび二面角φを変数として使用する、付記１１〜１４のいずれかに記載のタンパク質の立体構造探索装置。
【図面の簡単な説明】
【００６６】
【図１】タンパク質の立体構造既知部分を説明するための模式図である。
【図２】本発明の原理を説明するための模式図である。
【図３】ペプチド構造の二面角φ，ψ，ωを説明するための模式図である。
【図４】本発明に係るタンパク質の立体構造探索方法を実施するためのフローチャートである。
【図５】本発明に係るタンパク質の立体構造探索装置のブロック図である。
【図１】
2006171982.tif 000004
【図２】
2006171982.tif 000005
【図３】
2006171982.tif 000006
【図４】
2006171982.tif 000007
【図５】
2006171982.tif 000008