(57)【要約】【課題】質量スペクトル上の多数のピークから高速・高精度に糖鎖の配列を決定する。【解決手段】分析対象となる糖鎖および開裂施行によって得た複数の糖鎖断片に開裂施行を含む質量分析を行い、得られた複数の質量データとあらかじめ求められている各種糖鎖の構造または配列データとを照合し糖鎖配列を決定すること、もしくは質量スペクトル中の質量ピーク自身の組合せにより、糖鎖配列を直接的に決定することにおいて発生する膨大な配列決定演算を糖鎖の生合成における配列・構造特性ならびに開裂前の分析質量における制約を用いて減少させることで、実測の質量スペクトルが示す糖鎖配列を決定する。【特許請求の範囲】
【請求項１】
分析対象となる糖鎖および開裂施行により得た複数の糖鎖断片に対して開裂試行を含む質量分析を行い、得られた複数の質量スペクトルをあらかじめ求められている各種糖鎖の構造または配列データと照合し糖鎖配列を決定すること、もしくは糖鎖を構成する単糖の質量と質量スペクトル中の質量ピーク自身の組合せにより糖鎖配列を直接的に決定することにおいて発生する膨大な配列決定演算を糖鎖の生合成特性による制約を用いて減少させることを特徴とする実測の質量スペクトルが示す糖鎖配列を決定する方法。
【請求項２】
請求項１の膨大な配列決定演算を質量分析における開裂前の分析質量による制約を用いて減少させることを特徴とする実測の質量スペクトルが示す糖鎖配列を決定する方法。
【請求項３】
上記の制約の適用により、開裂スペクトルを測定するピークを限定する方法。
【発明の詳細な説明】
【技術分野】
【０００１】  本発明は、質量分析を用いた糖鎖配列解析の方法に関する。
【背景技術】
【０００２】  現在、生命現象の解明にはゲノムやタンパク質だけでなく、糖鎖まで含めた研究が必要であるとの認識が高まってきている。例えば、タンパク質が生体内で機能する上では糖鎖が介在してタンパク質の機能を調整していることが知られており（非特許文献１）、生体内におけるタンパク質の約５０％に糖鎖が修飾された状態で存在していると推測されている（非特許文献２）。【０００３】  糖鎖が生体内で重要かつ多岐に渡る機能を果たす理由として、糖鎖の構造多様性が挙げられる。例えば、核酸は４種類、アミノ酸は２０種類あるがその結合様式は１種類しかないので、核酸２塩基の配列は４×４＝１６通り、ジペプチドの配列は２０×２０＝４００通りである。それに対し糖鎖は、ヘキソース（Ｈｅｘ）ではグルコース（Ｇｌｃ）、ガラクトース（Ｇａｌ）、マンノース（Ｍａｎ）の３種しか生体内には存在しないが、α型、β型のアノマーが存在し、かつ各々が結合可能な水酸基を５つ有する。Ｈｅｘのみを検討しても考えられる結合様式も含めた糖鎖配列は多様であり、また、実際に糖鎖を構成する単糖はＨｅｘに限らない。更に、核酸、タンパク質は直鎖分子であるが、糖鎖は分岐構造を有するものも存在する為、糖鎖分析は核酸やタンパク質の研究と比べて難しく、解析に多くの時間が費やされている。【０００４】  糖鎖を開裂施行を伴う質量分析により測定する場合、質量分析器内においての開裂は不均一に起こるわけでなく、単糖間のグリコシル結合、もしくはピラノース環の特定位置において開裂が起き易いことが知られている（図１）。また、開裂位置と開裂が起こる確率がある程度均一であれば、質量分析器内で開裂を繰り返すことで、糖鎖を構成する単糖全てのフラグメント質量が検出される。例えば、図２に示したアスパラギン（Ａｓｎ）結合型糖鎖（Ｎ‐結合型糖タンパク質糖鎖）の開裂を含む質量分析（ＭＳｎ）の分析結果においては、初めの開裂（ＭＳ２）により糖鎖の各末端を含む糖鎖フラグメント、つまり、還元末端側のＮ‐アセチルグルコサミン（ＧｌｃＮＡｃ）を含んだ糖鎖配列とそれ以外の非還元末端側の糖鎖配列が検出され、更なる開裂（ＭＳ３）により糖鎖フラグメント自身の開裂を起こす。したがって、開裂により検出されるフラグメント質量に対応する単糖の組合せを求めること、または既知の糖鎖データベースと比較検討することで糖鎖配列の推定が可能となる。【０００５】  配列決定において糖鎖データベースと比較するべく、質量分析ではあらゆるピークに対してＭＳｎ分析を行うが、測定者による開裂対象ピークの選定に掛かる時間や無秩序に行われる多数の開裂を伴う質量分析が分析時間を増大させるばかりでなく、生成される分析データをも増大させ、配列解析における計算量が増大している。【０００６】  フラグメント質量に対応する糖鎖配列を単糖の組合せより求める方法ならびに、糖鎖データベースと比較する方法は、タンパク質の質量分析による同定ならびに配列決定法（特許文献１、２）におけるタンパク質を糖鎖として、アミノ酸を単糖として考えることで、同等の手法により比較評価が可能である。糖鎖を構成する単糖が３種類のＨｅｘからなる２糖の配列組合せであると考えるだけでも、２つのアノマーと５つの結合可能な部位（水酸基）より、（３（単糖）×２（アノマー）×５（結合部位））２＝９００通りもの結合様式も含めた糖鎖配列が想定され、実際の糖鎖配列において考えられる配列長と分岐構造を含めると、莫大な数の組合せが発生し、データベースとの比較検討が難解になる。【０００７】
【非特許文献１】Ｖａｒｋｉ  Ａ，ｅｔ  ａｌ．，Ｇｌｙｃｏｂｉｏｌｏｇｙ，３，９７‐１３０（１９９３）
【非特許文献２】Ａｐｗｅｉｌｅｒ  Ｒ，ｅｔ  ａｌ．，Ｂｉｏｃｈｉｍ  Ｂｉｏｐｈｙ  ｓＡｃｔａ，１４７３，４‐８（１９９９）
【特許文献１】特開２００６−１６２５５６
【特許文献２】特開平１１−２３７３８３
【発明の開示】
【発明が解決しようとする課題】
【０００８】  本発明が解決しようとする課題は、質量スペクトル用いた糖鎖配列決定において必要とする質量分析を限定し、かつ開裂試行を含む単糖の質量データもしくは糖鎖のデータベースとの比較検討の上で発生する莫大な糖鎖配列の組合せ計算を制限することで、高速・高精度に糖鎖配列を決定する方法を提供することである。
【課題を解決する為の手段】
【０００９】  上記課題を克服するため、本発明の糖鎖配列決定法では、糖鎖および開裂施行によって得た複数の糖鎖断片に複数回の開裂施行を含む質量分析を行い、得られた複数の質量スペクトルをあらかじめ求められている各種糖鎖配列データベースと照合し、対応する配列データを有する糖鎖に各質量スペクトル中の質量ピークを割り当てる。それとともに、割り当てられた各質量ピークの開裂施行によって得た質量スペクトル中の質量ピークと配列データとの一致状態を求める。複数の糖鎖配列データに同様の質量データの割当てを行い、開裂による質量データと配列データとの一致状態が高い糖鎖を選出する。
【００１０】  また本発明は、得られた複数の質量スペクトルを糖鎖配列データベースではなく、糖鎖を構成する単糖の質量の組合せと照合し、割り当てられた各質量ピークの開裂によるフラグメント質量のデータと配列データとの一致状態を求め、複数の単糖の質量の組合せに同様の質量データの割り当てを行い、開裂によるフラグメント質量のデータと配列データとの一致状態が高い糖鎖を選出することも含む。
【００１１】  得られた複数の質量スペクトルに対して、糖鎖配列データベースとの照合、または糖鎖を構成する単糖の質量の組合せと照合を問わず、質量スペクトルを導出したプレカーサーイオンから糖鎖配列を構成する単糖の組合せを予測する。
【００１２】  配列解析における質量分析ではあらゆるピークに対して開裂を含む質量分析を行うが、測定者による開裂対象ピーク選定に掛かる時間や無秩序に行われる開裂を伴う多数の質量分析が分析時間を増大させるばかりでなく、分析データを増大させ、配列解析における計算量を増大させている。その為、この単糖組合せによる配列予測を質量分析に反映させることで、開裂試行を用いて測定するべきピークを選択し、解析時間を増大させる質量分析を減少させる。
【００１３】  糖鎖および糖鎖断片の開裂による質量スペクトルの開裂前のプレカーサー質量は、糖鎖を構成するいくつかの単糖の組合せによって算出される質量もしくはそれにアグリコン（糖鎖成分以外の部分）の質量を加えたと一致する必要がある。単糖を組合せて、開裂前の質量を超える、もしくは不足する質量と組合せは予測する組合せに該当しない。この単糖の組合せを予測する具体的な手法は、ナップサック問題を代表とする制限付き組合せ問題の解法に該当する。
【００１４】  単糖の組合せの予測においては、糖鎖の生合成における制約を同時に用いる。コア構造は複合糖質の種類によって異なるが生体内での糖鎖のコア構造に用いられる糖鎖配列は、解析対象の糖鎖に含まれるという前提で予測する単糖の組合せを計算する。
【００１５】  生体内における糖鎖のコア構造に対し、決められた分岐もしくは決められた単糖の決められた結合が許されることを糖鎖の生合成における制限として単糖の組合せを予測する。制限によって単純に質量スペクトルを導いたプレカーサーイオンの質量を満たす単糖の組合せを網羅的に列挙する必要がなくなり、高速に予測する単糖の組合せが計算できる。
【００１６】  〔００１１〕から〔００１５〕によって予測された単糖の組合せより、糖鎖配列データベースから該当する糖鎖配列を選択すること、もしくは予測された単糖の組合せと質量スペクトルとを直接的に比較することで、膨大な照合における計算を削減し糖鎖配列を決定する。
【発明の効果】
【００１７】開裂試行を含んだ糖鎖の質量分析により導出される膨大な数の質量ピークと糖鎖データベースとの比較（照合）または質量データより単糖の質量の組合せを行う上で障害となる膨大な計算量を、本発明による糖鎖の生合成特性と質量による制約を用いることで大幅に削減し、糖鎖の配列決定を可能にする。
【００１８】  また、生合成特性と質量による制約による配列予測を測定にフィードバックすることにより、開裂試行を実施する必要のある質量ピークの選択的な制限を行い、配列決定において必要最小限となる質量スペクトルの測定ならびに計算を可能にし、糖鎖配列決定の高速化を実現する。
【発明を実施するための最良の形態】
【００１９】  以下、本発明を更に詳細に説明するが、以下の構成要件の説明は本発明の実施態様の代表例であり、本発明はこれらの内容のみに特定されるものではない。
本発明を実施する上で必要となる測定糖鎖は、あらかじめ還元末端に特定の修飾を行うことで、質量分析前に測定サンプル中のアグリコンの分子量を明確にする。Ａｓｎ結合型糖鎖（Ｎ‐結合型糖タンパク質糖鎖）を例に挙げると、還元末端側のＧｌｃＮＡｃのアノマー位にピリジルアミノ基等を導入する形で修飾を行い、測定サンプル中の糖鎖部分の質量とアグリコンの質量が把握できる状態で質量分析を実施。【００２０】  測定サンプルに対しては、ラベル修飾を行ったサンプル全体の質量電荷比、つまりプレカーサーイオンを測定すると共に、その測定サンプルに対して質量分析装置により定義されるエネルギーを加えることで測定サンプルを開裂させた質量スペクトル、またその開裂スペクトルを更に開裂させて得られる質量スペクトルといった他段階の開裂試行を含む質量分析の実施。
【実施例】
【００２１】  以下、実施例を挙げて本発明を詳細に説明するが、本発明の範囲はこれらの実施例により限定されるものではない。
図１に示されたＡｓｎ結合型糖鎖を例として、配列解析の方法を説明する。糖鎖を質量分析する場合、あらかじめ還元末端側のＧｌｃＮＡｃのアノマー位に対してピリジルアミノ基等の修飾によるラベル化を実施することが一般的であり、糖鎖と測定者があらかじめ知ることの出来るラベルが結合したサンプルを測定することとなる。本例においては、そのラベルを省略し、末端がＡｓｎである糖鎖を測定した場合の解析を示している。開裂を伴わない質量分析（ＭＳ）実施では、図１に示される分子全体の質量にナトリウムイオンが付加した質量電荷比（ｍ／ｚ）＝１７７７．７が得られ、その質量電荷比を質量に換算し、Ａｓｎの分子量を差し引くことで糖鎖全体の分子量が１６４０Ｄａであると決定できる。【００２２】  図２では、Ａｓｎ結合型糖鎖の糖鎖生合成のルールを制約として纏めている。この制約２により与えられる還元末端側の５つの単糖からなる糖鎖がＡｓｎ結合型糖鎖のコア構造（９１０．５Ｄａ）であり、その構造以外に図１に含まれる単糖の集合全体は１６４０．０Ｄａ−９１０．５Ｄａ＝７３０．４Ｄａの計算結果により、２つのＨｅｘと２つのＮ‐アセチルヘキソサミン（ＨｅｘＮＡｃ）が基本構造に付加された図３の配列が測定された糖鎖配列の候補であると判断できる。【００２３】  また、得られた図１の全体のｍ／ｚを示すピークに対して開裂を含む質量分析を行った（ＭＳ２）の結果である図４中のｍ／ｚ＝１４４２．７のピークは、制約２による基本構造から、図５に示す還元末端側に存在する２つのＧｌｃＮＡｃ間のグリコシル結合が開裂により、Ａｓｎとそれに結合するＧｌｃＮＡｃが欠如した非還元末端側のピークと考えられる。このピークは配列を求める上で明らかになっていない２つのＨｅｘと２つのＨｅｘＮＡｃが含まれる糖鎖断片であるため、この質量電荷比をもつピークを選択的に次の開裂試行を用いた質量分析（ＭＳ３）の対象として質量分析を実施する。【００２４】  図６はＭＳ３の分析結果を示している。ｍ／ｚ＝１４４２．７のピーク（Ｓ１１）は、前述のＡｓｎと還元末端側のＧｌｃＮＡｃが一つ欠失した図５の糖鎖を示すピークであり本開裂においてフラグメント化できなかったと考えられるため、開裂前の構造と同等である。Ｓ１はＨｅｘ２つ分を示しており、２つのＭａｎ、２つのＧａｌ、もしくはＭａｎとＧａｌの組合せであると考えられる。Ｓ３のピークはＨｅｘ３つ分を示しており、制約２のコア構造内にある３つのＭａｎの分岐構造、またＳ１はそのＭａｎのうちの連続する２つに由来するピークである可能性が示唆される。Ｓ３のピークからさらにＨｅｘ分を加えたピーク、つまりＨｅｘ４つ分の６４８Ｄａがイオン化したピークおよびＨｅｘ５つ分８１０Ｄａに相当するピークが存在しないため、コア構造のＭａｎの分岐に繋がる単糖はＨｅｘではなく、また、図２の制約３よりＧｌｃＮＡｃが繋がる複合型のＡｓｎ結合型糖鎖と考えられる。よって、図３ｃおよびｅは想定できなくなる。【００２５】  ここで、図３ａのＭａｎの分岐後に結合した、２つのＨｅｘと２つのＨｅｘＮＡｃを考える。ＨｅｘＮＡｃが連続した配列と連続しない配列が考えられるが、２つのＨｅｘＮＡｃに相当するフラグメントピークが得られないことから���者が、また、２つのＨｅｘと２つのＨｅｘＮＡｃに相当するフラグメントピークが無いことより後者の配列が想定できなくなる。さらに、図３ｄのＭａｎの分岐後に結合した、２つのＨｅｘと１つのＨｅｘＮＡｃに関しては、コア構造内の３つのＭａｎと２つのＨｅｘ、１つのＨｅｘＮＡｃを含んだピークが存在しない為、図３ｄの配列を想定することも難しい。【００２６】  よって図３ｃが配列解析候補として、単糖の配列としての順序を考える。前述の通り４つのＨｅｘによるピークが存在しない為、コア構造中の３つのＭａｎ以降はＨｅｘＮＡｃであり、制約３ｉｉよりＧｌｃＮＡｃと決定できる。ＧｌｃＮＡｃに続くＨｅｘは制約３ｉ，ｉｉｉよりＭａｎではなく、Ｇａｌであると考えられる。【００２７】  上記により予想される配列は、図７のようになる。データベースを参照することなく予想配列を一つに絞り込めたが、図７をデータベースと照らし合わせることで予想される配列の確定が行える。また配列の候補が複数ある場合は、候補配列のみをデータを検索することで配列決定は実施可能である。【００２８】  以上の一連の解析の流れを図８のフローチャートとして示す。
【図面の簡単な説明】
【００２９】
【図１】実施例において配列決定を行ったＡｓｎ結合型糖鎖を示した図である。
【図２】Ａｓｎ結合型糖鎖における生合成を用いた制約を示した図である。
【図３】制約の検討により考えられる配列の組合せを示した図である。
【図４】開裂を伴う質量分析（ＭＳ２）の結果を示した図である。
【図５】実施例において配列決定を行った糖鎖の単糖を用いた構成を示した図である。
【図６】更なる開裂を伴う質量分析（ＭＳ３）の結果を示した図である。
【図７】配列解析後に絞り込まれた糖鎖配列の候補を示した図である。
【図８】実施例で示した解析のフローを示した図である。
【図１】
2008175793.tif 000002
【図２】
2008175793.tif 000003
【図３】
2008175793.tif 000004
【図４】
2008175793.tif 000005
【図５】
2008175793.tif 000006
【図６】
2008175793.tif 000007
【図７】
2008175793.tif 000008
【図８】
2008175793.tif 000009