(57)【要約】【課題】  複数の化合物について第１の特性と第２の特性とを関連付けて効率的に表示できる化合物群表示装置であって、第１の特性として４つ以上の特性をも同時に取り扱うことのできる化合物群表示装置を提供する。【解決手段】  化合物群表示装置１は、化合物を第１の特性情報に基づいて複数の階層的なクラスターに分類する階層分類部１３と、これによって分類された階層的なクラスターに基づいた再帰的な入れ子構造となるように、化合物の２次元平面内における配置を決定する化合物配置部１４と、これによって決定された配置で、第２の特性情報に基づいた色によって化合物が描画された画像データを生成する画像データ生成部１５と、これによって生成された画像データに基づいた画像を表示する表示部１６とを備えている。【選択図】  図１
2006318048.tif 000002
【特許請求の範囲】
【請求項１】
複数の化合物について、当該各化合物の第１の特性を表す第１の特性情報と第２の特性を表す第２の特性情報とに基づいて、第１の特性と第２の特性との関係を表示する化合物群表示装置であって、
上記化合物を、第１の特性情報に基づいて複数の階層的なクラスターに分類する階層分類手段と、
上記階層分類手段によって分類された階層的なクラスターに基づいた再帰的な入れ子構造となるように、上記化合物のアイコンの２次元平面内における配置を決定する化合物配置手段と、
上記化合物配置手段によって決定された配置で、第２の特性情報に基づいた色、模様、及び／又は形状によって上記化合物のアイコンが描画された画像データを生成する画像データ生成手段と、
上記画像データ生成手段によって生成された画像データに基づいた画像を表示する表示部とを備えていることを特徴とする化合物群表示装置。
【請求項２】
第１の特性情報が、ｎ個（ただし、ｎは自然数）の量的変数からなるものであり、
上記階層分類手段が、上記化合物を、各量的変数を座標軸とするｎ次元空間における点として扱い、各化合物間又は各クラスター間の距離に基づいて、上記化合物を分類することを特徴とする請求項１に記載の化合物群表示装置。
【請求項３】
上記階層分類手段が、第１の特性情報を説明変数、第２の特性情報を目的変数として、決定木を作成することによって上記化合物を分類することを特徴とする請求項１に記載の化合物群表示装置。
【請求項４】
第２の特性情報に基づいて上記化合物を非階層的なクラスターに分類する非階層分類手段をさらに備え、
上記階層分類手段が、第１の特性情報を説明変数、上記非階層分類手段によって形成されたクラスターに対応するカテゴリ変数を目的変数として、決定木を作成することによって上記化合物を分類することを特徴とする請求項１に記載の化合物群表示装置。
【請求項５】
第２の特性情報に基づいて上記化合物を非階層的なクラスターに分類する非階層分類手段をさらに備え、
上記画像データ生成手段が、上記非階層分類手段によって形成された各クラスターに対応する色、模様、及び／又は形状によって上記化合物のアイコンが描画された画像データを生成することを特徴とする請求項１に記載の化合物群表示装置。
【請求項６】
第２の特性情報が複数の量的変数からなるものであり、
上記非階層分類手段が、上記第２の特性情報をパターンベクトルとした自己組織化マップ法によって、上記化合物を非階層的なクラスターに分類することを特徴とする請求項４又は５に記載の化合物群表示装置。
【請求項７】
第１の特性情報が、
(a)化合物の分子構造の構造フラグメントもしくは構造トポロジーを表す１つ以上の分子記述子；
(b)計算もしくは実験によって求められた物理化学的性質に対応する１つ以上の分子記述子、
の少なくとも何れか一方を含んでいることを特徴とする請求項１から６の何れか１項に記載の化合物群表示装置。
【請求項８】
第２の特性情報が、上記化合物の生物活性、毒性、物理化学的性質、薬物動態学的性質の少なくとも何れかのパラメータであることを特徴とする請求項１から７の何れか１項に記載の化合物群表示装置。
【請求項９】
第２の特性情報がカテゴリ変数からなることを特徴とする請求項１から８の何れか１項に記載の化合物群表示装置。
【請求項１０】
請求項１から９の何れか１項に記載の各手段として、コンピュータを動作させるためのプログラム。
【請求項１１】
請求項１０に記載のプログラムが記録されたコンピュータ読み取り可能な記録媒体。
【請求項１２】
複数の化合物について、当該各化合物の第１の特性を表す第１の特性情報と第２の特性を表す第２の特性情報とに基づいて、第１の特性と第２の特性との関係を表示する化合物群表示装置における化合物群表示方法であって、
上記化合物群表示装置が、階層分類手段、化合物配置手段、画像データ生成手段、及び表示部を備え、
上記階層分類手段が第１の特性情報に基づいて複数の階層的なクラスターに分類する階層分類工程と、
上記化合物配置手段が上記階層分類工程において分類された階層的なクラスターに基づいた再帰的な入れ子構造となるように、上記化合物のアイコンの２次元平面内における配置を決定する化合物配置工程と、
上記画像データ生成手段が上記化合物配置工程において決定された配置で、第２の特性情報に基づいた色、模様、及び／又は形状によって上記化合物のアイコンが描画された画像データを生成する画像データ生成工程と、
上記表示部が上記画像データ生成工程において生成された画像データに基づいた画像を表示する表示工程とを備えていることを特徴とする化合物群表示方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、複数の化合物を分類して表示する化合物群表示装置及び化合物群表示方法と、この装置を実現するためのプログラム、さらには、このプログラムを記録したコンピュータ読み取り可能な記録媒体に関するものである。
【背景技術】
【０００２】
医薬品、農薬、化粧品、及びその他の化学工業の分野では、活性や物性等の機能特性が優れた化合物を探索する研究が常に行われている。最近では、例えば、コンビナトリアルケミストリー技術とハイスループットスクリーニング法との組み合わせが頻繁に利用されている。これにより、多数の化合物を同時にかつ高速に合成し、評価することが可能となり、化合物の機能特性に関する大規模なスクリーニングデータが蓄積されるようになった。
【０００３】
このようにして蓄積されたスクリーニングデータを基に、さらに機能特性の優れた新規化合物を合成展開するためには、得られたスクリーニングデータを化合物の構造と対応付けて整理することが不可欠である。
【０００４】
化合物の構造とスクリーニングデータとの関係を整理するためには、まず、化合物の構造を数量的に表現する手法が必要である。この手法として、これまでに数多くの分子記述子が考案されている。その代表的なものとしては、構造キー（Daylight Chemical Infomation Systems社、MDL Information Systems社）、Molconn-Zディスクリプタ（Hall Associate Consulting社）、ClogP（Biobyte社）、ACD/LogD（ACD Labs社）等が挙げられる。
【０００５】
一方、これらの分子記述子を用いてスクリーニングデータを管理する方法についても、これまで多くの方法が提案されている。その方法は以下の二つに大別される。
【０００６】
一つは、スクリーニングデータを目的変数、分子記述子を説明変数として多変量解析を行い、それらの関係を定量的に表現するモデル式として管理する方法である（例えば、非特許文献１）。これにより、未知化合物についても、化学構造から分子記述子を計算し、これをモデル式に適用することによって機能特性を予測することが可能となる。
【０００７】
もう一つの方法は、各化合物の分子記述子とスクリーニングデータ（機能特性）とを対応付けたデータベースを構築するする方法である（例えば、非特許文献２）。構造的に類似すれば機能特性も類似することが予想されるため、未知化合物と構造的に類似する化合物をデータベースから検索すれば、検索した化合物のスクリーニングデータに基づいて未知化合物の機能特性を予測することが可能となる。この場合、類似性／非類似性の評価は、分子記述子に基づいて行われる。
【０００８】
上記の各方法は、何れも統計学又は推計学的なアプローチに基づくものであり、化合物における分子記述子と機能特性との関係を経験的に予測するものである。従って、これらの方法では、機構論に基づく理論的な予測に比して計算結果の妥当性が問われることが多く、妥当性を客観的に評価しうる表現方法が望まれている。これを解決する極めて有効な方法がデータの可視化である。データを可視化すれば、人間が経験によって培った高度な分析能力を活用して適切な判定を行うことができる。また、可視化技術により、大規模な情報リソースを共有することが容易になるため、プロジェクト研究において意思の統一を図り、プロジェクトを効率よく推進することも可能になる。
【０００９】
従来、化合物の機能特性を構造情報と関連付けて可視化する方法として、(a)単一の分子記述子、又は(b)複数の分子記述子に基づく総合特性値、を軸とする３次元以下の特性空間に化合物をマッピングする方法が採用されてきた（例えば、非特許文献３、特許文献１）。これにより、データベースに収載される化合物群について、機能特性と構造情報とを対応付けて表示することができる。そして、特性の評価を行いたい未知化合物があれば、構造情報に基づいて未知化合物を特性空間上にプロットすることによって、データベースに蓄積された化合物の機能特性から未知化合物の特性を視覚的に評価することができる。
【特許文献１】特公２００２−５３１８９４（平成１４年(2002)９月２４日公開）
【非特許文献１】C. Hansch et al., Chem-bioinformatics: comparative QSAR at the interface between chemistry and biology. Chem. Rev., 102: 783-812, 2002
【非特許文献２】R.P. Sheridan and S. K. Kearsley, Why do we need so many chemical similarity search methods? Drug Discovery Today, 7: 903-911, 2002
【非特許文献３】Y. Takahashi et al., MolSpace: a computer desktop tool for visualization of massive molecular data. J. Mol. Graph Model., 21: 333-339, 2003
【発明の開示】
【発明が解決しようとする課題】
【００１０】
しかしながら、上記従来の方法には以下の問題がある。
【００１１】
まず、従来の方法では、特性空間における各座標軸が間隔尺度を表すものであるため、プロットされた化合物の配置が離散的になったり、偏在したりする。これにより、ディスプレイの表示領域の一部に化合物の存在しない空白部分ができたり、それとは逆に、一部に化合物が密集して重複部分ができたりするおそれがある。このように、従来の方法では、化合物を特性空間内にプロットする際に、ディスプレイの表示領域を有効に利用できない。
【００１２】
また、ディスプレイの表示領域では、３次元の特性空間を表示するのが限界であり、４次元以上の特性空間を表示することはできない。これは、化合物の構造情報として、４つ以上の構造情報を同時に取り扱えないことを意味する。よって、構造情報を何らかの形で３つ以下に減少させる必要があり、このときに、化合物の構造情報の一部が失われてしまう。このように、従来の方法では、化合物について、多数の構造情報からなる構造特性と機能特性との関係を包括的に評価することができない。
【００１３】
本発明は上記課題に鑑みてなされたものであり、その目的は、複数の化合物について、例えば構造特性及び機能特性のような、第１の特性と第２の特性とを関連付けて効率的に表示できる化合物群表示装置であって、第１の特性として多数の特性情報をも同時に取り扱うことのできる化合物群表示装置を提供することにある。
【課題を解決するための手段】
【００１４】
上記課題を解決するために、本発明に係る化合物群表示装置は、複数の化合物について、当該各化合物の第１の特性を表す第１の特性情報と第２の特性を表す第２の特性情報とに基づいて、第１の特性と第２の特性との関係を表示する化合物群表示装置であって、上記化合物を、第１の特性情報に基づいて複数の階層的なクラスターに分類する階層分類手段と、上記階層分類手段によって分類された階層的なクラスターに基づいた再帰的な入れ子構造となるように、上記化合物のアイコンの２次元平面内における配置を決定する化合物配置手段と、上記化合物配置手段によって決定された配置で、第２の特性情報に基づいた色、模様、及び／又は形状によって上記化合物のアイコンが描画された画像データを生成する画像データ生成手段と、上記画像データ生成手段によって生成された画像データに基づいた画像を表示する表示部とを備えていることを特徴とする。
【００１５】
なお、第１の特性としては、例えば、化合物の構造特性等が挙げられ、第２の特性としては、例えば、化合物の機能特性等が挙げられる。
【００１６】
上記構成によれば、階層分類手段が化合物を第１の特性情報に基づいて階層的なクラスターに分類する。そして、化合物配置手段及び画像データ生成手段が、分類したクラスターに基づいた入れ子構造ととなるように化合物を２次元平面内に表示する。これにより、各化合物は、アイコンの含まれる入れ子構造によって第１の特性が表現され、アイコンの色、模様、及び／又は形状によって第２の特性が表現される。以上のように、本発明の化合物群表示装置は、化合物群を、第１の特性と第２の特性とを関連付けて表示することができる。なお、上記の色には、モノクロ及びグレースケールも含まれるものとする。
【００１７】
ここで、階層的なクラスターに分類するにあたって、第１の特性情報に含まれる変数の個数は複数であってもよい。すなわち、第１の特性情報が多次元ベクトルからなるものであっても、化合物を階層的なクラスターに分類することができる。よって、化合物群表示装置は、第１の特性として多数の特性情報をも同時に取り扱うことができる。
【００１８】
また、化合物は、第１の特性情報を座標軸とする特性空間における点によって表現されるのではなく、第１の特性情報に基づいた再帰的な入れ子構造におけるアイコンによって表現される。ここで、各アイコンの位置は、第１の特性情報の数値に基づいて一点に決定されるわけではなく、類似する化合物のアイコン同士は近隣に存在する、という風に相対的な位置関係が大まかに決定される。よって、アイコンの配置が離散的になったり、偏在したりしないように、なるべく一様な間隔になるように配置することができる。このように、ディスプレイの表示領域を有効に利用できる。
【００１９】
また、第１の特性情報は、ｎ個（ただし、ｎは自然数）の量的変数からなるものであり、上記階層分類手段が、上記化合物を、各量的変数を座標軸とするｎ次元空間における点として扱い、各化合物間又は各クラスター間の距離に基づいて、上記化合物を分類してもよい。
【００２０】
上記構成によれば、化合物の座標やクラスターの重心などの座標から算出した距離を非類似度の尺度として、化合物を分類することができる。
【００２１】
また、上記階層分類手段は、第１の特性情報を説明変数、第２の特性情報を目的変数として、決定木を作成することによって上記化合物を分類してもよい。
【００２２】
上記構成によれば、第１の特性情報に基づいて決定木を作成することによって分類されたクラスター内の化合物同士は、第２の特性情報が一様になる傾向にある。よって、第１の特性と第２の特性との対応関係が明確になり、未知の化合物がどのような第２の特性を有するかを、その化合物の第１の特性情報に基づいて容易に推定できるようになる。なお、決定木には回帰木も含まれるものとする。
【００２３】
また、上記化合物群表示装置は、第２の特性情報に基づいて上記化合物を非階層的なクラスターに分類する非階層分類手段をさらに備え、上記階層分類手段が、第１の特性情報を説明変数、上記非階層分類手段によって形成されたクラスターに対応するカテゴリ変数を目的変数として、決定木を作成することによって上記化合物を分類してもよい。
【００２４】
上記構成によれば、第２の特性情報が複数の変数からなるものであっても、非階層分類手段が化合物を第２の特性情報に基づいて非階層的なクラスターに分類するので、化合物が有する第２の特性情報を、クラスター数と同数の値を取り得るカテゴリ変数として扱うことができるようになる。これにより、決定木の各クラスターにおける化合物の第２の特性を明確にすることが可能になる。
【００２５】
また、上記化合物群表示装置は、第２の特性情報に基づいて上記化合物を非階層的なクラスターに分類する非階層分類手段をさらに備え、上記画像データ生成手段が、上記非階層分類手段によって形成された各クラスターに対応する色、模様、及び／又は形状によって上記化合物のアイコンが描画された画像データを生成するものであってもよい。
【００２６】
上記構成によれば、第２の特性情報が複数の変数からなるものであっても、非階層分類手段が化合物を第２の特性情報に基づいて非階層的なクラスターに分類するので、各化合物が有する第２の特性を、分類したクラスターに基づいた色、模様、及び／又は形状によって表現することができる。
【００２７】
また、第２の特性情報が複数の量的変数からなるものであり、上記非階層分類手段が、上記第２の特性情報をパターンベクトルとした自己組織化マップ法によって、上記化合物を非階層的なクラスターに分類してもよい。
【００２８】
自己組織化マップ法によれば、多次元のベクトルにより表されたデータを、その特徴を残し他のデータとの相互関係を保ったまま、２次元のマップに写像することが出来る。この２次元のマップ上では、類似する多次元データを有するノードが近接して配置され、データの関係が視覚的に容易に理解できる。
【００２９】
また、第１の特性情報が、(a)化合物の分子構造の構造フラグメントもしくは構造トポロジーを表す１つ以上の分子記述子；(b)計算もしくは実験によって求められた物理化学的性質に対応する１つ以上の分子記述子、の少なくとも何れか一方を含んでいてもよい。
【００３０】
また、第２の特性情報が、上記化合物の生物活性、毒性、物理化学的性質、薬物動態学的性質の少なくとも何れかのパラメータであってもよい。
【００３１】
また、第２の特性情報がカテゴリ変数からなるものであってもよい。
【００３２】
ところで、上記化合物群表示装置における各手段は、ハードウェアで実現してもよいし、プログラムをコンピュータに実行させることによって実現してもよい。具体的には、本発明に係るプログラムは、上記各手段としてコンピュータを動作させるプログラムであり、本発明に係る記録媒体は、当該プログラムが記録されたものである。
【００３３】
これらのプログラムがコンピュータによって実行されると、当該コンピュータは、上記化合物群表示装置の各手段として動作する。したがって、上記化合物群表示装置と同様の効果が得られる
上記課題を解決するために、本発明に係る化合物群表示方法は、複数の化合物について、当該各化合物の第１の特性を表す第１の特性情報と第２の特性を表す第２の特性情報とに基づいて、第１の特性と第２の特性との関係を表示する化合物群表示装置における化合物群表示方法であって、上記化合物群表示装置が、階層分類手段、化合物配置手段、画像データ生成手段、及び表示部を備え、上記階層分類手段が第１の特性情報に基づいて複数の階層的なクラスターに分類する階層分類工程と、上記化合物配置手段が上記階層分類工程において分類された階層的なクラスターに基づいた再帰的な入れ子構造となるように、上記化合物のアイコンの２次元平面内における配置を決定する化合物配置工程と、上記画像データ生成手段が上記化合物配置工程において決定された配置で、第２の特性情報に基づいた色、模様、及び／又は形状によって上記化合物のアイコンが描画された画像データを生成する画像データ生成工程と、上記表示部が上記画像データ生成工程において生成された画像データに基づいた画像を表示する表示工程とを備えていることを特徴とする。
【００３４】
なお、第１の特性としては、例えば、化合物の構造特性等が挙げられ、第２の特性としては、例えば、化合物の機能特性等が挙げられる。
【００３５】
上記構成によれば、階層分類工程において化合物を第１の特性情報に基づいて階層的なクラスターに分類する。そして、化合物配置工程及び画像データ生成工程において、分類したクラスターに基づいた入れ子構造ととなるように化合物を２次元平面内に表示する。これにより、各化合物は、アイコンの含まれる入れ子構造によって第１の特性が表現され、アイコンの色、模様、及び／又は形状によって第２の特性が表現される。以上のように、本発明の化合物群表示方法は、化合物群を、第１の特性と第２の特性とを関連付けて表示することができる。なお、上記の色には、モノクロ及びグレースケールも含まれるものとする。
【００３６】
ここで、階層的なクラスターに分類するにあたって、第１の特性情報に含まれる変数の個数は複数であってもよい。すなわち、第１の特性情報が多次元ベクトルからなるものであっても、化合物を階層的なクラスターに分類することができる。よって、本方法は、第１の特性として４つ以上の特性をも同時に取り扱うことができる。
【００３７】
また、化合物は、第１の特性情報を座標軸とする特性空間における点によって表現されるのではなく、第１の特性情報に基づいた再帰的な入れ子構造におけるアイコンによって表現される。ここで、各アイコンの位置は、第１の特性情報の数値に基づいて一点に決定されるわけではなく、類似する化合物のアイコン同士は近隣に存在する、という風に相対的な位置関係が大まかに決定される。よって、アイコンの配置が離散的になったり、偏在したりしないように、なるべく一様な間隔になるように配置することができる。このように、ディスプレイの表示領域を有効に利用できる。
【発明の効果】
【００３８】
本発明に係る化合物群表示装置は、以上のように、化合物を第１の特性情報に基づいて複数の階層的なクラスターに分類する階層分類手段と、階層分類手段によって分類された階層的なクラスターに基づいた再帰的な入れ子構造となるように、化合物のアイコンの２次元平面内における配置を決定する化合物配置手段と、化合物配置手段によって決定された配置で、第２の特性情報に基づいた色、模様、及び／又は形状によって化合物のアイコンが描画された画像データを生成する画像データ生成手段と、画像データ生成手段によって生成された画像データに基づいた画像を表示する表示部とを備えた構成となっている。
【００３９】
また、本発明に係る化合物群表示方法は、階層分類手段が第１の特性情報に基づいて複数の階層的なクラスターに分類する階層分類工程と、化合物配置手段が階層分類工程において分類された階層的なクラスターに基づいた再帰的な入れ子構造となるように、化合物のアイコンの２次元平面内における配置を決定する化合物配置工程と、画像データ生成手段が化合物配置工程において決定された配置で、第２の特性情報に基づいた色、模様、及び／又は形状によって上記化合物のアイコンが描画された画像データを生成する画像データ生成工程と、表示部が画像データ生成工程において生成された画像データに基づいた画像を表示する表示工程とを備えた構成となっている。
【００４０】
よって、上述したように、複数の化合物について第１の特性と第２の特性とを関連付けて効率的に表示でき、さらに、第１の特性として多数の特性情報をも同時に取り扱うこともできるという効果を奏する。
【発明を実施するための最良の形態】
【００４１】
〔実施形態１〕
本発明の一実施形態について図１から図５に基づいて説明すると以下の通りである。本実施形態の化合物群表示装置１は、複数の化合物について、第１の特性と第２の特性とを関連付けて表示するものである。第１の特性の種類としては特に限定されるものではないが、本実施形態では、一例として構造特性が用いられ、これを表す構造情報として分子記述子が用いられている。なお、構造情報は、単一の分子記述子であってもよいし、ｎ個の分子記述子、すなわち、各分子記述子の値を成分とするｎ次元ベクトルとして表現されていてもよい。また、分子記述子としては、化合物の分子構造の構造フラグメントもしくは構造トポロジーを表すものであってもよいし、計算もしくは実験によって求められた物理化学的性質に対応するものであってもよいし、これらの組み合わせであってもよい。
【００４２】
また、第２の特性についても、第１の特性情報と同様に特に限定されるものではなく、例えば、化合物の生物活性、毒性、物理化学的性質、又は薬物動態学的性質等の機能特性等を用いることができ、これを表す機能特性情報として各種パラメータを用いることができる。また、これらのパラメータを適宜組み合わせてもよい。本明細書では、これらのパラメータについて以下「機能特性情報」という。
【００４３】
なお、本実施形態において、機能特性情報は、１次元ベクトルで表される量的変数、カテゴリ変数の何れであってもよい。量的変数とは、例えば、化合物の活性値等の間隔尺度で測られる変数であり、離散変数などであってもよく、また、カテゴリ変数とは、例えば、化合物の活性の種類（Ａ活性、Ｂ活性、又はＣ活性）等の複数の値を取りうる変数をいう。
【００４４】
このような構成により、本実施形態の化合物群表示装置は、化合物の構造と機能特性とを関連付けて表示することができる。よって、ユーザは、未知の化合物の機能特性を知りたい場合、構造情報から未知の化合物がどこに表示されるかを知ることによって、構造が類似する化合物の機能特性情報から未知の化合物の機能特性を推測することができる。また、表示されている既知の化合物の分布から、構造と機能特性とがどのように関係しているかを、視覚的に判断することもできる。
【００４５】
図１は、本実施形態の化合物群表示装置１の機能ブロック図である。化合物群表示装置１は、図１に示すように、記述子入力部１１、機能特性入力部１２、階層分類部（階層分類手段）１３、化合物配置部（化合物配置手段）１４、画像データ生成部（画像データ生成手段）１５、及び表示部１６を備えている。
【００４６】
記述子入力部（第１の入力部）１１は、ユーザ又は他の装置等から、表示する化合物の分子記述子（第１の特性情報）を取得するためのものである。記述子入力部１１のハードウェア構成としては、例えばキーボード、マウス、タッチパネル、又はタブレット等のようにユーザからの入力を受け付けるものであってもよいし、各種入出力インターフェース等のように外部記憶装置や他のコンピュータからの入力を受け付けるものであってもよい。
【００４７】
また、機能特性入力部（第２の入力部）１２は、ユーザ又は他の装置等から、表示する化合物の機能特性情報（第２の特性情報）を取得するためのものである。機能特性入力部１２のハードウェア構成については、記述子入力部１１と同様である。
【００４８】
階層分類部１３は、化合物を分子記述子（第１の特性情報）に基づいて複数の階層的なクラスターに分類するためのものである。より具体的には、分子記述子に基づいて、各種アルゴリズムによって化合物を階層的に分類する。
【００４９】
化合物配置部１４は、階層分類部１３によって分類された階層的なクラスターに基づいた再帰的な入れ子構造となるように、化合物の２次元平面内における配置を決定する。より具体的には、化合物配置部１４は、階層的に分類された化合物を葉ノードとして、子ノードが親ノードに順次含まれた入れ子構造となるように、各化合物を２次元平面内に配置する。なお、このように配置するためのアルゴリズムとしては、例えば、平安京ビュー（「平安京ビュー  〜  階層型データを碁盤上に配置する視覚化手法」、伊藤貴之、小山田耕二、可視化情報学会第9回ビジュアリゼーションカンファレンス抄録、2003年発行）を用いることができる。平安京ビューでは、まず階層型データ中の葉ノード群を格子状に配列する。続いてこれらの葉ノードの親にあたる枝ノードを画面上に効率よく配置することで、限られた画面空間上に大量の情報を表現する。
【００５０】
画像データ生成部１５は、化合物配置部１４によって決定された配置で、第２の特性情報に基づいた色によって化合物のアイコンが描画された画像データを生成する。より具体的には、画像データ生成部１５は、化合物配置部１４によって決定された位置に化合物が配置された画像データを生成する。このとき、葉ノードに相当する各化合物は、例えば、円形、正方形、又は点等の簡単な図形で表現してもよいし、化合物の種類が分かるようなシンボルによって表現してもよい。以下、葉ノードを表現する図形又はシンボルをアイコンと称する。また、このとき、各化合物（葉ノード）に相当するアイコンを、第２の特性情報に基づいて色分けする。さらに、画像データ生成部１５は、葉ノードである各化合物が属する親ノードを枠で囲み、枠で囲んだ各ノードが属するさらなる親ノードを枠で囲み…、といった風に、順次、ノードを枠で囲んでいく。これにより、葉ノードである各化合物が、どのような階層構造のどのノードに属しているかが表現される。
【００５１】
なお、上記の色分けに用いられる色は必ずしもカラーである必要はなく、モノトーンの濃淡（グレースケール）等でもよい。また、化合物の機能特性を色によって表現する代わりに、網掛け等の模様や、アイコンの形状によって表現してもよい。もちろん、これらを適宜組み合わせて用いてもよい。
【００５２】
表示部１６は、画像データ生成部１５によって生成された画像データに基づいた画像を表示するためのものである。表示部１６のハードウェア構成としては、画像をカラー表示できるものであれば特に限定されるものではなく、例えば、カラー表示可能な液晶ディスプレイやＣＲＴ（Cathode Ray Tube）ディスプレイ等の各種ディスプレイを備えていてもよいし、カラープリンタを備えていてもよい。これにより、表示部１６は、画像データ生成部１５によって生成された画像データによって表現された画像を表示することができる。
【００５３】
なお、上記の階層分類部１３、化合物配置部１４、及び画像データ生成部１５は、ＣＰＵ等の演算装置がＲＯＭやＲＡＭ等の記憶装置に格納されたプログラムコードを実行することで実現される機能ブロックである。
【００５４】
次に、本実施形態の化合物群表示装置の動作について説明する。図２は、本実施形態の化合物群表示方法の工程図である。
【００５５】
まず、ステップＳ１０において、機能特性入力部１２が機能特性情報を取得する。取得する機能特性情報は、本実施形態では、１次元ベクトルで表される量的変数、カテゴリ変数の何れであってもよい。表１は、取得した機能特性情報の例を示すテーブルである。
【００５６】
【表１】
2006318048.tif 000003
【００５７】
取得した機能特性情報は、画像データ生成部１５に送信される。
【００５８】
次に、ステップＳ１１において、記述子入力部１１が分子記述子を取得する。取得する分子記述子は一種類のみでもよいが、本実施形態ではｎ種類とする。これらの分子記述子は、化合物と対応付けられている。表２は、取得した分子記述子の例を示すテーブルである。
【００５９】
【表２】
2006318048.tif 000004
【００６０】
なお、表２では、各化合物に対応付けられた分子記述子の種類は３つであるが、本発明はこれに限定されるものではなく、任意の個数の分子記述子を取り扱うことができる。
【００６１】
そして、取得した分子記述子は階層分類部１３に送信される。
【００６２】
次に、ステップＳ１２において、階層分類部１３が、受信した分子記述子に基づいて化合物を階層的なクラスターに分類する。なお、本ステップにおいて、特に記述しない場合は、階層分類部１３が各処理を行っているものとする。
【００６３】
説明変数が上述の分子記述子のような量的変数である場合、階層分類部１３は、各化合物をｎ次元空間内の点として取り扱う。すなわち、表１の場合、分子記述子の種類が３つであるので、化合物１は、３次元空間内の点（０．５，２，０）となる。また同様に、化合物２は、点（０．９，４，１）となる。なお、このとき、分子記述子に対して標準化処理を行ってもよい。具体的には、例えば、各分子記述子の化合物間における平均値が０となり、標準偏差値が１となるように標準化すること等が挙げられる。
【００６４】
また、構造情報が複数の分子記述子を要素とする多次元ベクトルからなる場合、構造情報について主成分分析を行い、構造情報の次元（変数の数）を縮約してもよい。
【００６５】
そして、各化合物がこのように配置される空間において、化合物間の距離から化合物間の非類似度を算出する。すなわち、化合物間の距離の遠いものほど化合物間の類似度が小さいものとして扱い、同一クラスターに分類されないようにする。なお、非類似度を算出する際の距離としては、例えば、ユークリッド距離、マンハッタン距離、及びべき乗距離等のうち好適な距離を利用することができる。
【００６６】
なお、変形例としては、距離に基づく非類似度の代わりに類似度を算出してもよい。類似度を算出する方法としては、例えば、ピアソン相関係数、Tanimoto係数、コサイン関数等を利用することができる。
【００６７】
また、距離に基づいて階層的なクラスターに分類するアルゴリズムとしては、例えば、最近隣法、群平均法、又はＷａｒｄ法等の公知のものを用いることができる。
【００６８】
以上の処理により、例えば図３のようなデンドログラムを表現できる各クラスター間の非類似度（距離）及び階層構造の情報が得られる。なお、図３における階層構造及び距離は正確なものではない。得られた情報のうち、階層構造の情報が、化合物配置部１４に送信される。表３は、階層構造を示すデータ（階層構造テーブル）の一例である。
【００６９】
【表３】
2006318048.tif 000005
【００７０】
表３に示すように、階層構造テーブルには、枝ノードのＩＤと、その枝ノードに含まれる子ノード及び葉ノードのＩＤが対応付けられて格納されている。
【００７１】
次に、ステップＳ１３において、化合物配置部１４が、上記階層構造テーブルによって表現される階層構造に基づいて、２次元平面内に各化合物（葉ノード）を配置する。このとき、同一の枝ノードに含まれる葉ノード同士が近接するように、かつ、隣接する葉ノード同士の間隔が葉ノード間の距離にかかわらず略一定の間隔になるように、各葉ノードを配置する。
【００７２】
図４は、２次元平面内における化合物の配置の一例を示す図である。なお、図中、各葉ノードは、黒塗りの正方形で表現されている。また、各葉ノードを囲む枠は、階層構造を視覚的に認識できるように付記したものであり、化合物配置部１４によって作成される配置データには含まれていない。すなわち、化合物配置部１４によって生成される配置データには、化合物（葉ノード）とその化合物の表示位置（例えば、図中の正方形における左上の座標）の情報とが対応付けて含まれている。生成された配置データは、画像データ生成部１５に送信される。
【００７３】
次に、ステップＳ１４において、画像データ生成部１５が、受信した配置データに基づいて、画像データを生成する。なお、生成された画像データでは、各化合物（葉ノード）が枝ノードを示す枠によって囲まれ、枠に囲まれた枝ノードがさらにその親ノードを示す枠によって囲まれ、再帰的な入れ子構造となっている。また、各化合物（葉ノード）のアイコンは、ステップＳ１０で取得した機能特性情報に基づいて、色分けされる。
【００７４】
例えば、入力された機能特性情報がカテゴリ変数である場合、画像データ生成部１５は、各カテゴリ変数と色とが対応付けられたテーブル等に基づいて、化合物の色を決定する。すなわち、画像データ生成部１５は、同じカテゴリ変数を有する化合物が同じ色となるように、各化合物の色を決定する。図５は、化合物１と化合物２、化合物３と化合物４が、それぞれ同じカテゴリ変数からなる機能特性情報を有している場合の色分けの一例である。
【００７５】
また、入力された機能特性情報が量的変数である場合、画像データ生成部１５は、各機能特性情報の値を、変換式によって色情報（例えば、Ｒ，Ｇ，Ｂの何れか１つの色における明度等）に変換して、アイコンの色が機能特性情報の値に応じて連続的に変化するように、化合物の色を決定してもよい。
【００７６】
このようにして生成された画像データは、表示部１６に送信される。なお、ステップＳ１４において、画像データ生成部１５は、平面内に配置された各化合物（葉ノード）を示すアイコンに対して、高さを与えてもよい。例えば、各化合物が１次元ベクトルで表現される第３の特性情報を有している場合、第３の特性情報の値を高さに変換し、各化合物のアイコンが高さを有するように立体的に描画された画像データを生成してもよい。これにより、化合物群表示装置１は、第１から第３の特性情報を関連付けて表示することができる。
【００７７】
最後に、ステップＳ１５において、表示部１６が、受信した画像データに対応する画像を表示する。図５は表示部１６による表示の一例である。なお、化合物群表示装置１では、表示部１６によって表示された画像における各化合物をポインタで指定すると、化合物の分子記述子や機能特性情報が表示されることが好ましい。
【００７８】
〔実施形態２〕
本発明の別の実施形態について図６及び図７に基づいて説明すると以下の通りである。なお、上述した実施形態１と同一の機能を有する部材については、実施形態１と同一の番号を付記し、説明を省略する。
【００７９】
本実施形態の化合物群表示装置２は、複数の化合物について、第１の特性情報と第２の特性情報とを関連付けて表示するものである。
【００８０】
第１の特性情報の種類としては特に限定されるものではなく、量的変数、カテゴリ変数の何れであってもよく、また、それらの変数の数は、単一であっても複数であってもよい。なお、本実施形態では、実施形態１と同様に複数の分子記述子が用いられるものとする。
【００８１】
また、第２の特性情報についても、第１の特性情報と同様に特に限定されるものではなく、例えば、化合物の生物活性、毒性、物理化学的性質、又は薬物動態学的性質等のパラメータを用いることができる。また、これらのパラメータを適宜組み合わせてもよい。本明細書では、これらのパラメータについて以下「機能特性情報」という。
【００８２】
なお、本実施形態の化合物群表示装置２は、化合物を階層的なクラスターに分類する際に、分子記述子を説明変数、機能特性情報を目的変数として、決定木を作成する点が実施形態１と異なっている。これにより、階層的に分類された化合物が、クラスター内において一様な機能特性を示すようになる。よって、構造と機能特性との関係が明確になり、未知の化合物がどのような機能特性を有するかを、その分子記述子に基づいて容易に推定できるようになる。
【００８３】
また、機能特性情報が、１次元ベクトルで表される量的変数、カテゴリ変数に加えて、ｎ次元ベクトル（ｎ≧２）で表される量的変数であってもよい点も実施形態１と異なっている。すなわち、本実施形態の化合物群表示装置２は、機能特性情報がｎ次元で表されるような化合物についても、化合物の構造と機能特性とを関連付けて表示することができる。
【００８４】
図６は、本実施形態の化合物群表示装置２の機能ブロック図である。化合物群表示装置２は、図６に示すように、記述子入力部１１、機能特性入力部１２、階層分類部（階層分類手段）２３、非階層分類部（非階層分類手段）２７、化合物配置部（化合物配置手段）１４、画像データ生成部（画像データ生成手段）１５、及び表示部１６を備えている。
【００８５】
機能特性入力部１２は、実施形態１と基本的に同様であるが、本実施形態では、ｎ個の特性値（量的変数）を取得するものとする。
【００８６】
非階層分類部２７は、機能特性入力部１２に入力されたｎ個の特性値（量的変数）に基づいて化合物を非階層的なクラスターに分類するためのものである。
【００８７】
階層分類部２３は、化合物を分子記述子（第１の特性情報）に基づいて複数の階層的なクラスターに分類するためのものである。より具体的には、分子記述子を説明変数、非階層分類部２７によって形成されたクラスターをカテゴリとして、決定木を作成する。
【００８８】
なお、第２の特性情報が単一の量的変数からなる場合は、階層分類部２３は、分子記述子を説明変数、第２の特性情報を目的変数として決定木（回帰木）を作成することによって、化合物を階層的に分類してもよい。
【００８９】
画像データ生成部２５は、基本的には実施形態１と同様であるが、色分けを行う際に、非階層分類部２７によって形成されたクラスターに基づいた色によって化合物のアイコンが描画された画像データを生成する点が、実施形態１と異なっている。
【００９０】
記述子入力部１１、化合物配置部１４、及び表示部１６については、実施形態１と同様である。
【００９１】
なお、上記の非階層分類部２７、階層分類部２３、化合物配置部１４、及び画像データ生成部２５は、ＣＰＵ等の演算装置がＲＯＭやＲＡＭ等の記憶装置に格納されたプログラムコードを実行することで実現される機能ブロックである。
【００９２】
次に、本実施形態の化合物群表示装置の動作について説明する。図７は、本実施形態の化合物群表示方法の工程図である。
【００９３】
まず、ステップＳ２０において、機能特性入力部１２が機能特性情報を取得する。取得する機能特性情報は、本実施形態では、一例としてｎ次元ベクトルで表される量的変数であるものとする。表４は、取得した機能特性情報の一例である。
【００９４】
【表４】
2006318048.tif 000006
【００９５】
なお、表４では、各化合物に対応付けられた機能特性情報（活性値）の種類は３つであるが、本発明はこれに限定されるものではなく、任意の個数の機能特性情報を取り扱うことができる。
【００９６】
そして、取得した機能特性情報は非階層分類部２７に送信される。
【００９７】
次に、ステップＳ２１において、記述子入力部１１が分子記述子を取得する。このステップについては、実施形態１のステップＳ１１と同様である。取得した分子記述子は階層分類部２３に送信される。
【００９８】
次に、ステップＳ２２において、非階層分類部２７が、受信した機能特性情報に基づいて化合物を非階層的なクラスターに分類する。以下、機能特性情報に基づいて分類したクラスターのことを機能クラスターともいう。ここで、非階層的なクラスターに分類する方法としては、公知のものを利用することができる。一般的なものとしては、ｋ−ｍｅａｎｓ法等が挙げられる。また、ニューラルネットワークモデル、例えば、コホネン（Kohonen）の自己組織化マップ法を用いてもよい。さらに、複数の方法を組み合わせてもよい。
【００９９】
これにより、各化合物の所属する機能クラスターと、各機能クラスターの平均的な機能特性とが得られる。換言すれば、各化合物が、どのような機能特性を有するどのクラスターに所属するのかが決定される。なお、各クラスターの機能特性は、例えばｋ−ｍｅａｎｓ法の場合は、クラスターの重心位置の座標として得られ、自己組織化マップ法の場合は、ノードのパターンベクトルとして得られる。
【０１００】
表５は、分類によって得られた、化合物の所属する機能クラスターが示されたテーブルの一例である。
【０１０１】
【表５】
2006318048.tif 000007
【０１０２】
表５に示すように、上記テーブルでは、化合物と機能クラスターのＩＤとが対応付けられている。
【０１０３】
分類によって得られた、化合物の所属する機能クラスターの情報は、階層分類部２３及び画像データ生成部２５に送信される。
【０１０４】
次に、ステップＳ２３において、階層分類部２３が、記述子入力部１１から受信した分子記述子に基づいて、化合物を階層的なクラスターに分類する。ここで、本実施形態では、クラスターに分類する際に、各クラスター内の化合物の有する機能特性がなるべく一様になるように分類する。具体的には、化合物の分子記述子を説明変数、非階層分類部２７から受信した化合物の所属する機能クラスターを目的変数として、決定木を作成する。このときの分岐ルールとしては、公知のものを用いることができ、例えば、クインラン（Quinlan）によるＣ５．０、Ｃ＆ＲＴ、ＱＵＥＳＴ、ＣＨＡＩＤ又は実施例に示すミッチェルらによる情報利得Ｇを用いる手法等が挙げられる。
【０１０５】
これにより、実施形態１の表３のような階層構造の情報が得られる。得られた階層構造の情報は、化合物配置部１４に送信される。
【０１０６】
次に、ステップＳ２４において、化合物配置部１４が、階層分類部２３によって形成された階層構造の情報に基づいて、２次元平面内に各化合物（葉ノード）を配置する。このステップは、実施形態１のステップＳ１３と同様である。
【０１０７】
次に、ステップＳ２５において、画像データ生成部２５が、受信した配置データに基づいて、画像データを生成する。このステップは、各化合物（葉ノード）がステップＳ２３で決定された機能クラスターに基づいて色分けされる点を除いて、実施形態１のステップＳ１４と同様である。
【０１０８】
最後に、ステップＳ２６において、表示部１６が、受信した画像データに対応する画像を表示する。このステップは、実施形態１のステップＳ１５と同様である。
【０１０９】
以上のように、本実施形態の化合物群表示装置２は、化合物を階層的なクラスターに分類する際に、分子記述子を説明変数、機能特性情報を目的変数として、決定木を作成する。これにより、階層的に分類された化合物が、クラスター内において一様な機能特性を示すようになる。よって、構造と機能特性との関係が明確になり、未知の化合物がどのような機能特性を有するかを、その分子記述子に基づいて容易に推定できるようになる。
【０１１０】
また、本実施形態の化合物群表示装置２は、各化合物に関する機能特性情報が複数個の特性値からなるような場合であっても、化合物を機能特性情報に基づいて非階層的なクラスター（機能クラスター）に分類し、所属する機能クラスターをその化合物の機能特性情報として扱う。そして、画像データ生成部２５がアイコンの色分けをする際に、この機能クラスターに応じて色分けすることによって、化合物の機能特性を表現している。
【０１１１】
なお、上記各実施形態では、化合物群表示装置を構成する階層分類部、非階層分類部、化合物配置部、及び画像データ生成部が、「ＣＰＵなどの演算装置がＲＯＭやＲＡＭなどの記憶装置に格納されたプログラムコードを実行することで実現される機能ブロックである」場合を例にして説明したが、同様の処理を行うハードウェアで実現してもよい。また、処理の一部を行うハードウェアと、当該ハードウェアの制御や残余の処理を行うプログラムコードを実行する上記演算手段とを組み合わせても実現することもできる。さらに、上記各部材のうち、ハードウェアとして説明した部材であっても、処理の一部を行うハードウェアと、当該ハードウェアの制御や残余の処理を行うプログラムコードを実行する上記演算手段とを組み合わせても実現することもできる。なお、上記演算手段は、単体であってもよいし、装置内部のバスや種々の通信路を介して接続された複数の演算手段が共同してプログラムコードを実行してもよい。
【０１１２】
上記演算手段によって直接実行可能なプログラムコード自体、または、後述する解凍などの処理によってプログラムコードを生成可能なデータとしてのプログラムは、当該プログラム（プログラムコードまたは上記データ）を記録媒体に格納し、当該記録媒体を配付したり、あるいは、上記プログラムを、有線または無線の通信路を介して伝送するための通信手段で送信したりして配付され、上記演算手段で実行される。
【０１１３】
なお、通信路を介して伝送する場合、通信路を構成する各伝送媒体が、プログラムを示す信号列を伝搬し合うことによって、当該通信路を介して、上記プログラムが伝送される。また、信号列を伝送する際、送信装置が、プログラムを示す信号列により搬送波を変調することによって、上記信号列を搬送波に重畳してもよい。この場合、受信装置が搬送波を復調することによって信号列が復元される。一方、上記信号列を伝送する際、送信装置が、デジタルデータ列としての信号列をパケット分割して伝送してもよい。この場合、受信装置は、受信したパケット群を連結して、上記信号列を復元する。また、送信装置が、信号列を送信する際、時分割／周波数分割／符号分割などの方法で、信号列を他の信号列と多重化して伝送してもよい。この場合、受信装置は、多重化された信号列から、個々の信号列を抽出して復元する。いずれの場合であっても、通信路を介してプログラムを伝送できれば、同様の効果が得られる。
【０１１４】
ここで、プログラムを配付する際の記録媒体は、取外し可能である方が好ましいが、プログラムを配付した後の記録媒体は、取外し可能か否かを問わない。また、上記記録媒体は、プログラムが記憶されていれば、書換え（書き込み）可能か否か、揮発性か否か、記録方法および形状を問わない。記録媒体の一例として、磁気テープやカセットテープなどのテープ、あるいは、フロッピー（登録商標）ディスクやハードディスクなどの磁気ディスク、または、ＣＤ−ＲＯＭや光磁気ディスク（ＭＯ）、ミニディスク（ＭＤ）やデジタルビデオディスク（ＤＶＤ）などのディスクが挙げられる。また、記録媒体は、ＩＣカードや光カードのようなカード、あるいは、マスクＲＯＭやＥＰＲＯＭ、ＥＥＰＲＯＭまたはフラッシュＲＯＭなどのような半導体メモリであってもよい。あるいは、ＣＰＵなどの演算手段内に形成されたメモリであってもよい。
【０１１５】
なお、上記プログラムコードは、上記各処理の全手順を上記演算手段へ指示するコードであってもよいし、所定の手順で呼び出すことで、上記各処理の一部または全部を実行可能な基本プログラム（例えば、オペレーティングシステムやライブラリなど）が既に存在していれば、当該基本プログラムの呼び出しを上記演算手段へ指示するコードやポインタなどで、上記全手順の一部または全部を置き換えてもよい。
【０１１６】
また、上記記録媒体にプログラムを格納する際の形式は、例えば、実メモリに配置した状態のように、演算手段がアクセスして実行可能な格納形式であってもよいし、実メモリに配置する前で、演算手段が常時アクセス可能なローカルな記録媒体（例えば、実メモリやハードディスクなど）にインストールした後の格納形式、あるいは、ネットワークや搬送可能な記録媒体などから上記ローカルな記録媒体にインストールする前の格納形式などであってもよい。また、プログラムは、コンパイル後のオブジェクトコードに限るものではなく、ソースコードや、インタプリトまたはコンパイルの途中で生成される中間コードとして格納されていてもよい。いずれの場合であっても、圧縮された情報の解凍、符号化された情報の復号、インタプリト、コンパイル、リンク、または、実メモリへの配置などの処理、あるいは、各処理の組み合わせによって、上記演算手段が実行可能な形式に変換可能であれば、プログラムを記録媒体に格納する際の形式に拘わらず、同様の効果を得ることができる。
【実施例】
【０１１７】
〔実施例１〕
本発明の一実施例について以下に説明する。なお、本発明はこれに限定されるものではない。本実施例は、化合物の水に対する溶解度の測定結果の表示に関するものである。具体的には、化合物の構造と溶解度とを関係付けて表示した。
【０１１８】
本実施例では、第２の特性情報（機能特性情報）として、文献（Y. Ran、 et al., Prediction of aqueous solubility of organic compounds by the general solubility equation (GSE). J. Chem. Inf. Comput. Sci., 41: 1208-1217、 2001; G. Yan, et al., Prediction of the aqueous solubility: comparison of the general solubility equation and the method using an amended solvation energy relationship. J. Pharm. Sci., 91: 517-533, 2002）より収集した化合物の溶解度データを用いた。収集した溶解度データの総数、すなわち、本実施例で扱った化合物の総数は908個であった。なお、溶解度のデータは、化合物と対応付けられたテーブルとして取得した。
【０１１９】
また、第１の特性情報（構造情報）としては、Lipinskiら（C.A. Lipinski, et al., Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings. Adv. Drug Deliv. Rev. 23: 3-25, 1997）が薬剤様化合物を識別するために提示している４つの分子記述子、具体的にはClogP、分子量、水素結合供与基数、及び水素結合受容基数を用いた。これらの分子記述子のデータは、化合物と対応付けられたテーブルとして取得した。この対応テーブルの一部を表６に示す。
【０１２０】
【表６】
2006318048.tif 000008
【０１２１】
そして、これらの各分子記述子について、化合物間で平均が０、かつ、標準偏差が１となるような標準化を行った。標準化を行った後の各分子記述子の対応テーブルを表７に示す。
【０１２２】
【表７】
2006318048.tif 000009
【０１２３】
次に、標準化後の分子記述子に基づいて、９０８個の化合物を階層的なクラスターに分類した。分類するアルゴリズムには、Ｗａｒｄ法（J.H. Ward, Hierarchical Grouping to Optimize an Objective Function. J. Am. Stat. Assoc., 58: 236-244、 1963）を用いた。Ｗａｒｄ法は、初期状態（未結合状態）における情報の欠損を０とし、クラスター合併における各結合のステップで情報欠損ΔＳpqの増加を最小にするようなクラスターの組を選ぶ方法であり、情報欠損ΔＳpqは次式で与えられる。
【０１２４】
【数１】
2006318048.tif 000010
【０１２５】
ただし、バーのついたｘpm、ｘqmは、それぞれ結合前のクラスターp、クラスターqの重心におけるｍ番目の分子記述子の値であり、ｎp、ｎqはクラスターｐ、ｑに属するクラスター又は化合物の数である。また、Ｌは、本実施例では分子記述子の種類数、すなわち４である。
【０１２６】
この方法に基づき、クラスター数が１になるまで、クラスターの合併を繰り返した。そして、情報欠損ΔＳpqの閾値を１．０に設定し、それ以下のものは同一階層とみなして９０８個の化合物を分類した。
【０１２７】
次に、階層的に分類した化合物について、階層的なクラスターに基づいた再帰的な入れ子構造となるように、２次元平面内における配置を決定した。この配置を決定するアルゴリズムには、平安京ビュー（「平安京ビュー  〜  階層型データを碁盤上に配置する視覚化手法」、伊藤貴之、小山田耕二、可視化情報学会第9回ビジュアリゼーションカンファレンス抄録、2003年発行）を用いた。
【０１２８】
そして、９０８個の化合物について、決定された配置で、溶解度の対数値を色相とした色でアイコンを描画した画像データを作成した。この画像データに基づく画像を図８に示す。なお、図中、化合物を示すアイコンはカラー表示されている。
【０１２９】
図８に示すように、比較的溶解度の高いもののみ、もしくは、低いもののみが集まっているクラスターが認められ、化学構造が類似する化合物は比較的良く似た溶解度を有することが示された。
【０１３０】
また、各化合物を示すアイコンに高さを与えたものを図９（ａ）、図９（ｂ）に示す。図９（ａ）は、各化合物のClogPをバーの高さによって表したものである。一方、図９（ｂ）は、ClogPの代わりに分子量をバーの高さによって表したものである。このように、９０８個もの化合物に関する非常に大規模なデータであるにもかかわらず、ClogPの値が小さいほど、また分子量が小さいほど水に対して溶けやすいことを視覚的に容易に捉えることができた。
【０１３１】
〔実施例２〕
本発明の一実施例について以下に説明する。本実施例は、薬物代謝酵素チトクロムＰ４５０（ＣＹＰ）による化合物（薬物）の代謝パターンの分類に関するものであり、化合物の構造と代謝パターンとを関連付けて表示した。
【０１３２】
なお、ＣＹＰは、生体の物質代謝や異物解毒に関わる重要な薬物代謝酵素群の一つであり、ヒトでは約２０種類の分子種が同定されている。生体にとって異物である薬物の多くはＣＹＰによる代謝を受けて消失するため、ＣＹＰによる代謝は薬物の有効性や安全性と深く関わっているといえる。ＣＹＰスーパーファミリーの中でも、CYP1A2、CYP2C9、CYP2C19、CYP2D6、CYP3A4の５つは薬物代謝に関わる主要な分子種であり、薬物代謝の約９０％がこれらによって説明できると言われている。したがって、これら５つのＣＹＰによる薬物の代謝を評価・予測することは創薬研究において非常に重要な課題となっている。
【０１３３】
本実施例では、第２の特性情報（機能特性情報）として、ボナブリーらが調査した、１６１個の薬物についての、上記の５つのＣＹＰによる代謝感受性のデータ（P. Bonnabry, et al., Quantitative drug interactions prediction system (Q-DIPS): a dynamic computer-based method to assist in the choice of clinically relevant in vivo studies. Clin. Pharmacokinet., 40: 631-640, 2001）を用いた。なお、ＣＹＰによる代謝感受性のデータは、３段階で評価されている。代謝感受性のデータは、化合物と対応付けられたテーブルとして取得した。実施例１と同様に標準化した後の代謝感受性のデータの一部を表８に示す。
【０１３４】
【表８】
2006318048.tif 000011
【０１３５】
また、第１の特性情報（構造情報）として、Molconn-Zによって計算される２２０個の分子記述子（Molconn-Zディスクリプタ）のうち、パス長が７以上のもの、及びデータセット内での分散が０となる分子記述子を除いた計１１５個の分子記述子を用いた。なお、各薬物のMolconn-Zディスクリプタは、Molconn-Z（商標、Hall Associate Consulting社製）を用いて化学構造から算出した。各化合物のMolconn-Zディスクリプタの一部を表９に示す。
【０１３６】
【表９】
2006318048.tif 000012
【０１３７】
次に、本実施例では、各化合物の機能特性情報（代謝パターン）が５つの量的変数（５次元ベクトル）からなるため、化合物を機能特性情報に基づいて非階層的な機能クラスターに分類した。機能クラスターに分類するアルゴリズムには、コホネンの自己組織化マップ法（T. Kohonen, “Self-Organizing Maps”, Springer, Berlin, Heidelberg, 1995）とｋ−ｍｅａｎｓ法とを組み合わせたものを用いた。これらについて以下に詳細に説明する。
【０１３８】
まず、各ＣＹＰによる代謝強度を属性とする５次元のパターンベクトルを考え、取得したＣＹＰによる代謝感受性のデータを用いて自己組織化マップ（７×９＝６３ノード）の学習を行った。学習によって得られた各ノードのパターンベクトルを示すテーブルの一部を表１０に示す。
【０１３９】
【表１０】
2006318048.tif 000013
【０１４０】
また、上記の表１０を可視化したものを図１０に示す。図１０では、自己組織化マップの各ノードが有するパターンベクトルの各成分を色相によって表示している。図１０において、CYP3A4に対応する要素については広範囲のノードが高い値を示すことから、多くの薬物がCYP3A4によって代謝されることが視覚的に捉えられた。また、CYP1A2とCYP2C9あるいはCYP2C19との間ではほとんどオーバーラップがないことから、CYP1A2はCYP2C9やCYP2C19とは大きく異なった基質選択性を示すことが明らかとなった。
【０１４１】
次に、自己組織化マップにおけるノード数は６３個であり、機能クラスターとしてそのまま用いるには数が多いため、ノードをｋ−ｍｅａｎｓ法によって最適なクラスター数にさらに分類することにした。ここで、最適なクラスター数を求める指標として、次式で表されるDavies-Bouldin Index（ＤＢ）を用いた。
【０１４２】
【数２】
2006318048.tif 000014
【０１４３】
ただし、Ｓn（Ｑi）はクラスターＱiにおけるクラスター中心と各クラスター要素とのユークリッド距離の平均値を表し、Ｓ（Ｑi,Ｑj）はクラスターＱiとＱjの中心間のユークリッド距離を表す。
【０１４４】
ＤＢにより、最適なクラスター数は６であることが求められた。
【０１４５】
次に、自己組織化マップにおける６３個のノードをｋ−ｍｅａｎｓ法によって６つのクラスターに分類し、これを機能クラスターとした。これにより、自己組織化マップにおける各ノードが６つのうちのどのクラスターに属するかを示すテーブルが得られた。そのテーブルの一部を表１１に示す。
【０１４６】
【表１１】
2006318048.tif 000015
【０１４７】
また、表１１を可視化したものを図１１に示す。図１１において、各ノードが所属するクラスターは、色相によって表現されている。
【０１４８】
次に、各薬物の代謝のパターンベクトルと自己組織化マップの各ノードのパターンベクトルとの距離から薬物の帰属するノード（ベストマッチノード）を決定し、そのノードが上記の６つのクラスターの何れに帰属するかに応じて、薬物を機能クラスターに分類した。表１２は各機能クラスターに分類された薬物のＣＹＰによる代謝感受性の平均スコアを示したものである。
【０１４９】
【表１２】
2006318048.tif 000016
【０１５０】
表１２から、機能クラスター１はCYP3A4及びCYP2C9による代謝が、機能クラスター２はCYP2C19による代謝が有意に高く、機能クラスター３はCYP2C9でのみ、機能クラスター４はCYP3A4でのみ、機能クラスター５はCYP2D6によって、機能クラスター６はCYP1A2によって代謝される薬物の集合であることがわかる。
【０１５１】
次に、構造情報に基づいて化合物を階層的に分類するために、二分木を作成した。二分木を作成するにあたって、説明変数としてMolconn-Zディスクリプタ、目的変数として機能クラスターを用い、分岐ルールには次に示す情報利得Ｇ（T.M. Mitchell, “Machine Learning”, McGraw-Hill, Singapore, 1997）を用いた。
【０１５２】
【数３】
2006318048.tif 000017
【０１５３】
ただし、Ｎは薬物を機能クラスターに分類した際のクラスター数（＝６）を示し、ｓは分割前の集合における薬物の総数を、ｃiは機能クラスターｉに属する薬物数を示し、ｓ1及びｓ2は二分割されたそれぞれの集合における薬物の総数を、ｃ1,i及びｃ2,iは二分割されたそれぞれの集合における、機能クラスターｉに属する薬物数を示す。
【０１５４】
ここで、二分木の各分割ステップにおいて、最大の情報利得Ｇが得られる条件（分割の判定に用いるMolconn-Zディスクリプタの種類及び閾値）を決定し、得られる最大の情報利得Ｇの値が０．２５以下になるまで再帰的に分割を繰り返した。二分木の作成によって得られたデータの一部を表１３に示す。
【０１５５】
【表１３】
2006318048.tif 000018
【０１５６】
表１３に示すように、ルートノード（ノード１）では、Molconn-ZディスクリプタのSHCsatsが4.80369未満か否かで薬物を２つの子ノード（ノード２及びノード３）に分割している。そして、同様に、ノード２では、Molconn-ZディスクリプタのSHHBAが31.98136未満か否かで薬物を２つの子ノード（ノード３０及びノード３１）に分割している。
【０１５７】
次に、階層的に分類した薬物について、階層的なクラスターに基づいた再帰的な入れ子構造となるように、２次元平面内における配置を決定した。この配置を決定するアルゴリズムには、平安京ビューを用いた。
【０１５８】
そして、１６１個の薬物について、決定された配置で、機能クラスターの種類を色相とした色でアイコンを描画した画像データを作成した。この画像データに基づく画像を図１２に示す。なお、画像中において、薬物を示すアイコンはカラー表示されており、アイコンの色は、代謝パターンに基づいて分類した際に帰属した機能クラスターに応じた色となっている。
【０１５９】
図１２に示すように、表示された画像では、共通の構造的特徴を有している薬物が同じ代謝パターンを示すように分類されている。このように、１６１個もの大規模な数の薬物について、構造と代謝パターンとの関係を効果的に視覚化できることが示された。
【０１６０】
また、図１３は、チトクロムＰ４５０の各分子種CYP1A2、CYP2C9、CYP2C19、CYP2D6、CYP3A4による代謝感受性を個別の画像で示したものである。画像中において、代謝感受性（３段階）はアイコンの色で表現されている。このように、表示された画像を一望するだけで各薬物代謝酵素の感受性を包括的に理解することができ、本発明の有用性が示された。
【０１６１】
なお、本発明は上述した実施形態及び実施例のみに限定されるものではなく、請求項に示した範囲で種々の変更が可能であり、異なる実施形態及び実施例にそれぞれ開示された技術的手段を適宜組み合わせて得られる実施形態についても本発明の技術的範囲に含まれる。
【産業上の利用可能性】
【０１６２】
本発明は、スクリーニングデータの全体像を効果的に視覚化できるので、候補化合物の選択を支援するとともに、機能的にさらに優れた化合物の合成展開を行う方向性を示すことができる。したがって、医薬品、農薬に限らず、その他の機能性化学物質の探索研究を行う化学工業分野一般において利用可能であり、非常に有用性が高いと考えられる。
【図面の簡単な説明】
【０１６３】
【図１】本発明の一実施形態を示すものであり、化合物群表示装置の機能ブロック図である。
【図２】本発明の一実施形態を示すものであり、化合物群表示装置の処理工程を示す工程図である。
【図３】階層的な分類によって得られたデンドログラムの一例を示す図である。
【図４】化合物配置部による化合物の配置の一例を示す図である。
【図５】画像データ生成部によって生成された画像データの一例を示す図である。
【図６】本発明の別の実施形態を示すものであり、化合物群表示装置の機能ブロック図である。
【図７】本発明の別の実施形態を示すものであり、化合物群表示装置の処理工程を示す工程図である。
【図８】本発明の一実施例を示すものであり、化合物群表示装置によって表示された画像を示す図である。
【図９】本発明の一実施例を示すものであり、化合物群表示装置によって表示された画像を示す図である。
【図１０】本発明の別の実施例を示すものであり、自己組織化マップにおける各ノードのパターンベクトルの各成分を色相で示した図である。
【図１１】本発明の別の実施例を示すものであり、自己組織化マップにおける各ノードが帰属する機能クラスターを色相で示した図である。
【図１２】本発明の別の実施例を示すものであり、化合物群表示装置によって表示された画像を示す図である。
【図１３】本発明の別の実施例を示すものであり、化合物群表示装置によって表示された画像を示す図である。
【符号の説明】
【０１６４】
１，２        化合物群表示装置
１３，２３    階層分類部（階層分類手段）
１４          化合物配置部（化合物配置手段）
１５，２５    画像データ生成部（画像データ生成手段）
１６          表示部
【図１】
2006318048.tif 000019
【図２】
2006318048.tif 000020
【図３】
2006318048.tif 000021
【図４】
2006318048.tif 000022
【図５】
2006318048.tif 000023
【図６】
2006318048.tif 000024
【図７】
2006318048.tif 000025
【図８】
2006318048.tif 000026
【図９】
2006318048.tif 000027
【図１０】
2006318048.tif 000028
【図１１】
2006318048.tif 000029
【図１２】
2006318048.tif 000030
【図１３】
2006318048.tif 000031