searching molecular databases with blast
Download
Skip this Video
Download Presentation
Searching Molecular Databases with BLAST

Loading in 2 Seconds...

play fullscreen
1 / 54

Searching Molecular Databases with BLAST - PowerPoint PPT Presentation


  • 296 Views
  • Uploaded on

Searching Molecular Databases with BLAST. Searching Molecular Databases with BLAST. B asic L ocal A lignment S earch T ool How BLAST works Interpreting search results The NCBI Web BLAST interface Demonstration and exercises. Why learn sequence database searching?. What have I cloned ?

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 'Searching Molecular Databases with BLAST' - MikeCarlo


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
searching molecular databases with blast2
Searching Molecular Databases with BLAST
  • Basic Local Alignment Search Tool
  • How BLAST works
  • Interpreting search results
  • The NCBI Web BLAST interface
  • Demonstration and exercises
why learn sequence database searching
Why learn sequence database searching?
  • What have I cloned ?
  • Is this really “my gene” ?
  • Has someone else already found it ?
  • What is this protein’s function ?
  • What is it related to ?
  • Can I get more sequence easily ?
search programs are sequence alignment programs
Search programs are sequence alignment programs
  • They try to find the best alignment between your probe sequence and every target sequence in the database
  • Finding optimal alignments is computationally a very resource intensive process
  • It is usually not necessary to find optimal alignments, particularly for large databases
  • Alignments are ranked and only top scores are reported
practical database search methods incorporate shortcuts
Practical database search methods incorporate shortcuts
  • The fastest sequence database searching programs use heuristic algorithms
  • The basic concept is to break the search and alignment process down into several steps
  • At each step, only a best scoring subset is retained for further analysis
what does heuristic mean
What does ‘HEURISTIC’ mean?
  • “a commonsense rule (or set of rules) intended to increase the probability of solving some problem”
  • Why consider every possible alignment once a reasonably good alignment is found?
heuristic programs find approximate alignments
Heuristic programs find approximate alignments
  • They are less sensitive than “dynamic programming” algorithms such as Smith-Waterman for detecting weak similarity
  • In practice, they run much faster and are usually adequate
  • The BLAST program developed by Stephen Altschul and coworkers at the NCBI is the most widely used heuristic program
blast is a collection of five programs for different combinations of query and database sequences
BLAST is a collection of five programs for different combinations of query and database sequences
why blast is great
Why BLAST is great
  • Very fast and can be used to search extremely large databases
  • Sufficiently sensitive and selective for most purposes
  • Robust - the default parameters can usually be used
blast scores are reported in two columns
BLAST scores are reported in two columns
  • Raw values based on the specific scoring matrix employed
  • As bits, which are matrix independent normalized values (bigger = better)
  • Significance is represented by E values (smaller = better)
the expect e threshold is used to control score reporting
The EXPECT (E) threshold is used to control score reporting
  • A match will only be reported if its E value falls below the threshold set
  • The default value for E is 10, which means that 10 matches with scores this high are expected to be found by chance
  • Lower EXPECT thresholds are more stringent, and report fewer matches
interpreting blast scores
Interpreting BLAST scores
  • Score interpretation is based on context
    • What is the question?
    • What else do you know about the sequences?
    • Scoring is highly dependent on probe length
  • Exact matches will usually have the highest scores (and lowest E values)
    • Short exact matches may score lower than longer partial matches
interpreting blast scores15
Interpreting BLAST scores
  • Short exact matches are expected to occur at random.
  • Partial matches over the entire length of a query are stronger evidence for homology than are short exact matches.
homology vs identity
Homology vs Identity
  • Homologous sequences are descended from a common ancestral sequence.
  • Homology is either true or false. It can never be partial! Saying two sequences are 45% homologous is a misuse of the term.
  • Sequence identity and similarity can be described as a percentage and are used as evidence of homology.
blast example
BLAST Example

Is this sequence known? What does it encode?

search strategy
Search Strategy
  • Choose the BLAST program:
    • nucleotide query vs. nucleotide db
    • megabalst: optimized to find identical sequences
    • blastn: will find identical and similar sequences
  • Choose the Database
    • nr (non-redundant) – everything
    • genome specific
blastn options
blastn Options

Paste Query

Sequence

HERE

Choose Database

HERE

Choose search program

HERE

slide21
Each line is a hit

in the database

sorted vertically

by E value

Colored rectangles along the X axis show where in the query sequence

a similarity in the database has been found. Color indicates degree of similarity

blastp databases31
blastp Databases
  • nr - All non-redundant GenBank CDS translations + PDB + SwissProt+PIR
  • swissprot - the last major release of the SWISS-PROT protein sequence database
  • pat - patented sequences
  • pdb - Sequences derived from the 3-dimensional structure Protein Data Bank
  • env_nr - Non-redundant environmental samples
blastp output
BLASTP Output

Conserved Domain Search

Conserved domains are shown

graphically. Link to explanation

of the domain.

protein scoring matrices
Protein Scoring Matrices

Blosom 62 is the default BLASTP scoring matrix

other blast programs psi blast
Other BLAST Programs:Psi-BLAST

4.6 PSI-BLAST is designed for more sensitive protein-protein similarity searches.

Position-Specific Iterated (PSI)-BLAST is the most sensitive BLAST program, making it useful for finding very distantly related proteins or new members of a protein family. Use PSI-BLAST when your standard protein-protein BLAST search either failed to find significant hits, or returned hits with descriptions such as "hypothetical protein" or "similar to...".

other blast programs phi blast
Other BLAST Programs:Phi-BLAST

4.7 PHI-BLAST can do a restricted protein pattern search.

Pattern-Hit Initiated (PHI)-BLAST is designed to search for proteins that contain a pattern specified by the user AND are similar to the query sequence in the vicinity of the pattern. This dual requirement is intended to reduce the number of database hits that contain the pattern, but are likely to have no true homology to the query.

sequence filters
Sequence filters
  • Since only a limited number of matches are reported, hits to simple repeats and other low complexity sequences can obscure other more biologically meaningful similarities
  • Filters are used to remove low complexity sequences from the probe
  • Low Complexity, human repeats (blastn)
blastn vs blastp
BLASTN vs BLASTP
  • Protein sequences have much higher information content than nucleotide sequence
  • To find evidence for sequence homology, use BLASTP and search protein sequences
  • Is my sequence already in the database?
  • To find identical sequences, search nucleotide databases
translated blast searches
Translated BLAST Searches
  • translations use all 6 frames
  • computationally intensive
  • tblastx searches can be very slow with some large databases
  • must specify genetic code
blast tutorial
BLAST Tutorial
  • BLAST tutorial on Biocomp Web page
  • Goal: demonstrate utility and difference between BLASTN and BLASTP searches
  • BLASTN: is my DNA sequence in the database?
  • BLASTP: are there related (homologs) proteins in the database?
ad