評価対象化合物の活性度評価方法
【課題】高分子化合物（タンパク質、酵素）に対する評価対象化合物の活性度評価を、従来より短時間で行うことができる活性度評価方法の提供。
【解決手段】高分子化合物（例えば、アデニル酸シクラーゼ）に対する評価対象化合物の活性度を、当該評価対象化合物（例えば、フォルスコリン誘導体）を構成する原子の内、一または複数を解析ポイントとし、当該解析ポイントについて数値計算により求めた福井関数の値を変数とする所定の評価式に基づいて算出する活性度評価方法。
【解決手段】高分子化合物（例えば、アデニル酸シクラーゼ）に対する評価対象化合物の活性度を、当該評価対象化合物（例えば、フォルスコリン誘導体）を構成する原子の内、一または複数を解析ポイントとし、当該解析ポイントについて数値計算により求めた福井関数の値を変数とする所定の評価式に基づいて算出する活性度評価方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、評価対象化合物の活性度評価方法に関する。
【背景技術】
【０００２】
薬の開発では、生体を構成する高分子化合物の1つであるタンパク質に対する薬候補剤の活性度が評価される。すなわち、タンパク質に対する活性度の低いものは薬候補剤として適切ではなく、そのような低活性度のものは薬候補剤から除外される。しかるに、タンパク質には数十万以上の種類が存在し、機能が類似するにも関わらず臓器発現・生化学的特性の異なるタンパク質も存在する（タンパク質のサブタイプ）。例えば、下記先行技術文献には、薬候補剤としてフォルスコリンを側鎖修飾したフォルスコリン誘導体を各種生成し、タンパク質であるアデニル酸シクラーゼ（II型、III型、Ｖ型サブタイプ）に対する上記フォルスコリン誘導体の活性度を実験的に求めた実験結果が開示されている。
【非特許文献１】Takeshi Onda et al., J. Bio. Chem., 276(2001) 47785-47793
【発明の開示】
【発明が解決しようとする課題】
【０００３】
ところで、薬の開発においては、長い開発期間が大きな問題となっており、当該開発期間を短縮することによって競合他社に対して優位性を確保すると共に、開発コストの大幅な削減を行うことができる。しかしながら、上述したように、タンパク質、例えばアデニル酸シクラーゼに対するフォルスコリン誘導体の活性度を評価し、なおかつ複数存在するアデニル酸シクラーゼのサブタイプ毎に評価するためには、当該活性度を実験によって求める他なく、この実験には膨大な労力且つ時間と共に特殊技術を必要とするため、従来の活性度評価方法では非常に効率が悪いという問題があった。
【０００４】
本発明は、上述した事情に鑑みたものであり、従来より短時間で高分子化合物に対する評価対象化合物の活性度評価を行うことを目的とする。
【課題を解決するための手段】
【０００５】
上記目的を達成するために、本発明では、第１の解決手段として、高分子化合物に対する評価対象化合物の活性度を、当該評価対象化合物を構成する原子の内、一または複数を解析ポイントとし、当該解析ポイントについて数値計算により求めた福井関数の値を変数とする所定の評価式に基づいて算出する、という手段を採用する。
【０００６】
また、本発明では、第２の解決手段として、上記第１の解決手段において、前記解析ポイントｋ毎に求核的反応を示す第１の福井関数の値ｆk＋と、求電子的反応を示す第２の福井関数の値ｆk−とを数値計算により求め、前記第１の福井関数の値ｆk＋に係る第１係数Ｄk＋、前記第２の福井関数の値ｆk−に係る第２係数Ｄk−から成る下記評価式（５）に基づいて前記高分子化合物に対する評価対象化合物の活性度Ｐを算出することを特徴とする。
【０００７】
【数１】
【０００８】
また、本発明では、第３の解決手段として、上記第２の解決手段において、評価対象化合物の複数の誘導体に関する活性度実験値Ｐｅと、前記各誘導体に関する第１の福井関数の値ｆk＋及び第２の福井関数の値ｆk−とに最小二乗法を適用することにより、前記第１の福井関数の値ｆk＋及び第２の福井関数の値ｆk−と活性度Ｐとの関係を示す第１係数Ｄk＋及び第２係数Ｄk−を演算することを特徴とする。
【０００９】
また、本発明では、第４の解決手段として、上記第２の解決手段において、評価対象化合物の複数の誘導体に関する活性度実験値Ｐｅと、前記各誘導体に関する第１の福井関数の値ｆk＋及び第２の福井関数の値ｆk−とを前記評価式に代入して得られる複数の方程式を連立方程式として解くことにより、前記第１係数Ｄk＋及び第２係数Ｄk−を演算することを特徴とする。
【００１０】
また、本発明では、第５の解決手段として、上記第３または４の解決手段において、前記評価対象化合物の各誘導体を構成する原子の内、当該各誘導体に共通して存在している原子を解析ポイントとすることを特徴とする。
【００１１】
また、本発明では、第６の解決手段として、上記第１〜５いずれかの解決手段において、前記高分子化合物はタンパク質であることを特徴とする。
【００１２】
また、本発明では、第７の解決手段として、上記第１〜６いずれかの解決手段において、前記評価対象化合物は薬候補剤であることを特徴とする。
【発明の効果】
【００１３】
本発明によれば、高分子化合物に対する評価対象化合物の活性度を、当該評価対象化合物を構成する原子の内、一または複数を解析ポイントとし、当該解析ポイントについて数値計算により求めた福井関数の値を変数とする所定の評価式に基づいて算出するので、当該数値計算によって求めた福井関数の値を評価式に代入するだけで活性度を算出することができるので、従来のように実験に基づく活性度評価方法と比較して、非常に短時間で活性度の評価を行うことができる。
【発明を実施するための最良の形態】
【００１４】
以下、図面を参照して、本発明の一実施形態について説明する。
なお、本実施形態は、高分子化合物としてのタンパク質、より具体的にはアデニル酸シクラーゼと、評価対象化合物としての薬候補剤、より具体的にはフォルスコリン誘導体との活性度評価方法に関するものである。
【００１５】
図１は、本活性度評価方法の処理手順を示すフローチャートである。まず、評価対象となるフォルスコリン誘導体の選定を行う（ステップＳ１）。図２は、フォルスコリン誘導体の分子構造を示す模式図である。この図において、Ｒ６、Ｒ７及びＲ１３はフォルスコリン誘導体の側鎖を修飾するための原子または分子が結合する位置を示しており、当該位置にどのような原子または分子が結合するかによってフォルスコリン誘導体の物性が変化する。具体的には、アデニル酸シクラーゼに対する活性度が変化する。図３にフォルスコリン誘導体の側鎖を修飾する原子または分子の一例を示す。この図において、Ｒ６にＨ、Ｒ７にＣＨ３、Ｒ１３にＣＨ＝ＣＨ２が結合している「Forskolin」は、自然界に唯一存在するフォルスコリンであり、その他の「ＦＤ１」〜「ＦＤ６」は、側鎖の構造を変えることで人工的に生成されたフォルスコリンである。このように人工的に側鎖の構造を変えられたフォルスコリンをフォルスコリン誘導体という。本実施形態では、図３に示すフォルスコリン誘導体「ＦＤ５」を評価対象として選定し、以下説明する。なお、以下のステップＳ２〜Ｓ１３は、第１原理計算法に基づいたコンピュータシミュレーションによって、評価対象として選定したフォルスコリン誘導体のアデニル酸シクラーゼに対する活性度を求める工程である。よって、以下説明する処理は、コンピュータシミュレーション上で行われるものである。
【００１６】
さて、上記のように評価対象となるフォルスコリン誘導体「ＦＤ５」を選定すると、次に、フォルスコリン誘導体「ＦＤ５」の化学式をコンピュータシミュレーションに入力する（ステップＳ２）。そして、フォルスコリン誘導体「ＦＤ５」を構成するベンゼン環におけるＣ（炭素）の内、いずれかを解析ポイントとして決定（入力）する（ステップＳ３）。本実施形態では、図２に示すようにＣ１〜Ｃ１３までを解析ポイントとして決定する。なお、以下では、解析ポイントｋ（ｋ＝１〜１３）に該当するＣをＣｋと記載する。
【００１７】
上記のように、ベンゼン環におけるＣｋのように、各誘導体に共通して存在する、つまり安定的な位置に存在する原子を解析ポイントとして決定することで、以下に述べる活性度を精度良く求めることができる。また、解析ポイントｋの数が多い程、活性度を精度良く求めることができる。
【００１８】
続いて、フォルスコリン誘導体「ＦＤ５」の化学式に基づいて波動関数Φ（ｒ）、有効ポテンシャルＶ（ｒ）及び電荷密度ρ（ｒ）の初期値を設定し入力する（ステップＳ４）。なお、これら波動関数Φ（ｒ）、有効ポテンシャルＶ（ｒ）及び電荷密度ρ（ｒ）において、ｒは３次元空間上の座標を示す変数である。
【００１９】
波動関数Φ（ｒ）の初期値は、フォルスコリン誘導体「ＦＤ５」を構成する各原子が孤立原子として３次元空間上に存在する場合の波動関数を各原子毎に求め、当該各原子毎に求めた波動関数を全て加算したものである。同様に、有効ポテンシャルＶ（ｒ）の初期値も、フォルスコリン誘導体「ＦＤ５」を構成する各原子が孤立原子として３次元空間上に存在する場合の波動関数に基づいて有効ポテンシャルを各原子毎に求め、当該各原子毎に求めた有効ポテンシャルを全て加算したものである。また、電荷密度ρ（ｒ）の初期値は、上記波動関数Φ（ｒ）の初期値を下記演算式（１）に代入することによって求める。なお、下記演算式（１）において、Φ＊（ｒ）は、波動関数Φ（ｒ）の共役複素数である。
【００２０】
【数２】
【００２１】
次に、上記有効ポテンシャルＶ（ｒ）の初期値及び電荷密度ρ（ｒ）の初期値に基づいて下記Kohn Sham方程式（２）を解くことにより、フォルスコリン誘導体「ＦＤ５」の波動関数Φ（ｒ）とエネルギ固有値εとを算出する（ステップＳ５）。
【００２２】
【数３】
【００２３】
そして、ステップＳ４で求めたフォルスコリン誘導体「ＦＤ５」の波動関数Φ（ｒ）に基づいてフォルスコリン誘導体「ＦＤ５」の電荷密度ρ（ｒ）及び有効ポテンシャルＶ（ｒ）を算出し（ステップＳ６）、この電荷密度ρ（ｒ）と当該電荷密度ρ（ｒ）の前回値、つまりここでは初期値とが等しいか否かを判断する（ステップＳ７）。このステップＳ７において、「NO」、すなわち電荷密度ρ（ｒ）の前回値（初期値）と、ステップＳ６で求めた電荷密度ρ（ｒ）の今回値とが等しくないと判断された場合、ステップＳ６で求めた有効ポテンシャルＶ（ｒ）及び電荷密度ρ（ｒ）を新たな初期値と設定し（ステップＳ８）、ステップＳ５に移行して、再度Kohn Sham方程式（２）を解くことにより、新たな波動関数Φ（ｒ）とエネルギ固有値εとを算出する。すなわち、ステップＳ７において、電荷密度ρ（ｒ）の前回値と今回値とが等しくなるまでステップＳ５〜Ｓ８の処理を繰り返すことにより、Kohn Sham方程式（２）を満足する波動関数Φ（ｒ）及びエネルギ固有値εが求まる。
【００２４】
一方、ステップＳ７において、「YES」、すなわち電荷密度ρ（ｒ）の前回値と今回値とが等しいと判断された場合、上記のようにKohn Sham方程式（２）を満足する波動関数Φ（ｒ）及びエネルギ固有値εが求まるので、この波動関数Φ（ｒ）及びエネルギ固有値εに基づいて各原子に働く原子間力を算出すると共に、フォルスコリン誘導体「ＦＤ５」の構造の最適化を行う（ステップＳ９）。つまり、ステップＳ５〜Ｓ８の繰り返しで求めた波動関数Φ（ｒ）等は、あくまで図２に示すような２次元平面上のモデルにおいて最適な値であって、実際には３次元空間上におけるフォルスコリン誘導体「ＦＤ５」の構造を考慮する必要がある。
【００２５】
具体的には、ステップＳ９では、フォルスコリン誘導体「ＦＤ５」を構成する各原子を、３次元空間上において波動関数Φ（ｒ）から推測される最適な方向に所定の距離だけ移動させ、その時の各原子に働く原子間力を算出する。この時の原子間力が０になり、各原子が動かなくなった場合にフォルスコリン誘導体「ＦＤ５」の構造が最適化されたと判断できる。よって、移動後の各原子に働く原子間力を算出し、該原子間力が０になったか否かを判断する（ステップＳ１０）。このステップＳ１０において、「NO」、すなわち原子間力が０ではなく、構造が最適化されていない場合、各原子の移動後の構造における波動関数Φ（ｒ）を求め、当該波動関数Φ（ｒ）から求めた有効ポテンシャルＶ（ｒ）及び電荷密度ρ（ｒ）を新たな初期値と設定し（ステップＳ１１）、ステップＳ５〜Ｓ９の処理を繰り返す。ここで、ステップＳ５に戻るのは、各原子の移動後の構造変化によって波動関数Φ（ｒ）が変化するためである。
【００２６】
一方、このステップＳ１０において、「YES」、すなわち各原子に働く原子間力が０になり、フォルスコリン誘導体「ＦＤ５」の構造が最適化された場合、その最適化された構造における波動関数Φ（ｒ）に基づいて各解析ポイントＣｋにおける電荷密度ρ（ｒ）を求め、この電荷密度ρ(ｒ)を以下の求核的反応を示す第１の福井関数（３），求電子的反応を示す第２の福井関数（４）に代入することにより、各解析ポイントＣｋにおける第１の福井関数（３）の値としてのＬＵＭＯｆk＋及び第２の福井関数（４）の値としてのＨＯＭＯｆk−を算出する（ステップＳ１２）。
【００２７】
【数４】
【００２８】
そして、上記のように各解析ポイントＣｋ毎のＬＵＭＯｆk＋及びＨＯＭＯｆk−が求まると、下記活性度評価式（５）に代入することによってアデニル酸シクラーゼに対するフォルスコリン誘導体「ＦＤ５」の活性度Ｐを算出する（ステップＳ１３）。なお、アデニル酸シクラーゼにはサブタイプが存在するので、実際には、このサブタイプｉ固有の活性度評価式（６）に基づいて活性度Ｐｉを算出することになる。
【００２９】
【数５】
【００３０】
なお、活性度評価式（６）において、Ｄi,k＋はＬＵＭＯｆk＋に係る第１係数、Ｄi,k−はＨＯＭＯｆk−に係る第２係数、ｎは解析ポイントｋの総数、ｉはアデニル酸シクラーゼのサブタイプである。例えば、ｉ＝ＶならばＶ型アデニル酸シクラーゼに対する活性度Ｐｖ、第１係数Ｄｖ,k＋、第２係数Ｄｖ,k−を示し、ｉ＝IIならばII型アデニル酸シクラーゼに対する活性度ＰII等を示す。つまり、アデニル酸シクラーゼのサブタイプｉによって第１係数Ｄi,k＋及び第２係数Ｄi,k−は変化するので、アデニル酸シクラーゼのサブタイプｉ固有の活性度評価式（６）が存在する。
【００３１】
ここで、第１係数Ｄi,k＋及び第２係数Ｄi,k−の設定方法について説明する。図４は、フォルスコリン「Forskolin」及びフォルスコリン誘導体「ＦＤ１」〜「ＦＤ６」における解析ポイントＣ１〜Ｃ１３について数値計算により算出したＬＵＭＯｆk＋の値である。また、図５は、同様にフォルスコリン「Forskolin」及びフォルスコリン誘導体「ＦＤ１」〜「ＦＤ６」における解析ポイントＣ１〜Ｃ１３について数値計算により算出したＨＯＭＯｆk−の値である。図６は、サブタイプｉ＝Ｖ、つまりＶ型アデニル酸シクラーゼに対するフォルスコリン「Forskolin」及びフォルスコリン誘導体「ＦＤ１」〜「ＦＤ６」の活性度を示す実験値（図６において白色の棒グラフ）である。なお、図６において縦軸は、自然界に存在するフォルスコリン「Forskolin」の活性度を１００％とした時の活性度を示している。
【００３２】
例えば、フォルスコリン誘導体「ＦＤ５」に着目し、図４及び図５からフォルスコリン誘導体「ＦＤ５」に関するＬＵＭＯｆk＋及びＨＯＭＯｆk−を抽出し、図６からフォルスコリン誘導体「ＦＤ５」に関する活性度の実験値を抽出して活性度評価式（６）に代入すると下記方程式（７）が得られる。
【００３３】
【数６】
【００３４】
上記方程式（７）からわかるように、合計２６個の未知数である第１係数Ｄｖ,k＋及び第２係数Ｄｖ,k−が与えられる。すなわち、上記方程式（７）を２６種類のフォルスコリン誘導体について求め、これら２６個の方程式による連立一次方程式を解くことにより、Ｖ型アデニル酸シクラーゼに対する第１係数Ｄｖ,k＋及び第２係数Ｄｖ,k−を求めることができる。または、図４及び図５に示すＬＵＭＯｆk＋及びＨＯＭＯｆk−の値と、図６に示す活性度の実験値とに多変数最小二乗法を適用することにより、第１係数Ｄｖ,k＋及び第２係数Ｄｖ,k−を求めることもできる。なお、各サブタイプのアデニル酸シクラーゼに対する各種フォルスコリン誘導体の活性度に関する実験は、過去に様々な機関で行われており、上記実験値は既知のものを使用している。
【００３５】
上記のように、第１係数Ｄｖ,k＋及び第２係数Ｄｖ,k−が求まれば、上記活性度評価式（６）に代入することで、Ｖ型アデニル酸シクラーゼに対する各種フォルスコリン誘導体の活性度Ｐｖについての活性度評価式（８）が得られる。
【００３６】
【数７】
【００３７】
図６において、斜線の棒グラフは上記活性度評価式（８）を用いて算出したＶ型アデニル酸シクラーゼに対する各種フォルスコリン誘導体の活性度Ｐｖの計算値である。この図に示すように、活性度評価式（８）を用いて算出した活性度Ｐｖは、実験値と相関があることがわかる。
【００３８】
さらに、図７はサブタイプｉ＝II、つまりII型アデニル酸シクラーゼに対するフォルスコリン「Forskolin」及びフォルスコリン誘導体「ＦＤ１」〜「ＦＤ６」の活性度を示すものであり、下記活性度評価式（９）は、図７に示す活性度の実験値と、図４及び図５に示すＬＵＭＯｆk＋及びＨＯＭＯｆk−とに基づいて第１係数ＤII、k＋及び第２係数ＤII、k−を求めることによって導出したII型アデニル酸シクラーゼに関する活性度評価式である。
【００３９】
【数８】
【００４０】
図７に示すように、活性度評価式（９）を用いて算出した活性度ＰIIは、実験値と相関があることがわかる。すなわち、活性度評価式（８）、（９）のように、アデニル酸シクラーゼのサブタイプ毎に予め活性度評価式を求めておき、ステップＳ１〜Ｓ１３までの処理によって、評価対象となるフォルスコリン誘導体の各解析ポイントｋにおけるＬＵＭＯｆk＋及びＨＯＭＯｆk−を算出し、アデニル酸シクラーゼのサブタイプ毎に用意した活性度評価式に代入することで、評価対象として選定したフォルスコリン誘導体が、どのサブタイプのアデニル酸シクラーゼに対して、どの程度の活性度を示すかを計算によって求めることができるのである。
【００４１】
本願発明者は、他のサブタイプのアデニル酸シクラーゼに関する活性度評価式を求め、フォルスコリン誘導体「ＦＤ１」〜「ＦＤ６」以外の複数のフォルスコリン誘導体について当該活性度評価式に基づいて活性度を求め、実験値との相関を評価し、いずれも実験値と計算値とでは高い相関を示すことを確認した。つまり、上記活性度評価式（６）に基づく活性度評価は、非常に信頼性が高い評価方法といえる。
【００４２】
また、フォルスコリン誘導体としては理論的に無数の誘導体が合成可能であるが、上述したように、所定のアデニル酸シクラーゼに対して少なくとも解析ポイント数の２倍の種類のフォルスコリン誘導体の活性度についての実験値があれば、活性度評価式を求めることができるので、他のフォルスコリン誘導体についての活性度評価は上記活性度評価式に基づいて算出すれば良いことになり、従来と比べて非常に短時間で活性度評価を行うことができる。
【００４３】
さらに、上記活性度評価式（６）の特徴は、評価対象となるフォルスコリン誘導体における所定の解析ポイントｋについてＬＵＭＯｆk＋及びＨＯＭＯｆk−を数値計算により算出するだけでアデニル酸シクラーゼのサブタイプ毎に活性度を算出することができる、すなわち活性度を算出する過程において、フォルスコリン誘導体の結合相手であるアデニル酸シクラーゼの分子構造を全く考慮せずに活性度を算出可能であるということである。従来では、例えば、フォルスコリン誘導体とアデニル酸シクラーゼとの結合性を評価する場合、フォルスコリン誘導体の分子構造のみならず、高分子化合物であるアデニル酸シクラーゼの分子構造をも考慮する必要があり、非常に複雑で長時間を要する数値計算によって評価を行わざるを得なかった。しかしながら、本活性度評価方法では、上述したように予め実験で求められている既知の活性度を用いて、活性度評価式（６）を各サブタイプのアデニル酸シクラーゼに対して求めておけば、フォルスコリン誘導体の分子構造のみを考慮して数値計算を行うことによって短時間で活性度を算出することができる。
【００４４】
以上のように、本活性度評価方法によれば、飛躍的に活性度の評価効率を向上することができる。従って、薬の開発期間を大幅に短縮することによって競合他社に対して優位性を確保すると共に、開発コストの大幅な削減を行うことができる。
【００４５】
なお、本発明は、上記実施形態に限定されるものではなく、例えば以下のような変形例が考えられる。
【００４６】
（１） 上記実施形態では、高分子化合物としてアデニル酸シクラーゼ、評価対象化合物
としてフォルスコリン誘導体を用いて説明したが、これに限らず、高分子化合物としては他のタンパク質でも良く、また、評価対象化合物としてはフォルスコリン誘導体に限定されない。
【００４７】
（２）上記実施形態では、福井関数として求核的反応を示すＬＵＭＯｆk＋と、求電子的反応を示すＨＯＭＯｆk−とを用いて活性度評価式を構成したが、これに限らず、ラジカル反応を示す福井関数、すなわち、上記求核的反応と求電子的反応との両方を示す福井関数の値を用いて活性度評価式を構成しても良い。
【００４８】
（３）上記実施形態では、解析ポイントｋをフォルスコリン誘導体を構成するベンゼン環のＣとしたが、これに限らず、フォルスコリン誘導体を構成する他の原子でも良い。
【００４９】
（４）上記実施形態では、解析ポイントｋを１３個としたが、これに限定されない。しかしながら、活性度の計算値の精度を上げるためには、解析ポイントｋは多数の方が好ましい。
【図面の簡単な説明】
【００５０】
【図１】本発明の一実施形態に係る活性度評価方法の処理手順を示すフローチャートである。
【図２】本発明の一実施形態におけるフォルスコリン誘導体の分子構造モデル図である。
【図３】本発明の一実施形態におけるフォルスコリン誘導体に結合する側鎖の種類を示すものである。
【図４】本発明の一実施形態における各種フォルスコリン誘導体の各解析ポイントにおけるＬＵＭＯｆk＋の計算値である。
【図５】本発明の一実施形態における各種フォルスコリン誘導体の各解析ポイントにおけるＨＯＭＯｆk−の計算値である。
【図６】本発明の一実施形態におけるＶ型アデニル酸シクラーゼに対する各種フォルスコリン誘導体の活性度に関する実験値及び計算値を示す棒グラフである。
【図７】本発明の一実施形態におけるII型アデニル酸シクラーゼに対する各種フォルスコリン誘導体の活性度に関する実験値および計算値を示す棒グラフである。
【技術分野】
【０００１】
本発明は、評価対象化合物の活性度評価方法に関する。
【背景技術】
【０００２】
薬の開発では、生体を構成する高分子化合物の1つであるタンパク質に対する薬候補剤の活性度が評価される。すなわち、タンパク質に対する活性度の低いものは薬候補剤として適切ではなく、そのような低活性度のものは薬候補剤から除外される。しかるに、タンパク質には数十万以上の種類が存在し、機能が類似するにも関わらず臓器発現・生化学的特性の異なるタンパク質も存在する（タンパク質のサブタイプ）。例えば、下記先行技術文献には、薬候補剤としてフォルスコリンを側鎖修飾したフォルスコリン誘導体を各種生成し、タンパク質であるアデニル酸シクラーゼ（II型、III型、Ｖ型サブタイプ）に対する上記フォルスコリン誘導体の活性度を実験的に求めた実験結果が開示されている。
【非特許文献１】Takeshi Onda et al., J. Bio. Chem., 276(2001) 47785-47793
【発明の開示】
【発明が解決しようとする課題】
【０００３】
ところで、薬の開発においては、長い開発期間が大きな問題となっており、当該開発期間を短縮することによって競合他社に対して優位性を確保すると共に、開発コストの大幅な削減を行うことができる。しかしながら、上述したように、タンパク質、例えばアデニル酸シクラーゼに対するフォルスコリン誘導体の活性度を評価し、なおかつ複数存在するアデニル酸シクラーゼのサブタイプ毎に評価するためには、当該活性度を実験によって求める他なく、この実験には膨大な労力且つ時間と共に特殊技術を必要とするため、従来の活性度評価方法では非常に効率が悪いという問題があった。
【０００４】
本発明は、上述した事情に鑑みたものであり、従来より短時間で高分子化合物に対する評価対象化合物の活性度評価を行うことを目的とする。
【課題を解決するための手段】
【０００５】
上記目的を達成するために、本発明では、第１の解決手段として、高分子化合物に対する評価対象化合物の活性度を、当該評価対象化合物を構成する原子の内、一または複数を解析ポイントとし、当該解析ポイントについて数値計算により求めた福井関数の値を変数とする所定の評価式に基づいて算出する、という手段を採用する。
【０００６】
また、本発明では、第２の解決手段として、上記第１の解決手段において、前記解析ポイントｋ毎に求核的反応を示す第１の福井関数の値ｆk＋と、求電子的反応を示す第２の福井関数の値ｆk−とを数値計算により求め、前記第１の福井関数の値ｆk＋に係る第１係数Ｄk＋、前記第２の福井関数の値ｆk−に係る第２係数Ｄk−から成る下記評価式（５）に基づいて前記高分子化合物に対する評価対象化合物の活性度Ｐを算出することを特徴とする。
【０００７】
【数１】
【０００８】
また、本発明では、第３の解決手段として、上記第２の解決手段において、評価対象化合物の複数の誘導体に関する活性度実験値Ｐｅと、前記各誘導体に関する第１の福井関数の値ｆk＋及び第２の福井関数の値ｆk−とに最小二乗法を適用することにより、前記第１の福井関数の値ｆk＋及び第２の福井関数の値ｆk−と活性度Ｐとの関係を示す第１係数Ｄk＋及び第２係数Ｄk−を演算することを特徴とする。
【０００９】
また、本発明では、第４の解決手段として、上記第２の解決手段において、評価対象化合物の複数の誘導体に関する活性度実験値Ｐｅと、前記各誘導体に関する第１の福井関数の値ｆk＋及び第２の福井関数の値ｆk−とを前記評価式に代入して得られる複数の方程式を連立方程式として解くことにより、前記第１係数Ｄk＋及び第２係数Ｄk−を演算することを特徴とする。
【００１０】
また、本発明では、第５の解決手段として、上記第３または４の解決手段において、前記評価対象化合物の各誘導体を構成する原子の内、当該各誘導体に共通して存在している原子を解析ポイントとすることを特徴とする。
【００１１】
また、本発明では、第６の解決手段として、上記第１〜５いずれかの解決手段において、前記高分子化合物はタンパク質であることを特徴とする。
【００１２】
また、本発明では、第７の解決手段として、上記第１〜６いずれかの解決手段において、前記評価対象化合物は薬候補剤であることを特徴とする。
【発明の効果】
【００１３】
本発明によれば、高分子化合物に対する評価対象化合物の活性度を、当該評価対象化合物を構成する原子の内、一または複数を解析ポイントとし、当該解析ポイントについて数値計算により求めた福井関数の値を変数とする所定の評価式に基づいて算出するので、当該数値計算によって求めた福井関数の値を評価式に代入するだけで活性度を算出することができるので、従来のように実験に基づく活性度評価方法と比較して、非常に短時間で活性度の評価を行うことができる。
【発明を実施するための最良の形態】
【００１４】
以下、図面を参照して、本発明の一実施形態について説明する。
なお、本実施形態は、高分子化合物としてのタンパク質、より具体的にはアデニル酸シクラーゼと、評価対象化合物としての薬候補剤、より具体的にはフォルスコリン誘導体との活性度評価方法に関するものである。
【００１５】
図１は、本活性度評価方法の処理手順を示すフローチャートである。まず、評価対象となるフォルスコリン誘導体の選定を行う（ステップＳ１）。図２は、フォルスコリン誘導体の分子構造を示す模式図である。この図において、Ｒ６、Ｒ７及びＲ１３はフォルスコリン誘導体の側鎖を修飾するための原子または分子が結合する位置を示しており、当該位置にどのような原子または分子が結合するかによってフォルスコリン誘導体の物性が変化する。具体的には、アデニル酸シクラーゼに対する活性度が変化する。図３にフォルスコリン誘導体の側鎖を修飾する原子または分子の一例を示す。この図において、Ｒ６にＨ、Ｒ７にＣＨ３、Ｒ１３にＣＨ＝ＣＨ２が結合している「Forskolin」は、自然界に唯一存在するフォルスコリンであり、その他の「ＦＤ１」〜「ＦＤ６」は、側鎖の構造を変えることで人工的に生成されたフォルスコリンである。このように人工的に側鎖の構造を変えられたフォルスコリンをフォルスコリン誘導体という。本実施形態では、図３に示すフォルスコリン誘導体「ＦＤ５」を評価対象として選定し、以下説明する。なお、以下のステップＳ２〜Ｓ１３は、第１原理計算法に基づいたコンピュータシミュレーションによって、評価対象として選定したフォルスコリン誘導体のアデニル酸シクラーゼに対する活性度を求める工程である。よって、以下説明する処理は、コンピュータシミュレーション上で行われるものである。
【００１６】
さて、上記のように評価対象となるフォルスコリン誘導体「ＦＤ５」を選定すると、次に、フォルスコリン誘導体「ＦＤ５」の化学式をコンピュータシミュレーションに入力する（ステップＳ２）。そして、フォルスコリン誘導体「ＦＤ５」を構成するベンゼン環におけるＣ（炭素）の内、いずれかを解析ポイントとして決定（入力）する（ステップＳ３）。本実施形態では、図２に示すようにＣ１〜Ｃ１３までを解析ポイントとして決定する。なお、以下では、解析ポイントｋ（ｋ＝１〜１３）に該当するＣをＣｋと記載する。
【００１７】
上記のように、ベンゼン環におけるＣｋのように、各誘導体に共通して存在する、つまり安定的な位置に存在する原子を解析ポイントとして決定することで、以下に述べる活性度を精度良く求めることができる。また、解析ポイントｋの数が多い程、活性度を精度良く求めることができる。
【００１８】
続いて、フォルスコリン誘導体「ＦＤ５」の化学式に基づいて波動関数Φ（ｒ）、有効ポテンシャルＶ（ｒ）及び電荷密度ρ（ｒ）の初期値を設定し入力する（ステップＳ４）。なお、これら波動関数Φ（ｒ）、有効ポテンシャルＶ（ｒ）及び電荷密度ρ（ｒ）において、ｒは３次元空間上の座標を示す変数である。
【００１９】
波動関数Φ（ｒ）の初期値は、フォルスコリン誘導体「ＦＤ５」を構成する各原子が孤立原子として３次元空間上に存在する場合の波動関数を各原子毎に求め、当該各原子毎に求めた波動関数を全て加算したものである。同様に、有効ポテンシャルＶ（ｒ）の初期値も、フォルスコリン誘導体「ＦＤ５」を構成する各原子が孤立原子として３次元空間上に存在する場合の波動関数に基づいて有効ポテンシャルを各原子毎に求め、当該各原子毎に求めた有効ポテンシャルを全て加算したものである。また、電荷密度ρ（ｒ）の初期値は、上記波動関数Φ（ｒ）の初期値を下記演算式（１）に代入することによって求める。なお、下記演算式（１）において、Φ＊（ｒ）は、波動関数Φ（ｒ）の共役複素数である。
【００２０】
【数２】
【００２１】
次に、上記有効ポテンシャルＶ（ｒ）の初期値及び電荷密度ρ（ｒ）の初期値に基づいて下記Kohn Sham方程式（２）を解くことにより、フォルスコリン誘導体「ＦＤ５」の波動関数Φ（ｒ）とエネルギ固有値εとを算出する（ステップＳ５）。
【００２２】
【数３】
【００２３】
そして、ステップＳ４で求めたフォルスコリン誘導体「ＦＤ５」の波動関数Φ（ｒ）に基づいてフォルスコリン誘導体「ＦＤ５」の電荷密度ρ（ｒ）及び有効ポテンシャルＶ（ｒ）を算出し（ステップＳ６）、この電荷密度ρ（ｒ）と当該電荷密度ρ（ｒ）の前回値、つまりここでは初期値とが等しいか否かを判断する（ステップＳ７）。このステップＳ７において、「NO」、すなわち電荷密度ρ（ｒ）の前回値（初期値）と、ステップＳ６で求めた電荷密度ρ（ｒ）の今回値とが等しくないと判断された場合、ステップＳ６で求めた有効ポテンシャルＶ（ｒ）及び電荷密度ρ（ｒ）を新たな初期値と設定し（ステップＳ８）、ステップＳ５に移行して、再度Kohn Sham方程式（２）を解くことにより、新たな波動関数Φ（ｒ）とエネルギ固有値εとを算出する。すなわち、ステップＳ７において、電荷密度ρ（ｒ）の前回値と今回値とが等しくなるまでステップＳ５〜Ｓ８の処理を繰り返すことにより、Kohn Sham方程式（２）を満足する波動関数Φ（ｒ）及びエネルギ固有値εが求まる。
【００２４】
一方、ステップＳ７において、「YES」、すなわち電荷密度ρ（ｒ）の前回値と今回値とが等しいと判断された場合、上記のようにKohn Sham方程式（２）を満足する波動関数Φ（ｒ）及びエネルギ固有値εが求まるので、この波動関数Φ（ｒ）及びエネルギ固有値εに基づいて各原子に働く原子間力を算出すると共に、フォルスコリン誘導体「ＦＤ５」の構造の最適化を行う（ステップＳ９）。つまり、ステップＳ５〜Ｓ８の繰り返しで求めた波動関数Φ（ｒ）等は、あくまで図２に示すような２次元平面上のモデルにおいて最適な値であって、実際には３次元空間上におけるフォルスコリン誘導体「ＦＤ５」の構造を考慮する必要がある。
【００２５】
具体的には、ステップＳ９では、フォルスコリン誘導体「ＦＤ５」を構成する各原子を、３次元空間上において波動関数Φ（ｒ）から推測される最適な方向に所定の距離だけ移動させ、その時の各原子に働く原子間力を算出する。この時の原子間力が０になり、各原子が動かなくなった場合にフォルスコリン誘導体「ＦＤ５」の構造が最適化されたと判断できる。よって、移動後の各原子に働く原子間力を算出し、該原子間力が０になったか否かを判断する（ステップＳ１０）。このステップＳ１０において、「NO」、すなわち原子間力が０ではなく、構造が最適化されていない場合、各原子の移動後の構造における波動関数Φ（ｒ）を求め、当該波動関数Φ（ｒ）から求めた有効ポテンシャルＶ（ｒ）及び電荷密度ρ（ｒ）を新たな初期値と設定し（ステップＳ１１）、ステップＳ５〜Ｓ９の処理を繰り返す。ここで、ステップＳ５に戻るのは、各原子の移動後の構造変化によって波動関数Φ（ｒ）が変化するためである。
【００２６】
一方、このステップＳ１０において、「YES」、すなわち各原子に働く原子間力が０になり、フォルスコリン誘導体「ＦＤ５」の構造が最適化された場合、その最適化された構造における波動関数Φ（ｒ）に基づいて各解析ポイントＣｋにおける電荷密度ρ（ｒ）を求め、この電荷密度ρ(ｒ)を以下の求核的反応を示す第１の福井関数（３），求電子的反応を示す第２の福井関数（４）に代入することにより、各解析ポイントＣｋにおける第１の福井関数（３）の値としてのＬＵＭＯｆk＋及び第２の福井関数（４）の値としてのＨＯＭＯｆk−を算出する（ステップＳ１２）。
【００２７】
【数４】
【００２８】
そして、上記のように各解析ポイントＣｋ毎のＬＵＭＯｆk＋及びＨＯＭＯｆk−が求まると、下記活性度評価式（５）に代入することによってアデニル酸シクラーゼに対するフォルスコリン誘導体「ＦＤ５」の活性度Ｐを算出する（ステップＳ１３）。なお、アデニル酸シクラーゼにはサブタイプが存在するので、実際には、このサブタイプｉ固有の活性度評価式（６）に基づいて活性度Ｐｉを算出することになる。
【００２９】
【数５】
【００３０】
なお、活性度評価式（６）において、Ｄi,k＋はＬＵＭＯｆk＋に係る第１係数、Ｄi,k−はＨＯＭＯｆk−に係る第２係数、ｎは解析ポイントｋの総数、ｉはアデニル酸シクラーゼのサブタイプである。例えば、ｉ＝ＶならばＶ型アデニル酸シクラーゼに対する活性度Ｐｖ、第１係数Ｄｖ,k＋、第２係数Ｄｖ,k−を示し、ｉ＝IIならばII型アデニル酸シクラーゼに対する活性度ＰII等を示す。つまり、アデニル酸シクラーゼのサブタイプｉによって第１係数Ｄi,k＋及び第２係数Ｄi,k−は変化するので、アデニル酸シクラーゼのサブタイプｉ固有の活性度評価式（６）が存在する。
【００３１】
ここで、第１係数Ｄi,k＋及び第２係数Ｄi,k−の設定方法について説明する。図４は、フォルスコリン「Forskolin」及びフォルスコリン誘導体「ＦＤ１」〜「ＦＤ６」における解析ポイントＣ１〜Ｃ１３について数値計算により算出したＬＵＭＯｆk＋の値である。また、図５は、同様にフォルスコリン「Forskolin」及びフォルスコリン誘導体「ＦＤ１」〜「ＦＤ６」における解析ポイントＣ１〜Ｃ１３について数値計算により算出したＨＯＭＯｆk−の値である。図６は、サブタイプｉ＝Ｖ、つまりＶ型アデニル酸シクラーゼに対するフォルスコリン「Forskolin」及びフォルスコリン誘導体「ＦＤ１」〜「ＦＤ６」の活性度を示す実験値（図６において白色の棒グラフ）である。なお、図６において縦軸は、自然界に存在するフォルスコリン「Forskolin」の活性度を１００％とした時の活性度を示している。
【００３２】
例えば、フォルスコリン誘導体「ＦＤ５」に着目し、図４及び図５からフォルスコリン誘導体「ＦＤ５」に関するＬＵＭＯｆk＋及びＨＯＭＯｆk−を抽出し、図６からフォルスコリン誘導体「ＦＤ５」に関する活性度の実験値を抽出して活性度評価式（６）に代入すると下記方程式（７）が得られる。
【００３３】
【数６】
【００３４】
上記方程式（７）からわかるように、合計２６個の未知数である第１係数Ｄｖ,k＋及び第２係数Ｄｖ,k−が与えられる。すなわち、上記方程式（７）を２６種類のフォルスコリン誘導体について求め、これら２６個の方程式による連立一次方程式を解くことにより、Ｖ型アデニル酸シクラーゼに対する第１係数Ｄｖ,k＋及び第２係数Ｄｖ,k−を求めることができる。または、図４及び図５に示すＬＵＭＯｆk＋及びＨＯＭＯｆk−の値と、図６に示す活性度の実験値とに多変数最小二乗法を適用することにより、第１係数Ｄｖ,k＋及び第２係数Ｄｖ,k−を求めることもできる。なお、各サブタイプのアデニル酸シクラーゼに対する各種フォルスコリン誘導体の活性度に関する実験は、過去に様々な機関で行われており、上記実験値は既知のものを使用している。
【００３５】
上記のように、第１係数Ｄｖ,k＋及び第２係数Ｄｖ,k−が求まれば、上記活性度評価式（６）に代入することで、Ｖ型アデニル酸シクラーゼに対する各種フォルスコリン誘導体の活性度Ｐｖについての活性度評価式（８）が得られる。
【００３６】
【数７】
【００３７】
図６において、斜線の棒グラフは上記活性度評価式（８）を用いて算出したＶ型アデニル酸シクラーゼに対する各種フォルスコリン誘導体の活性度Ｐｖの計算値である。この図に示すように、活性度評価式（８）を用いて算出した活性度Ｐｖは、実験値と相関があることがわかる。
【００３８】
さらに、図７はサブタイプｉ＝II、つまりII型アデニル酸シクラーゼに対するフォルスコリン「Forskolin」及びフォルスコリン誘導体「ＦＤ１」〜「ＦＤ６」の活性度を示すものであり、下記活性度評価式（９）は、図７に示す活性度の実験値と、図４及び図５に示すＬＵＭＯｆk＋及びＨＯＭＯｆk−とに基づいて第１係数ＤII、k＋及び第２係数ＤII、k−を求めることによって導出したII型アデニル酸シクラーゼに関する活性度評価式である。
【００３９】
【数８】
【００４０】
図７に示すように、活性度評価式（９）を用いて算出した活性度ＰIIは、実験値と相関があることがわかる。すなわち、活性度評価式（８）、（９）のように、アデニル酸シクラーゼのサブタイプ毎に予め活性度評価式を求めておき、ステップＳ１〜Ｓ１３までの処理によって、評価対象となるフォルスコリン誘導体の各解析ポイントｋにおけるＬＵＭＯｆk＋及びＨＯＭＯｆk−を算出し、アデニル酸シクラーゼのサブタイプ毎に用意した活性度評価式に代入することで、評価対象として選定したフォルスコリン誘導体が、どのサブタイプのアデニル酸シクラーゼに対して、どの程度の活性度を示すかを計算によって求めることができるのである。
【００４１】
本願発明者は、他のサブタイプのアデニル酸シクラーゼに関する活性度評価式を求め、フォルスコリン誘導体「ＦＤ１」〜「ＦＤ６」以外の複数のフォルスコリン誘導体について当該活性度評価式に基づいて活性度を求め、実験値との相関を評価し、いずれも実験値と計算値とでは高い相関を示すことを確認した。つまり、上記活性度評価式（６）に基づく活性度評価は、非常に信頼性が高い評価方法といえる。
【００４２】
また、フォルスコリン誘導体としては理論的に無数の誘導体が合成可能であるが、上述したように、所定のアデニル酸シクラーゼに対して少なくとも解析ポイント数の２倍の種類のフォルスコリン誘導体の活性度についての実験値があれば、活性度評価式を求めることができるので、他のフォルスコリン誘導体についての活性度評価は上記活性度評価式に基づいて算出すれば良いことになり、従来と比べて非常に短時間で活性度評価を行うことができる。
【００４３】
さらに、上記活性度評価式（６）の特徴は、評価対象となるフォルスコリン誘導体における所定の解析ポイントｋについてＬＵＭＯｆk＋及びＨＯＭＯｆk−を数値計算により算出するだけでアデニル酸シクラーゼのサブタイプ毎に活性度を算出することができる、すなわち活性度を算出する過程において、フォルスコリン誘導体の結合相手であるアデニル酸シクラーゼの分子構造を全く考慮せずに活性度を算出可能であるということである。従来では、例えば、フォルスコリン誘導体とアデニル酸シクラーゼとの結合性を評価する場合、フォルスコリン誘導体の分子構造のみならず、高分子化合物であるアデニル酸シクラーゼの分子構造をも考慮する必要があり、非常に複雑で長時間を要する数値計算によって評価を行わざるを得なかった。しかしながら、本活性度評価方法では、上述したように予め実験で求められている既知の活性度を用いて、活性度評価式（６）を各サブタイプのアデニル酸シクラーゼに対して求めておけば、フォルスコリン誘導体の分子構造のみを考慮して数値計算を行うことによって短時間で活性度を算出することができる。
【００４４】
以上のように、本活性度評価方法によれば、飛躍的に活性度の評価効率を向上することができる。従って、薬の開発期間を大幅に短縮することによって競合他社に対して優位性を確保すると共に、開発コストの大幅な削減を行うことができる。
【００４５】
なお、本発明は、上記実施形態に限定されるものではなく、例えば以下のような変形例が考えられる。
【００４６】
（１） 上記実施形態では、高分子化合物としてアデニル酸シクラーゼ、評価対象化合物
としてフォルスコリン誘導体を用いて説明したが、これに限らず、高分子化合物としては他のタンパク質でも良く、また、評価対象化合物としてはフォルスコリン誘導体に限定されない。
【００４７】
（２）上記実施形態では、福井関数として求核的反応を示すＬＵＭＯｆk＋と、求電子的反応を示すＨＯＭＯｆk−とを用いて活性度評価式を構成したが、これに限らず、ラジカル反応を示す福井関数、すなわち、上記求核的反応と求電子的反応との両方を示す福井関数の値を用いて活性度評価式を構成しても良い。
【００４８】
（３）上記実施形態では、解析ポイントｋをフォルスコリン誘導体を構成するベンゼン環のＣとしたが、これに限らず、フォルスコリン誘導体を構成する他の原子でも良い。
【００４９】
（４）上記実施形態では、解析ポイントｋを１３個としたが、これに限定されない。しかしながら、活性度の計算値の精度を上げるためには、解析ポイントｋは多数の方が好ましい。
【図面の簡単な説明】
【００５０】
【図１】本発明の一実施形態に係る活性度評価方法の処理手順を示すフローチャートである。
【図２】本発明の一実施形態におけるフォルスコリン誘導体の分子構造モデル図である。
【図３】本発明の一実施形態におけるフォルスコリン誘導体に結合する側鎖の種類を示すものである。
【図４】本発明の一実施形態における各種フォルスコリン誘導体の各解析ポイントにおけるＬＵＭＯｆk＋の計算値である。
【図５】本発明の一実施形態における各種フォルスコリン誘導体の各解析ポイントにおけるＨＯＭＯｆk−の計算値である。
【図６】本発明の一実施形態におけるＶ型アデニル酸シクラーゼに対する各種フォルスコリン誘導体の活性度に関する実験値及び計算値を示す棒グラフである。
【図７】本発明の一実施形態におけるII型アデニル酸シクラーゼに対する各種フォルスコリン誘導体の活性度に関する実験値および計算値を示す棒グラフである。
【特許請求の範囲】
【請求項１】
高分子化合物に対する評価対象化合物の活性度を、当該評価対象化合物を構成する原子の内、一または複数を解析ポイントとし、当該解析ポイントについて数値計算により求めた福井関数の値を変数とする所定の評価式に基づいて算出することを特徴とする評価対象化合物の活性度評価方法。
【請求項２】
前記解析ポイントｋ毎に求核的反応を示す第１の福井関数の値ｆk＋と、求電子的反応を示す第２の福井関数の値ｆk−とを数値計算により求め、前記第１の福井関数の値ｆk＋に係る第１係数Ｄk＋、前記第２の福井関数の値ｆk−に係る第２係数Ｄk−から成る下記評価式（５）に基づいて前記高分子化合物に対する評価対象化合物の活性度Ｐを算出することを特徴とする請求項１記載の評価対象化合物の活性度評価方法。
【数１】
【請求項３】
評価対象化合物の複数の誘導体に関する活性度実験値Ｐｅと、前記各誘導体に関する第１の福井関数の値ｆk＋及び第２の福井関数の値ｆk−とに最小二乗法を適用することにより、前記第１の福井関数の値ｆk＋及び第２の福井関数の値ｆk−と活性度Ｐとの関係を示す第１係数Ｄk＋及び第２係数Ｄk−を演算することを特徴とする請求項２記載の評価対象化合物の活性度評価方法。
【請求項４】
評価対象化合物の複数の誘導体に関する活性度実験値Ｐｅと、前記各誘導体に関する第１の福井関数の値ｆk＋及び第２の福井関数の値ｆk−とを前記評価式に代入して得られる複数の方程式を連立方程式として解くことにより、前記第１係数Ｄk＋及び第２係数Ｄk−を演算することを特徴とする請求項２記載の評価対象化合物の活性度評価方法。
【請求項５】
前記評価対象化合物の各誘導体を構成する原子の内、当該各誘導体に共通して存在している原子を解析ポイントとすることを特徴とする請求項３または４記載の評価対象化合物の活性度評価方法。
【請求項６】
前記高分子化合物はタンパク質であることを特徴とする請求項1〜５いずれかに記載の評価対象化合物の活性度評価方法。
【請求項７】
前記評価対象化合物は薬候補剤であることを特徴とする請求項1〜６いずれかに記載の評価対象化合物の活性度評価方法。
【請求項１】
高分子化合物に対する評価対象化合物の活性度を、当該評価対象化合物を構成する原子の内、一または複数を解析ポイントとし、当該解析ポイントについて数値計算により求めた福井関数の値を変数とする所定の評価式に基づいて算出することを特徴とする評価対象化合物の活性度評価方法。
【請求項２】
前記解析ポイントｋ毎に求核的反応を示す第１の福井関数の値ｆk＋と、求電子的反応を示す第２の福井関数の値ｆk−とを数値計算により求め、前記第１の福井関数の値ｆk＋に係る第１係数Ｄk＋、前記第２の福井関数の値ｆk−に係る第２係数Ｄk−から成る下記評価式（５）に基づいて前記高分子化合物に対する評価対象化合物の活性度Ｐを算出することを特徴とする請求項１記載の評価対象化合物の活性度評価方法。
【数１】
【請求項３】
評価対象化合物の複数の誘導体に関する活性度実験値Ｐｅと、前記各誘導体に関する第１の福井関数の値ｆk＋及び第２の福井関数の値ｆk−とに最小二乗法を適用することにより、前記第１の福井関数の値ｆk＋及び第２の福井関数の値ｆk−と活性度Ｐとの関係を示す第１係数Ｄk＋及び第２係数Ｄk−を演算することを特徴とする請求項２記載の評価対象化合物の活性度評価方法。
【請求項４】
評価対象化合物の複数の誘導体に関する活性度実験値Ｐｅと、前記各誘導体に関する第１の福井関数の値ｆk＋及び第２の福井関数の値ｆk−とを前記評価式に代入して得られる複数の方程式を連立方程式として解くことにより、前記第１係数Ｄk＋及び第２係数Ｄk−を演算することを特徴とする請求項２記載の評価対象化合物の活性度評価方法。
【請求項５】
前記評価対象化合物の各誘導体を構成する原子の内、当該各誘導体に共通して存在している原子を解析ポイントとすることを特徴とする請求項３または４記載の評価対象化合物の活性度評価方法。
【請求項６】
前記高分子化合物はタンパク質であることを特徴とする請求項1〜５いずれかに記載の評価対象化合物の活性度評価方法。
【請求項７】
前記評価対象化合物は薬候補剤であることを特徴とする請求項1〜６いずれかに記載の評価対象化合物の活性度評価方法。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【公開番号】特開２００７−６４７９５（Ｐ２００７−６４７９５Ａ）
【公開日】平成１９年３月１５日（２００７．３．１５）
【国際特許分類】
【出願番号】特願２００５−２５１１８９（Ｐ２００５−２５１１８９）
【出願日】平成１７年８月３１日（２００５．８．３１）
【出願人】（０００００００９９）石川島播磨重工業株式会社 (5,014)
【出願人】（５０５３２８６８３）
【Ｆターム（参考）】
【公開日】平成１９年３月１５日（２００７．３．１５）
【国際特許分類】
【出願日】平成１７年８月３１日（２００５．８．３１）
【出願人】（０００００００９９）石川島播磨重工業株式会社 (5,014)
【出願人】（５０５３２８６８３）
【Ｆターム（参考）】
[ Back to top ]