<p>It sounds like <a href="https://blast.ncbi.nlm.nih.gov/smartblast/" rel="noreferrer">Smart BLAST</a> might do what you want. Here's the description of it's goal:</p>

<blockquote>
  <p>SmartBLAST is a new and experimental NCBI tool that makes it easier to
  complete common sequence analysis tasks, such as finding a candidate
  protein name for a sequence, locating regions of high sequence
  conservation, or identifying regions covered by database sequences but
  missing from the query.</p>
  
  <p>To do this, SmartBLAST performs the following tasks in much less time
  than it takes to run a typical BLASTp search:</p>
  
  <ul>
  <li>a BLASTp comparison of the query with the closest matching sequences available;</li>
  <li>a parallel BLASTp search to find the closest matches to high quality sequences from model organisms;</li>
  <li>a multiple alignment between the query and five of the closest matching sequences (usually including two high quality sequences);</li>
  <li>an analysis that produces a phylogenetic tree from the multiple sequence alignment.</li>
  </ul>
</blockquote>

<p>[from <a href="https://ncbiinsights.ncbi.nlm.nih.gov/2015/07/29/smartblast/" rel="noreferrer">NCBI Insights</a>]</p>

CC BY-SA 3.0
2017-06-19T11:02:35.177Z
SmartBLAST is a very good idea as it gives so much more divergent organism, but is there any way to give its result to PSI-BLAST for more iterations?
Sara
CC BY-SA 3.0
2017-06-19T11:25:42.060Z
@Sara just note the accession numbers of the sequences found by Smart BLAST and make sure to add them to your alignment in the psi-blast. They should be found as hits there too.
terdon
CC BY-SA 3.0
2017-06-19T04:39:10.010Z
2017-06-19T04:39:10.010Z
gringer
<p>It depends on what you're looking for. If you're just looking for sequence homology, then you can simply pick the best hits from a blast search. If, however, you are referring to <em>functional</em> homology, if you are looking for the protein which has the same functions as your query, then it's more complicated. </p>

<p>Sequence homology is not enough to infer functional homology. For example, you can have cases of gene duplication and subsequent functional divergence. Such <em>paralogs</em> are still <em>homologs</em> (<a href="https://biology.stackexchange.com/a/4964/1306">paralogs are a subset of homologs</a>), but they don't necessarily have the same function. It is also often the case that the homolog (be it orthologous or paralogous) of a protein in species B has a completely different function than its homolog in species A despite a high level of sequence similarity. This is usually very hard to determine <em>in silico</em>.</p>

<p>To find the functionally homologous protein(s), you would ideally need to identify the essential residues that allow your protein to perform its function. This could be done using something like <a href="http://pfam.xfam.org/" rel="noreferrer">PFam</a> which will identify protein domains. You can then check whether the homologs you find also have this domain. </p>

<p>This is essentially what PSI-blast does. Although it doesn't take domains into account, each successive iteration is used to build a model of your proptein. The model is built under the reasonable assumption that highly conserved residues are important. So it will consider more diverged sequences as homologous if those residues are conserved. </p>

<p>If you know how your protein works and what residues are important, you can use that knowledge to refine the results of your PSI-blast. If you don't, you'll have to use only "good" hits to make the model. One way to do this, for well studied proteins, is to only add proteins that are already annotated as homologs of what you are looking for to build your model, then use that model to search in un-annotated species. </p>

<p>If you don't know, you could try looking for recognizable protein domains in your query protein (use PFam) and then use the HMM (hidden markov model) of the domain to identify important residues. For example, <a href="http://pfam.xfam.org/family/PF01754#tabview=tab4" rel="noreferrer">this is</a> the HMM logo for the zf-A20 zinc finger domain:</p>

<p><a href="https://i.sstatic.net/90nPc.png" rel="noreferrer"><img src="https://i.sstatic.net/90nPc.png" alt="zinc finger sequence logo"></a></p>

<p>The huge cysteine (C) residues are shown that size because they are very conserved across proteins carrying this domain and, presumably, are functionally important for the domain. So, if you pass your protein through PFam and identify domains, find the important residues and make sure all your homologs have those conserved. If using PSI-blast, only include sequences where those residues are conserved in the results you keep.</p>

<p>Finally, another useful tool that works in the same way is <a href="http://hmmer.org/" rel="noreferrer">HMMER</a>. This takes a protein alignment as input, like PSI-blast builds an HMM model from it and then can use this model to query a protein database for more hits. Methods like HMMER and PSI-blast are far better than simple sequence similarity approaches when looking for homologs.  </p>

CC BY-SA 3.0
2017-06-19T14:02:38.123Z
Great tool HMMER, is as fast as BLAST and works pretty well
llrs
CC BY-SA 3.0
2017-06-19T16:27:44.687Z
Dear terdon, your answer was the best answer someone could have given to this question and it was exactly the answer I was looking for. Thank you very much for your excellent answer
Sara
CC BY-SA 3.0
2017-06-19T16:28:09.320Z
@Sara you're very welcome :)
terdon
CC BY-SA 3.0
2017-06-19T13:54:59.920Z
2017-06-19T13:54:59.920Z
terdon
<p>I have a refseq ID of a protein from E.coli and I want to find homologs of this protein. I ran Blast against refseq database but I got a lot of sequences most of which were from Ecoli again. I decided to run PSI-Blast to get more divergent species, but I do not exactly know if my result are real homologs or false positives. what is your idea for finding homologous protein sequences from more divergent species? And what can I do for selecting the real (but not false positive) hits?</p>

