結合評価方法
【課題】高分子化合物と評価対象化合物との結合部位を従来よりも効率良く解明すると供に、高分子化合物と評価対象化合物との結合位置の信憑性を従来よりも向上させる。
【解決手段】高分子化合物に対する評価対象分子の結合を数値計算に基づいて評価する方法であって、評価対象分子を構成する各原子の電子受容性及び電子供与性の度合いに基づいて高分子化合物のどの原子と評価対象分子のどの原子とが結合するかを特定する。
【解決手段】高分子化合物に対する評価対象分子の結合を数値計算に基づいて評価する方法であって、評価対象分子を構成する各原子の電子受容性及び電子供与性の度合いに基づいて高分子化合物のどの原子と評価対象分子のどの原子とが結合するかを特定する。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、結合評価方法に関する。
【背景技術】
【０００２】
薬の開発では、生体を構成する高分子化合物の1つであるタンパク質と薬候補剤との結合が評価される。すなわち、タンパク質に対する結合性の低いものは薬候補剤として適当ではなく、そのような低結合性のものは薬候補剤から除外される。また、上記結合性の評価手法としては、Ｘ線構造解析や古典分子動力学法が一般的に用いられている。
【０００３】
例えば下記の先行技術文献には、古典分子動力学法を用いた結合性の解析事例として、タンパク質であるアデニル酸シクラーゼ（タイプII型）と薬候補剤としてのフォルスコリンとの結合メカニズムの検討結果が開示されている。この先行技術文献では、アデニル酸シクラーゼの活性部位の中心にフォルスコリン等の各種薬候補剤を置いた解析モデル（系）の内部エネルギー（＝位置エネルギー＋運動エネルギー）を古典分子動力学法を用いてそれぞれ算出することにより当該内部エネルギーが比較的小さいもの（つまり系をエネルギー的に安定化するもの）をアデニル酸シクラーゼに対して結合性の優れた薬候補剤として抽出すると共に、このようにして抽出された薬候補剤の１つであるフォルスコリンについてアデニル酸シクラーゼとの結合部位を検討している。
【先行技術文献】
【非特許文献】
【０００４】
【非特許文献１】Takeshi Onda et al., J. Bio. Chem., 276(2001) 47785-47793
【発明の概要】
【発明が解決しようとする課題】
【０００５】
ところで、上記古典分子動力学法はアデニル酸シクラーゼの活性部位と薬候補剤との結合性について一応の評価結果を与えてくれるものの、それのみによってアデニル酸シクラーゼと薬候補剤との結合部位とその結合状態とを明らかにすることはできない。例えば上記先行技術文献では、フォルスコリンの側鎖に結合している原子あるいは分子を種々の原子や分子と置換することによってフォルスコリンとアデニル酸シクラーゼの活性部位との結合部位を明らかにしているが、この作業には長期間と膨大な労力を必要とする実験に基づくものであり、よって非効率であり実用的ではない。
【０００６】
また、上記先行技術文献の結合性評価手法はアデニル酸シクラーゼの活性部位の中心に薬候補剤を置いた解析モデルの内部エネルギーを評価するものであるが、信憑性に問題がありより信憑性のある評価手法が望まれている。すなわち、上記先行技術文献の結合性評価手法では、アデニル酸シクラーゼの活性部位に対する薬候補剤の位置が内部エネルギーに及ぼす影響を考慮していないので、信憑性が低い。本願の発明者は、フォルスコリンの位置をアデニル酸シクラーゼの活性部位の中心からずらした解析モデルの方がフォルスコリンをアデニル酸シクラーゼの活性部位の中心に置いた解析モデルよりも内部エネルギーが小さくなるという計算結果を得た。
【０００７】
本発明は、上述した事情に鑑みたものであり、以下の点を目的とするものである。
（１）高分子化合物と評価対象化合物との結合部位を従来よりも効率良く解明する。
（２）高分子化合物と評価対象化合物との結合位置の信憑性を従来よりも向上させる。
【課題を解決するための手段】
【０００８】
上記目的を達成するために、本発明では、高分子化合物に対する評価対象分子の結合を数値計算に基づいて評価する方法であって、前記高分子化合物の活性部位に対する前記評価対象分子の結合位置を複数設定して得られた、前記高分子化合物と前記評価対象分子とから成る各解析モデルのそれぞれについて、古典分子動力学法を用いて位置エネルギーと運動エネルギーとの和である内部エネルギーを数値計算し、該内部エネルギーが最小値となる結合位置を前記高分子化合物の活性部位に対する前記評価対象分子の結合位置として特定する第１ステップと、前記第１ステップで特定された結合位置にて結合する前記高分子化合物と前記評価対象分子とから成る解析モデルについて、前記評価対象分子を構成する各原子の原子核を中心として波動関数Φ（ｒ）が「０」に漸近する距離ｒを半径とする３次元空間を限定領域として設定し、当該限定領域についての波動関数Φ（ｒ）と有効ポテンシャルＶ（ｒ）の初期値を設定する第２ステップと、前記第２ステップで設定した波動関数Φ（ｒ）の初期値を下記（１）式に代入することで、前記限定領域の電荷密度ρ（ｒ）の初期値を算出する第３ステップと、前記第２ステップで設定した有効ポテンシャルＶ（ｒ）の初期値と、前記第３ステップで算出した電荷密度ρ（ｒ）の初期値とに基づいて、下記（２）式で表されるシュレディンガーの波動方程式を解くことにより、前記限定領域の波動関数Φ（ｒ）を算出する第４ステップと、前記第４ステップで算出した波動関数Φ（ｒ）を下記（１）式に代入することで、前記限定領域の電荷密度ρ（ｒ）を算出する第５ステップと、 前記第５ステップで算出した電荷密度ρ（ｒ）と前回算出した電荷密度ρ（ｒ）とが等しいか否かを判断し、両者が等しくない場合には、前記シュレディンガーの波動方程式を解くことで得られた波動関数Φ（ｒ）及び有効ポテンシャルＶ（ｒ）を新たな初期値として第４ステップに戻る一方、両者が等しい場合には、次ステップに移行する第６ステップと、前記第６ステップで等しいと判断された電荷密度ρ（ｒ）を下記（３）及び（４）式で表される福井関数に代入することで、前記限定領域における電子受容性の度合いを示すＬＵＭＯｆ＋（ｒ）及び電子供与性の度合いを示すＨＯＭＯｆ−（ｒ）を算出する第７ステップと、を有し、前記第７ステップによるＬＵＭＯｆ＋（ｒ）及びＨＯＭＯｆ−（ｒ）の算出結果に基づいて、前記限定領域において前記高分子化合物のどの原子と前記評価対象分子のどの原子とが結合するかを特定する、という解決手段を採用する。
【０００９】
【数１】
【発明の効果】
【００１０】
本発明によれば、高分子化合物と評価対象化合物との結合部位を、評価対象分子及び高分子化合物の評価対象分子に隣接する部分を含む所定限定領域における電子状態の不安定な部位として数値計算によって検出するので、電子数が極めて多い高分子化合物を含む解析対象を効率的に数値解析することが可能であり、したがって従来の実験に基づく結合部位の解析手法に比較して結合部位の解析効率を飛躍的に向上させることができる。
【図面の簡単な説明】
【００１１】
【図１】本発明の一実施形態に係わる結合評価方法の処理手順を示すフローチャートである。
【図２】本発明の一実施形態においてアデニル酸シクラーゼＡの活性部位に対するフォルスコリンＢの結合位置を特定するための解析モデルの模式図である。
【図３】本発明の一実施形態において、動径方向の波動関数の変化特性の一例として、水素原子に属する１ｓ軌道の電子の波動関数の変化を示すグラフである。
【図４】本発明の一実施形態において、フォルスコリン近傍領域におけるＬＵＭＯｆ＋(ｒ)及びＨＯＭＯｆ−(ｒ)をIsosurfaceとして示す特性図である。
【発明を実施するための形態】
【００１２】
以下、図面を参照して、本発明の一実施形態について説明する。
なお、本実施形態は、高分子化合物としてのタンパク質、より具体的にはアデニル酸シクラーゼＡの活性部位と評価対象化合物としての薬候補剤、より具体的にはフォルスコリンＢとの結合評価方法に関するものである。
【００１３】
図１は、本結合評価方法の処理手順を示すフローチャートである。この処理手順において、ステップＳ1は前処理としての結合位置検出工程を示し、ステップＳ2〜Ｓ9は本処理としての結合部位検出工程を示している。
【００１４】
さて、結合位置検出工程では、アデニル酸シクラーゼＡとフォルスコリンＢとから成る解析モデルの内部エネルギー（＝位置エネルギー＋運動エネルギー）をコンピュータを用いた数値解析に基づいて評価することにより結合位置を特定する（ステップＳ1）。
【００１５】
図２は、結合位置を特定するための解析モデルを模式的に示した模式図である。フォルスコリンＢとアデニル酸シクラーゼＡとが結合する場合、フォルスコリンＢは、破線で示すアデニル酸シクラーゼＡの活性部位のいずれかの位置に結合することになるが、本結合評価方法では、フォルスコリンＢの活性部位に対する位置を複数設定した状態における解析対象モデルの各内部エネルギーを古典分子動力学法を用いて各々数値計算し、各位置に対応する内部エネルギーのうち最小値となるものに対応するフォルスコリンＢの位置をアデニル酸シクラーゼＡの活性部位におけるフォルスコリンＢの結合位置として特定する。
【００１６】
例えば、図２に示すフォルスコリンＢの位置ｐ1,ｐ2，ｐ3は代表としてピックアップした３箇所であるが、フォルスコリンＢを各位置ｐ1,ｐ2，ｐ3に位置させた各解析対象モデルの内部エネルギーは以下のようになる。
位置ｐ1：−１２２４７７．４３３kcal/mol
位置ｐ2：−１２０８８９．４６３kcal/mol
位置ｐ3：−１２４１５１．４１３kcal/mol
【００１７】
これら内部エネルギーの相互比較から判るように、３つの位置ｐ1,ｐ2，ｐ3のうちフォルスコリンＢを位置ｐ3に配置した場合の内部エネルギーが最も小さい。これら３つの位置ｐ1,ｐ2，ｐ3の他に複数の位置について内部エネルギーを求めたが、結果的に位置ｐ3に配置した場合の内部エネルギーが最小値となる。すなわち、フォルスコリンＢは、アデニル酸シクラーゼＡの活性部位のうち位置ｐ3に置かれた状態のときに内部エネルギー的に最も安定した状態で活性部位と結合するので、フォルスコリンＢの活性部位に対する結合位置は位置ｐ3であると特定される。
【００１８】
このようにしてフォルスコリンＢの結合位置が特定されると、当該結合位置におけるフォルスコリンＢとアデニル酸シクラーゼＡとの結合部位がステップＳ2〜Ｓ9の処理に基づいて特定される。このうちステップＳ2〜Ｓ6は、フォルスコリンＢ及びアデニル酸シクラーゼＡのフォルスコリンＢに隣接する部分を含む所定の限定領域Ｇについて波動関数Φ(ｒ)を求めることにより当該限定領域Ｇの電荷密度ρ(ｒ)をコンピュータを用いた数値計算によって求める工程（電荷密度算出工程）である。なお、波動関数Φ(ｒ)及び電荷密度ρ(ｒ)において、ｒは３次元空間上の座標を示す変数である。
【００１９】
ここで、上記限定領域Ｇの設定方法について図３をも参照して説明する。図３は、水素原子に属する電子（１ｓ軌道の電子）の波動関数Ｒ(ｋ)を原子核を基点とした動径方向の距離ｋを変数として数値計算したグラフである。なお、このグラフでは、便宜上波動関数Ｒ(ｋ)に距離ｋを乗算した値を縦軸とし、また横軸は原子核からの距離を単位a.u.（１a.u.＝０．５２９オングストローム）で示している。
【００２０】
波動関数Ｒ(ｋ)は、周知のように電子の存在確率を示すものであり、このグラフでも原子核の位置では「０」であるが、距離ｋが大きくなるに従って大きくなり、原子核に比較的近いある距離において最大値を取ると、その後徐々に減少して「０」に漸近する。すなわち、波動関数Ｒ(ｋ)が有意な値を取り得る距離範囲は原子核近傍に限定されており、この範囲よりも遠方の距離範囲では殆ど「０」とみなすことができる。このような波動関数Ｒ(ｋ)の変化特性は、水素原子に属する電子だけではなく、他の原子、他の軌道に属する電子についても同様である。
【００２１】
本実施形態では、このような波動関数の変化特徴を考慮することにより、波動関数が有意な値を取り得る距離範囲を上記限定領域Ｇとし、つまり３次元位置を表す変数ｒの範囲を波動関数が十分に大きな値を取り得る領域に限定し、この限定領域Ｇの波動関数Φ(ｒ)のみを求めることにより当該限定領域Ｇの電荷密度ρ(ｒ)を算出する。波動関数Φ(ｒ)が有意な値を取り得る距離範囲は原子の種類及び電子軌道によって多少異なるが、原子核から１０〜２０オングストローム離れた距離では、殆どの原子及び電子軌道における電子の波動関数Φ(ｒ)は「０」に漸近していると考えることができる。
【００２２】
以上の観点から本実施形態では、限定領域Ｇを規定するしきい値として１０〜２０オングストロームを採用し、これによってフォルスコリンＢを構成する各原子の原子核を中心として半径１０〜２０オングストローム内の３次元空間を限定領域Ｇに設定している。この限定領域Ｇは、フォルスコリンＢとアデニル酸シクラーゼＡのフォルスコリンＢに隣接する部分とを含む３次元空間である。このような限定領域Ｇを設定することにより、以下に説明する電荷密度算出工程における計算量が限定領域Ｇを設けない場合に比較して大幅に減少する。
【００２３】
さて、電荷密度算出工程では、最初に上記限定領域Ｇについての波動関数Φ(ｒ)と有効ポテンシャルＶ(ｒ)の各初期値を設定し（ステップＳ2）、波動関数Φ(ｒ)の初期値を次式（１）に代入することにより限定領域Ｇの全ての電子に基づく電荷密度ρ(ｒ)の初期値
を算出する（ステップＳ3）。
【００２４】
【数２】
【００２５】
そして、上記有効ポテンシャルＶ(ｒ)の初期値及び電荷密度ρ(ｒ)の初期値に基づいて以下のシュレディンガーの波動方程式（２）を解くことにより、上記限定領域Ｇの波動関数Φ(ｒ)を算出する（ステップＳ4）。なお、このシュレディンガーの波動方程式（２）において、εは内部エネルギーの期待値である。
【００２６】
【数３】
【００２７】
さらに、シュレディンガーの波動方程式（２）から求められた波動関数Φ(ｒ)を式（１）に代入することによって電荷密度ρ(ｒ)を求め（ステップＳ5）、この電荷密度ρ(ｒ)を前回の計算で求められた電荷密度ρ(ｒ)とが等しいか否かを判断し（ステップＳ6）、両者が等しくない場合（ＮＯの場合）は、ステップＳ5におけるシュレディンガーの波動方程式（２）を解くことによって求められた波動関数Φ(ｒ)及び有効ポテンシャルＶ(ｒ)を新たな初期値としてステップＳ4，Ｓ5の処理を繰り返すことにより、新たな波動関数Φ(ｒ)及び電荷密度ρ(ｒ)を求める。
【００２８】
すなわち、今回の計算で求められた電荷密度ρ(ｒ)が前回の計算で求められた電荷密度ρ(ｒ)と等しくなるまでシュレディンガーの波動方程式（２）を繰り返し解くことにより、シュレディンガーの波動方程式（２）を満足する、つまり前回の電荷密度ρ(ｒ)と等しい電荷密度ρ(ｒ)を求める。
【００２９】
このような電荷密度ρ(ｒ)が求まると、ステップＳ6の判断はＹＥＳとなり、この電荷密度ρ(ｒ)を以下の福井関数（３），（４）に代入することにより、上記限定領域Ｇにおける福井関数（３），（４）の値としてのＬＵＭＯｆ＋(ｒ)及びＨＯＭＯｆ−(ｒ)が算出される（ステップＳ7）。
【００３０】
【数４】
【００３１】
【数５】
【００３２】
このＬＵＭＯｆ＋(ｒ)及びＨＯＭＯｆ−(ｒ)は、限定領域Ｇにおける電子状態の不安定さの度合いを示すものであり、一方のＬＵＭＯｆ＋(ｒ)は電子受容性（電子を受け取り易い性質）の度合いを示し、他方のＨＯＭＯｆ−(ｒ)は電子供与性（電子を放出し易い性質）を示している。このようなＬＵＭＯｆ＋(ｒ)及びＨＯＭＯｆ−(ｒ)を画像表示し（ステップＳ8）、この表示結果に基づいて限定領域ＧにおけるフォルスコリンＢとアデニル酸シクラーゼＡとの結合部位、つまりフォルスコリンＢを構成するどの原子とアデニル酸シクラーゼＡを構成するどの原子とが結合するかを特定する（ステップＳ9）。
【００３３】
図４は、上記限定領域Ｇ内のフォルスコリン近傍領域におけるＬＵＭＯｆ＋(ｒ)及びＨＯＭＯｆ−(ｒ)を「Isosurface」として示す特性図である。すなわち、この特性図では、このフォルスコリン近傍領域における最大値に対して９０％に相当するＬＵＭＯｆ＋(ｒ)及びＨＯＭＯｆ−(ｒ)に対応する部分（空間領域）を「Isosurface」とし、この「Isosurface」をフォルスコリンＢの分子構造に対応させて表示している。
【００３４】
図４（ａ）はＬＵＭＯｆ＋(ｒ)に関する上記「Isosurface」を示しているが、この図４（ａ）とフォルスコリンＢの分子構造を示す図４（ｃ）とを対比させると解るように、図４（ａ）において比較的大きな空間として「Isosurface」が表示されている矢印の部分は、フォルスコリンＢの側鎖の１つである酸素原子（Ｏ）に該当する部分である。
【００３５】
一方、図４（ｂ）はＨＯＭＯｆ−(ｒ)に関する「Isosurface」を示しているが、この図４（ｂ）とフォルスコリンＢの分子構造を示す図４（ｃ）とを対比させると解るように、図４（ｂ）において比較的大きな空間として「Isosurface」が表示されている矢印の部分は、フォルスコリンＢの側鎖の１つであるＯＨ基（ＯＨ）に該当する部分である。
【００３６】
この図４から、フォルスコリンＢの側鎖である酸素原子（Ｏ）は、フォルスコリンＢを構成する各原子の中でも電子受容性が著しく高く、またフォルスコリンＢの側鎖であるＯＨ基（ＯＨ）は、フォルスコリンＢを構成する各原子の中でも電子供与性が著しく高いことが分かる。フォルスコリンＢとアデニル酸シクラーゼＡの活性部位との結合は両者に各々属する原子の共有結合によって行われるので、フォルスコリンＢの上記酸素原子（Ｏ）及びＯＨ基（ＯＨ）は、アデニル酸シクラーゼＡの活性部位と結合する部位（結合部位）と考えられる。
【００３７】
そして、このフォルスコリンＢにおける酸素原子（Ｏ）及びＯＨ基（ＯＨ）がアデニル酸シクラーゼＡの活性部位との結合部位であるという本結合評価方法の結果は、従来技術として示した先行技術文献の結果（実験結果）と完全に一致する。
【００３８】
一方、上記限定領域Ｇにおいてアデニル酸シクラーゼＡのフォルスコリンＢに隣接する部分に関するＬＵＭＯｆ＋(ｒ)及びＨＯＭＯｆ−(ｒ)を上記福井関数（２），（３）に基づいて算出し、かつ「Isosurface」として表示することにより、アデニル酸シクラーゼＡのフォルスコリンＢに隣接する部分について電子受容性及び電子供与性が高い部位（原子）を特定する。この結果も、従来技術として示した先行技術文献の結果（実験結果）と完全に一致する。なお、限定領域Ｇにおいてアデニル酸シクラーゼＡのフォルスコリンＢに隣接する部分に関するＬＵＭＯｆ＋(ｒ)及びＨＯＭＯｆ−(ｒ)を「Isosurface」として示す特性図については、原子が入り組んでいて図示し難いので省略する。
【００３９】
すなわち、限定領域ＧにおけるフォルスコリンＢ及びアデニル酸シクラーゼＡのフォルスコリンＢに隣接する部分に関するＬＵＭＯｆ＋(ｒ)及びＨＯＭＯｆ−(ｒ)に基づいてアデニル酸シクラーゼＡの活性部位とフォルスコリンＢとの結合部位を特定する本結合評価方法は実験結果と完全に一致する十分に有効な結合部位評価手法である。
【００４０】
また、本結合評価方法は、限定領域Ｇ内の原子のみを数値解析の対象とするので、フォルスコリンＢ及びアデニル酸シクラーゼＡの全構成原子を数値解析の対象とする場合に比較して波動関数Φ(ｒ)及び電荷密度ρ(ｒ)の算出に要する時間を大幅に短縮することが可能であり、よって従来の実験に基づく結合部位の解析手法に比較して結合部位の解析効率を飛躍的に向上させることが可能である。
【００４１】
なお、本発明は上記実施形態に限定されるものではなく、例えば以下のような変形が考えられる。
（１）上記実施形態は評価対象化合物としてのフォルスコリンＢと高分子化合物の１つであるアデニル酸シクラーゼＡとの結合に関するものであるが、本発明はこれに限定されるものではない。すなわち、高分子化合物はアデニル酸シクラーゼＡに限定されず、他のタンパク質あるいは他の高分子化合物であっても良く、また評価対象化合物はフォルスコリンＢに限定されない。
【００４２】
（２）上記実施形態では、限定領域Ｇを規定するしきい値として１０〜２０オングストロームを採用することによりフォルスコリンＢを構成する各原子の原子核を中心として半径１０〜２０オングストローム内の３次元空間を限定領域Ｇに設定し、このような限定領域Ｇを設定することによりＬＵＭＯｆ＋(ｒ)及びＨＯＭＯｆ−(ｒ)の算出に必要となる電荷密度ρ(ｒ)の算出における計算量を大幅に削減して解析効率を飛躍的に向上させている。そして、このような限定領域Ｇは、フォルスコリンの構成原子に属する電子の波動関数Φ(ｒ)が原子核から１０〜２０オングストローム離れた距離において「０」に漸近していると考えることができるという知見に基づいている。
【００４３】
しかしながら、限定領域Ｇを規定する上記しきい値、つまり波動関数Φ(ｒ)が有意な値を取るとみなす距離は１０〜２０オングストロームに限定されるものではない。また、フォルスコリンの構成原子の種類や電子軌道によって波動関数Φ(ｒ)の変化特性は異なるので、波動関数Φ(ｒ)が有意な値を取る距離は構成原子の種類や電子軌道に応じて異なる。したがって、構成原子の種類や電子軌道毎に限定領域Ｇを規定するしきい値を設定するようにしても良い。
【００４４】
（３）上記実施形態では電子状態の不安定な部位を福井関数（３），（４）を用いて特定するが、電子状態の不安定な部位の特定方法は福井関数（３），（４）に限定されるものではない。他の評価関数を用いて電子状態の不安定な部位を特定しても良い。
【符号の説明】
【００４５】
Ａ…アデニル酸シクラーゼ、Ｂ…フォルスコリン、Ｇ…限定領域
【技術分野】
【０００１】
本発明は、結合評価方法に関する。
【背景技術】
【０００２】
薬の開発では、生体を構成する高分子化合物の1つであるタンパク質と薬候補剤との結合が評価される。すなわち、タンパク質に対する結合性の低いものは薬候補剤として適当ではなく、そのような低結合性のものは薬候補剤から除外される。また、上記結合性の評価手法としては、Ｘ線構造解析や古典分子動力学法が一般的に用いられている。
【０００３】
例えば下記の先行技術文献には、古典分子動力学法を用いた結合性の解析事例として、タンパク質であるアデニル酸シクラーゼ（タイプII型）と薬候補剤としてのフォルスコリンとの結合メカニズムの検討結果が開示されている。この先行技術文献では、アデニル酸シクラーゼの活性部位の中心にフォルスコリン等の各種薬候補剤を置いた解析モデル（系）の内部エネルギー（＝位置エネルギー＋運動エネルギー）を古典分子動力学法を用いてそれぞれ算出することにより当該内部エネルギーが比較的小さいもの（つまり系をエネルギー的に安定化するもの）をアデニル酸シクラーゼに対して結合性の優れた薬候補剤として抽出すると共に、このようにして抽出された薬候補剤の１つであるフォルスコリンについてアデニル酸シクラーゼとの結合部位を検討している。
【先行技術文献】
【非特許文献】
【０００４】
【非特許文献１】Takeshi Onda et al., J. Bio. Chem., 276(2001) 47785-47793
【発明の概要】
【発明が解決しようとする課題】
【０００５】
ところで、上記古典分子動力学法はアデニル酸シクラーゼの活性部位と薬候補剤との結合性について一応の評価結果を与えてくれるものの、それのみによってアデニル酸シクラーゼと薬候補剤との結合部位とその結合状態とを明らかにすることはできない。例えば上記先行技術文献では、フォルスコリンの側鎖に結合している原子あるいは分子を種々の原子や分子と置換することによってフォルスコリンとアデニル酸シクラーゼの活性部位との結合部位を明らかにしているが、この作業には長期間と膨大な労力を必要とする実験に基づくものであり、よって非効率であり実用的ではない。
【０００６】
また、上記先行技術文献の結合性評価手法はアデニル酸シクラーゼの活性部位の中心に薬候補剤を置いた解析モデルの内部エネルギーを評価するものであるが、信憑性に問題がありより信憑性のある評価手法が望まれている。すなわち、上記先行技術文献の結合性評価手法では、アデニル酸シクラーゼの活性部位に対する薬候補剤の位置が内部エネルギーに及ぼす影響を考慮していないので、信憑性が低い。本願の発明者は、フォルスコリンの位置をアデニル酸シクラーゼの活性部位の中心からずらした解析モデルの方がフォルスコリンをアデニル酸シクラーゼの活性部位の中心に置いた解析モデルよりも内部エネルギーが小さくなるという計算結果を得た。
【０００７】
本発明は、上述した事情に鑑みたものであり、以下の点を目的とするものである。
（１）高分子化合物と評価対象化合物との結合部位を従来よりも効率良く解明する。
（２）高分子化合物と評価対象化合物との結合位置の信憑性を従来よりも向上させる。
【課題を解決するための手段】
【０００８】
上記目的を達成するために、本発明では、高分子化合物に対する評価対象分子の結合を数値計算に基づいて評価する方法であって、前記高分子化合物の活性部位に対する前記評価対象分子の結合位置を複数設定して得られた、前記高分子化合物と前記評価対象分子とから成る各解析モデルのそれぞれについて、古典分子動力学法を用いて位置エネルギーと運動エネルギーとの和である内部エネルギーを数値計算し、該内部エネルギーが最小値となる結合位置を前記高分子化合物の活性部位に対する前記評価対象分子の結合位置として特定する第１ステップと、前記第１ステップで特定された結合位置にて結合する前記高分子化合物と前記評価対象分子とから成る解析モデルについて、前記評価対象分子を構成する各原子の原子核を中心として波動関数Φ（ｒ）が「０」に漸近する距離ｒを半径とする３次元空間を限定領域として設定し、当該限定領域についての波動関数Φ（ｒ）と有効ポテンシャルＶ（ｒ）の初期値を設定する第２ステップと、前記第２ステップで設定した波動関数Φ（ｒ）の初期値を下記（１）式に代入することで、前記限定領域の電荷密度ρ（ｒ）の初期値を算出する第３ステップと、前記第２ステップで設定した有効ポテンシャルＶ（ｒ）の初期値と、前記第３ステップで算出した電荷密度ρ（ｒ）の初期値とに基づいて、下記（２）式で表されるシュレディンガーの波動方程式を解くことにより、前記限定領域の波動関数Φ（ｒ）を算出する第４ステップと、前記第４ステップで算出した波動関数Φ（ｒ）を下記（１）式に代入することで、前記限定領域の電荷密度ρ（ｒ）を算出する第５ステップと、 前記第５ステップで算出した電荷密度ρ（ｒ）と前回算出した電荷密度ρ（ｒ）とが等しいか否かを判断し、両者が等しくない場合には、前記シュレディンガーの波動方程式を解くことで得られた波動関数Φ（ｒ）及び有効ポテンシャルＶ（ｒ）を新たな初期値として第４ステップに戻る一方、両者が等しい場合には、次ステップに移行する第６ステップと、前記第６ステップで等しいと判断された電荷密度ρ（ｒ）を下記（３）及び（４）式で表される福井関数に代入することで、前記限定領域における電子受容性の度合いを示すＬＵＭＯｆ＋（ｒ）及び電子供与性の度合いを示すＨＯＭＯｆ−（ｒ）を算出する第７ステップと、を有し、前記第７ステップによるＬＵＭＯｆ＋（ｒ）及びＨＯＭＯｆ−（ｒ）の算出結果に基づいて、前記限定領域において前記高分子化合物のどの原子と前記評価対象分子のどの原子とが結合するかを特定する、という解決手段を採用する。
【０００９】
【数１】
【発明の効果】
【００１０】
本発明によれば、高分子化合物と評価対象化合物との結合部位を、評価対象分子及び高分子化合物の評価対象分子に隣接する部分を含む所定限定領域における電子状態の不安定な部位として数値計算によって検出するので、電子数が極めて多い高分子化合物を含む解析対象を効率的に数値解析することが可能であり、したがって従来の実験に基づく結合部位の解析手法に比較して結合部位の解析効率を飛躍的に向上させることができる。
【図面の簡単な説明】
【００１１】
【図１】本発明の一実施形態に係わる結合評価方法の処理手順を示すフローチャートである。
【図２】本発明の一実施形態においてアデニル酸シクラーゼＡの活性部位に対するフォルスコリンＢの結合位置を特定するための解析モデルの模式図である。
【図３】本発明の一実施形態において、動径方向の波動関数の変化特性の一例として、水素原子に属する１ｓ軌道の電子の波動関数の変化を示すグラフである。
【図４】本発明の一実施形態において、フォルスコリン近傍領域におけるＬＵＭＯｆ＋(ｒ)及びＨＯＭＯｆ−(ｒ)をIsosurfaceとして示す特性図である。
【発明を実施するための形態】
【００１２】
以下、図面を参照して、本発明の一実施形態について説明する。
なお、本実施形態は、高分子化合物としてのタンパク質、より具体的にはアデニル酸シクラーゼＡの活性部位と評価対象化合物としての薬候補剤、より具体的にはフォルスコリンＢとの結合評価方法に関するものである。
【００１３】
図１は、本結合評価方法の処理手順を示すフローチャートである。この処理手順において、ステップＳ1は前処理としての結合位置検出工程を示し、ステップＳ2〜Ｓ9は本処理としての結合部位検出工程を示している。
【００１４】
さて、結合位置検出工程では、アデニル酸シクラーゼＡとフォルスコリンＢとから成る解析モデルの内部エネルギー（＝位置エネルギー＋運動エネルギー）をコンピュータを用いた数値解析に基づいて評価することにより結合位置を特定する（ステップＳ1）。
【００１５】
図２は、結合位置を特定するための解析モデルを模式的に示した模式図である。フォルスコリンＢとアデニル酸シクラーゼＡとが結合する場合、フォルスコリンＢは、破線で示すアデニル酸シクラーゼＡの活性部位のいずれかの位置に結合することになるが、本結合評価方法では、フォルスコリンＢの活性部位に対する位置を複数設定した状態における解析対象モデルの各内部エネルギーを古典分子動力学法を用いて各々数値計算し、各位置に対応する内部エネルギーのうち最小値となるものに対応するフォルスコリンＢの位置をアデニル酸シクラーゼＡの活性部位におけるフォルスコリンＢの結合位置として特定する。
【００１６】
例えば、図２に示すフォルスコリンＢの位置ｐ1,ｐ2，ｐ3は代表としてピックアップした３箇所であるが、フォルスコリンＢを各位置ｐ1,ｐ2，ｐ3に位置させた各解析対象モデルの内部エネルギーは以下のようになる。
位置ｐ1：−１２２４７７．４３３kcal/mol
位置ｐ2：−１２０８８９．４６３kcal/mol
位置ｐ3：−１２４１５１．４１３kcal/mol
【００１７】
これら内部エネルギーの相互比較から判るように、３つの位置ｐ1,ｐ2，ｐ3のうちフォルスコリンＢを位置ｐ3に配置した場合の内部エネルギーが最も小さい。これら３つの位置ｐ1,ｐ2，ｐ3の他に複数の位置について内部エネルギーを求めたが、結果的に位置ｐ3に配置した場合の内部エネルギーが最小値となる。すなわち、フォルスコリンＢは、アデニル酸シクラーゼＡの活性部位のうち位置ｐ3に置かれた状態のときに内部エネルギー的に最も安定した状態で活性部位と結合するので、フォルスコリンＢの活性部位に対する結合位置は位置ｐ3であると特定される。
【００１８】
このようにしてフォルスコリンＢの結合位置が特定されると、当該結合位置におけるフォルスコリンＢとアデニル酸シクラーゼＡとの結合部位がステップＳ2〜Ｓ9の処理に基づいて特定される。このうちステップＳ2〜Ｓ6は、フォルスコリンＢ及びアデニル酸シクラーゼＡのフォルスコリンＢに隣接する部分を含む所定の限定領域Ｇについて波動関数Φ(ｒ)を求めることにより当該限定領域Ｇの電荷密度ρ(ｒ)をコンピュータを用いた数値計算によって求める工程（電荷密度算出工程）である。なお、波動関数Φ(ｒ)及び電荷密度ρ(ｒ)において、ｒは３次元空間上の座標を示す変数である。
【００１９】
ここで、上記限定領域Ｇの設定方法について図３をも参照して説明する。図３は、水素原子に属する電子（１ｓ軌道の電子）の波動関数Ｒ(ｋ)を原子核を基点とした動径方向の距離ｋを変数として数値計算したグラフである。なお、このグラフでは、便宜上波動関数Ｒ(ｋ)に距離ｋを乗算した値を縦軸とし、また横軸は原子核からの距離を単位a.u.（１a.u.＝０．５２９オングストローム）で示している。
【００２０】
波動関数Ｒ(ｋ)は、周知のように電子の存在確率を示すものであり、このグラフでも原子核の位置では「０」であるが、距離ｋが大きくなるに従って大きくなり、原子核に比較的近いある距離において最大値を取ると、その後徐々に減少して「０」に漸近する。すなわち、波動関数Ｒ(ｋ)が有意な値を取り得る距離範囲は原子核近傍に限定されており、この範囲よりも遠方の距離範囲では殆ど「０」とみなすことができる。このような波動関数Ｒ(ｋ)の変化特性は、水素原子に属する電子だけではなく、他の原子、他の軌道に属する電子についても同様である。
【００２１】
本実施形態では、このような波動関数の変化特徴を考慮することにより、波動関数が有意な値を取り得る距離範囲を上記限定領域Ｇとし、つまり３次元位置を表す変数ｒの範囲を波動関数が十分に大きな値を取り得る領域に限定し、この限定領域Ｇの波動関数Φ(ｒ)のみを求めることにより当該限定領域Ｇの電荷密度ρ(ｒ)を算出する。波動関数Φ(ｒ)が有意な値を取り得る距離範囲は原子の種類及び電子軌道によって多少異なるが、原子核から１０〜２０オングストローム離れた距離では、殆どの原子及び電子軌道における電子の波動関数Φ(ｒ)は「０」に漸近していると考えることができる。
【００２２】
以上の観点から本実施形態では、限定領域Ｇを規定するしきい値として１０〜２０オングストロームを採用し、これによってフォルスコリンＢを構成する各原子の原子核を中心として半径１０〜２０オングストローム内の３次元空間を限定領域Ｇに設定している。この限定領域Ｇは、フォルスコリンＢとアデニル酸シクラーゼＡのフォルスコリンＢに隣接する部分とを含む３次元空間である。このような限定領域Ｇを設定することにより、以下に説明する電荷密度算出工程における計算量が限定領域Ｇを設けない場合に比較して大幅に減少する。
【００２３】
さて、電荷密度算出工程では、最初に上記限定領域Ｇについての波動関数Φ(ｒ)と有効ポテンシャルＶ(ｒ)の各初期値を設定し（ステップＳ2）、波動関数Φ(ｒ)の初期値を次式（１）に代入することにより限定領域Ｇの全ての電子に基づく電荷密度ρ(ｒ)の初期値
を算出する（ステップＳ3）。
【００２４】
【数２】
【００２５】
そして、上記有効ポテンシャルＶ(ｒ)の初期値及び電荷密度ρ(ｒ)の初期値に基づいて以下のシュレディンガーの波動方程式（２）を解くことにより、上記限定領域Ｇの波動関数Φ(ｒ)を算出する（ステップＳ4）。なお、このシュレディンガーの波動方程式（２）において、εは内部エネルギーの期待値である。
【００２６】
【数３】
【００２７】
さらに、シュレディンガーの波動方程式（２）から求められた波動関数Φ(ｒ)を式（１）に代入することによって電荷密度ρ(ｒ)を求め（ステップＳ5）、この電荷密度ρ(ｒ)を前回の計算で求められた電荷密度ρ(ｒ)とが等しいか否かを判断し（ステップＳ6）、両者が等しくない場合（ＮＯの場合）は、ステップＳ5におけるシュレディンガーの波動方程式（２）を解くことによって求められた波動関数Φ(ｒ)及び有効ポテンシャルＶ(ｒ)を新たな初期値としてステップＳ4，Ｓ5の処理を繰り返すことにより、新たな波動関数Φ(ｒ)及び電荷密度ρ(ｒ)を求める。
【００２８】
すなわち、今回の計算で求められた電荷密度ρ(ｒ)が前回の計算で求められた電荷密度ρ(ｒ)と等しくなるまでシュレディンガーの波動方程式（２）を繰り返し解くことにより、シュレディンガーの波動方程式（２）を満足する、つまり前回の電荷密度ρ(ｒ)と等しい電荷密度ρ(ｒ)を求める。
【００２９】
このような電荷密度ρ(ｒ)が求まると、ステップＳ6の判断はＹＥＳとなり、この電荷密度ρ(ｒ)を以下の福井関数（３），（４）に代入することにより、上記限定領域Ｇにおける福井関数（３），（４）の値としてのＬＵＭＯｆ＋(ｒ)及びＨＯＭＯｆ−(ｒ)が算出される（ステップＳ7）。
【００３０】
【数４】
【００３１】
【数５】
【００３２】
このＬＵＭＯｆ＋(ｒ)及びＨＯＭＯｆ−(ｒ)は、限定領域Ｇにおける電子状態の不安定さの度合いを示すものであり、一方のＬＵＭＯｆ＋(ｒ)は電子受容性（電子を受け取り易い性質）の度合いを示し、他方のＨＯＭＯｆ−(ｒ)は電子供与性（電子を放出し易い性質）を示している。このようなＬＵＭＯｆ＋(ｒ)及びＨＯＭＯｆ−(ｒ)を画像表示し（ステップＳ8）、この表示結果に基づいて限定領域ＧにおけるフォルスコリンＢとアデニル酸シクラーゼＡとの結合部位、つまりフォルスコリンＢを構成するどの原子とアデニル酸シクラーゼＡを構成するどの原子とが結合するかを特定する（ステップＳ9）。
【００３３】
図４は、上記限定領域Ｇ内のフォルスコリン近傍領域におけるＬＵＭＯｆ＋(ｒ)及びＨＯＭＯｆ−(ｒ)を「Isosurface」として示す特性図である。すなわち、この特性図では、このフォルスコリン近傍領域における最大値に対して９０％に相当するＬＵＭＯｆ＋(ｒ)及びＨＯＭＯｆ−(ｒ)に対応する部分（空間領域）を「Isosurface」とし、この「Isosurface」をフォルスコリンＢの分子構造に対応させて表示している。
【００３４】
図４（ａ）はＬＵＭＯｆ＋(ｒ)に関する上記「Isosurface」を示しているが、この図４（ａ）とフォルスコリンＢの分子構造を示す図４（ｃ）とを対比させると解るように、図４（ａ）において比較的大きな空間として「Isosurface」が表示されている矢印の部分は、フォルスコリンＢの側鎖の１つである酸素原子（Ｏ）に該当する部分である。
【００３５】
一方、図４（ｂ）はＨＯＭＯｆ−(ｒ)に関する「Isosurface」を示しているが、この図４（ｂ）とフォルスコリンＢの分子構造を示す図４（ｃ）とを対比させると解るように、図４（ｂ）において比較的大きな空間として「Isosurface」が表示されている矢印の部分は、フォルスコリンＢの側鎖の１つであるＯＨ基（ＯＨ）に該当する部分である。
【００３６】
この図４から、フォルスコリンＢの側鎖である酸素原子（Ｏ）は、フォルスコリンＢを構成する各原子の中でも電子受容性が著しく高く、またフォルスコリンＢの側鎖であるＯＨ基（ＯＨ）は、フォルスコリンＢを構成する各原子の中でも電子供与性が著しく高いことが分かる。フォルスコリンＢとアデニル酸シクラーゼＡの活性部位との結合は両者に各々属する原子の共有結合によって行われるので、フォルスコリンＢの上記酸素原子（Ｏ）及びＯＨ基（ＯＨ）は、アデニル酸シクラーゼＡの活性部位と結合する部位（結合部位）と考えられる。
【００３７】
そして、このフォルスコリンＢにおける酸素原子（Ｏ）及びＯＨ基（ＯＨ）がアデニル酸シクラーゼＡの活性部位との結合部位であるという本結合評価方法の結果は、従来技術として示した先行技術文献の結果（実験結果）と完全に一致する。
【００３８】
一方、上記限定領域Ｇにおいてアデニル酸シクラーゼＡのフォルスコリンＢに隣接する部分に関するＬＵＭＯｆ＋(ｒ)及びＨＯＭＯｆ−(ｒ)を上記福井関数（２），（３）に基づいて算出し、かつ「Isosurface」として表示することにより、アデニル酸シクラーゼＡのフォルスコリンＢに隣接する部分について電子受容性及び電子供与性が高い部位（原子）を特定する。この結果も、従来技術として示した先行技術文献の結果（実験結果）と完全に一致する。なお、限定領域Ｇにおいてアデニル酸シクラーゼＡのフォルスコリンＢに隣接する部分に関するＬＵＭＯｆ＋(ｒ)及びＨＯＭＯｆ−(ｒ)を「Isosurface」として示す特性図については、原子が入り組んでいて図示し難いので省略する。
【００３９】
すなわち、限定領域ＧにおけるフォルスコリンＢ及びアデニル酸シクラーゼＡのフォルスコリンＢに隣接する部分に関するＬＵＭＯｆ＋(ｒ)及びＨＯＭＯｆ−(ｒ)に基づいてアデニル酸シクラーゼＡの活性部位とフォルスコリンＢとの結合部位を特定する本結合評価方法は実験結果と完全に一致する十分に有効な結合部位評価手法である。
【００４０】
また、本結合評価方法は、限定領域Ｇ内の原子のみを数値解析の対象とするので、フォルスコリンＢ及びアデニル酸シクラーゼＡの全構成原子を数値解析の対象とする場合に比較して波動関数Φ(ｒ)及び電荷密度ρ(ｒ)の算出に要する時間を大幅に短縮することが可能であり、よって従来の実験に基づく結合部位の解析手法に比較して結合部位の解析効率を飛躍的に向上させることが可能である。
【００４１】
なお、本発明は上記実施形態に限定されるものではなく、例えば以下のような変形が考えられる。
（１）上記実施形態は評価対象化合物としてのフォルスコリンＢと高分子化合物の１つであるアデニル酸シクラーゼＡとの結合に関するものであるが、本発明はこれに限定されるものではない。すなわち、高分子化合物はアデニル酸シクラーゼＡに限定されず、他のタンパク質あるいは他の高分子化合物であっても良く、また評価対象化合物はフォルスコリンＢに限定されない。
【００４２】
（２）上記実施形態では、限定領域Ｇを規定するしきい値として１０〜２０オングストロームを採用することによりフォルスコリンＢを構成する各原子の原子核を中心として半径１０〜２０オングストローム内の３次元空間を限定領域Ｇに設定し、このような限定領域Ｇを設定することによりＬＵＭＯｆ＋(ｒ)及びＨＯＭＯｆ−(ｒ)の算出に必要となる電荷密度ρ(ｒ)の算出における計算量を大幅に削減して解析効率を飛躍的に向上させている。そして、このような限定領域Ｇは、フォルスコリンの構成原子に属する電子の波動関数Φ(ｒ)が原子核から１０〜２０オングストローム離れた距離において「０」に漸近していると考えることができるという知見に基づいている。
【００４３】
しかしながら、限定領域Ｇを規定する上記しきい値、つまり波動関数Φ(ｒ)が有意な値を取るとみなす距離は１０〜２０オングストロームに限定されるものではない。また、フォルスコリンの構成原子の種類や電子軌道によって波動関数Φ(ｒ)の変化特性は異なるので、波動関数Φ(ｒ)が有意な値を取る距離は構成原子の種類や電子軌道に応じて異なる。したがって、構成原子の種類や電子軌道毎に限定領域Ｇを規定するしきい値を設定するようにしても良い。
【００４４】
（３）上記実施形態では電子状態の不安定な部位を福井関数（３），（４）を用いて特定するが、電子状態の不安定な部位の特定方法は福井関数（３），（４）に限定されるものではない。他の評価関数を用いて電子状態の不安定な部位を特定しても良い。
【符号の説明】
【００４５】
Ａ…アデニル酸シクラーゼ、Ｂ…フォルスコリン、Ｇ…限定領域
【特許請求の範囲】
【請求項１】
高分子化合物に対する評価対象分子の結合を数値計算に基づいて評価する方法であって、
前記高分子化合物の活性部位に対する前記評価対象分子の結合位置を複数設定して得られた、前記高分子化合物と前記評価対象分子とから成る各解析モデルのそれぞれについて、古典分子動力学法を用いて位置エネルギーと運動エネルギーとの和である内部エネルギーを数値計算し、該内部エネルギーが最小値となる結合位置を前記高分子化合物の活性部位に対する前記評価対象分子の結合位置として特定する第１ステップと、
前記第１ステップで特定された結合位置にて結合する前記高分子化合物と前記評価対象分子とから成る解析モデルについて、前記評価対象分子を構成する各原子の原子核を中心として波動関数Φ（ｒ）が「０」に漸近する距離ｒを半径とする３次元空間を限定領域として設定し、当該限定領域についての波動関数Φ（ｒ）と有効ポテンシャルＶ（ｒ）の初期値を設定する第２ステップと、
前記第２ステップで設定した波動関数Φ（ｒ）の初期値を下記（１）式に代入することで、前記限定領域の電荷密度ρ（ｒ）の初期値を算出する第３ステップと、
前記第２ステップで設定した有効ポテンシャルＶ（ｒ）の初期値と、前記第３ステップで算出した電荷密度ρ（ｒ）の初期値とに基づいて、下記（２）式で表されるシュレディンガーの波動方程式を解くことにより、前記限定領域の波動関数Φ（ｒ）を算出する第４ステップと、
前記第４ステップで算出した波動関数Φ（ｒ）を下記（１）式に代入することで、前記限定領域の電荷密度ρ（ｒ）を算出する第５ステップと、
前記第５ステップで算出した電荷密度ρ（ｒ）と前回算出した電荷密度ρ（ｒ）とが等しいか否かを判断し、両者が等しくない場合には、前記シュレディンガーの波動方程式を解くことで得られた波動関数Φ（ｒ）及び有効ポテンシャルＶ（ｒ）を新たな初期値として第４ステップに戻る一方、両者が等しい場合には、次ステップに移行する第６ステップと、
前記第６ステップで等しいと判断された電荷密度ρ（ｒ）を下記（３）及び（４）式で表される福井関数に代入することで、前記限定領域における電子受容性の度合いを示すＬＵＭＯｆ＋（ｒ）及び電子供与性の度合いを示すＨＯＭＯｆ−（ｒ）を算出する第７ステップと、を有し、
前記第７ステップによるＬＵＭＯｆ＋（ｒ）及びＨＯＭＯｆ−（ｒ）の算出結果に基づいて、前記限定領域において前記高分子化合物のどの原子と前記評価対象分子のどの原子とが結合するかを特定することを特徴とする結合評価方法。
【数１】
【請求項２】
前記第７ステップで算出したＬＵＭＯｆ＋（ｒ）及びＨＯＭＯｆ−（ｒ）を画像表示する第８ステップと、
前記第８ステップによる画像表示結果に基づいて、前記限定領域において前記高分子化合物のどの原子と前記評価対象分子のどの原子とが結合するかを特定する第９ステップと、
をさらに有することを特徴とする請求項１記載の結合評価方法。
【請求項３】
前記高分子化合物は、タンパク質であることを特徴とする請求項１または２記載の結合評価方法。
【請求項４】
評価対象分子は薬候補剤であることを特徴とする請求項１〜３のいずれか一項に記載の結合評価方法。
【請求項５】
前記タンパク質は、アデニル酸シクラーゼであることを特徴とする請求項３記載の結合評価方法。
【請求項６】
前記薬候補剤は、フォルスコリンであることを特徴とする請求項４記載の結合評価方法。
【請求項１】
高分子化合物に対する評価対象分子の結合を数値計算に基づいて評価する方法であって、
前記高分子化合物の活性部位に対する前記評価対象分子の結合位置を複数設定して得られた、前記高分子化合物と前記評価対象分子とから成る各解析モデルのそれぞれについて、古典分子動力学法を用いて位置エネルギーと運動エネルギーとの和である内部エネルギーを数値計算し、該内部エネルギーが最小値となる結合位置を前記高分子化合物の活性部位に対する前記評価対象分子の結合位置として特定する第１ステップと、
前記第１ステップで特定された結合位置にて結合する前記高分子化合物と前記評価対象分子とから成る解析モデルについて、前記評価対象分子を構成する各原子の原子核を中心として波動関数Φ（ｒ）が「０」に漸近する距離ｒを半径とする３次元空間を限定領域として設定し、当該限定領域についての波動関数Φ（ｒ）と有効ポテンシャルＶ（ｒ）の初期値を設定する第２ステップと、
前記第２ステップで設定した波動関数Φ（ｒ）の初期値を下記（１）式に代入することで、前記限定領域の電荷密度ρ（ｒ）の初期値を算出する第３ステップと、
前記第２ステップで設定した有効ポテンシャルＶ（ｒ）の初期値と、前記第３ステップで算出した電荷密度ρ（ｒ）の初期値とに基づいて、下記（２）式で表されるシュレディンガーの波動方程式を解くことにより、前記限定領域の波動関数Φ（ｒ）を算出する第４ステップと、
前記第４ステップで算出した波動関数Φ（ｒ）を下記（１）式に代入することで、前記限定領域の電荷密度ρ（ｒ）を算出する第５ステップと、
前記第５ステップで算出した電荷密度ρ（ｒ）と前回算出した電荷密度ρ（ｒ）とが等しいか否かを判断し、両者が等しくない場合には、前記シュレディンガーの波動方程式を解くことで得られた波動関数Φ（ｒ）及び有効ポテンシャルＶ（ｒ）を新たな初期値として第４ステップに戻る一方、両者が等しい場合には、次ステップに移行する第６ステップと、
前記第６ステップで等しいと判断された電荷密度ρ（ｒ）を下記（３）及び（４）式で表される福井関数に代入することで、前記限定領域における電子受容性の度合いを示すＬＵＭＯｆ＋（ｒ）及び電子供与性の度合いを示すＨＯＭＯｆ−（ｒ）を算出する第７ステップと、を有し、
前記第７ステップによるＬＵＭＯｆ＋（ｒ）及びＨＯＭＯｆ−（ｒ）の算出結果に基づいて、前記限定領域において前記高分子化合物のどの原子と前記評価対象分子のどの原子とが結合するかを特定することを特徴とする結合評価方法。
【数１】
【請求項２】
前記第７ステップで算出したＬＵＭＯｆ＋（ｒ）及びＨＯＭＯｆ−（ｒ）を画像表示する第８ステップと、
前記第８ステップによる画像表示結果に基づいて、前記限定領域において前記高分子化合物のどの原子と前記評価対象分子のどの原子とが結合するかを特定する第９ステップと、
をさらに有することを特徴とする請求項１記載の結合評価方法。
【請求項３】
前記高分子化合物は、タンパク質であることを特徴とする請求項１または２記載の結合評価方法。
【請求項４】
評価対象分子は薬候補剤であることを特徴とする請求項１〜３のいずれか一項に記載の結合評価方法。
【請求項５】
前記タンパク質は、アデニル酸シクラーゼであることを特徴とする請求項３記載の結合評価方法。
【請求項６】
前記薬候補剤は、フォルスコリンであることを特徴とする請求項４記載の結合評価方法。
【図１】
【図２】
【図３】
【図４】
【図２】
【図３】
【図４】
【公開番号】特開２０１３−９１６４０（Ｐ２０１３−９１６４０Ａ）
【公開日】平成２５年５月１６日（２０１３．５．１６）
【国際特許分類】
【出願番号】特願２０１２−２４３９９９（Ｐ２０１２−２４３９９９）
【出願日】平成２４年１１月５日（２０１２．１１．５）
【分割の表示】特願２００４−２０５７４６（Ｐ２００４−２０５７４６）の分割
【原出願日】平成１６年７月１３日（２００４．７．１３）
【出願人】（０００００００９９）株式会社ＩＨＩ (5,014)
【Ｆターム（参考）】
【公開日】平成２５年５月１６日（２０１３．５．１６）
【国際特許分類】
【出願日】平成２４年１１月５日（２０１２．１１．５）
【分割の表示】特願２００４−２０５７４６（Ｐ２００４−２０５７４６）の分割
【原出願日】平成１６年７月１３日（２００４．７．１３）
【出願人】（０００００００９９）株式会社ＩＨＩ (5,014)
【Ｆターム（参考）】
[ Back to top ]