計算生命科学の基礎10
2023年10⽉11⽇

タンパク質の機能の進化とその情報解析
藤博幸

関⻄学院⼤学
⽣命環境学部⽣物科学科

本⽇の予定1. タンパク質2. タンパク質の機能
3. 情報解析に使⽤するタンパク質のデータ
4. 配列解析5. 構造解析
6. ネットワーク解析

図1.2 DNA、RNA、タンパク質の生合成とセントラルドグマ

5’…ATGGAACATCTGACT…3’
3’…TACCTTGTAGACTGA…5’DNA
転写

5’…AUGGAACAUCUGACU…3’

逆転写

RNA

リボソーム

翻訳

タンパク質
M E H V T…

DNA複製

RNA複製

http://www.wdic.org/w/SCI/ペプチド結合

タンパク質は20種類のアミノ酸がペプチド結合でつながった紐状の分⼦
各アミノ酸を１⽂字で表現すると、⼀つのタンパク質は、20種類のアルファベットで構成された⽂字列として表すことができる。

アミノ酸
Asp
(D)
Glu (E)
アスパラギン酸
グルタミン酸

側鎖

H
N
H
アミノ基

R

主鎖

C

OH

H

C

Tyr
(Y)
チロシン

Ser
(S)
セリン

解離性アミノ酸
Lys
(K)
リジン

His (H)
ヒスチジン

Thr (T)
スレオニン

Leu(L) イソロイシン
Ile (I)
ロイシン

親水性
アミノ酸

Arg (R)
アルギニン

O

カルボキシ基

基本構造

Gly
(G)
グリシン

Ala
(A)
アラニン

Val
(V)
バリン

Asn (N)
アスパラギン

アミノ酸は
”大文字”
で表記する

疎水性アミノ酸
Met
(M)
メチオニン

p.4, 3

Cys
(C)
システイン

Pro
(P)
プロリン

Phe (F)
フェニルアラニン

Trp (W)
トリプトファン

Gln
(Q)
グルタミン

http://www.wdic.org/w/SCI/ペプチド結合

タンパク質は20種類のアミノ酸がペプチド結合でつながった紐状の分⼦
各アミノ酸を１⽂字で表現すると、⼀つのタンパク質は、20種類のアルファベットで構成された⽂字列として表すことができる。

リゾチームのアミノ酸配列>LYC_HUMAN ヒト Lysozyme C
MKALIVLGLVLLSVTVQGKVFERCELARTLKRLGMDGYRGISLANWMCLAKWESGYNTRATNYNAGDRSTDYGIFQINSRYWCNDGKTPGAVNACHLSCSALLQDNIADAVACAKRVVRDPQGIRAWVAWRNRCQNRDVR
QYVQGCGV>LYC1_BOVIN ウシ Lysozyme C 1
MKALIILGFLFLSVAVQGKVFERCELARTLKKLGLDGYKGVSLANWLCLTKWESSYNTKATNYNPGSEST
DYGIFQINSKWWCNDGKTPNAVDGCHVSCSELMENDIAKAVACAKQIVSEQGITAWVAWKSHCRDHDVSSYVEGCTL>LYC_CHICK ニワトリ Lysozyme C
MRSLLILVLCFLPLAALGKVFGRCELAAAMKRHGLDNYRGYSLGNWVCAAKFESNFNTQATNRNTDGSTDYGILQINSRWWCNDGRTPGSRNLCNIPCSALLSSDITASVNCAKKIVSDGNGMNAWVAWRNRCKGTDVQAWIRGCRL>LYC2_ONCMY マス Lysozyme C IIMRAVVVLLLVAVASAKVYDRCELARALKASGMDGYAGNSLPNWVCLSKWESSYNTQATNRNTDGSTDYGIFQINSRYWCDDGRTPGAKNVCGIRCSQLLTADLTVAIRCAKRVVLDPNGIGAWVAWRLHCQNQDLRSYVA
GCGV>LYC_BOMMO カイコ Lysozyme
MQKLIIFALVVLCVGSEAKTFTRCGLVHELRKHGFEENLMRNWVCLVEHESSRDTSKTNTNRNGSKDYGL
FQINDRYWCSKGASPGKDCNVKCSDLLTDDITKAAKCAKKIYKRHRFDAWYGWKNHCQGSLPDISSC>LYSP_DROME ハエ Lysozyme P
MKAFLVICALTLTAVATQARTMDRCSLAREMSKLGVPRDQLAKWTCIAQHESSFRTGVVGPANSNGSNDYGIFQINNKYWCKPADGRFSYNECGLSCNALLTDDITNSVKCARKIQRQQGWTAWSTWKYCSGSLPSINSCFFASTA フォーマット
>LYC_HUMAN ヒト Lysozyme C
MKALIVLGLVLLSVTVQGKVFERCELARTLKRLGMDGYRGISLANWMCLAKWESGYNTRATNYNAGDRSTDYGIFQINSRYWCNDGKTPGAVNACHLSCSALLQDNIADAVACAKRVVRDPQGIRAWVAWRNRCQNRDVRQYVQGCGV“>” で始まる注釈⾏（通常、配列名）と改⾏後に、１⽂字表記で塩基あるいはアミノ酸の配列を記す書式を
FASTAフォーマットといい、分⼦系統解析ばかりでなく、バイオインフォ
マティクス分野では配列の書式としてよく使われている。
前ページのように、⼀つのファイルに複数のFASTAフォーマットの配列がおさめられている場合は、マルチFASTAフォーマットとよぶ。
p.9,81

配列データベース
塩基配列データベース
INSDC (International Nucleotide Sequence Database Collaboration)
⽶国 NCBIのGenBank
欧州 EMBL-EBIのENA
⽇本 国⽴遺伝学のDDBJ3者が連携して国際的に維持管理

タンパク質のデータベース
UniProt
以下の３者のコンソーシアムとして統合
スイスバイオインフォマティクス研究所(SIB)のSwiss-Prot
EBIのTrEMBL
ジョージタウン⼤学のPIR

