(57)【要約】  計算機を用いて、標的分子と、複数の薬候補分子の各々との相互作用エネルギーの各々を算出する、相互作用エネルギーの算出方法であって、  前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの少なくともいずれかに影響を与えるが、前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの全てには影響を与えない計算対象水分子を抽出する工程と、  前記相互作用エネルギーの各々の算出において陽に扱う水分子として前記計算対象水分子のみを考慮し、かつ、前記計算対象水分子が影響を与える相互作用エネルギーの算出においてのみ前記計算対象水分子を考慮して、前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの各々を算出する工程と、を含む相互作用エネルギーの算出方法である。
2016103336.tif 000002
【特許請求の範囲】
【請求項１】
計算機を用いて、標的分子と、複数の薬候補分子の各々との相互作用エネルギーの各々を算出する、相互作用エネルギーの算出方法であって、
前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの少なくともいずれかに影響を与えるが、前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの全てには影響を与えない計算対象水分子を抽出する工程と、
前記相互作用エネルギーの各々の算出において陽に扱う水分子として前記計算対象水分子のみを考慮し、かつ、前記計算対象水分子が影響を与える相互作用エネルギーの算出においてのみ前記計算対象水分子を考慮して、前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの各々を算出する工程と、
を含むことを特徴とする相互作用エネルギーの算出方法。
【請求項２】
前記計算対象水分子を抽出する工程が、
前記標的分子と、前記複数の薬候補分子における一の薬候補分子とを含む一の相互作用算出対象空間における、
前記標的分子、及び前記一の薬候補分子の双方と相互作用する相互作用水分子を抽出する第１の処理と、
前記相互作用水分子のうちから、前記相互作用水分子の座標の揺らぎ幅、及び前記相互作用水分子の存在頻度の少なくともいずれかに基づいて、計算候補水分子を抽出する第２の処理と、
を含み、
更に、前記計算対象水分子を抽出する工程が、
前記標的分子と、前記複数の薬候補分子の各々とを含む全ての相互作用算出対象空間のうちの少なくとも前記計算候補水分子が存在する複数の相互作用算出対象空間を対比して、前記計算候補水分子のうちから、前記全ての相互作用算出対象空間において同じ座標にはない水分子を抽出する処理、
を含む請求項１に記載の相互作用エネルギーの算出方法。
【請求項３】
前記第２の処理が、前記相互作用水分子のうちから、前記相互作用水分子の座標の揺らぎ幅、及び前記相互作用水分子の存在頻度に基づいて、前記計算候補水分子を抽出する処理である、請求項２に記載の相互作用エネルギーの算出方法。
【請求項４】
前記相互作用水分子の座標の揺らぎ幅、及び前記相互作用水分子の存在頻度が、分子動力学シミュレーションにより求められる、請求項２から３のいずれかに記載の相互作用エネルギーの算出方法。
【請求項５】
前記相互作用エネルギーの各々を算出する工程が、
前記標的分子と、前記複数の薬候補分子における一の薬候補分子との一の相互作用エネルギーを算出する際に、
前記標的分子と、前記一の薬候補分子との複合体における、前記標的分子と、前記一の薬候補分子とのペア相互作用エネルギーを算出する処理Ａと、
前記複合体における前記標的分子の、前記標的分子単体に対する分極相互作用エネルギーを算出する処理Ｂと、
前記複合体における前記一の薬候補分子の、前記一の薬候補分子単体に対する分極相互作用エネルギーを算出する処理Ｃと、
を含み、
更に、前記標的分子と、前記一の薬候補分子との前記一の相互作用エネルギーを算出する際に、前記計算対象水分子を考慮するときには、前記計算対象水分子と、前記一の薬候補分子とのペア相互作用エネルギーを算出する処理Ｄを含む、
請求項１から４のいずれかに記載の相互作用エネルギーの算出方法。
【請求項６】
前記標的分子と、前記一の薬候補分子との前記相互作用エネルギーを算出する際に、前記計算対象水分子を考慮するときには、前記計算対象水分子が、前記処理Ａ、前記処理Ｂ、及び前記処理Ｃにおいて考慮される、請求項５に記載の相互作用エネルギーの算出方法。
【請求項７】
コンピュータに、
標的分子と、複数の薬候補分子の各々との相互作用エネルギーの少なくともいずれかに影響を与えるが、前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの全てには影響を与えない計算対象水分子を抽出する工程と、
前記相互作用エネルギーの各々の算出において陽に扱う水分子として前記計算対象水分子のみを考慮し、かつ、前記計算対象水分子が影響を与える相互作用エネルギーの算出においてのみ前記計算対象水分子を考慮して、前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの各々を算出する工程と、
を実行させることを特徴とするプログラム。
【請求項８】
前記計算対象水分子を抽出する工程が、
前記標的分子と、前記複数の薬候補分子における一の薬候補分子とを含む一の相互作用算出対象空間における、
前記標的分子、及び前記一の薬候補分子の双方と相互作用する相互作用水分子を抽出する第１の処理と、
前記相互作用水分子のうちから、前記相互作用水分子の座標の揺らぎ幅、及び前記相互作用水分子の存在頻度の少なくともいずれかに基づいて、計算候補水分子を抽出する第２の処理と、
を含み、
更に、前記計算対象水分子を抽出する工程が、
前記標的分子と、前記複数の薬候補分子の各々とを含む全ての相互作用算出対象空間のうちの少なくとも前記計算候補水分子が存在する複数の相互作用算出対象空間を対比して、前記計算候補水分子のうちから、前記全ての相互作用算出対象空間において同じ座標にはない水分子を抽出する処理、
を含む請求項７に記載のプログラム。
【請求項９】
前記第２の処理が、前記相互作用水分子のうちから、前記相互作用水分子の座標の揺らぎ幅、及び前記相互作用水分子の存在頻度に基づいて、前記計算候補水分子を抽出する処理である、請求項８に記載のプログラム。
【請求項１０】
前記相互作用水分子の座標の揺らぎ幅、及び前記相互作用水分子の存在頻度が、分子動力学シミュレーションにより求められる、請求項７から８のいずれかに記載のプログラム。
【請求項１１】
前記相互作用エネルギーの各々を算出する工程が、
前記標的分子と、前記複数の薬候補分子における一の薬候補分子との一の相互作用エネルギーを算出する際に、
前記標的分子と、前記一の薬候補分子との複合体における、前記標的分子と、前記一の薬候補分子とのペア相互作用エネルギーを算出する処理Ａと、
前記複合体における前記標的分子の、前記標的分子単体に対する分極相互作用エネルギーを算出する処理Ｂと、
前記複合体における前記一の薬候補分子の、前記一の薬候補分子単体に対する分極相互作用エネルギーを算出する処理Ｃと、
を含み、
更に、前記標的分子と、前記一の薬候補分子との前記一の相互作用エネルギーを算出する際に、前記計算対象水分子を考慮するときには、前記計算対象水分子と、前記一の薬候補分子とのペア相互作用エネルギーを算出する処理Ｄを含む、
請求項７から１０のいずれかに記載のプログラム。
【請求項１２】
前記標的分子と、前記一の薬候補分子との前記相互作用エネルギーを算出する際に、前記計算対象水分子を考慮するときには、前記計算対象水分子が、前記処理Ａ、前記処理Ｂ、及び前記処理Ｃにおいて考慮される、請求項１１に記載のプログラム。
【請求項１３】
標的分子と、複数の薬候補分子の各々との相互作用エネルギーの少なくともいずれかに影響を与えるが、前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの全てには影響を与えない計算対象水分子を抽出する抽出部と、
前記相互作用エネルギーの各々の算出において陽に扱う水分子として前記計算対象水分子のみを考慮し、かつ、前記計算対象水分子が影響を与える相互作用エネルギーの算出においてのみ前記計算対象水分子を考慮して、前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの各々を算出する算出部と、
を備えることを特徴とする相互作用エネルギーの算出装置。
【請求項１４】
前記抽出部が、
前記標的分子と、前記複数の薬候補分子における一の薬候補分子とを含む一の相互作用算出対象空間における、
前記標的分子、及び前記一の薬候補分子の双方と相互作用する相互作用水分子を抽出する第１の処理と、
前記相互作用水分子のうちから、前記相互作用水分子の座標の揺らぎ幅、及び前記相互作用水分子の存在頻度の少なくともいずれかに基づいて、計算候補水分子を抽出する第２の処理と、
を行い、
更に、前記抽出部が、
前記標的分子と、前記複数の薬候補分子の各々とを含む全ての相互作用算出対象空間のうちの少なくとも前記計算候補水分子が存在する複数の相互作用算出対象空間を対比して、前記計算候補水分子のうちから、前記全ての相互作用算出対象空間において同じ座標にはない水分子を抽出する処理、
を行う請求項１３に記載の相互作用エネルギーの算出装置。
【請求項１５】
前記第２の処理が、前記相互作用水分子のうちから、前記相互作用水分子の座標の揺らぎ幅、及び前記相互作用水分子の存在頻度に基づいて、前記計算候補水分子を抽出する処理である、請求項１４に記載の相互作用エネルギーの算出装置。
【請求項１６】
前記相互作用水分子の座標の揺らぎ幅、及び前記相互作用水分子の存在頻度が、分子動力学シミュレーションにより求められる、請求項１４から１５のいずれかに記載の相互作用エネルギーの算出装置。
【請求項１７】
前記算出部における前記相互作用エネルギーの各々の算出が、
前記標的分子と、前記複数の薬候補分子における一の薬候補分子との一の相互作用エネルギーを算出する際に、
前記標的分子と、前記一の薬候補分子との複合体における、前記標的分子と、前記一の薬候補分子とのペア相互作用エネルギーを算出する処理Ａと、
前記複合体における前記標的分子の、前記標的分子単体に対する分極相互作用エネルギーを算出する処理Ｂと、
前記複合体における前記一の薬候補分子の、前記一の薬候補分子単体に対する分極相互作用エネルギーを算出する処理Ｃと、
を含み、
更に、前記標的分子と、前記一の薬候補分子との前記一の相互作用エネルギーを算出する際に、前記計算対象水分子を考慮するときには、前記計算対象水分子と、前記一の薬候補分子とのペア相互作用エネルギーを算出する処理Ｄを含む、
請求項１３から１６のいずれかに記載の相互作用エネルギーの算出装置。
【請求項１８】
前記標的分子と、前記一の薬候補分子との前記相互作用エネルギーを算出する際に、前記計算対象水分子を考慮するときには、前記計算対象水分子が、前記処理Ａ、前記処理Ｂ、及び前記処理Ｃにおいて考慮される、請求項１７に記載の相互作用エネルギーの算出装置。
【請求項１９】
標的分子と、複数の薬候補分子の各々との相互作用エネルギーの少なくともいずれかに影響を与えるが、前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの全てには影響を与えず、
前記相互作用エネルギーの各々の算出において陽に扱う水分子として唯一考慮され、かつ、前記計算対象水分子が影響を与える相互作用エネルギーの算出においてのみ考慮され、
前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの各々を算出するのに用いられる、
計算対象水分子、
を抽出する抽出部を備えることを特徴とする相互作用エネルギーの算出装置。
【請求項２０】
標的分子と、複数の薬候補分子の各々との相互作用エネルギーの少なくともいずれかに影響を与えるが、前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの全てには影響を与えない、抽出された計算対象水分子のみを、前記相互作用エネルギーの各々の算出において陽に扱う水分子として考慮し、かつ、前記計算対象水分子が影響を与える相互作用エネルギーの算出においてのみ前記計算対象水分子を考慮して、前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの各々を算出する算出部を備えることを特徴とする相互作用エネルギーの算出装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本件は、標的分子と、複数の薬候補分子の各々との相互作用エネルギーの算出方法、及び算出装置、並びに前記算出方法を実行するプログラムに関する。
【背景技術】
【０００２】
あるタンパク質が体に悪影響を及ぼす機能部位を持つ場合、前記タンパク質をターゲットとする創薬では、前記タンパク質の前記機能部位と安定に結合するリガンドを設計することが必要である。前記リガンドが前記タンパク質に安定に結合することより、前記タンパク質の前記機能部位が塞がれる。その結果として、前記タンパク質による体への悪影響が抑制される。【０００３】
前記タンパク質と前記リガンドとの結合が安定かどうかの判定には、一般的に、前記タンパク質と前記リガンドとの相互作用エネルギーを量子力学（ＱＭ；ｑｕａｎｔｕｍ  ｍｅｃｈａｎｉｃｓ）に基づく手法で高精度に求める技術が用いられている（例えば、非特許文献１参照）。
前記ＱＭに基づく手法には、計算対象の原子の数が増えると、計算時間が指数関数的に増大してしまうという問題がある。タンパク質とリガンドとの複合体のみであれば、原子数は数千であっても、周囲に存在する水分子を含めると原子数は数万にもなり得る。そのため、水分子を含めたＱＭ計算を実施することには、計算時間の点で問題がある。【０００４】
そこで、前記相互作用エネルギー計算の効率化のために、水分子を一様な連続体としての溶媒で近似する技術が提案されている（例えば、非特許文献２参照）。この提案の技術では、原子数が増大しないため、計算時間の指数関数的な増大は回避できる。
しかし、この提案の技術では、前記タンパク質及び前記リガンドの複合体と水とが水素結合する場合に連続体近似が破綻してしまう。そのため、適切な相互作用エネルギーを求めるには、水を明示的に考慮する必要があり、計算時間の効率化の点では、十分ではない。【０００５】
複数の薬候補分子（リガンド）のうちから、タンパク質などの標的分子との相互作用エネルギーが大きい薬候補分子を予測することがある。その場合には、水分子の存在下で、前記標的分子と、前記複数の薬候補分子の各々との相互作用エネルギーの各々を求める。そして、求められた複数の相互作用エネルギーを対比する。前記複数の相互作用エネルギーを対比する際には、相互作用エネルギーの絶対値よりも、求められた前記複数の相互作用エネルギーの値の相対的な関係の信頼性が重要である。そのため、相対的な関係の信頼性が高い複数の相互作用エネルギーを、効率的に算出することが求められる。しかし、前述の提案の技術では、そのような要求に対して、十分であるとはいえない。
【先行技術文献】
【非特許文献】
【０００６】
【非特許文献１】Kａｚｕｏ  Ｋｉｔａｕｒａ，  Ｅｉｊｉ  Ｉｋｅｏ，  Ｔｏｓｈｉｏ  Ａｓａｄａ，  Ｔａｔｓｕｙａ  Ｎａｋａｎｏ，  Ｍａｓａｍｉ  Ｕｅｂａｙａｓｉ，  Ｃｈｅｍｉｃａｌ  Ｐｈｙｓｉｃｓ  Ｌｅｔｔｅｒｓ，  ３１３，  (１９９９)，  ７０１−７０６
【非特許文献２】Ｄ．  Ｇ．  Ｆｅｄｏｒｏｖ，  Ｋ．  Ｋｉｔａｕｒａ，  Ｈ．  Ｌｉ，  Ｊ．  Ｈ．  Ｊｅｎｓｅｎ，  Ｍ．  Ｓ．  Ｇｏｒｄｏｎ，  Ｊ．  Ｃｏｍｐｕｔ．  Ｃｈｅｍ．  ２７，  ９７６−９８５  （２００６）
【発明の概要】
【発明が解決しようとする課題】
【０００７】
本件は、従来における前記諸問題を解決し、以下の目的を達成することを課題とする。即ち、本件は、相対的な関係の信頼性が高い複数の相互作用エネルギーを、効率的に算出することができる相互作用エネルギーの算出方法、及び算出装置、並びに前記算出方法を実行するプログラムを提供することを目的とする。
【課題を解決するための手段】
【０００８】
前記課題を解決するための手段としては、以下の通りである。即ち、
開示の相互作用エネルギーの算出方法は、
計算機を用いて、標的分子と、複数の薬候補分子の各々との相互作用エネルギーの各々を算出する、相互作用エネルギーの算出方法であって、
前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの少なくともいずれかに影響を与えるが、前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの全てには影響を与えない計算対象水分子を抽出する工程と、
前記相互作用エネルギーの各々の算出において陽に扱う水分子として前記計算対象水分子のみを考慮し、かつ、前記計算対象水分子が影響を与える相互作用エネルギーの算出においてのみ前記計算対象水分子を考慮して、前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの各々を算出する工程と、
を含む。【０００９】
開示のプログラムは、コンピュータに、
標的分子と、複数の薬候補分子の各々との相互作用エネルギーの少なくともいずれかに影響を与えるが、前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの全てには影響を与えない計算対象水分子を抽出する工程と、
前記相互作用エネルギーの各々の算出において陽に扱う水分子として前記計算対象水分子のみを考慮し、かつ、前記計算対象水分子が影響を与える相互作用エネルギーの算出においてのみ前記計算対象水分子を考慮して、前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの各々を算出する工程と、
を実行させる。【００１０】
開示の算出装置は、
標的分子と、複数の薬候補分子の各々との相互作用エネルギーの少なくともいずれかに影響を与えるが、前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの全てには影響を与えない計算対象水分子を抽出する抽出部と、
前記相互作用エネルギーの各々の算出において陽に扱う水分子として前記計算対象水分子のみを考慮し、かつ、前記計算対象水分子が影響を与える相互作用エネルギーの算出においてのみ前記計算対象水分子を考慮して、前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの各々を算出する算出部と、
を備える。
【発明の効果】
【００１１】
開示の相互作用エネルギーの算出方法によると、従来における前記諸問題を解決し、前記目的を達成することができ、相対的な関係の信頼性が高い複数の相互作用エネルギーを、効率的に算出することができる。
開示のプログラムによると、従来における前記諸問題を解決し、前記目的を達成することができ、相対的な関係の信頼性が高い複数の相互作用エネルギーを、効率的に算出することができる。
開示の算出装置によると、従来における前記諸問題を解決し、前記目的を達成することができ、相対的な関係の信頼性が高い複数の相互作用エネルギーを、効率的に算出することができる。
【図面の簡単な説明】
【００１２】
【図１】図１は、開示の相互作用エネルギーの算出方法を説明するためのフローチャートである。
【図２】図２は、標的分子Ｐと薬候補分子Ｌとの複合体及び前記複合体を取り巻く水分子Ｗの模式図である。
【図３】図３は、第１の処理を説明するための模式図である。
【図４】図４は、第２の処理を説明するための模式図である。
【図５】図５は、第２の処理を説明するための模式図である。
【図６Ａ】図６Ａは、非重複水分子抽出処理を説明するための模式図である。
【図６Ｂ】図６Ｂは、非重複水分子抽出処理を説明するための模式図である。
【図７】図７は、計算対象水分子を抽出する工程の一例のフローチャートである。
【図８】図８は、相互作用エネルギーの各々を算出する工程の一例のフローチャートである。
【図９】図９は、開示の算出装置の構成例である。
【図１０】図１０は、開示の算出装置の他の構成例である。
【図１１】図１１は、開示の算出装置の他の構成例である。
【図１２】図１２は、実施例１の結果を示すグラフである。
【図１３】図１３は、比較例１の結果を示すグラフである。
【発明を実施するための形態】
【００１３】
（相互作用エネルギーの算出方法）
開示の相互作用エネルギーの算出方法は、計算機を用いて、標的分子と、複数の薬候補分子の各々との相互作用エネルギーの各々を算出する、相互作用エネルギーの算出方法である。
前記相互作用エネルギーの算出方法は、計算対象水分子を抽出する工程と、相互作用エネルギーの各々を算出する工程とを含み、更に必要に応じて、その他の工程を含む。
前記相互作用エネルギーの算出方法は、計算機を用いて行われる。前記相互作用エネルギーの算出方法に使用される前記計算機は、１つであってもよいし、複数であってもよい。例えば、複数の計算機に前記相互作用エネルギーの算出方法を分散させて実行させてもよい。
【００１４】
図１は、前記相互作用エネルギーの算出方法を説明するためのフローチャートである。
まず、前記計算対象水分子を抽出する工程を行う。この工程では、標的分子と、複数の薬候補分子の各々との相互作用エネルギーの少なくともいずれかに影響を与えるが、前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの全てには影響を与えない計算対象水分子を抽出する。
次に、前記相互作用エネルギーの各々を算出する工程を行う。この工程では、前記相互作用エネルギーの各々の算出において陽に扱う水分子として前記計算対象水分子のみを考慮し、かつ、前記計算対象水分子が影響を与える相互作用エネルギーの算出においてのみ前記計算対象水分子を考慮して、前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの各々を算出する。
【００１５】
標的分子と、薬候補分子との相互作用エネルギーの算出において、前記標的分子、及び前記薬候補分子とから形成される複合体を取り巻く水分子を陽に扱い、かつそれら水分子の全てを考慮すると、計算量が膨大になり、効率的な計算はできない。
一方、前記相互作用エネルギーの算出において、水分子を考慮しないと、計算の精度が大きく低下する。
【００１６】
開示の前記相互作用エネルギーの算出方法では、標的分子と、複数の薬候補分子の各々との相互作用エネルギーの少なくともいずれかに影響を与えるが、それらの相互作用エネルギーの全てには影響を与えない水分子を抽出し、その水分子のみを相互作用エネルギーの算出において、陽に扱う水分子として考慮する。そのような水分子は、相互作用エネルギーの算出において、考慮するか、考慮しないかにより、算出される複数の相互作用エネルギーの相対的な関係に変化を与える。そして、そのような水分子のみを相互作用エネルギーの算出において陽に扱う水分子として考慮すること、言い換えれば、そのような水分子以外の水分子を前記相互作用エネルギーの算出において陽に扱う水分子として考慮しないことにより、前記相互作用エネルギーの算出における計算量の増大を抑えることができる。また、そのような水分子を相互作用エネルギーの算出において陽に扱う水分子として考慮することで、複数の相互作用エネルギーの算出における相対的な関係の信頼性を高くできる。
【００１７】
＜計算対象水分子を抽出する工程＞
前記計算対象水分子を抽出する工程では、標的分子と、複数の薬候補分子の各々との相互作用エネルギーの少なくともいずれかに影響を与えるが、前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの全てには影響を与えない水分子を抽出する。本明細書、及び請求の範囲では、前記工程において抽出される前記水分子を、計算対象水分子と称する。
【００１８】
前記計算対象水分子を抽出する工程により抽出される前記計算対象水分子は、１つであってもよいし、複数であってもよい。
【００１９】
前記標的分子としては、特に制限はなく、目的に応じて適宜選択することができ、例えば、タンパク質、ＲＮＡ（リボ核酸、ｒｉｂｏｎｕｃｌｅｉｃ  ａｃｉｄ）、ＤＮＡ（デオキシリボ核酸、ｄｅｏｘｙｒｉｂｏｎｕｃｌｅｉｃ  ａｃｉｄ）などが挙げられる。
【００２０】
ここで、相互作用エネルギーに影響を与える水分子とは、前記相互作用エネルギーの算出において、その水分子を考慮する場合と、考慮しない場合とで、前記相互作用エネルギーの値が変わる水分子を意味する。
【００２１】
前記計算対象水分子を抽出する工程は、後述する第１の処理と、後述する第２の処理と、後述する非重複水分子抽出処理とを含むことが、算出される複数の相互作用エネルギーにおける相対的な関係の信頼性がより高くなる点で、好ましい。
【００２２】
なお、前記第１の処理、及び前記第２の処理は、前記標的分子と、前記複数の薬候補分子における一の薬候補分子とを含む一の相互作用算出対象空間における処理である。
前記非重複水分子抽出処理は、前記標的分子と、前記複数の薬候補分子の各々とを含む全ての相互作用算出対象空間のうちの少なくとも計算候補水分子が存在する複数の相互作用算出対象空間を対比して行う処理である。
ここで、前記相互作用算出対象空間とは、一の相互作用エネルギーを算出する際の、標的分子と、一の薬候補分子とを含む空間を意味する。
【００２３】
＜＜第１の処理＞＞
前記第１の処理では、前記標的分子、及び一の薬候補分子の双方と相互作用する水分子を抽出する。本明細書、及び請求の範囲では、前記第１の処理において抽出される前記水分子を、相互作用水分子と称する。
ここでの相互作用は、非結合性の相互作用を意味し、例えば、静電相互作用、ファンデルワールス相互作用などが挙げられる。
【００２４】
前記第１の処理で抽出される前記相互作用水分子は、通常、前記相互作用算出対象空間に、複数存在する。
【００２５】
前記第１の処理において前記相互作用水分子を抽出する方法としては、例えば、前記標的分子、及び前記一の薬候補分子の双方との距離がしきい値以下にある水分子を、前記相互作用水分子として抽出する方法などが挙げられる。
前記しきい値は、予め計算機内にデフォルト値として定められた値を用いてもよいし、使用者が任意の値を設定してもよい。
前記しきい値としては、例えば、６Åなどが挙げられる。即ち、前記相互作用水分子を抽出する方法としては、例えば、前記標的分子、及び前記一の薬候補分子の双方との距離が６Å以下にある水分子を、前記相互作用水分子として抽出する方法などが挙げられる。
前記標的分子と、前記水分子との距離とは、例えば、前記標的分子の表面と、前記水分子の酸素原子との距離である。
前記薬候補分子と、前記水分子との距離とは、例えば、前記薬候補分子の表面と、前記水分子の酸素原子との距離である。
【００２６】
＜＜第２の処理＞＞
前記第２の処理では、前記第１の処理で抽出された前記相互作用水分子のうちから、前記相互作用水分子の座標の揺らぎ幅、及び前記相互作用水分子の存在頻度の少なくともいずれかに基づいて、水分子を抽出する。本明細書、及び請求の範囲では、前記第２の処理において抽出される前記水分子を、計算候補水分子と称する。
前記第２の処理では、前記第１の処理で抽出された前記相互作用水分子のうちから、前記水分子の座標の揺らぎ幅、及び前記水分子の存在頻度に基づいて、前記計算候補水分子を抽出することが、算出される複数の相互作用エネルギーにおける相対的な関係の信頼性がより高くなる点で、好ましい。
【００２７】
前記計算候補水分子の座標の揺らぎ幅、及び前記計算候補水分子の存在頻度は、分子動力学シミュレーションにより求めることが好ましい。
【００２８】
前記計算候補水分子の座標の揺らぎ幅は、例えば、前記標的分子、前記一の薬候補分子、及び前記相互作用水分子について、分子動力学シミュレーションを行い、前記相互作用水分子の座標を解析し、根平均二乗揺らぎ（Ｒｏｏｔ  Ｍｅａｎ  Ｓｑｕａｒｅ  Ｆｌｕｃｔｕａｔｉｏｎ：ＲＭＳＦ）により評価することで求めることができる。ここで、前記ＲＭＳＦを求める式を以下に示す。
【化１】
2016103336.tif 000003
ここで、前記式中、Ｗ（ｔ）は、時刻ｔにおける相互作用水分子の座標を表し、ＷＲＥＦは、リファレンスとなる相互作用水分子の座標を表し、Ｔは、シミュレーション時間を表す。
前記ＲＭＳＦ（Ｗ）が小さい相互作用水分子は、標的分子及び一の薬候補分子の双方から近距離の座標に存在しつづけるため、前記標的分子と前記一の薬候補分子との相互作用への影響が大きい。
【００２９】
前記計算候補水分子の座標の揺らぎ幅に基づいて前記計算候補水分子を抽出する際には、しきい値を設けて抽出を行ってもよい。前記計算候補水分子の座標の揺らぎ幅に基づいて計算候補水分子を抽出する方法としては、例えば、前記ＲＭＳＦが２Å以下の相互作用水分子を抽出する方法などが挙げられる。前記しきい値は、予め計算機内にデフォルト値として定められた値を用いてもよいし、使用者が任意の値を設定してもよい。
【００３０】
前記相互作用水分子の存在頻度は、例えば、前記標的分子、前記一の薬候補分子、及び前記相互作用水分子について、分子動力学シミュレーションを行い、前記分子動力学シミュレーションの全時間における全ての相互作用水分子を１つの空間に配置することで得られる相互作用水分子の密度分布から求めることができる。相互作用水分子の存在頻度が高い座標には、特定の水分子ではないが、頻繁に相互作用水分子が存在する。そのため、相互作用水分子の存在頻度が高い座標には、仮想的な水分子が常に存在すると仮定できる。相互作用水分子の存在頻度が高い座標に存在する仮想的な水分子は、座標の揺らぎ幅が小さい相互作用水分子と同様に、前記標的分子及び前記一の薬候補分子の双方から近距離の座標に存在しつづけるため、前記標的分子と前記一の薬候補分子との相互作用への影響が大きい。
なお、上記のとおり、前記相互作用水分子の存在頻度に基づいて抽出される相互作用水分子は、特定の水分子ではなく、仮想的な水分子である。
【００３１】
相互作用水分子の存在頻度が高い座標の選択方法としては、例えば、相互作用水分子の密度分布における任意の低密度領域の密度よりも高い密度（例えば、１０倍以上の密度）の座標を選択する方法などが挙げられる。
【００３２】
前記分子動力学シミュレーションは、分子動力学計算プログラムを用いて行うことができる。前記分子動力学計算プログラムとしては、例えば、ＡＭＢＥＲ、ＣＨＡＲＭｍ、ＧＲＯＭＡＣＳ、ＧＲＯＭＯＳ、ＮＡＭＤ、ｍｙＰｒｅｓｔｏなどが挙げられる。
【００３３】
＜＜非重複水分子抽出処理＞＞
前記非重複水分子抽出処理では、前記標的分子と、前記複数の薬候補分子の各々とを含む全ての相互作用算出対象空間のうちの少なくとも前記計算候補水分子が存在する複数の相互作用算出対象空間を対比して、前記計算候補水分子のうちから、前記全ての相互作用算出対象空間において同じ座標にはない水分子を抽出する。
【００３４】
前記非重複水分子抽出処理は、前記標的分子と、前記複数の薬候補分子の各々とを含む全ての相互作用算出対象空間のうちの少なくとも前記計算候補水分子が存在する複数の相互作用算出対象空間を対比して、前記計算候補水分子のうちから、前記全ての相互作用算出対象空間において同じ座標にある水分子を排除する処理であってもよい。
【００３５】
なお、前記第２の処理で抽出された前記計算候補水分子のうちで、前記標的分子と、前記複数の薬候補分子の各々との複合体の全てにおいて同じ座標にある水分子は、標的分子と薬候補分子との相互作用エネルギー自体に影響を与えても、複数の相互作用エネルギーの相対的な関係には影響を与えない。そのため、そのような水分子を相互作用エネルギーの算出において考慮しても、計算量の増大になるだけである。
【００３６】
ここで、前記計算対象水分子を抽出する工程の一例である、前記第１の処理、前記第２の処理、及び前記非重複水分子抽出処理を、図２、図３、図４、図５、並びに図６Ａ、及び図６Ｂを用いて説明する。
図２は、標的分子Ｐと薬候補分子Ｌとの複合体及び前記複合体を取り巻く水分子Ｗの模式図である。前記複合体は、一の相互作用算出対象空間内に存在している。前記複合体の周囲には、多数の水分子Ｗが存在している。
まず、前記第１の処理では、図３に示すように、多数の水分子Ｗのうちから、標的分子Ｐ、及び薬候補分子Ｌの双方と相互作用する水分子Ｗを抽出する。通常、抽出される水分子Ｗは、標的分子Ｐ、及び薬候補分子Ｌの双方に近距離にある水分子Ｗである。
次に、前記第２の処理では、図４及び図５に示すように、前記第１の処理で抽出された水分子Ｗのうちから、水分子Ｗの座標の揺らぎ幅（図４）、及び水分子Ｗの存在頻度（図５）の少なくともいずれかに基づいて、更に水分子を抽出する。水分子Ｗの存在頻度に基づいて抽出される水分子は、仮想的な水分子Ｗ’である。この仮想的な水分子Ｗ’の座標には、頻繁に水分子Ｗが存在する。
次に、前記非重複水分子抽出処理では、全ての相互作用算出対象空間のうちの少なくとも計算候補水分子が存在する複数の相互作用算出対象空間を対比して、前記第２の処理で抽出された水分子のうちから、前記全ての相互作用算出対象空間において同じ座標にはない水分子Ｗを抽出する。２つの薬候補分子を開示の前記相互作用エネルギーの算出方法に用いる場合の、前記非重複水分子抽出処理の一例を以下に説明する。図６Ａに示すように、標的分子Ｐと一の薬候補分子Ｌ１とを含む相互作用算出対象空間において、前記第２の処理により、水分子Ｗ１、Ｗ２が抽出される。一方、図６Ｂに示すように、標的分子Ｐと他の薬候補分子Ｌ２とを含む相互作用算出対象空間において、前記第２の処理により、水分子Ｗ１が抽出される。この場合、水分子Ｗ１は、２つの相互作用算出対象空間において同じ座標に存在する。一方、水分子Ｗ２は、２つの相互作用算出対象空間において同じ座標に存在しない。したがって、前記非重複水分子抽出処理では、水分子Ｗ２が抽出される。
【００３７】
次に、前記計算対象水分子を抽出する工程の一例を説明する。
図７は、前記計算対象水分子を抽出する工程の一例のフローチャートである。
まず、一の薬候補分子を選択する。
次に、標的分子、及び前記一の薬候補分子の双方と相互作用する水分子を抽出する（第１の処理）。
次に、前記第１の処理で抽出された前記水分子のうちから、前記水分子の座標の揺らぎ幅、及び前記水分子の存在頻度の少なくともいずれかに基づいて、更に水分子を抽出する（第２の処理）。
次に、他の薬候補分子について、前記第１の処理、及び前記第２の処理を行うために、薬候補分子の変更を行う。
次に、前記他の薬候補分子について、前記第１の処理、及び前記第２の処理を行う。
前記標的分子との相互作用エネルギーを調べたい複数の薬候補分子について前記第１の処理、及び前記第２の処理が終われば、非重複水分子抽出処理に進む。
前記非重複水分子抽出処理では、複数の相互作用算出対象空間を対比して、前記第２の処理で抽出された前記水分子のうちから、前記標的分子と、前記複数の薬候補分子の各々との複合体の全てにおいて同じ座標にはない水分子を抽出する。
以上により、前記計算対象水分子を抽出する工程の一例が終了する。
【００３８】
＜相互作用エネルギーの各々を算出する工程＞
前記相互作用エネルギーの各々を算出する工程では、前記相互作用エネルギーの各々の算出において陽に扱う水分子として前記計算対象水分子のみを考慮し、かつ、前記計算対象水分子が影響を与える相互作用エネルギーの算出においてのみ前記計算対象水分子を考慮して、前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの各々を算出する。
【００３９】
ここで、陽に扱うとは、分子として明示的に考慮することを意味する。なお、陽に扱わない方法としては、連続体近似によって陰に扱う行う方法などが挙げられる。
【００４０】
前記相互作用エネルギーの各々を算出する工程は、前記標的分子と、前記複数の薬候補分子における一の薬候補分子との一の相互作用エネルギーを算出する際に、後述する処理Ａと、後述する処理Ｂと、後述する処理Ｃとを含み、更に、前記標的分子と、前記一の薬候補分子との前記一の相互作用エネルギーを算出する際に、前記薬候補水分子を考慮するときには、後述する処理Ｄとを含むことが、算出される複数の相互作用エネルギーにおける相対的な関係の信頼性がより高くなる点で、好ましい。
なお、これらの処理の順序は問わない。
【００４１】
＜＜処理Ａ＞＞
前記処理Ａでは、前記標的分子と、前記一の薬候補分子との複合体における、前記標的分子と、前記一の薬候補分子とのペア相互作用エネルギーを算出する。
【００４２】
＜＜処理Ｂ＞＞
前記処理Ｂでは、前記複合体における前記標的分子の、前記標的分子単体に対する分極相互作用エネルギーを算出する。
【００４３】
前記分極相互作用エネルギーとは、前記標的分子単体の状態から、前記標的分子が前記一の薬候補分子とともに前記複合体を形成した際に生じる、前記標的分子内の分極のエネルギーである。
【００４４】
前記分極相互作用エネルギーは、例えば、以下の方法で求めることができる。
前記標的分子単体の状態のエネルギー（ＥＰ０）を算出する。この際、前記標的分子の周囲には、仮想的な水分子を配置する。仮想的な水分子の配置は、例えば、連続体近似により行うことができる。
前記複合体における前記標的分子のエネルギー（ＥＰ１）を算出する。この際、前記複合体の周囲には、仮想的な水分子を配置する。仮想的な水分子の配置は、例えば、連続体近似により行うことができる。
そして、前記エネルギー（ＥＰ０）と前記エネルギー（ＥＰ１）との差から、前記分極相互作用エネルギーが求められる。
【００４５】
＜＜処理Ｃ＞＞
前記処理Ｃでは、前記複合体における前記一の薬候補分子の、前記一の薬候補分子単体に対する分極相互作用エネルギーを算出する。
【００４６】
前記分極相互作用エネルギーとは、前記一の薬候補分子単体の状態から、前記一の薬候補分子が前記標的分子とともに前記複合体を形成した際に生じる、前記一の薬候補分子内の分極のエネルギーである。
【００４７】
前記分極相互作用エネルギーは、例えば、以下の方法で求めることができる。
前記一の薬候補分子が単体の状態のエネルギー（ＥＬ０）を算出する。この際、前記一の薬候補分子の周囲には、仮想的な水分子を配置する。仮想的な水分子の配置は、例えば、連続体近似により行うことができる。
前記複合体における前記一の薬候補分子のエネルギー（ＥＬ１）を算出する。この際、前記複合体の周囲には、仮想的な水分子を配置する。仮想的な水分子の配置は、例えば、連続体近似により行うことができる。
そして、前記エネルギー（ＥＬ０）と前記エネルギー（ＥＬ１）との差から、前記分極相互作用エネルギーが求められる。
【００４８】
＜＜処理Ｄ＞＞
前記処理Ｄでは、前記計算対象水分子水分子と、前記一の薬候補分子とのペア相互作用エネルギーを算出する。
前記処理Ｄは、前記標的分子と、前記一の薬候補分子との前記一の相互作用エネルギーを算出する際に、前記計算対象水分子を考慮するときに行われる。
【００４９】
前記標的分子と、前記一の薬候補分子とを含む一の相互作用算出対象空間において、前記計算対象水分子が存在しない場合には、前記標的分子と、前記一の薬候補分子との前記一の相互作用エネルギーを算出する際に、前記処理Ｄは行われない。
【００５０】
前記処理Ｄにおいては、前記計算対象水分子は、前記標的分子の一部とみなされることが、算出される前記相互作用エネルギーの精度が向上する点で好ましい。
【００５１】
前記標的分子と、前記一の薬候補分子との前記一の相互作用エネルギーを算出する際に、前記計算対象水分子を考慮するときには、前記計算対象水分子は、前記処理Ａ、前記処理Ｂ、及び前記処理Ｃにおいて考慮されることが、算出される前記相互作用エネルギーの精度が向上する点で好ましい。
前記計算対象水分子を考慮する際には、前記計算対象水分子を前記複合体の一部とみなすことが好ましい。
【００５２】
前記相互作用エネルギーの各々を算出する工程は、例えば、量子力学（ＱＭ）計算により求められる。前記量子力学計算としては、例えば、分子軌道法による分子軌道計算が挙げられる。前記分子軌道計算としては、例えば、非経験的分子軌道計算（ａｂ  ｉｎｉｔｉｏ分子軌道計算）、半経験的分子軌道計算などが挙げられる。
前記非経験的分子軌道計算の方法論としては、例えば、ハートリー−フォック法、電子相関法などが挙げられる。
前記半経験的分子軌道計算の方法論としては、例えば、ＣＮＤＯ、ＩＮＤＯ、ＡＭ１、ＰＭ３などが挙げられる。
前記非経験的分子軌道計算のプログラムとしては、例えば、Ｇａｕｓｓｉａｎ０３、ＧＡＭＥＳＳ、ＡＢＩＮＩＴ−ＭＰ、Ｐｒｏｔｅｉｎ  ＤＦなどが挙げられる。
前記半経験的分子軌道計算のプログラムとしては、例えば、ＭＯＰＡＣなどが挙げられる。
【００５３】
次に、前記相互作用エネルギーの各々を算出する工程の一例を説明する。
図８は、前記相互作用エネルギーの各々を算出する工程の一例のフローチャートである。
まず、前記標的分子と、一の薬候補分子との複合体における、前記標的分子と、前記一の薬候補分子とのペア相互作用エネルギー（ＰＩＥ）を算出する（処理Ａ）。
次に、前記標的分子の分極相互作用エネルギー（分極ＩＥ）を算出する（処理Ｂ）。
次に、前記一の薬候補分子の分極相互作用エネルギー（分極ＩＥ）を算出する（処理Ｃ）。
そして、前記標的分子と、前記一の薬候補分子とを含む一の相互作用算出対象空間に前記計算対象水分子が存在する場合には、前記計算対象水分子と前記一の薬候補分子とのペア相互作用エネルギー（ＰＩＥ）を算出する（処理Ｄ）。一方、前記一の相互作用算出対象空間に前記計算対象水分子が存在しない場合には、前記処理Ｄは行わない。
次に、他の薬候補分子について、前記処理Ａ、前記処理Ｂ、及び前記処理Ｃ、並びに、前記標的分子と、前記他の薬候補分子とを含む一の相互作用算出対象空間に前記計算対象水分子が存在する場合には、前記処理Ｄを行うために、薬候補分子の変更を行う。
次に、前記他の薬候補分子について、前記処理Ａ、前記処理Ｂ、及び前記処理Ｃ、並びに、前記標的分子と、前記他の薬候補分子とを含む一の相互作用算出対象空間に前記計算対象水分子が存在する場合には、前記第処理Ｄを行う。
そして、前記標的分子との相互作用エネルギーを調べたい複数の薬候補分子について、前記処理Ａ、前記処理Ｂ、前記処理Ｃ、及び前述の条件に応じて前記処理Ｄが終われば、本工程が終了する。
【００５４】
前記相互作用エネルギーの算出方法は、例えば、ＣＰＵ（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）、ＲＡＭ（Ｒａｎｄｏｍ  Ａｃｃｅｓｓ  Ｍｅｍｏｒｙ）、ハードディスク、各種周辺機器等を備えた通常のコンピュータシステム（例えば、各種ネットワークサーバ、ワークステーション、パーソナルコンピュータ等）を用いることによって実現することができる。
【００５５】
（プログラム）
開示のプログラムは、コンピュータに、開示の前記相互作用エネルギーの算出方法を実行させるプログラムである。
前記相互作用エネルギーの算出方法の実行における好適な態様は、前記相互作用エネルギーの算出方法における好適な態様と同じである。
【００５６】
前記プログラムは、使用するコンピュータシステムの構成及びオペレーティングシステムの種類・バージョンなどに応じて、公知の各種のプログラム言語を用いて作成することができる。
【００５７】
前記プログラムは、内蔵ハードディスク、外付けハードディスクなどの記録媒体に記録しておいてもよいし、ＣＤ−ＲＯＭ（Ｃｏｍｐａｃｔ  Ｄｉｓｃ  Ｒｅａｄ  Ｏｎｌｙ  Ｍｅｍｏｒｙ）、ＤＶＤ−ＲＯＭ（Ｄｉｇｉｔａｌ  Ｖｅｒｓａｔｉｌｅ  Ｄｉｓｋ  Ｒｅａｄ  Ｏｎｌｙ  Ｍｅｍｏｒｙ）、ＭＯディスク（Ｍａｇｎｅｔｏ−Ｏｐｔｉｃａｌ  ｄｉｓｋ）、ＵＳＢメモリ〔ＵＳＢ（Ｕｎｉｖｅｒｓａｌ  Ｓｅｒｉａｌ  Ｂｕｓ）  ｆｌａｓｈ  ｄｒｉｖｅ〕などの記録媒体に記録しておいてもよい。前記プログラムをＣＤ−ＲＯＭ、ＤＶＤ−ＲＯＭ、ＭＯディスク、ＵＳＢメモリなどの記録媒体に記録する場合には、必要に応じて随時、コンピュータシステムが有する記録媒体読取装置を通じて、これを直接、又はハードディスクにインストールして使用することができる。また、コンピュータシステムから情報通信ネットワークを通じてアクセス可能な外部記憶領域（他のコンピュータ等）に前記プログラムを記録しておき、必要に応じて随時、前記外部記憶領域から情報通信ネットワークを通じてこれを直接、又はハードディスクにインストールして使用することもできる。
前記プログラムは、複数の記録媒体に、任意の処理毎に分割されて記録されていてもよい。
【００５８】
（コンピュータが読み取り可能な記録媒体）
開示のコンピュータが読み取り可能な記録媒体は、開示の前記プログラムを記録してなる。
前記コンピュータが読み取り可能な記録媒体としては、特に制限はなく、目的に応じて適宜選択することができ、例えば、内蔵ハードディスク、外付けハードディスク、ＣＤ−ＲＯＭ、ＤＶＤ−ＲＯＭ、ＭＯディスク、ＵＳＢメモリなどが挙げられる。
前記記録媒体は、前記プログラムが任意の処理毎に分割されて記録された複数の記録媒体であってもよい。
【００５９】
（相互作用エネルギーの算出装置）
開示の相互作用エネルギーの算出装置の一態様は、抽出部と、算出部とを少なくとも備え、更に必要に応じて、その他の部を備える。
前記抽出部は、標的分子と、複数の薬候補分子の各々との相互作用エネルギーの少なくともいずれかに影響を与えるが、前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの全てには影響を与えない計算対象水分子を抽出する。
前記算出部は、前記相互作用エネルギーの各々の算出において陽に扱う水分子として前記計算対象水分子のみを考慮し、かつ、前記計算対象水分子が影響を与える相互作用エネルギーの算出においてのみ前記計算対象水分子を考慮して、前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの各々を算出する。
【００６０】
開示の相互作用エネルギーの算出装置の他の一態様は、抽出部を備え、更に必要に応じて、その他の部を備える。
この態様における前記抽出部は、計算対象水分子を抽出する。
この態様における前記計算対象水分子は、標的分子と、複数の薬候補分子の各々との相互作用エネルギーの少なくともいずれかに影響を与えるが、前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの全てには影響を与えない。
この態様における前記計算対象水分子は、前記相互作用エネルギーの各々の算出において陽に扱う水分子として唯一考慮され、かつ、前記計算対象水分子が影響を与える相互作用エネルギーの算出においてのみ考慮される。
この態様における前記計算対象水分子は、前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの各々を算出するのに用いられる。
【００６１】
開示の相互作用エネルギーの算出装置の他の一態様は、算出部を備え、更に必要に応じて、その他の部を備える。
この態様における前記算出部は、標的分子と、複数の薬候補分子の各々との相互作用エネルギーの少なくともいずれかに影響を与えるが、前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの全てには影響を与えない、抽出された計算対象水分子のみを、前記相互作用エネルギーの各々の算出において陽に扱う水分子として考慮し、かつ、前記計算対象水分子が影響を与える相互作用エネルギーの算出においてのみ前記計算対象水分子を考慮して、前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの各々を算出する。
【００６２】
前記相互作用エネルギーの算出装置における各部の処理方法の好適な態様は、前記相互作用エネルギーの算出方法における各工程の好適な態様と同じである。
前記算出装置は、前記プログラムが任意の処理毎に分割されて記録された複数の記録媒体をそれぞれに備える複数の算出装置であってもよい。
【００６３】
図９に、開示の算出装置の構成例を示す。
算出装置１０は、例えば、ＣＰＵ１１、メモリ１２、記憶部１３、表示部１４、入力部１５、出力部１６、Ｉ／Ｏインターフェース部１７等がシステムバス１８を介して接続されて構成される。
【００６４】
ＣＰＵ（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）１１は、演算（四則演算、比較演算等）、ハードウエア及びソフトウエアの動作制御などを行う。
【００６５】
メモリ１２は、ＲＡＭ（Ｒａｎｄｏｍ  Ａｃｃｅｓｓ  Ｍｅｍｏｒｙ）、ＲＯＭ（Ｒｅａｄ  Ｏｎｌｙ  Ｍｅｍｏｒｙ）などのメモリである。前記ＲＡＭは、前記ＲＯＭ及び記憶部１３から読み出されたＯＳ（Ｏｐｅｒａｔｉｎｇ  Ｓｙｓｔｅｍ）及びアプリケーションプログラムなどを記憶し、ＣＰＵ１１の主メモリ及びワークエリアとして機能する。
【００６６】
記憶部１３は、各種プログラム及びデータを記憶する装置であり、例えば、ハードディスクである。記憶部１３には、ＣＰＵ１１が実行するプログラム、プログラム実行に必要なデータ、ＯＳなどが格納される。
前記プログラムは、記憶部１３に格納され、メモリ１２のＲＡＭ（主メモリ）にロードされ、ＣＰＵ１１により実行される。
【００６７】
表示部１４は、表示装置であり、例えば、ＣＲＴモニタ、液晶パネル等のディスプレイ装置である。
入力部１５は、各種データの入力装置であり、例えば、キーボード、ポインティングデバイス（例えば、マウス等）などである。
出力部１６は、各種データの出力装置であり、例えば、プリンタである。
Ｉ／Ｏインターフェース部１７は、各種の外部装置を接続するためのインターフェースである。例えば、ＣＤ−ＲＯＭ、ＤＶＤ−ＲＯＭ、ＭＯディスク、ＵＳＢメモリなどのデータの入出力を可能にする。
【００６８】
図１０に、開示の算出装置の他の構成例を示す。
図１０の構成例は、クラウド型の構成例であり、ＣＰＵ１１が、記憶部１３等とは独立している。この構成例では、ネットワークインターフェース部１９、２０を介して、記憶部１３等を格納するコンピュータ３０と、ＣＰＵ１１を格納するコンピュータ４０とが接続される。
ネットワークインターフェース部１９、２０は、インターネットを利用して、通信を行うハードウェアである。
【００６９】
図１１に、開示の算出装置の他の構成例を示す。
図１１の構成例は、クラウド型の構成例であり、記憶部１３が、ＣＰＵ１１等とは独立している。この構成例では、ネットワークインターフェース部１９、２０を介して、ＣＰＵ１１等を格納するコンピュータ３０と、記憶部１３を格納するコンピュータ４０とが接続される。
【実施例】
【００７０】
以下、開示の技術について説明するが、開示の技術は下記実施例に何ら限定されるものではない。【００７１】
（実施例１）
開示の相互作用エネルギーの算出方法を、図１、図７、及び図８に示すフローに従って、実施した。
標的分子としては、タンパク質である血液凝固因子Ｘａを用いた。
薬候補分子としては、４種類のリガンド（ＲＲＰ，  ＲＴＲ，  ＲＤＲ，  ＲＲＲ）を用いた。
タンパク質である血液凝固因子Ｘａと、４種類のリガンド（ＲＲＰ，  ＲＴＲ，  ＲＤＲ，  ＲＲＲ）とは、それぞれの複合体構造が、ＰＤＢ結晶構造（それぞれ１ＮＦＵ，  １ＮＦＹ，  １ＮＦＸ，  １ＮＦＷ）として知られている。【００７２】
水分子の座標の揺らぎ幅、及び水分子の存在頻度を求める際には、前記ＰＤＢ結晶構造を初期値とする分子動力学シミュレーションを１０ｎｓ実施した。
なお、前記第１の処理においては、距離のしきい値を３Åとし、前記標的分子、及び前記薬候補分子との距離が３Å以内の水分子を抽出した。
また、前記第２の処理の、水分子の座標の揺らぎ幅に基づく水分子の抽出においては、ＲＭＳＦの値のしきい値を２Åとし、ＲＭＳＦの値が２Å以内の水分子を抽出した。
また、前記第２の処理の、水分子の存在頻度は、前記分子動力学シミュレーションの終わりの２ｎｓにおける全ての水分子を１つの空間に配置することで得られる水分子の密度分布から求めた。そして、水分子の密度分布における任意の低密度領域の密度より密度が１０倍高い座標を選択し、前記座標に仮想的な水分子が存在するものとした。【００７３】
図７のフローに従って水分子の抽出を行ったところ、ＲＲＲ近傍にのみ存在する水分子が２つ抽出できた。【００７４】
次に、図８のフローに従って、血液凝固因子Ｘａと、４種類のリガンド（ＲＲＰ，  ＲＴＲ，  ＲＤＲ，  ＲＲＲ）の各々との相互作用エネルギーを算出した。ここで、血液凝固因子ＸａとリガンドＲＲＲとの相互作用エネルギーの算出においてのみ、抽出された２つの水分子を考慮した。【００７５】
得られた４つの相互作用エネルギーと、血液凝固因子Ｘａと４種類のリガンド（ＲＲＰ，  ＲＴＲ，  ＲＤＲ，  ＲＲＲ）の各々との相互作用エネルギーの実測値とを対比した。結果を図１２に示した。横軸の数値が実測値であり、縦軸の数値が本方法で算出した相互作用エネルギー値（算出値）である。実測値と、算出値との間で線形の相関関係が得られた。【００７６】
（比較例１）
実施例１において、水分子の抽出を行わなかった以外は、実施例１と同様にして、血液凝固因子Ｘａと、４種類のリガンド（ＲＲＰ，  ＲＴＲ，  ＲＤＲ，  ＲＲＲ）の各々との相互作用エネルギーを算出した。【００７７】
得られた４つの相互作用エネルギーと、血液凝固因子Ｘａと４種類のリガンド（ＲＲＰ，  ＲＴＲ，  ＲＤＲ，  ＲＲＲ）の各々との相互作用エネルギーの実測値とを対比した。結果を図１３に示した。横軸の数値が実測値であり、縦軸の数値が本方法で算出した相互作用エネルギー値（算出値）である。実測値と、算出値との間で線形の相関関係が得られなかった。【００７８】
実施例１は、比較例１と比べて、相互作用エネルギーの計算時間の増大はほとんどなかった。
【符号の説明】
【００７９】
１０    算出装置
１１    ＣＰＵ
１２    メモリ
１３    記憶部
１４    表示部
１５    入力部
１６    出力部
１７    Ｉ／Ｏインターフェース部
１８    システムバス
１９    ネットワークインターフェース部
２０    ネットワークインターフェース部
３０    コンピュータ
４０    コンピュータ
Ｌ      薬候補分子
Ｌ１    一の薬候補分子
Ｌ２    他の薬候補分子
Ｐ      標的分子
Ｗ      水分子
Ｗ’    仮想的な水分子
【図１】
2016103336.tif 000004
【図２】
2016103336.tif 000005
【図３】
2016103336.tif 000006
【図４】
2016103336.tif 000007
【図５】
2016103336.tif 000008
【図６Ａ】
2016103336.tif 000009
【図６Ｂ】
2016103336.tif 000010
【図７】
2016103336.tif 000011
【図８】
2016103336.tif 000012
【図９】
2016103336.tif 000013
【図１０】
2016103336.tif 000014
【図１１】
2016103336.tif 000015
【図１２】
2016103336.tif 000016
【図１３】
2016103336.tif 000017
【手続補正書】
【提出日】平成29年6月8日(2017.6.8)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
計算機を用いて、標的分子と、複数の薬候補分子の各々との相互作用エネルギーの各々を算出する、相互作用エネルギーの算出方法であって、
前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの少なくともいずれかに影響を与えるが、前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの全てには影響を与えない計算対象水分子を抽出する工程と、
前記相互作用エネルギーの各々の算出において陽に扱う水分子として前記計算対象水分子のみを考慮し、かつ、前記計算対象水分子が影響を与える相互作用エネルギーの算出においてのみ前記計算対象水分子を考慮して、前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの各々を算出する工程と、
を含むことを特徴とする相互作用エネルギーの算出方法。
【請求項２】
前記計算対象水分子を抽出する工程が、
前記標的分子と、前記複数の薬候補分子における一の薬候補分子とを含む一の相互作用算出対象空間における、
前記標的分子、及び前記一の薬候補分子の双方と相互作用する相互作用水分子を抽出する第１の処理と、
前記相互作用水分子のうちから、前記相互作用水分子の座標の揺らぎ幅、及び前記相互作用水分子の存在頻度の少なくともいずれかに基づいて、計算候補水分子を抽出する第２の処理と、
を含み、
更に、前記計算対象水分子を抽出する工程が、
前記標的分子と、前記複数の薬候補分子の各々とを含む全ての相互作用算出対象空間のうちの少なくとも前記計算候補水分子が存在する複数の相互作用算出対象空間を対比して、前記計算候補水分子のうちから、前記全ての相互作用算出対象空間において同じ座標にはない水分子を抽出する処理、
を含む請求項１に記載の相互作用エネルギーの算出方法。
【請求項３】
前記第２の処理が、前記相互作用水分子のうちから、前記相互作用水分子の座標の揺らぎ幅、及び前記相互作用水分子の存在頻度に基づいて、前記計算候補水分子を抽出する処理である、請求項２に記載の相互作用エネルギーの算出方法。
【請求項４】
前記相互作用水分子の座標の揺らぎ幅、及び前記相互作用水分子の存在頻度が、分子動力学シミュレーションにより求められる、請求項２から３のいずれかに記載の相互作用エネルギーの算出方法。
【請求項５】
前記相互作用エネルギーの各々を算出する工程が、
前記標的分子と、前記複数の薬候補分子における一の薬候補分子との一の相互作用エネルギーを算出する際に、
前記標的分子と、前記一の薬候補分子との複合体における、前記標的分子と、前記一の薬候補分子とのペア相互作用エネルギーを算出する処理Ａと、
前記複合体における前記標的分子の、前記標的分子単体に対する分極相互作用エネルギーを算出する処理Ｂと、
前記複合体における前記一の薬候補分子の、前記一の薬候補分子単体に対する分極相互作用エネルギーを算出する処理Ｃと、
を含み、
更に、前記標的分子と、前記一の薬候補分子との前記一の相互作用エネルギーを算出する際に、前記計算対象水分子を考慮するときには、前記計算対象水分子と、前記一の薬候補分子とのペア相互作用エネルギーを算出する処理Ｄを含む、
請求項１から４のいずれかに記載の相互作用エネルギーの算出方法。
【請求項６】
前記標的分子と、前記一の薬候補分子との前記相互作用エネルギーを算出する際に、前記計算対象水分子を考慮するときには、前記計算対象水分子が、前記処理Ａ、前記処理Ｂ、及び前記処理Ｃにおいて考慮される、請求項５に記載の相互作用エネルギーの算出方法。
【請求項７】
コンピュータに、
標的分子と、複数の薬候補分子の各々との相互作用エネルギーの少なくともいずれかに影響を与えるが、前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの全てには影響を与えない計算対象水分子を抽出する工程と、
前記相互作用エネルギーの各々の算出において陽に扱う水分子として前記計算対象水分子のみを考慮し、かつ、前記計算対象水分子が影響を与える相互作用エネルギーの算出においてのみ前記計算対象水分子を考慮して、前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの各々を算出する工程と、
を実行させることを特徴とするプログラム。
【請求項８】
標的分子と、複数の薬候補分子の各々との相互作用エネルギーの少なくともいずれかに影響を与えるが、前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの全てには影響を与えない計算対象水分子を抽出する抽出部と、
前記相互作用エネルギーの各々の算出において陽に扱う水分子として前記計算対象水分子のみを考慮し、かつ、前記計算対象水分子が影響を与える相互作用エネルギーの算出においてのみ前記計算対象水分子を考慮して、前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの各々を算出する算出部と、
を備えることを特徴とする相互作用エネルギーの算出装置。
【請求項９】
標的分子と、複数の薬候補分子の各々との相互作用エネルギーの少なくともいずれかに影響を与えるが、前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの全てには影響を与えず、
前記相互作用エネルギーの各々の算出において陽に扱う水分子として唯一考慮され、かつ、前記計算対象水分子が影響を与える相互作用エネルギーの算出においてのみ考慮され、
前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの各々を算出するのに用いられる、
計算対象水分子、
を抽出する抽出部を備えることを特徴とする相互作用エネルギーの算出装置。
【請求項１０】
標的分子と、複数の薬候補分子の各々との相互作用エネルギーの少なくともいずれかに影響を与えるが、前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの全てには影響を与えない、抽出された計算対象水分子のみを、前記相互作用エネルギーの各々の算出において陽に扱う水分子として考慮し、かつ、前記計算対象水分子が影響を与える相互作用エネルギーの算出においてのみ前記計算対象水分子を考慮して、前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの各々を算出する算出部を備えることを特徴とする相互作用エネルギーの算出装置。
【国際調査報告】
2016103336.tif 000018
2016103336.tif 000019
2016103336.tif 000020
2016103336.tif 000021