Question: Building Hidden Markov Model (HMM) for proteins
0
gravatar for saadleeshehreen
11 days ago by
saadleeshehreen30 wrote:

Hi,

I have 460 amino acid sequences for a specific protein family X. I want to make a HMM model for those sequences. I will use that HMM to search homologs of X / neighbour of X in some bacterial genomes. How can I make a HMM model ? Is there any software available? What is the step to do that?

Thanks in advance.

sequence • 107 views
ADD COMMENTlink modified 11 days ago by toheitka160 • written 11 days ago by saadleeshehreen30
1
gravatar for toheitka
11 days ago by
toheitka160
Germany/Dresden
toheitka160 wrote:

The HMMER software (which is well documented) can be used to produce HMMs from alignments.

As for searching DNA with protein HMMs:

  • nhmmer can be used to build nHMMs from nucleotide alignments. They can be used to query DNA.
  • hmmer can be used to build HMMs from amino acid alignments. They can be used to query proteins.

The current HMMER manual (2015, version 3.1) states: "Still missing: Translated comparisons. We’d of course love to have the HMM equivalents of BLASTX, TBLASTN, and TBLASTX. They’ll come."

(Of course, you could always translate your genome in all frames , chop it up, and then screen it using protein HMMs. It is a bit ugly, you might run into some frameshift trouble, but maybe it works?)

ADD COMMENTlink written 11 days ago by toheitka160

To give a bit more details, to build a HMM for a set of proteins, the steps are:
- build a multiple sequence alignment with e.g. Clustal
- run hmmbuild (from the HMMER package) with the multiple sequence alignment as input

ADD REPLYlink modified 11 days ago • written 11 days ago by Jean-Karim Heriche15k

Thanks. From their manual I came to know, I need to give .sto file as input to get a .hmm file. Does Clustal give .sto file as output? Or I have to use different software to convert my file after clustal?

Cheers

ADD REPLYlink written 10 days ago by saadleeshehreen30

As far as I remember, hmmbuild can read alignments in several formats. Just check the docs.

ADD REPLYlink written 8 days ago by Jean-Karim Heriche15k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 727 users visited in the last hour