氨基酸保守性的着色

网上有vmd进行原子数值着色方法的教程，格式要求较为严格，譬如，原子序号与数值所有对应，数值只能是一列数据。其实，能够没必要要求这么严格的格式，只要有原子与数值对应关系便可。这里采用数值替换氨基酸的B-factor数据的例子，而后在chimera中采用render by attribute显示。
该例子中氨基酸的数值是用生物信息学方法计算获得的保守打分，采用的是SCA5（statistics coupling analysis）算法获得。计算结果由Matlab导出，最终转到文本可编辑的数据，数据格式是
序列号 TAB 数值
例如：
1  0.58
2  0.34
...
采用bash脚本处理
[code=Bash width=600px]
#!/bin/bash
#usage: chBfactor.sh cons_score.file pdb.file > new_pdb.file

#reading conservation from file
declare -a cons_score_data
while read residue cons_score; do
   cons_score_data[$residue]=$cons_score
done < $1

#reading pdb file
cat $2 | grep '^ATOM' > tmp.pdb
while read LINE; do
   res_number=${LINE:22:4}
   echo "${LINE:0:60}" ${cons_score_data[`echo $res_number`]}"${LINE:67:13}"
done < tmp.pdb

rm tmp.pdb
[/code]

pymol中也是能够替换b factor来数值着色，但是一样在pymol中替换b factor也须要通过处理。该脚本简单易用，输出的new_pdb.file能够在pymol, chimera, vmd等各个软件中从新着色。

分享http://bioms.org/forum.php?mod=viewthread&tid=970&highlight=%E6%B0%A8%E5%9F%BA%E9%85%B8%E6%95%B0%E5%80%BC%E7%9D%80%E8%89%B2