JP6409880B2 - 相互作用エネルギーの算出方法、及び算出装置、並びにプログラム - Google Patents

相互作用エネルギーの算出方法、及び算出装置、並びにプログラム Download PDF

Info

Publication number
JP6409880B2
JP6409880B2 JP2016565634A JP2016565634A JP6409880B2 JP 6409880 B2 JP6409880 B2 JP 6409880B2 JP 2016565634 A JP2016565634 A JP 2016565634A JP 2016565634 A JP2016565634 A JP 2016565634A JP 6409880 B2 JP6409880 B2 JP 6409880B2
Authority
JP
Japan
Prior art keywords
molecule
calculation
target
drug candidate
interaction
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2016565634A
Other languages
English (en)
Other versions
JPWO2016103336A1 (ja
Inventor
佐藤　博之
博之 佐藤
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fujitsu Ltd
Original Assignee
Fujitsu Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Fujitsu Ltd filed Critical Fujitsu Ltd
Publication of JPWO2016103336A1 publication Critical patent/JPWO2016103336A1/ja
Application granted granted Critical
Publication of JP6409880B2 publication Critical patent/JP6409880B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B15/00—ICT specially adapted for analysing two-dimensional or three-dimensional molecular structures, e.g. structural or functional relations or structure alignment
    • G16B15/30—Drug targeting using structural data; Docking or binding prediction
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C10/00—Computational theoretical chemistry, i.e. ICT specially adapted for theoretical aspects of quantum chemistry, molecular mechanics, molecular dynamics or the like
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B15/00—ICT specially adapted for analysing two-dimensional or three-dimensional molecular structures, e.g. structural or functional relations or structure alignment
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B5/00—ICT specially adapted for modelling or simulations in systems biology, e.g. gene-regulatory networks, protein interaction networks or metabolic networks

Description

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

Claims (8)

  1. 計算機を用いて、標的分子と、複数の薬候補分子の各々との相互作用エネルギーの各々を算出する、相互作用エネルギーの算出方法であって、
    前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの少なくともいずれかに影響を与えるが、前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの全てには影響を与えない計算対象水分子を抽出する工程と、
    前記相互作用エネルギーの各々の算出において陽に扱う水分子として前記計算対象水分子のみを考慮し、かつ、前記計算対象水分子が影響を与える相互作用エネルギーの算出においてのみ前記計算対象水分子を考慮して、前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの各々を算出する工程と、
    を含み、
    前記計算対象水分子を抽出する工程が、
    前記標的分子と、前記複数の薬候補分子における一の薬候補分子とを含む一の相互作用算出対象空間における、
    前記標的分子、及び前記一の薬候補分子の双方と相互作用する相互作用水分子を抽出する第１の処理と、
    前記相互作用水分子のうちから、前記相互作用水分子の座標の揺らぎ幅、及び前記相互作用水分子の存在頻度の少なくともいずれかに基づいて、計算候補水分子を抽出する第２の処理と、
    を含み、
    更に、前記計算対象水分子を抽出する工程が、
    前記標的分子と、前記複数の薬候補分子の各々とを含む全ての相互作用算出対象空間のうちの少なくとも前記計算候補水分子が存在する複数の相互作用算出対象空間を対比して、前記計算候補水分子のうちから、前記全ての相互作用算出対象空間において同じ座標にはない水分子を抽出する処理、
    を含み、
    前記相互作用エネルギーの各々を算出する工程が、
    前記標的分子と、前記複数の薬候補分子における一の薬候補分子との一の相互作用エネルギーを算出する際に、
    前記標的分子と、前記一の薬候補分子との複合体における、前記標的分子と、前記一の薬候補分子とのペア相互作用エネルギーを算出する処理Ａと、
    前記複合体における前記標的分子の、前記標的分子単体に対する分極相互作用エネルギーを算出する処理Ｂと、
    前記複合体における前記一の薬候補分子の、前記一の薬候補分子単体に対する分極相互作用エネルギーを算出する処理Ｃと、
    を含み、
    更に、前記標的分子と、前記一の薬候補分子との前記一の相互作用エネルギーを算出する際に、前記計算対象水分子を考慮するときには、前記計算対象水分子と、前記一の薬候補分子とのペア相互作用エネルギーを算出する処理Ｄを含む、
    ことを特徴とする相互作用エネルギーの算出方法。
  2. 前記第２の処理が、前記相互作用水分子のうちから、前記相互作用水分子の座標の揺らぎ幅、及び前記相互作用水分子の存在頻度に基づいて、前記計算候補水分子を抽出する処理である、請求項１に記載の相互作用エネルギーの算出方法。
  3. 前記相互作用水分子の座標の揺らぎ幅、及び前記相互作用水分子の存在頻度が、分子動力学シミュレーションにより求められる、請求項１から２のいずれかに記載の相互作用エネルギーの算出方法。
  4. 前記標的分子と、前記一の薬候補分子との前記相互作用エネルギーを算出する際に、前記計算対象水分子を考慮するときには、前記計算対象水分子が、前記処理Ａ、前記処理Ｂ、及び前記処理Ｃにおいて考慮される、請求項１から３のいずれかに記載の相互作用エネルギーの算出方法。
  5. コンピュータに、
    標的分子と、複数の薬候補分子の各々との相互作用エネルギーの少なくともいずれかに影響を与えるが、前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの全てには影響を与えない計算対象水分子を抽出する工程と、
    前記相互作用エネルギーの各々の算出において陽に扱う水分子として前記計算対象水分子のみを考慮し、かつ、前記計算対象水分子が影響を与える相互作用エネルギーの算出においてのみ前記計算対象水分子を考慮して、前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの各々を算出する工程と、
    を実行させるプログラムであって、
    前記計算対象水分子を抽出する工程が、
    前記標的分子と、前記複数の薬候補分子における一の薬候補分子とを含む一の相互作用算出対象空間における、
    前記標的分子、及び前記一の薬候補分子の双方と相互作用する相互作用水分子を抽出する第１の処理と、
    前記相互作用水分子のうちから、前記相互作用水分子の座標の揺らぎ幅、及び前記相互作用水分子の存在頻度の少なくともいずれかに基づいて、計算候補水分子を抽出する第２の処理と、
    を含み、
    更に、前記計算対象水分子を抽出する工程が、
    前記標的分子と、前記複数の薬候補分子の各々とを含む全ての相互作用算出対象空間のうちの少なくとも前記計算候補水分子が存在する複数の相互作用算出対象空間を対比して、前記計算候補水分子のうちから、前記全ての相互作用算出対象空間において同じ座標にはない水分子を抽出する処理、
    を含み、
    前記相互作用エネルギーの各々を算出する工程が、
    前記標的分子と、前記複数の薬候補分子における一の薬候補分子との一の相互作用エネルギーを算出する際に、
    前記標的分子と、前記一の薬候補分子との複合体における、前記標的分子と、前記一の薬候補分子とのペア相互作用エネルギーを算出する処理Ａと、
    前記複合体における前記標的分子の、前記標的分子単体に対する分極相互作用エネルギーを算出する処理Ｂと、
    前記複合体における前記一の薬候補分子の、前記一の薬候補分子単体に対する分極相互作用エネルギーを算出する処理Ｃと、
    を含み、
    更に、前記標的分子と、前記一の薬候補分子との前記一の相互作用エネルギーを算出する際に、前記計算対象水分子を考慮するときには、前記計算対象水分子と、前記一の薬候補分子とのペア相互作用エネルギーを算出する処理Ｄを含む、
    ことを特徴とするプログラム。
  6. 標的分子と、複数の薬候補分子の各々との相互作用エネルギーの少なくともいずれかに影響を与えるが、前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの全てには影響を与えない計算対象水分子を抽出する抽出部と、
    前記相互作用エネルギーの各々の算出において陽に扱う水分子として前記計算対象水分子のみを考慮し、かつ、前記計算対象水分子が影響を与える相互作用エネルギーの算出においてのみ前記計算対象水分子を考慮して、前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの各々を算出する算出部と、
    を備え、
    前記抽出部が、
    前記標的分子と、前記複数の薬候補分子における一の薬候補分子とを含む一の相互作用算出対象空間における、
    前記標的分子、及び前記一の薬候補分子の双方と相互作用する相互作用水分子を抽出する第１の処理と、
    前記相互作用水分子のうちから、前記相互作用水分子の座標の揺らぎ幅、及び前記相互作用水分子の存在頻度の少なくともいずれかに基づいて、計算候補水分子を抽出する第２の処理と、
    を行い、
    更に、前記抽出部が、
    前記標的分子と、前記複数の薬候補分子の各々とを含む全ての相互作用算出対象空間のうちの少なくとも前記計算候補水分子が存在する複数の相互作用算出対象空間を対比して、前記計算候補水分子のうちから、前記全ての相互作用算出対象空間において同じ座標にはない水分子を抽出する処理、
    を行い、
    前記算出部における前記相互作用エネルギーの各々の算出が、
    前記標的分子と、前記複数の薬候補分子における一の薬候補分子との一の相互作用エネルギーを算出する際に、
    前記標的分子と、前記一の薬候補分子との複合体における、前記標的分子と、前記一の薬候補分子とのペア相互作用エネルギーを算出する処理Ａと、
    前記複合体における前記標的分子の、前記標的分子単体に対する分極相互作用エネルギーを算出する処理Ｂと、
    前記複合体における前記一の薬候補分子の、前記一の薬候補分子単体に対する分極相互作用エネルギーを算出する処理Ｃと、
    を含み、
    更に、前記標的分子と、前記一の薬候補分子との前記一の相互作用エネルギーを算出する際に、前記計算対象水分子を考慮するときには、前記計算対象水分子と、前記一の薬候補分子とのペア相互作用エネルギーを算出する処理Ｄを含む、
    ことを特徴とする相互作用エネルギーの算出装置。
  7. 標的分子と、複数の薬候補分子の各々との相互作用エネルギーの少なくともいずれかに影響を与えるが、前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの全てには影響を与えず、
    前記相互作用エネルギーの各々の算出において陽に扱う水分子として唯一考慮され、かつ、計算対象水分子が影響を与える相互作用エネルギーの算出においてのみ考慮され、
    前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの各々を算出するのに用いられる、
    計算対象水分子、
    を抽出する抽出部を備え、
    前記抽出部が、
    前記標的分子と、前記複数の薬候補分子における一の薬候補分子とを含む一の相互作用算出対象空間における、
    前記標的分子、及び前記一の薬候補分子の双方と相互作用する相互作用水分子を抽出する第１の処理と、
    前記相互作用水分子のうちから、前記相互作用水分子の座標の揺らぎ幅、及び前記相互作用水分子の存在頻度の少なくともいずれかに基づいて、計算候補水分子を抽出する第２の処理と、
    を行い、
    更に、前記抽出部が、
    前記標的分子と、前記複数の薬候補分子の各々とを含む全ての相互作用算出対象空間のうちの少なくとも前記計算候補水分子が存在する複数の相互作用算出対象空間を対比して、前記計算候補水分子のうちから、前記全ての相互作用算出対象空間において同じ座標にはない水分子を抽出する処理、
    を行い、
    前記相互作用エネルギーの各々の算出が、
    前記標的分子と、前記複数の薬候補分子における一の薬候補分子との一の相互作用エネルギーを算出する際に、
    前記標的分子と、前記一の薬候補分子との複合体における、前記標的分子と、前記一の薬候補分子とのペア相互作用エネルギーを算出する処理Ａと、
    前記複合体における前記標的分子の、前記標的分子単体に対する分極相互作用エネルギーを算出する処理Ｂと、
    前記複合体における前記一の薬候補分子の、前記一の薬候補分子単体に対する分極相互作用エネルギーを算出する処理Ｃと、
    を含み、
    更に、前記標的分子と、前記一の薬候補分子との前記一の相互作用エネルギーを算出する際に、前記計算対象水分子を考慮するときには、前記計算対象水分子と、前記一の薬候補分子とのペア相互作用エネルギーを算出する処理Ｄを含む、
    ことを特徴とする相互作用エネルギーの算出装置。
  8. 標的分子と、複数の薬候補分子の各々との相互作用エネルギーの少なくともいずれかに影響を与えるが、前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの全てには影響を与えない、抽出された計算対象水分子のみを、前記相互作用エネルギーの各々の算出において陽に扱う水分子として考慮し、かつ、前記計算対象水分子が影響を与える相互作用エネルギーの算出においてのみ前記計算対象水分子を考慮して、前記標的分子と、前記複数の薬候補分子の各々との前記相互作用エネルギーの各々を算出する算出部を備え、
    前記算出部における前記相互作用エネルギーの各々の算出が、
    前記標的分子と、前記複数の薬候補分子における一の薬候補分子との一の相互作用エネルギーを算出する際に、
    前記標的分子と、前記一の薬候補分子との複合体における、前記標的分子と、前記一の薬候補分子とのペア相互作用エネルギーを算出する処理Ａと、
    前記複合体における前記標的分子の、前記標的分子単体に対する分極相互作用エネルギーを算出する処理Ｂと、
    前記複合体における前記一の薬候補分子の、前記一の薬候補分子単体に対する分極相互作用エネルギーを算出する処理Ｃと、
    を含み、
    更に、前記標的分子と、前記一の薬候補分子との前記一の相互作用エネルギーを算出する際に、前記計算対象水分子を考慮するときには、前記計算対象水分子と、前記一の薬候補分子とのペア相互作用エネルギーを算出する処理Ｄを含み、
    前記計算対象水分子の抽出が、
    前記標的分子と、前記複数の薬候補分子における一の薬候補分子とを含む一の相互作用算出対象空間における、
    前記標的分子、及び前記一の薬候補分子の双方と相互作用する相互作用水分子を抽出する第１の処理と、
    前記相互作用水分子のうちから、前記相互作用水分子の座標の揺らぎ幅、及び前記相互作用水分子の存在頻度の少なくともいずれかに基づいて、計算候補水分子を抽出する第２の処理と、
    を含み、
    更に、前記計算対象水分子の抽出が、
    前記標的分子と、前記複数の薬候補分子の各々とを含む全ての相互作用算出対象空間のうちの少なくとも前記計算候補水分子が存在する複数の相互作用算出対象空間を対比して、前記計算候補水分子のうちから、前記全ての相互作用算出対象空間において同じ座標にはない水分子を抽出する処理、
    を含む、
    ことを特徴とする相互作用エネルギーの算出装置。
JP2016565634A 2014-12-24 2014-12-24 相互作用エネルギーの算出方法、及び算出装置、並びにプログラム Active JP6409880B2 (ja)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2014/084022 WO2016103336A1 (ja) 2014-12-24 2014-12-24 相互作用エネルギーの算出方法、及び算出装置、並びにプログラム

Publications (2)

Publication Number Publication Date
JPWO2016103336A1 JPWO2016103336A1 (ja) 2017-09-28
JP6409880B2 true JP6409880B2 (ja) 2018-10-24

Family

ID=56149438

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2016565634A Active JP6409880B2 (ja) 2014-12-24 2014-12-24 相互作用エネルギーの算出方法、及び算出装置、並びにプログラム

Country Status (4)

Country Link
US (1) US10650913B2 (ja)
EP (1) EP3239874B1 (ja)
JP (1) JP6409880B2 (ja)
WO (1) WO2016103336A1 (ja)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3646250A1 (en) * 2017-05-30 2020-05-06 GTN Ltd Tensor network machine learning system

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4315960B2 (ja) 2001-01-19 2009-08-19 株式会社インシリコサイエンス タンパク質のリガンド結合部位の特定方法およびタンパク質−リガンド複合体の立体構造構築方法
US7577553B2 (en) * 2002-07-10 2009-08-18 Numerate, Inc. Method and apparatus for molecular mechanics analysis of molecular systems
RU2265244C2 (ru) * 2003-09-22 2005-11-27 Алгодайн Ллс Способ моделирования и прогнозирования связывания молекул-лигандов с молекулами-мишенями методами квантовой механики с учетом влияния растворителя
US7467046B2 (en) * 2003-10-01 2008-12-16 Duke University Geranylgeranyl transferase type I (GGTase-I) structure and uses thereof
US7756674B2 (en) * 2007-08-03 2010-07-13 The Trustees Of Columbia University In The City Of New York Methods of calculating differences of binding affinities between congeneric pairs of ligands by way of a displaced solvent functional
JP5348681B2 (ja) * 2008-11-14 2013-11-20 Ｎｅｃソフト株式会社 水溶性複合体中における、分子間相互作用エネルギーの評価方法
WO2011006158A2 (en) * 2009-07-10 2011-01-13 University Of Maryland, Baltimore Targeting nad biosynthesis in bacterial pathogens
WO2013092922A2 (en) * 2011-12-22 2013-06-27 Universitat De Barcelona Method of binding site and binding energy determination by mixed explicit solvent simulations
JP6079411B2 (ja) * 2013-04-25 2017-02-15 富士通株式会社 座標データの変換方法、相互作用の計算方法、プログラム、記録媒体及び装置

Also Published As

Publication number Publication date
EP3239874A4 (en) 2018-01-03
WO2016103336A1 (ja) 2016-06-30
US10650913B2 (en) 2020-05-12
EP3239874A1 (en) 2017-11-01
EP3239874B1 (en) 2021-04-07
US20170270281A1 (en) 2017-09-21
JPWO2016103336A1 (ja) 2017-09-28

Similar Documents

Publication Publication Date Title
Lee et al. Toward fast and accurate binding affinity prediction with pmemdGTI: An efficient implementation of GPU-accelerated thermodynamic integration
DeLuca et al. Fully flexible docking of medium sized ligand libraries with RosettaLigand
Brandsdal et al. Free energy calculations and ligand binding
Li et al. istar: A web platform for large-scale protein-ligand docking
Roe et al. Evaluation of enhanced sampling provided by accelerated molecular dynamics with Hamiltonian replica exchange methods
Gill et al. Emerging role of bioinformatics tools and software in evolution of clinical research
Cao et al. Power analysis of transcriptome-wide association study: implications for practical protocol choice
Kirsopp et al. Quantum computational quantification of protein–ligand interactions
JP6186785B2 (ja) 結合自由エネルギーの算出方法、及び結合自由エネルギーの算出装置、プログラム、並びに化合物のスクリーニング方法
Jiao et al. Computational drug repositioning for cancer therapeutics
WO2017013802A1 (ja) 結合自由エネルギーの算出方法、及び算出装置、並びにプログラム
Anandakrishnan et al. An n log n generalized Born approximation
JP6409880B2 (ja) 相互作用エネルギーの算出方法、及び算出装置、並びにプログラム
JP7379810B2 (ja) 結合自由エネルギーの算出方法、及び算出装置、並びにプログラム
JP6610182B2 (ja) 結合自由エネルギー計算の前処理方法、結合自由エネルギーの算出方法、及び装置、並びにプログラム
Huang et al. Mean-force scoring functions for protein–ligand binding
JP6623697B2 (ja) 相互作用エネルギーの算出方法、及び算出装置、並びにプログラム
JP6974766B2 (ja) 安定結合構造の算出方法、及び算出装置、並びにプログラム
WO2017199279A1 (ja) 結合自由エネルギーの算出方法、及び算出装置、並びにプログラム
JP6940752B2 (ja) プローブ分子の配置方法、及び配置装置、標的分子の結合サイトの探索方法、及び探索装置、並びにプログラム
JP6311320B2 (ja) 結合構造の算出方法、及び算出装置、プログラム、並びに記録媒体
JP6079411B2 (ja) 座標データの変換方法、相互作用の計算方法、プログラム、記録媒体及び装置
JP6488728B2 (ja) アンカー点の決定方法、結合自由エネルギーの算出方法、及び算出装置、並びにプログラム
Yamasaki et al. ToGo-WF: prediction of RNA tertiary structures and RNA–RNA/protein interactions using the KNIME workflow
WO2016072027A1 (ja) 安定結合構造の算出方法、及び算出装置、並びにプログラム

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20170608

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20170608

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20180501

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20180524

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20180828

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20180910

R150 Certificate of patent or registration of utility model

Ref document number: 6409880

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150