(57)【特許請求の範囲】
【請求項１】
（ａ）蛋白質の構成要素である２０個のアミノ酸それぞれについて、疎水性指標Ｈ、両親媒性指標Ａ、及び非電荷指標Ｃを割り当て、（ｂ）５〜９個の何れかの数の一連の窓を有するウィンドウＷを用いて、アミノ酸配列が決定された被判別蛋白質の該アミノ酸配列から抽出される前記ウィンドウＷと対応する全てのアミノ酸配列について、該アミノ酸配列の疎水性コア指標Icoreを、前記疎水性指標を用いて計算した一回平均疎水性［Ｈ］を介して算出し、（ｃ）得られたIcore値を縦軸とし、被判別蛋白質のアミノ酸鎖を横軸にしてプロットし、（ｄ）得られたグラフが一定の閾値を超える領域の一連のアミノ酸残基の数が１５以上となったとき、その領域を構成する一連のアミノ酸残基からなるセグメントを一次型膜貫通ヘリックスの候補とし、（ｅ）得られた各候補について二回平均疎水性＜Ｈ＞及び平均両親媒性＜Ａ＞を算出し、（ｆ）被判別蛋白質のアミノ酸鎖の長さ（アミノ酸の数）、上記二回平均疎水性、及び、平均両親媒性を用いた一定の判別式に従って一次型膜貫通ヘリックスの有無を判別し、（ｇ）一次型膜貫通ヘリックスが存在しない場合には被判別蛋白質を水溶性蛋白質と認定し、（ｈ）一次型膜貫通ヘリックスが存在する場合には、少なくとも、前記（ｄ）における閾値を変更すること以外は、前記（ｂ）〜（ｄ）と全く同様にしてセグメントを抽出し、該抽出したセグメントの集合から既に一次型膜貫通ヘリックスと認定されたセグメントを除去して二次型膜貫通ヘリックスの候補セグメントを得、（ｉ）次いで前記（ｅ）と同様にして各候補セグメントについて＜Ｈ＞’と＜Ａ＞’を算出し、（ｊ）前記二次型膜貫通ヘリックスの候補セグメントについてその中心が４５残基以下である場合と４５残基以上である場合に分類し、前記＜Ｈ＞’及び＜Ａ＞’を用いた新たな判別式によって二次型膜貫通ヘリックスの有無を判別すると共に、（ｋ）全く新しく、前記（ｂ）で抽出された全てのアミノ酸配列の単位に対して非電荷指標Ｃを割り当て、該非電荷指標Ｃが１であるアミノ酸残基が連続して１９残基以上存在するセグメントを抽出し、該セグメントについて前記＜Ｈ＞’及び＜Ａ＞’を用いた判別式によって二次型膜貫通ヘリックスの有無を判別する方法であって、前記＜Ｈ＞’及び＜Ａ＞’の演算の基本データとして（ａ）で各アミノ酸に割り当てた各指標及びその値を使用することを特徴とする、蛋白質の判別方法。
【請求項２】
前記(a)-(k)から成る判別によって膜蛋白質であると認定された場合に、得られた膜貫通ヘリックスの本数を、予め記録された、膜貫通ヘリックスの本数と膜蛋白質の機能との対応表に照し合わせることにより、判別に係る蛋白質が有する機能を更に推定することを特徴とする、請求項１に記載されたアミノ酸配列の決定された蛋白質の判別方法。
【請求項３】
前記(a)-(k)から成る判別によって膜蛋白質であると認定された場合に、膜貫通ヘリックスであると認定されたアミノ酸配列について、（１）前後７残基領域の３．６周期のスペクトル密度を最大エントロピー法で計算すると共にその領域の平均疎水性値を求め、（２）膜貫通ヘリックスの中心からＮ端及びＣ端側に向けて、それぞれ１乃至１０迄の数値ｊを下記一般式で表される膜環境フィルターに代入しながら（１）で計算したスペクトル強度を掛け合わせ、（３）上記（２）で得られた値と平均疎水性値の和の値が連続して設定された閾値を超えた時、その閾値の存在個所によって更に膜貫通領域を特定することを特徴とする、請求項１に記載されたアミノ酸配列の決定された蛋白質の判別方法。
0003722420.tif 000009
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、アミノ酸配列が決定された蛋白質について、その性質や機能あるいは構造を判別する方法に関し、特にコンピュータを用いて、迅速且つ高い精度で膜蛋白質であるか否か、膜蛋白質である場合には、膜蛋白質であるとの情報を与えるセグメントについての情報と共に、該膜蛋白質がいかなる機能を発揮するかを推定する蛋白質の判別方法に関する。
【０００２】
【従来技術】
膜蛋白質の二次構造の予測は、一般にアミノ酸配列に疎水性指標を割り当てたプロフィールから簡単な計算によって予測することができる[Kyte, j. and Doolittle, R. F., J. Mol. Biol., 157(1982)103-132.]。更に、ニュートラルネットワーク[Rost, B., Casadio, R., Fariseli, P. and Sander, C., Prot. Sci., 4(1995)521-533]や統計的手法[Jones, D. T., Taylor, W. R. and Thornton, J. M., Biochemistry, 33(1994)3038-3049.]による方法も報告されている。しかしながら、水溶性蛋白質に比べて一見精度よく見える膜蛋白質二次構造予測にも、次のようないくつかの問題点がある。
【０００３】
（１）ある種の蛋白質では、従来の方法では予測が困難な二次構造（膜貫通ヘリックス）が存在する。
（２）膜貫通ヘリックスの領域予測が十分でない。
（３）膜貫通ヘリックスに対して成績の良い方法でも、実際には水溶性蛋白質である多くのものを膜貫通ヘリックスとして誤認予測してしまう。
【０００４】
膜貫通ヘリックスは、基本的には疎水性相互作用によって脂質二重膜の疎水性領域に入り、安定となる。しかしながら、ある疎水性インデックスの値を閾値とし、それを越える一次構造領域を膜貫通ヘリックスとするような従来の単純な方法では、本来単なる水溶性蛋白質であり、膜貫通ヘリックスではないアミノ酸配列の解析からも、多くの膜貫通ヘリックスが誤って予測されるという欠点があった。一方、膜貫通ヘリックスにはしばしば極性残基が含まれており、特に細胞内側には正電荷を持つ残基が多いことなどが知られている[von Heijne, G., J. Mol. Biol., 225(1992)487-494.]。このことは膜貫通ヘリックス形成にも極性相互作用の存在（極性と疎水性のバランス）を考慮することが必要であることを示している。
【０００５】
一方、膜蛋白質ができるプロセスを考えてみると、アミノ酸配列にはいくつかの役割があるということが分かる。つまり、あるアミノ酸配列からなるセグメントを想定した場合、（イ）そのセグメント自体が膜貫通ヘリックスを作るべく設計された疎水性の高い配列である場合や、（ロ）膜貫通ヘリックスになったときに一定の機能を持つように設計されている場合などに分類することができる。後者の場合、膜蛋白質の機能に本質的に重要なアミノ酸配列は疎水性の低いアミノ酸残基を含む場合が多く、本来膜貫通ヘリックスを作ることに不利に働くが、近傍の疎水性の高いセグメントによって膜内に引き込まれ、膜貫通ヘリックスになると考えられる。
【０００６】
つまり、膜貫通ヘリックスには次の２つのタイプがあると考えられる。
１）アミノ酸配列が本来持つ高い疎水性によって、そのセグメントが自発的に膜貫通ヘリックスを形成するタイプ
２）膜蛋白質の機能的に重要なセグメントなど比較的親水的なセグメントが、上記1）のタイプの膜貫通ヘリックスによって膜内に引き込まれ、膜貫通ヘリックスを形成するタイプ
これらの膜貫通ヘリックスが周囲のヘリックス同士で相互作用することによって強固にパッキングされ、膜蛋白質の立体構造が安定化されるものと考えられる。以上のような膜蛋白質の立体構造形成についての考察から、本発明においては主として上記１）のタイプのものを一次型膜貫通ヘリックス、２）のものを二次型膜貫通ヘリックスと分類する。
【０００７】
そこで本発明者は、膜貫通ヘリックスをより高い精度で予測することについて鋭意検討した結果、Kyte-Doolittleの疎水性指標（以下、K-D疎水性指標とする）と、新しく定義した両親媒性指標および非電荷指標を用いて膜貫通候補セグメントを抽出・分類し、さらにタンパク質の大きさを考慮して予測した場合には、格段に高精度の予測が可能となることを見出し、本発明に到達した。
【０００８】
【発明が解決しようとする課題】
従って本発明の第１の目的は、アミノ酸の配列が決定された蛋白質が膜蛋白質であるか否かを高い精度で判別するための、蛋白質の判別方法を提供することにある。
本発明の第２の目的は、膜蛋白質について、該膜蛋白質の機能を推定することのできる蛋白質の判別方法を提供することにある。
更に本発明の第３の目的は、膜蛋白質における膜貫通ヘリックスの領域を高い精度で予測することのできる、蛋白質の判別方法を提供することにある。
【０００９】
【課題を解決するための手段】
本発明の上記の諸目的は、
（ａ）蛋白質の構成要素である２０個のアミノ酸それぞれについて、疎水性指標Ｈ、両親媒性指標Ａ、及び非電荷指標Ｃを割り当て、（ｂ）５〜９個の何れかの数の一連の窓を有するウィンドウＷを用いて、アミノ酸配列が決定された被判別蛋白質の該アミノ酸配列から抽出される前記ウィンドウＷと対応する全てのアミノ酸配列について、該アミノ酸配列の疎水性コア指標Icoreを、前記疎水性指標を用いて計算した一回平均疎水性［Ｈ］を介して算出し、（ｃ）得られたIcore値を縦軸とし、被判別蛋白質のアミノ酸鎖を横軸にしてプロットし、（ｄ）得られたグラフが一定の閾値を超える領域の一連のアミノ酸残基の数が１５以上となったとき、その領域を構成する一連のアミノ酸残基からなるセグメントを一次型膜貫通ヘリックスの候補とし、（ｅ）得られた各候補について二回平均疎水性＜Ｈ＞及び平均両親媒性＜Ａ＞を算出し、（ｆ）被判別蛋白質のアミノ酸鎖の長さ（アミノ酸の数）、上記二回平均疎水性、及び、平均両親媒性を用いた一定の判別式に従って一次型膜貫通ヘリックスの有無を判別し、（ｇ）一次型膜貫通ヘリックスが存在しない場合には被判別蛋白質を水溶性蛋白質と認定し、（ｈ）一次型膜貫通ヘリックスが存在する場合には、少なくとも、前記（ｄ）における閾値を変更すること以外は、前記（ｂ）〜（ｄ）と全く同様にしてセグメントを抽出し、該抽出したセグメントの集合から既に一次型膜貫通ヘリックスと認定されたセグメントを除去して二次型膜貫通ヘリックスの候補セグメントを得、（ｉ）次いで前記（ｅ）と同様にして各候補セグメントについて＜Ｈ＞’と＜Ａ＞’を算出し、（ｊ）前記二次型膜貫通ヘリックスの候補セグメントについてその中心が４５残基以下である場合と４５残基以上である場合に分類し、前記＜Ｈ＞’及び＜Ａ＞’を用いた新たな判別式によって二次型膜貫通ヘリックスの有無を判別すると共に、（ｋ）全く新しく、前記（ｂ）で抽出された全てのアミノ酸配列の単位に対して非電荷指標Ｃを割り当て、該非電荷指標Ｃが１であるアミノ酸残基が連続して１９残基以上存在するセグメントを抽出し、該セグメントについて前記＜Ｈ＞’及び＜Ａ＞’を用いた判別式によって二次型膜貫通ヘリックスの有無を判別する方法であって、前記＜Ｈ＞’及び＜Ａ＞’の演算の基本データとして（ａ）で各アミノ酸に割り当てた各指標及びその値を使用し、その結果膜蛋白質であると認定された場合に、更にその膜蛋白質の機能を推定し、及び/又は、その膜蛋白質における膜貫通ヘリックスの領域を予測することのできる、蛋白質の判別方法によって達成された。
【００１０】
【発明の実施の形態】
一般的に、一次型膜貫通ヘリックスは、水溶性蛋白質や膜蛋白質の水溶性ドメインには見られないような特徴を持っていて、特に疎水性が高い。しかしながら、疎水性だけではそれらを判断することができない。即ち、脂質からなる二重膜自体が、非常に疎水性の高い炭化水素鎖領域と表面に分布する極性基でできているために、膜貫通ヘリックスも、中心の疎水性部分だけで安定化しているのではなく、両端の極性部位による構造安定化がかなり重要であると考えられる。従って、中心部分に対する疎水性パラメーターである疎水性指標Ｈと共に、もう１つのパラメータとして膜貫通ヘリックスの両端にある極性残基の分布を考慮することが妥当であると考えられる。
【００１１】
また、一般に極性基があると共に炭化水素部分が大きい場合には、分子の両親媒性が高く、ミセルやラメラ構造などの超分子構造を作りやすい。これと同様に、側鎖自体の両親媒性が高いと膜貫通ヘリックスの両端の構造を安定化し易いのではないかと考えられる。本発明者は、このことを用いて、新たに、下記表１に示す両親媒性指標Ａ１，Ａ２を作成した。以後これを両親媒性指標Ａとする。尚、Ｃは後述する非電荷指標であり、（）内のアルファベットは、三文字表記のアミノ酸を一文字表記する場合の記号である。
【００１２】
【表１】
0003722420.tif 000002
【００１３】
以下、本発明を具体的手順に従って説明する。
本発明においては、先ず判別しようとする蛋白質を構成する各アミノ酸に、表１によって予め設定されているＫ−Ｄ疎水性指標Ｈを割り当てる。次に、５〜９箇から選択されるアミノ酸残基の数、好ましくは７箇の連続するアミノ酸残基に当てはめることのできるウィンドウＷを用いて、下記（１）式によって一回平均疎水性［Ｈ］を計算する。
0003722420.tif 000003
但し、ｉはウィンドウの中心のアミノ酸の座標である。
上記の計算は、例えば７残基について行う場合（Ｗ＝７、Ｗ／２＝３とする）には蛋白質を構成するアミノ酸鎖の左端から７残基ウィンドウを１残基ずつずらしながら、５残基（Ｗ＝５、Ｗ／２＝２とする）について行う場合には５残基ウィンドウを１残基ずつずらしながら、全ての単位について計算する。
【００１４】
例えば、下記（５）式によって、疎水性コア指標であるＩｃｏｒｅ（ｉ）を計算する。
Ｉｃｏｒｅ（ｉ）＝Positive（［Ｈｍａｘ（ｉ）］+Ｓ（ｉ）・［Ｈｍｉｎ（ｉ）］）  （５）
但し、［Ｈｍａｘ（ｉ）］、［Ｈｍｉｎ（ｉ）］及びＳ（ｉ）は下記６〜８式で表される。
［Ｈｍａｘ（ｉ）］＝Max（［Ｈ（ｊ）］：ｊ＝ｉ−ｎ…，０，…，ｉ＋ｎ）  （６）
［Ｈｍｉｎ（ｉ）］＝Min（［Ｈ（ｊ）］：ｊ＝ｉ−ｎ，…，０，…ｉ＋ｎ）  （７）
Ｓ（ｉ）＝Ｓｉｇｎ（［Ｈｍａｘ（ｉ）］・［Ｈｍｉｎ（ｉ）］）                （８）
【００１５】
上記のＩｃｏｒｅ（ｉ）は正負の間を変動する関数の、正の部分（Positive）を表す。尚、Positive関数は、負の値の場合は０とする。この正の値が特定の閾値ｋより大きくなる部分を抽出し、その部分の連続するアミノ数の数がＬより大であるものを、膜蛋白質か否かを判別するためのデータを保有する第一次セグメントの候補Ｓ（ｉ）とし、この候補セグメントの集合をＳ１とする。但しＬは、１５〜２５から選択された何れかのアミノ酸残基の数、即ち長さであり、好ましくは１８〜２０から選択される整数である。尚、Ｌを１９とした場合には、前記ｋとして、１．１５を設定することが好ましい。また、Ｍａｘ及びＭｉｎは、領域の中から最大値及び最小値を選び出す関数であり、Ｓｉｇｎは（）内が正のときＳ（ｉ）が１となり、（）内が負のときＳ（ｉ）が−１となるように符号を返す関数である。
【００１６】
上記の如くして抽出された候補Ｓ（ｉ）の各セグメントについて、下記の式により、二回平均疎水性＜Ｈ＞と平均両親媒性＜Ａ＞を計算する。下記式は、ウィンドウＷを７とした場合について示したものである。尚、Ａは表１に示された両親媒性指標であり、[Ｈ]は前記（１）式で求めたものである。但し、ｉｃはセグメントの中心、ｍはセグメント長の１／２である。
0003722420.tif 000004
0003722420.tif 000005
0003722420.tif 000006
【００１７】
次に、被判別蛋白質の長さＮが約１７５（ここでは１７５と設定）より長い場合と短い場合に分けて、Ｓ（ｉ）の集合Ｓ１に含まれる全てのセグメントについて次式の判別式を実行する。但し、通常、Ｎは９０〜２５０の整数である。
（ア）：Ｎ≧１７５の場合：＜Ａ＞≧-４＜Ｈ＞+４．８及び＜Ａ＞≧-１＜Ｈ＞+１．５の両判別式を同時に満足するセグメントは一次型膜貫通ヘリックスである。
（イ）：Ｎ＜１７５の場合：判別式＜Ａ＞≧-１．２５＜Ｈ＞+１．６２５を満足するセグメントは一次型膜貫通ヘリックスである。
尚、上記判別式は、一般式＜Ａ＞≧-ａ＜Ｈ＞+ｂで表され、このａ及びｂは、既知の多数の蛋白質に対する＜Ａ＞及び＜Ｈ＞を計算し、＜Ａ＞及び＜Ｈ＞を夫々Ｘ軸又はＹ軸としてグラフ化した散布図に基づいて選択される。従って、上記判別式における定数ａ、ｂは、夫々上記の値±１０％の範囲で異なった値に設定される場合がある。
【００１８】
上記判別式を満たすセグメントが前記候補セグメントの中に１つも存在しない場合には、被判別蛋白質は水溶性蛋白質であって、膜蛋白質ではないと判定され、判定操作（プログラム）は終了する。従って、この蛋白質については以降の演算はされず、別の被判別蛋白質について、初めから本判定操作（プログラム）が実行される。
一方、上記判別式を満たすセグメントが抽出されたとき、そのセグメントを一次型膜貫通ヘリックスであると判定する。
尚、上記の判別式は、既知の蛋白質について、平均両親媒性＜Ａ＞の二回平均疎水性＜Ｈ＞依存性を表す散布図を作成することによって求めることができる。
【００１９】
次に、類似の方法によって二次型膜貫通ヘリックスの存在の有無を判定する。この場合の候補セグメントとして、少なくともＳ１を求めたときの閾値ｋの値を変えること以外は、Ｓ１の場合と同様にして候補セグメントを抽出し、この候補セグメントの集合の中から、前記一次型膜貫通ヘリックスを除去した残りの全セグメント（その集合をＳ１’とする）、及び、新たに、非電荷指標Ｃが１であるアミノ酸が、１５〜２５から選択された何れかの長さＬ’以上に連なった全てのセグメント（その集合をＳ２とする）の集合を含める。Ｓ１’を求めるとき、ウィンドウＷや閾値ｋより大きくなる部分の連続するアミノ酸の数Ｌ等も、ｋの値と共に変更しても良い。
【００２０】
次に、一次型膜貫通ヘリックスの判定の場合と同様にして、次式によって二次型膜貫通ヘリックスの有無を判定する。
Ｓ１’の集合に含まれる各セグメントに対する判別式は下記の通りである。
Ｉｃｏｒｅで決まる領域の中心が４５残基以下の場合：＜Ａ＞≧-２＜Ｈ＞+１．６
Ｉｃｏｒｅで決まる領域の中心がＮ端から４５残基以上、且つ前記中心から前後６５残基以内に一次型膜貫通ヘリックスの中心が存在する場合：＜Ａ＞≧-２．７１＜Ｈ＞+２．５７
Ｓ２の集合に含まれる各セグメントについては、Ｓ１の場合と全く同様にして＜Ｈ＞及び＜Ａ＞を求める。次に各候補セグメントについて、中心から前後５０残基以内に一次型膜貫通ヘリックス又は上記判別式によって判別されたセグメントが存在することを条件として＜Ａ＞≧-５＜Ｈ＞+１．９の判別式を満たすものを抽出する。
但し、これらの判別式における勾配や定数も±１０％の範囲で変更する場合もある。
Ｓ１’及びＳ２について上記判別式を満足するセグメントがない場合には二次膜貫通ヘリックスは存在しないことになり、判別操作（プログラム）は終了する。
上記判別式を満足するセグメントが存在する場合には、そのセグメントを二次型膜貫通ヘリックスであると判定する。
【００２１】
以上の判定によって得られた膜貫通ヘリックスの末端の最終的な領域を決定するために、本発明においては、最大エントロピー法によって周期のスペクトル密度を求め、それを下記ｆ（ｊ）の式で表される膜環境フィルタと組み合わせて膜貫通ヘリックスの末端を決定する。
0003722420.tif 000007
但し、ｊはコアの中心から左右へ伸長する値である。
膜貫通ヘリックスの末端は、以下の操作によって決定する。
１．膜貫通ヘリックスのアミノ酸配列について、その前後７残基領域の３．６周期のスペクトル密度を最大エントロピー法で計算し、同時にその領域の平均疎水性値を求める。
２．膜貫通ヘリックスコアの中心からＮ端、Ｃ端側に向けて１，２，…１０までの数値ｊを膜環境フィルターに代入しながら、（１）で計算したスペクトル強度を掛け合わせる。
３．（２）で得られた値と平均疎水性値を足しあわせてその値が連続してある閾値を超えれば、その領域を膜貫通ヘリックス領域とする。
尚、前記最大エントロピー法を用いる計算の詳細は、Biophysical Chemistry, 34(1989)の第７０−７１頁、「２．２  Ｓｐｅｃｔｒａｌ  ａｎａｌｙｓｉｓ」の項に記載されている。
【００２２】
被判別蛋白質のアミノ酸配列について以上のデータ入力と演算を行わせ、不等式の判別式を満足する解がない場合には、判別操作（プログラム）はそこで終了する。本発明においては、各判別段階によって得られた全ての膜貫通ヘリックスの本数とそれらの領域についての情報を、モニター及び／又はプリンターによって出力する。勿論、各判別段階における、それらの情報を、個別に出力するようにすることも、全体情報と併せて各個別データを付加的に出力することもできる。
【００２３】
上記の判別を実施する為のコンピュータプログラムのフローチャートは、図１に示される通りである。本発明のコンピュータプログラムは、Ｃ言語等を用いて記載することができる。以下に更に詳述する。
【００２４】
本発明のプログラムを格納した記録手段を有するコンピュータを起動し、被判別蛋白質を形成するＮ個のアミノ酸配列を入力すると、各アミノ酸にＮ個の連続した整数が割り当てられ（ＳＴＰ１００）、入力されたアミノ酸残基の数がＮ個であることが認識される。次に、入力された全てのアミノ酸残基に、予めメモリ部（記録手段）に保存されている前記表１のデータのうち、該当する疎水性指標Ｈの値を割り当てる（ＳＴＰ１０１）。
【００２５】
所定のウィンドウ（例えば７つの連続した窓を有するウィンドウ）を、端から１残基ずつずらしながら、各ウィンドウに対応するアミノ酸列を抽出し、抽出された全てのアミノ酸列に対して、前記（１）式に従って一回平均疎水性［Ｈ］を計算する（ＳＴＰ１０２）。次に前記（５）式に従って、上記の全てのアミノ酸列に対してIcoreを計算しこれを被判別蛋白質のアミノ酸配列の順に対応させて整列させる（Icoreのグラフ化：ＳＴＰ１０３）。
【００２６】
Icoreがｋ（例えば１．１５）より大きくなる領域を抽出し、その領域のうちＬ（例えば１９）より大きい数のアミノ酸残基が連続する部分のアミノ酸配列を候補セグメントとして抽出し、その集合をＳ１とする（ＳＴＰ１０４）。次いで、この各候補セグメントについて、表１のデータを用い、前記した（２）〜（４）式に従って二回平均疎水性＜Ｈ＞及び平均両親媒性＜Ａ＞を計算する（ＳＴＰ１０５）。
【００２７】
被判別蛋白質を形成するアミノ酸残基の数Ｎが１７５より大きい場合には、ＳＴＰ１０４で得られた集合Ｓ１に含まれる全ての候補セグメントについて、＜Ａ＞≧-４＜Ｈ＞+４．８及び＜Ａ＞≧-１＜Ｈ＞+１．５を同時に満たすか否か判別し、これらの判別式を満たすセグメントが抽出された場合にそれを一次型膜貫通ヘリックスと認定する。
一次型膜貫通ヘリックスが存在しない場合には、プログラムは終了し、被判別蛋白質が水溶性蛋白質であるとの結果を、モニター及び/又はプリンターに出力する（ＳＴＰ１０６）。
【００２８】
一次型膜貫通ヘリックスが存在した場合には、次に二次型膜貫通ヘリックスの有無の判断に移る。
少なくとも一次形膜貫通ヘリックスの判別操作で使用した閾値ｋをｋ’に変更し、必要に応じて、ウィンドウＷや、候補セグメントを抽出するための閾値ｋ’を超える領域として設定する一連のアミノ酸残基の数を更に変更すること以外は、Ｓ１の場合と同様にして候補セグメントを抽出し、この中から一次型膜貫通ヘリックスに該当するセグメントを除去し、この集合Ｓ１’を二次型膜貫通ヘリックスの候補セグメントの部分集合と認定する（ＳＴＰ１０７）。次に、この集合Ｓ１’に含まれる各セグメントについて、Ｓ１の場合と同様にして＜Ｈ＞及び＜Ａ＞を求める（ＳＴＰ１０８）。
【００２９】
セグメントの中心が４５残基以下であるものに対しては判別式＜Ａ＞≧-２＜Ｈ＞+１．６を満足するものを二次型膜貫通ヘリックスと認定し、４５残基以上である場合には、更に、該セグメントの前後６５残基以内に一次型膜貫通ヘリックスが存在することを条件に、判別式＜Ａ＞≧-２．７１＜Ｈ＞+２．５７を満足するものを抽出し、これも二次型膜貫通ヘリックスと認定する（ＳＴＰ１０９）。
【００３０】
次に、全く新しく、被蛋白質の全アミノ酸残基に、表１に基いてＣの値を割り当てる。その結果Ｃ＝１のアミノ酸残基が連続して１９残基を超える場合、その領域を二次型膜貫通ヘリックスの候補セグメントとし、その集合をＳ２とする。該候補セグメントのうち、その中心から前後５０残基以内に一次型膜貫通ヘリックス又は【００２７】で認定した二次型膜貫通ヘリックスが存在するものに対しては、判別式＜Ａ＞≧-５＜Ｈ＞+１．９を満足するものを二次型膜貫通ヘリックスと認定する（ＳＴＰ１１０）。
これらの全ての演算が完了すると、一次型膜貫通ヘリックスの本数と領域、及び二次型膜貫通ヘリックスの本数と領域を出力してプログラムは、終了する。
【００３１】
膜貫通ヘリックスの本数と膜蛋白質の機能とは一定の関係を有するので、その対応関係を表１の場合と同様に、予めコンピュータの記憶部に記憶させておくと、被判別蛋白質から抽出された全膜貫通ヘリックスの本数を該データと比較対照して、被判別蛋白質の機能をモニター及び/又はプリンターに出力することができる。
また、得られた膜貫通ヘリックスのアミノ酸配列について、（１）該アミノ酸配列の前後７残基領域の３．６周期のスペクトル密度を最大エントロピー法で計算すると共にその領域の平均疎水性値を求め、（２）膜貫通ヘリックスの中心からＮ端及びＣ端側に向けて、それぞれ１乃至１０迄の数値ｊを下記一般式で表される膜環境フィルターに代入しながら、（１）で計算したスペクトル強度を掛け合わせ、（３）上記（２）で得られた値と平均疎水性値の和の値が連続して設定された閾値を超えた時、その閾値の存在個所によって特定された領域を、膜貫通ヘリックスの領域として抽出し、モニター及び/又はプリンターに出力することができる。
0003722420.tif 000008
【００３２】
【発明の効果】
本発明によれば、アミノ酸の一次配列が解明された蛋白質について、９７％以上の正答率で、迅速に、被判別蛋白質が膜蛋白質であるか、水溶性蛋白質であるかを判定することができる。
【図面の簡単な説明】
【図１】本発明のコンピュータプログラムのフローチャートの全体図の１例である。
【図２】図１におけるＳＴＰ１０９のフローチャートである。
【図３】図１におけるＳＴＰ１１０のフローチャートである。
【図１】
0003722420.tif 000010
【図２】
0003722420.tif 000011
【図３】
0003722420.tif 000012