塩基配列データベースのエントリ : GenBank形式 (タンパク質の場合: GenPept形式 ）LOCUS
DEFINITION

AH003092
6907 bpDNAlinear
ROD 01-AUG-2016
Nannospalax ehrenbergi clone pSCr-1 crystallin (CRY-alpha-A) gene,complete cds, alternatively spliced.
ACCESSION
AH003092 M17247 M17248 M17249 Y00464
アクセッション番号
VERSION
AH003092.2
キーワード
KEYWORDS
alpha-crystallin; crystallin.
SOURCE
Nannospalax ehrenbergi (Ehrenberg's mole-rat)
ORGANISM Nannospalax ehrenbergi
Eukaryota; Metazoa; Chordata; Craniata; Vertebrata; Euteleostomi;
Mammalia; Eutheria; Euarchontoglires; Glires; Rodentia; Myomorpha;
Muroidea; Spalacidae; Spalacinae; Nannospalax.
REFERENCE
1 (bases 2161 to 2880; 3841 to 4080; 5592 to 6431)
AUTHORS
Hendriks,W., Leunissen,J., Nevo,E., Bloemendal,H. and de Jong,W.W.
TITLE
The lens protein alpha A-crystallin of the blind mole rat, Spalax
ehrenbergi: evolutionary change and functional constraints
JOURNAL
Proc. Natl. Acad. Sci. U.S.A. 84 (15), 5320-5324 (1987)
PUBMED3474658:
FEATURES
Location/Qualifiers
source1..6907/organism="Nannospalax ehrenbergi"/mol_type="genomic DNA”:
CDS
join(2349..2537,2736..2804,3859..3981,5715..5924)
コード領域/gene="CRY-alpha-A"/codon_start=1/product="alpha-A-ins-crystallin"/protein_id="AAA66165.1"/translation="MDVTIQHPWFKHALGPFYPSRLFDQFFGQGLFEYDLLPFLSSTISPYYRQTLLRTVLDSCISELMTHRWFVPHQPHAGNPENNPIKVRSDRDKFVIFLDVKH
アミノ酸配列FSPEDLTVKVLEDFVEIHGKHNERQDDHGYISREFHRRYRLPSSVDQSALSCSLSADG（ins あり）
MLTFSGPKVQSGLDAGHSERAIPVSQEEKPSSAPLF"
CDS
join(2349..2537,3859..3981,5715..5924)
コード領域/gene="CRY-alpha-A"/codon_start=1/product="alpha A-crystallin"/protein_id="AAA66166.1"/translation="MDVTIQHPWFKHALGPFYPSRLFDQFFGQGLFEYDLLPFLSSTISPYYRQTLLRTVLDSCISEVRSDRDKFVIFLDVKHFSPEDLTVKVLEDFVEIHGKHNE
アミノ酸配列RQDDHGYISREFHRRYRLPSSVDQSALSCSLSADGMLTFSGPKVQSGLDAGHSERAIP（ins なし）
VSQEEKPSSAPLF”:
ORIGIN
1 bp upstream of EcoRI site.
1 gaattctctc ttgtggagca ggccttaaat ccaagcagaa agtggtttct cacccccata:
2281 gctccccact ccagcggccc ctgcccgact caccaccagc cttcagtccc tgcgttgcca
2341 cgtggaacat ggacgtcacc atccagcacc cctggttcaa gcacgccctg ggtcccttct
2401 atcccagccg gctgttcgac cagtttttcg gccagggcct ctttgagtac gatctgctgc
2461 ccttcctgtc ctccaccatc agcccctatt accgccagac cctcctccgc acagtgctgg
2521 actcctgcat ctctgaggta agatgcggct ggtcagaggg gctcagggcc ttggagaggt
2581 cagctaggga ttcctgaggg gccacctgcc tcgtccctgc ctctcccagc ctcggaccct
イントロン１
2641 ggaccctgct gcggctgggt ggaagtccgc ccgcccagcc gactgggaat ggggtttgaa
2701 gagaaggccc gactgctgtc ttttcctccc atcagctcat gacccatagg tggtttgtac
2761 cgcaccaacc acatgctgga aaccccgaga acaaccccat caaggcaagt ttcgtggcaa
2821 aggcgagcgt cctggacccg ctgctccctg tgcttttctg ccgggtgcgg gtgcggtgat:
イントロン２3781 accatgctgt tgcccctccc acccccgttc ctggagctgt tagcctctgg gagctccaca
3841 tgagggtttg ggtttcaggt ccgatccgac cgggacaaat ttgtcatctt cctggatgtg
3901 aagcacttct ctcctgagga cctcaccgtg aaggtgctgg aggacttcgt ggagatccat
3961 ggcaagcaca acgagaggca ggtgagctcg agtgggacct cacactgcgg gcatcctgta
4021 tcccttgaga ggaagcaagg ggaggtggca ctggctcagc tgctcctgct agaaggagcc:
イントロン３5641 ggtctttttc cagtgtggga tggcagtcct gtgtctcatc atgcctgacc atgctcttcc
5701 ctggttgtcc ccaggacgac cacggctaca tttcccgtga gttccaccgc cgataccgcc
5761 tgccttccag tgtggaccag tctgcccttt cctgctccct gtctgcagat ggcatgctga
5821 ccttctctgg ccccaaagtc caatctggcc tggatgctgg ccacagtgaa agggccattc
5881 ctgtgtctca ggaggagaag cccagctcgg cacccttgtt ctaagcaggc ctcaccttag
5941 ctgccccctg aggccactga cccatcaagc tcagggacca cagcaaaggg tcttccttcc:
6841 acacacacgc atacacatga cacatatact actgcgcatg cacagcacac acacactgat
6901 caccaca
//

配列の定義
（学名、遺伝子名）

学名（生物資料の由来）
分類

文献（一部省略）（一部省略）（一部省略）（一部省略）

イントロン
１-ins-２（一部省略）（一部省略）（一部省略）

アミノ酸配列

ヒト・リゾチームのアミノ酸配列（一文字表記）>LYC_HUMANMKALIVLGLVLLSVTVQGKVFERCELARTLKRLGMDGYRGISLANWMCLAKWESGYNTRATNYNAGDRSTDYGIFQINSRYWCNDGKTPG
AVNACHLSCSALLQDNIADAVACAKRVVRD
PQGIRAWVAWRNRCQNRDVRQYVQGCGV

148 残基ヒト・リゾチームの立体構造
立体構造

タンパク質⽴体構造の階層性
立体構造をとる二本以上のペプチド鎖から形成される構造
四次構造（具体例：ヘモグロビンは４つのペプチド鎖が集まり機能 (αサブユニット×２＋βサブユニット×２）

三次構造

立体構造

立体構造

ドメイン構造

二次構造

一次構造

αへリックス

βシート

アミノ酸配列

ターン

ループ

⼆次構造

規則的な水素結合のパターンを持つ

αへリックスβシート

4残基先のアミノ酸と水素結合を作る

配列上離れた位置にアミノ酸と水素結合を作る『MOLECULAR CELL BIOLOGY』より引用

三次構造

内側 :
疎水性残基が集まる
外側 :
親水性残基が多い

ドメイン :
ヘモグロビンα鎖
独自の安定な構造をとる

・比較的安定な構造をとる単位・三次構造より（たいていの場合）小さな構造・ドメインは、様々な三次構造中に見られる

二次構造の組み合わせで、ある形に折りたたまれたもの。
一本のポリペプチドからなる。

四次構造
サブユニット

多量体の呼び方
１本 ： モノマー
２本 ： ダイマー
３本 ： トリマー
４本 ： テトラマー
５本 ： ペンタマー
６本 ： ヘキサマー

ヘモグロビン（2W6V）
生体内で機能しているタンパク質の多くは、四次構造を形成している

複数のポリペプチドによって形成される立体構造

⽣体⾼分⼦の⽴体構造を扱ったデータベース• 網羅的に立体構造情報を収集しているもの
ヨーロッパ
欧州バイオインフォマティクス研究所(EBI)
⾼分⼦構造データベース (MSD-EBI)PDBeRSCB PDB
PDBj

⽇本
⼤阪⼤学蛋⽩質研究所

wwPDB

アメリカ
構造バイオインフォマティ
クス
研究共同体 (RCSB)

⽣体⾼分⼦の⽴体構造を扱ったデータベース• PDBファイルの中身
HEADER
TRANSFERASE
30-AUG-04
1XBB
ボディ
TITLE
CRYSTAL STRUCTURE OFヘッダー
THE SYK TYROSINE KINASE DOMAIN WITH
TITLE
2 GLEEVECHEADER
物質の分類名、日付、ファイル名
COMPND
MOL_ID: 1;
COMPND KINASE SYK;
タンパク質名COMPND
2 MOLECULE: TYROSINE-PROTEIN
COMPND
3 CHAIN: A;
SOURCE
物質が由来する生物名
COMPND
4 SYNONYM: SPLEEN TYROSINE KINASE;
COMPND
5 EC: 2.7.1.112;AUTHOR
座標を作成した著者名
COMPND
6 ENGINEERED: YESREVDAT
データが登録された日
SOURCE
MOL_ID: 1;
SOURCE
2 ORGANISM_SCIENTIFIC:JRNLHOMO SAPIENS; 文献情報
SOURCE
3 ORGANISM_COMMON: HUMAN;
REMARK
実験条件などに関する情報
SOURCE
4 ORGANISM_TAXID: 9606;
SOURCE
5 GENE: SYK;SEQRES
アミノ酸配列SOURCE
6 EXPRESSION_SYSTEM: SPODOPTERA FRUGIPERDA;
HET
アミノ酸以外の原子
SOURCE
7 EXPRESSION_SYSTEM_COMMON:
FALL ARMYWORM;
SOURCE
8 EXPRESSION_SYSTEM_TAXID:
HELIX 7108;
αーヘリックスに関する情報
SOURCE
9 EXPRESSION_SYSTEM_CELL_LINE: SF9;SHEET
βーシートに関する情報
SOURCE 10 EXPRESSION_SYSTEM_VECTOR_TYPE:BACULOVIRUS
KEYWDSGLEEVEC, STI-571,
IMATINIB, SYK, SPLEEN TYPROSINE KINASE,
TURN
βーターンに関する情報
KEYWDS
2 ACTIVE CONFORMATION, STRUCTURAL GENOMICS, STRUCTURALSSBOND
ジスルフィド結合に関する情報
KEYWDS
3 GENOMIX, TRANSFERASE
EXPDTA
X-RAY DIFFRACTION
ATOM
個々のアミノ酸の原子の座標
AUTHOR
V.L.NIENABER,S.ATWELL,J.M.ADAMS,J.BADGER,M.D.BUCHANAN,HETATM
AUTHOR
2 I.K.FEIL,K.J.FRONING,X.GAO,J.HENDLE,K.KEEGAN,B.C.LEON,
アミノ酸以外の原子の座標
AUTHOR
3 H.J.MULLER-DEICKMANN,B.W.NOLAND,K.POST,K.R.RAJASHANKAR,CONTACT
原子の結合に関する情報
AUTHOR
4 A.RAMOS,M.RUSSELL,S.K.BURLEY,S.G.BUCHANAN（省略）
END
ファイルの終了

⽣体⾼分⼦の⽴体構造を扱ったデータベース• 座標データ
ATOM
ATOM
ATOM
ATOM
ATOM
ATOM
ATOM
ATOM
ATOM
ATOM
ATOM
ATOM
ATOM
ATOM
ATOM
ATOM
ATOM

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17

N VAL A 363
CA VAL A 363
C VAL A 363
O VAL A 363
CB VAL A 363
N TYR A 364
CA TYR A 364
C TYR A 364
O TYR A 364
CB TYR A 364
CG TYR A 364
CD1 TYR A 364
CD2 TYR A 364
CE1 TYR A 364
CE2 TYR A 364
CZ TYR A 364
OH TYR A 364

原⼦の名前
アミノ酸の名前チェインの名前1原子22.74121.55720.95419.73721.88321.79821.31020.92921.73522.34922.61921.72523.77221.95924.01923.10123.330-1.397 11.729 1.00 33.32-0.831 11.024 1.00 32.13-1.757 9.943 1.00 31.73-1.906 9.845 1.00 30.94
0.552 10.391 1.00 33.45-2.389 9.135 1.00 29.77-3.035 7.928 1.00 27.96-4.485 8.157 1.00 26.43-5.274 8.650 1.00 29.31-2.901 6.808 1.00 28.13-1.461 6.442 1.00 28.30-0.751 5.638 1.00 30.33-0.815 6.893 1.00 30.80
0.588 5.285 1.00 31.59
0.523 6.541 1.00 32.22
1.208 5.733 1.00 32.48
2.524 5.389 1.00 36.29

温度因⼦

y座標
x座標
アミノ酸残基番
号

N
C
C
O
C
N
C
C
O
C
C
C
C
C
C
C
O

占有率
z座標

１残基

１残基

原⼦

mmCIF形式
現在のPDBのフォーマット
旧PDBフォーマットは1⾏80⽂字の固定幅であったため、原⼦番号、
残基番号、チェイン番号などの⽂字数に制限があり、100個以上チェインからなるような巨⼤なタンパク質複合体の構造を取り扱えなくなった。

可視化ツールが必要• 分子グラフィクス
RasMolhttp://www.umass.edu/microbio/rasmol/
PyMol
http://pymol.sourceforge.net/
Chime
http://www.mdli.com/downloads/downloadable/index.jsp
PDBj viewer http://www.pdbj.org/PDBjViewer/index_j.html
MolScript
http://www.avatar.se/molscript/• 粗視化コンタクトマップ
ラマチャンドランプロット
TOPS
http://www.tops.leeds.ac.uk/
LIGPLOT
http://www.biochem.ucl.ac.uk/bsm/ligplot/ligplot.html

チロシンキナーゼの構造がPyMolのウィンドウに表示される
赤いドットは水タンパク質の機能
⽣化学的機能vs
Biochemical Function
酵素活性、リガンド結合能

⽣物学的機能Biological Function
記憶、発⽣、⾏動XX

相互作⽤データベース
STRING
IntActBioGRID

STRING
https://string-db.org

クリック

BioGRID

https://thebiogrid.org

ダウンロードされるファイルには⼆項関係の形で相互作⽤を記述
タンパク質A タンパク質Bタンパク質A タンパク質Fタンパク質C タンパク質D...

進化的視点からのタンパク質の機能解析
相同タンパク質の⽐較解析相同性 (homology)
共通の起源から分岐してきた遺伝⼦の間の関係性

祖先遺伝⼦

分岐

相同

アデロバシレウス （Adelobasileus cromptoni

２種類の分岐（１）種分化 (speciation)）

2億2500万年前

遺伝⼦

A

7500万年前
遺伝⼦ A”
遺伝⼦ Aʼ

（２）遺伝⼦重複 (gene duplication)

重複した時には同じコピー

分⼦進化 (molecular evolution)
進化の過程で配列が変化すること向上、複雑化などの意味は持たない（分⼦進化だけでなく、進化⾃体も同じ）

分子進化の基本ステップ
点変異･挿入変異･欠失変異（a） 点変異(point mutation)

5’...ATGGTGCTGTC...3’5’...ATGGTACTGTC...3’（b） 挿入変異(insertion)（c） 欠失変異(deletion)

5’...ATGGTGCTGTC...3’5’...ATGGTGCTGTC...3’5’...ATGGTGAGCATACTGTC...3’5’...ATGTC...3’

種分化に伴う相同タンパク質の形成

種分化
10
20
30
40
50
human VLSPADKTNVKAAWGKVGAHAGEYGAEALERMFLSFPTTKTYFPHFDLSH
祖先型の生物a::: ::.:.::::::.:.:..::::::::::: :::::::::::::.::
mouse VLSGEDKSNIKAAWGKIGGHGAEYGAEALERMFASFPTTKTYFPHFDVSH
10
20
30
40
50
60
70
80
90
100
human GSAQVKGHGKKVADALTNAVAHVDDMPNALSALSDLHAHKLRVDPVNFKL
突然変異::::::::::::::::..:..:.::.:.::::::::::::::::::::::
mouse GSAQVKGHGKKVADALASAAGHLDDLPGALSALSDLHAHKLRVDPVNFKL
60
70
80
90
100

種分化
110
120
130
140
human LSHCLLVTLAAHLPAEFTPAVHASLDKFLASVSTVLTSKYR::::::::::.: ::.:::::::::::::::::::::::::
mouse LSHCLLVTLASHHPADFTPAVHASLDKFLASVSTVLTSKYR
110
120
130
140
a

進化

20/141a

生物の種分化を反映したタンパク質の増加

突然変異(mutation)と置換（substitution)

突然変異は集団中の個体に⽣じる※ 進化に寄与するのは体細胞ではなく⽣殖系列の細胞に⽣じる突然変異

突然変異(mutation)と置換（substitution)

有害な突然変異は、集団から除去される（負の選択 or 純化淘汰）

突然変異(mutation)と置換（substitution)

