BMutate
XYZ 形式,PDB 形式,あるいは独自形式である BCL 形式のファイルを開いて,ペプチドに変異を 1 回導入します.
配列上の位置とアミノ酸を指定,乱数を発生させて位置とアミノ酸を設定,乱数に BLOSUM62 で重み付けをして位置とアミノ酸を設定するというオプションがあります.
変異したペプチドは(複数のペプチドが含まれる場合は変異したペプチドのみ),libopenbabel でエネルギー極小化してから出力します.
GUI 版分子モデリングソフト Builcule と比べて,多様な条件で変異を導入できます.
例えば,BLOSUM62 を利用した重み付けは,手持ちの材料で完全にランダムな変異よりも現実を反映した処理ができないか考えてみた試論です.
これを第 1 版としておきます.
ビルドとインストール
ソースファイル
- libbuilcule-2:ライブラリのページへのリンクです.先にライブラリをインストールしてください
- bmutate-1.2.0.tar.gz
ビルドとインストール
ダウンロードした bmutate-1.2.0.tar.gz を解凍し,生成するディレクトリに移動します.
~$ tar -zxvf bmutate-1.2.0.tar.gz
~$ cd bmutate-1.2.0
Makefile 作成ツールには,Automake を使っています.Automake の常法どおり,
~/bmutate-1.2.0$ ./configure
~/bmutate-1.2.0$ make
~/bmutate-1.2.0$ sudo make install
とすれば,/usr/local/bin に bmutate がインストールされます.
アンインストール
実行ファイルがインストールされるだけなので,これを削除してください.
~/$ sudo rm /usr/local/bin/bmutate
不具合
入力ファイルに複数の分子が含まれており,これらが空間的に近い位置にある場合,嵩高いアミノ酸への変異により分子間の共有結合が生成する場合があります.
これは,原子間の距離で共有結合を検知しているためです.
改定履歴
2023 年 10 月 23 日 : 1.2.0 リリース.変異の位置を乱数で,アミノ酸の種類を引数で指定するオプションを追加しました.
2023 年 10 月 19 日 : 1.1.0 リリース.ファイルに複数の分子が含まれる場合に,変異した分子しか出力されないバグを修正しました.
2023 年 10 月 19 日 : 1.0.0 リリース.
マニュアル
分子ファイルを読んで,検知できたペプチド中のアミノ酸を,下記の方法で 1 個変異します.
アミノ酸は一文字記号しか使えません.大文字小文字の区別しません.
[1] オプションで分子,ペプチド,配列番号,アミノ酸を指定
[2] 乱数で,分子,ペプチド,配列番号,アミノ酸を設定
[3] 乱数で,分子,ペプチド,配列番号.オプションでアミノ酸を指定
[4] 乱数を BLOSUM62 で重み付けし,分子,ペプチド,配列番号,アミノ酸を設定
ファイル形式
- 入力ファイルは,XYZ 形式,PDB 形式,および 独自形式である BCL 形式に対応しています
- 出力ファイルの形式は,OpenBabel の PDB 形式としています
- ファイル形式は,拡張子で判断します
オプション
ショートオプション(1 ハイフン + 1 文字)のみ使えます.
bpeptgen i:o:sm:p:a:c:rbhv
- i 入力ファイルを指定してください.相対パスでも絶対パスでも OK です
- o 出力ファイルを指定してください.相対パスでも絶対パスでも OK です
- s 位置とアミノ酸を指定して変異を導入します
- r 乱数で,分子,ペプチド,配列番号を設定します,アミノ酸については,乱数で設定することもオプションで指定することも可能です
- b BLOSUM62 で重み付けした乱数で,分子,ペプチド,配列番号,アミノ酸を設定します
- m s オプション使用時に分子番号を指定してください
- p s オプション使用時にペプチド番号を指定してください
- a s オプション使用時にアミノ酸番号を指定してください
- c アミノ酸を指定するオプションです.s オプション使用時,または アミノ酸を指定して r オプションを使用する場合に利用します
- h ヘルプを表示します
- v バージョンを表示します
生成するペプチドについて
生成するペプチドは,OpenBabel でエネルギー極小化しています.
これは,側鎖どうしの空間的な重なり合いを排除するためです.
使用例
[1] オプションで分子,ペプチド,配列番号,アミノ酸を指定
bmutate -i 入力ファイル -o 出力力ファイル -s -m 分子番号 -p ペプチド番号 -a 配列番号 -c 変異後のアミノ酸(一文字記号)
例:
~$ bmutate -i/home/User/wk/input.pdb -o/home/User/wk/output.pdb -s -m0 -p0 -a1 -cq
とすると,/home/User/wk/input.pdb の第 0 番めの分子の第 0 番めのペプチドの第一番めのアミノ酸を Q に変異した /home/User/wk/output.pdb というファイルが生成します.
[2] 乱数で,分子,ペプチド,配列番号,アミノ酸を設定
bmutate -i 入力ファイル -o 出力力ファイル -r
例:
~$ bmutate -i/home/User/wk/input.pdb -o/home/User/wk/output.pdb -r
とすると,/home/User/wk/input.pdb のどれかのアミノ酸が変異した /home/User/wk/output.pdb というファイルが生成します.
[3] 乱数で,分子,ペプチド,配列番号.オプションでアミノ酸を指定
bmutate -i 入力ファイル -o 出力力ファイル -r -c 変異後のアミノ酸
例:
~$ bmutate -i/home/User/wk/input.pdb -o/home/User/wk/output.pdb -r -cq
とすると,/home/User/wk/input.pdb のどれかのアミノ酸が Q に変異した /home/User/wk/output.pdb というファイルが生成します.
入力されるペプチドがすべて c オプションで指定されるアミノ酸である場合は,エラーとしています.
[4] 乱数を BLOSUM62 で重み付けし,分子,ペプチド,配列番号,アミノ酸を設定
bmutate -i 入力ファイル -o 出力力ファイル -r -c 変異後のアミノ酸
例:
~$ bmutate -i/home/User/wk/input.pdb -o/home/User/wk/output.pdb -b
とすると,/home/User/wk/input.pdb のどれかのアミノ酸が変異した /home/User/wk/output.pdb というファイルが生成します.
変異は,乱数を BLOSUM62 で重み付けして決定されます.
BLOSUM62 を利用した重み付け
考え方
乱数で変異を起こすよりも多少なりとも現実を反映するように,BLOSUM62 行列を次のように変形してみました.
これを第 1 版としておきます.
- 行列の各スコアについて,2(スコア - 1)を計算してスコアをオッズ比に戻す
- オッズ比を,「アミノ酸 X」 から「アミノ酸 Y」 への変異しやすさ(変異性向)と解釈する.なお,コーディング上の都合,オッズ比を 100000 倍している.
- 同じアミノ酸どうしの変異性向は 0 とする.これはプログラム上の無駄を省くため
- 行の合計を,あるアミノ酸が他のアミノ酸に変異する変異性向と解釈する
- 行の合計の合計(== 行列の総計 == オッズ比の総計)を,アミノ酸全体の変異性向と解釈する
計算結果 1.あるアミノ酸の変異確率
行の合計を行列の総計で割った値を,あるアミノ酸の変異確率と解釈します.
言い換えると,ペプチドに一回変異が起こるとしたら,各アミノ酸が何 % の確率で変異するかということです.
その値を表に示します(数値は %).
Cys や Pro は低い値となりました.
C | S | T | P | A | G | N | D | E | Q | H | R | K | M | I | L | V | F | Y | W |
2.2 | 5.6 | 4.2 | 2.1 | 4.3 | 2.3 | 5.1 | 4.5 | 6.4 | 6.1 | 4.9 | 4.8 | 5.4 | 5.7 | 6.9 | 5.6 | 6.6 | 6.0 | 7.8 | 3.4 |
計算結果 2.Cys が他のアミノ酸に変異する際の変異確率
Cys を例にして考えてみます.
Cys 行の各変異性向を行の総和で割った値を,Cys の変異確率と解釈します.
これは,Cys が変異するとしたら,何 % の確率でそれぞれのアミノ酸に変異するかということです.
その値を表に示します(数値は %).
Cys が変異するとすれば,Ala に変異する確率が最大で,続いて Met,Ile,Leu,Val となりました.
C | S | T | P | A | G | N | D | E | Q | H | R | K | M | I | L | V | F | Y | W |
0.0 | 8.6 | 8.6 | 2.2 | 17.2 | 2.2 | 2.2 | 2.2 | 1.1 | 2.2 | 2.2 | 2.2 | 2.2 | 8.6 | 8.6 | 8.6 | 8.6 | 4.3 | 4.3 | 4.3 |
すべてのアミノ酸について計算してみると…
F → Y が 50.2% だとか,V → I が 45.7% だとか,なんとなくそれらしい確率になりました.
ただし,逆方向の変異は異なる値になっています.
- | C | S | T | P | A | G | N | D | E | Q | H | R | K | M | I | L | V | F | Y | W |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
C | 0.0 | 8.6 | 8.6 | 2.2 | 17.2 | 2.2 | 2.2 | 2.2 | 1.1 | 2.2 | 2.2 | 2.2 | 2.2 | 8.6 | 8.6 | 8.6 | 8.6 | 4.3 | 4.3 | 4.3 |
S | 3.4 | 0.0 | 13.4 | 3.4 | 13.4 | 6.7 | 13.4 | 6.7 | 6.7 | 6.7 | 3.4 | 3.4 | 6.7 | 3.4 | 1.7 | 1.7 | 1.7 | 1.7 | 1.7 | 0.8 |
T | 4.4 | 17.8 | 0.0 | 4.4 | 8.9 | 2.2 | 8.9 | 4.4 | 4.4 | 4.4 | 2.2 | 4.4 | 4.4 | 4.4 | 4.4 | 4.4 | 8.9 | 2.2 | 2.2 | 2.2 |
P | 2.2 | 8.9 | 8.9 | 0.0 | 8.9 | 4.4 | 4.4 | 8.9 | 8.9 | 8.9 | 4.4 | 4.4 | 8.9 | 4.4 | 2.2 | 2.2 | 4.4 | 1.1 | 2.2 | 1.1 |
A | 8.8 | 17.6 | 8.8 | 4.4 | 0.0 | 8.8 | 2.2 | 2.2 | 4.4 | 4.4 | 2.2 | 4.4 | 4.4 | 4.4 | 4.4 | 4.4 | 8.8 | 2.2 | 2.2 | 1.1 |
G | 2.0 | 16.0 | 4.0 | 4.0 | 16.0 | 0.0 | 16.0 | 8.0 | 4.0 | 4.0 | 4.0 | 4.0 | 4.0 | 2.0 | 1.0 | 1.0 | 2.0 | 2.0 | 2.0 | 4.0 |
N | 0.9 | 14.6 | 7.3 | 1.8 | 1.8 | 7.3 | 0.0 | 14.6 | 7.3 | 7.3 | 14.6 | 7.3 | 7.3 | 1.8 | 0.9 | 0.9 | 0.9 | 0.9 | 1.8 | 0.5 |
D | 1.1 | 8.4 | 4.2 | 4.2 | 2.1 | 4.2 | 16.8 | 0.0 | 33.7 | 8.4 | 4.2 | 2.1 | 4.2 | 1.1 | 1.1 | 0.5 | 1.1 | 1.1 | 1.1 | 0.5 |
E | 0.4 | 5.9 | 2.9 | 2.9 | 2.9 | 1.5 | 5.9 | 23.4 | 0.0 | 23.4 | 5.9 | 5.9 | 11.7 | 1.5 | 0.7 | 0.7 | 1.5 | 0.7 | 1.5 | 0.7 |
Q | 0.8 | 6.1 | 3.1 | 3.1 | 3.1 | 1.5 | 6.1 | 6.1 | 24.4 | 0.0 | 6.1 | 12.2 | 12.2 | 6.1 | 0.8 | 1.5 | 1.5 | 0.8 | 3.1 | 1.5 |
H | 1.0 | 3.8 | 1.9 | 1.9 | 1.9 | 1.9 | 15.4 | 3.8 | 7.7 | 7.7 | 0.0 | 7.7 | 3.8 | 1.9 | 1.0 | 1.0 | 1.0 | 3.8 | 30.8 | 1.9 |
R | 1.0 | 3.9 | 3.9 | 1.9 | 3.9 | 1.9 | 7.8 | 1.9 | 7.8 | 15.5 | 7.8 | 0.0 | 31.1 | 3.9 | 1.0 | 1.9 | 1.0 | 1.0 | 1.9 | 1.0 |
K | 0.9 | 6.9 | 3.4 | 3.4 | 3.4 | 1.7 | 6.9 | 3.4 | 13.8 | 13.8 | 3.4 | 27.6 | 0.0 | 3.4 | 0.9 | 1.7 | 1.7 | 0.9 | 1.7 | 0.9 |
M | 3.3 | 3.3 | 3.3 | 1.6 | 3.3 | 0.8 | 1.6 | 0.8 | 1.6 | 6.6 | 1.6 | 3.3 | 3.3 | 0.0 | 13.1 | 26.2 | 13.1 | 6.6 | 3.3 | 3.3 |
I | 2.7 | 1.4 | 2.7 | 0.7 | 2.7 | 0.3 | 0.7 | 0.7 | 0.7 | 0.7 | 0.7 | 0.7 | 0.7 | 10.8 | 0.0 | 21.7 | 43.4 | 5.4 | 2.7 | 0.7 |
L | 3.4 | 1.7 | 3.4 | 0.8 | 3.4 | 0.4 | 0.8 | 0.4 | 0.8 | 1.7 | 0.8 | 1.7 | 1.7 | 26.9 | 26.9 | 0.0 | 13.4 | 6.7 | 3.4 | 1.7 |
V | 2.9 | 1.4 | 5.7 | 1.4 | 5.7 | 0.7 | 0.7 | 0.7 | 1.4 | 1.4 | 0.7 | 0.7 | 1.4 | 11.4 | 45.7 | 11.4 | 0.0 | 2.9 | 2.9 | 0.7 |
F | 1.6 | 1.6 | 1.6 | 0.4 | 1.6 | 0.8 | 0.8 | 0.8 | 0.8 | 0.8 | 3.1 | 0.8 | 0.8 | 6.3 | 6.3 | 6.3 | 3.1 | 0.0 | 50.2 | 12.5 |
Y | 1.2 | 1.2 | 1.2 | 0.6 | 1.2 | 0.6 | 1.2 | 0.6 | 1.2 | 2.4 | 19.2 | 1.2 | 1.2 | 2.4 | 2.4 | 2.4 | 2.4 | 38.3 | 0.0 | 19.2 |
W | 2.8 | 1.4 | 2.8 | 0.7 | 1.4 | 2.8 | 0.7 | 0.7 | 1.4 | 2.8 | 2.8 | 1.4 | 1.4 | 5.5 | 1.4 | 2.8 | 1.4 | 22.1 | 44.1 | 0.0 |