CC BY-SA 3.0
2017-06-19T09:28:51.850Z
False positive can be judged by the E and P-values so remember to use them.
Kristoffer Vitting-Seerup
CC BY-SA 3.0
2017-06-19T10:57:19.390Z
I would suggest that if you can robustly obtain the same topology using a variety of phylogenetic methods and models on your proteins, there is a chance that the homology is true. Unfortunately, this might be difficult to obtain when dealing with quite distant species, even if the homology is true.
bli
<sequence-homology>
Finding homologs of a protein sequence
# Finding homologs of a protein sequence

(Asked by: username_0 on 2017-06-19)

I have a refseq ID of a protein from E.coli and I want to find homologs of this protein. I ran Blast against refseq database but I got a lot of sequences most of which were from Ecoli again. I decided to run PSI-Blast to get more divergent species, but I do not exactly know if my result are real homologs or false positives. what is your idea for finding homologous protein sequences from more divergent species? And what can I do for selecting the real (but not false positive) hits?

**Tags:** sequence-homology
> **Comments:**
> **username_1**: False positive can be judged by the E and P-values so remember to use them.  
> **username_2**: I would suggest that if you can robustly obtain the same topology using a variety of phylogenetic methods and models on your proteins, there is a chance that the homology is true. Unfortunately, this might be difficult to obtain when dealing with quite distant species, even if the homology is true.  

---

## 2 answers

### Answer by username_4 on 2017-06-19. Score: 5 | Accepted answer

It depends on what you're looking for. If you're just looking for sequence homology, then you can simply pick the best hits from a blast search. If, however, you are referring to *functional* homology, if you are looking for the protein which has the same functions as your query, then it's more complicated.

Sequence homology is not enough to infer functional homology. For example, you can have cases of gene duplication and subsequent functional divergence. Such *paralogs* are still *homologs* ([paralogs are a subset of homologs](https://biology.stackexchange.com/a/4964/1306)), but they don't necessarily have the same function. It is also often the case that the homolog (be it orthologous or paralogous) of a protein in species B has a completely different function than its homolog in species A despite a high level of sequence similarity. This is usually very hard to determine *in silico*.

To find the functionally homologous protein(s), you would ideally need to identify the essential residues that allow your protein to perform its function. This could be done using something like [PFam](http://pfam.xfam.org/) which will identify protein domains. You can then check whether the homologs you find also have this domain.

This is essentially what PSI-blast does. Although it doesn't take domains into account, each successive iteration is used to build a model of your proptein. The model is built under the reasonable assumption that highly conserved residues are important. So it will consider more diverged sequences as homologous if those residues are conserved.

If you know how your protein works and what residues are important, you can use that knowledge to refine the results of your PSI-blast. If you don't, you'll have to use only "good" hits to make the model. One way to do this, for well studied proteins, is to only add proteins that are already annotated as homologs of what you are looking for to build your model, then use that model to search in un-annotated species.

If you don't know, you could try looking for recognizable protein domains in your query protein (use PFam) and then use the HMM (hidden markov model) of the domain to identify important residues. For example, [this is](http://pfam.xfam.org/family/PF01754#tabview=tab4) the HMM logo for the zf-A20 zinc finger domain:

[![zinc finger sequence logo](https://i.sstatic.net/90nPc.png)](https://i.sstatic.net/90nPc.png)

The huge cysteine (C) residues are shown that size because they are very conserved across proteins carrying this domain and, presumably, are functionally important for the domain. So, if you pass your protein through PFam and identify domains, find the important residues and make sure all your homologs have those conserved. If using PSI-blast, only include sequences where those residues are conserved in the results you keep.

Finally, another useful tool that works in the same way is [HMMER](http://hmmer.org/). This takes a protein alignment as input, like PSI-blast builds an HMM model from it and then can use this model to query a protein database for more hits. Methods like HMMER and PSI-blast are far better than simple sequence similarity approaches when looking for homologs.

> **Comments:**
> **username_5**: Great tool HMMER, is as fast as BLAST and works pretty well  
> **username_0**: Dear terdon, your answer was the best answer someone could have given to this question and it was exactly the answer I was looking for. Thank you very much for your excellent answer  
> **username_4**: @Sara you're very welcome :)  

---

### Answer by username_3 on 2017-06-19. Score: 5

It sounds like [Smart BLAST](https://blast.ncbi.nlm.nih.gov/smartblast/) might do what you want. Here's the description of it's goal:

> SmartBLAST is a new and experimental NCBI tool that makes it easier to
> complete common sequence analysis tasks, such as finding a candidate
> protein name for a sequence, locating regions of high sequence
> conservation, or identifying regions covered by database sequences but
> missing from the query.
>
> To do this, SmartBLAST performs the following tasks in much less time
> than it takes to run a typical BLASTp search:
>
> * a BLASTp comparison of the query with the closest matching sequences available;
> * a parallel BLASTp search to find the closest matches to high quality sequences from model organisms;
> * a multiple alignment between the query and five of the closest matching sequences (usually including two high quality sequences);
> * an analysis that produces a phylogenetic tree from the multiple sequence alignment.

[from [NCBI Insights](https://ncbiinsights.ncbi.nlm.nih.gov/2015/07/29/smartblast/)]

> **Comments:**
> **username_0**: SmartBLAST is a very good idea as it gives so much more divergent organism, but is there any way to give its result to PSI-BLAST for more iterations?  
> **username_4**: @Sara just note the accession numbers of the sequences found by Smart BLAST and make sure to add them to your alignment in the psi-blast. They should be found as hits there too.  

---

