質量分析を用いたペプチド同定方法
【課題】ペプチドの開裂に伴い生じる複数種のフラグメントイオン強度の予測値を高い精度で且つ効率良く求める。
【解決手段】フラグメントイオンを構成するアミノ酸残基のそれぞれについて原子数にスケーリング因子を乗じた有効原子数を考え、各アミノ酸残基の有効原子数の和をフラグメントイオンの有効原子数とし（Ｓ１）、この有効原子数を所定の対数式に代入してプロトン親和力を求める（Ｓ２）。アミノ酸配列の特異性などに応じてスケーリング因子を変えることで、これにより、たとえアミノ酸残基の構成が同じでもプロトン親和力が変化する場合がある。同じフラグメント化反応で生じる複数のフラグメントイオンのプロトン親和力に応じて、活性化エネルギーから計算されるフラグメントイオン総量を分配し、より精度の高いフラグメントイオン強度を求める。
【解決手段】フラグメントイオンを構成するアミノ酸残基のそれぞれについて原子数にスケーリング因子を乗じた有効原子数を考え、各アミノ酸残基の有効原子数の和をフラグメントイオンの有効原子数とし（Ｓ１）、この有効原子数を所定の対数式に代入してプロトン親和力を求める（Ｓ２）。アミノ酸配列の特異性などに応じてスケーリング因子を変えることで、これにより、たとえアミノ酸残基の構成が同じでもプロトン親和力が変化する場合がある。同じフラグメント化反応で生じる複数のフラグメントイオンのプロトン親和力に応じて、活性化エネルギーから計算されるフラグメントイオン総量を分配し、より精度の高いフラグメントイオン強度を求める。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ペプチド混合物を含む試料を質量分析し、これにより得られた質量スペクトルデータを用いて各ペプチドのアミノ酸配列を同定するための方法に関する。
【背景技術】
【０００２】
近年、ポストゲノム研究としてタンパク質の構造や機能の解析が急速に進められている。このようなタンパク質の構造・機能解析手法（プロテオーム解析）の一つとして、質量分析装置を用いたタンパク質の発現解析や一次構造解析が広く行われるようになってきており、四重極型イオントラップや衝突誘起分解（ＣＩＤ）などによって特定のピークの捕捉と開裂を行う、いわゆるＭＳｎ分析（ｎは２以上の整数）が威力を発揮している。一般にＭＳ２（＝ＭＳ／ＭＳ）分析では、まず、分析対象物から特定の質量電荷比（ｍ／ｚ）を有するイオンをプリカーサイオンとして選別し、該プリカーサイオンをＣＩＤによって開裂させる。その後、開裂によって生成したイオン（プロダクトイオン）を質量分析することによって、目的とするイオンの質量や化学構造についての情報を得ることができる。
【０００３】
上記のようなＭＳｎ分析によってタンパク質のアミノ酸配列を決定する場合には、まず、タンパク質を適当な酵素で消化してペプチド断片の混合物としてから、該ペプチド混合物を質量分析する。このとき、各ペプチドを構成する元素には質量の異なる安定同位体が存在するため、同一のアミノ酸配列から成るペプチドであっても、その同位体組成の違いによって質量電荷比の異なる複数のピークを生じる。該複数のピークは、天然存在比が最大の同位体のみで構成されたイオン（主イオン）のピークと、それ以外の同位体を含むイオン（同位体イオン）のピークから成り、これらは１Ｄａから数Ｄａ間隔で並んだ複数本のピークから成るピーク群を形成する。
【０００４】
続いて、上記のようなペプチド混合物の質量スペクトルデータの中から、単一のペプチドに由来する一組の同位体ピーク群をプリカーサイオンとして選択し、該プリカーサイオンを開裂させて得られたイオン（プロダクトイオン）の質量分析（ＭＳ２分析）を行う。
【０００５】
以上のようにして得られたプロダクトイオンのスペクトルパターンや、上記プリカーサイオンのスペクトルパターンをデータベース検索に供することによって、被検ペプチドのアミノ酸配列を決定することができる。このように質量分析により得られた質量スペクトルデータに基づいてペプチドを同定するための標準的なソフトウエアが従来より利用されている。例えば代表的なプロテオーム解析用データベース検索エンジンとして、マトリックスサイエンス社が提供しているマスコット（MASCOT）がよく知られている。
【０００６】
MASCOTを用いた解析では、まず得られた質量スペクトルに出現しているフラグメントイオンの質量をキーとして、データベースに登録されているペプチドの中から質量がマッチするペプチドを検索し、検索により得られた複数の候補ペプチド名をマッチ度順にランク付けしてリストアップする。分析者はこのリストを参考にしてペプチドを推定する。しかしながら、実際には候補ペプチドのリスト上で最高ランクのものが正解の物質でないことがよくあり、上記ソフトウエアによる検索は必ずしも信頼性が高いものではない。その理由として、上記ソフトウエアではフラグメントイオンの質量のマッチングは考慮されているものの、フラグメントイオンのピーク強度とのマッチングは全く又は殆ど考慮されていないことに拠るものと考えられる。
【０００７】
また、MASCOTとは別にサーモ・エレクトロン社が提供するシークエスト（SEQUEST）というデータベース検索エンジンも利用されている。このソフトウエアでは、フラグメントイオン質量とそのピーク強度との両方のマッチングが考慮されている。しかしながら、基本的にはMASCOTと同様に、データベースとしては質量に関するものしか備えておらず、強度に関しては非常に簡易的な強度予測（具体的にはb/yイオン強度は全て50、aイオンや脱水・脱塩イオン強度は10）の下で候補ペプチドのフラグメントイオン強度を与えて測定により得られた質量スペクトルデータとのマッチングを評価しているだけである。したがって、強度マッチングに関してはあまり信頼性を期待できない。
【０００８】
前述のような従来の問題に鑑み、本発明者は、候補ペプチドとして与えられたアミノ酸配列に対して質量スペクトルパターンを予測し、実測のフラグメントイオン強度との相関を計算することで、候補ペプチドの信頼性を評価してペプチドを同定するという手法を提案している（非特許文献１など参照）。こうした手法では、質量スペクトルパターンを予測するためにペプチドのフラグメント強度を高い精度で予測することが重要である。従来のペプチドのフラグメント強度の基本的な予測方法を説明する。
【０００９】
ペプチド分子イオンの開裂の過程は単分子反応と呼ばれる化学反応として知られており、その反応速度は、或る温度での化学反応の予測式であるアレニウス（Arrhenius）の式を用いて与えられる。
フラグメントイオン強度∝反応定数＝Ａ・exp（−Ｅａ／Ｒ・Ｔ） …(1)
ここでＥａは見かけ上の活性化エネルギー、Ｒは気体定数、Ｔは絶対温度、Ａは頻度因子と呼ばれる値である。或いは、活性化エネルギーＥａの代わりにギブスの自由エネルギーΔＧａを用い、
フラグメントイオン強度∝Ａ０・exp（−ΔＧａ／Ｒ・Ｔ） …(2)
として表すこともできる。ここではＡ０が頻度因子である。(2)式は(1)式の変形であって、この変形によって形式上ではあるが頻度因子が異なってくる。
【００１０】
このようなフラグメント化反応は一種類ではなく、ペプチド分子イオンが幾つかのフラグメントイオンに開裂する際には、それぞれの開裂に対応したフラグメント化反応が存在する。その中には、ｂフラグメントイオンやｙフラグメントイオンのように、同じ結合の開裂によるフラグメント化反応で生成されるフラグメントイオンもある。
【００１１】
多くのフラグメント化反応では、アミノ酸配列の長さこそ違え、前述したようにｙタイプとｂタイプのフラグメントイオンが同時に発生することが多い。そして、このフラグメントイオンの総量から個々のフラグメントイオンの強度、即ち、上記例におけるｙフラグメントイオンとｂフラグメントイオンの強度に分配するために、カイネティック法（Kinetic Method)と呼ばれる分子運動論的な方法で分配する方法が知られている（非特許文献２参照）。
【００１２】
一方、最近、数十個のパラメータから成る数学的な算出モデルを仮定し、この仮定の下でフラグメントイオン強度のスペクトルパターンを予測する手段を与えることで、フラグメントイオンの質量と強度とを合わせたマッチングを行う手法が、非特許文献３において提案されている。この文献に記載されている計算モデルにおいては、上記アレニウスの式、及びカイネティック法に形式上合ったモデル式を採用しており、基本的に以下のような手法でペプチドの質量スペクトル強度を求めることができる。
【００１３】
基本式は次の(3)式である。
〔フラグメントイオン強度〕＝〔関係するフラグメント化反応で生じるフラグメントイオンの総量〕×〔フラグメントイオンへの気相塩基性度に応じた配分〕 …(3)
ここで、気相塩基性度（Gas-phase basicity、以下ＧＢと略す）は気相での塩基の強さを表す指標である。上記(3)式の右辺の２つの項は次のような手順で求める。
【００１４】
〔第１ステップ：関係する結合位置でのフラグメント化反応で生じるフラグメントイオンの総量の計算〕
関係するフラグメント化反応で生じるフラグメントイオンの総量は、着目するフラグメント化反応が起こる前提となる分子配置（以下、反応配置という）となる確率に、反応位置から活性化状態へ励起される割合を乗じることにより求める。反応位置となる確率は、フラグメント化反応で開裂するペプチド結合に対する気相塩基性度から求める。一方、反応位置から活性化状態へ励起される割合は、活性化エネルギーをアレニウスの式に適用することで求める。そして、得られた割合の値を上記反応位置となる確率に乗じることで、着目する結合位置で開裂したフラグメントイオン総量を求める。
【００１５】
〔第２ステップ：フラグメントイオンへの気相塩基性度に応じた配分の計算〕
着目しているフラグメント化反応で発生した複数のフラグメントイオンの発生量を個別に求めるために、上記第１ステップで求めたフラグメントイオン総量をカイネティック法に基づく手法、即ち、フラグメントイオンの気相塩基性度の大きさに応じて、exp｛ＧＢ／（Ｒ・Ｔeff）｝（但し、Ｔeffは有効温度と呼ばれ、衝突ガスの種類や圧力、ペプチド分子の大きさ、分解前のイオンの電荷数等に依存する値）の重みを付与して各フラグメントイオンに配分する。この際、フラグメントイオンの気相塩基性度はアミノ酸残基の気相塩基性度を用い、次の(4)式を仮定して求められるものとしている。
exp｛ＧＢ／（Ｒ・Ｔeff）｝＝Σexp｛ＧＢｉ／（Ｒ・Ｔeff）｝ …(4)
ここでＧＢは求めようとするフラグメントイオンの気相塩基性度で、ｉはいま想定しているペプチド分子を構成するアミノ酸残基に付与した連番である。つまり、ＧＢｉはｉ番目のアミノ酸残基の気相塩基性度を意味する。以下の説明においても、便宜上、連番はＮ末端側のアミノ酸残基から順番に付与されているものとする。
【００１６】
以上のように非特許文献３に記載の手法では、気相塩基性度や活性化エネルギー、有効温度などの物理的な量を用いてフラグメントイオン強度を計算している。そして、簡単な計算式や評価式で上記物理量を算出できるようにし、さらに、仮定した計算式にこれら物理量を代入することによりペプチド分子の質量スペクトル強度を求めることができるようにしている。この際、これら物理量を求めるための簡単な計算式や評価式に用いるパラメータは予め、幾つかのフラグメントイオン強度の観測事例に矛盾しないようにパラメータフィッティングなどの手法で決められる。また、仮定した計算式の下でできるだけ多くの事例に対応できるように、一つの物理量を二つ以上のパラメータに分割し、パラメータフィッティングによってそれらパラメータ値を決めるようにしている。
【００１７】
しかしながら、上述した従来のフラグメントイオン強度の算出方法では、計算式や評価式を立てる上での幾つかの仮定や簡単化に無理があり、それ故に正確性を損ねるおそれがあると思われる。具体的には次のようなことが挙げられる。
〔１〕単一アミノ酸の気相塩基性度の評価値は標準的な文献データと矛盾する。例えば、非特許文献３に基づけば、グリシン、アラニンやイソロイシンの気相塩基性度の評価値はお互いの差が１kJ/mol程度でしかないが、ＮＩＳＴで公開されている標準データ（非特許文献４参照）で比べるとその差がグリシンとアラニンとで１５kJ/mol程度、グリシンとイソロイシンとで３０kJ/mol程度、と大きく相違する。
〔２〕フラグメントイオンの気相塩基性度のフラグメントイオンの種類への依存性について例を挙げて言えば、同じアミノ酸配列であれば、ｂタイプのフラグメントイオンとｙフラグメントイオンの気相塩基性度は末端のＮＨ2とＣＯＯＨの寄与について若干考慮されているのみで、ｂタイプのフラグメントイオンの持つ特徴的なＣ末端構造による寄与を十分には反映できていない。
【００１８】
〔３〕(4)式はフラグメントイオンを構成するアミノ酸残基の種類の依存性のみを考慮し、アミノ酸残基の配列の依存性についてはあまり考慮していない。即ち、Ｃ末端側に隣り合うアミノ酸残基の影響を僅かに考慮しているだけである。しかしながら、これは明らかに一般性に欠けるものである。例えばプロリン（ピロリジン−２−カルボン酸）のように、他のアミノ酸と大きく異なり環の中にアミノ基の窒素を含む形態の分子構造を持つような特徴的な分子結合を示すアミノ酸残基が含まれている場合には、そのアミノ酸残基（例えばプロリン）がアミノ酸配列の中でどこに位置するのかによって、イオン化のエネルギー状態は異なる筈である。したがって、少なくとも上記のような特異なアミノ酸残基を含む場合には、アミノ酸残基の配列によって気相塩基性度も変わると考えるのが妥当である。
【００１９】
表１は、プロリン（Ｐ）を含む各種ジペプチドについて、プロリンがＮ末端側に存在する場合（Ｐ−Ｘのアミノ酸配列）とＣ末端側に存在する場合（Ｘ−Ｐのアミノ酸配列）のプロトン親和力（プロトンを引きつける力）ＰＡをそれぞれ分子動力学法による計算で求めた結果を示している。
【表１】
この表により、同じアミノ酸の組み合わせであってもプロリンの位置によってプロトン親和力が相違することは明らかであり、同様のことは気相塩基性度についても予想される。なお、この表において、Ｐ−ＸとＸ−ＰとのＰＡの値の開きは大きな場合で１０〜１５kcal/mol即ち４０〜６０kJ/molほどあり、このような大きな相違は非特許文献３による評価方法では説明できない。また非特許文献３に基づいて、Ｇ−Ｐ、Ａ−Ｐ、Ｉ−Ｐ、Ｌ−ＰやＶ−Ｐなどの気相塩基性度を評価しても、互いの評価値には殆ど差がなく１kJ/mol程度であるが、表１に示すＰＡの計算結果ではその差が最大１３kcal/mol即ち５０kJ/molほどにも及ぶ。
【００２０】
【非特許文献１】梶原、石田、竹内、「第一原理計算を利用したペプチドフラグメント予測法と定性分析への応用」、第５４回質量分析総合討論会、大阪、２００６
【非特許文献２】ラスキン(Julia Laskin)ほか１名、「ザ・セオリティカル・ベイシス・オブ・ザ・カイネマティック・メソッド・フロム・ザ・ポイント・オブ・ビュー・オブ・フィニット・ヒート・バス・セオリー(The Theoretical Basis of the Kinematic Method from the Point of the View of Finite heat bath Theory)」、ジャーナル・オブ・フィジックス・ケミストリー・エー(J. Phys. Chem. A)、2000、pp.8829-8837
【非特許文献３】ジャン(Z.Zhang)、「プレディケイション・オブ・ロー−エナジー・コリジョン−インデュースド・ディソシエイション・スペクトラ・オブ・ペプチドズ(Predication of Low-Energy Collision-Induced Dissociation Spectra of Peptides)」、アナリティカル・ケミストリー(Analytical Chemistry)、 Vol.76、No.14、July 15, 2004、 pp.3908-3922
【非特許文献４】「ナショナル・インスティテュート・オブ・スタンダーズ・アンド・テクノロジー（National Institute of Standards and Technology）」、ＮＩＳＴ、[平成１８年１１月１５日検索]、インターネット＜URL : http://webbook.nist.gov/＞
【発明の開示】
【発明が解決しようとする課題】
【００２１】
本発明は上記課題を解決するために成されたものであり、その目的とするところは、与えられたペプチドのフラグメントイオン強度を従来よりも正確に予測することにより、質量スペクトルデータに基づくタンパク質やペプチドのアミノ酸配列の決定精度を一層高めることができるペプチド同定方法を提供することである。
【課題を解決するための手段】
【００２２】
上記課題を解決するために成された本発明は、ＭＳｎ分析により取得された質量スペクトルデータに基づいて被検試料中のペプチド混合物のアミノ酸配列を決定するペプチド同定方法であって、候補ペプチドのフラグメントイオン強度を予測して質量スペクトルパターンを求め、これと実測の質量スペクトルとの比較により候補ペプチドの妥当性を判断するペプチド同定方法において、ペプチドのフラグメント化反応で生じる複数種のフラグメントイオン毎のフラグメントイオン強度を予測するために、
a)フラグメント化反応で生じるフラグメントイオンの強度を求める強度計算ステップと、
b)該フラグメントイオン強度をフラグメントイオン種類毎に与えられるプロトン親和力又は気相塩基性度に応じて分配する分配ステップと、
とを含み、フラグメントイオン種類毎のプロトン親和力又は気相塩基性度を計算する際に、フラグメントイオンを構成するアミノ酸残基のそれぞれについて実際の原子数にアミノ酸配列の特異性を反映した補正係数を乗じた有効原子数を求め、アミノ酸残基毎の有効原子数の和から求まる該フラグメントイオンの有効原子数を、所定の近似計算式に適用してプロトン親和力又は気相塩基性度を算出するようにしたことを特徴としている。
【００２３】
本発明に係るペプチド同定方法の一態様として、前記所定の近似計算式は、有効原子数とプロトン親和力又は気相塩基性度との非線形関係を対数で近似した式とすることができる。
【発明の効果】
【００２４】
或るアミノ酸配列のペプチド分子についてのフラグメント化反応で生じるフラグメントイオン種類（例えばｙタイプとｂタイプのフラグメントイオン）毎のフラグメントイオン強度を算出する際に、各フラグメントイオンのプロトン親和力（又は気相塩基性度）が必要となるが、従来、プロトン親和力や気相塩基性度といった値はかなり大まかにしか与えられなかった。これに対し本発明に係るペプチド同定方法では、プロトン親和力や気相塩基性度を従来よりも高い精度で、且つ実用的な計算量や計算時間で求める方法を提供する。
【００２５】
即ち、アミノ酸残基毎に原子数に適当な補正係数を乗じた有効原子数なるものを考え、アミノ酸配列に含まれるアミノ酸残基の種類がたとえ同じであっても、その配列順序によってプロトン親和力や気相塩基性度が影響を受けるような場合には、補正係数をそれぞれに合わせて修正することで有効原子数を変え、それによってプロトン親和力や気相塩基性度も変化するようにした。具体的には、例えば、アミノ酸配列中に塩基性アミノ酸が存在する場合に、該塩基性アミノ酸がＮ末端に存在するか否かで異なる補正係数を与えるようにするとよい。また、アミノ酸配列中にプロリンが存在する場合に、該プロリンと隣接する他のアミノ酸残基との組み合わせで該プロリンがＮ末端に存在するか否かで異なる補正係数を与えるようにするとよい。
【００２６】
また同一のフラグメント化反応から生じるフラグメントイオンでもそのタイプによって有効原子数を変えることで、プロトン親和力や気相塩基性度の精度をさらに上げることができる。例えばｂ−ｙフラグメント化反応ではｙフラグメントイオンとｂフラグメントイオンとが生成されるが、同じアミノ酸配列であっても後者は前者より原子数が３個少ないためにそれを有効原子数に反映させる。
【００２７】
従来、(4)式にあるように、気相塩基性度（プロトン親和力も同様）はフラグメントイオンの種類やアミノ酸残基の配列の影響を考慮していないものであったが、本発明に係るペプチド同定方法によれば、プロトン親和力や気相塩基性度をフラグメントイオンの種類やアミノ酸残基の配列の影響を考慮してより高い精度で求めることができる。それにより、ペプチドのフラグメント化反応で生じるフラグメントイオン種類毎にフラグメントイオン強度を高い精度で計算することができ、それ故に、これに基づいたペプチドの同定の信頼性を向上させることができる。
【００２８】
また、上記のような補正係数は特にプロトン親和力や気相塩基性度に対する影響が大きいものについてのみ異なる値を用意しておけばよく、さらに同定処理の際に一々計算するのではなく、予め計算しておいた値をデータベース化して同定処理の際にはデータベース検索により適当な値を得るようにすることができる。したがって、上記のような精度向上を図りながら、コンピュータ上での計算量や計算時間が膨大になることを回避でき、効率的な同定処理を実行してスループットの向上を図ることができる。
【発明を実施するための最良の形態】
【００２９】
本発明の一実施形態によるペプチド同定方法の全体の手順をまず説明する。図１は本実施形態のペプチド同定方法における解析処理の概略フローチャートである。実際には、この解析処理は、新規のプロテオーム解析用ソフトウエアをコンピュータ上で動作させることにより実行される。
【００３０】
まず分析対象の物質（ペプチド）を例えばイオントラップ飛行時間型質量分析装置で質量分析（ＭＳｎ分析）して質量スペクトルデータを収集する（ステップＳ１）。次いで、この質量スペクトルデータに対し、既存のプロテオーム解析用データベース検索ソフトウエア（例えばMASCOT）を利用して解析処理を実行する。即ち、質量スペクトルに出現しているフラグメントイオンの質量（厳密には質量電荷比、ｍ／ｚ値）をキーとして、データベースに登録されているペプチドの中から質量がマッチするペプチドを検索する（ステップＳ２）。そして、検索により得られた複数のペプチド分子をマッチ度順にランク付けしてリストアップする（ステップＳ３）。従来は、このリストアップされたペプチド分子の中から、分析者自らの判断（例えば最もスコアが高いものを選択する等）により同定を行う必要があったが、本発明に係る解析処理によれば、ステップＳ４以降の処理によりさらに候補の絞り込み（換言すればより確度の高い判断材料の提供）を行うことができる。
【００３１】
即ち、リストアップされた候補ペプチドの中で他と比べて明らかにスコアが高い複数の候補ペプチドを選定し、各候補のそれぞれのアミノ酸配列についてフラグメントイオン強度を計算することで質量スペクトルパターンを予測する。そのために、まずフラグメント化反応で生じるフラグメントイオン総量を計算し（ステップＳ４）、さらに計算により求まるプロトン親和力（又は気相塩基性度）に応じて同一のフラグメント化反応で生じる複数のフラグメントイオンへの強度の配分を行ってフラグメントイオン強度を求める（ステップＳ５）。こうして求められたフラグメントイオン強度から各候補ペプチドの質量スペクトルパターンを予測し、実測により得られた質量スペクトルとの相関による強度マッチング評価を行う（ステップＳ６）。そして、その相関性の評価結果を比べて最も可能性の高い候補ペプチドを選定して例えばモニタの表示画面上等に提示する（ステップＳ７）。このステップＳ４〜Ｓ７による新規の機能が、従来の標準的なデータベース検索ソフトウエアによる機能と組み合わせられ、これにより従来よりも信頼度の高い同定結果が期待できる。
【００３２】
本実施形態によるペプチド同定方法は、特にステップＳ４及びＳ５の処理に特徴を有している。これらについて以下に詳しく説明する。
【００３３】
（１）フラグメント化反応で生じるフラグメントイオン総量の計算方法
図３はフラグメントイオン強度の計算方法をエネルギー準位で示した原理説明図であり、（ａ）はプロトン親和力をベースにした、（ｂ）は気相塩基性度をベースにした原理説明図である。表現は異なるが、いずれも(1)式で示したアレニウスの反応式に基づいたフラグメントイオン強度計算式を利用している。
【００３４】
ＣＩＤ等のフラグメント化反応で生じるフラグメントイオンの総量は、(1)式に示したアレニウスの反応式を用いて求める。その際、パラメータの１つとして関連フラグメントの活性化エネルギーＥａを与える必要がある、活性化エネルギーは、ペプチドの初期状態から、着目しているフラグメントイオンの活性化状態、即ち遷移状態までのエネルギーの差として直接計算で求めることも可能であるが、ここでは二つの状態エネルギー差の和で考える。即ち、例えば図３（ａ）に示すように、初期状態であるプリカーサイオン（ＸＹ−Ｈ＋）の状態から開裂するペプチド結合位置（反応分子位置と呼ぶ）までプロトン（Ｈ＋）が移動するための移動エネルギーＥＨと、反応分子位置（Ｘ−Ｈ＋−Ｙ）から遷移状態（Ｘ−Ｈ＋−Ｙ＊）に移行するためのエネルギー変化量（以下、活性化エネルギーＥａ’と呼ぶ）とに分けて計算を行う。
【００３５】
このようにエネルギーを分割した上で、近似計算として、活性化エネルギーＥａ’は開裂するペプチド結合の両側に存在する２又は３個のアミノ酸配列分子モデルに簡単化して求める。このように簡単化できる理由は、フラグメント化反応がペプチド結合付近の局所的な反応であるとみて構わないためである。また、プロトンの移動エネルギーＥＨに関しては、移動の初めと最後にプロトンの存在するアミノ酸を組み合わせてできるジペプチドに分子モデルを簡単化して求める。このように分子モデルを簡単化できるのは、プロトンはそれ自身近くの原子と強く結合していて、遠くに存在する原子や分子との結合力は無視できるものとし、プロトンとの結合に直接関与しない部分を省略しても構わないためである。
【００３６】
上述したように２又は３個のアミノ酸配列分子モデルに簡単化して活性化エネルギーＥａ’やプロトンの移動エネルギーＥＨを求めるためには、予めアミノ酸配列の組み合わせ毎に活性化エネルギーＥａ’、プロトン移動エネルギーＥＨを計算してデータベース化しておき、結合を与えるとそのデータベースから対応したエネルギーＥａ’、ＥＨが導出されるようにしておくとよい。共通アミノ酸の種類はたかだか２０種類であるが、例えば５個のアミノ酸配列の組み合わせを考えただけでも、ｂフラグメントイオン、ｙフラグメントイオンの２種を考えると、エネルギーの計算対象は２０5×２＝６４×１０5と膨大な数となり、これを全て計算するのには膨大な時間を要する。これに対し、上記のように２又は３個のアミノ酸配列分子モデルに簡単化することで、計算対象の数を大幅に減らすことができ、コンピュータを用いて比較的短時間で計算することができるので、活性化エネルギーＥａ’、プロトン移動エネルギーＥＨを呼び出すためのデータベースを比較的容易に構築することができる。
【００３７】
なお、過去の研究により、フラグメント化反応は大別して、モバイルプロトンモデル又はリモートプロトンモデルと呼ばれる２つの反応メカニズムのいずれかで行われることが知られており、リモートプロトンモデルではフラグメント化反応においてプロトンが移動しない場合がある（それらの反応メカニズムの詳細については例えば特願２００６−４５５３号など参照）。この場合にはプロトンの移動エネルギーＥＨの計算は行わずに直接、活性化エネルギーＥａのみを計算するが、その場合でも計算方法はＥａ’の計算方法に準じればよい。
【００３８】
また、フラグメント化反応の過程をプロトン親和力でなく気相塩基性度で考える場合には、図３（ａ）と対比した図３（ｂ）に示すように、プロトンの移動エネルギーＥＨに代えたプロトンの移動自由エネルギー（ΔＧＨ）、Ｅａ’に代えた反応分子位置（Ｘ−Ｈ＋−Ｙ）から遷移状態（Ｘ−Ｈ＋−Ｙ＊）を経て分解するときの活性化自由エネルギー（ΔＧａ’）の和で活性化自由エネルギー（ΔＧａ）を表せばよい。
【００３９】
以上のようにして活性化エネルギーＥａ（又は活性化自由エネルギーΔＧａ）が求まれば、他のパラメータ、気体定数Ｒ、温度Ｔ、頻度因子Ａ（又はＡ０）は既知であるから、アレニウスの反応式によりフラグメントイオン強度、つまりはフラグメント化反応で同時に生じるＸイオンとＹイオンのフラグメントイオンの総量を算出することができる。
【００４０】
（２）複数種のフラグメントイオンへのプロトン親和力（又は気相塩基性度）に応じたフラグメントイオン強度の配分計算方法
質量スペクトルパターンを求めるには、着目しているフラグメント化反応で発生した複数種類のフラグメントイオンの発生量を種類毎（例えばｂフラグメントイオンとｙフラグメントイオン）に個別に求める必要がある。そこで、ステップＳ４で求めたフラグメントイオン総量をカイネティック法に基づく手法、即ち、フラグメントイオンのプロトン親和力ＰＡ（又は気相塩基性度ＧＢ）の大きさに応じて、exp｛ＰＡ／（Ｒ・Ｔeff）｝（又はexp｛ＧＢ／（Ｒ・Ｔeff）｝）の重みでもって各フラグメントイオンに配分する。
【００４１】
例えば代表的なフラグメント化反応プロセスであるｂ−ｙフラグメント化反応では、ｂフラグメントイオンとｙフラグメントイオンとが発生するが、ｂフラグメントイオン、ｙフラグメントイオンのプロトン親和力（又は気相塩基性度）をそれぞれＰＡｂ、ＰＡｙ（又はＧＢｂ、ＧＢｙ）としたとき、配分比は、
ｂフラグメントイオン配分比＝expＰＡｂ／（expＰＡｂ＋expＰＡｙ） …(5)
ｙフラグメントイオン配分比＝expＰＡｙ／（expＰＡｂ＋expＰＡｙ） …(6)
となる。気相塩基性度ＧＢを用いる場合には、上記(5)、(6)式においてＰＡｂ、ＰＡｙの代わりにＧＢｂ、ＧＢｙを使用すればよい。この配分方法自体はカイネテッィク法そのものであるが、この計算に用いるフラグメントイオンのプロトン親和力（又は気相塩基性度）を次のような特徴的な方法で求める。
【００４２】
過去の研究（ジョン・ホーメス、ほか２名(John L. Hohmes, Christiane Aubry, and Paul M. Mayer)、「プロトン・アフィニティーズ・オブ・プライマリー・アルカノールズ：アン・アプレイザル・オブ・ザ・カイネマティック・メソッド(Proton Affinities of Primary Alkanols: An Appraisal of the Kinematic Method)」、ジャーナル・フィジックス・ケミストリー(J. Phys. Chem.）、A,103, pp.705-709(1999)参照）によれば、同一系列の一級アルカノール(Primary Alkanols) 化合物について、原子数とプロトン親和力との間には一つの曲線関係が存在すると仮定し、この曲線を利用して同一系列に属する未知のアルカノールのプロトン親和力を予測している。
【００４３】
本実施形態によるペプチド同定方法でも基本的にはこの原理を利用してフラグメントイオンの原子数からプロトン親和力を求めるが、単なる原子数ではなくアミノ酸残基の種類や配列などを考慮した有効原子数という新たな概念を導入する。即ち、フラグメントイオンを構成するアミノ酸残基のそれぞれについて、原子数に補正係数（以下、スケーリング因子と呼ぶ）を乗じた有効原子数を与え、各アミノ酸残基の有効原子数の和をフラグメントイオンの有効原子数とする。また、フラグメントイオンのプロトン親和力（又は気相塩基性度）と有効原子数との非線形関係を対数式で表しておき、この対数式に有効原子数を代入することでプロトン親和力（又は気相塩基性度）を簡単に求められるようにする。
【００４４】
具体的には、例えばｙフラグメントイオンの場合について考えると、ｙフラグメントイオンは一般にｎ個のアミノ酸残基から成るペプチド分子イオンの構造を有しているので、次の(7)式により有効原子数Ｎを表す。
Ｎ＝Ｐ１・Ｌ１＋Ｐ２・Ｌ２＋…＋Ｐｎ・Ｌｎ＋３ …(7)
ここでＬｉ及びＰｉはそれぞれ、アミノ酸配列のＮ末端側から数えてｉ番目のアミノ酸残基の原子数及びスケーリング因子である。全てのスケーリング因子Ｐｉが１ならば、有効原子数Ｎはｙフラグメントイオンの原子数と同じとなる。ここで、(7)式の右辺右端の「３」は、Ｃ末端に結合する水酸基（−ＯＨ）とＮ末端に結合する水素基（−Ｈ）の原子数である。そして(7)式で定義した有効原子数Ｎを用いて、ｙフラグメントイオンのプロトン親和力を次の(8)式で表す。
ＰＡ＝ａ・ln｛ｂ・（Ｎ−３）｝＝ａ・ln｛ｂ・（Ｐ１・Ｌ１＋Ｐ２・Ｌ２＋…＋Ｐｎ・Ｌｎ）｝ …(8)
ここでａ、ｂは定数であり、実際のプロトン親和力の値とできるだけマッチングするように予めパラメータフィッティングで決められる。
【００４５】
なお、(8)式に示した対数式において、有効原子数Ｎをそのまま用いずに（Ｎ−３）を変数としているが、これはＮとしても構わない。但し、その場合には定数ａ、ｂを新たなパラメータフィッティングで決めればよい。
【００４６】
一方、ｂフラグメントイオンの場合には、次の(9)式によりプロトン親和力ＰＡｂを計算する。
ＰＡｂ＝ａ・ln｛ｂ・（Ｎ−６）｝＝ａ・ln｛ｂ・（Ｐ１・Ｌ１＋Ｐ２・Ｌ２＋…＋Ｐｎ・Ｌｎ−３）｝ …(9)
同じアミノ酸残基から成るｙフラグメントイオンとｂフラグメントイオンとでは、ｂフラグメントイオンのほうが原子数が３個少ないため、(9)式ではこれを考慮した式としている。
【００４７】
上記(8)、(9)式においては、フラグメントイオンを構成するアミノ酸残基の種類には依存するものの、その序列には依存しない形式としてプロトン親和力を定義している。但し、プロトン親和力の値がアミノ酸残基の序列にも多少の影響を受けることは、従来から知られている。特にアミノ酸配列の両端であるＮ末端、Ｃ末端は分子構造が途切れることからプロトン親和力が影響を受け易いと考えられ、その中でもＮ末端はモバイルプロトンモデルにおいては、フラグメント化反応の初期にプロトンのいる分子位置であることから、このＮ末端にプロリンや塩基性アミノ酸などの特定のアミノ酸が存在するとプロトン親和力が無視できない程度に変化することが予想される。
【００４８】
そこで、特に分子構造の途切れるＮ末端での影響の大きさを考慮し、Ｎ末端に塩基性アミノ酸が存在する場合にはＮ末端にそうしたアミノ酸が無い場合とは異なるスケーリング因子を割り当てるようにする。また、プロリンのような他のアミノ酸とは異なる特徴的な分子構造を持つアミノ酸については、隣接するアミノ酸との組み合わせによってもプロトン親和力が影響を受けると考えられる。そこで、プロリンのような特異なアミノ酸が存在する場合には、そのプロリンと隣接するアミノ酸との組み合わせ一つ一つについて異なるスケーリング因子を与えるものとする。このようにスケーリング因子が変わると、仮に元の原子数は同じであっても有効原子数は変わり、それに伴いプロトン親和力も変化することになる。
【００４９】
気相塩基性度ＧＢに対してもプロトン親和力ＰＡと同様に表すことができる。即ち、ｙフラグメントイオン及びｂフラグメントイオンについてそれぞれ次式で気相塩基性度ＧＢｙ、ＧＢｂを表すものとする。
ＧＢｙ＝ａ’・ln｛ｂ’・（Ｐ１’・Ｌ１＋Ｐ２’・Ｌ２＋…＋Ｐｎ’・Ｌｎ）｝ …(10)
ＧＢｂ＝ａ’・ln｛ｂ’・（Ｐ１’・Ｌ１＋Ｐ２’・Ｌ２＋…＋Ｐｎ’・Ｌｎ−３）｝ …(11)
ここでａ’、ｂ’は定数であり、実際の気相塩基性度とできるだけマッチするようにカーブフィッティングで予め決めておく。Ｐｉ’、Ｌｉはそれぞれｉ番目のアミノ酸残基の原子数、スケーリング因子である。そしてプロトン親和力と同様に、特定のアミノ酸がＮ末端にある場合とそうでない場合とで異なるスケーリング因子を割り当てて気相塩基性度が変化するようにする。またプロリンについても同様である。
【００５０】
具体的には、(8)〜(11)式で用いられるスケーリング因子Ｐｉ（又はＰｉ’）は定数ａ、ｂ（又はａ’、ｂ’）とともに実験データとマッチングするように決められるべき定数であるから、これらの値は多くの観測データと比較して最小二乗法により予め求めておくことができる。そうした求めた値を例えばテーブル形式でデータベース化しておき、与えられたアミノ酸配列に応じて適当な値を読み出して来て上記のような計算式を用いて、フラグメントイオン毎にプロトン親和力（又は気相塩基性度）を算出することができる。
【００５１】
また、別の方法でスケーリング因子Ｐｉを決めることもできる。即ち、ｎ＝１の場合、ｙフラグメントイオンはアミノ酸イオンであることに着目すれば、アミノ酸のプロトン親和力はよく知られているので、逆にアミノ酸のプロトン親和力からスケーリング因子を求めることができる。この場合は、(8)式においてｉ番目のアミノ酸残基一つのみから成るｙフラグメントイオン（＝アミノ酸イオン）を想定し、そのアミノ酸イオンのプロトン親和力をＰＡｉと置けば、(8)式は
ＰＡｉ＝ａ・ln（ｂ・Ｐｉ・Ｌｉ） …(12)
つまり、
Ｐｉ＝｛１／（ｂ・Ｌｉ）｝exp（ＰＡｉ／ａ） …(13)
と変形することができる。したがって、実験や理論的な考察などから得られるアミノ酸イオンのプロトン親和力ＰＡｉの値を(13)式に代入すれば、スケーリング因子Ｐｉを求めることができる。
【００５２】
また(8)式を変形すると、
exp（ＰＡ／ａ）＝ｂ・（Ｐ１・Ｌ１＋Ｐ２・Ｌ２＋…＋Ｐｎ・Ｌｎ） …(14)
のようになる。この式に上記(13)式を代入すれば、
exp（ＰＡ／ａ）＝exp（ＰＡ１／ａ）＋exp（ＰＡ２／ａ）＋…＋exp（ＰＡｎ／ａ） …(15)
となる。この式は形式上は(4)式と似ているが、意味するところは全く相違する。即ち、(15)式の指数因子の分母ａは定数であり、一方、(4)式の指数因子の分母Ｒ・Ｔeffで有効温度Ｔeffはペプチド分子の質量電荷比、ｍ／ｚ値等に依存する値であって定数ではない（上記非特許文献３参照）。
【００５３】
なお、必ずしも(13)式を用いてスケーリング因子Ｐｉを決める必要はなく、できるだけ実際のフラグメントイオンのプロトン親和力（又は気相塩基性度）にマッチするようにパラメータフィッティング等で決めればよい。
【００５４】
上述の方法でフラグメントイオン毎のプロトン親和力を算出する際の具体例を説明する。表２はポリグリシン（Ｇ、ＧＧ、ＧＧＧ、ＧＧＧＧ、ＧＧＧＧＧ）やポリアラニン（Ａ、ＡＡ、ＡＡＡ、ＡＡＡＡ）など、分子量の小さなアミノ酸の配列から構成されるペプチドのプロトン親和力を、Ｂ３ＬＹＰ／６−３１１＋Ｇ（２ｄ，ｐ）のレベルでの密度汎関数法を用いた分子動力学計算で求めた結果を示している。
【表２】
【００５５】
この結果を利用して、(8)式にあるパラメータａ、ｂを、横軸を原子数Ｎ−３、縦軸をプロトン親和力とした曲線のカーブフィッティングから決めることができる。即ち、グリシン（Ｇ）やアラニン（Ａ）などのアミノ酸のスケーリング因子Ｐｉをともに１であると仮定してカーブフィッティングを行うと、その結果は、
ＰＡ[kcal/mol]＝１１．３９５・ln（Ｎ−３）＋１９０．８６ …(16)
と得られる。これを(8)式の形に変形すれば、定数ａ、ｂは、
ａ＝１１．３９５[kcal/mol]、ｂ＝１．８８０１６×１０７
と求まる。
【００５６】
表３はグリシン、アラニン以外の各種アミノ酸に対して上記定数ａ、ｂの値を用い、アミノ酸のプロトン親和力の値からスケーリング因子を求めた例である。ここで示していないアミノ酸についても同様にしてスケーリング因子を求めることができる。
【表３】
このように定数ａ、ｂ及びスケーリング因子Ｐを決めておけば、(8)式を用いて任意のｙフラグメントイオンのプロトン親和力を、(9)式を用いて任意のｂフラグメントイオンのプロトン親和力をそれぞれ予測値として計算することができる。
【００５７】
表４は上記結果を利用して２個のアミノ酸から成るｙフラグメントイオン、換言すればジペプチドイオンに対するプロトン親和力ＰＡの予測値を算出した例である。この例では、アラニン（Ａ）を含むジペプチドのプロトン親和力を予測しているが、Ｂ３ＬＹＰ／６−３１１＋Ｇ（２ｄ，ｐ）のレベルでの密度汎関数法を用いた分子動力学計算から求めた結果と比較的良好に一致しており、上述のような本実施形態におけるプロトン親和力の予測方法が実用的に有効であることを確認することができる。
【表４】
【００５８】
表５は、(9)式を利用して２個のグリシン（Ｇ）から成るｂフラグメントイオンのプロトン親和力を予測し、それを分子動力学による計算値と比較した結果である。この場合にも両者は比較的良好に一致していることが分かる。
【表５】
【００５９】
前述のように表１はプロリン（Ｐ）とその他のアミノ酸の組み合わせによるジペプチドについてプロトン親和力の分子動力学的計算値を示したものであるが、この表１の結果から求められたプロリンのスケーリング因子を表６に示す。
【表６】
プロリン以外のアミノ酸残基のスケーリング因子は既に決定されているものとし、表１に示したジペプチドのプロトン親和力を(8)式に代入することにより、プロリンのスケーリング因子を決定することができる。
【００６０】
Ｎ末端にプロリンが存在する場合のプロリンのスケーリング因子は、表１中のＰ−Ｘ態様のジペプチドのプロトン親和力の値を(8)式に代入して得ることができる。そして、表６中の左側に示したスケーリング因子が、プロリンと結合するアミノ酸残基の種類に応じて与えられる。またプロリンがＮ末端に無い場合のスケーリング因子は、表１中のＸ−Ｐ態様のジペプチドのプロトン親和力の値を(8)式に代入して得ることができ、表６中の右側に示したスケーリング因子が、プロリンと結合するアミノ酸残基の種類に応じて与えられる。
【００６１】
以上例示したように、アミノ酸配列に応じて適宜のスケーリング因子を決定することができる。図２は実際にコンピュータでフラグメントイオン強度の配分計算を行う際の処理手順を示すフローチャートである。即ち、評価対象のペプチド分子が与えられると、そのアミノ酸配列に応じて所定のアルゴリズムに従ってスケーリング因子が決定される（ステップＳ１１）。具体的には、例えばアミノ酸配列が入力されると、例えば上述したようにスケーリング因子に影響を与える条件が抽出され、その条件に応じたデータベース検索によりアミノ酸残基毎のスケーリング因子が求まる。次に、各アミノ酸残基のスケーリング因子と原子数とから、(8)、(9)式に基づいて各フラグメントイオンのプロトン親和力を計算する（ステップＳ１２）。或いは、(10)、(11)式に基づいて各フラグメントイオンの気相塩基性度を計算してもよい。
【００６２】
そうして各フラグメントイオンのプロトン親和力又は気相塩基性度が求まったならば、例えばｂ−ｙフラグメントであれば(5)、(6)式に基づいてフラグメントイオン強度を配分する（ステップＳ１３）。もちろん、ｂ−ｙフラグメントでなくても同様に配分できることは当然である。
【００６３】
以上説明したように本発明に係るペプチド同定方法によれば、よく知られた化学反応の標準的な理論、即ちアレニウスの反応速度式やカイネテッィク法に沿う形で、フラグメントイオン強度を簡単なモデル式により近似計算で求めることができる。こうしたモデル式を用いず直接コンピュータでフラグメントイオン強度を求めようとすると、一般にペプチドのアミノ酸配列のサイズは大きいため、計算量が膨大になって実用的な時間での処理は困難になる。これに対し、上記実施形態で採用しているような方法によれば、計算に必要な活性化エネルギーやプロトン親和力（又は気相塩基性度）をペプチド分子に対して効率的に、つまり実用的な時間の範囲で求めることが可能となる。特にプロトン親和力（又は気相塩基性度）はアミノ酸配列のサイズが大きくなるほど大きくなる傾向にあり、一定の値に収束する傾向はない。そのため、プロトン親和力や気相塩基性度を簡単な評価式から少ない計算量で求めることができることは非常に有効であると言える。
【００６４】
また、各フラグメントイオンのプロトン親和力を計算する際に、そのフラグメントイオンの種類やアミノ酸配列の特異性（例えばプロリンの存在や存在位置など）を考慮して値が修正されるので、より高い精度でもってフラグメントイオン強度を求めることができ、ひいてはペプチドの同定の正確性の向上が期待できる。
【００６５】
なお、上記実施形態は本発明の一例にすぎず、本発明の趣旨の範囲で適宜変形、修正、追加等を行っても本願特許請求の範囲に包含されることは当然である。
【図面の簡単な説明】
【００６６】
【図１】本発明の一実施形態によるペプチド同定方法における解析処理の概略フローチャート。
【図２】フラグメントイオン強度の配分計算を行う際の処理手順を示すフローチャート。
【図３】フラグメントイオン強度の計算方法を反応経路に沿ったエネルギープロファイルで示した原理説明図。
【技術分野】
【０００１】
本発明は、ペプチド混合物を含む試料を質量分析し、これにより得られた質量スペクトルデータを用いて各ペプチドのアミノ酸配列を同定するための方法に関する。
【背景技術】
【０００２】
近年、ポストゲノム研究としてタンパク質の構造や機能の解析が急速に進められている。このようなタンパク質の構造・機能解析手法（プロテオーム解析）の一つとして、質量分析装置を用いたタンパク質の発現解析や一次構造解析が広く行われるようになってきており、四重極型イオントラップや衝突誘起分解（ＣＩＤ）などによって特定のピークの捕捉と開裂を行う、いわゆるＭＳｎ分析（ｎは２以上の整数）が威力を発揮している。一般にＭＳ２（＝ＭＳ／ＭＳ）分析では、まず、分析対象物から特定の質量電荷比（ｍ／ｚ）を有するイオンをプリカーサイオンとして選別し、該プリカーサイオンをＣＩＤによって開裂させる。その後、開裂によって生成したイオン（プロダクトイオン）を質量分析することによって、目的とするイオンの質量や化学構造についての情報を得ることができる。
【０００３】
上記のようなＭＳｎ分析によってタンパク質のアミノ酸配列を決定する場合には、まず、タンパク質を適当な酵素で消化してペプチド断片の混合物としてから、該ペプチド混合物を質量分析する。このとき、各ペプチドを構成する元素には質量の異なる安定同位体が存在するため、同一のアミノ酸配列から成るペプチドであっても、その同位体組成の違いによって質量電荷比の異なる複数のピークを生じる。該複数のピークは、天然存在比が最大の同位体のみで構成されたイオン（主イオン）のピークと、それ以外の同位体を含むイオン（同位体イオン）のピークから成り、これらは１Ｄａから数Ｄａ間隔で並んだ複数本のピークから成るピーク群を形成する。
【０００４】
続いて、上記のようなペプチド混合物の質量スペクトルデータの中から、単一のペプチドに由来する一組の同位体ピーク群をプリカーサイオンとして選択し、該プリカーサイオンを開裂させて得られたイオン（プロダクトイオン）の質量分析（ＭＳ２分析）を行う。
【０００５】
以上のようにして得られたプロダクトイオンのスペクトルパターンや、上記プリカーサイオンのスペクトルパターンをデータベース検索に供することによって、被検ペプチドのアミノ酸配列を決定することができる。このように質量分析により得られた質量スペクトルデータに基づいてペプチドを同定するための標準的なソフトウエアが従来より利用されている。例えば代表的なプロテオーム解析用データベース検索エンジンとして、マトリックスサイエンス社が提供しているマスコット（MASCOT）がよく知られている。
【０００６】
MASCOTを用いた解析では、まず得られた質量スペクトルに出現しているフラグメントイオンの質量をキーとして、データベースに登録されているペプチドの中から質量がマッチするペプチドを検索し、検索により得られた複数の候補ペプチド名をマッチ度順にランク付けしてリストアップする。分析者はこのリストを参考にしてペプチドを推定する。しかしながら、実際には候補ペプチドのリスト上で最高ランクのものが正解の物質でないことがよくあり、上記ソフトウエアによる検索は必ずしも信頼性が高いものではない。その理由として、上記ソフトウエアではフラグメントイオンの質量のマッチングは考慮されているものの、フラグメントイオンのピーク強度とのマッチングは全く又は殆ど考慮されていないことに拠るものと考えられる。
【０００７】
また、MASCOTとは別にサーモ・エレクトロン社が提供するシークエスト（SEQUEST）というデータベース検索エンジンも利用されている。このソフトウエアでは、フラグメントイオン質量とそのピーク強度との両方のマッチングが考慮されている。しかしながら、基本的にはMASCOTと同様に、データベースとしては質量に関するものしか備えておらず、強度に関しては非常に簡易的な強度予測（具体的にはb/yイオン強度は全て50、aイオンや脱水・脱塩イオン強度は10）の下で候補ペプチドのフラグメントイオン強度を与えて測定により得られた質量スペクトルデータとのマッチングを評価しているだけである。したがって、強度マッチングに関してはあまり信頼性を期待できない。
【０００８】
前述のような従来の問題に鑑み、本発明者は、候補ペプチドとして与えられたアミノ酸配列に対して質量スペクトルパターンを予測し、実測のフラグメントイオン強度との相関を計算することで、候補ペプチドの信頼性を評価してペプチドを同定するという手法を提案している（非特許文献１など参照）。こうした手法では、質量スペクトルパターンを予測するためにペプチドのフラグメント強度を高い精度で予測することが重要である。従来のペプチドのフラグメント強度の基本的な予測方法を説明する。
【０００９】
ペプチド分子イオンの開裂の過程は単分子反応と呼ばれる化学反応として知られており、その反応速度は、或る温度での化学反応の予測式であるアレニウス（Arrhenius）の式を用いて与えられる。
フラグメントイオン強度∝反応定数＝Ａ・exp（−Ｅａ／Ｒ・Ｔ） …(1)
ここでＥａは見かけ上の活性化エネルギー、Ｒは気体定数、Ｔは絶対温度、Ａは頻度因子と呼ばれる値である。或いは、活性化エネルギーＥａの代わりにギブスの自由エネルギーΔＧａを用い、
フラグメントイオン強度∝Ａ０・exp（−ΔＧａ／Ｒ・Ｔ） …(2)
として表すこともできる。ここではＡ０が頻度因子である。(2)式は(1)式の変形であって、この変形によって形式上ではあるが頻度因子が異なってくる。
【００１０】
このようなフラグメント化反応は一種類ではなく、ペプチド分子イオンが幾つかのフラグメントイオンに開裂する際には、それぞれの開裂に対応したフラグメント化反応が存在する。その中には、ｂフラグメントイオンやｙフラグメントイオンのように、同じ結合の開裂によるフラグメント化反応で生成されるフラグメントイオンもある。
【００１１】
多くのフラグメント化反応では、アミノ酸配列の長さこそ違え、前述したようにｙタイプとｂタイプのフラグメントイオンが同時に発生することが多い。そして、このフラグメントイオンの総量から個々のフラグメントイオンの強度、即ち、上記例におけるｙフラグメントイオンとｂフラグメントイオンの強度に分配するために、カイネティック法（Kinetic Method)と呼ばれる分子運動論的な方法で分配する方法が知られている（非特許文献２参照）。
【００１２】
一方、最近、数十個のパラメータから成る数学的な算出モデルを仮定し、この仮定の下でフラグメントイオン強度のスペクトルパターンを予測する手段を与えることで、フラグメントイオンの質量と強度とを合わせたマッチングを行う手法が、非特許文献３において提案されている。この文献に記載されている計算モデルにおいては、上記アレニウスの式、及びカイネティック法に形式上合ったモデル式を採用しており、基本的に以下のような手法でペプチドの質量スペクトル強度を求めることができる。
【００１３】
基本式は次の(3)式である。
〔フラグメントイオン強度〕＝〔関係するフラグメント化反応で生じるフラグメントイオンの総量〕×〔フラグメントイオンへの気相塩基性度に応じた配分〕 …(3)
ここで、気相塩基性度（Gas-phase basicity、以下ＧＢと略す）は気相での塩基の強さを表す指標である。上記(3)式の右辺の２つの項は次のような手順で求める。
【００１４】
〔第１ステップ：関係する結合位置でのフラグメント化反応で生じるフラグメントイオンの総量の計算〕
関係するフラグメント化反応で生じるフラグメントイオンの総量は、着目するフラグメント化反応が起こる前提となる分子配置（以下、反応配置という）となる確率に、反応位置から活性化状態へ励起される割合を乗じることにより求める。反応位置となる確率は、フラグメント化反応で開裂するペプチド結合に対する気相塩基性度から求める。一方、反応位置から活性化状態へ励起される割合は、活性化エネルギーをアレニウスの式に適用することで求める。そして、得られた割合の値を上記反応位置となる確率に乗じることで、着目する結合位置で開裂したフラグメントイオン総量を求める。
【００１５】
〔第２ステップ：フラグメントイオンへの気相塩基性度に応じた配分の計算〕
着目しているフラグメント化反応で発生した複数のフラグメントイオンの発生量を個別に求めるために、上記第１ステップで求めたフラグメントイオン総量をカイネティック法に基づく手法、即ち、フラグメントイオンの気相塩基性度の大きさに応じて、exp｛ＧＢ／（Ｒ・Ｔeff）｝（但し、Ｔeffは有効温度と呼ばれ、衝突ガスの種類や圧力、ペプチド分子の大きさ、分解前のイオンの電荷数等に依存する値）の重みを付与して各フラグメントイオンに配分する。この際、フラグメントイオンの気相塩基性度はアミノ酸残基の気相塩基性度を用い、次の(4)式を仮定して求められるものとしている。
exp｛ＧＢ／（Ｒ・Ｔeff）｝＝Σexp｛ＧＢｉ／（Ｒ・Ｔeff）｝ …(4)
ここでＧＢは求めようとするフラグメントイオンの気相塩基性度で、ｉはいま想定しているペプチド分子を構成するアミノ酸残基に付与した連番である。つまり、ＧＢｉはｉ番目のアミノ酸残基の気相塩基性度を意味する。以下の説明においても、便宜上、連番はＮ末端側のアミノ酸残基から順番に付与されているものとする。
【００１６】
以上のように非特許文献３に記載の手法では、気相塩基性度や活性化エネルギー、有効温度などの物理的な量を用いてフラグメントイオン強度を計算している。そして、簡単な計算式や評価式で上記物理量を算出できるようにし、さらに、仮定した計算式にこれら物理量を代入することによりペプチド分子の質量スペクトル強度を求めることができるようにしている。この際、これら物理量を求めるための簡単な計算式や評価式に用いるパラメータは予め、幾つかのフラグメントイオン強度の観測事例に矛盾しないようにパラメータフィッティングなどの手法で決められる。また、仮定した計算式の下でできるだけ多くの事例に対応できるように、一つの物理量を二つ以上のパラメータに分割し、パラメータフィッティングによってそれらパラメータ値を決めるようにしている。
【００１７】
しかしながら、上述した従来のフラグメントイオン強度の算出方法では、計算式や評価式を立てる上での幾つかの仮定や簡単化に無理があり、それ故に正確性を損ねるおそれがあると思われる。具体的には次のようなことが挙げられる。
〔１〕単一アミノ酸の気相塩基性度の評価値は標準的な文献データと矛盾する。例えば、非特許文献３に基づけば、グリシン、アラニンやイソロイシンの気相塩基性度の評価値はお互いの差が１kJ/mol程度でしかないが、ＮＩＳＴで公開されている標準データ（非特許文献４参照）で比べるとその差がグリシンとアラニンとで１５kJ/mol程度、グリシンとイソロイシンとで３０kJ/mol程度、と大きく相違する。
〔２〕フラグメントイオンの気相塩基性度のフラグメントイオンの種類への依存性について例を挙げて言えば、同じアミノ酸配列であれば、ｂタイプのフラグメントイオンとｙフラグメントイオンの気相塩基性度は末端のＮＨ2とＣＯＯＨの寄与について若干考慮されているのみで、ｂタイプのフラグメントイオンの持つ特徴的なＣ末端構造による寄与を十分には反映できていない。
【００１８】
〔３〕(4)式はフラグメントイオンを構成するアミノ酸残基の種類の依存性のみを考慮し、アミノ酸残基の配列の依存性についてはあまり考慮していない。即ち、Ｃ末端側に隣り合うアミノ酸残基の影響を僅かに考慮しているだけである。しかしながら、これは明らかに一般性に欠けるものである。例えばプロリン（ピロリジン−２−カルボン酸）のように、他のアミノ酸と大きく異なり環の中にアミノ基の窒素を含む形態の分子構造を持つような特徴的な分子結合を示すアミノ酸残基が含まれている場合には、そのアミノ酸残基（例えばプロリン）がアミノ酸配列の中でどこに位置するのかによって、イオン化のエネルギー状態は異なる筈である。したがって、少なくとも上記のような特異なアミノ酸残基を含む場合には、アミノ酸残基の配列によって気相塩基性度も変わると考えるのが妥当である。
【００１９】
表１は、プロリン（Ｐ）を含む各種ジペプチドについて、プロリンがＮ末端側に存在する場合（Ｐ−Ｘのアミノ酸配列）とＣ末端側に存在する場合（Ｘ−Ｐのアミノ酸配列）のプロトン親和力（プロトンを引きつける力）ＰＡをそれぞれ分子動力学法による計算で求めた結果を示している。
【表１】
この表により、同じアミノ酸の組み合わせであってもプロリンの位置によってプロトン親和力が相違することは明らかであり、同様のことは気相塩基性度についても予想される。なお、この表において、Ｐ−ＸとＸ−ＰとのＰＡの値の開きは大きな場合で１０〜１５kcal/mol即ち４０〜６０kJ/molほどあり、このような大きな相違は非特許文献３による評価方法では説明できない。また非特許文献３に基づいて、Ｇ−Ｐ、Ａ−Ｐ、Ｉ−Ｐ、Ｌ−ＰやＶ−Ｐなどの気相塩基性度を評価しても、互いの評価値には殆ど差がなく１kJ/mol程度であるが、表１に示すＰＡの計算結果ではその差が最大１３kcal/mol即ち５０kJ/molほどにも及ぶ。
【００２０】
【非特許文献１】梶原、石田、竹内、「第一原理計算を利用したペプチドフラグメント予測法と定性分析への応用」、第５４回質量分析総合討論会、大阪、２００６
【非特許文献２】ラスキン(Julia Laskin)ほか１名、「ザ・セオリティカル・ベイシス・オブ・ザ・カイネマティック・メソッド・フロム・ザ・ポイント・オブ・ビュー・オブ・フィニット・ヒート・バス・セオリー(The Theoretical Basis of the Kinematic Method from the Point of the View of Finite heat bath Theory)」、ジャーナル・オブ・フィジックス・ケミストリー・エー(J. Phys. Chem. A)、2000、pp.8829-8837
【非特許文献３】ジャン(Z.Zhang)、「プレディケイション・オブ・ロー−エナジー・コリジョン−インデュースド・ディソシエイション・スペクトラ・オブ・ペプチドズ(Predication of Low-Energy Collision-Induced Dissociation Spectra of Peptides)」、アナリティカル・ケミストリー(Analytical Chemistry)、 Vol.76、No.14、July 15, 2004、 pp.3908-3922
【非特許文献４】「ナショナル・インスティテュート・オブ・スタンダーズ・アンド・テクノロジー（National Institute of Standards and Technology）」、ＮＩＳＴ、[平成１８年１１月１５日検索]、インターネット＜URL : http://webbook.nist.gov/＞
【発明の開示】
【発明が解決しようとする課題】
【００２１】
本発明は上記課題を解決するために成されたものであり、その目的とするところは、与えられたペプチドのフラグメントイオン強度を従来よりも正確に予測することにより、質量スペクトルデータに基づくタンパク質やペプチドのアミノ酸配列の決定精度を一層高めることができるペプチド同定方法を提供することである。
【課題を解決するための手段】
【００２２】
上記課題を解決するために成された本発明は、ＭＳｎ分析により取得された質量スペクトルデータに基づいて被検試料中のペプチド混合物のアミノ酸配列を決定するペプチド同定方法であって、候補ペプチドのフラグメントイオン強度を予測して質量スペクトルパターンを求め、これと実測の質量スペクトルとの比較により候補ペプチドの妥当性を判断するペプチド同定方法において、ペプチドのフラグメント化反応で生じる複数種のフラグメントイオン毎のフラグメントイオン強度を予測するために、
a)フラグメント化反応で生じるフラグメントイオンの強度を求める強度計算ステップと、
b)該フラグメントイオン強度をフラグメントイオン種類毎に与えられるプロトン親和力又は気相塩基性度に応じて分配する分配ステップと、
とを含み、フラグメントイオン種類毎のプロトン親和力又は気相塩基性度を計算する際に、フラグメントイオンを構成するアミノ酸残基のそれぞれについて実際の原子数にアミノ酸配列の特異性を反映した補正係数を乗じた有効原子数を求め、アミノ酸残基毎の有効原子数の和から求まる該フラグメントイオンの有効原子数を、所定の近似計算式に適用してプロトン親和力又は気相塩基性度を算出するようにしたことを特徴としている。
【００２３】
本発明に係るペプチド同定方法の一態様として、前記所定の近似計算式は、有効原子数とプロトン親和力又は気相塩基性度との非線形関係を対数で近似した式とすることができる。
【発明の効果】
【００２４】
或るアミノ酸配列のペプチド分子についてのフラグメント化反応で生じるフラグメントイオン種類（例えばｙタイプとｂタイプのフラグメントイオン）毎のフラグメントイオン強度を算出する際に、各フラグメントイオンのプロトン親和力（又は気相塩基性度）が必要となるが、従来、プロトン親和力や気相塩基性度といった値はかなり大まかにしか与えられなかった。これに対し本発明に係るペプチド同定方法では、プロトン親和力や気相塩基性度を従来よりも高い精度で、且つ実用的な計算量や計算時間で求める方法を提供する。
【００２５】
即ち、アミノ酸残基毎に原子数に適当な補正係数を乗じた有効原子数なるものを考え、アミノ酸配列に含まれるアミノ酸残基の種類がたとえ同じであっても、その配列順序によってプロトン親和力や気相塩基性度が影響を受けるような場合には、補正係数をそれぞれに合わせて修正することで有効原子数を変え、それによってプロトン親和力や気相塩基性度も変化するようにした。具体的には、例えば、アミノ酸配列中に塩基性アミノ酸が存在する場合に、該塩基性アミノ酸がＮ末端に存在するか否かで異なる補正係数を与えるようにするとよい。また、アミノ酸配列中にプロリンが存在する場合に、該プロリンと隣接する他のアミノ酸残基との組み合わせで該プロリンがＮ末端に存在するか否かで異なる補正係数を与えるようにするとよい。
【００２６】
また同一のフラグメント化反応から生じるフラグメントイオンでもそのタイプによって有効原子数を変えることで、プロトン親和力や気相塩基性度の精度をさらに上げることができる。例えばｂ−ｙフラグメント化反応ではｙフラグメントイオンとｂフラグメントイオンとが生成されるが、同じアミノ酸配列であっても後者は前者より原子数が３個少ないためにそれを有効原子数に反映させる。
【００２７】
従来、(4)式にあるように、気相塩基性度（プロトン親和力も同様）はフラグメントイオンの種類やアミノ酸残基の配列の影響を考慮していないものであったが、本発明に係るペプチド同定方法によれば、プロトン親和力や気相塩基性度をフラグメントイオンの種類やアミノ酸残基の配列の影響を考慮してより高い精度で求めることができる。それにより、ペプチドのフラグメント化反応で生じるフラグメントイオン種類毎にフラグメントイオン強度を高い精度で計算することができ、それ故に、これに基づいたペプチドの同定の信頼性を向上させることができる。
【００２８】
また、上記のような補正係数は特にプロトン親和力や気相塩基性度に対する影響が大きいものについてのみ異なる値を用意しておけばよく、さらに同定処理の際に一々計算するのではなく、予め計算しておいた値をデータベース化して同定処理の際にはデータベース検索により適当な値を得るようにすることができる。したがって、上記のような精度向上を図りながら、コンピュータ上での計算量や計算時間が膨大になることを回避でき、効率的な同定処理を実行してスループットの向上を図ることができる。
【発明を実施するための最良の形態】
【００２９】
本発明の一実施形態によるペプチド同定方法の全体の手順をまず説明する。図１は本実施形態のペプチド同定方法における解析処理の概略フローチャートである。実際には、この解析処理は、新規のプロテオーム解析用ソフトウエアをコンピュータ上で動作させることにより実行される。
【００３０】
まず分析対象の物質（ペプチド）を例えばイオントラップ飛行時間型質量分析装置で質量分析（ＭＳｎ分析）して質量スペクトルデータを収集する（ステップＳ１）。次いで、この質量スペクトルデータに対し、既存のプロテオーム解析用データベース検索ソフトウエア（例えばMASCOT）を利用して解析処理を実行する。即ち、質量スペクトルに出現しているフラグメントイオンの質量（厳密には質量電荷比、ｍ／ｚ値）をキーとして、データベースに登録されているペプチドの中から質量がマッチするペプチドを検索する（ステップＳ２）。そして、検索により得られた複数のペプチド分子をマッチ度順にランク付けしてリストアップする（ステップＳ３）。従来は、このリストアップされたペプチド分子の中から、分析者自らの判断（例えば最もスコアが高いものを選択する等）により同定を行う必要があったが、本発明に係る解析処理によれば、ステップＳ４以降の処理によりさらに候補の絞り込み（換言すればより確度の高い判断材料の提供）を行うことができる。
【００３１】
即ち、リストアップされた候補ペプチドの中で他と比べて明らかにスコアが高い複数の候補ペプチドを選定し、各候補のそれぞれのアミノ酸配列についてフラグメントイオン強度を計算することで質量スペクトルパターンを予測する。そのために、まずフラグメント化反応で生じるフラグメントイオン総量を計算し（ステップＳ４）、さらに計算により求まるプロトン親和力（又は気相塩基性度）に応じて同一のフラグメント化反応で生じる複数のフラグメントイオンへの強度の配分を行ってフラグメントイオン強度を求める（ステップＳ５）。こうして求められたフラグメントイオン強度から各候補ペプチドの質量スペクトルパターンを予測し、実測により得られた質量スペクトルとの相関による強度マッチング評価を行う（ステップＳ６）。そして、その相関性の評価結果を比べて最も可能性の高い候補ペプチドを選定して例えばモニタの表示画面上等に提示する（ステップＳ７）。このステップＳ４〜Ｓ７による新規の機能が、従来の標準的なデータベース検索ソフトウエアによる機能と組み合わせられ、これにより従来よりも信頼度の高い同定結果が期待できる。
【００３２】
本実施形態によるペプチド同定方法は、特にステップＳ４及びＳ５の処理に特徴を有している。これらについて以下に詳しく説明する。
【００３３】
（１）フラグメント化反応で生じるフラグメントイオン総量の計算方法
図３はフラグメントイオン強度の計算方法をエネルギー準位で示した原理説明図であり、（ａ）はプロトン親和力をベースにした、（ｂ）は気相塩基性度をベースにした原理説明図である。表現は異なるが、いずれも(1)式で示したアレニウスの反応式に基づいたフラグメントイオン強度計算式を利用している。
【００３４】
ＣＩＤ等のフラグメント化反応で生じるフラグメントイオンの総量は、(1)式に示したアレニウスの反応式を用いて求める。その際、パラメータの１つとして関連フラグメントの活性化エネルギーＥａを与える必要がある、活性化エネルギーは、ペプチドの初期状態から、着目しているフラグメントイオンの活性化状態、即ち遷移状態までのエネルギーの差として直接計算で求めることも可能であるが、ここでは二つの状態エネルギー差の和で考える。即ち、例えば図３（ａ）に示すように、初期状態であるプリカーサイオン（ＸＹ−Ｈ＋）の状態から開裂するペプチド結合位置（反応分子位置と呼ぶ）までプロトン（Ｈ＋）が移動するための移動エネルギーＥＨと、反応分子位置（Ｘ−Ｈ＋−Ｙ）から遷移状態（Ｘ−Ｈ＋−Ｙ＊）に移行するためのエネルギー変化量（以下、活性化エネルギーＥａ’と呼ぶ）とに分けて計算を行う。
【００３５】
このようにエネルギーを分割した上で、近似計算として、活性化エネルギーＥａ’は開裂するペプチド結合の両側に存在する２又は３個のアミノ酸配列分子モデルに簡単化して求める。このように簡単化できる理由は、フラグメント化反応がペプチド結合付近の局所的な反応であるとみて構わないためである。また、プロトンの移動エネルギーＥＨに関しては、移動の初めと最後にプロトンの存在するアミノ酸を組み合わせてできるジペプチドに分子モデルを簡単化して求める。このように分子モデルを簡単化できるのは、プロトンはそれ自身近くの原子と強く結合していて、遠くに存在する原子や分子との結合力は無視できるものとし、プロトンとの結合に直接関与しない部分を省略しても構わないためである。
【００３６】
上述したように２又は３個のアミノ酸配列分子モデルに簡単化して活性化エネルギーＥａ’やプロトンの移動エネルギーＥＨを求めるためには、予めアミノ酸配列の組み合わせ毎に活性化エネルギーＥａ’、プロトン移動エネルギーＥＨを計算してデータベース化しておき、結合を与えるとそのデータベースから対応したエネルギーＥａ’、ＥＨが導出されるようにしておくとよい。共通アミノ酸の種類はたかだか２０種類であるが、例えば５個のアミノ酸配列の組み合わせを考えただけでも、ｂフラグメントイオン、ｙフラグメントイオンの２種を考えると、エネルギーの計算対象は２０5×２＝６４×１０5と膨大な数となり、これを全て計算するのには膨大な時間を要する。これに対し、上記のように２又は３個のアミノ酸配列分子モデルに簡単化することで、計算対象の数を大幅に減らすことができ、コンピュータを用いて比較的短時間で計算することができるので、活性化エネルギーＥａ’、プロトン移動エネルギーＥＨを呼び出すためのデータベースを比較的容易に構築することができる。
【００３７】
なお、過去の研究により、フラグメント化反応は大別して、モバイルプロトンモデル又はリモートプロトンモデルと呼ばれる２つの反応メカニズムのいずれかで行われることが知られており、リモートプロトンモデルではフラグメント化反応においてプロトンが移動しない場合がある（それらの反応メカニズムの詳細については例えば特願２００６−４５５３号など参照）。この場合にはプロトンの移動エネルギーＥＨの計算は行わずに直接、活性化エネルギーＥａのみを計算するが、その場合でも計算方法はＥａ’の計算方法に準じればよい。
【００３８】
また、フラグメント化反応の過程をプロトン親和力でなく気相塩基性度で考える場合には、図３（ａ）と対比した図３（ｂ）に示すように、プロトンの移動エネルギーＥＨに代えたプロトンの移動自由エネルギー（ΔＧＨ）、Ｅａ’に代えた反応分子位置（Ｘ−Ｈ＋−Ｙ）から遷移状態（Ｘ−Ｈ＋−Ｙ＊）を経て分解するときの活性化自由エネルギー（ΔＧａ’）の和で活性化自由エネルギー（ΔＧａ）を表せばよい。
【００３９】
以上のようにして活性化エネルギーＥａ（又は活性化自由エネルギーΔＧａ）が求まれば、他のパラメータ、気体定数Ｒ、温度Ｔ、頻度因子Ａ（又はＡ０）は既知であるから、アレニウスの反応式によりフラグメントイオン強度、つまりはフラグメント化反応で同時に生じるＸイオンとＹイオンのフラグメントイオンの総量を算出することができる。
【００４０】
（２）複数種のフラグメントイオンへのプロトン親和力（又は気相塩基性度）に応じたフラグメントイオン強度の配分計算方法
質量スペクトルパターンを求めるには、着目しているフラグメント化反応で発生した複数種類のフラグメントイオンの発生量を種類毎（例えばｂフラグメントイオンとｙフラグメントイオン）に個別に求める必要がある。そこで、ステップＳ４で求めたフラグメントイオン総量をカイネティック法に基づく手法、即ち、フラグメントイオンのプロトン親和力ＰＡ（又は気相塩基性度ＧＢ）の大きさに応じて、exp｛ＰＡ／（Ｒ・Ｔeff）｝（又はexp｛ＧＢ／（Ｒ・Ｔeff）｝）の重みでもって各フラグメントイオンに配分する。
【００４１】
例えば代表的なフラグメント化反応プロセスであるｂ−ｙフラグメント化反応では、ｂフラグメントイオンとｙフラグメントイオンとが発生するが、ｂフラグメントイオン、ｙフラグメントイオンのプロトン親和力（又は気相塩基性度）をそれぞれＰＡｂ、ＰＡｙ（又はＧＢｂ、ＧＢｙ）としたとき、配分比は、
ｂフラグメントイオン配分比＝expＰＡｂ／（expＰＡｂ＋expＰＡｙ） …(5)
ｙフラグメントイオン配分比＝expＰＡｙ／（expＰＡｂ＋expＰＡｙ） …(6)
となる。気相塩基性度ＧＢを用いる場合には、上記(5)、(6)式においてＰＡｂ、ＰＡｙの代わりにＧＢｂ、ＧＢｙを使用すればよい。この配分方法自体はカイネテッィク法そのものであるが、この計算に用いるフラグメントイオンのプロトン親和力（又は気相塩基性度）を次のような特徴的な方法で求める。
【００４２】
過去の研究（ジョン・ホーメス、ほか２名(John L. Hohmes, Christiane Aubry, and Paul M. Mayer)、「プロトン・アフィニティーズ・オブ・プライマリー・アルカノールズ：アン・アプレイザル・オブ・ザ・カイネマティック・メソッド(Proton Affinities of Primary Alkanols: An Appraisal of the Kinematic Method)」、ジャーナル・フィジックス・ケミストリー(J. Phys. Chem.）、A,103, pp.705-709(1999)参照）によれば、同一系列の一級アルカノール(Primary Alkanols) 化合物について、原子数とプロトン親和力との間には一つの曲線関係が存在すると仮定し、この曲線を利用して同一系列に属する未知のアルカノールのプロトン親和力を予測している。
【００４３】
本実施形態によるペプチド同定方法でも基本的にはこの原理を利用してフラグメントイオンの原子数からプロトン親和力を求めるが、単なる原子数ではなくアミノ酸残基の種類や配列などを考慮した有効原子数という新たな概念を導入する。即ち、フラグメントイオンを構成するアミノ酸残基のそれぞれについて、原子数に補正係数（以下、スケーリング因子と呼ぶ）を乗じた有効原子数を与え、各アミノ酸残基の有効原子数の和をフラグメントイオンの有効原子数とする。また、フラグメントイオンのプロトン親和力（又は気相塩基性度）と有効原子数との非線形関係を対数式で表しておき、この対数式に有効原子数を代入することでプロトン親和力（又は気相塩基性度）を簡単に求められるようにする。
【００４４】
具体的には、例えばｙフラグメントイオンの場合について考えると、ｙフラグメントイオンは一般にｎ個のアミノ酸残基から成るペプチド分子イオンの構造を有しているので、次の(7)式により有効原子数Ｎを表す。
Ｎ＝Ｐ１・Ｌ１＋Ｐ２・Ｌ２＋…＋Ｐｎ・Ｌｎ＋３ …(7)
ここでＬｉ及びＰｉはそれぞれ、アミノ酸配列のＮ末端側から数えてｉ番目のアミノ酸残基の原子数及びスケーリング因子である。全てのスケーリング因子Ｐｉが１ならば、有効原子数Ｎはｙフラグメントイオンの原子数と同じとなる。ここで、(7)式の右辺右端の「３」は、Ｃ末端に結合する水酸基（−ＯＨ）とＮ末端に結合する水素基（−Ｈ）の原子数である。そして(7)式で定義した有効原子数Ｎを用いて、ｙフラグメントイオンのプロトン親和力を次の(8)式で表す。
ＰＡ＝ａ・ln｛ｂ・（Ｎ−３）｝＝ａ・ln｛ｂ・（Ｐ１・Ｌ１＋Ｐ２・Ｌ２＋…＋Ｐｎ・Ｌｎ）｝ …(8)
ここでａ、ｂは定数であり、実際のプロトン親和力の値とできるだけマッチングするように予めパラメータフィッティングで決められる。
【００４５】
なお、(8)式に示した対数式において、有効原子数Ｎをそのまま用いずに（Ｎ−３）を変数としているが、これはＮとしても構わない。但し、その場合には定数ａ、ｂを新たなパラメータフィッティングで決めればよい。
【００４６】
一方、ｂフラグメントイオンの場合には、次の(9)式によりプロトン親和力ＰＡｂを計算する。
ＰＡｂ＝ａ・ln｛ｂ・（Ｎ−６）｝＝ａ・ln｛ｂ・（Ｐ１・Ｌ１＋Ｐ２・Ｌ２＋…＋Ｐｎ・Ｌｎ−３）｝ …(9)
同じアミノ酸残基から成るｙフラグメントイオンとｂフラグメントイオンとでは、ｂフラグメントイオンのほうが原子数が３個少ないため、(9)式ではこれを考慮した式としている。
【００４７】
上記(8)、(9)式においては、フラグメントイオンを構成するアミノ酸残基の種類には依存するものの、その序列には依存しない形式としてプロトン親和力を定義している。但し、プロトン親和力の値がアミノ酸残基の序列にも多少の影響を受けることは、従来から知られている。特にアミノ酸配列の両端であるＮ末端、Ｃ末端は分子構造が途切れることからプロトン親和力が影響を受け易いと考えられ、その中でもＮ末端はモバイルプロトンモデルにおいては、フラグメント化反応の初期にプロトンのいる分子位置であることから、このＮ末端にプロリンや塩基性アミノ酸などの特定のアミノ酸が存在するとプロトン親和力が無視できない程度に変化することが予想される。
【００４８】
そこで、特に分子構造の途切れるＮ末端での影響の大きさを考慮し、Ｎ末端に塩基性アミノ酸が存在する場合にはＮ末端にそうしたアミノ酸が無い場合とは異なるスケーリング因子を割り当てるようにする。また、プロリンのような他のアミノ酸とは異なる特徴的な分子構造を持つアミノ酸については、隣接するアミノ酸との組み合わせによってもプロトン親和力が影響を受けると考えられる。そこで、プロリンのような特異なアミノ酸が存在する場合には、そのプロリンと隣接するアミノ酸との組み合わせ一つ一つについて異なるスケーリング因子を与えるものとする。このようにスケーリング因子が変わると、仮に元の原子数は同じであっても有効原子数は変わり、それに伴いプロトン親和力も変化することになる。
【００４９】
気相塩基性度ＧＢに対してもプロトン親和力ＰＡと同様に表すことができる。即ち、ｙフラグメントイオン及びｂフラグメントイオンについてそれぞれ次式で気相塩基性度ＧＢｙ、ＧＢｂを表すものとする。
ＧＢｙ＝ａ’・ln｛ｂ’・（Ｐ１’・Ｌ１＋Ｐ２’・Ｌ２＋…＋Ｐｎ’・Ｌｎ）｝ …(10)
ＧＢｂ＝ａ’・ln｛ｂ’・（Ｐ１’・Ｌ１＋Ｐ２’・Ｌ２＋…＋Ｐｎ’・Ｌｎ−３）｝ …(11)
ここでａ’、ｂ’は定数であり、実際の気相塩基性度とできるだけマッチするようにカーブフィッティングで予め決めておく。Ｐｉ’、Ｌｉはそれぞれｉ番目のアミノ酸残基の原子数、スケーリング因子である。そしてプロトン親和力と同様に、特定のアミノ酸がＮ末端にある場合とそうでない場合とで異なるスケーリング因子を割り当てて気相塩基性度が変化するようにする。またプロリンについても同様である。
【００５０】
具体的には、(8)〜(11)式で用いられるスケーリング因子Ｐｉ（又はＰｉ’）は定数ａ、ｂ（又はａ’、ｂ’）とともに実験データとマッチングするように決められるべき定数であるから、これらの値は多くの観測データと比較して最小二乗法により予め求めておくことができる。そうした求めた値を例えばテーブル形式でデータベース化しておき、与えられたアミノ酸配列に応じて適当な値を読み出して来て上記のような計算式を用いて、フラグメントイオン毎にプロトン親和力（又は気相塩基性度）を算出することができる。
【００５１】
また、別の方法でスケーリング因子Ｐｉを決めることもできる。即ち、ｎ＝１の場合、ｙフラグメントイオンはアミノ酸イオンであることに着目すれば、アミノ酸のプロトン親和力はよく知られているので、逆にアミノ酸のプロトン親和力からスケーリング因子を求めることができる。この場合は、(8)式においてｉ番目のアミノ酸残基一つのみから成るｙフラグメントイオン（＝アミノ酸イオン）を想定し、そのアミノ酸イオンのプロトン親和力をＰＡｉと置けば、(8)式は
ＰＡｉ＝ａ・ln（ｂ・Ｐｉ・Ｌｉ） …(12)
つまり、
Ｐｉ＝｛１／（ｂ・Ｌｉ）｝exp（ＰＡｉ／ａ） …(13)
と変形することができる。したがって、実験や理論的な考察などから得られるアミノ酸イオンのプロトン親和力ＰＡｉの値を(13)式に代入すれば、スケーリング因子Ｐｉを求めることができる。
【００５２】
また(8)式を変形すると、
exp（ＰＡ／ａ）＝ｂ・（Ｐ１・Ｌ１＋Ｐ２・Ｌ２＋…＋Ｐｎ・Ｌｎ） …(14)
のようになる。この式に上記(13)式を代入すれば、
exp（ＰＡ／ａ）＝exp（ＰＡ１／ａ）＋exp（ＰＡ２／ａ）＋…＋exp（ＰＡｎ／ａ） …(15)
となる。この式は形式上は(4)式と似ているが、意味するところは全く相違する。即ち、(15)式の指数因子の分母ａは定数であり、一方、(4)式の指数因子の分母Ｒ・Ｔeffで有効温度Ｔeffはペプチド分子の質量電荷比、ｍ／ｚ値等に依存する値であって定数ではない（上記非特許文献３参照）。
【００５３】
なお、必ずしも(13)式を用いてスケーリング因子Ｐｉを決める必要はなく、できるだけ実際のフラグメントイオンのプロトン親和力（又は気相塩基性度）にマッチするようにパラメータフィッティング等で決めればよい。
【００５４】
上述の方法でフラグメントイオン毎のプロトン親和力を算出する際の具体例を説明する。表２はポリグリシン（Ｇ、ＧＧ、ＧＧＧ、ＧＧＧＧ、ＧＧＧＧＧ）やポリアラニン（Ａ、ＡＡ、ＡＡＡ、ＡＡＡＡ）など、分子量の小さなアミノ酸の配列から構成されるペプチドのプロトン親和力を、Ｂ３ＬＹＰ／６−３１１＋Ｇ（２ｄ，ｐ）のレベルでの密度汎関数法を用いた分子動力学計算で求めた結果を示している。
【表２】
【００５５】
この結果を利用して、(8)式にあるパラメータａ、ｂを、横軸を原子数Ｎ−３、縦軸をプロトン親和力とした曲線のカーブフィッティングから決めることができる。即ち、グリシン（Ｇ）やアラニン（Ａ）などのアミノ酸のスケーリング因子Ｐｉをともに１であると仮定してカーブフィッティングを行うと、その結果は、
ＰＡ[kcal/mol]＝１１．３９５・ln（Ｎ−３）＋１９０．８６ …(16)
と得られる。これを(8)式の形に変形すれば、定数ａ、ｂは、
ａ＝１１．３９５[kcal/mol]、ｂ＝１．８８０１６×１０７
と求まる。
【００５６】
表３はグリシン、アラニン以外の各種アミノ酸に対して上記定数ａ、ｂの値を用い、アミノ酸のプロトン親和力の値からスケーリング因子を求めた例である。ここで示していないアミノ酸についても同様にしてスケーリング因子を求めることができる。
【表３】
このように定数ａ、ｂ及びスケーリング因子Ｐを決めておけば、(8)式を用いて任意のｙフラグメントイオンのプロトン親和力を、(9)式を用いて任意のｂフラグメントイオンのプロトン親和力をそれぞれ予測値として計算することができる。
【００５７】
表４は上記結果を利用して２個のアミノ酸から成るｙフラグメントイオン、換言すればジペプチドイオンに対するプロトン親和力ＰＡの予測値を算出した例である。この例では、アラニン（Ａ）を含むジペプチドのプロトン親和力を予測しているが、Ｂ３ＬＹＰ／６−３１１＋Ｇ（２ｄ，ｐ）のレベルでの密度汎関数法を用いた分子動力学計算から求めた結果と比較的良好に一致しており、上述のような本実施形態におけるプロトン親和力の予測方法が実用的に有効であることを確認することができる。
【表４】
【００５８】
表５は、(9)式を利用して２個のグリシン（Ｇ）から成るｂフラグメントイオンのプロトン親和力を予測し、それを分子動力学による計算値と比較した結果である。この場合にも両者は比較的良好に一致していることが分かる。
【表５】
【００５９】
前述のように表１はプロリン（Ｐ）とその他のアミノ酸の組み合わせによるジペプチドについてプロトン親和力の分子動力学的計算値を示したものであるが、この表１の結果から求められたプロリンのスケーリング因子を表６に示す。
【表６】
プロリン以外のアミノ酸残基のスケーリング因子は既に決定されているものとし、表１に示したジペプチドのプロトン親和力を(8)式に代入することにより、プロリンのスケーリング因子を決定することができる。
【００６０】
Ｎ末端にプロリンが存在する場合のプロリンのスケーリング因子は、表１中のＰ−Ｘ態様のジペプチドのプロトン親和力の値を(8)式に代入して得ることができる。そして、表６中の左側に示したスケーリング因子が、プロリンと結合するアミノ酸残基の種類に応じて与えられる。またプロリンがＮ末端に無い場合のスケーリング因子は、表１中のＸ−Ｐ態様のジペプチドのプロトン親和力の値を(8)式に代入して得ることができ、表６中の右側に示したスケーリング因子が、プロリンと結合するアミノ酸残基の種類に応じて与えられる。
【００６１】
以上例示したように、アミノ酸配列に応じて適宜のスケーリング因子を決定することができる。図２は実際にコンピュータでフラグメントイオン強度の配分計算を行う際の処理手順を示すフローチャートである。即ち、評価対象のペプチド分子が与えられると、そのアミノ酸配列に応じて所定のアルゴリズムに従ってスケーリング因子が決定される（ステップＳ１１）。具体的には、例えばアミノ酸配列が入力されると、例えば上述したようにスケーリング因子に影響を与える条件が抽出され、その条件に応じたデータベース検索によりアミノ酸残基毎のスケーリング因子が求まる。次に、各アミノ酸残基のスケーリング因子と原子数とから、(8)、(9)式に基づいて各フラグメントイオンのプロトン親和力を計算する（ステップＳ１２）。或いは、(10)、(11)式に基づいて各フラグメントイオンの気相塩基性度を計算してもよい。
【００６２】
そうして各フラグメントイオンのプロトン親和力又は気相塩基性度が求まったならば、例えばｂ−ｙフラグメントであれば(5)、(6)式に基づいてフラグメントイオン強度を配分する（ステップＳ１３）。もちろん、ｂ−ｙフラグメントでなくても同様に配分できることは当然である。
【００６３】
以上説明したように本発明に係るペプチド同定方法によれば、よく知られた化学反応の標準的な理論、即ちアレニウスの反応速度式やカイネテッィク法に沿う形で、フラグメントイオン強度を簡単なモデル式により近似計算で求めることができる。こうしたモデル式を用いず直接コンピュータでフラグメントイオン強度を求めようとすると、一般にペプチドのアミノ酸配列のサイズは大きいため、計算量が膨大になって実用的な時間での処理は困難になる。これに対し、上記実施形態で採用しているような方法によれば、計算に必要な活性化エネルギーやプロトン親和力（又は気相塩基性度）をペプチド分子に対して効率的に、つまり実用的な時間の範囲で求めることが可能となる。特にプロトン親和力（又は気相塩基性度）はアミノ酸配列のサイズが大きくなるほど大きくなる傾向にあり、一定の値に収束する傾向はない。そのため、プロトン親和力や気相塩基性度を簡単な評価式から少ない計算量で求めることができることは非常に有効であると言える。
【００６４】
また、各フラグメントイオンのプロトン親和力を計算する際に、そのフラグメントイオンの種類やアミノ酸配列の特異性（例えばプロリンの存在や存在位置など）を考慮して値が修正されるので、より高い精度でもってフラグメントイオン強度を求めることができ、ひいてはペプチドの同定の正確性の向上が期待できる。
【００６５】
なお、上記実施形態は本発明の一例にすぎず、本発明の趣旨の範囲で適宜変形、修正、追加等を行っても本願特許請求の範囲に包含されることは当然である。
【図面の簡単な説明】
【００６６】
【図１】本発明の一実施形態によるペプチド同定方法における解析処理の概略フローチャート。
【図２】フラグメントイオン強度の配分計算を行う際の処理手順を示すフローチャート。
【図３】フラグメントイオン強度の計算方法を反応経路に沿ったエネルギープロファイルで示した原理説明図。
【特許請求の範囲】
【請求項１】
ＭＳｎ分析により取得された質量スペクトルデータに基づいて被検試料中のペプチド混合物のアミノ酸配列を決定するペプチド同定方法であって、候補ペプチドのフラグメントイオン強度を予測して質量スペクトルパターンを求め、これと実測の質量スペクトルとの比較により候補ペプチドの妥当性を判断するペプチド同定方法において、ペプチドのフラグメント化反応で生じる複数種のフラグメントイオン毎のフラグメントイオン強度を予測するために、
a)フラグメント化反応で生じるフラグメントイオンの強度を求める強度計算ステップと、
b)該フラグメントイオン強度をフラグメントイオン種類毎に与えられるプロトン親和力又は気相塩基性度に応じて分配する分配ステップと、
を含み、フラグメントイオン種類毎のプロトン親和力又は気相塩基性度を計算する際に、フラグメントイオンを構成するアミノ酸残基のそれぞれについて実際の原子数にアミノ酸配列の特異性を反映した補正係数を乗じた有効原子数を求め、アミノ酸残基毎の有効原子数の和から求まる該フラグメントイオンの有効原子数を、所定の近似計算式に適用してプロトン親和力又は気相塩基性度を算出するようにしたことを特徴とする、質量分析を用いたペプチド同定方法。
【請求項２】
前記所定の近似計算式は、有効原子数とプロトン親和力又は気相塩基性度との非線形関係を対数で近似した式であることを特徴とする、請求項１に記載の質量分析を用いたペプチド同定方法。
【請求項３】
アミノ酸配列中に塩基性アミノ酸が存在する場合に、該塩基性アミノ酸がＮ末端に存在するか否かで異なる補正係数を与えることを特徴とする、請求項１又は２に記載の質量分析を用いたペプチド同定方法。
【請求項４】
アミノ酸配列中にプロリンが存在する場合に、該プロリンと隣接する他のアミノ酸残基との組み合わせで該プロリンがＮ末端に存在するか否かで異なる補正係数を与えることを特徴とする、請求項１又は２に記載の質量分析を用いたペプチド同定方法。
【請求項１】
ＭＳｎ分析により取得された質量スペクトルデータに基づいて被検試料中のペプチド混合物のアミノ酸配列を決定するペプチド同定方法であって、候補ペプチドのフラグメントイオン強度を予測して質量スペクトルパターンを求め、これと実測の質量スペクトルとの比較により候補ペプチドの妥当性を判断するペプチド同定方法において、ペプチドのフラグメント化反応で生じる複数種のフラグメントイオン毎のフラグメントイオン強度を予測するために、
a)フラグメント化反応で生じるフラグメントイオンの強度を求める強度計算ステップと、
b)該フラグメントイオン強度をフラグメントイオン種類毎に与えられるプロトン親和力又は気相塩基性度に応じて分配する分配ステップと、
を含み、フラグメントイオン種類毎のプロトン親和力又は気相塩基性度を計算する際に、フラグメントイオンを構成するアミノ酸残基のそれぞれについて実際の原子数にアミノ酸配列の特異性を反映した補正係数を乗じた有効原子数を求め、アミノ酸残基毎の有効原子数の和から求まる該フラグメントイオンの有効原子数を、所定の近似計算式に適用してプロトン親和力又は気相塩基性度を算出するようにしたことを特徴とする、質量分析を用いたペプチド同定方法。
【請求項２】
前記所定の近似計算式は、有効原子数とプロトン親和力又は気相塩基性度との非線形関係を対数で近似した式であることを特徴とする、請求項１に記載の質量分析を用いたペプチド同定方法。
【請求項３】
アミノ酸配列中に塩基性アミノ酸が存在する場合に、該塩基性アミノ酸がＮ末端に存在するか否かで異なる補正係数を与えることを特徴とする、請求項１又は２に記載の質量分析を用いたペプチド同定方法。
【請求項４】
アミノ酸配列中にプロリンが存在する場合に、該プロリンと隣接する他のアミノ酸残基との組み合わせで該プロリンがＮ末端に存在するか否かで異なる補正係数を与えることを特徴とする、請求項１又は２に記載の質量分析を用いたペプチド同定方法。
【図１】
【図２】
【図３】
【図２】
【図３】
【公開番号】特開２００８−１２８８４９（Ｐ２００８−１２８８４９Ａ）
【公開日】平成２０年６月５日（２００８．６．５）
【国際特許分類】
【出願番号】特願２００６−３１４９２２（Ｐ２００６−３１４９２２）
【出願日】平成１８年１１月２２日（２００６．１１．２２）
【出願人】（０００００１９９３）株式会社島津製作所 (3,708)
【出願人】（５０５１９５３８４）国立大学法人奈良女子大学 (15)
【Ｆターム（参考）】
【公開日】平成２０年６月５日（２００８．６．５）
【国際特許分類】
【出願日】平成１８年１１月２２日（２００６．１１．２２）
【出願人】（０００００１９９３）株式会社島津製作所 (3,708)
【出願人】（５０５１９５３８４）国立大学法人奈良女子大学 (15)
【Ｆターム（参考）】
[ Back to top ]