blast l.
Download
Skip this Video
Download Presentation
BLAST

Loading in 2 Seconds...

play fullscreen
1 / 48

BLAST - PowerPoint PPT Presentation


  • 170 Views
  • Uploaded on

BLAST. BLAST. B asic L ocal A lignment S earch T ool Developed in 1990 and 1997 (S. Altschul) A heuristic method for performing local alignments through searches of high scoring segment pairs (HSP’s) 1st to use statistics to predict significance of initial matches - saves on false leads

loader
I am the owner, or an agent authorized to act on behalf of the owner, of the copyrighted work described.
capcha
Download Presentation

PowerPoint Slideshow about 'BLAST' - morey


Download Now An Image/Link below is provided (as is) to download presentation

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.


- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -
Presentation Transcript
blast2
BLAST
  • Basic Local Alignment Search Tool
  • Developed in 1990 and 1997 (S. Altschul)
  • A heuristic method for performing local alignments through searches of high scoring segment pairs (HSP’s)
  • 1st to use statistics to predict significance of initial matches - saves on false leads
  • Offers both sensitivity and speed
blast3
BLAST
  • Looks for clusters of nearby or locally dense “similar or homologous” k-tuples
  • Uses “look-up” tables to shorten search time
  • Uses larger “word size” than FASTA to accelerate the search process
  • Performs both Global and Local alignment
  • Fastest and most frequently used sequence alignment tool -- THE STANDARD
blast access
BLAST Access
  • NCBI BLAST
  • http://www.ncbi.nlm.nih.gov/BLAST/
  • Canadian Bioinformatics Resource BLAST
  • http://cbr-rbc.nrc-cnrc.gc.ca/blast/
  • European Bioinformatics Institute BLAST
  • http://www.ebi.ac.uk/blastall/
  • http://www.ebi.ac.uk/blast2/
different flavours of blast
Different Flavours of BLAST
  • BLASTP - protein query against protein DB
  • BLASTN - DNA/RNA query against GenBank (DNA)
  • BLASTX - 6 frame trans. DNA query against proteinDB
  • TBLASTN - protein query against 6 frame GB transl.
  • TBLASTX - 6 frame DNA query to 6 frame GB transl.
  • PSI-BLAST - protein ‘profile’ query against protein DB
  • PHI-BLAST - protein pattern against protein DB
other blast services
Other BLAST Services
  • MEGABLAST - for comparison of large sets of long DNA sequences
  • RPS-BLAST - Conserved Domain Detection
  • BLAST 2 Sequences - for performing pairwise alignments for 2 chosen sequences
  • Genomic BLAST - for alignments against select human, microbial or malarial genomes
  • VecScreen - for detecting cloning vector contamination in sequenced data
mt0895
MT0895
  • MMKIQIYGTGCANCQMLEKNAREAVKELGIDAEFEKIKEMDQILEAGLTALPGLAVDGELKIMGRVASKEEIKKILS
running ncbi blast12
Running NCBI BLAST
  • Paste in sequence (FASTA format, raw sequence or type in GI or accession number)

>Mysequence MT0895 KIQIYGTGCANCQMLEKNAREAVKELGIDAEFEKIKEMDQILEAGLTALPGLAVDGELKIDS

OR

> KIQIYGTGCANCQMLEKNAREAVKELGIDAEFEKIKEMDQILEAGLTALPGLAVDGELKIDS

OR

KIQIYGTGCANCQMLEKNAREAVKELGIDAEFEKIKEMDQILEAGLTALPGLAVDGELKIDS

running ncbi blast13
Running NCBI BLAST
  • Choose a range of interest in the sequence “set subsequences” (not usually used)
  • Select the database from pull-down menu (usually choose nr = non-redundant)
  • Keep CD Search “check box” on
  • Leave “Options” unchanged (use defaults)
  • Go to “Format” menu and adjust Number of descriptions and alignments as desired
running ncbi blast14
Running NCBI BLAST

Select Database

conserved domain database
Conserved Domain Database
  • Contains a collection of pre-identified functional or structural domains
  • Derived from Pfam and Smart databases as well as other sources
  • Uses Reverse Position Specific BLAST (RPS-BLAST) to perform search
  • Query sequence is compared to a PSSM derived from each of the aligned domains
running ncbi blast16
Running NCBI BLAST

Click BLAST!

blast parameters
BLAST Parameters
  • Identities - No. & % exact residue matches
  • Positives - No. and % similar & ID matches
  • Gaps - No. & % gaps introduced
  • Score - Summed HSP score (S)
  • Bit Score - a normalized score (S’)
  • Expect (E) - Expected # of chance HSP aligns
  • P - Probability of getting a score > X
  • T - Minimum word or k-tuple score (Threshold)
