BLAST Basic Local Alignment Search Tool BLAST BLAST

BLAST Basic Local Alignment Search Tool

BLAST החכה BLAST (Basic Local Alignment Search Tool) allows rapid sequence comparison of a query sequence [[ (רצף שאילתא nucleotides or amino acids) הפיתיון בחכה against a database הים הגדול. לצורך דיג מוצלח פיתיון ומקווה מים בהתאם , יש לבחור חכה . לשאלה הביולוגית

Comparing the query sequence to known sequences in databases is fundamental to understanding the relatedness of any query sequence to other known proteins or DNA sequences. Applications include: • Identifying shared similarities with sequences already deposited in the databanks (orthologs and paralogs? ) • Discovering new genes or proteins (ascertaining existence of a putative ORF) • Discovering variants of genes or proteins • Identifying functional motifs shared with other proteins. • Investigating expressed sequence tags (ESTs) • Exploring protein structure and function

Why use local alignment for database searches? Local alignment is a useful approach to DB searching because many query sequences have domains, active sites or other motifs that have local but not global regions of similarity to other sequences.

BLAST (1) for the query, find the list of high scoring words of length w Query Sequence of length L For each word from the query sequence find the list of words that will score at least T when scored using a pair-score matrix (e. g. PAM 250, BLOSUM)

BLAST (cont. ) (2) Compare the word list to the database and identify exact matches database sequence Word List Exact matches of words from word lists (3) For each word match, extend the alignment in both directions to find alignments that score greater than a threshold of value S maximal segment pairs (MSPs)

Raw score "S" of the alignment is usually calculated by summing the scores for matches, mismatches and gaps in the alignment. Normalized score (bits) - bit scores from different alignments, even those employing different scoring matrices can be compared. The higher the score the better the alignment, BUT the significance of an alignment can not be deduced from the score alone.

E-value (Expectation value) • Expect value of 10 for a match means, in a database of current size, one might expect to see 10 matches with a similar or better score, simply by chance alone • E-value is the most commonly used threshold in database searches. Only those hits with E-values smaller than the set threshold will be reported in the output • Increasing the E-value enables you to see biologically related sequences but statistically insignificant

To evaluate the alignment • Examine statistical parameters: �� Normalized score �� E value �� % identity �� % similarity �� % gaps • Examine the alignment itself. • Use biological common sense. Don’t rely on statistical significance!!!

• Limit DB • Limit organism ספירת האפשרויות השונות • Filter reported entries by keyword • (Limit to a specific domain) • Change matrix and/or gap penalties • Change E-value • Add filter for low complexity

What can we do if there are hardly any matches?

• Check choice of DB • Check choice of organism • Remove filter for low complexity • Change matrix or gap penalties • Increase E-value

DNA vs. Protein searches If we have a nucleotide sequence, should we search the DNA databases only? Or should we translate it to protein and search protein databases? Query: DNA Protein Database: DNA Protein Translating causes loss of information but protein sequence is more conserved than DNA sequence It is therefore advisable to translate a nucleotide sequence to protein and search protein databases for homology

Why use a nucleotide sequence after all? n No ORF found. No similar protein sequences were found Specific DNA databases are available (EST) n To find duplicated genes in a genome n To find pseudogenes n To find the location of non-protein coding genes n n in the genome (si. RNA etc. )

Blast flavors Query: DB: n n n DNA Protein Blast. N - nt versus nt database Blast. P - protein versus protein database Blast. X - translated nt (6 frames) versus protein database t. Blast. N - protein versus translated nt database (6 frames) t. Blast. X - translated nt versus translated nt database (both 6 frames)

Uses of BLAST programs BLASTx – compares a nucleotide query seq translated in all reading frames against a prot seq db. DNA protein If you have a DNA seq and you want to now what protein (if any) it encodes, you can perform BLASTx search.

t. BLASTn – compares a protein query seq against a nucleotide seq db which is translated in all reading frames. Protein DNA You can use this program to ask whether a DNA or ESTs db contains a nuc seq encoding a protein that matches your protein of interest.

t. BLASTx – translates DNA from query and compares it to db of DNA seqs all translated to all reading frames DNA (nr db cannot be used, because it’s too large) Used to determine whether an entire DNA db contains genes that encodes proteins similar to your query. (If blastx or tblastn fail)

E-value