(57)【要約】  【課題】  経験的力場関数を用いたエネルギー最適化を基礎とすることなく、活性未知な化合物や、機能未知または立体構造が未知であるタンパク質に対しても適用可能な、タンパク質と化合物との相互作用に関するデータ処理装置を提供する。  【解決手段】  ２以上の既知タンパク質相互の距離を示すプロファイルもしくはカーネル、および、２以上の既知化合物相互の距離を示すプロファイルもしくはカーネルを、それぞれクラスタリングするクラスタリング処理部と、前記クラスタリング処理部による既知タンパク質のクラスタリング結果および既知化合物のクラスタリング結果に従って、前記既知タンパク質と既知化合物との既知相互作用の情報に基づき、前記既知タンパク質と既知化合物に対する前記既知相互作用の分布データを作成する分布データ作成部とを備えるデータ処理装置とする。【選択図】  図１
2007011752.tif 000002
【特許請求の範囲】
【請求項１】
２以上の既知タンパク質相互の距離を示すプロファイルもしくはカーネルおよび２以上の既知化合物相互の距離を示すプロファイルもしくはカーネルを、それぞれクラスタリングするクラスタリング処理部と、
前記クラスタリング処理部による既知タンパク質のクラスタリング結果および既知化合物のクラスタリング結果に従って、前記既知タンパク質と既知化合物との既知相互作用の情報に基づき、前記既知タンパク質と既知化合物に対する前記既知相互作用の分布データを作成する分布データ作成部とを備えるデータ処理装置。
【請求項２】
さらに、前記既知タンパク質のクラスタリング結果および既知化合物のクラスタリング結果に従って、前記既知タンパク質と既知化合物とを二軸座標に表し、前記二軸座標に、前記既知相互作用の情報をマッピングし、前記分布データを表示するマップ表示部を備える、請求項１記載のデータ処理装置。
【請求項３】
さらに、２以上の既知タンパク質の配列情報、２以上の既知化合物の情報および前記既知タンパク質と既知化合物との既知相互作用の情報を入力する第１の情報入力部と、
前記第１の情報入力部より入力された前記既知タンパク質の配列情報および前記既知化合物の情報を、それぞれ前記タンパク質相互の距離を示すプロファイルまたはカーネル、および、前記既知化合物相互の距離を示すプロファイルまたはカーネルに変換し、前記タンパク質相互の距離を示すプロファイルまたはカーネル、および、前記既知化合物相互の距離を示すプロファイルまたはカーネルを、前記クラスタリング処理部に出力するプロファイル変換部とを備える、請求項１または２記載のデータ処理装置。
【請求項４】
さらに、前記既知タンパク質相互の距離を示すプロファイルもしくはカーネル、前記既知化合物相互の距離を示すプロファイルもしくはカーネル、および前記既知タンパク質と既知化合物との既知相互作用の情報を入力する第２の情報入力部を備える、請求項１または２記載のデータ処理装置。
【請求項５】
さらに、前記既知タンパク質のプロファイルもしくはカーネル、前記既知化合物のプロファイルもしくはカーネル、および前記分布データ作成部による前記既知相互作用の分布データを教師データとして、タンパク質と化合物との相互作用に関する統計的法則を機械学習し、前記機械学習による統計的法則から得られる結果より、目的タンパク質に対して相互作用を示す化合物または目的化合物が相互作用を示すタンパク質を予測する予測部を備える、請求項１〜４のいずれか一項に記載のデータ処理装置。
【請求項６】
前記予測部が、機械学習により統計的法則から得られる結果をパラメータとして算出し、前記パラメータに基づいて予測を行う、請求項５記載のデータ処理装置。
【請求項７】
さらに、前記分布データ作成部による前記既知相互作用の分布データより、タンパク質と化合物との相互作用が適切な範囲を選択し、これを教師データとして前記予測部に出力するデータ選択部を備える、請求項５または６記載のデータ処理装置。
【請求項８】
さらに、予測部による予測結果を出力する予測結果出力部を備える、請求項５〜７のいずれか一項に記載のデータ処理装置。
【請求項９】
さらに、２以上の既知タンパク質の配列情報、２以上の既知化合物の情報および前記既知タンパク質と既知化合物との既知相互作用の情報を入力する第１の情報入力部と、
前記第１の情報入力部より入力された前記既知タンパク質の配列情報および前記既知化合物の情報を、それぞれ前記タンパク質相互の距離を示すプロファイルまたはカーネル、および、前記既知化合物相互の距離を示すプロファイルまたはカーネルに変換し、前記タンパク質相互の距離を示すプロファイルまたはカーネル、および、前記既知化合物相互の距離を示すプロファイルまたはカーネルを、前記クラスタリング処理部に出力するプロファイル変換部と、
目的タンパク質の配列情報または目的化合物の情報を入力する第３の情報入力部を備え、
前記プロファイル変換部が、さらに、前記第３の情報入力部より入力された目的タンパク質の配列情報を、前記第１の情報入力部より入力された既知タンパク質の配列情報に基づき、前記目的タンパク質の前記既知タンパク質に対する距離を示すプロファイルもしくはカーネルに変換し、または、前記第３の情報入力部により入力された目的化合物の情報を、前記第１の情報入力部により入力された既知化合物の情報に基づき、前記目的化合物の前記既知化合物に対する距離を示すプロファイルもしくはカーネルに変換し、変換された前記目的タンパク質の前記既知タンパク質に対する距離を示すプロファイルもしくはカーネル、または、前記目的化合物の前記既知化合物に対する距離を示すプロファイルもしくはカーネルを前記予測部に出力し、
前記予測部が、前記統計的法則から得られる結果と、前記プロファイル変換部による前記目的タンパク質の前記既知タンパク質に対する距離を示すプロファイルもしくはカーネル、または、前記目的化合物の前記既知化合物に対する距離を示すプロファイルもしくはカーネルとから、目的タンパク質に対して相互作用を示す化合物または前記目的化合物が相互作用を示すタンパク質を予測する、請求項５〜８のいずれか一項に記載のデータ処理装置。
【請求項１０】
さらに、前記既知タンパク質相互の距離を示すプロファイルもしくはカーネル、前記既知化合物相互の距離を示すプロファイルもしくはカーネル、および前記既知タンパク質と既知化合物との既知相互作用の情報を入力する第２の情報入力部と、
目的タンパク質の前記既知タンパク質に対する距離を示すプロファイルもしくはカーネル、または、目的化合物の前記既知化合物に対する距離を示すプロファイルもしくはカーネルを入力する第４の情報入力部とを備え、
前記予測部において、前記統計的法則から得られる結果と、前記第４の情報入力部より入力された前記目的タンパク質の前記既知タンパク質に対する距離を示すプロファイルもしくはカーネル、または、目的化合物の前記既知化合物に対する距離を示すプロファイルもしくはカーネルとから、目的タンパク質に対して相互作用を示す化合物または前記目的化合物が相互作用を示すタンパク質を予測する、請求項５〜８のいずれか一項に記載のデータ処理装置。
【請求項１１】
前記化合物相互の距離を示すプロファイルもしくはカーネルが、化合物の化学構造および物理化学特性の少なくとも一方の情報から変換されたプロファイルもしくはカーネルである、請求項１〜１０のいずれか一項に記載のデータ処理装置。
【請求項１２】
前記タンパク質相互の距離を示すプロファイルもしくはカーネルが、タンパク質の塩基配列、アミノ酸配列、塩基配列ならびにアミノ酸配列の少なくとも一方から決定されたモチーフ、疎水性指標および二次構造からなる群から選択された少なくとも一つの情報から変換されたプロファイルもしくはカーネルである、請求項１〜１１のいずれか一項に記載のデータ処理装置。
【請求項１３】
前記タンパク質が、立体構造不明のタンパク質である請求項１〜１２のいずれか一項に記載のデータ処理装置。
【請求項１４】
前記タンパク質が、Ｇタンパク質共役型受容体である請求項１〜１３のいずれか一項に記載のデータ処理装置。
【請求項１５】
２以上の既知タンパク質相互の距離を示すプロファイルもしくはカーネル、および、２以上の既知化合物相互の距離を示すプロファイルもしくはカーネルを、それぞれクラスタリングするクラスタリング処理ステップと、
前記既知タンパク質のクラスタリング結果および前記既知化合物のクラスタリング結果に従って、前記既知タンパク質と既知化合物との既知相互作用の情報に基づき、前記既知タンパク質と既知化合物に対する前記既知相互作用の分布データを作成する分布データ作成ステップとを、コンピュータに実行させるデータ処理プログラム。
【請求項１６】
さらに、前記既知タンパク質のクラスタリング結果および既知化合物のクラスタリング結果に従って、前記既知タンパク質と既知化合物とを二軸座標に表し、前記二軸座標に、前記既知相互作用の情報をマッピングし、前記分布データを表示するマップ表示ステップをコンピュータに実行させる、請求項１５記載のデータ処理プログラム。
【請求項１７】
さらに、２以上の既知タンパク質の配列情報、２以上の既知化合物の情報および前記既知タンパク質と既知化合物との既知相互作用の情報の入力を受け付ける第１の入力ステップと、
入力された前記既知タンパク質の配列情報および前記既知化合物の情報を、それぞれ前記タンパク質相互の距離を示すプロファイルもしくはカーネル、および、前記既知化合物相互の距離を示すプロファイルもしくはカーネルに変換し、前記タンパク質相互の距離を示すプロファイルもしくはカーネル、および、前記既知化合物相互の距離を示すプロファイルもしくはカーネルを、前記クラスタリング処理部に出力する第１のプロファイル変換ステップとを、コンピュータに実行させる、請求項１５または１６記載のデータ処理プログラム。
【請求項１８】
さらに、前記既知タンパク質相互の距離を示すプロファイルもしくはカーネル、前記既知化合物相互の距離を示すプロファイルもしくはカーネル、および前記既知タンパク質と既知化合物との既知相互作用の情報の入力を受け付ける第２の入力ステップを、コンピュータに実行させる請求項１５または１６記載のデータ処理プログラム。
【請求項１９】
さらに、前記既知タンパク質のプロファイルもしくはカーネル、前記既知化合物のプロファイルもしくはカーネル、および前記分布データ作成ステップにおける前記既知相互作用の分布データを教師データとして、タンパク質と化合物との相互作用に関する統計的法則を機械学習し、前記機械学習による統計的法則から得られる結果より、目的タンパク質に対して相互作用を示す化合物または目的化合物が相互作用を示すタンパク質を予測する予測ステップを、コンピュータに実行させる請求項１５〜１８のいずれか一項に記載のデータ処理プログラム。
【請求項２０】
前記予測ステップが、前記機械学習により統計的法則から得られる結果をパラメータとして算出し、前記パラメータに基づいて予測を行うステップである、請求項１９記載のデータ処理プログラム。
【請求項２１】
さらに、前記分布データ作成ステップにおける前記既知相互作用の分布データより、タンパク質と化合物との相互作用が適切な範囲を選択し、これを教師データとして前記予測ステップに出力するデータ選択ステップを、コンピュータに実行させる請求項１９または２０記載のデータ処理プログラム。
【請求項２２】
さらに、前記予測ステップにおける予測結果を出力する予測結果出力ステップを、コンピュータに実行させる請求項１９〜２１のいずれか一項に記載のデータ処理プログラム。
【請求項２３】
さらに、２以上の既知タンパク質の配列情報、２以上の既知化合物の情報および前記既知タンパク質と既知化合物との既知相互作用の情報の入力を受け付ける第１の情報入力ステップと、
入力された前記既知タンパク質の配列情報および前記既知化合物の情報を、それぞれ前記タンパク質相互の距離を示すプロファイルもしくはカーネル、および、前記既知化合物相互の距離を示すプロファイルもしくはカーネルに変換し、前記タンパク質相互の距離を示すプロファイルもしくはカーネル、および、前記既知化合物相互の距離を示すプロファイルもしくはカーネルを、前記クラスタリング処理ステップに出力する第１のプロファイル変換ステップと、
目的タンパク質の配列情報または目的化合物の情報の入力を受け付ける第３の情報入力ステップと、
入力された目的タンパク質の配列情報を、入力された既知タンパク質の配列情報に基づき、前記目的タンパク質の前記既知タンパク質に対する距離を示すプロファイルもしくはカーネルに変換し、または、入力された目的化合物の情報を、入力された既知化合物の情報に基づき、前記目的化合物の前記既知化合物に対する距離を示すプロファイルもしくはカーネルに変換し、変換された前記目的タンパク質の前記既知タンパク質に対する距離を示すプロファイルもしくはカーネル、または、前記目的化合物の前記既知化合物に対する距離を示すプロファイルもしくはカーネルを前記予測ステップに出力する第２のプロファイル変換ステップと、
前記統計的法則から得られる結果と、前記目的タンパク質の前記既知タンパク質に対する距離を示すプロファイルもしくはカーネル、または、前記目的化合物の前記既知化合物に対する距離を示すプロファイルもしくはカーネルとから、目的タンパク質に対して相互作用を示す化合物または前記目的化合物が相互作用を示すタンパク質を予測する前記予測ステップとを、コンピュータに実行させる請求項１９〜２１のいずれか一項に記載のデータ処理プログラム。
【請求項２４】
さらに、前記既知タンパク質相互の距離を示すプロファイルもしくはカーネル、前記既知化合物相互の距離を示すプロファイルもしくはカーネル、および前記既知タンパク質と既知化合物との既知相互作用の情報の入力を受け付ける第２の情報入力ステップと、
目的タンパク質の前記既知タンパク質に対する距離を示すプロファイルもしくはカーネル、または、目的化合物の前記既知化合物に対する距離を示すプロファイルもしくはカーネルを入力する第４の情報入力ステップと、
前記統計的法則から得られる結果と、入力された前記目的タンパク質の前記既知タンパク質に対する距離を示すプロファイルもしくはカーネル、または目的化合物の前記既知化合物に対する距離を示すプロファイルもしくはカーネルとから、目的タンパク質に対して相互作用を示す化合物または前記目的化合物が相互作用を示すタンパク質を予測する前記予測ステップとを、コンピュータに実行させる、請求項１９〜２２のいずれか一項に記載のデータ処理プログラム。
【請求項２５】
請求項１５〜２４のいずれか一項に記載のデータ処理プログラムを格納したコンピュータ読み取り可能な記録媒体。
【請求項２６】
タンパク質と化合物との相互作用に関するデータ処理方法であって、
２以上の既知タンパク質相互の距離を示すプロファイルもしくはカーネル、および、２以上の既知化合物相互の距離を示すプロファイルもしくはカーネルを、それぞれクラスタリングするクラスタリング処理工程と、
前記既知タンパク質のクラスタリング結果および前記既知化合物のクラスタリング結果に従って、前記既知タンパク質と既知化合物との既知相互作用の情報に基づき、前記既知タンパク質と既知化合物に対する前記既知相互作用の分布データを作成する分布データ作成工程とを含むデータ処理方法。
【請求項２７】
さらに、前記既知タンパク質のクラスタリング結果および既知化合物のクラスタリング結果に従って、前記既知タンパク質と既知化合物とを二軸座標に表し、前記二軸座標に、前記既知相互作用の情報をマッピングし、前記分布データを表示するマップ表示工程を含む、請求項２６記載のデータ処理方法。
【請求項２８】
さらに、２以上の既知タンパク質の配列情報および２以上の既知化合物の情報から、前記タンパク質相互の距離を示すプロファイルもしくはカーネル、および、前記既知化合物相互の距離を示すプロファイルもしくはカーネルに変換する工程を含む請求項２６または２７記載のデータ処理方法。
【請求項２９】
さらに、前記既知タンパク質のプロファイルもしくはカーネル、前記既知化合物のプロファイルもしくはカーネル、および前記分布データ作成工程における前記既知相互作用の分布データを教師データとして、タンパク質と化合物との相互作用に関する統計的法則を機械学習し、前記機械学習による統計的法則から得られる結果より、目的タンパク質に対して相互作用を示す化合物または目的化合物が相互作用を示すタンパク質を予測する予測工程を含む、請求項２６〜２８のいずれか一項に記載のデータ処理方法。
【請求項３０】
前記予測工程において、機械学習により統計的法則から得られる結果をパラメータとして算出し、前記パラメータに基づいて予測を行う、請求項２９記載のデータ処理方法。
【請求項３１】
さらに、前記分布データ作成工程における前記既知相互作用の分布データより、タンパク質と化合物との相互作用が適切な範囲を選択し、これを教師データとするデータ選択工程を含む、請求項２９または３０記載のデータ処理方法。
【請求項３２】
さらに、予測工程における予測結果を出力する予測結果出力工程を含む、請求項２９〜３１のいずれか一項に記載のデータ処理方法。
【請求項３３】
前記予測工程が、前記統計的法則から得られる結果と前記目的タンパク質の前記既知タンパク質に対する距離を示すプロファイルもしくはカーネル、または、前記目的化合物の前記既知化合物に対する距離を示すプロファイルもしくはカーネルとから、前記目的タンパク質に対して相互作用を示す化合物または前記目的化合物が相互作用を示すタンパク質を予測する工程である、請求項２９〜３２のいずれか一項に記載のデータ処理方法。
【請求項３４】
既知タンパク質と既知化合物との既知相互作用の分布データのデータ構造であって、
２以上の既知タンパク質相互の距離を示すプロファイルもしくはカーネルをクラスタリングして得られた既知タンパク質プロファイルもしくはカーネルと、
２以上の既知化合物相互の距離を示すプロファイルもしくはカーネルをクラスタリングして得られた既知化合物プロファイルもしくはカーネルと、
前記２以上の既知タンパク質のうち１の既知タンパク質に対応する既知タンパク質プロファイルもしくはカーネルと、前記２以上の既知化合物のうち１の既知化合物に対応する既知化合物プロファイルもしくはカーネルとに対応づけて、当該１の既知タンパク質と当該１の既知化合物との相互作用の情報が含まれていることを特徴とするデータ構造。
【請求項３５】
請求項３４に記載のデータ構造を有するデータを記録した記録媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、タンパク質と化合物との相互作用に関するデータ処理装置、データ処理プログラム、それを格納したコンピュータ読み取り可能な記録媒体、およびデータ処理方法に関する。
【背景技術】
【０００２】
医薬品や生理性物質の作用は、標的となるタンパク質や遺伝子との結合によって活性発現されることから、タンパク質または遺伝子（以下、「タンパク質等」という）と化合物との相互作用の予測や、メカニズムの解明は、医薬品開発やバイオ産業に多大な貢献を果たしている。
【０００３】
タンパク質等と化合物との相互作用を予測する従来技術としては、一般に、立体構造モデルを用いた予測システムが知られている。このシステムは、生体高分子の立体構造の情報に基づき、リガンドとの安定な複合体構造およびその結合の強さを推定する方法であり、このようなアプローチはドッキングスタディと呼ばれている。そして、グラフィクス画面において、マニュアル操作によって安定しそうなリガンド分子との複合体構造を探す対話的手法に対して、近年、結合様式や配座についての自由度を自動的に考慮しながら存在可能性の高い複合体構造を探索する自動ドッキング法の開発が進んでいる（非特許文献１、非特許文献２、特許文献１参照）。さらに、自動ドッキング法で得られる各化合物の安定複合体構造がエネルギー値やその他の条件を満たすか否かによって、膨大な化合物の３次元構造データベースから少数の候補化合物を選択し、実験的に活性を確認することによって新規リガンドを見出す方法の開発も進んでいる（特許文献２参照）。
【０００４】
しかしながら、これらのアプローチは、いずれも経験に基づいて定められる力場関数を用いたエネルギー最適化を基礎としている。このため、信頼性の高い相互作用予測は、その対象が、使用する経験的力場パラメータの適用範囲内である化合物群やタンパク質群に限定されてしまう。したがって、例えば、活性未知の化合物や機能未知のタンパク質に対する予測には適用が困難な状態である。また、前述のドッキング手順からわかるように、立体構造が未知のタンパク質や化合物に対しても適用が困難という問題がある。
【０００５】
また、このような研究を行うためには、タンパク質と化合物に関するデータベースの入手が不可欠であり、代表的なデータベースとして、例えば、ＫＥＧＧ（京都大学）、Ｂｉｏｃｈｅｍｉｃａｌ Ｐａｔｈｗａｙｓ(Boehringer Mannheim)、ＷＩＴ(Russian Academy of Sciences)、Ｂｉｏｆｒｏｎｔｉｒｅ(呉羽化学)、Ｐｒｏｔｅｉｎ Ｐａｔｈｗａｙ(AxCell)、ｂｉｏＳＣＯＵＴ(LION)、ＥｃｏＣｙｃ(DoubleTWist)、ＵＭ−ＢＢＤ(ミネソタ大学)、ＫｅｙＭｏｌｎｅｔ(医薬分子設計研)、ＧＰＣＲＤＢ(CMBI)、ＩＵＰＨＡＲなどがあげられる。しかしながら、これらのデータベースは、今日までに実験的に決定された代謝経路やシグナル伝達経路に関する情報をテキスト形式やＧＩＦ画像形式で表現しているため、情報提供にすぎず、相互作用予測に使用するには限界があった。
【非特許文献１】”DOCK” R. L. DesJarlais et al, J. Med. Chem. 29, p.2149, 1986
【非特許文献２】“FlexX” Protesin, 37, p.228, 1999
【特許文献１】米国特許5,642,292 “ADAM”
【特許文献２】国際公開第96/13785号パンフレット
【発明の開示】
【発明が解決しようとする課題】
【０００６】
したがって、本発明は、例えば、従来のように経験的力場関数を用いたエネルギー最適化を基礎とすることなく、活性未知な化合物や、機能未知または立体構造が未知であるタンパク質に対しても適用可能な、タンパク質と化合物との相互作用に関するデータ処理装置、データ処理プログラム、データ処理方法等の提供を目的とする。
【課題を解決するための手段】
【０００７】
前記目的を達成するため、本発明のデータ処理装置は、２以上の既知タンパク質相互の距離を示すプロファイルもしくはカーネル、および、２以上の既知化合物相互の距離を示すプロファイルもしくはカーネルを、それぞれクラスタリングするクラスタリング処理部と、前記クラスタリング処理部による既知タンパク質のクラスタリング結果および既知化合物のクラスタリング結果に従って、前記既知タンパク質と既知化合物との既知相互作用の情報に基づき、前記既知タンパク質と既知化合物に対する前記既知相互作用の分布データを作成する分布データ作成部とを備えることを特徴とする。
【０００８】
また、本発明のデータ処理プログラムは、２以上の既知タンパク質相互の距離を示すプロファイルもしくはカーネル、および、２以上の既知化合物相互の距離を示すプロファイルもしくはカーネルを、それぞれクラスタリングするクラスタリング処理ステップと、前記既知タンパク質のクラスタリング結果および前記既知化合物のクラスタリング結果に従って、前記既知タンパク質と既知化合物との既知相互作用の情報に基づき、前記既知タンパク質と既知化合物に対する前記既知相互作用の分布データを作成する分布データ作成ステップとを、コンピュータに実行させることを特徴とする。また、本発明の記録媒体は、本発明のデータ処理プログラムを格納したコンピュータ読み取り可能な記録媒体である。
【０００９】
また、本発明のデータ処理方法は、タンパク質と化合物との相互作用に関するデータ処理方法であって、２以上の既知タンパク質相互の距離を示すプロファイルもしくはカーネル、および、２以上の既知化合物相互の距離を示すプロファイルもしくはカーネルを、それぞれクラスタリングするクラスタリング処理工程と、前記既知タンパク質のクラスタリング結果および前記既知化合物のクラスタリング結果に従って、前記既知タンパク質と既知化合物との既知相互作用の情報に基づき、前記既知タンパク質と既知化合物に対する前記既知相互作用の分布データを作成する分布データ作成工程とを含むことを特徴とする。
【００１０】
また、本発明のデータ構造は、既知タンパク質と既知化合物との既知相互作用の分布データのデータ構造であって、２以上の既知タンパク質相互の距離を示すプロファイルもしくはカーネルをクラスタリングして得られた既知タンパク質プロファイルもしくはカーネルと、２以上の既知化合物相互の距離を示すプロファイルもしくはカーネルをクラスタリングして得られた既知化合物プロファイルもしくはカーネルと、前記２以上の既知タンパク質のうち１の既知タンパク質に対応する既知タンパク質プロファイルもしくはカーネルと、前記２以上の既知化合物のうち１の既知化合物に対応する既知化合物プロファイルもしくはカーネルとに対応づけて、当該１の既知タンパク質と当該１の既知化合物との相互作用の情報が含まれていることを特徴とする。また、本発明の記録媒体は、本発明のデータ構造を有するデータを記録した記録媒体である。
【発明の効果】
【００１１】
本発明によれば、既知タンパク質のクラスタリング結果および既知化合物のクラスタリング結果に従って、既知タンパク質と既知化合物との既知相互作用の情報に基づき、前記既知相互作用の分布データが作成される。具体的には、例えば、タンパク質側のデータと既知相互作用の情報との対比もしくは化合物側のデータと既知相互作用の情報との対比という、いずれか一方の対比のみではなく、タンパク質相互の対比、化合物相互の対比に加え、さらにタンパク質と化合物の両方に前記既知相互作用の情報を対比させている。このような分布データに基づけば、例えば、目的タンパク質に相互作用する化合物または目的化合物が相互作用するタンパク質の予測を、従来よりも優れた精度で行うことができる。また、従来の既存のデータベースは、タンパク質と化合物間の相互作用を様々な観点から把握する機能はなく、機械学習を利用した相互作用の予測への適用には適していないことからも、本発明は、タンパク質と化合物との間における相互作用に関する極めて新しい手法であるといえる。
【００１２】
さらに、本発明で使用するタンパク質のデータは、タンパク質相互のクラスタリングを行うことができるデータ（タンパク質相互の距離を示すプロファイルもしくはカーネル）であればよく、従来のような機能や立体構造等の情報は必須ではない。また、同様に使用する化合物のデータも、化合物相互のクラスタリングを行うことができるデータ（化合物相互の距離を示すプロファイルもしくはカーネル）であればよく、従来のような活性等に関する情報は必須ではない。このため、従来法では予測できなかった機能や立体構造が未知のタンパク質、活性未知の化合物についても、本発明により既知相互作用に関する分布データを作成すれば、容易に相互作用を示す化合物やタンパク質の予測が可能となる。したがって、本発明のデータ処理装置は、例えば、前述のような医薬品開発やバイオ産業において極めて有用な装置と言える。
【発明を実施するための最良の形態】
【００１３】
本発明のデータ処理装置は、さらに、前記既知タンパク質のクラスタリング結果および既知化合物のクラスタリング結果に従って、前記既知タンパク質と既知化合物とを二軸座標に表し、前記二軸座標に、前記既知相互作用の情報をマッピングし、前記既知相互作用の分布データを表示するマップ表示部を備えることが好ましい。
【００１４】
本発明のデータ処理装置は、既知タンパク質の配列情報や既知化合物の情報をそのまま使用する場合、さらに、２以上の既知タンパク質の配列情報、２以上の既知化合物の情報および前記既知タンパク質と既知化合物との既知相互作用の情報を入力する第１の情報入力部と、前記第１の情報入力部より入力された前記既知タンパク質の配列情報および前記既知化合物の情報を、それぞれ前記タンパク質相互の距離を示すプロファイル（以下、「既知タンパク質プロファイル」ともいう）もしくはカーネル、および前記既知化合物相互の距離を示すプロファイル（以下、「既知化合物プロファイル」ともいう）もしくはカーネルに変換し、前記既知タンパク質プロファイルもしくはカーネル、および、前記既知化合物プロファイルもしくはカーネルを、前記クラスタリング処理部に出力するプロファイル変換部とを備えることが好ましい。一方、前記既知タンパク質プロファイルもしくはカーネル、および、前記化合物プロファイルもしくはカーネルを使用する場合には、さらに、前記既知タンパク質プロファイルもしくはカーネル、前記既知化合物プロファイルもしくはカーネル、および、前記既知タンパク質と既知化合物との既知相互作用の情報を入力する第２の情報入力部を備えることが好ましい。
【００１５】
本発明のデータ処理装置は、さらに、前記既知タンパク質のプロファイルもしくはカーネル、前記既知化合物のプロファイルもしくはカーネル、および、前記分布データ作成部による前記既知相互作用の分布データを教師データとして、タンパク質と化合物との相互作用に関する統計的法則を機械学習し、前記機械学習による統計的法則から得られる結果より、目的タンパク質に対して相互作用を示す化合物または目的化合物が相互作用を示すタンパク質を予測する予測部を備えることが好ましい。このようにタンパク質と化合物との相互作用に関する統計的法則を機械学習すれば、統計的法則から得られる結果により、従来よりも精度よく目的タンパク質に対して相互作用を示す化合物の予測や目的化合物に対して相互作用を示す化合物の予測を行うことができる。
【００１６】
本発明のデータ処理装置において、前記予測部が、機械学習により統計的法則から得られる結果をパラメータとして算出し、前記パラメータに基づいて予測を行うことが好ましい。パラメータの決定により、例えば、新規化合物や新規タンパク質との相互作用を容易に予測できるシステムを構築できる。
【００１７】
本発明のデータ処理装置において、さらに、前記分布データ作成部による前記既知相互作用の分布データより、タンパク質と化合物との相互作用が適切な範囲を選択し、これを教師データとして前記予測部に出力するデータ選択部を備えることが好ましい。このように分布データから適切な範囲を選択することによって、より精度良く予測を行うことができる。
【００１８】
本発明のデータ処理装置において、さらに、予測部による予測結果を出力する予測結果出力部を備えることが好ましい。
【００１９】
本発明のデータ処理装置において、タンパク質の配列情報や化合物の情報をそのまま使用する場合には、さらに、２以上の既知タンパク質の配列情報、２以上の既知化合物の情報および前記既知タンパク質と既知化合物との既知相互作用の情報を入力する第１の情報入力部と、前記第１の情報入力部より入力された前記既知タンパク質の配列情報および前記既知化合物の情報を、それぞれ前記タンパク質相互の距離を示すプロファイルもしくはカーネル、および、前記既知化合物相互の距離を示すプロファイルもしくはカーネルに変換し、前記タンパク質相互の距離を示すプロファイルもしくはカーネル、および、前記既知化合物相互の距離を示すプロファイルもしくはカーネルを、前記クラスタリング処理部に出力するプロファイル変換部と、目的タンパク質の配列情報または目的化合物の情報を入力する第３の情報入力部を備え、前記プロファイル変換部が、さらに、前記第３の情報入力部より入力された目的タンパク質の配列情報を、前記第１の情報入力部より入力された既知タンパク質の配列情報に基づき、前記目的タンパク質の前記既知タンパク質に対する距離を示すプロファイルもしくはカーネルに変換し、または、前記第３の情報入力部により入力された目的化合物の情報を、前記第１の情報入力部により入力された既知化合物の情報に基づき、前記目的化合物の前記既知化合物に対する距離を示すプロファイルもしくはカーネルに変換し、変換された前記目的タンパク質の前記既知タンパク質に対する距離を示すプロファイルもしくはカーネル、または、前記目的化合物の前記既知化合物に対する距離を示すプロファイルもしくはカーネルを前記予測部に出力し、
前記予測部が、前記統計的法則から得られる結果と、前記プロファイル変換部による前記目的タンパク質の前記既知タンパク質に対する距離を示すプロファイルもしくはカーネル、または、前記目的化合物の前記既知化合物に対する距離を示すプロファイルもしくはカーネルとから、目的タンパク質に対して相互作用を示す化合物または前記目的化合物が相互作用を示すタンパク質を予測することが好ましい。このように目的タンパク質を既知タンパク質と同様にプロファイルもしくはカーネルに変換することによって、目的タンパク質または目的化合物についての相互作用を容易に予測することができる。
【００２０】
また、本発明のデータ処理装置において、タンパク質または化合物のプロファイルもしくはカーネルを使用する場合には、さらに、前記既知タンパク質相互の距離を示すプロファイルもしくはカーネル、前記既知化合物相互の距離を示すプロファイルもしくはカーネル、および、前記既知タンパク質と既知化合物との既知相互作用の情報を入力する第２の情報入力部と、目的タンパク質の前記既知タンパク質に対する距離を示すプロファイルもしくはカーネル、または、目的化合物の前記既知化合物に対する距離を示すプロファイルもしくはカーネルを入力する第４の情報入力部とを備え、前記予測部において、前記統計的法則から得られる結果と、前記第４の情報入力部より入力された前記目的タンパク質の前記既知タンパク質に対する距離を示すプロファイルもしくはカーネル、または、目的化合物の前記既知化合物に対する距離を示すプロファイルもしくはカーネルとから、目的タンパク質に対して相互作用を示す化合物または前記目的化合物が相互作用を示すタンパク質を予測することが好ましい。このように目的タンパク質を既知タンパク質と同様に変換したプロファイルもしくはカーネルを使用することによって、目的タンパク質または目的化合物についての相互作用を容易に予測することができる。
【００２１】
本発明のデータ処理プログラムは、２以上の既知タンパク質相互の距離を示すプロファイルもしくはカーネル、および、２以上の既知化合物相互の距離を示すプロファイルもしくはカーネルを、それぞれクラスタリングするクラスタリング処理ステップと、前記既知タンパク質のクラスタリング結果および前記既知化合物のクラスタリング結果に従って、前記既知タンパク質と既知化合物との既知相互作用の情報に基づき、前記既知タンパク質と既知化合物に対する前記既知相互作用の分布データを作成する分布データ作成ステップとを、コンピュータに実行させることを特徴とする。
【００２２】
本発明のデータ処理プログラムは、さらに、前記既知タンパク質のクラスタリング結果および既知化合物のクラスタリング結果に従って、前記既知タンパク質と既知化合物とを二軸座標に表し、前記二軸座標に、前記既知相互作用の情報をマッピングし、前記分布データを表示するマップ表示ステップをコンピュータに実行させることが好ましい。
【００２３】
本発明のデータ処理プログラムは、さらに、２以上の既知タンパク質の配列情報、２以上の既知化合物の情報および前記既知タンパク質と既知化合物との既知相互作用の情報の入力を受け付ける第１の入力ステップと、入力された前記既知タンパク質の配列情報および前記既知化合物の情報を、それぞれ前記タンパク質相互の距離を示すプロファイルもしくはカーネル、および、前記既知化合物相互の距離を示すプロファイルもしくはカーネルに変換するプロファイル変換し、前記タンパク質相互の距離を示すプロファイルもしくはカーネル、および、前記既知化合物相互の距離を示すプロファイルもしくはカーネルを、前記クラスタリング処理部に出力する第１のプロファイル変換ステップとを、コンピュータに実行させることが好ましい。
【００２４】
本発明のデータ処理プログラムは、さらに、前記既知タンパク質相互の距離を示すプロファイルもしくはカーネル、前記既知化合物相互の距離を示すプロファイルもしくはカーネル、および、前記既知タンパク質と既知化合物との既知相互作用の情報の入力を受け付ける第２の入力ステップを、コンピュータに実行させることが好ましい。
【００２５】
本発明のデータ処理プログラムは、さらに、前記既知タンパク質のプロファイルもしくはカーネル、前記既知化合物のプロファイルもしくはカーネル、および、前記分布データ作成ステップにおける前記既知相互作用の分布データを教師データとして、タンパク質と化合物との相互作用に関する統計的法則を機械学習し、前記機械学習による統計的法則から得られる結果より、目的タンパク質に対して相互作用を示す化合物または目的化合物が相互作用を示すタンパク質を予測する予測ステップを、コンピュータに実行させることが好ましい。
【００２６】
本発明のデータ処理プログラムは、前記予測ステップが、前記機械学習により統計的法則から得られる結果をパラメータとして算出し、前記パラメータに基づいて予測を行うステップであることが好ましい。
【００２７】
本発明のデータ処理プログラムは、さらに、前記分布データ作成ステップにおける前記既知相互作用の分布データより、タンパク質と化合物との相互作用が適切な範囲を選択し、これを教師データとして前記予測ステップに出力するデータ選択ステップを、コンピュータに実行させることが好ましい。
【００２８】
本発明のデータ処理プログラムは、さらに、前記予測ステップにおける予測結果を出力する予測結果出力ステップを、コンピュータに実行させることが好ましい。
【００２９】
本発明のデータ処理プログラムは、さらに、２以上の既知タンパク質の配列情報、２以上の既知化合物の情報および前記既知タンパク質と既知化合物との既知相互作用の情報の入力を受け付ける第１の情報入力ステップと、
入力された前記既知タンパク質の配列情報および前記既知化合物の情報を、それぞれ前記タンパク質相互の距離を示すプロファイルもしくはカーネル、および、前記既知化合物相互の距離を示すプロファイルもしくはカーネルに変換し、前記タンパク質相互の距離を示すプロファイルもしくはカーネル、および、前記既知化合物相互の距離を示すプロファイルもしくはカーネルを、前記クラスタリング処理ステップに出力する第１のプロファイル変換ステップと、目的タンパク質の配列情報または目的化合物の情報の入力を受け付ける第３の情報入力ステップと、入力された目的タンパク質の配列情報を、入力された既知タンパク質の配列情報に基づき、前記目的タンパク質の前記既知タンパク質に対する距離を示すプロファイルもしくはカーネルに変換し、または、入力された目的化合物の情報を、入力された既知化合物の情報に基づき、前記目的化合物の前記既知化合物に対する距離を示すプロファイルもしくはカーネルに変換し、変換された前記目的タンパク質の前記既知タンパク質に対する距離を示すプロファイルもしくはカーネル、または、前記目的化合物の前記既知化合物に対する距離を示すプロファイルもしくはカーネルを前記予測ステップに出力する第２のプロファイル変換ステップと、前記統計的法則から得られる結果と、前記目的タンパク質の前記既知タンパク質に対する距離を示すプロファイルもしくはカーネル、または、前記目的化合物の前記既知化合物に対する距離を示すプロファイルもしくはカーネルとから、目的タンパク質に対して相互作用を示す化合物または前記目的化合物が相互作用を示すタンパク質を予測する前記予測ステップとを、コンピュータに実行させることが好ましい。
【００３０】
本発明のデータ処理プログラムは、さらに、前記既知タンパク質相互の距離を示すプロファイルもしくはカーネル、前記既知化合物相互の距離を示すプロファイルもしくはカーネル、および、前記既知タンパク質と既知化合物との既知相互作用の情報の入力を受け付ける第２の情報入力ステップと、目的タンパク質の前記既知タンパク質に対する距離を示す、または、目的化合物の前記既知化合物に対する距離を示すプロファイルもしくはカーネルを入力する第４の情報入力ステップと、前記統計的法則から得られる結果と、入力された前記目的タンパク質の前記既知タンパク質に対する距離を示すプロファイルもしくはカーネル、または、目的化合物の前記既知化合物に対する距離を示すプロファイルもしくはカーネルとから、目的タンパク質に対して相互作用を示す化合物または前記目的化合物が相互作用を示すタンパク質を予測する前記予測ステップとを、コンピュータに実行させることが好ましい。
【００３１】
本発明の記録媒体は、本発明のデータ処理プログラムを格納したコンピュータ読み取り可能な記録媒体である。
【００３２】
本発明のデータ処理方法は、タンパク質と化合物との相互作用に関するデータ処理方法であって、２以上の既知タンパク質相互の距離を示すプロファイルもしくはカーネル、および、２以上の既知化合物相互の距離を示すプロファイルもしくはカーネルを、それぞれクラスタリングするクラスタリング処理工程と、前記既知タンパク質のクラスタリング結果および前記既知化合物のクラスタリング結果に従って、前記既知タンパク質と既知化合物との既知相互作用の情報に基づき、前記既知タンパク質と既知化合物に対する前記既知相互作用の分布データを作成する分布データ作成工程とを含むことを特徴とする。
【００３３】
本発明のデータ処理方法は、さらに、前記既知タンパク質のクラスタリング結果および既知化合物のクラスタリング結果に従って、前記既知タンパク質と既知化合物とを二軸座標に表し、前記二軸座標に、前記既知相互作用の情報をマッピングし、前記分布データを表示するマップ表示工程を含むことが好ましい。
【００３４】
本発明のデータ処理方法は、さらに、２以上の既知タンパク質の配列情報および２以上の既知化合物の情報から、前記タンパク質相互の距離を示すプロファイルもしくはカーネル、および、前記既知化合物相互の距離を示すプロファイルもしくはカーネルに変換する工程を含むことが好ましい。
【００３５】
本発明のデータ処理方法は、さらに、前記既知タンパク質のプロファイルもしくはカーネル、前記既知化合物のプロファイルもしくはカーネル、および、前記分布データ作成工程における前記既知相互作用の分布データを教師データとして、タンパク質と化合物との相互作用に関する統計的法則を機械学習し、前記機械学習による統計的法則から得られる結果より、目的タンパク質に対して相互作用を示す化合物または目的化合物が相互作用を示すタンパク質を予測する予測工程を含むことが好ましい。
【００３６】
本発明のデータ処理方法は、前記予測工程において、機械学習により統計的法則から得られる結果をパラメータとして算出し、前記パラメータに基づいて予測を行うことが好ましい。
【００３７】
本発明のデータ処理方法は、さらに、前記分布データ作成工程における前記既知相互作用の分布データより、タンパク質と化合物との相互作用が適切な範囲を選択し、これを教師データとするデータ選択工程を含むことが好ましい。
【００３８】
本発明のデータ処理方法は、さらに、予測工程における予測結果を出力する予測結果出力工程を含むことが好ましい。
【００３９】
本発明のデータ処理方法は、前記予測工程が、前記統計的法則から得られる結果と前記目的タンパク質の前記既知タンパク質に対する距離を示すプロファイルもしくはカーネル、または、前記目的化合物の前記既知化合物に対する距離を示すプロファイルもしくはカーネルとから、前記目的タンパク質に対して相互作用を示す化合物または前記目的化合物が相互作用を示すタンパク質を予測する工程であることが好ましい。
【００４０】
（第１の実施形態）
本発明のデータ処理装置の実施形態の一例について以下に説明する。図１は、本実施形態におけるデータ処理装置の概略構成を示すブロック図である。なお、以下の実施形態において「距離を示すプロファイル」は「カーネル」であってもよく、他の実施形態においても同様である。
【００４１】
図１に示すように、本実施形態のデータ処理装置は、情報入力部１０１、プロファイル変換部１０２、クラスタリング処理部１０３、分布データ作成部１０４およびマップ表示部１０５を備えている。
【００４２】
このデータ処理装置は、コンピュータによって構築され、プロファイル変換部１０２、クラスタリング処理１０３、分布データ作成部１０４およびマップ表示部１０５は、前記コンピュータのプロセッサが所定のプログラムを実行することによって実現される機能的ブロックである。すなわち、これらの各部は、必ずしもハードウェアとして実装されなくともよい。
【００４３】
情報入力部１０１は、２以上の既知タンパク質の配列情報、２以上の既知化合物の情報および前記既知タンパク質と既知化合物との既知相互作用の情報の入力を受け付ける機能を有する。この情報の入力は、例えば、キーボード、テンキー等の入力デバイスによって行われてもよいし、また、情報が格納されたＣＤ−ＲＯＭ等の記録媒体や、インターネット等の外部通信網にアクセスして入手してもよい。外部通信網から前述の情報を入手する場合、例えば、データベースの情報を入力部にダウンロードしてもよいし、通信回路を接続した状態で前述の情報を利用することもできる。情報元となるデータベース等の例は後述する。
【００４４】
タンパク質の種類は特に制限されないが、前述のように、機能や立体構造等の情報は必須でないことから、本発明は、特に、機能未知のタンパク質や立体構造が未知のタンパク質のデータ処理に有効である。前記立体構造が未知のタンパク質としては、例えば、Ｇタンパク質共役型受容体（ＧＰＣＲ）があげられる。
【００４５】
タンパク質の情報としては、後述するようにタンパク質相互の距離を示すプロファイルへの変換が可能なものであれば特に制限されないが、前述のように立体構造等の情報である必要はなく、例えば、塩基配列やアミノ酸配列（一次配列）等の配列情報があげられる。また、塩基配列やアミノ酸配列から決定されたモチーフ、疎水性指標または二次構造等でもよく、これらの情報は、従来公知のインフォマティクス技術から得られる。情報の種類は、例えば、いずれか一種類の配列情報でもよいし、二以上の配列情報であってもよい。なお、このようなタンパク質情報を後述するようにプロファイル表現することは、これまでの多くのバイオインフォマティクス研究により、タンパク質の構造、機能、物理化学特性を十分に表現し得ることが立証されている。
【００４６】
化合物の情報としては、後述するように化合物相互の距離を示すプロファイルへの変換が可能なものであれば特に制限されず、例えば、化学構造や物理化学特性の情報があげられる。前記化学構造の情報としては、例えば、化合物全体の情報（例えば、化学式、構造式）であってもよいし、部分構造の情報であってもよい。具体的な情報としては、例えば、軌道原子の軌道状態、官能基の種類や数、各種原子（炭素、窒素、硫黄、酸素等）の種類や数、骨格となる炭素の数、環構造の有無や種類、芳香族性の有無等があげられる。また、物理化学特性としては、例えば、疎水性指標、沸点、融点、分子量等があげられ、これらの情報は、実際に測定された情報であってもよいし、化合物の構造から経験的計算によって算出された情報であってもよい。化合物の情報の種類は、例えば、いずれか一種類でもよいが、二種類以上であることが好ましい。
【００４７】
なお、化合物の化学構造を部分構造に分解して表現する手法は、分子構造を計算機上で表現する際に特に有効であり、分子構造の高速構造検索に適用されていることも報告されている。また、化合物の物理化学特性についてのプロファイルも、例えば、薬物の構造活性相関（ＱＳＡＲ）等にも利用され、その実用性が確認されており、化合物と生体系との相互作用を数値化する有効な手法である。
【００４８】
タンパク質の配列情報、化合物の情報、ならびにタンパク質と化合物との相互作用の情報は、各種データベース、文献、辞書等から収集でき、例えば、ＰｕｂＭｅｄ、Ｍｅｒｃｋ ｉｎｄｅｘ、ＭｅＳＨ、ＧＰＣＲＤＢ、ＩＵＰＨＡＲ、ＫＥＧＧ、ＭＤＬ、ＩＳＩＳ等があげられる。
【００４９】
プロファイル変換部１０２は、入力された既知タンパク質の配列情報および既知化合物の情報を、既知タンパク質プロファイルおよび既知化合物プロファイルに変換し、クラスタリング処理部に出力する機能を有する。なお、既知タンパク質プロファイルへの変換は、特に制限されないが、例えば、既存のプログラムであるＢＬＡＳＴ（http://ncbi.nlm.nih.gov/BLAST/）等により実現できる。また、既知化合物プロファイルへの変換は、例えば、分子内の原子タイプを６８種類に分類する方法（Hattori et al. J.AM.CHEM.SOC, Vol.125, No. 39, 2003）に基づき、各化合物について、構造内における前記原子タイプの有無や個数を抽出し、化合物相互の類似性を示す距離計算を行うプログラムによって実現できる。
【００５０】
クラスタリング処理部１０３は、既知タンパク質プロファイルおよび既知化合物プロファイルをクラスタリングする機能を持ち、例えば、Ｒ（フリーソフト）等の一般的な統計解析ツールによって実現できる。
【００５１】
分布データ作成部１０４は、既知タンパク質のクラスタリング結果および既知化合物のクラスタリング結果に従って、前記既知タンパク質と既知化合物との既知相互作用の情報に基づき、前記既知タンパク質と既知化合物に対する前記既知相互作用の分布データを作成する機能を有する。分布データ作成部１０４の機能は、例えば、一般的な統計解析ツールであるＲ（フリーソフト）等によって実現できる。
【００５２】
マップ表示部１０５は、前記既知タンパク質のクラスタリング結果および既知化合物のクラスタリング結果に従って、前記既知タンパク質と既知化合物とを二軸座標に表し、前記二軸座標に、前記既知相互作用の情報をマッピングし、前記分布データを表示する機能を有する。マップ表示部の機能は、例えば、一般的な統計解析ツールであるＲ（フリーソフト）等によって実現できる。前記既知相互作用の情報のマッピングは、二軸座標に対して相互作用の有無を、例えば、色等でプロットすることにより行える。
【００５３】
次に、図２に、本実施形態のデータ処理装置における処理の流れをフローチャートで示す。まず、２以上の既知タンパク質の配列情報、２以上の既知化合物の情報および前記既知タンパク質と既知化合物との既知相互作用の情報の入力（Ｓ１１）から開始される。そして、入力された情報のうち、既知タンパク質の配列情報と既知化合物の情報を、それぞれ前記既知タンパク質のプロファイルおよび前記既知化合物のプロファイルに変換し、得られた既知タンパク質プロファイルと既知化合物プロファイルとの出力を行い（Ｓ１２）、前記既知タンパク質プロファイルおよび既知化合物プロファイルのクラスタリング処理が行われる（Ｓ１３）。続いて、得られた既知タンパク質のクラスタリング結果、既知化合物のクラスタリング結果に従って、前記既知タンパク質と既知化合物との既知相互作用の情報に基づき、前記既知タンパク質と既知化合物に対する前記既知相互作用の分布データを作成する処理（Ｓ１４）を行う。さらに、前記既知タンパク質のクラスタリング結果および既知化合物のクラスタリング結果に従って、前記既知タンパク質と既知化合物とを二軸座標に表し、前記二軸座標に、前記既知相互作用の情報をマッピングし、前記分布データを表示するマップ表示処理（Ｓ１５）を行う。
【００５４】
このように分布データのマップにより、既知タンパク質相互の関係、既知化合物相互の関係、既知タンパク質と既知相互作用の関係、および既知化合物と既知相互作用の関係が表示されるため、肉眼でも容易に結果を把握でき、タンパク質と化合物間との相互作用に関する統計的法則を、視覚を通じた判断により抽出することができる。従って、統計的法則の結果に基づけば、例えば、目的タンパク質に対して相互作用を示す化合物、目的化合物に対して相互作用を示すタンパク質を予測することができる。
【００５５】
なお、既知タンパク質の配列情報や既知化合物の情報を直接入力するのではなく、予め前記各情報から既知タンパク質プロファイルと既知化合物プロファイルを算出しておき、これらのプロファイルを直接入力してもよい。この場合、情報入力部１０１およびプロファイル変換部１０２に代えて、既知タンパク質プロファイル、既知化合物プロファイルおよび既知相互作用の情報を入力する情報入力部を備えてもよい。この場合、例えば、Ｓ１１およびＳ１２を省き、前記プロファイルを入力するステップから、クラスタリング処理のステップに移ってもよい。
【００５６】
また、マップ表示部１０５は必須の構成ではなく、例えば、分布データ作成部により得られた分布データを、そのまま表示してもよいし、前記分布データを外部に出力してもよい。
【００５７】
既知タンパク質の配列情報、既知化合物の情報および前記既知相互作用の情報を入力することによって、これらの情報を格納したデータベースを、本実施形態のデータ処理装置内に構築してもよい。また、このデータベースに格納した情報は、例えば、データ処理装置が表示ユニットを備える場合、その表示画面に表示できることが好ましい。このため、これらの情報、特に化合物の構造に関する情報は、ディスプレイ用のイメージ形式のファイル（例えば、ＧＩＦ、ＪＰＧ等）をコンピュータが認識できるように数値化したファイル（例えば、ＭＯＬ、ＳＭＩＬＥ等）として格納することが好ましい。
【００５８】
さらに、本実施形態のデータ処理装置において、既知タンパク質の情報、既知化合物の情報および既知相互作用は、格納されてプロファイルへの変換に使用されるだけでなく、例えば、情報を検索し、その結果を表示できるように設計されていることが好ましい。このようにすれば、例えば、格納した既知タンパク質情報、既知化合物情報および既知相互作用情報を様々な方法で検索し、いわゆるケモゲノミクスに有用な結果を読み出すことができる。この場合、これらの情報が、データベースにおいてタンパク質カテゴリと化合物カテゴリの２つのカテゴリから構成され、既知相互作用の情報を、前記両カテゴリ間のクロスリンクとして表示されることが好ましい。このような設定であれば、タンパク質-化合物間のクロス検索も可能となり、タンパク質側からも化合物側からもより簡単にダイナミックな検索が可能となる。また、あわせて分布データも表示できるように設計されることが好ましい。
【００５９】
なお、化合物に関する入力する情報としては、前述のような構造や物理化学的特性等の他に、例えば、ＣＡＳ登録番号、タンパク質リンク情報、そのタンパク質との生物活性情報（アゴニスト、アンタゴニスト、インヒビター等）、データベースリンク情報等があげられる。また、タンパク質に関する入力情報としては、前述のような配列情報の他に、化合物リンク情報、その化合物とのＩＤ、名称、データベースリンク情報等もあわせて入力してもよい。
【００６０】
（第２の実施形態）
本発明のデータ処理装置の実施形態のさらに他の例について以下に説明する。図３は、本実施形態のデータ処理装置の概略構成を示すブロック図であり、前記第１の実施形態における構成と同様の機能を有する構成については、前記第１の実施形態と同じ符号を付している。
【００６１】
本実施形態のデータ処理装置は、第１の実施形態におけるデータ処理装置に加えて、さらに予測部１０６を備えており、既知タンパク質プロファイル、既知化合物プロファイルおよび分布データに基づき、目的タンパク質に対して相互作用を示す化合物または目的化合物が相互作用を示すタンパク質を予測する装置の形態である。
【００６２】
予測部１０６は、プロファイル変換部１０２における前記既知タンパク質プロファイルならびに前記既知化合物プロファイル、および、分布データを教師データとして、タンパク質と化合物との相互作用に関する統計的法則を機械学習し、前記機械学習による統計的法則から得られる結果により、目的タンパク質に対して相互作用を示す化合物または目的化合物が相互作用を示すタンパク質を予測する機能を有する。機械学習は、特に制限されないが、一般的なサポートベクタマシン、アスペクトモデル等の分類機によって実現できる。
【００６３】
また、予測部１０６は、機械学習により統計的法則から得られる結果をパラメータとして算出し、前記パラメータに基づいて予測を行う機能を有してもよい。
【００６４】
また、本実施形態のデータ処理装置は、さらに、プロファイル変換部１０２による既知タンパク質プロファイルならびに既知化合物プロファイルおよび分布データ作成部１０４による既知相互作用の分布データより、タンパク質と化合物との相互作用が適切な範囲を選択し、これを教師データとして予測部６０１に出力するデータ選択部（図示せず）を備えてもよい。これによって、より適切な教師データが与えられるため、より信頼性が向上した機械学習による統計的法則の結果を得ることができる。また、既知タンパク質プロファイル、既知化合物プロファイルおよび分布データからの教師データの選択は、マップの結果に基づいて行ってもよい。
【００６５】
また、本実施形態のデータ処理装置は、さらに、予測部１０６による予測結果を出力する予測結果出力部を備えてもよい。出力した結果は、例えば、装置が表示ユニットを備える場合には、その表示画面に表示されてもよいし、プリントアウト等により外部に表示されてもよく、また、各種記録媒体に格納されてもよい。
【００６６】
次に、図４に、本実施形態のデータ処理装置における処理の流れをフローチャートで示す。なお、前述の実施形態と同一の処理工程については説明を省く。
【００６７】
既知相互作用の分布データを教師データとして、タンパク質と化合物との相互作用に関する統計的法則を、教師あり機械学習処理し、機械学習から得られる結果より、目的タンパク質に対して相互作用を示す化合物または目的化合物の相互作用を示すタンパク質を予測する予測処理（Ｓ１６）を行う。前記予測処理は、例えば、前記機械学習により統計的法則から得られる結果をパラメータとして算出し、前記パラメータに基づいて行ってもよい。また、得られた既知タンパク質プロファイル、既知化合物プロファイルおよび前記既知相互作用の分布データより、タンパク質と化合物との相互作用が適切な範囲を選択し、これを教師データとして予測処理してもよい。さらに、得られた予測結果を出力してもよい。
【００６８】
（第３の実施形態）
本発明のデータ処理装置の実施形態のさらに他の例について以下に説明する。図５は、本実施形態のデータ処理装置の概略構成を示すブロック図であり、前述の実施形態における構成と同様の機能を有する構成については、同じ符号を付している。
【００６９】
本実施形態のデータ処理装置は、第２の実施形態におけるデータ処理装置と同様の構成であるが、情報入力部１０１は、さらに、目的タンパク質の配列情報または目的化合物の情報の入力を受け付ける機能を有し、プロファイル変換部１０２は、さらに、入力された目的タンパク質の配列情報を、情報入力部１０１より入力された既知タンパク質の配列情報に基づき、前記目的タンパク質の前記既知タンパク質に対する距離を示すプロファイル（以下、「目的タンパク質プロファイル」という）に変換し、または、情報入力部１０１により入力された目的化合物の情報を、情報入力部１０１により入力された既知化合物の情報に基づき、前記目的化合物の前記既知化合物に対する距離を示すプロファイル（以下、「目的化合物プロファイル」という）に変換し、変換された前記目的タンパク質プロファイルまたは前記目的化合物プロファイルを予測部１０６に出力する機能を有する。また、予測部１０６は、前記統計的法則から得られる結果と、前記目的タンパク質プロファイルまたは前記目的化合物プロファイルとから、目的タンパク質に対して相互作用を示す化合物または前記目的化合物が相互作用を示すタンパク質を予測する機能を有する。
【００７０】
次に、図６に、本実施形態のデータ処理装置における処理の流れをフローチャートで示し、その具体例を図７〜１１を用いて説明する。
【００７１】
まず、２以上の既知タンパク質の配列情報、２以上の既知化合物の情報および前記既知タンパク質と既知化合物との既知相互作用の情報の入力（Ｓ１１）から開始される。入力によって、これらの情報を保存したデータベースを、本実施形態のデータ処理装置内に構築してもよい。図７に、このデータベースの構築に関する概略を示す。同図に示すように、例えば、タンパク質情報、化合物情報、相互作用の実験情報を、それぞれ既存のデータベースや文献等から抽出し、これらの入力によってデータベース（ＤＢ）を構築する。なお、同図において「Ｐ００００１、Ｐ００００２、・・・」とはタンパク質の種類を示し、「Ｌ００００１、Ｌ００００２、・・・」とは化合物の種類を示し、タンパク質−化合物間の矢印は相互作用の関係がある組合せを示す。
【００７２】
そして、入力された情報のうち、既知タンパク質の配列情報と既知化合物の情報を、それぞれ前記既知タンパク質のプロファイルおよび前記既知化合物のプロファイルに変換し、得られた既知タンパク質プロファイルと既知化合物プロファイルとの出力を行う（Ｓ１２）。図８に、データベースの既知タンパク質情報からプロファイル変換を行うステップの概略を示す。同図に示すように、データベースのタンパク質データから、例えば、遺伝子配列データを抽出し、アミノ酸配列に置き換えてから、タンパク質相互の相同性を既知プログラム（例えば、ＢＬＡＳＴ）により、プロファイルに変換する。前記ＢＬＡＳＴを利用する場合、通常、配列の一致度が高いほど、相互の関係に高いスコアが与えられる。このようにスコアを算出することによって、同図に示すように、タンパク質間における相同性を示す距離行列が得られる。なお、得られたスコアは、例えば、距離に変換してもよい。
【００７３】
また、図９に、データベースの化合物情報からプロファイル変換を行うステップの概略を示す。同図に示すように、データベースの化合物データから、例えば、化学構造および物理化学特性を抽出し、数値化する。化学構造や物理化学特性に関しては、例えば、ＭＤＬデータベース等のＭＯＬ形式からプロファイルへの変換を行う。この変換は、例えば、前述の方法（Hattori et al. J.AM.CHEM.SOC, Vol.125, No. 39, 2003）に基づいて行うことができる。この数値化によって、化合物は、ｎ個の数値化された特徴（Chemical Descriptor）からなるｎ次ベクトルで表現できる。そして、このベクトル空間における距離を算出することによって、各化合物間の類似度を示す距離行列が算出できる。
【００７４】
続いて、既知タンパク質プロファイルおよび既知化合物プロファイルを、それぞれクラスタリング処理する（Ｓ１３）。このクラスタリング処理は、例えば、前述のような距離行列の結果に基づいて行うことができる。
【００７５】
続いて、得られた既知タンパク質のクラスタリング結果、既知化合物のクラスタリング結果に従って、前記既知タンパク質と既知化合物との既知相互作用の情報に基づき、前記既知タンパク質と既知化合物に対する前記既知相互作用の分布データを作成する処理（Ｓ１４）を行う。さらに、前記既知タンパク質のクラスタリング結果および既知化合物のクラスタリング結果に従って、前記既知タンパク質と既知化合物とを二軸座標に表し、前記二軸座標に、前記既知相互作用の情報をマッピングし、前記分布データを表示するマップ表示処理（Ｓ１５）を行ってもよい。
【００７６】
マッピングした分布データの概略を図１０に示す。同図に示すように、既知タンパク質のプロファイル（距離行列）および既知化合物のプロファイル（距離行列）をクラスタリング処理し、それぞれの系統樹を横軸（タンパク質）、縦軸（化合物）に表した二軸座標のマップを作成する。そして、このマップに、前述の既知相互作用を示す既知タンパク質と既知化合物との組合せについてプロットを行う。同図において、黒塗り部分が、相互作用が既知である座標を示す。
【００７７】
他方、相互作用を予測する目的タンパク質もしくは目的化合物について、タンパク質の配列情報または目的化合物の情報を入力する処理を行う（Ｓ２１）。これらの情報に基づいて、次のステップにおいて既知タンパク質または既知化合物に対するプロファイルへの変換を行うことから、目的タンパク質の配列情報と目的化合物の情報とは、すでに入力した既知タンパク質や既知化合物と同じ種類の情報であることが望ましい。
【００７８】
続いて、入力された目的タンパク質の配列情報を、入力された既知タンパク質の配列情報に基づき、前記目的タンパク質の前記既知タンパク質に対する距離を示すプロファイルに変換し、または、入力された目的化合物の情報を、入力された既知化合物の情報に基づき、前記目的化合物の前記既知化合物に対する距離を示すプロファイルに変換する（Ｓ２２）。
【００７９】
そして、前記既知タンパク質プロファイル、既知化合物プロファイルおよび前述の分布データと、前記目的タンパク質プロファイルまたは目的化合物プロファイルとから、目的タンパク質に対して相互作用を示す化合物または前記目的化合物が相互作用を示すタンパク質を予測する（Ｓ１６）。すなわち、目的タンパク質プロファイルまたは目的化合物プロファイルと、既知タンパク質プロファイル、既知化合物プロファイルおよび分布データを教師データとする機械学習による相互作用に関する統計的法則の結果とから、相互作用を示す化合物またはタンパク質を予測できる。
【００８０】
図１１に、マップを用いた予測処理の概略を示す。例えば、目的タンパク質「Ｐ００？？？」のプロファイルを前記二軸座標マップに加えると、例えば、系統樹における目的タンパク質の位置決定、目的タンパク質と相同性のある既知タンパク質群の選択、前記タンパク質群と相互作用を示す既知化合物の選択、前記相互作用を示す既知化合物と類似性を示す既知化合物群の選択等が行われる。そして、これらの結果と、統計的法則から得られる結果から、目的タンパク質「Ｐ００？？？」と相互作用を示す化合物が予測される。また、目的化合物の場合、目的化合物「Ｌ００？？？」のプロファイルを前記二軸座標マップに加えると、例えば、系統樹における目的化合物の位置決定、目的化合物と類似性を示す既知化合物群の選択、前記既知化合物群が相互作用を示す既知タンパク質の選択、前記相互作用を示す既知タンパク質と類似性を示す既知タンパク質群の選択等が行われる。そして、これらの結果と統計的法則から得られる結果から、目的化合物「Ｐ００？？？」が相互作用を示すタンパク質が予測される。
【００８１】
なお、本発明は、これらの実施形態には何ら制限されず、例えば、前述のような処理をコンピュータに実行させるプログラムや、前記プログラムを格納したコンピュータ読み取り可能な記録媒体、データ構造も、本発明の対象である。
【実施例１】
【００８２】
１．情報の抽出
既知タンパク質として、ＧＰＣＲｓを選択した。ＧＰＣＲのアミノ酸配列情報は、ＧＰＣＲＤＢ、ＩＵＰＨＡＲ、ＫＥＧＧ、Ｌｏｃｕｓ ｌｉｎｋ等のリンク情報より収集した。一方、各種ＧＰＣＲｓのリガンド化合物の情報は、ＰｕｂＭｅｄ、Ｍｅｒｃｋ ｉｎｄｅｘ、ＭｅＳＨ、ＧＰＣＲＤＢ、ＩＵＰＨＡＲ、ＫＥＧＧ等の文献、辞書、公的データベースから収集した。化合物の情報は、化学構造をＧＩＦイメージファイルとして、また、構造の数値化をＭＯＬ形式ファイルとして格納した。これらの抽出情報について、化合物の名称、分子量、ＣＡＳ登録番号、化学式（formula）を格納した化合物情報テーブルと、ＧＰＣＲｓの前記リンク情報を格納したＧＰＣＲｓ情報テーブルを作成した。そして、ＧＰＣＲとリガンド化合物の相互作用の情報は、相互に関係があるＧＰＣＲと化合物とを整列した関係情報テーブルに格納した。なお、本発明におけるデータベースとして好ましい構成と機能との概略を図１２に示す。
【００８３】
２．データベースの基本構造の決定
前記３種類の情報テーブルと化合物の構造ファイルとを利用して、ＧＰＣＲとリガンド化合物のデータベースの基本構造を構築した。図１３は、前記データベースを構成する２つ情報テーブルが関係情報テーブルによって連結されている状態を表す。収集した各ＧＰＣＲｓ情報と化合物情報は、それぞれ固有の結果ページを持つ。ＧＰＣＲの結果ページとして、ＧＰＣＲ情報テーブルには、ファミリーの情報や公的データベースへのリンク情報がディスプレイされ、関係情報テーブルには、そのＧＰＣＲと相互作用するリガンド化合物とその情報がディスプレイされる。この結果ページに表示されるリガンド化合物とＧＰＣＲの職別番号は、それぞれ個々の結果ページにハイパリンク（hyperlink）されているため、相互関係があるＧＰＣＲとリガンド化合物間のダイナミックな移動と検索が可能である。また、図１４は、本発明により制作したテータベースの結果を示す一例であり、検索で読み出された化合物の結果ページと、この化合物と相互作用するタンパク質の一つの結果ページとのリンク関係を表している。
【００８４】
３．化合物とＧＰＣＲｓの分類テーブルの作成
このデータベースの検索は、例えば、キーワード検索とともに、ＧＰＣＲの分類テーブルやリガンドの分類テーブルから行うこともできる。この分類テーブルは、プロファイルに従って、化合物またはＧＰＣＲｓをグルーピングしたものである。この分類は、既存のものを利用してもよいし、既存の分類を利用して、対象として好ましい性質の単位で再度分類したものを使用することもできる。
【００８５】
４．化合物のプロファイルによる分類テーブル
分子内の原子タイプを６８種類に分類する方法（Hattori et al. J.AM.CHEM.SOC, Vol.125, No. 39, 2003）に基づき、各化合物について、構造内における前記原子タイプの有無や個数を抽出し、化合物相互の類似性を示す距離計算を行った。この距離結果に基づいて各化合物をクラスタリングし、クラスタリング結果を化合物の系統樹（tree）として画面に表示して、その末端を各化合物の結果ページにリンクした。
【００８６】
５．ＧＰＣＲｓの分類テーブル
ＧＰＣＲＤＢ（http://www.gpcr.org/7tm）の分類テーブルを導入してＧＰＣＲｓの分類テーブルを作成し、クラスタリング結果をＧＰＣＲの系統樹（tree）として画面に表示し、その末端を各ＧＰＣＲの結果ページにリンクした。
【００８７】
６．類似化合物リストの表示
前記化合物の結果ページに、前記距離計算による類似度が近い化合物のリストを、その情報（例えば、構造イメージ等を含む）とともに表示した。
【００８８】
７．タンパク質と化合物の相互作用の表示
前記化合物の結果ページに、前記距離計算により距離が近い２０個の化合物を横軸、これらの化合物と相互作用の関係を示す４０個のＧＰＣＲｓを縦軸として、二軸座標でマッピングし、グラフィクで表示した。化合物に関する横軸には、前記分類方法に基づいてクラスタリングしたtreeを、ＧＰＣＲｓに関する縦軸には、ＢＬＡＳＴを用いて計算した各ＧＰＣＲの配列類似性によりクラスタリングしたtreeを表示した。そして、この二軸座標を示すマップには、化合物とＧＰＣＲとの相互作用の有無を色で表示し、肉眼で容易に把握できるようにした。なお、化合物間およびＧＰＣＲ間の距離計算ならびにとマッピングは、統計解析ツールの一種であるＲ（フリーソフト）を使用した。
【００８９】
８．教師データの決定
タンパク質とリガンド化合物との様々な組合せについて、前記方法により相互作用の結果を分析した。特に、相互作用を示す化合物とＧＰＣＲの中で、化合物間の距離とＧＰＣＲ間の距離にある傾向が存在すると思われるデータの組合せのみを選び、これを教師データとした。
【００９０】
９．機械学習によるパラメータの決定
機械学習の分類ツールであるサポートベクタマシンを用い、前段階で選択した教師データについて、教師あり機械学習を行い、学習パラメータを決定した。
【産業上の利用可能性】
【００９１】
このように本発明によれば、例えば、従来のように経験的力場関数を用いたエネルギー最適化を基礎とすることなく、活性未知な化合物や、機能未知または立体構造が未知であるタンパク質に対しても適用可能な、タンパク質と化合物との相互作用に関するデータ処理を行うことができる。また、既知タンパク質の情報、既知化合物の情報ならびに既知タンパク質と既知化合物との既知相互作用の情報に基づいて相関関係に関するデータを処理するため、目的タンパク質や目的化合物と相互作用を示す化合物やタンパク質の予測を精度良く行うことができる。このため、例えば、医薬品開発やバイオ産業において極めて有用な装置と言える。
【図面の簡単な説明】
【００９２】
【図１】図１は、本発明の実施形態におけるデータ処理装置の概略構成を示すブロック図である。
【図２】図２は、本発明の実施形態のデータ処理装置における処理の流れを示すフローチャートである。
【図３】図３は、本発明のその他の実施形態のデータ処理装置の概略構成を示すブロック図である。
【図４】図４は、本発明のその他の実施形態のデータ処理装置における処理の流れを示すフローチャートである。
【図５】図５は、本発明のその他の実施形態のデータ処理装置の概略構成を示すブロック図である。
【図６】図６は、本発明のその他の実施形態のデータ処理装置における処理の流れを示すフローチャートである。
【図７】図７は、本発明のその他の実施形態におけるデータベースの構築を示す概略図である。
【図８】図８は、本発明のその他の実施形態における、データベースの既知タンパク質情報からプロファイル変換を行うステップの概略を示す図である。
【図９】図９は、本発明のその他の実施形態における、データベースの化合物情報からプロファイル変換を行うステップの概略を示す図である。
【図１０】図１０は、本発明のその他の実施形態における、マッピングした分布データの概略図である。
【図１１】図１１は、本発明のその他の実施形態における、マップを用いた予測処理の概略を示す図である。
【図１２】図１２は、本発明の実施例における、データベースの構成と機能とを示す概略図である。
【図１３】図１３は、本発明の実施例における、前記データベースを構成する２つの情報テーブルと関係情報テーブルとの連結状態を示す概略図である。
【図１４】図１４は、本発明の実施例における、テータベースの結果の一例を示す図である。
【図１】
2007011752.tif 000003
【図２】
2007011752.tif 000004
【図３】
2007011752.tif 000005
【図４】
2007011752.tif 000006
【図５】
2007011752.tif 000007
【図６】
2007011752.tif 000008
【図７】
2007011752.tif 000009
【図８】
2007011752.tif 000010
【図９】
2007011752.tif 000011
【図１０】
2007011752.tif 000012
【図１１】
2007011752.tif 000013
【図１２】
2007011752.tif 000014
【図１３】
2007011752.tif 000015
【図１４】
2007011752.tif 000016