このページはまだ書きかけです.

はじめに

ここでは,HMMERを使って相同性検索を行います.

相同性検索

相同性検索については次のページで説明しています.

HMMERによる相同性検索を行う前に,BLASTによる相同性検索を先に理解しておくことをおすすめします.

HMMER

HMMERモチーフをクエリーとして相同性検索を行うツールです. (反対に配列をクエリーとしてモチーフを検索することもできますが,ここでは説明しません.)

モチーフとは,タンパク質のドメインを構造化したものです.

タンパク質のドメインは,タンパク質のアミノ酸配列の部分配列のうち,タンパ質の機能を発現するのに必要な部分です. 同じ機能をもつ(異なる種の)タンパク質には,同じ(あるいは,とても良く似た)ドメインが含まれています. ドメイン以外の部分は大きく異なっていても,ドメインが同じであれば(あるいは,とても良く似ていれば)同じ機能を発現します.

HMMERでは,モチーフはプロファイルHMMで表現されます.

プロファイルHMM

プロファイルHMMは,多樹配列アラインメントの結果を一つにまとめたものであるプロファイル隠れマルコフモデル (HMM, hidden Markov model) で表したものです.

状態遷移確率が直前の状態のみに依存するというマルコフ性を仮定して状態遷移確率をモデル化したものがマルコフ・モデルであり,マルコフ・モデルにおいて状態が観測できず(隠れていて)出力だけが観測できるものがHMMです.

たとえば,長さ3のドメインのHMMは次のようなものになります.

hmm.png

Sが開始状態,Eが終端状態,Mが一致状態,Iが挿入状態,Dが削除状態を表します. 出力は,一致状態と挿入状態ではアミノ酸のいずれか,削除状態では出力なし(ギャップだけど明示的には観測できない)となります.

トップ   新規 一覧 単語検索 最終更新   ヘルプ   最終更新のRSS