Hi,
I have a set of proteins and I need to search homologous partners for them. I wanted to automate the process of searching and I wrote a perl script for that. Now, the question is what should be the % of identity (both min and max) should I use in searching the homologous sequences. Also, the database that I using for blastp is PDB , since I need the structures only .. Kindly help out with this problem ??
Regards, BHARAT
•
written
by
bharat.85.monu •
1 answer
I recalls that 30% is an empirical cutoff in term of protein sequence similarity.
If you use BLAST, then E-value serves as a better indicator of homology, comparing to identity. Because E-value takes into account the lengths of query and subject sequences. For example, a short protein is more likely to be somehow similar to another random guy simply by chance, in which case, a high E-value speaks stronger than a high identity. As I know there's no standard cutoff for E-value. You can try from 1e-2 (IMG's protocol) to 1e-10 or even lower.
On the other hand, there are a bunch of programs to help you identify orthologs, using more sophisticated algorithms, such as OrthoMCL. You can try those...
Log in to answer this question.
-
Finding homologous sequences that have crystal structures
written by Hansen_869 •
-
Identifying genes from contigs (obtained from transcriptome of a fungi)
written by bharat.85.monu •
-
Quantifying Structural Conservation From Msa
written by bharat.85.monu •
-
Finding Structures Of Protein Complexes!
written by kandoigaurav
-
N-Cap And C-Cap Residues In Strands In Beta-Hairpins
written by bharat.85.monu •
-
How To Study Protein Structure Evolution
written by bharat.85.monu •
-
Tool For Evolutionary Conservation Of Protein Folds
written by bharat.85.monu •
-
Missing Pdb Enteries In Scop And Cath
written by bharat.85.monu •
-
Ncbi Standalone Blast -Outfmt
written by bharat.85.monu •
-
Generating All Interacting Partners Of A Protein
written by Reyhaneh
The PID that you choose varies for different purposes. Give us more details about your problem for answering your question.
My query is that I want to obtain homologous protein structures for my dataset ?? . Since the dataset is large I need some sort of cutoff values to identify the homologous sequences
Please take care with the word "homologous". It has a very specific meaning which you should understand.
You can use Dali to make the identity scores from structural alignments.