ペプチドの in silico コンフォメーション実験
ペプチドがどのようなコンフォメーションを取るか in silico 実験した.
ペプチドは乱数を使って生成し,libopenbabel を使ってコンフォメーション解析をおこなった.
コンフォメーション解析後のペプチドについて,分子の中心からの平均距離と分散,および X,Y,Z 座標に関する分散共分散行列の固有値を計算した.
コンフォメーションに電荷が大きく影響しているであろうことが確認できた.
目次(ページ内リンク)
実験
結果
実験
ペプチドの生成
ペプチドは自作ソフトウェア BCMD シリーズの BPeptGen で作成した.
鎖長は 10 で,ペプチド数は 20 である.アミノ酸は乱数を発生させて決定した.
コンフォメーション探索
OpenBabel のライブラリ版 libopenbabel を利用した簡単なソフトウェアを使っておこなった.
もちろんこの計算は,実行ファイル版 OpenBabel でも可能である.
ペプチド情報の取得
ペプチドの情報は,libbuilcule を利用したソフトウェアを作成して計算した.
項目は,電荷,分子の中心からの平均距離と分散,および X,Y,Z 座標に関する分散共分散行列の固有値である.
正電荷は,Lys,Arg,His の残基数,および N-末分 1 個である.
負電荷は,Asp,Glu の残基数,および C-末分 1 個である.
アミノ酸をランダムに発生すると,統計的に正電荷が多くなる.
ここでの固有値の意味は,分子を全体として楕円体と見たときの 3 軸の長さである.
結果
名前 | 配列 | 正電荷 | 負電荷 | 距離の平均 | 距離の分散 | 固有値0 | 固有値1 | 固有値2 | 固有値2 / 固有値0 |
---|---|---|---|---|---|---|---|---|---|
00.pdh | WCYQPPLPAA | 1 | -1 | 7.33861 | 9.43939 | 5.60652 | 9.55527 | 48.5359 | 8.6570 |
01.pdh | TSLSMKMTTF | 2 | -1 | 7.78562 | 8.65488 | 6.24088 | 7.40032 | 56.0572 | 8.9823 |
02.pdh | YDMRNIIEFF | 2 | -3 | 6.32376 | 6.4817 | 6.46609 | 6.98172 | 33.2764 | 5.1463 |
03.pdh | MDVDTLHVRH | 4 | -3 | 6.72893 | 5.09099 | 5.59977 | 13.5443 | 31.5235 | 5.6294 |
04.pdh | WMRIVFIIPT | 2 | -1 | 8.65772 | 14.371 | 6.7992 | 12.8336 | 70.1669 | 10.3199 |
05.pdh | NNMCGRAGVF | 2 | -1 | 6.59993 | 4.78737 | 4.77211 | 11.5099 | 32.4026 | 6.7900 |
06.pdh | CKIWMSYTPQ | 2 | -1 | 6.90592 | 7.12461 | 6.87802 | 11.6818 | 36.5734 | 5.3174 |
07.pdh | SMSLLPQLYV | 1 | -1 | 7.55705 | 10.0228 | 8.08338 | 8.61016 | 50.8402 | 6.2895 |
08.pdh | VDLISHCTNY | 2 | -2 | 9.21513 | 12.8115 | 6.02248 | 12.1875 | 80.1387 | 13.3066 |
09.pdh | APYIDLDDVI | 1 | -4 | 8.31937 | 12.9252 | 3.37628 | 7.84797 | 71.436 | 21.1582 |
10.pdh | IPVEQYSHKI | 3 | -2 | 7.9604 | 10.3516 | 4.03182 | 9.52373 | 60.5853 | 15.0268 |
11.pdh | FGLGFGNHQS | 2 | -1 | 7.09994 | 7.10804 | 3.55309 | 12.6514 | 41.7148 | 11.7404 |
12.pdh | VQITQDVIFR | 2 | -2 | 8.70787 | 11.9557 | 5.25832 | 8.60649 | 74.4139 | 14.1516 |
13.pdh | RNYFRPEDNE | 3 | -4 | 8.44553 | 6.00883 | 8.38621 | 27.9233 | 41.4656 | 4.9445 |
14.pdh | LIAMPFISDA | 1 | -2 | 7.32482 | 7.09511 | 4.89306 | 12.7505 | 43.499 | 8.8899 |
15.pdh | DYHDINFFMK | 3 | -3 | 6.6087 | 5.37234 | 5.82924 | 14.8909 | 28.6026 | 4.9067 |
16.pdh | NGNMCWRTAG | 2 | -1 | 7.75978 | 11.5624 | 7.2114 | 7.42933 | 57.6308 | 7.9916 |
17.pdh | YCFYEKENFH | 3 | -3 | 7.54838 | 7.38951 | 6.32253 | 10.937 | 47.4657 | 7.5074 |
18.pdh | EHVCKGPFLS | 3 | -2 | 8.6758 | 15.0929 | 4.7563 | 6.66874 | 79.5166 | 16.7182 |
19.pdh | PPPTRQRLLH | 4 | -1 | 7.37071 | 7.47655 | 3.35442 | 20.5102 | 38.2827 | 11.4126 |
固有値とコンフォメーションについて
細長い分子では,固有値2 / 固有値0 が大きくなる.逆に丸い分子では 1 に近くなる.
今回の実験で最大となったのは,09.pdh であった.分子の正電荷は 1 個で,負電荷は 4 個である.
逆に,固有値2 / 固有値0 が最小となった 15.pdh の正電荷は 3,負電荷は -3 である.
画像下側の分子は,09.pdh の,上側は 15.pdh の棒球モデルである.
09.pdh は負電荷の反発力で分子が伸びているようだ.
一方,15.pdh は正電荷と負電荷の引力で分子が折りたたまれているようだ.
ただし,15pdh は,すべての正電荷と負電荷がペアを形成しているわけではないようである.
これらの結果は,電荷の相互作用がコンフォメーションに影響していることと矛盾しない.
分子の概形は,ここで計算している固有値で表現できそうである.
シミュレーション実験で,コンフォメーションに電荷が大きく影響しているであろうことが確認できた.
今後,さらに研究を進めたい.
参考書の検索
- Amazon の「科学・テクノロジー」本カテゴリーでの,「タンパク質 構造」の検索結果です
- Amazon の「科学・テクノロジー」本カテゴリーでの,「生化学」の検索結果です