HMMER学习笔记
简单的讲HMMER就是一款比对软件，可以类似于blast等，但是它提供的比对结果要精确于blast,相应的速度也要慢。与其类似的比对软件有PSI-blast、SAM、PFTOOLS等。
用途一：对于一条未知的序列在蛋白质数据库中寻找比对，HMMER完全可以取代目前常用的BLASTP 和 PSI-BLAST。
用途二：可以自动注释蛋白质结构域。包括pfam、SMART。
用途三：在数据库中寻找已有家族的相似序列。
HMMER的程序如下，注意每个程序都是输出文件在前面，输入文件在后面：
hmmbuild （建立参考数据的隐马尔可科夫模型）
输入：多重序列比对的文件（Stockholm file）
输出：建立的这些多重序列比对的隐马尔可科夫模型，
例子： hmmbuild globins4.hmm tutorial/globins4.sto
hmmalign （多重序列比对）输入：输入是fasta文件输出：Stockholm 格式的多重序列比对文件例子：hmmalign globins4.hmm tutorial/globins45.fa
hmmsearch（在数据库中寻找已经建立好的模体）
输入：建立好的参考数据模型、要搜索的数据库
输出：输出结果
例子： hmmsearch globins4.hmm uniprot sprot.fasta > globins4.out
====================================================
phmmer（对于只有一条未知蛋白序列的情况下，就省去多重序列比对和建立模型的过程直接一步到位，期间使用了BLOSUM62 scores打分矩阵，类似于BLASTP-like）
输入:tutorial/HBB HUMAN 是你要检测的fasta格式的序列，uniprot sprot.fasta你要搜索的数据库
1 | phmmer tutorial/HBB HUMAN uniprot sprot.fasta |
nhmmer（类似于blastn，输入的query可以是fasta格式的DNA或者RNA序列,但是只能是一条。如果是多条请先做多重序列比对（hmmalign），然后使用建立模型(hmmbuild)，然后生成的hmm作为输入）
1 | nhmmer MADE1.hmm dna target.fa > MADE1.out |
=================================================
如果你有未知的序列去搜寻已知的模型数据库，例如：Pfam, SMART, or TIGRFams.这个时候搜寻的过程就使用hmmscan,相应的核酸序列就使用nhmmscan.
首先还是要建立模型数据库，如果你要合并多个可以直接:
1 | cat globins4.hmm fn3.hmm Pkinase.hmm > minifam |
如果你仅仅有的是模型数据库的多重序列比对格式文件例如来自pfam数据库的（Stockholm file）Pfam-A.seed文件：你可以利用hmmbuild 建立模型数据库：
1 | hmmbuild Pfam-A.hmm Pfam-A.see |
生成的文件为pfam-A.hmm由于你一般建立的模型数据库比较大，所以需要建立index,运行：hmmpress minifam（感觉有点像formatdb）
最后运行：
1 | mmscan minifam tutorial/7LESS DROME |
原文来自：http://blog.sina.com.cn/s/blog_83f77c940101rzlo.html
相关推荐：
- NCBI在线BLAST使用方法与结果详解 2941
- 神经网络术语：Epoch、Batch Size和迭代 527
- Consed的安装与使用教程 465
- 陈连福的NGS生物信息学培训教材V2.1 277
- WGCNA分析使用教程 272
-
-
***1Levin 2015-03-31 10:40#2
我想问，那个多重序列比对的文件，以sto命名的文件是怎么来的，怎么可以得到，如果用muscle的话，结果问价也不是sto这种格式的，怎么才能得到Stockholm 格式的多重序列比对文件
-
***2xuyk1993 2015-04-01 10:16#3
it’s from Userguide of Hmmer , maybe it can help you ,sorry ,i can’t used the chinese now. “Most popular alignment formats are similar block-based formats, and can be turned into Stockholm format with a little editing or scripting. Don’t forget the # STOCKHOLM 1.0 line at the start of the alignment, nor the //at the end. Stockholm alignments can be concatenated to create an alignment database flatfile containing many alignments” 出现问题中文输入法
-
***3Levin 2015-04-02 10:46#4
谢谢，原来他们的格式还真的很相似，我已经改过来了，这里我还有个问题，我用muscle软件做alignment，选用-msf参数（Write output in GCG MSF format）和-clw参数（Write output in CLUSTALW format）的结果有什么区别？因为我接下来要用HMMbuild，所以我不知道该选哪个参数？或者还有其他的alignment软件推荐使用的吗？
-
***2潇潇 2016-03-10 21:54#5
在Pfam里面可以下载sto格式的
最新创建圈子
-
原料药研发及国内外注册申报
2019-01-25 10:41圈主：caolianhui 帖子：33 -
制药工程交流
2019-01-25 10:40圈主：polysciences 帖子：30 -
健康管理
2019-01-25 10:40圈主：neuromics 帖子：20 -
发酵技术
2019-01-25 10:39圈主：fitzgerald 帖子：17 -
医学肿瘤学临床试验
2019-01-25 10:39圈主：bma 帖子：58
你有联系方式吗？我想学习这个软件