(57)【要約】【課題】分子の構造の反応特性に対する相関を解析する際のユーザの負担を軽減してその解析の精度を向上させること。【解決手段】少なくとも２つの分子を所定の規則に基づいて重畳する分子重畳部(11)と、前記重畳された分子に対して、着目する原子或いは基毎に複数の領域に分割する領域分割部(12)と、前記分子重畳部で重畳された分子をメッシュに分割し、メッシュの各格子点における少なくとも静電及び立体ポテンシャルエネルギーとを含むポテンシャルエネルギーを計算するポテンシャルエネルギー計算部(13)と、前記複数の領域ごとに平均化ポテンシャルエネルギーを計算し、前記平均化ポテンシャルエネルギーと反応特性値との相関を計算する相関計算部(14)と、前記計算結果に基づいて、反応特性値を、前記ポテンシャルエネルギーを独立変数とする関数（一次式）として提示する結果提示部(15)とを備えた。【選択図】      図１
2007106708.tif 000002
【特許請求の範囲】
【請求項１】
少なくとも２つの分子を所定の規則に基づいて重畳する分子重畳部と、
前記分子重畳部で重畳された分子に対して、着目する原子或いは基毎に複数の領域に分割する領域分割部と、
前記分子重畳部で重畳された分子をメッシュに分割し、メッシュの各格子点における少なくとも静電ポテンシャルエネルギーと立体ポテンシャルエネルギーとを含むポテンシャルエネルギーを計算するポテンシャルエネルギー計算部と、
前記格子点が属する前記複数の領域ごとに平均化ポテンシャルエネルギーを計算し、前記平均化ポテンシャルエネルギーと反応特性値との相関を計算する相関計算部と、
前記相関計算部の計算結果に基づいて、反応特性値を、前記ポテンシャルエネルギーを独立変数とする一次式として提示する結果提示部とを具備することを特徴とする構造反応特性相関解析装置。
【請求項２】
請求項１に記載の構造反応特性相関解析装置において、前記分子重畳部で重畳された分子の重ね合わせの状態を記憶するデータ記憶部を更に具備することを特徴とする構造反応特性相関解析装置。
【請求項３】
請求項１又は請求項２に記載の構造反応特性相関解析装置において、前記分子重畳部における分子の重畳の前後における状況を表示する表示部と、
ユーザの入力を受け付ける入力部とを更に具備し、
前記分子重畳部は、前記入力部に入力された操作情報に基づいて、前記表示部に表示された重畳後の分子について再重畳を行うことを特徴とする構造反応特性相関解析装置。
【請求項４】
請求項１から請求項３のいずれか１項に記載の構造反応特性相関解析装置において、前記領域分割部は、前記重畳後の分子の各原子又は基を中心とするボロノイ（Voronoi）分割を用いて前記分子を複数の領域に分割することを特徴とする構造反応特性相関解析装置。
【請求項５】
請求項１から請求項４のいずれか１項に記載の構造反応特性相関解析装置において、前記相関計算部は、サポートベクタマシンを用いて、前記反応特性値と前記ポテンシャルエネルギーとの相関を計算することを特徴とする構造反応特性相関解析装置。
【請求項６】
請求項１から請求項５のいずれか１項に記載の構造反応特性相関解析装置において、前記ポテンシャルエネルギー計算部は、各分子のポテンシャルエネルギーを計算し、前記結果提示部は、前記重畳される分子の類似性を、各分子の前記ポテンシャルエネルギーに基づいて判定することを特徴とする構造反応特性相関解析装置。
【請求項７】
請求項１から請求項６のいずれか１項に記載の構造反応特性相関解析装置において、前記表示部は、前記ポテンシャルエネルギーの値の分布を色相や、明暗により明示的に表示することを特徴とする構造反応特性相関解析装置。
【請求項８】
少なくとも２つの分子を所定の規則に基づいて重畳する分子重畳部と、
前記分子重畳部で重畳された分子に対して、着目する原子或いは基毎に複数の領域に分割する領域分割部と、
前記分子重畳部で重畳された分子をメッシュに分割し、メッシュの各格子点における少なくとも静電ポテンシャルエネルギーと立体ポテンシャルエネルギーとを含むポテンシャルエネルギーを計算するポテンシャルエネルギー計算部と、
前記格子点が属する前記複数の領域ごとに平均化ポテンシャルエネルギーを計算し、前記平均化ポテンシャルエネルギーと反応特性値との相関を、サポートベクタマシンを用いて計算する相関計算部と、
前記相関計算部の計算結果に基づいて、反応特性値を提示する結果提示部とを具備することを特徴とする構造反応特性相関解析装置。
【請求項９】
分子の反応特性値と分子のポテンシャルエネルギーとの相関を計算するためのコンピュータで実行可能な構造反応特性相関解析プログラムにおいて、
少なくとも２つの分子を所定の規則に基づいて重畳する分子重畳手段と、
前記分子重畳手段で重畳された分子に対して、着目する原子或いは基毎に複数の領域に分割する領域分割手段と、
前記分子重畳手段で重畳された分子をメッシュに分割し、メッシュの各格子点における少なくとも静電ポテンシャルエネルギーと立体ポテンシャルエネルギーとを含むポテンシャルエネルギーを計算するポテンシャルエネルギー計算手段と、
前記格子点が属する前記複数の領域ごとに平均化ポテンシャルエネルギーを計算し、前記平均化ポテンシャルエネルギーと反応特性値との相関を計算する相関計算手段と、
前記相関計算手段の計算結果に基づいて、反応特性値を、前記ポテンシャルエネルギーを独立変数とする一次式として提示する結果提示手段とを具備することを特徴とする構造反応特性相関解析プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、分子の反応特性値と所定の対象分子に対する各ポテンシャルエネルギーとの間の相関を解析する構造反応特性相関解析装置および構造反応特性相関解析プログラムに関する。
【背景技術】
【０００２】
研究開発において分子の反応特性の情報を必要とする分野の１つに薬物開発がある。薬物開発では、種々の薬物が生体に及ぼす薬物活性を調べ、より高い活性を有する薬物分子が設計される。薬物活性は、薬物がある強さで生体に対して効力を示す濃度Ｃの逆数の対数値である活性値（ｌｏｇ（１／Ｃ））によって評価される。
【０００３】
薬物活性は、分子の反応特性の１つであり、生体の細胞膜などに存在する受容体分子に対する薬物分子の相互作用ポテンシャルエネルギーによって決定されると考えられている。言い換えれば、相互作用ポテンシャルエネルギーと薬物活性の間に相関があると考えられている。但し、その物理化学的な相互作用がどのような効果をもたらすかは、大規模分子である受容体およびその周辺の環境に依存するため直接説明することはできない。そこで、指定された薬物の分子に類似した構造を有する参照分子群の各分子を空間的に重ね合わせて、それらの各分子の３次元的な相互作用ポテンシャルエネルギーと活性値との相関係数を計算することにより薬物の構造と活性値との関係を解析する３次元定量的構造活性相関解析（３Ｄ−ＱＳＡＲ）という計算手法が開発されている。
【０００４】
この計算手法の１つとして、本発明者の１人により開発されたＶＦＡ（Voronoi Field Analysis）が知られている（例えば、非特許文献１参照）。ＶＦＡは、分子を囲む空間をその分子を構成する各原子を囲むボロノイ（Voronoi）領域に分割して、相互作用ポテンシャルエネルギーをこのボロノイ領域ごとに平均化して取り扱うことを特徴とするものである。一方、それ以前に知られていた３Ｄ−ＱＳＡＲ計算手法であるＣｏＭＦＡ（Comparative Molecular Field Analysis：特許文献１参照）は、空間を１−２オングストローム間隔のグリッドで分割し、その各点での相互作用ポテンシャルエネルギーの値を取り扱う方法である。ＶＦＡは、ＣｏＭＦＡの、そのグリッドの数の増加に伴って上述の相関係数が増加しやすくその計算量が膨大になるという問題、後述する立体ポテンシャルエネルギーの計算においてエネルギーカットオフが必要であることにより解析結果に任意性が生ずる問題、３０年あまりにわたり解析例の膨大な蓄積があるＨａｎｓｈ−Ｆｕｊｉｔａ法（例えば、非特許文献２参照）との直接的な比較対応が困難でありその資産を生かしにくいといった問題などを解決している。
【特許文献１】米国特許第５，３０７，２８７号
【非特許文献１】H. Chuman他, Quant. Struct. -Act. Relat., 17, 313-326 (1998)
【非特許文献２】C. Hansch, T. Fujita, J. Am. Chem.Soc., 86, 1616-1626 (1963)
【発明の開示】
【発明が解決しようとする課題】
【０００５】
しかし、ＶＦＡにも、参照分子群の重ね合わせおよび相関係数決定の際のパラメータ選択など、調整および最適化の操作が多く、それらの操作がユーザの負担となり、ＶＦＡが本来もたらすべき解析結果の精度を十分に引き出せないという問題を有する。
【０００６】
本発明は、上記の課題を解決するためになされたものであって、計算量を少なくすると共に、分子の構造の反応特性に対する相関を解析する際のユーザの負担を軽減してその解析の精度を向上させる構造反応特性相関解析装置および構造反応特性相関解析プログラムを提供することを目的とする。
【課題を解決するための手段】
【０００７】
本発明の局面に係る発明は、少なくとも２つの分子を所定の規則に基づいて重畳する分子重畳部と、前記分子重畳部で重畳された分子に対して、着目する原子或いは基毎に複数の領域に分割する領域分割部と、前記分子重畳部で重畳された分子をメッシュに分割し、メッシュの各格子点における少なくとも静電ポテンシャルエネルギーと立体ポテンシャルエネルギーとを含むポテンシャルエネルギーを計算するポテンシャルエネルギー計算部と、前記格子点が属する前記複数の領域ごとに平均化ポテンシャルエネルギーを計算し、前記平均化ポテンシャルエネルギーと反応特性値との相関を計算する相関計算部と、前記相関計算部の計算結果に基づいて、反応特性値を、前記ポテンシャルエネルギーを独立変数とする一次式として提示する結果提示部とを具備することを特徴とする。なお、本発明は、物（装置）の発明以外に、装置で実現する方法の発明としても成立するし、更には、当該方法を実現するためのプログラムの発明としても成立する。
【発明の効果】
【０００８】
本発明によれば、相関計算に対する計算量が減少すると共に、分子の構造の反応特性に対する相関を解析する際のユーザの負担を軽減してその解析の精度を向上させることができる。
【発明を実施するための最良の形態】
【０００９】
図面を参照して本発明の実施の形態を説明する。図１は、本発明の一実施形態に係る構造反応特性相関解析装置の概略構成を示すブロック図である。なお、図１のブロック図において、矢印は、データの流れを示している。なお、本発明の一実施形態に係る構造反応特性相関解析装置は、汎用のコンピュータを用いて実現しても良いし、専用のコンピュータを用いて実現しても良く、本発明で要求する相関解析を実行可能であれば、その構成はそのような構成であっても良い。  
本発明の一実施形態に係る構造反応特性相関解析装置は、少なくとも２つの分子の重畳を行って反応特性の相関解析を行う構造反応特性相関解析部１０と、種々のデータの入力や制御データの入力を行う入力部２０と、分子データや重畳データなどを記憶するデータ記憶部３０と、反応特性を取得する反応特性値取得部４０と、構造反応特性相関解析部１０の相関解析の結果などを表示する表示部５０とを備えている。  
構造反応特性相関解析部１０は、複数の分子を入力して、所定の方法で重畳する分子重畳部１１と、特異点を中心として重畳後の分子を所定の領域毎に分割する領域分割部１２と、重畳した分子を所定のメッシュに分割して、各格子点におけるポテンシャルエネルギーを計算するポテンシャルエネルギー計算部１３と、分割された領域におけるポテンシャルエネルギーと反応特性との相関を計算する相関計算部１４と、相関計算部１４で算出された結果を提示する結果提示部１５とを備えている。
【００１０】
上記のように構成された構造反応特性相関解析装置の具体的な動作について説明する。
まず、ユーザは入力部２０を介して、データ記憶部３０から少なくとも２つの所望の相互に類似する分子を選択する。ここで、「類似」とは、分子の主要部が同一である分子をいい、例えば、置換基（なお、本明細書では、置換基に限らず、分子として一つの作用を有するものを「基」として取り扱う）を取り除いた場合の分子がまったく同一の分子を云う。なお、以下の説明では、分子は２つ以上でもよいが、説明の便宜上、２つの分子が選択されたものとして説明する。この場合において、例えば、通信回線からダウンロードしたものや、記憶媒体を介して入手した分子が予めデータ記憶部３０に記憶されているものとする。また、データ記憶部３０に記憶された分子は、ユーザによってその一部が置換基で置換されたような分子であっても良いし、置換基を除いた部分や、置換基のみのデータを含んでいても良い。
【００１１】
入力部２０で選択された分子は、構造反応特性相関解析部１０の分子重畳部１１に入力して、次のように２つの分子が重畳される。分子の重畳は、対応する原子を演算によって自動的に求めても良いし、ユーザが表示部５０に表示された２つの分子について、それぞれ対応する原子として、重ね合わせられるべき（または重ね合わせるべきと考えられる）原子の一部或いは全部を指定しても良い。また、自動的に重ね合わせた分子をユーザが対応する部分（例えば原子）を指定して再重畳しても良い。なお、重ねあわせに対して、２つの分子Ａ及び分子Ｂにおいて、それぞれ対応する原子（或いは原子団）について、それらの位置ベクトルについて、分子の内部自由度を考えないものとすると、ｒAi、ｒBi（ｉ＝１、２、３、・・・、ｎ）であるときに、  
Ｆ＝Σ｜ｒAi−ｒBi｜2  
が最小になるときに、最適な重ね合わせが得られるものとする。このように、対応する原子を基にして、２つの分子の原子同士の位置がなるべく一致するように、原子間の結合部分を回転させることによって（これを「配座の変更」と称する）、２つの分子を重畳する。この重ね合わせの様子を図２に示す。まず、図２（ａ）に示すように、選択された異なる２つの分子Ａと分子Ｂとが表示部５０に表示されている。これらの２つの分子Ａ，Ｂについて、その類似した部分を重ね合わせて（すなわち、重ね合わせるのに妥当な原子同士を一致させて）、重畳結果として表示部５０に表示する（図２（ｂ）参照）。なお、この重ね合わせるのに際して、適宜、それぞれの分子の原子の位置がより一致するように重畳するために、分子を回転させたり、結合部を回転させたりすることが好ましい。なお、ある化合物に置換基を重畳させる場合には、置換基が結合する原子同士を重ね合わせることにより、重畳すればよい。例えば、ベンゼン環に、メチル基を重畳させる場合には、メチル基の炭素がベンゼン環の１つの炭素に結合するので、メチル基の炭素原子にダミーの原子を付加しておき、その原子とベンゼン環の１つの炭素と重ね合わせるようにすれば良い。
【００１２】
上記のように重畳された分子重畳部１１で重畳された分子は、重畳の都度、領域分割部１２と、ポテンシャルエネルギー計算部１３とに出力されると共に、データ記憶部３０に出力されてその都度記憶される。なお、この場合において、データ記憶部３０は、重畳後のデータと共に、重畳前のそれぞれの分子のデータも記憶しておくことが好ましい。ここで、分子が占める空間は計算上閉空間であることが必要である。そのため、この閉空間を次のように定め、その閉空間内で計算を行うものとする。この閉空間は、各原子のファンデルワールス半径をつなぎ合わせた空間よりも数オングストローム大きい空間で分子全体を覆うような空間とする。
【００１３】
領域分割部１２は、詳細は後述する相関計算部１４で相関計算を行う際の変数を少なくするために、分子中の原子或いは基に注目して分割を行う。具体的には、ボロノイ分割を行うものとし、この場合において、分子中の原子を基本として原子ごとに領域を分割する場合には、各原子の作用の及ぶ範囲が隣接する原子の中間点までと定義して、その中間点を結ぶ平面で領域を分割する。図３は、分割後の分割領域を示す図である。図３に示すように、各原子の作用の及ぶ範囲としてボロノイ領域が定められていることがわかる。なお、領域分割の基本として原子を中心とした分割を行うことが最も適切であるが、置換基を一単位として、その置換基を一原子として捉えることによってボロノイ分割を行っても良い。
【００１４】
ポテンシャルエネルギー計算部１３は、まず、分子全体を、領域分割部１２で分割された領域より小さいメッシュである３次元メッシュに分割する。このメッシュの間隔は、計算機の能力と精度のとの兼ね合いにより定まるものであるので、特に定めないが、前述の非特許文献１に記載のように０．３オングストローム程度でも良い。そして、３次元メッシュの全ての格子点についてポテンシャルエネルギー（静電ポテンシャルエネルギーと立体ポテンシャルエネルギー）を求める。ここで、静電ポテンシャルエネルギーは、Ｈ＋（水素原子核）を格子点に置いたときの分子との静電ポテンシャルエネルギーであり、立体ポテンシャルエネルギーは、格子点がファンデルワールス表面上に存在するか否かを示すものである。なお、このポテンシャルエネルギーは、分子の類似性を判定するための指標として使用しても良い。「類似分子」とは、図４（ａ）に示すような構造を有する分子において、ＸｎとＲ１、Ｒ２以外は共通の構造を有する分子において、ＸｎとＲ１、Ｒ２とを除いたものが同一であるような構造を有する分子をいう。例えば、図４（ａ）の分子では、図４（ｂ）に示すように、Ｘｎは、Ｈ、２−Ｃｌ、・・・、Ｒ１は、Ｈ、Ｍｅ又はＥｔ、Ｒ２はＨ又はＭｅと置換可能であることを示している。
【００１５】
上記の領域分割部１２での分割結果と、ポテンシャルエネルギー計算部１３で計算されたポテンシャルエネルギーは、相関計算部１４に出力される。相関計算部１４には、反応特性値取得部４０で取得された反応特性値も出力される。反応特性値取得部４０で取得される反応特性値は、分子の活性値（すなわち、反応のし易さ）を示すものであって、例えば、ある薬物が効き始める濃度をＣとした場合に、この逆数の対数値（＝ｌｏｇ（１／Ｃ））で表すことができる。従って、以下の説明において、反応特性値を活性値と証する場合もある。
【００１６】
相関計算部１４は、図５に示すようなフローチャートに従って、相関計算を行う。ポテンシャルエネルギー計算部１３でメッシュの格子点が領域分割部のどの領域に属するかを判定する（ステップＳ１）。例えば、領域分割部１２で分子が２つの領域Ａ及び領域Ｂに分割された場合に、メッシュの格子点がそれぞれどの領域に属するかを決定する。そして、領域毎に平均化ポテンシャルエネルギーを求める（ステップＳ２）。ここで、平均化ポテンシャルエネルギーは、算術平均を用いる。なお、この算術平均として、単純平均を求めても良いし、原子に近いほど重みを付ける加重平均としても良いし、その他、様々な平均方法が考えられる。そして、特性値と平均化ポテンシャルエネルギーとの相関をサポートベクタマシン（ＳＶＭ）による近似法で求める（ステップＳ３）。このサポートベクタマシンは、数多くの応用が考えられており、テキストのカテゴリ化や、画像認識や、文字認識などにおけるそれらの分類に応用されている。サポートベクタマシンでは、空間中に配置された２つの集合を最適に分離する境界を、その集合に属する点の確率分布モデルを考えることなく求めることができる。ここで、関数近似に対応した損失関数を  
ｍａｘ｛０，｜ｙ−ｆ（ｘ）｜−ε｝  
というε不感応関数で定義することにより、識別問題の場合と似たような凸２次計画問題に帰着することが知られている。このアルゴリズムを用いて、近似式が得られる。この場合において、サポートベクタマシンを用いることにより、大域的な最適化を行うことができ、他の最適化法のような局所最適解を最適解として出力することはない。この結果は、例えば、標的分子中のそれぞれの原子或いは分子（置換基）に対する活性値（アクティビティ：Activity）、すなわち反応特性値を示す係数を持つ線型方程式（一次式）として与えられる。
【００１７】
上記の結果は結果提示部１５によって、例えば、一次近似式として提示され、  
Activity＝ａ1ｘ1＋ａ2ｘ2＋ａ3ｘ3＋・・・
＋ｂ1ｙ1＋ｂ2ｙ2＋ｂ3ｙ3＋・・・＋ｃ  
といった形式で、表示部５０に表示される。ここで、例えば、ｘは静電ポテンシャルエネルギーによるものであり、ｙは立体ポテンシャルエネルギーによるものである。また、ｃは定数である。この一次式により、どの原子或いは分子（置換基）を他の原子或いは分子（置換基）に置き換えたときに全体の活性値がどの程度大きく又は小さくなるかが、当該原子或いは分子の係数ａ、ｂにより決定されるので、直感的に、どの原子或いは分子を置換すれば活性値が高くなるのかの目安を得ることができる。このため、分子の一部を他の原子や分子に置換して、新たな化合物を生成した場合においても、その化合物の活性値が定量的に予測できるので、新たな化合物の生成を検討するのに有効である。
【００１８】
そして、提示された結果に基づいて、所望の部分における原子或いは分子の置換（すなわち、再重畳）を行い、上記のステップを繰り返すことにより、化合物の活性値の推定が可能になる。
【００１９】
上記のような活性値の推定により、所望の反応特性値を有すると思われる化合物が決定した時点で、通常実験により正確な反応特性値が求められる。しかし、本発明の実施形態では、相関計算をサポートベクタマシンを利用して求めていることから、カーネル関数を変更して非線形解（すなわち反応特性値）を求めるようにすれば、高精度な解を求めることができる。すなわち、本発明では、サポートベクタマシンを利用して最適解の計算を行っているので、異なるアルゴリズムを用いることなく、線形方程式と非線形解とが得られる。従って、本発明の一実施形態において、化合物の実際の反応特性値を実験で求める前に、相関計算部１４でサポートベクタマシンを使用して非線形のより高精度な解を求めるようにすることにより更に正確な反応特性値（すなわち活性値）の予測値を予め得ることができるので、実験の手間を減少させることができる。
【００２０】
また、上記では、サポートベクタマシンによって得られた一次近似式を用いてさまざまな化合物の反応特性値を推定した後に、更にサポートベクタマシンによって非線形のより精度の高い反応特性値を計算するようにしたが、これに限らず、化合物を入力部２０から入力（或いは入力部２０で選択）したら、その化合物について、サポートベクタマシンを用いて、そのまま非線形の最適化計算によって精度の高い反応特性値の予測値を計算するようにしても良い。
【００２１】
上記のように、本発明によれば、ボロノイ分割により分子を複数の領域に分割して、その領域内の平均化ポテンシャルエネルギーと反応特性値との相関を計算するようにしたので、相関計算に係る変数の数を格段に減らすことができる。従って、計算の速度が速くなり、計算機にかかる負担が軽減される。また、ポテンシャルエネルギーを独立変数とした一次式で相関を定義するようにしたので、一次式の係数を参照して最も反応特性値に影響のありそうな因子（例えば、置換基や分子、或いは原子）を容易に推定することができる。また、分子の類似性をポテンシャルエネルギーに基づいて提示するようにしたので、ユーザによる分子の重ね合わせをより高精度で行うことができる。
【００２２】
更に、分子の重ね合わせの状態をデータ記憶部に記憶するようにしたので、重ね合わせの様子が容易に再現できると共に、一部を変更して再重畳するのも容易にできる。また、これにより、原子同士の対応や結合周りの回転などの設定パラメータ数が多い大きな分子の重ね合わせに対して、ユーザは分子の様々な重ね合わせ状態を比較検討しながら容易に適切な重ね合わせ状態の決定を行うことができる。更に、ポテンシャルエネルギーの値の分布を色相や、明暗により明示的に表示するようにしたので、ユーザは、各（ボロノイ）分割領域毎の平均値を計算する以前にそのポテンシャルエネルギーの空間分布を確認しながら分子の重ね合わせをより適切に行うことができるようになる。
【００２３】
また、相関計算をサポートベクタマシンを使用して実行させることにより、変数（パラメータ）が多くなっても高速な計算を実行させることができる。また、大域的な最適解を提示することができ、局所最小解（Local Minimum）の問題は生じない。更に、サポートベクタマシンによる解は、他の方法で求めた解よりも一般性（精度）が高いため、新たなデータを用意して、そのデータに対する反応特性相関を予測する場合に、他の方法よりも妥当な予測が可能である。また、サポートベクタマシンのみで、一次式による表現と非線形のより高い精度の予測計算が、カーネル関数の変更のみできるので、ユーザの要求に応じたさまざまな計算が１台の計算機で行うことが可能になる。
【００２４】
本発明は、上記各実施の形態に限ることなく、その他、実施段階ではその要旨を逸脱しない範囲で種々の変形を実施し得ることが可能である。さらに、上記各実施形態には、種々の段階の発明が含まれており、開示される複数の構成要件における適宜な組合せにより種々の発明が抽出され得る。
【００２５】
また、例えば各実施形態に示される全構成要件から幾つかの構成要件が削除されても、発明が解決しようとする課題の欄で述べた課題が解決でき、発明の効果で述べられている効果が得られる場合には、この構成要件が削除された構成が発明として抽出され得る。
【図面の簡単な説明】
【００２６】
【図１】本発明の一実施形態に係る構造反応特性相関解析装置の概略構成を示すブロック図。
【図２】２つの異なる分子を重ね合わせる様子を示す図。
【図３】ボロノイ分割による分割後の分割領域を示す図。
【図４】類似分子を説明するための図。
【図５】相関計算を行うためのフローチャート。
【符号の説明】
【００２７】
１０…構造反応特性相関解析部
１１…分子重畳部
１２…領域分割部
１３…ポテンシャルエネルギー計算部
１４…相関計算部
１５…結果提示部
２０…入力部
３０…データ記憶部
４０…反応特性値取得部
５０…表示部
【図１】
2007106708.tif 000003
【図２】
2007106708.tif 000004
【図３】
2007106708.tif 000005
【図４】
2007106708.tif 000006
【図５】
2007106708.tif 000007