有利な突然変異は急速に集団中に広まり集団全体がその突然変異遺伝⼦で置き換えられる （正の選択）突然変異(mutation)と置換（substitution)

中⽴な突然変異の頻度はランダムに浮動し確率的に集団中に固定 （選択的に中⽴）固定にいたる途中の過程

多型(polymorphism)タンパク質をコードする遺伝⼦における分⼦レベルの変化の⼤部分は中⽴的

中⽴な変化

⾃然選択

遺伝⼦重複に伴う相同タンパク質の形成

遺伝⼦重複古典的な遺伝⼦重複の描像

重複した時には同じコピー

多くの変異は、偽遺伝⼦を⽣じ、まれに新機能を持った遺伝⼦が⽣じる
タンパク質が新たな機能を獲得する進化機構

2R WGD 仮説
2 round whole genome duplication
脊椎動物のゲノムは進化の初期の段階で２回の全ゲノム重複を起こしたという
仮説
⼤野乾により「Evolution by Gene Duplication」(1970)で提唱された。

Susumu Ohno
(1928 - 2000)

Evolution by Gene Duplication (1970)

DDC (duplication, degeneration, complementation) model
３つのシナリオよりなる
Gene A

Gene A'

遺伝⼦重複 (gene duplication)

duplication
Non-functionalization