blast rules of thumb
BLAST - Rules of Thumb
  • Expect (E-value) is equal to the number of BLAST alignments with a given Score that are expected to be seen simply due to chance
  • Don’t trust a BLAST alignment with an Expect score > 0.01 (Grey zone is between 0.01 - 1)
  • Expect and Score are related, but Expect contains more information. Note that %Identies is more useful than the bit Score
  • Recall Doolittle’s Curve (%ID vs. Length, next slide) %ID > 30 - numres/50
  • If uncertain about a hit, perform a PSI-BLAST search
doolittle s curve
Doolittle’s Curve

Twilight Zone

blast options30
BLAST Options
  • Composition-based statistics (Yes)
  • Sequence Complexity Filter (Yes)
  • Expect (E) value (10)
  • Word Size (3)
  • Substitution or Scoring Matrix (Blosum62)
  • Gap Insertion Penalty (11)
  • Gap Extension Penalty (1)
composition statistics
Composition Statistics
  • Recent addition to BLAST algorithm
  • Permits calculated E (Expect) values to account for amino acid composition of queries and database hits
  • Improves accuracy and reduces false positives
  • Effectively conducts a different scoring procedure for each sequence in database
lcr s low complexity
LCR’s (low complexity)
  • Watch out for…
    • transmembrane or signal peptide regions
    • coil-coil regions
    • short amino acid repeats (collagen, elastin)
    • homopolymeric repeats
  • BLAST uses SEG to mask amino acids
  • BLAST uses DUST to mask bases
scoring matrices
Scoring Matrices
  • BLOSUM Matrices
    • Developed by Henikoff & Henikoff (1992)
    • BLOcks SUbstitution Matrix
    • Derived from the BLOCKS database
  • PAM Matrices
    • Developed by Schwarz and Dayhoff (1978)
    • Point Accepted Mutation
    • Derived from manual alignments of closely related proteins
how to make your own matrix
How to Make Your Own Matrix

A

C D ...

ACDEFGH..

ACDEFGK..

AADEFGH..

GCDEFGH..

ACAEYGK..

ACAEFAH..

#Aobs

f(A,A) =

A

0.8 -- --

#Aexp

C

0.2 0.8 --

D

0.0 0.3 1.0

#C/Aobs

f(C,A) =

E

-- -- --

#Aexp

+

#Cexp

Perform Calculate Fill Sub

Alignment Frequencies Matrix

pam versus blosum
First useful scoring matrix for protein

Assumed a Markov Model of evolution (I.e. all sites equally mutable and independent)

Derived from small, closely related proteins with ~15% divergence

Much later entry to matrix “sweepstakes”

No evolutionary model is assumed

Built from PROSITE derived sequence blocks

Uses much larger, more diverse set of protein sequences (30% - 90% ID)

PAM versus BLOSUM
pam versus blosum36
Higher PAM numbers to detect more remote sequence similarities

Lower PAM numbers to detect high similarities

1 PAM ~ 1 million years of divergence

Errors in PAM 1 are scaled 250X in PAM 250

Lower BLOSUM numbers to detect more remote sequence similarities

Higher BLOSUM numbers to detect high similarities

Sensitive to structural and functional subsitution

Errors in BLOSUM arise from errors in alignment

PAM versus BLOSUM
pam matricies
PAM Matricies
  • PAM 40 - prepared by multiplying PAM 1 by itself a total of 40 times best for short alignments with high similarity
  • PAM 120 - prepared by multiplying PAM 1 by itself a total of 120 times best for general alignment
  • PAM 250 - prepared by multiplying PAM 1 by itself a total of 250 times best for detecting distant sequence similarity
blosum matricies
BLOSUM Matricies
  • BLOSUM 90 - prepared from BLOCKS sequences with >90% sequence ID best for short alignments with high similarity
  • BLOSUM 62 - prepared from BLOCKS sequences with >62% sequence ID best for general alignment (default)
  • BLOSUM 30 - prepared from BLOCKS sequences with >30% sequence ID best for detecting weak local alignments
psi blast algorithm
PSI-BLAST Algorithm
  • Perform initial alignment with BLAST using BLOSUM 62 substitution matrix
  • Construct a multiple alignment from matches
  • Prepare position specific scoring matrix
  • Use PSSM profile as the scoring matrix for a second BLAST run against database
  • Repeat steps 3-5 until convergence
psi blast43
PSI-BLAST

PresSIterate!

psi blast45
PSI-BLAST
  • For Protein Sequences ONLY
  • Much more sensitive than BLAST
  • Slower (iterative process)
  • Often yields results that are as good as many common threading methods
  • SHOULD BE YOUR FIRST CHOICE IN ANALYZING A NEW SEQUENCE
still confused
Still Confused?

http://www.ncbi.nlm.nih.gov/Education/BLASTinfo/information3.html

conclusions
Conclusions
  • BLAST is the most important program in bioinformatics (maybe all of biology)
  • BLAST is based on sound statistical principles (key to its speed and sensitivity)
  • A basic understanding of its principles is key for using/interpreting BLAST output
  • Use NBLAST or MEGABLAST for DNA
  • Use PSI-BLAST for protein searches
ad