计算模拟,有说没有用,只能算作湿实验的锦上添花,有说他很有用,可以节省成本,缩短周期。然而我认为至少他可以做出许许多多美丽的图片,It is an art.下面就看看精美的序列比对图如何做的吧。
如这篇文章的多序列比对图:
SH Sternberg et al. Nature 1-4 (2015) doi:10.1038/nature15544
进行氨基酸序列比对时,10条以内序列用经典的Clustal比对,如果序列更多可采用MUSCLE或MAFFT,网址链接里也有在线版本,相比clustal优点是速度快,准确度高。比对完成后生成序列比对文件,复制或下载,命名为…….aln。在Espript网站提交aln文件,以及结构文件,比对图中可生成二级结构。特别注意一点网页允许弹窗。
提交pdb文件,可以在序列比对结果中生成相应的二级结构比对。所有设置完成后点击SUBMIT,会弹出窗口,得到序列比对图。特别注意一点网页允许弹窗。在线版本在断网情况下,就没法做了,推荐一个本地安装版Aline,在win 7 64下没能成功安装,主要是没有成功安装perl 的Tk模块,linux下比较容易安装,这部分请参考
http://mp.weixin.qq.com/s?__biz=MzIwODEzMjU3MQ==&mid=2666445566&idx=1&sn=2344f92a01d5458743e8e50c964d7766&scene=1&srcid=0906qWZzhDk9mM6ZbMYXZYi7#rd
Fig.2 Sequence alignment was prepared with the program Aline
http://bondxray.org/software/aline.html
如果研究进化关系,有成百上千条序列呢?怎么用完美的图形表示呢?如下如所示:
YuanS, Filipek S, Vogel H. A Gating Mechanismof the Serotonin 5-HT 3 Receptor[J]. Structure, 2016, 24(5): 816-825.
所采用的工具呢,就是WebLogo http://weblogo.threeplusone.com/
这段区域是高度保守的区域,根据结构分析,同样属于催化活性中心,可以把催化三联体之一的Ser40突变为Thr40,活性中心侧链多一个甲基会有什么影响?可以定点突变试一下。
来源:公众号(生物研究生)
棒!