Gene A(pseudogenization)
Neofunctionalization

Subfunctionalization

Non-processed
pseudogene A（偽遺伝子）
Function A

Gene A

Gene B

Gene A'

Gene A''

古典的モデルと同じ
Functions A + B

Functions A' + A'' = A

Phase I of DDC モデル

⼩さなボックスが発現制御領域を表す。⼤きなボックスはコード領域を表す。
Phase Iにおいて、subfunctionalizationでは、制御領域のいくつか壊れ、重複したコピーは相補的な発現パターンを⽰し、1コピーの時の機能を維持できる。
Neofunctionalizationにおいても、制御領域が変化することが知られてている。

Phase II of DDC モデル

Force et al. (1999)
Neofunctionzalitionまたsubfunctinalizationで⽣じたコピーの持つ制御領域の中の冗⻑な領域が壊れ、発現パターンの分化が進む。
Pseufogenizationは、Phase IIには関係しない

祖先遺伝子遺伝子重複による
α、βの形成ヒト、ラットの種分化

paralogousヒト

ラット ヒトヘモグロビンαラット

ヘモグロビンβorthologousオーソロガス (orthologous)
種分化に伴い分岐した相同な関係

相同(homologous)

パラロガス(paralogous)

遺伝⼦重複によって分岐した相同な関係

何故、相同タンパク質の解析をするのか？
“A major objective of the analysis of protein sequences is to discover which share
a common ancestor. Determining an evolutionary relationship is useful because it
generally implies a common or similar function.”
Hubbard,T.J.P. (1997) Current Opinion in Str. Biol. 7, 190-193.

相同配列の⽐較解析の⼿順
第⼀ステップ：相同配列の収集

第⼆ステップ：相同配列のマルチプルアラインメント

第三ステップ：アラインメントからの機能・構造・進化的情報の抽出

相同配列の⽐較解析の⼿順
第⼀ステップ：相同配列の収集

第⼆ステップ：相同配列のマルチプルアラインメント

第三ステップ：アラインメントからの機能・構造・進化的
情報の抽出

