(57)【特許請求の範囲】
【請求項１】
アミノ酸配列の決定された水溶性蛋白質が、ダンベル型蛋白質であるか否かを判別する方法であって、該方法が、（ａ）蛋白質の構成要素である20種のアミノ酸それぞれについて電荷指標Cを割り当てて該蛋白質の総電荷指標の絶対値｜CT｜を算出し、得られた|CT|が一定の閾値を超える蛋白質については、（ｂ）蛋白質全体の電荷密度の絶対値|D|を算出し、該|D|が一定の閾値を超える蛋白質について、（ｃ）蛋白質のアミノ酸配列を半分に分割した後、分割されたそれぞれの総電荷指標CN及びCCを算出し、CN×CCが一定の閾値以上となる蛋白質をダンベル型蛋白質の候補とし、（ｄ）該蛋白質を構成する各アミノ酸に疎水性指標Hと両親媒性指標Aを割り当て、（ｅ）次いで5〜9個の何れかの数の一連の窓を有するウィンドウWを用いて被判別蛋白質のアミノ酸配列から抽出される前記ウィンドウWと対応する全てのセグメントについて一回平均疎水性[H]及び一回平均両親媒性[A]を計算し、（ｆ）前記アミノ酸の連続10残基以上が、前記一回平均疎水性[H]として-0.25以下をとる領域があるか否かを判別し、該領域があればその領域を前記アミノ酸配列のCore領域として抽出し、（ｇ）得られたCore領域について二回平均疎水性<H>及び二回平均両親媒性<A>を算出し、これらがそれぞれ一定の閾値を超える場合に、（ｈ）Core領域の両端残基を除いた領域の[H]の最大値Maxと最小値Minを算出し、Max−Minが２．０以下となる場合に、（ｉ）中央ヘリックス候補セグメントを、前記Core領域外両端の一回平均疎水性[H]が-0.30以下になる直前のアミノ酸残基を中央へリックスの末端とすることによって抽出し、（ｊ）得られた候補セグメントを除いたN端側のアミノ酸残基数NN（アミノ酸の長さ）とC端側のアミノ酸残基数NC（アミノ酸の長さ）を算出し、該NNとNCが、    NC≦（3/2）×NN    かつ    NC≧（2/3）×NNの関係を有する場合に、（ｋ）候補セグメントの総電荷指標CHと候補セグメントを除いたN端側の総電荷指標CNDとC端側の総電荷指標CCDを算出し、これらの総電荷指標が全て正であるか全て負である場合にその蛋白質をダンベル型水溶性蛋白質と判別することを特徴とするダンベル型水溶性蛋白質の判別方法。
【請求項２】
蛋白質の構成要素である２０種のアミノ酸のそれぞれに対して割り当てられた疎水性指標Ｈ、両親媒性指標Ａ，及び電荷指標Ｃを用い、入力されたアミノ酸配列について、一回平均疎水性[H]、二回平均疎水性<H>、二回平均両親媒性<A>、総電荷指標|CT|及び電荷密度|D|を計算し、これらの値を利用して、アミノ酸配列の決定された蛋白質がダンベル型蛋白質であるか否かを判別するコンピュータプログラムであって、被判別蛋白質の全てのアミノ酸配列を入力し、判別プログラムを実行させると、（ａ）予め記憶部に記録された、前記２０種のアミノ酸に設定された前記疎水性指標Ｈ、両親媒性指標Ａ及び電荷指標Ｃのデータのうち、電荷指標Ｃのデータを参照して、入力した全てのアミノ酸残基のそれぞれに前記電荷指標Ｃを割り当て、該蛋白質の総電荷指標の絶対値｜CT｜を算出し、得られた|CT|が２０以上であるか否かを判別し、２０未満であればプログラムを終了し、（ｂ）２０以上となる蛋白質については蛋白質全体の電荷密度の絶対値|D|を算出し、該|D|が０．１４以上であるかを判別する；０．１４未満であればプログラムを終了し、（ｃ）０．１４以上を満足する蛋白質について、全アミノ酸配列を半分に分割した後、分割されたそれぞれの総電荷指標CN及びCCを算出し、CN×CCが１００以上である蛋白質をダンベル型蛋白質の候補とし、該蛋白質を構成する各アミノ酸に疎水性指標Hと両親媒性指標Aを割り当て、（ｄ）次いで5〜9個の何れかの数の一連の窓を有するウィンドウWを用いて被判別蛋白質のアミノ酸配列から抽出される前記ウィンドウWと対応する全てのセグメントについて一回平均疎水性[H]及び一回平均両親媒性[A]を計算する；次に（ｅ）前記一回平均疎水性[H]が連続して１０残基以上−０．２５以下となる領域を探索し、該当する領域が無ければプログラムを終了し、該当する領域がある場合には、その領域をＣｏｒｅ領域として定義する；（ｆ）得られたCore領域について二回平均疎水性<H>及び二回平均両親媒性<A>を算出し、<H>の値が−１．１５以下であると共に、<A>の値が０．９０以上となるか否かを判別し、この条件が満たされない場合にはプログラムを終了する；（ｇ）前記（ｆ）の条件が満たされた場合に、前記Core領域の両端３残基を除いた領域の一回平均疎水性[H]の最大値Maxと最小値Minを算出し、Max−Minが２．０以下となるかを判別し、２．０を超える場合にはプログラムを終了する；（ｈ）前記（ｇ）の条件が満たされた場合、Core領域外の両端に向けて一回平均疎水性[H]を探索し、その値が−０．３０以下になる直前のアミノ酸残基を、それぞれ中央ヘリックス候補セグメントの両末端とする；次に（ｉ）得られた候補セグメントを除いたN端側のアミノ酸残基数NN（アミノ酸の長さ）とC端側のアミノ酸残基数NC（アミノ酸の長さ）を算出し、NC≦（３／２）NN且つNC≧（２／３）NNの判別条件を実行し、これらの条件が満たされない場合にはプログラムを終了する；（ｊ）上記判別条件を満たすセグメントが存在した場合、該セグメントを中央ヘリックスの候補セグメントとし、（ｋ）該候補セグメントの総電荷指標CHと候補セグメントを除いたN端側の総電荷指標CNDとC端側の総電荷指標CCDを算出し、これらの総電荷指標が全て正であるか全て負である場合にその蛋白質をダンベル型水溶性蛋白質と判���し、その中央ヘリックスを出力して、全プログラムを終了することを特徴とするダンベル型水溶性蛋白質を判別するためのコンピュータプログラム。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、アミノ酸配列が決定された水溶性蛋白質のコンピュータを用いた判別方法に関し、特に、コンピュータを用いて迅速且つ高い精度でダンベル型水溶性蛋白質であるか否かを判別する方法、及びダンベル型水溶性蛋白質の中央ヘリックスについての情報を出力する方法、並びにそれを実施する為のコンピュータプログラムに関する。
【０００２】
【従来技術】
蛋白質の立体構造の分類方法及び構造予測方法は既に数多く知られているが、水溶性蛋白質の二次構造予測の精度はあまり高くないというのが現状である。その原因の一つとして、二次構造の分類がうまくなされていないということが挙げられる。即ち、統計的・経験的手法による水溶性蛋白質の二次構造予測の正当率は7割を越えていないが、これは、二次構造が同じであれば、それだけで一様に分類するということが一因と考えられる。
【０００３】
例えば、水溶性蛋白質に含まれるヘリックスでの結合には、水相にむき出しになっているタイプのものもあれば、内部コアで形成されているタイプのものや極性と非極性の界面付近でうまく安定しているタイプのものもある。つまり、構造的には同じヘリックスでも環境が異なる場合には蛋白質の立体構造が同一であるとすることはできないのであるが、統計的・経験的手法ではこれらを同じものとして分類せざるを得ないので予測精度が下がる。従って、このようなヘリックスの環境が異なるものを、異なるものとして分類することが予測精度の向上につながると推定される。
【０００４】
本発明が対象としている、水相にむき出しになったヘリックスをもつダンベル型水溶性蛋白質についても状況は同様である。尚、ここでダンベル型蛋白質とは、二つの水溶性ドメインが、一本の長い中央ヘリックスによって結ばれた構造をしている蛋白質である。
一般に、水溶性蛋白質の大部分は球状の形態をとることによって安定化しているが、ダンベル型蛋白質の場合には、二つの水溶性ドメインが同程度の大きな電荷を帯びて反発し合うことにより安定化していると考えられる。また、水溶性ドメインを結ぶ中央ヘリックスは、両親媒性残基が多く出現しそれらがヘリックスの安定化に寄与していると考えられる。
【０００５】
そこで本発明者は、ダンベル型水溶性蛋白質と中央ヘリックスを高い精度で予測することについて鋭意検討した結果、蛋白質の大きさとその蛋白質を構成する各アミノ酸固有の電荷指標を用いて蛋白質を選別し、さらにKyte-Doolittleの疎水性指標（以下、K-D疎水性指標とする）と、新しく各アミノ酸に対して定義した両親媒性指標を用いて中央ヘリックス領域を予測したところ、極めて高精度の予測が可能となることを見出し、本発明に到達した。
【０００６】
【発明が解決しようとする課題】
従って本発明の第１の目的は、アミノ酸配列が決定された水溶性蛋白質がダンベル型水溶性蛋白質であるか否かを迅速且つ高い精度で判別する方法を提供することにある。
本発明の第２の目的は、ダンベル型水溶性蛋白質の中央ヘリックスを、高い精度で判別するための方法を提供することにある。
更に本発明の第３の目的は、水溶性蛋白質の一次構造情報を用いて、その蛋白質がダンベル型水溶性蛋白質であるか否か、ダンベル型水溶性蛋白質である場合には更にその中央ヘリックスを、高い精度で判別するためのコンピュータプログラムを提供することにある。
【０００７】
【課題を解決するための手段】
本発明の上記の諸目的は、アミノ酸配列の決定された水溶性蛋白質が、ダンベル型蛋白質であるか否かを判別する方法であって、該方法が、（ａ）蛋白質の構成要素である20種のアミノ酸それぞれについて電荷指標Cを割り当てて該蛋白質の総電荷指標の絶対値｜CT｜を算出し、得られた|CT|が一定の閾値を超える蛋白質については、（ｂ）蛋白質全体の電荷密度の絶対値|D|を算出し、該|D|が一定の閾値を超える蛋白質について、（ｃ）蛋白質のアミノ酸配列を半分に分割した後、分割されたそれぞれの総電荷指標CN及びCCを算出し、CN×CCが一定の閾値以上となる蛋白質をダンベル型蛋白質の候補とし、（ｄ）該蛋白質を構成する各アミノ酸に疎水性指標Hと両親媒性指標Aを割り当て、（ｅ）次いで5〜9個の何れかの数の一連の窓を有するウィンドウWを用いて被判別蛋白質のアミノ酸配列から抽出される前記ウィンドウWと対応する全てのセグメントについて一回平均疎水性[H]及び一回平均両親媒性[A]を計算し、（ｆ）前記アミノ酸の連続10残基以上が、前記一回平均疎水性[H]として-0.25以下をとる領域があるか否かを判別し、該領域があればその領域を前記アミノ酸配列のCore領域として抽出し、（ｇ）得られたCore領域について二回平均疎水性<H>及び二回平均両親媒性<A>を算出し、これらがそれぞれ一定の閾値を超える場合に、（ｈ）Core領域の両端残基を除いた領域の[H]の最大値Maxと最小値Minを算出し、Max−Minが２．０以下となる場合に、（ｉ）中央ヘリックス候補セグメントを、前記Core領域外両端の一回平均疎水性[H]が-0.30以下になる直前のアミノ酸残基を中央へリックスの末端とすることによって抽出し、（ｊ）得られた候補セグメントを除いたN端側のアミノ酸残基数NN（アミノ酸の長さ）とC端側のアミノ酸残基数NC（アミノ酸の長さ）を算出し、該NNとNCが、    NC≦（3/2）×NN    かつ    NC≧（2/3）×NNの関係を有する場合に、（ｋ）候補セグメントの総電荷指標CHと候補セグメントを除いたN端側の総電荷指標CNDとC端側の総電荷指標CCDを算出し、これらの総電荷指標が全て正であるか全て負である場合にその蛋白質をダンベル型水溶性蛋白質と判別することを特徴とするダンベル型水溶性蛋白質の判別方法、及び、その方法を実施するコンピュータプログラムによって達成された。
【０００８】
【発明の実施の形態】
一般に、水溶性蛋白質の大部分は球状の形態をとることによって安定化しているが、ダンベル型蛋白質は球状の形態をとっていない。一方、ダンベル型蛋白質は、蛋白質の総電荷指標が大きく、電荷密度が高いという特徴と共に二つのドメインが同程度の大きな電荷を帯びているという特徴がある。従って、ダンベル型蛋白質の構造安定化には、蛋白質の電荷を考慮することが重要である。
【０００９】
また、極性基があると共に炭化水素部分が大きい場合には、一般に分子の両親媒性が高く、これによって極性と非極性の界面が安定化される。これと同様に、側鎖自体の両親媒性が高いと水に露出した中央ヘリックスの構造を安定化し易いと考えられる。本発明者は、このことを用いて、新たに下記表1に示す両親媒性指標Aを作成した。なお、Hは疎水性指標、Cは電荷指標であり、（）内のアルファベットは、三文字表記のアミノ酸を一文字表記する場合の記号である。
【００１０】
【表１】
0004021153.tif 000002
【００１１】
以下、本発明を具体的手順に従って説明する。
本発明においては、ダンベル型蛋白質を次のように定義する（図１参照）。
（１）PDBに記述されているヘリックス領域のアミノ酸残基[HN]が２０残基以上である。
（２）（１）のヘリックス両端のアミノ酸残基[NN][CN]が５０残基以上である。
（３）（１）のヘリックスの中央残基を通り、ヘリックス末端２残基のCα原子を結ぶベクトルを法線ベクトルとする平面を考えたとき、ヘリックスN(C)端側の全Cα原子がN(C)端側の空間に存在する割合[PN]（[PC]）を計算し、それが９９％以上である。
（４）蛋白質のサイズが大きいとドメインがいくつかあることが考えられるので、総残基数が５００残基以内である。
【００１２】
従って、本発明においては、蛋白質を構成するアミノ酸残基の数（アミノ酸鎖の長さ）Nが１００以上５００以下の蛋白質を判別の対象とする。尚、１００残基以上としたのは、ダンベル型蛋白質の定義が１２０残基以上であること、及び、１００残基未満の場合にはペプチドや蛋白質断片が多く含まれるからである。
本発明においては、先ず判別しようとする蛋白質を構成する各アミノ酸に、表1によって予め設定されている電荷指標Cを割り当て、下記（1）式によって総電荷指標の絶対値|CT|を計算する。
0004021153.tif 000003
（１）
但し、Nはアミノ酸残基の数（アミノ酸鎖の長さ）である。
【００１３】
この総電荷指標の絶対値|CT|が２０以上であるか否かを判別し、|CT|が２０以上である場合には下記（2）式によって蛋白質全体の電荷密度の絶対値|D|を計算する。
0004021153.tif 000004
（２）
尚、閾値を２０とした理由は、ダンベル型蛋白質は総電荷指標の絶対値が２０以上の大きな電荷を帯びており、総電荷指標の絶対値が２０以上の蛋白質の全体に占める割合は１．５％程度である。
【００１４】
次に、この蛋白質全体の電荷密度の絶対値|D|が０．１４以上であるか否かを判別し、上記|D|が０．１４以上である場合には、対象としているアミノ酸配列を半分に分割し、それぞれの総電荷指標CN、CCを、下記（3-1）及び（3-2）式によって計算する。
0004021153.tif 000005
（３−１）
0004021153.tif 000006
（３−２）
但し、C(i)はi番目のアミノ酸残基の電荷、Nは全アミノ酸残基の数（アミノ酸鎖の長さ）である。Nが奇数のときはN/2の小数を四捨五入する。尚、既知のダンベル型蛋白質の解析から、負の単位残基当たりの電荷は０．１４以上と大きく、０．１４以上の蛋白質の全体に占める割合は０．４％程度であることが分かっている。
【００１５】
上記のCN、CCよりCN×CCを計算し、その値が１００以上であるか否かを判別し、CN×CCが１００以上である場合には、対象としている蛋白質を構成する各アミノ酸に、表1によって予め設定されている疎水性指標H及び両親媒性指標Aを割り当てる。尚、既知のダンベル型蛋白質の解析から、負の単位残基当たりの電荷は１０以上と大きく、１０以上の蛋白質の全体に占める割合は０．３％程度であることが分かっている。
【００１６】
次に、5〜9箇から選択されるアミノ酸残基の数、好ましくは7箇の連続するアミノ酸残基に当てはめることのできるウィンドウWを用い、下記（4-1）及び（4-2）式によって一回平均疎水性 [H] 及び一回平均両親媒性  [A] を計算する。
0004021153.tif 000007
（４−１）
0004021153.tif 000008
（４−２）
但し、iはウィンドウの中心のアミノ酸の座標である。
【００１７】
上記の計算は、例えば7残基について行う場合（W=7、W/2=3とする）には蛋白質を構成するアミノ酸鎖の左端から7残基ウィンドウを1残基ずつずらしながら、5残基（W=5、W/2=2とする）について行う場合には5残基ウィンドウを1残基ずつずらしながら、全ての単位について計算する。次に、連続10残基以上が、一回平均疎水性 [H] として-0.25以下をとる領域（Core領域とする）があるか否かを判別し、[H]が連続１０残基以上で-0.25以下であるCore領域があればその領域を抽出する。尚、連続１０残基以上とした理由は、中央ヘリックスの疎水性の低い領域が１０残基より長いからである。
【００１８】
上記の如くして抽出されたCore領域に対して、下記（5-1）、（5-2）式によって、二回平均疎水性<H> 及び 二回平均両親媒性 <A> を計算する。
0004021153.tif 000009
（５−１）
但し、[H] は一回平均疎水性、S及びEはCore領域の開始残基と終始残基である。
また、NcoreはCore領域の残基数（E-S+1）である。
0004021153.tif 000010
（５−２）
但し、[A] は一回平均両親媒性、S及びEはCore領域の開始残基と終始残基である。
また、NcoreはCore領域の残基数（E-S+1）である。
【００１９】
次に、このCore領域の二回平均疎水性 <H> の値が-1.15以下であると共に、二回平均両親媒性 <A> の値が0.90以上であるか否かを判別する。尚、これらの閾値は既知のダンベル型蛋白質の解析に基づいて設定したものである。上記の２条件が満たされた場合には、抽出されたCore領域に対して、両端3残基を除いた領域の一回平均疎水性 [H] に対して最大値Maxと最小値Minを求める。次いで、Max-Minを計算し、その値が2.0以下となるか否かを判断する。
【００２０】
Max-Minの値が２．０以下であると判断されたCore領域に対して、その領域外両端に一回平均疎水性 [H] の値を探索していき、-0.30以下になる直前のアミノ酸残基をそれぞれ中央ヘリックス候補セグメントの末端とする。
上記の中央ヘリックス候補セグメントに対して、そのセグメントを除いたN端側の総残基数NNとC端側の総残基数NCを計算し、次式の判別式を実行する。尚、Ｎ端はアミノ基末端、Ｃ端はカルボキシル基末端を意味する。
NC ≦（3/2）×NN    かつ    N C ≧（2/3）×NN
上記判別式を満たすセグメントを抽出する。
【００２１】
上記の如くして抽出された中央ヘリックス候補セグメントに対して、下記（6-1）式によって総電荷指標CHを計算する。また、セグメントを除いたN端側の総電荷指標CNDとC端側の総電荷指標CCDを下記（6-2）、（6-3）式によって計算する。
0004021153.tif 000011
（６−１）
0004021153.tif 000012
（６−２）
0004021153.tif 000013
（６−３）
但し、C(i)はアミノ酸残基iの電荷指標、S及びEは中央ヘリックス候補セグメントの開始残基と終始残基である。Nはアミノ酸鎖の長さ（アミノ酸の長さ）である。
【００２２】
次に、次式の判別式を実行し、判別式を満たすセグメントを抽出する。
CH ＞ 0    かつ    CND ＞ 0    かつ    CCD ＞ 0  あるいは
CH ＜ 0    かつ    CND ＜ 0    かつ    CCD ＜ 0
上記の如くして抽出された中央ヘリックス候補セグメントを中央ヘリックス領域とし、その領域を含む蛋白質をダンベル型水溶性蛋白質と判別する。
【００２３】
上記の判別を実施する為のコンピュータプログラムのフローチャートは図2〜4に示される通りである。本発明のコンピュータプログラムはC言語などを用いて記載することができる。以下に更に詳述する。
【００２４】
本発明のプログラムを起動し、被判別蛋白質を形成するＮ個のアミノ酸列を入力すると、記憶部に貯えられている前記表１のデータと照らし合わされ、各アミノ酸に、それぞれ固有の電荷指標が割り振られる（ＳＴＰ１００）。次いで、前記（１）式によって総電荷指標の絶対値|CT|が計算される（ＳＴＰ１０１）。次に、得られた|CT|が２０以上であるか否かを判別し（ＳＴＰ１０２）、２０未満であればプログラムを終了し、必要に応じて新しい被判別蛋白質について、新たに操作を開始する。
【００２５】
|CT|が２０以上である場合には、（２）式によって蛋白質全体の電荷密度の絶対値|D|を計算する（ＳＴＰ１０３）。但し、（２）式中のＮは、被判別蛋白質の全アミノ酸残基数（アミノ酸鎖の長さ）である。次に|D|が０．１４以上であるか否かを判別する（ＳＴＰ１０４）。|D|が０．１４未満である場合にはプログラムは終了する。|D|が０．１４以上である場合には、Ｎ個からなる全アミノ酸配列を半分に分割し、分割されたそれぞれについて、総電荷指標CN（Ｎ端側）及びCC（Ｃ端側）を計算する（ＳＴＰ１０５）。
【００２６】
次に、上記CNとCCの積CN×CCが１００以上か否かを判別し（ＳＴＰ１０６）、１００未満であればプログラムを終了する。CN×CCが１００以上の場合には、表１に照らし合わせて、Ｎ個の全アミノ酸にそれぞれ固有のＫ−Ｄ疎水性指標Ｈと両親媒性指標Ａを割り当てる（ＳＴＰ１０７）。これらの割り当ては、ＳＴＰ１００で予め行っておいても良い。７残基ウィンドウを用い、前記式（４−１）によって一回平均疎水性[H]を計算する（ＳＴＰ１０８）。
【００２７】
次に、[H]の値が連続１０残基以上−０．２５以下となるＣｏｒｅ領域が存在するか否かを判別し（ＳＴＰ１０９）、存在しなければプログラムを終了する。Ｃｏｒｅ領域が存在する場合には、該Ｃｏｒｅ領域に対して、前記式（５−１）によって二回平均疎水性<H>を計算すると共に、式（５−２）によって二回平均両親媒性<A>を計算する（ＳＴＰ１１０）。尚、ＳＴＰ１０８における[A]の計算は、ＳＴＰ１０８では行わずにＳＴＰ１１０の<A>の計算の前に行うようにし、[A]から<A>迄の計算を一貫して行わせても良い。
【００２８】
次に、<H>≦−１．１５且つ<A>≧０．９０が満足されるか否かを判別し、満足されない場合にはプログラムを終了する。満足される場合には、Ｃｏｒｅ領域の両端３残基を除いた領域について、[H]の最大値Maxと最小値Minを求める（ＳＴＰ１１２）。次いでMax−Minが２．以下であるか否かを判別し（ＳＴＰ１１３）、２．０を超えればプログラムを終了する。
【００２９】
Max−Minが２．０以下となる場合には、Ｃｏｒｅ領域外の両端に向って[H]の値を探索していき、−０．３０以下になる直前の残基を候補セグメント末端と決定する。（ＳＴＰ１１４）。次に候補セグメントを除いたＮ端側総残基数NNとＣ端側総残基数NCを計算し（ＳＴＰ１１５）、NC≦３／２×NNかつNC≦２／３×NNを満足するか否かを判別する（ＳＴＰ１１６）。上記の条件が満たされない場合にはプログラムを終了する。
【００３０】
上記の条件が満足される場合、候補セグメントの総電荷指標CH、候補セグメントを除いたＮ端側の総電荷指標CND、及び、Ｃ端側の総電荷指標CCDを計算する（ＳＴＰ１１７）。次に、CH、CND及びCCDが何れも正又は負となるか判断する（ＳＴＰ１１８）。これらの指標の符号が一致して正又は負とならない場合には、プログラムは終了する。上記の条件が満足される場合には、候補セグメントをダンベル型水溶性蛋白質の中央ヘリックスと判断し（ＳＴＰ１１９）、ダンベル型水溶性蛋白質であるとの判定と中央ヘリックス領域を出力し、プログラムを終了する。
尚、本プログラム中で行った判断の各結果を同時に出力するようにできることは当然である。
【００３１】
下記表２に、１２個のダンベル型蛋白質について、本発明のコンピュータプログラムを用いて解析した結果を示す。
【表２】
0004021153.tif 000014
【００３２】
【発明の効果】
本発明はダンベル型水溶液蛋白質の二次構造の特徴を考慮し、新たに各アミノ酸に対して設定した、電化荷指標、K-D疎水性指標、及び両親媒性指標を用いて、特にダンベル型水溶性蛋白質に絞って判別を行うので、従来になく、迅速且つ精度良く判別することができる。
【図面の簡単な説明】
【図１】ダンベル型蛋白質の定義を表す図である。
【図２】本発明のコンピュータプログラムのフローチャートの内、ダンベル型蛋白質の候補を抽出するための部分フローチャートである。
【図３】本発明のコンピュータプログラムのフローチャートの内、ダンベル型蛋白質の候補について、中央ヘリックスの候補となるセグメントを抽出するための部分フローチャートである。
【図４】本発明のコンピュータプログラムのフローチャートの内、中央ヘリックスとダンベル型蛋白質についての最終判断をするための部分フローチャートである。
【図１】
0004021153.tif 000015
【図２】
0004021153.tif 000016
【図３】
0004021153.tif 000017
【図４】
0004021153.tif 000018