抗原と抗体の複合体立体構造の予測構造である抗原抗体複合体候補構造の適正を評価する抗原抗体複合体候補構造評価装置であって、前記抗原と前記抗体の反応速度定数である結合速度定数あるいは解離速度定数の少なくともどちらか一方を取得する反応速度情報取得部と、前記抗原抗体複合体候補構造における抗原抗体間の静電的相互作用の強さを表す静電作用スコアを取得する構造的作用スコア取得部と、前記解離速度定数または前記結合速度定数と、前記静電作用スコアとを用いて、前記抗原抗体複合体候補構造の適正を総合評価する複合体評価スコアを算出する総合スコア算出部とを備える抗原抗体複合体複合体候補構造評価装置。前記反応速度情報取得部が前記結合速度定数と前記結合速度定数を取得し、前記総合スコア算出部が前記結合速度定数と前記結合速度定数と、前記構造的作用スコア取得部から得られた前記静電作用スコアとを用いて、前記抗原抗体複合体候補構造の適正を総合評価する複合体評価スコアを算出する請求項1記載の抗原抗体複合体候補構造評価装置。前記抗原を構成するアミノ酸残基の数である抗原アミノ酸残基数を取得する単体構造情報取得部を備え、前記総合スコア算出部が、前記解離速度定数と前記抗原アミノ酸残基数から前記抗原と前記抗体が結合した複合体立体構造の結合部位面積の予測値である面積予測スコアを算出する請求項1または請求項2記載の抗原抗体複合体候補構造評価装置。前記抗原を構成するアミノ酸残基の数である抗原アミノ酸残基数と、前記抗体のCDR領域超可変領域の極性状態を示す抗体CDR極性スコアを取得する単体構造情報取得部を備え、前記総合スコア算出部が、前記結合速度定数と前記抗原アミノ酸残基数と前記抗体CDR極性スコアから、前記抗原と前記抗体の間の静電的な相互作用の強さを予測した値である静電予測スコアを算出する請求項1または請求項2記載の抗原抗体複合体候補構造評価装置。前記総合スコア算出部が、前記面積予測スコアと前記静電予測スコアを積分して、前記抗原と前記抗体の静電的な相互作用状態を表す実験スコアを算出し、前記静電作用スコアが前記実験スコアより大きい場合に前記静電作用スコアを前記実験スコアに置換えた補正後静電作用スコアを算出する請求項3および請求項4記載の抗原抗体複合体候補構造評価装置。前記構造的作用スコア取得部が、抗原抗体複合体候補構造における抗原抗体間の形状相補性による相互作用の強さを表す形状相補作用スコアを取得し、前記総合スコア算出部が、前記補正後静電作用スコアと形状相補作用スコアを用いて前記複合体評価スコアを算出する請求項5記載の抗原抗体複合体候補構造評価装置。コンピュータにより、抗原と抗体の複合体立体構造の予測構造である抗原抗体複合体候補構造の適正を評価する抗原抗体複合体候補構造評価方法であって、前記抗原と前記抗体の反応速度定数である結合速度定数と解離速度定数を取得し、前記抗原抗体複合体候補構造における抗原抗体間の静電的相互作用の強さを表す静電作用スコアを取得し、前記抗原抗体複合体候補構造における抗原抗体間の形状相補性による相互作用の強さを表す形状相補作用スコアを取得し、前記抗原を構成するアミノ酸残基の数である抗原アミノ酸残基数と、前記抗体のCDR領域超可変領域の極性状態を示す抗体CDR極性スコアを取得し、前記解離速度定数と前記抗原アミノ酸残基数から前記抗原と前記抗体が結合した複合体立体構造の結合部位面積の予測値である面積予測スコアを算出し、前記結合速度定数と前記抗原アミノ酸残基数と前記抗体CDR極性スコアから、前記抗原と前記抗体の間の静電的な相互作用の強さを予測した値である静電予測スコアを算出し、前記面積予測スコアと前記静電予測スコアを積分して、前記抗原と前記抗体の静電的な相互作用状態を表す実験スコアを算出し、前記静電作用スコアが前記実験スコアより大きい場合に前記静電作用スコアを前記実験スコアに置換えた補正後静電作用スコアを算出し、前記補正後静電作用スコアと形状相補作用スコアを用いて前記複合体評価スコアを算出する抗原抗体複合体候補構造評価方法。請求項7に記載の相互作用力変化予測方法をコンピュータに実行させるためのプログラム。しかしながら、アミノ酸残基やアミノ酸残基ペア単位での指標を生成する際にはアミノ酸残基単位へのまるめこみによる評価値誤差が生じ、立体構造情報の原子配置に基づく指標においては各タンパク質の構造精度や構造的な相補性評価方法による誤差などが生じ、既知のタンパク質間の結合情報に基づく指標においては統計データの偏りやばらつきによる誤差が生じる。前記従来の構成では、各指標で生じている誤差を補うためにこれらの指標を組み合わせることで精度向上を図っているものの、評価精度はまだ十分であるとはいえない。本発明は、前記従来の課題を解決するもので、評価する複合体の構造情報や既知のタンパク質の結合情報である統計情報だけではなく、結合反応状態を正しく示している実験値を用いることで、評価精度の高い複合体候補構造の適正評価方法を提供することを目的とする。前記従来の課題を解決するために、本発明の抗原抗体複合体候補構造の評価方法は、前記抗原抗体複合体を構成する抗原と抗体の反応速度定数である結合速度定数あるいは解離速度定数の少なくともどちらか一方を取得する反応速度情報取得部と、前記抗原抗体複合体候補構造における抗原抗体間の静電的相互作用の強さを表す静電作用スコアを取得する構造的作用スコア取得部と、前記反応速度情報取得部で得られた前記解離速度定数または前記結合速度定数と、前記構造的作用スコア取得部から得られた前記静電作用スコアとを用いて、前記抗原抗体複合体候補構造の適正を総合評価する複合体評価スコアを算出する総合スコア算出部を備える。本構成によって、結合反応状態を正しく示している実験値を用いて、構造的に評価を行った評価スコアの補正や再評価を行うことで、評価精度の高い複合体候補構造の適正評価方法を提供することができる。本発明の抗原抗体複合体候補構造の評価方法によれば、評価する複合体候補構造や統計情報だけでなく、結合反応状態を誤りなく示している実験値を用いることで、精度の高い複合体候補構造の適正評価方法を提供することができる。以下本発明の実施の形態について、図面を参照しながら説明する。実施の形態1図1は、本発明の実施の形態1における抗原抗体複合体候補構造の評価装置の全体構成を示す図である。図1において、抗原抗体複合体候補構造情報502は、複合体を構成する抗原と抗体の各立体構造を入力情報として、2つのタンパク質が相互作用したときの複合体立体構造を予測するドッキング予測プログラムの出力結果として得られる複合体候補構造情報である。例えば、ドッキング予測プログラムとして、非特許文献1の複合体立体構造予測の手法と、非特許文献2の複合体候補構造スコア評価手法に基づくsurFitがある。ドッキング予測プログラムsurFitのURLを示す。http://sysimm.ifrec.osakau.ac.jp/surFit/抗原抗体反応速度情報501は、SPR表面プラズモン共鳴装置などを用いて得られた抗原と抗体の反応速度定数である結合速度定数と解離速度定数の実験値である。抗原抗体複合体候補構造情報502は、ドッキング予測プログラムなどにより予測された抗原と抗体の複合体の予測構造である。複合体候補構造評価装置201は、反応速度情報取得部211と構造的作用スコア取得部212と単体構造情報取得部213と統合スコア算出部221とを備える。反応速度情報取得部211は、抗原抗体反応速度情報501から抗原と抗体の結合反応実験により得られた結合速度定数や解離速度定数を取得する。構造的作用スコア取得部212は、抗原抗体複合体候補構造情報502から、抗原と抗体間の結合部位における静電的な相互作用状態を表す静電作用スコアと、抗原と抗体間の結合部位の形状的な相補性の高さを表す形状相補作用スコアを取得する。具体的には、まず、抗原抗体複合体候補構造情報502から、非特許文献2においてthecomplementarityscorefortheelectrostaticpotentialとして示されている、抗原と抗体間の静電的な相互作用の強さを表す静電ポテンシャルスコアと、非特許文献2においてthecomplementarityscorefortheshapeとして示されている、抗原と抗体間の形状相補的な相互作用の強さを表す形状相補ポテンシャルスコアと、非特許文献2においてはvertex-pairsexistingintheinterfacesで示される頂点で囲まれた面の面積となる、結合部位面積とを取得する。次に、取得した静電ポテンシャルスコアと結合部位面積とを積算して静電作用スコアを算出し、形状相補ポテンシャルスコアと結合部位面積とを積算して形状相補作用スコアを算出する。単体構造情報取得部213は、抗原抗体複合体候補構造情報502から得られる抗原のアミノ酸配列情報から抗原アミノ酸残基数を取得する。さらに、抗体アミノ酸配列情報から抗原結合領域であるCDR領域超可変領域のアミノ酸配列情報を取り出し、アミノ酸残基タイプ毎の極性を示す値として、例えば非特許文献3に示される極性スコアを用いて、CDR領域を構成する各アミノ酸残基の極性スコアの平均値を求め、抗体CDRの極性状態を示す抗体CDR極性スコアを算出する。次に、統合スコア算出部221が、反応速度情報取得部211から得られた結合速度定数と解離速度定数、構造的作用スコア取得部212から得られた静電作用スコアと形状相補作用スコア、個別構造情報取得部213から得られた抗原アミノ酸残基数と抗体CDR極性スコアを用いて複合体評価スコア503を算出する方法について、図2を用いて説明する。図2は総合スコア算出部が実行する処理フローを示す図である。総合スコア算出部221は結合速度定数508と抗原アミノ酸残基数513と抗体CDR極性スコア514から、静電ポテンシャルスコアの予測値である静電予測スコアを、以下の式1に基づいて算出するS20。C11、C12、C13、C14は、静電的相互作用の強さに対しての結合速度定数508と抗原アミノ酸残基数513と抗体CDR極性スコア514のそれぞれの貢献度に応じた係数と定数である。これらの係数と定数は以下のように算出できる。PDBProteinDataBankに抗原抗体複合体構造が登録され、かつ結合速度定数の実験値が求められているデータを複数抽出する。抽出した抗原抗体複合体構造について非特許文献2に示される手法を用いて静電ポテンシャルスコアを算出し、抗体CDRのアミノ酸配列情報から抗体CDR極性スコアを算出する。算出した静電ポテンシャルスコアを目的変数、結合速度定数実験値と抗原アミノ酸残基数、算出した抗体CDR極性スコアを説明変数として回帰分析を行い、回帰分析の結果はC11=9.9、C12=0.24、C13=1.8、C14=-27.4、であった。結合が早い場合は静電ポテンシャルが大きく引きあう力が強くなっていると考えられ、抗原アミノ酸残基数が大きいほど抗原全体の電荷が強くなって静電的に引き合う力が強くなり、抗体CDR極性スコアが大きいほど静電的作用の強さが大きいと考えられる。なお、各係数は異なる値でもよい。次に、解離速度定数509と抗原アミノ酸残基数513から、結合部位面積の予測値である面積予測スコアを、以下の式2に基づいて算出するS21。C21、C22、C23は、結合部位面積に対しての解離速度定数509と抗原アミノ酸残基数513のそれぞれの貢献度に応じた係数と定数である。これらの係数と定数は以下のように算出できる。PDBに抗原抗体複合体構造が登録され、かつ解離速度定数の実験値が求められているデータを複数抽出する。抽出した抗原抗体複合体構造について非特許文献2に示される手法を用いて結合部位面積を算出する。算出した結合部位面積を目的変数、解離速度定数実験値と抗原アミノ酸残基数を説明変数として回帰分析を行い、回帰分析の結果はC21=-90.2、C22=1.3、C23=366、であった。解離が遅い場合は結合部位の面積が大きくなることで静電的あるいは疎水的な相互作用が大きくなっている状態と考えられ、抗原アミノ酸残基数が大きいほど結合部位面積が大きい傾向となる。なお、各係数は異なる値でもよい。次に、これらの算出した静電予測スコアと面積予測スコアを積分して、実験値を用いて求めた抗原抗体の静電的な相互作用状態を表す実験スコアの算出を行うS22。次に、静電作用スコア511を、算出した実験スコアを用いて補正して補正後静電作用スコアを算出するS23。補正は、静電作用スコア511が実験スコアより大きい場合は実験スコアの値に置換え、それ以外の場合はそのままの値とした。式3なお、静電作用スコアの補正では、補正前の静電作用スコアと実験スコアとの差分に応じて静電作用スコアの値を変更してもよく、複数の抗原抗体複合体構造から求められた実験スコアの分布を反映した値で補正してもよい。次に、補正後静電作用スコアと形状相補作用スコア512から、以下の式4に基づいて複合体評価スコアの算出を行うS24。C41、C42は、複合体構造の適正の高さに対しての静電作用スコア511と形状相補作用スコア512のそれぞれの貢献度に応じた係数である。これらの係数は以下のように算出できる。PDBから抗原抗体複合体構造を複数抽出し、それらの複合体を構成する抗原と抗体を入力情報としてドッキング予測プログラムsurFitにより複合体候補構造の予測を行い、複合体候補構造とその構造に対応する静電ポテンシャルスコア、形状相補ポテンシャルスコア、結合部位面積を得て、静電ポテンシャルスコアと結合部位面積とを積算して静電作用スコアを算出し、形状相補ポテンシャルスコアと結合部位面積とを積算して形状相補作用スコアを算出する。また、この複合体候補構造の適正の高さとして非特許文献4のfnatに示される値である結合残基ペア予測度を算出し、算出した予測精度を目的変数、静電作用スコアと形状相補作用スコアの抗原抗体毎の標準化値を説明変数として回帰分析を行う。なお、抗原抗体毎に複数の複合体候補構造が得られるが、抗原抗体毎に静電作用スコア、形状相補作用スコアの値の範囲にばらつきがあるため、静電作用スコアと形状相補作用スコア、の値は標準化して回帰分析を行い、その結果はC41=9、C42=14、であった。静電作用スコア511、形状相補作用スコア512共に、高いほうが複合体候補構造の適正が高くなる傾向にある。なお、各係数は異なる値でもよい。次に、本実施の形態に示した複合体候補構造評価装置201による抗原抗体複合体候補構造の適正評価の正当性について考察する。非特許文献5に結合速度定数kassocと解離速度定数kdissocの実験値データが示され、かつPDBに複合体結晶構造が登録されている、PDB-IDPDB登録番号が1C08、1MLC、1NDG、1NDMの抗原抗体複合体について、取得した抗原抗体複合体の抗原と抗体の立体構造情報を入力情報としてsurFitを実行し、抗原抗体複合体候補構造予測と各候補構造の静電ポテンシャルスコア、形状相補ポテンシャルスコアと結合部位面積を取得した。取得した各候補構造について、本実施の形態による抗原抗体複合体候補構造評価方法を用いた複合体評価スコアの取得を行い、さらに非特許文献4のfnatに示される値である結合残基ペア予測度を算出して適正が最も高い複合体候補構造を抽出し、この適正が最も高い複合体候補構造が複合体評価スコアの高い順における順位を算出する。また、本実施の形態ではS23において静電作用スコアの補正を行ったが、この補正を行わずに結合速度定数と解離速度定数を用いることなく、静電作用スコア511と形状相補作用スコア512から式4を用いて複合体評価スコアの算出を行い、適正が最も高い複合体候補構造が複合体評価スコアの高い順における順位を算出する。本実施の形態による結合速度定数と解離速度定数を用いた補正を行った補正ありの場合と、結合速度定数と解離速度定数による補正を行わなかった補正なしの場合の、PDB-ID毎の複合体評価スコアの順位を表1に示す。結合速度定数、解離速度定数による補正を行った場合は、補正なしの場合に比べると、特に補正なしで順位の低かった1MLC、1NDMの抗原抗体複合体構造において順位が高くなっており、補正を行うことにより、より精度よく複合体構造の適正が評価可能となっていることが示されている。かかる構成によれば、抗原と抗体の反応速度定数である抗原抗体反応速度情報501を用いて、抗原抗体複合体候補構造情報502から得られるスコアを補正または再評価することにより、抗原抗体の結合反応状態を正しく示している実験値を用いることで、精度の高い抗原抗体複合体の複合体候補構造の適正評価方法を提供することができる。以上、本発明の実施の形態に係る抗原抗体複合体候補構造評価装置201について説明したが、本発明は、この実施の形態に限定されるものではない。なお、本実施の形態において、反応速度情報取得部211において結合速度定数と解離速度定数を取得し、統合スコア算出部221において結合速度定数と解離速度定数の両方を用いて複合体評価スコアの算出を行ったが、結合速度定数と解離速度定数のどちらか一方のみを用いて複合体評価スコアの算出を行ってもよい。また、本実施の形態において、構造的相互作用スコア取得部212は、抗原抗体複合体候補構造情報502から静電ポテンシャルスコアや形状相補ポテンシャルスコア、結合部位面積を取得したが、複合体の立体構造情報から静電ポテンシャルスコアや形状相補ポテンシャルスコア、結合部位面積を算出する複合体候補構造スコア評価部を備えてもよい。また、構造的相互作用スコア取得部212は、静電ポテンシャルスコアと形状相補ポテンシャルスコアと結合部位面積から静電作用スコアと形状相補作用スコアを算出したが、抗原抗体複合体候補構造情報502からこれらの情報を直接取得してもよい。また、本実施の形態において、構造的相互作用スコア取得部212は静電作用スコアと形状相補作用スコアを取得し、総合スコア算出部221は静電作用スコアと形状相補作用スコアを用いて複合体評価スコアを算出したが、構造的相互作用スコア取得部212が取得する値は、静電相互作用スコア、形状相補作用スコア以外にも、疎水的な相互作用状態を表す疎水作用スコアや結合部位面積などの複合体立体構造情報から算出される抗原と抗体間の相互作用状態を示す値を含むことができる。また、静電相互作用スコア、形状相補作用スコア、疎水作用スコア、結合部位面積などの複合体立体構造情報から算出される抗原と抗体間の相互作用状態を示すいずれかの値を含めばよい。