データベース検索
配列A
配列B (機能既知）
配列DB

問い合わせ配列

配列C
配列D (構造既知）

問い合わせ配列に類似した配列を配列DBの中から検索
類似配列の中で、構造や機能既知のものがあれば、それらと同様の構造
あるいは機能を有するものと推測発癌遺伝⼦ v-sis

配列データベース

⾎⼩板由来成⻑因⼦(PDGF)の配列決定

v-sis gene
Simian Sarcoma Virus

サル⾁腫ウイルス
v-sis は、成⻑因⼦としての情報を送り続けることで癌化を引き起こす。

GPCRsMonomer function
PeptidesAminesNucleic acids

Coupling withtrimeric G-protein
Conformation
changeSecondarymessengerbaG-protein
activation

g

l n
a
n tio
g
i
S uc
d
s
n
a
r
T

GPCRs• Membrane proteins• Bindneurotransmitters(physiologically
activepeptides,amines, nucleic acids, etc).
• Ligand binding to GPCRs causes theirconformation changes.
• It
leads
to
severalsignaltransductions
conjugatedwithtrimeric G-proteins.

GPCRs• About 1000 genes in human genome• Target for ~45% of clinically marketed drugs• Divided into 5 classes based on sequencesimilarity (Class A-E, the other)• Atomically resolved structure in class A
GPCR: Bovine Rhodopsin

既知のGPCRの配列

ゲノムにコード
されているタンパク質との類似性を検索

ヒト・ゲノム
データベース

既知のGPCRの配列

創薬ターゲットとしての検討新規のGPCRの配列

GPCR以外のターゲット探索の例
中外製薬(株）中外分⼦医学研究所
ドラッグデザインのターゲットとしての新規サイトカイン及び新規サイトカイン受容体のゲノム配列、EST配列からの探索「ゲノム創薬」 個別化医療とゲノムデータマイニング
野村仁 著 サイエンス社 (2005)

BLAST (Basic Local Alignment Search Tool)
https://blast.ncbi.nlm.nih.gov/Blast.cgi

相同配列の⽐較解析の⼿順
第⼀ステップ：相同配列の収集

第⼆ステップ：相同配列のマルチプルアラインメント

第三ステップ：アラインメントからの機能・構造・進化的情報の抽出

2. 配列アラインメント (sequence alignment)
進化の過程でのINDELを考慮しながら、相同な配列の間の対応する塩基（あるいはアミノ酸）を同じ位置に並べる操作
あるいは、その操作によってできたもの。
INDELに対応してギャップ(gap)とよばれる空記号を挿⼊し位置をずらして、塩基やアミノ酸を対応づける。
通常、動的計画法(dynamic programing algorithm)や、そのバリエーションを⽤いて、配列間の類似度が⾼くとなるように
アラインメントが構築される。
アラインメントの原理は次回

マルチプルアラインメント(multiple alignment)

分⼦系統樹を構築するには、まず複数本の相同配列についてのマルチプルアラインメントを作成する。

リゾチームのアミノ酸配列>LYC_HUMAN ヒト Lysozyme C
MKALIVLGLVLLSVTVQGKVFERCELARTLKRLGMDGYRGISLANWMCLAKWESGYNTRATNYNAGDRSTDYGIFQINSRYWCNDGKTPGAVNACHLSCSALLQDNIADAVACAKRVVRDPQGIRAWVAWRNRCQNRDVR
QYVQGCGV>LYC1_BOVIN ウシ Lysozyme C 1
MKALIILGFLFLSVAVQGKVFERCELARTLKKLGLDGYKGVSLANWLCLTKWESSYNTKATNYNPGSEST
DYGIFQINSKWWCNDGKTPNAVDGCHVSCSELMENDIAKAVACAKQIVSEQGITAWVAWKSHCRDHDVSSYVEGCTL>LYC_CHICK ニワトリ Lysozyme C
MRSLLILVLCFLPLAALGKVFGRCELAAAMKRHGLDNYRGYSLGNWVCAAKFESNFNTQATNRNTDGSTDYGILQINSRWWCNDGRTPGSRNLCNIPCSALLSSDITASVNCAKKIVSDGNGMNAWVAWRNRCKGTDVQAWIRGCRL>LYC2_ONCMY マス Lysozyme C IIMRAVVVLLLVAVASAKVYDRCELARALKASGMDGYAGNSLPNWVCLSKWESSYNTQATNRNTDGSTDYGIFQINSRYWCDDGRTPGAKNVCGIRCSQLLTADLTVAIRCAKRVVLDPNGIGAWVAWRLHCQNQDLRSYVA
GCGV>LYC_BOMMO カイコ Lysozyme
MQKLIIFALVVLCVGSEAKTFTRCGLVHELRKHGFEENLMRNWVCLVEHESSRDTSKTNTNRNGSKDYGL
FQINDRYWCSKGASPGKDCNVKCSDLLTDDITKAAKCAKKIYKRHRFDAWYGWKNHCQGSLPDISSC>LYSP_DROME ハエ Lysozyme P
MKAFLVICALTLTAVATQARTMDRCSLAREMSKLGVPRDQLAKWTCIAQHESSFRTGVVGPANSNGSNDYGIFQINNKYWCKPADGRFSYNECGLSCNALLTDDITNSVKCARKIQRQQGWTAWSTWKYCSGSLPSINSC
F

マルチfasta形式>gi|443546|pdb|7HVP|A Chain A, Protease Of Human Immunodeficiency Virus 1
PQITLWQRPLVTIRIGGQLKEALLDTGADDTVLEEMNLPGKWKPKMIGGIGGFIKVRQYDQIPVEIXGHKAIGTVLVGPTPVNIIGRNLLTQIGXTLNF> HIV2 protease (NP_663784 REGION: 523..603)
vtayiedqpv evlldtgadd sivagielgd nytpkivggi ggfintkeyk nveikvlnkr
vratimtgdt pinifgrnil t> simian immunodeficiency virus SIV-mnd 2 (NP_758887 REGION: 72..166)
slwnrpttvv eiegqkveal ldtgaddtvi kdldlkgnwk pqiiggiggs invkqffnck
vtiagkttha svlvgptpvn ivgrnvlkkl gctln>gi|4389337|pdb|1BAI|A Chain A, Rous Sarcoma Virus Protease
LAMTMEHKDRPLVRVILTNTGSHPVKQRSVYITALLDTGADDTVISEEDWPTDWPVMEAANPQIHGIGGGIPVRKSRDMIELGVINRDGSLERPLLLFPLVAMTPVNILGRDCLQGLGLRLTNL>gi|224443|prf||1104339A MoMuLV protease
TLDDQGGQGQEPPPEPRITLKVGGQPVTFLVDTGAQHSVLTQNPGPLSDKSAWVQGATGGKRYRWTTDRKVHLATGKVTHSFLHVPDCPYPLLGRDLLTKLKAQIHFEGSGAQVMGPMGQPLQVL

1. 配列は、スペースがあいていても構わないことに注意
2. Mafftは⼤⽂字も、⼩⽂字も区別せずに配列を処理してくれるが、ソフトによっては⼤⽂字あるいは⼩⽂字しか処理できない場合があるかもしれない

MAFFTによるマルチプルアラインメント
mafftは宮⽥研究室で開発され、加藤和貴によって継続的に開発されているマルチプルアラインメントのフリーソフトウェア
海外の多くの研究機関で利⽤されている。
Web上でのアラインメントサービスに加え、ダウンロードして⾃⾝のPC上で利⽤できる。Mac, Windows, Linuxなど様々なOSに対応している
ここでは、web serviceとして公開されているmafftを利⽤する。
https://mafft.cbrc.jp/alignment/software/リゾチームのアミノ酸配列 マルチプル・アラインメント： Clustal形式

CLUSTAL W (1.81) multiple sequence alignment
ヒト LYC_HUMAN
ウシ LYC1_BOVIN
マス LYC2_ONCMY
ニワトリ LYC_CHICK
カイコ LYC_BOMMO
ハエ LYSP_DROME-MKALIVLGLVLLSVTVQGKVFERCELARTLKRLGMDGYRGISLANWMCLAKWESGYNTR-MKALIILGFLFLSVAVQGKVFERCELARTLKKLGLDGYKGVSLANWLCLTKWESSYNTK----MRAVVVLLLVAVASAKVYDRCELARALKASGMDGYAGNSLPNWVCLSKWESSYNTQ-MRSLLILVLCFLPLAALGKVFGRCELAAAMKRHGLDNYRGYSLGNWVCAAKFESNFNTQ-MQKLIIFALVVLCVGSEAKTFTRCGLVHELRKHGFEENL---MRNWVCLVEHESSRDTSMKAFLVICALTLTAVATQARTMDRCSLAREMSKLGVPRDQ---LAKWTCIAQHESSFRTG:. .
.:. ** *. :*.
: :* * : **. *LYC_HUMAN
LYC1_BOVIN
LYC2_ONCMY
LYC_CHICKLYC_BOMMO
LYSP_DROMEATNYNAGDRSTDYGIFQINSRYWCNDGKTPGAVNACHLSCSALLQDNIADAVACAKRVVRATNYNPGSESTDYGIFQINSKWWCNDGKTPNAVDGCHVSCSELMENDIAKAVACAKQIVSATNRN-TDGSTDYGIFQINSRYWCDDGRTPGAKNVCGIRCSQLLTADLTVAIRCAKRVVLATNRN-TDGSTDYGILQINSRWWCNDGRTPGSRNLCNIPCSALLSSDITASVNCAKKIVS
KTNTN-RNGSKDYGLFQINDRYWCS--KGASPGKDCNVKCSDLLTDDITKAAKCAKKIYKVVGPANSNGSNDYGIFQINNKYWCKPADGRFSYNECGLSCNALLTDDITNSVKCARKIQR..
. *.***::***.::**.
. . * : *. *: ::: : **:::

LYC_HUMAN
LYC1_BOVIN
LYC2_ONCMY
LYC_CHICKLYC_BOMMO
LYSP_DROMEDPQGIRAWVAWRNRCQNRDVRQYVQGCGVE-QGITAWVAWKSHCRDHDVSSYVEGCTLDPNGIGAWVAWRLHCQNQDLRSYVAGCGVDGNGMNAWVAWRNRCKGTDVQAWIRGCRL
R-HRFDAWYGWKNHCQGS--LPDISSC-Q-QGWTAWSTWK-YCSGS--LPSINSCF:** *: * .
: .*

相同配列の⽐較解析の⼿順
第⼀ステップ：相同配列の収集

第⼆ステップ：相同配列のマルチプルアラインメント

第三ステップ：アラインメントからの機能・構造・進化的情報の抽出

第三ステップ：アラインメントからの機能・構造・進化的情報の抽出
モチーフの抽出：機能部位予測
重複遺伝⼦の機能差を決定しているサイトの予測
マルチプルアラインメント

分⼦系統解析
⽴体構造モデリング
機械学習の特徴ベクトル...

リゾチームのアミノ酸配列 マルチプル・アラインメント： Clustal形式

CLUSTAL W (1.81) multiple sequence alignment
ヒト LYC_HUMAN
ウシ LYC1_BOVIN
マス LYC2_ONCMY
ニワトリ LYC_CHICK
カイコ LYC_BOMMO
ハエ LYSP_DROME-MKALIVLGLVLLSVTVQGKVFERCELARTLKRLGMDGYRGISLANWMCLAKWESGYNTR-MKALIILGFLFLSVAVQGKVFERCELARTLKKLGLDGYKGVSLANWLCLTKWESSYNTK----MRAVVVLLLVAVASAKVYDRCELARALKASGMDGYAGNSLPNWVCLSKWESSYNTQ-MRSLLILVLCFLPLAALGKVFGRCELAAAMKRHGLDNYRGYSLGNWVCAAKFESNFNTQ-MQKLIIFALVVLCVGSEAKTFTRCGLVHELRKHGFEENL---MRNWVCLVEHESSRDTSMKAFLVICALTLTAVATQARTMDRCSLAREMSKLGVPRDQ---LAKWTCIAQHESSFRTG:. .
.:. ** *. :*.
: :* * : **. *LYC_HUMAN
LYC1_BOVIN
LYC2_ONCMY
LYC_CHICKLYC_BOMMO
LYSP_DROMEATNYNAGDRSTDYGIFQINSRYWCNDGKTPGAVNACHLSCSALLQDNIADAVACAKRVVRATNYNPGSESTDYGIFQINSKWWCNDGKTPNAVDGCHVSCSELMENDIAKAVACAKQIVSATNRN-TDGSTDYGIFQINSRYWCDDGRTPGAKNVCGIRCSQLLTADLTVAIRCAKRVVLATNRN-TDGSTDYGILQINSRWWCNDGRTPGSRNLCNIPCSALLSSDITASVNCAKKIVS
KTNTN-RNGSKDYGLFQINDRYWCS--KGASPGKDCNVKCSDLLTDDITKAAKCAKKIYKVVGPANSNGSNDYGIFQINNKYWCKPADGRFSYNECGLSCNALLTDDITNSVKCARKIQR..
. *.***::***.::**.
. . * : *. *: ::: : **:::

LYC_HUMAN
LYC1_BOVIN
LYC2_ONCMY
LYC_CHICKLYC_BOMMO
LYSP_DROMEDPQGIRAWVAWRNRCQNRDVRQYVQGCGVE-QGITAWVAWKSHCRDHDVSSYVEGCTLDPNGIGAWVAWRLHCQNQDLRSYVAGCGVDGNGMNAWVAWRNRCKGTDVQAWIRGCRL
R-HRFDAWYGWKNHCQGS--LPDISSC-Q-QGWTAWSTWK-YCSGS--LPSINSCF:** *: * .
: .*Clustal形式アラインメント下段のシンボルの意味“*”では，完全に保存“:”では，強い物理化学的類似性のあるグループで保存“.”では，弱い類似性のあるグル―プで保存強い弱いの基準は，PAM250 ⾏列において，アミノ酸間のスコアが0.5よ
り⼤きいか，0.5以下かで分けているスコア・テーブル多くのアミノ酸置換は保存的
Dayhoffのスコアテーブル
G
A
S
T
P
L
I
M
V
D
N
E
Q
F
Y
W
K
R
H
C

G
5
1
1
0
-1
-4
-3
-3
-1
1
0
0
-1
-5
-5
-7
-2
-3
-2
-3

A

S

T

P

L

I

M

V

D

N

E

Q

F

Y

W

K

R

H

C

2
1
1
1
-2
-1
-1
0
0
0
0
0
-4
-3
-6
-1
-2
-1
-2

2
1
1
-3
-1
-2
-1
0
1
0
-1
-3
-3
-2
0
0
-1
0

3
0
-2
0
-1
0
0
0
0
-1
-3
-3
-5
0
-1
-1
-2

6
-3
-2
-2
-1
-1
-1
-1
0
-5
-5
-6
-1
0
0
-3

6
2
4
2
-4
-3
-3
-2
2
-1
-2
-3
-3
-2
-6

5
2
4
-2
-2
-2
-2
1
-1
-5
-2
-2
-2
-2

6
2
-3
-2
-2
-1
0
-2
-4
0
0
-2
-5

4
-2
-2
-2
-2
-1
-2
-6
-2
-2
-2
-2

4
2
3
2
-6
-4
-7
0
-1
1
-5

2
1
1
-4
-2
-4
1
0
2
-4

4
2
-5
-4
-7
0
-1
1
-5

4
-5
-4
-5
1
1
3
-5

9
7
0
-5
-4
-2
-4

10
0
-4
-4
0
0

17
-3
2
-3
-8

5
3
0
-5

6
2
-4

6
-3

12

(1) G,A,S, T, P: small hydrophilic residues(2) L, I, M, V: hydrophobic residues
Dayhoffコード (3) D, N, E, Q: nagatively charged residues and the relatives(4) F, Y, W: aromatic residues(5) K, R, H: positively charged residues(6) C: Cys

Multiple Alignment of Lysozyme Amino Acid Sequences
CLUSTAL W (1.81) multiple sequence alignment
Human
BovineTrout
ChikcenSilkworm
Fruitfly

LYC_HUMAN
LYC1_BOVIN
LYC2_ONCMY
LYC_CHICKLYC_BOMMO
LYSP_DROME-MKALIVLGLVLLSVTVQGKVFERCELARTLKRLGMDGYRGISLANWMCLAKWESGYNTR-MKALIILGFLFLSVAVQGKVFERCELARTLKKLGLDGYKGVSLANWLCLTKWESSYNTK----MRAVVVLLLVAVASAKVYDRCELARALKASGMDGYAGNSLPNWVCLSKWESSYNTQ-MRSLLILVLCFLPLAALGKVFGRCELAAAMKRHGLDNYRGYSLGNWVCAAKFESNFNTQ-MQKLIIFALVVLCVGSEAKTFTRCGLVHELRKHGFEENL---MRNWVCLVEHESSRDTSMKAFLVICALTLTAVATQARTMDRCSLAREMSKLGVPRDQ---LAKWTCIAQHESSFRTG:. .
.:. ** *. :*.
: :* * : **. *LYC_HUMAN
LYC1_BOVIN
LYC2_ONCMY
LYC_CHICKLYC_BOMMO
LYSP_DROMEATNYNAGDRSTDYGIFQINSRYWCNDGKTPGAVNACHLSCSALLQDNIADAVACAKRVVRATNYNPGSESTDYGIFQINSKWWCNDGKTPNAVDGCHVSCSELMENDIAKAVACAKQIVSATNRN-TDGSTDYGIFQINSRYWCDDGRTPGAKNVCGIRCSQLLTADLTVAIRCAKRVVLATNRN-TDGSTDYGILQINSRWWCNDGRTPGSRNLCNIPCSALLSSDITASVNCAKKIVS
KTNTN-RNGSKDYGLFQINDRYWCS--KGASPGKDCNVKCSDLLTDDITKAAKCAKKIYKVVGPANSNGSNDYGIFQINNKYWCKPADGRFSYNECGLSCNALLTDDITNSVKCARKIQR..
. *.***::***.::**.
. . * : *. *: ::: : **:::

LYC_HUMAN
LYC1_BOVIN
LYC2_ONCMY
LYC_CHICKLYC_BOMMO
LYSP_DROMEDPQGIRAWVAWRNRCQNRDVRQYVQGCGVE-QGITAWVAWKSHCRDHDVSSYVEGCTLDPNGIGAWVAWRLHCQNQDLRSYVAGCGVDGNGMNAWVAWRNRCKGTDVQAWIRGCRL
R-HRFDAWYGWKNHCQGS--LPDISSC-Q-QGWTAWSTWK-YCSGS--LPSINSCF:** *: * .
: .*モチーフ: short stretch of conserved amino acid sequences
When several motifs are found by the sequence comparison,
the sequences are closely located in the 3 dimensional structure.

A

B

A
C

C

B

C
B
A

A

B
C

(a)

(b)

(c)

catalytic center of enzyme

Interface of complexhydrophobic coreCharacterization of Retroviral Proteases- Development of HIV Protease Inhibitor -CLUSTAL format alignment by MAFFT L-INS-i (v7.130b)
gi|443546|pdb|7 PQITLW----------QRPLVTIRIGGQL----------KEALLDTGADDTVLEEMNLPG
HIV2--------------------VTAYIEDQP----------VEVLLDTGADDSIVAGIELGD
simian---SLW----------NRPTTVVEIEGQK----------VEALLDTGADDTVIKDLDLKG
gi|4389337|pdb| LAMTMEHK--------DRPLVRVILTNTGSHPVKQRSVYITALLDTGADDTVISEEDWPT
gi|224443|prf|| ---TLDDQGGQGQEPPPEPRITLKVGGQP----------VTFLVDTGAQHSVLTQNPGPL: .
*:****:.:::
gi|443546|pdb|7 KW------KPKMIGGIGGFIKVRQ---YDQIPVEIXGHKAIGTVL----VGPTPVNIIGR
HIV2
NY------TPKIVGGIGGFINTKE---YKNVEIKVLNKRVRATIM----TGDTPINIFGR
simianNW------KPQIIGGIGGSINVKQ---FFNCKVTIAGKTTHASVL----VGPTPVNIVGR
gi|4389337|pdb| DWPVMEAANPQ-IHGIGGGIPVRKSRDMIELGVINRDGSLERPLLLFPLVAMTPVNILGR
gi|224443|prf|| SD------KSAWVQGATGGKRYRW---TTDRKVHLATGKVTHSFLH---VPDCPYPLLGR.
.. : * *
:
: :
..:
.
* :.**
gi|443546|pdb|7 NLLTQIGXTLN------------------F
HIV2
NILT-------------------------simian
NVLKKLGCTLN------------------gi|4389337|pdb| DCLQGLGLRLT-----------------NL
gi|224443|prf|| DLLTKLKAQIHFEGSGAQVMGPMGQPLQVL: *レトロウイルス・プロテアーゼのコンセンサス配列(Hydrophobic)-Asp-Thr-Gly-(Small Hydrophilic)Proteases with the Sequences
Similar to Retroviral Protease Motif

Input: RetroviralProtease Motif

Amino Acid Sequence
Database

Retroviral Proteasecaliflower mosaic virus
copi-like element 17.6 fruit fly
Moloney murine leulkaemia virus
human T-cell leukaemia virus
Rous sarcoma virus

LHCFV-DTGASLCIASLKCLI-DTGSTVNMTSVTFLV-DTGAQHSVLTIEALL-DTGADMTVLP
ITALL-DSGADITIIS

Aspartic Acid Protease (C-terminal domain)
pepsinogen（human）
CQAIV-DTGTSLLTGP
pepsinogen（pig）
CQAIV-DTGTSLLTGP
prochymosin（cow）
CQAIL-DTGTSKLVGP
penisilopepsin
FSGIA-DTGTTLLLLP
renin (mouse）
CEVVV-DTGSSFISAP
renin（human）
CLALV-DTGASYISGS
Aspartic Acid Protease (N-terminal domain)
pepsinogen（human）
DFTVVFDTGSSNLWVP
pepsinogen（pig）
DFTVIFDTGSSNLWVP
prochymosin (cow）
EFTVLFDTGSSDFWVP
penusilopepsin
TLNLNFDTGSSDFWVP
renin (mouse）
TFKVMFDTGSANLWVP
renin (human)
TFKVVFDTGSSNVWVP

モチーフ配列は、酸性プロテアーゼの活性中⼼に類似

N末ドメインとC末ドメインが
相同N

DD

ホモダイマーDC

予測

レトロウイルスのプロテアーゼは酸性プロテアーゼである。
H. Toh et al. EMBO J. 4, 1267 (1985)
H. Toh et al. Nature 315, 691 (1985)

実験[1] 既知酸性プロテアーゼの立体構造を鋳型としたHIV プロテアーゼのホモロジー・モデリン
グ L.H. Pearl and W.R. Taylor, Nature 329, 351 (1987).
[2] 酸性プロテーゼ阻害剤（ペプスタチン）によるHIV プロテアーゼの阻害
R.F. Nutt et al. Proc. Natl. Acad. Sci. USA 85, 7129 (1988).
P.L. Darke et al. J. Biol. Chem. 264, 2307 (1989).

検証
X 線結晶構造解析からHIV プロテアーゼは酸性プロテーゼと類似の構造である。
M. Miller et al. Science 246, 1149 (1989).

Coffin, JM (1997)

N末ドメインとC末ドメインが
相同N

DD

ホモダイマーDCmafftによる配列アラインメント
CLUSTAL format alignment by MAFFT L-INS-i (v7.402)

2APRN
2APRC

AGVGTVPMTDYGNDIEYYGQVTIGTPGKKFNLDFDTGSSDLWIASTLCTNCGSGQTKYDP---------------------------------YD-------------------STKFKG
:*
.**:.

2APRN
2APRC

NQSSTYQADGRT-WSISYGDGSSASGILAKDNVNLGGLLIKGQTI-----ELAKREAASFSLTTVPIDNSRGWWGITVDRATVGTSTVAS---SFDGILDTGTTLLILPNNIAASVARAY. ::.
:.* *.*: . .: .:. :*.
.:.*:* .* *:
::*
* ::

2APRN
2APRC

ASGPN-DGLLGLGFDT------ITTVRGV--KTPMDNLISQGLISRPIFGVYLGK----GASDNGDGTYTISCDTSAFKPLVFSINGASFQVSPDSLVFEEFQGQCIAGFGYGNWGFAI.:. * **:. **: ::.*. :.. *.*: : : .: * *. *:

2APRN
2APRC

-----AKNGGGGEYI-------------FGG
IGDTFLKN----NYVVFNQGVPEVQIAPVAE
**
:*:
..

進化的関係が遠いと配列情報だけでは正しいアラインメントが作れない
入力ファイルは 2aprNC.fasta

構造を⽐較することの意味
アミノ酸配列に⽐べ、⽴体構造は保存しやすい

アミノ酸配列は類似していない、つまり配列レベルの検索で検出できないような類似性(遠い進化的関係）を、
構造⽐較により検出できる。（構造データベース検索）
また、⽴体構造⽐較から、配列レベルでは検出できない機能部位の同定も可能（構造アラインメント）

図４

Δ＝0.40 exp(1.87H)
Δ：common core 主鎖の重ねあわせのRMSD
H：配列の相違度
Hが小さい時：表面のみの変化
Hが大きい時：表面と内部の変化
内部の変化は表面に比べ大きな構造変化をもたらす

構造⽐較のサイト

http://ekhidna2.biocenter.helsinki.fi/dali/87

10/7/2388

タンパク質の機能
⽣化学的機能vs
Biochemical Function
酵素活性、リガンド結合能

⽣物学的機能Biological Function
記憶、発⽣、⾏動X
X

ネットワークの解析のツール(1) R

igraph(2) Cytoscapeタンパク質間相互作⽤ネットワーク解析
1. ネットワーク可視化
2. 次数分布3. 中⼼性解析
4. コミュニティ抽出
5. …

ネットワークの可視化
Igraphによる可視化の例

次数分布
次数：各ノードに接続する辺の数

degree(f)gene0 gene1 gene2 gene3 gene4 gene5
2
3
2
1
2
4

次数分布を調べる

次数分布
log-log plot

log-log プロットをとると直線関係---à スケールフリー
であることを⽰している
少数のハブ(hub =次数の⾼いノード）と多数の次数の低いノードから構成されているスケールフリーネットワーク
少数のハブと多数の相互作⽤の数の少ないノードからなる。
ハブ
多数の相互作⽤相⼿を持つ
従来の鍵-鍵⽳モデルやinduced fitモデルでは説明できない
天然変性領域がハブの役割を担っていることが多いスケールフリー性は、社会学をはじめとするこれまでの研究により、現実世界のネットワークで幅広く観察されている。
・ ⼈々の持っている知⼈関係の数をみると、⼀部の⼈は⾮常にたくさんの知⼈を持っているが、⼤多数の⼈々の知⼈の数は限られている。
・WWWで はごく少数の有名サイトが数百万単位のリンクを集めているが、⼤多数のサイトはわずかなリンク先からしかリンクされていない。
・⽣体内の相互作⽤でも、ごく ⼀部のタンパク質が多数のタンパク質と相互作⽤する構造になっている。
数学的には、スケールフリー性は頂点が次数 k を持つ確率 p(k) の確率分布がp(k) ∝ k-γ のべき乗則になると表現されるlog( p(k) ) ∝ γ log(k-)
スケールフリーグラフが持つ注⽬すべき特性として、ネットワーク障害など「ランダムな故障や攻撃」に対して頑強性が⾼いことがあげられる。
ただし、スケールフリーなネットワークの場合、特定の重要なハブをピンポイントで狙った攻撃に対しては脆弱

https://ja.wikipedia.org/wiki/複雑ネットワーク

より

中⼼性
ノードがネットワーク中でどの程度中⼼的であるかを評価する指標
ネットワークの中の各ノードの重要性の評価に⽤いられる
様々な中⼼性の定義がある
次数中⼼性
近接中⼼性
媒介中⼼性

次数の⾼いノードが重要
他のノードに短い距離で繋がるノードが重要
ノード間最短距離を考えた時によく通過する
ノードが重要固有ベクトル中⼼性 隣接するノードの中⼼性も考慮
次数中⼼性の拡張版 (Google Page rankで使⽤）タンパク質間相互作⽤ネットワークの中⼼性解析
ecoli_proteins_essentiality_Baba2006MSB.txt
Baba et al. (2006) Mol Syst Biol 2, 0008

gene essential. . .
b0023 Nnon essential
b0024 uunknownb0025 Nb0026 Eessential. . .

中⼼性が⾼いものほどessentialなものが多い重複可能性(duplicability)とネットワークの関係(2012)進化速度と中⼼性の相関は⼀般には低いが、特定の機能カテゴリーについては、発現レベルと進化速度の間の相関と同程度かそれ以上の相関を⽰す。
ネットワーク内で相互作⽤するタンパク質をコードする遺伝⼦は⽐較的類似した進化速度を持つ
シロイヌナズナでは重複している遺伝⼦由来のタンパク質の⽅が、シングルトン（重複のない）遺伝⼦
由来のタンパク質の⽅が⾼度に相互作⽤が多く、ヒトでも同じ傾向。⼀⽅、E. coli, 酵⺟、線⾍、ハエでは重複遺伝⼦はネットワークの周辺に位置する傾向がある。
---à 中⼼性と重複可能性(duplicability)の関係が、真核⽣物の進化の過程で少なくとも２回逆転した

(2012)

中⼼性と重複可能性との関係の⽅向性は普遍的ではない：
Escherichia coli、酵⺟、ハエ、および線⾍では、重複した遺伝⼦はネットワークの周辺に存在する傾向
ヒトでは重複しが遺伝⼦はネットワークの中⼼的な位置を占める傾向
霊⻑類の異なる系統の各枝で起きた重複について調査
種系統樹の13の枝のうち、10個の枝で重複した遺伝⼦の⽅が、重複のない遺伝⼦よりも次数が⾼い
そのうち、５つの枝でその差は統計的に有意3個の枝（ヒトへ⾄る外部枝を含む）で逆のパターン。そのうち２つ（ヒトへ⾄る外部枝含む）でその差は統計的に有意次数以外の中⼼性でも同様の傾向が⾒られる。
中⼼性と重複可能性には関連はあるが、その関係性は霊⻑類の進化の過程で変化している。
相互作⽤するタンパク質の遺伝⼦は、同じ枝で同時に重複する傾向がある。

Conclusion
⽣化学的機能の解析:
アミノ酸配列データ、⽴体構造データ
⽣物学的機能の解析
相互作⽤データタンパク質の情報解析には、その他にも使えるデータがある（塩基配列、ゲノムデータ、etc)
異なるタイプのデータを組み合わせることで、個別のデータのみを⽤いるより豊かな情報を得ることができる