Introduction to Computation & Pairwise Alignment
|
|
- Cody Douglas
- 6 years ago
- Views:
Transcription
1 Introduction to Computation & Pairwise Alignment Eunok Paek
2 Algorithm what you already know about programming Pan-Fried Fish with Spicy Dipping Sauce This spicy fish dish is quick to prepare and cooks in about 8 minutes. Ingredients: ½ c mayonnaise ½ t salt ½ t cayenne pepper ¼ t ground black pepper 2 T lemon juice 2 eggs, beaten 4 white fish fillets (6 oz.) 1 c bread crumbs 3 T vegetable oil Directions: In a small bowl whisk together mayonnaise, cayenne pepper and lemon juice; set aside. Season fish fillets with salt and pepper to taste. Dip in beaten egg and coat evenly with bread crumbs. Heat a large, nonstick skillet over medium-high heat. Add oil and when hot, but not smoking, saute fish until golden brown and thoroughly cooked, about 4 minutes per side. Serve warm with reserved spicy dipping sauce.
3 Algorithm what you already know about programming Recipes have to be refined - A new recipe is rarely right on the first attempt. - Modifications are made as necessary. - Trying the recipe on the intended audience may yield further modifications. - The recipe can be adapted for new ingredients. Writing a program is a lot like writing a recipe.
4 Algorithm Definition An algorithm is a finite set of precise instructions for performing a computation or for solving a problem Example: find a maximum value in a finite sequence of integers 1. Set the temporary maximum equal to the first integer in the sequence. 2. Compare the next integer in the sequence to the temporary maximum equal to this integer. 3. Repeat the previous step if there are more integers in the sequence. 4. Stop when there are no integers left in the sequence. The temporary maximum at this point is the largest integer in the sequence.
5 Algorithm Pseudo code assignment variable = value repetition (iteration) function name arguments function max (a 1, a 2,, a n : integer) max = a 1 ; for i = 2 to n if max < a i then max = a i ; return max; data type function value
6 Algorithm Pseudo code function binary search (x: integer, a 1, a 2,, a n : increasing integers) i = 1; j = n; repetition while i < j begin m = (i + j / 2); if x > a m then i = m + 1 else j = m; end if x = a i then location = i else location = 0; return location;
7 Algorithm Pseudo code function n_choose_k (n, k: integers) calling another function return Factorial(n) / (Factorial(n k) * Factorial(k)); function Factorial (n: integer) temp = 1; for i = 2 to n temp = temp * i; return temp;
8 Algorithm Recursion function fibonacci (n: nonnegative integer) if n = 0 then return 0 else if n = 1 then return 1 else return fibonacci(n-1) + fibonacci(n-2); recursive call F 4 F 4 F 3 F 2 F 2 F 3 F 2 F 1 F 1 F 0 F 1 F 0 F 0 F 1
9 Algorithm Iteration & Memory function fibonacci (n: nonnegative integer) if n = 0 then return 0 else begin what if fn_2 = 0; n = 1? fn_1 = 1; for i = 1 to n-1; begin fn = fn_1 + fn_2; fn_2 = fn_1; fn_1 = fn; end end return fn_1;
10 Computation Running Time Two ways to measure relative efficiency of an algorithm Mathematical analysis Empirical analysis Mathematical analysis of the running time Running time is measured by the number of basic steps (e.g., the number of python statements) that the algorithm makes. Running time is described as a function of input size, Tn ( ) We are usually interested in the worst case running time or average case running time. n
11 Computation Big-Oh(O) Notation Example: T(n) = 13n n 2 + 2nlogn + 4n As n grows larger, n 3 is MUCH larger than n 2, nlogn, and n, so it dominates T(n) The constant factor 13 can be ignored since it is affected by the compiler used or machine speed, etc. The running time grows roughly on the order of n 3 Notationally, T(n)=O(n 3 )
12 Computation Complexity g(n) n n 3 /2 5n 2 100n n
13 Computation Complexity Function Approximate Values n nlogn n 3 1,000 1,000, n n x x n logn x x10 29 n! 3,628, x
14 Computation Complexity Function Size of Instance Solved in One Day Size of Instance Solved in a Computer 10 Times Faster n nlogn 0.948x x10 12 n x10 6 n x n n n n logn n! 14 15
15 Take Home Message There can be many ways to solve the same problem. Running time can often be estimated mathematically, using parameter of input size n. What matter is the order of growth in computational time.
16 Sequence Alignment A C C T G A G A G A C G T G G C A G mismatch 70% identical indel
17 Sequence Alignment Eye of the tiger * In 1994 Walter Gehring et alum (Un. Basel) turn the gene eyeless on in various places on Drosophila melanogaster * Result: on multiple places eyes are formed * eyeless is a master regulatory gene that controls +/ other genes * eyeless on induces formation of an eye
18 Sequence Alignment Eyeless Drosophila
19 Sequence Alignment
20 Sequence Alignment Homeoboxes & Master regulatory genes
21 Sequence Alignment Homeoboxes & Master regulatory genes HOMEO BOX A homeobox is a DNA sequence found within genes that are involved in the regulation of development (morphogenesis) of animals, fungi and plants.
22 Sequence Alignment Sequence alignment is important for: * prediction of function * database searching * gene finding * sequence divergence * sequence assembly 22
23 Growth of GenBank and WGS
24 Pairwise Alignment Dot matrix Dynamic programming Needleman-Wunsch optimal global alignment Smith-Waterman optimal local alignment
25 Pairwise Alignment Types of Sequence Alignment Dot matrix Number of sequences pairwise alignment: compare two sequences multiple alignment: compare multiple sequences Portion of sequences aligned global alignment: align sequences over their entire lengths local alignment: find the longest/best subsequence pairs that give maximum similarity Algorithmic approach optimal methods: Needleman-Wunsch, Smith-Waterman heuristic methods: FASTA, BLAST
26 Pairwise Alignment Dot Matrix Dot Matrix A visual depiction of relationship between 2 sequences Reveals insertion/deletion Finds direct or inverted repeats Steps create a 2D matrix one sequence along the top the other along the left side for each cell of the matrix, place a dot if the two corresponding residues match
27 Pairwise Alignment Dot Matrix Running Time of Dot Matrix Lengths of sequences: m, n O(mn)
28 Pairwise Alignment Dot Matrix DNA sequences protein sequences
29 Pairwise Alignment Dot Matrix Random Matches in Dot Matrix When comparing DNA sequences, random matches occur with probability 1/4 When comparing protein sequences, 1/20 Thus, for comparisons of protein coding DNA sequences, we should translate them to amino acids first
30 Pairwise Alignment Dot Matrix To Reduce Random Noise in Dot Matrix Specify a window size, w Take w residues from each of the two sequences Among the w pairs of residues, count how many pairs are matches Specify a stringency
31 Pairwise Alignment Dot Matrix Simple dot matrix, Window size 1 P V I L E P M M K V T I E M P P V 1 1 I 1 1 L 1 E 1 1 P I 1 1 M R V 1 1 E 1 1 V 1 1 T 1 T 1 P 1 1 1
32 Pairwise Alignment Dot Matrix Window size is 3 P V I L E P M M K V T I E M P P V I L E P I M R V E V T T P
33 Pairwise Alignment Dot Matrix Window size is 3; Stringency is 2 P V I L E P M M K V T I E M P P 3 V 3 I 3 L 3 E 2 P 2 I M 2 R V E 2 V 2 T T 2 2 P 3
34 Pairwise Alignment Dot Matrix DNA Sequences single residue identity 16 out of 23 identical
35 Pairwise Alignment Dot Matrix Protein Sequences single residue identity 6 out of 23 identical
36 Pairwise Alignment Dot Matrix Insertion/Deletion, Inversion
37 Pairwise Alignment Dot Matrix ABCDEFGEFGHIJKLMNO tandem duplication compared to no duplication tandem duplication compared to self
38 Pairwise Alignment Dot Matrix What Is This? 5 GGCGG 3 Palindrome (Intrastrand)
39 Pairwise Alignment Global Alignment Optimal Alignment Consider two sequences, both of length n If no gaps are allowed, there is only one alignment, which is optimal If n gaps are allowed, there are possible alignments How to find the optimal ones? n (2n)! ( n ) 2 ( n!) 2n 2 2 n
40 Pairwise Alignment Global Alignment First, Define Optimality Scoring scheme a scoring matrix and gap penalties Examples of scoring schemes amino acids: PAM250, or BLOSUM62; -13 for gap opening, -2 for gap extension nucleotides: the matrix to the right; -8 for gap opening, -6 for gap extension A C G T A C G T
41 Pairwise Alignment Global Alignment Intuition of Dynamic Programming If we already have the optimal solution to: XY AB then we know the next pair of characters will either be: XYZ or XY- or XYZ ABC ABC AB- (where - indicates a gap). So we can extend the match by determining which of these has the highest score.
42 Pairwise Alignment Global Alignment Recursive Definition of Dynamic Programming Notations: F(i,j): the accumulated score of aligning x 1, x 2,, x i to y 1,, y j s(x,y): the score of matching residue x to residue y, from the scoring matrix (k): the penalty for a gap of length k F ( i, j) max F F F ( i 1, j 1) ( k, j) ( i ( i, k ) ( j s( x k ), k ), i, y k k j ), 0,..., 0,..., i 1, j 1.
43 Pairwise Alignment Global Alignment Illustration of Dynamic Programming X Y Z U V W
44 Pairwise Alignment Global Alignment Dynamic Programming: Units of Operations Y 1 Y 2 Y 3 Y 4 Y n total X n X n+1 (n+4)(n-1)/2+1 = (n 2 +3n-4)/2+1 X n+2 (n+6)(n-1)/2+1 = (n 2 +5n-6)/2+1 X n+3 (n+8)(n-1)/2+1 = (n 2 +7n-8)/2+1 X n 1 n+1 n+2 n+3 2n-1 (n+2n)(n-1)/2+1 = (n 2 +(2n-1)n-2n)/2+1 [n 2 (n-1)+n(n+1)(n-1)-(n+2)(n-1)]/2+2n-1 = [2n 3-3n 2 -n+2]/2 +2n -1 O(n 3 ) units of operations
45 Pairwise Alignment Global Alignment The Needleman-Wunsch Algorithm The method described in the previous slides is the Needleman- Wunsch (1970) algorithm It computes the optimal global alignment between two sequences The optimality is defined in terms of a scoring scheme (a scoring matrix plus gap penalties) The running time is O(n 3 )
46 Pairwise Alignment Global Alignment Needleman-Wunsch Implementation Details F( i, F( i 1, j 1) s( xi, j) maxf( k, j) ( i k), F( i, k) ( j k), At each cell of the matrix, keep track of how the maximum is arrived at After the entire matrix is filled, do a traceback from the bottom right corner to the top left corner y k j ), 0,..., i 1, k 0,..., j 1. A B C ABCDEFG-HIJ I J
47 Pairwise Alignment Global Alignment Gap Penalties Above, the function of gap penalties can take any form Below, using a simple gap penalty (-d for each gap position), we can speed up the alignment algorithm 1. 0,..., ), ( ), ( 1, 0,..., ), ( ), ( ),, ( 1) 1, ( max ), ( j k k j k i F i k k i j k F y x s j i F j i F j i. 1), (, ) 1, ( ),, ( 1) 1, ( max ), ( d j i F d j i F y x s j i F j i F j i
48 Pairwise Alignment Global Alignment Illustration of Gotoh s Algorithm X Y Z 0 -d -2d -3d U -d V -2d W -3d
49 Pairwise Alignment Global Alignment Example: match 1, mismatch -1, gap -1 A C G T A G C T
50 Pairwise Alignment Global Alignment Gotoh s Algorithm: Units of Operations O(n 2 ) units of operations to fill the matrix O(n) units to trace back Y 1 Y 2 Y 3 Y 4 Y n total n+1 X n+1 X n+1 X n+1 X n+1 X n n+1 3n 2 +2n+1
51 Pairwise Alignment Global Alignment Affine Gap Penalties -d for gap opening -e for gap extension (k) = -d - e (k-1) Running time is still O(n 2 ) Described in Gotoh (1982) Optimal global alignment F( i, R( i, C( i, F( i 1, j) maxr( i 1, C( i 1, j 1) F( i 1, j) d, j) max R( i 1, j) e. F( i, j 1) d, j) max C( i, j 1) e. s( x i j 1) s( xi, y j 1) s( x, y i, y j j j ), ), ).
52 Pairwise Alignment Local Alignment Smith-Waterman Running time is O(n 2 ) Described in Smith and Waterman (1981) Optimal local alignment Traceback is different F( i, R( i, C( i, F( i 1, j 1) R( i 1, j 1) j) max C( i 1, j 1) 0. F( i 1, j) d, j) max R( i 1, j) e. F( i, j 1) d, j) max C( i, j 1) e. s( x s( x s( x i i i,,, y y y j j j ), ), ),
53 Pairwise Alignment Local Alignment Global versus Local Alignments LGPSSKQTGKGS-SRIWDN (Global) LN-ITKSAGKGAIMRLGDA TGKG (Local) AGKG
54 Pairwise Alignment Local Alignment Smith-Waterman Traceback H E A G A W G H E E P A W H E A E
55 Pairwise Alignment Significance of Alignment Probability of Random Alignments Suppose we have a tetrahedron-shaped die whose four faces are labeled with A, C, G, T. Throw the die twice, and record the labels facing down. Probability of getting an identical pair: ¼*¼. There are 4 possible identical pairs: 4*¼*¼ = ¼. 6 identical pairs = (1/4)^6 = 2.4E-4. Probability of getting a mismatch: 1 ¼ = ¾. 6 mismatched pairs is (3/4)^6 =
56 Pairwise Alignment Significance of Alignment If A, C, G, T are not of Equal Proportions Probability of drawing an identical pair is given by: p p 2 A p 2 C p 2 G p 2 T p x is proportion of nucleotide x Probability of drawing a mismatch is 1 - p
57 Pairwise Alignment Significance of Alignment Longest Run of Heads in Coin Toss HTTHHHTHHTHHHTTTHHHHHHHTTTHHT Probability of head is p. We are looking at a sequence of length n. At a random position, probability of seeing a run of 5 heads p 5 There are n 4such positions Frequency of observing such a run is p 5 (n 4). In general, p K (n (K 1)). (Erdos-Renyi law, 1970) For large n, K = log 1/p n. Expected length of the longest run of heads: If p=0.5, after 100 tosses, the longest run is log = 6.65
58 Pairwise Alignment Significance of Alignment M: Longest Run in Random Alignment Sequence lengths: m, n p: probability of match q: 1 p γ: Euler s number, E(M) log 1/p (mn) + log 1/p (q) + γlog(e) ½, for large m, n If a local alignment is longer than E(M), then it is significant How significant?
59 Pairwise Alignment Significance of Alignment Significance of Local Alignment In biological experiments, after a set of values of an entity is obtained, we usually calculate the mean and variance Assume data follows the normal distribution The mean and variance are of interest For example, is the mean not equal to zero at the significance level of 0.05? This is not what we want in local alignment We want the significance of the highest scores not the mean score
60 Pairwise Alignment Significance of Alignment Distribution of Scores The scores of a pair of sequences are compared to those of two random sequences of the same length and composition The distribution of random sequence scores follows the Gumbel extreme value distribution Similar to the normal distribution, with a positively skewed tail The score must be greater than expected from a normal distribution to achieve the same level of significance
61 Pairwise Alignment Significance of Alignment Normal Distribution versus Extreme Value Distribution 0.4 Normal Extreme Value Normal distribution: y = exp(-x 2 /2) / sqrt(2π) Extreme value distribution: y = exp(-x exp(-x)) x
62 Pairwise Alignment Substitution Matrices DNA PAM 1 Matrix PAM 1 corresponds to 1% mutations, 99% conservation. Assume 4 nucleotides are present at equal frequencies Assume all mutations from any nucleotide to any other are equally likely A C G T A C G T A uniform model M
63 Pairwise Alignment Substitution Matrices Transitions and Transversions Purines: A and G Pyrimidines: C and T Transitions: more often purine to purine pyrimidine to pyrimidine Transversions: less often from purine to pyrimidine from pyrimidine to purine
64 Pairwise Alignment Substitution Matrices Another DNA PAM 1 Matrix Assume 4 nucleotides are present at equal frequencies Assume transitions are 3 times more often than transversions A biased model A C G T A C G T
65 Pairwise Alignment Substitution Matrices The Meaning of the Score of an Alignment Assume ACGT is aligned to CCGT Given a model (matrix) M Want: odds ratio Pr(A C) Pr(C C) Pr(G G) Pr(T T) given the model (P A M AC ) (P C M CC )(P G M GG )(P T M TT ) Divided by Pr(A C) Pr(C C) Pr(G G) Pr(T T) happened by chance (P A P C ) (P C P C )(P G P G )(P T P T ) Compute: Let S XY = log 2 (P X M XY / P X P Y ) S = S AC + S CC + S GG + S TT, log odds ratio 2 S is what we want (odds ratio)
66 Pairwise Alignment Substitution Matrices From PAM 1 Mutation Probability Matrix to PAM1 Log Odds Ratio Matrix A C G T A C G T A A C C G G T T
67 Pairwise Alignment Substitution Matrices From Another PAM 1 Mutation Probability Matrix to PAM1 Log Odds Ratio Matrix A C G T A C G T A A C C G G T T
68 Pairwise Alignment Substitution Matrices From PAM 1 to PAM 2 PAM 2 = PAM 1 * PAM 1 = (PAM 1 ) 2 PAM 2 (A C): PAM 1 (A A)*PAM 1 (A C) + PAM 1 (A C)*PAM 1 (C C) + PAM 1 (A G)*PAM 1 (G C) + PAM 1 (A T)*PAM 1 (T C) Markov process: the probability of change from nucleotide A to nucleotide C is the same, regardless of previous changes at the site or the position of the site in the sequence
69 Pairwise Alignment Substitution Matrices Amino Acid PAM Matrices Percent Accepted Mutation Dayhoff (1978), 1572 changes in 71 families of proteins, at least 85% similar For each amino acid, count 20 numbers For example, how many F (phenylalanine) stay the same, how many change to the other 19 amino acids Normalize: divide each of these 20 numbers by (sum of 20 numbers) PAM1: 1% probability of change
70 Pairwise Alignment Substitution Matrices The Column/Row of F in PAM1 F to A: F to R: F to N: F to D: F to C: F to Q: F to E: F to G: F to H: F to I: F to L: F to K: F to M: F to F: F to P: F to S: F to T: F to W: F to Y: F to V:
71 Pairwise Alignment Substitution Matrices Compute PAM250 PAM 2 = PAM 1 * PAM 1 = (PAM 1 ) 2 PAM 250 = (PAM1) 250 Convert to log odds: PAM 250 (F Y) = 0.15 Divide by the frequency of F, /0.04 = 3.75 log 10 (3.75) = 0.57 Similarly for Y F: log 10 (0.2/0.03) = 0.83 So PAM250(F Y) = 10*( )/2
72 Pairwise Alignment Substitution Matrices BLOSUM BLOcks of amino acid SUbstitution Matrices Start with highly-conserved patterns (blocks) in a large set of closely related proteins Use the likelihood of substitutions found in those sequences to create a substitution probability matrix BLOSUM-n means that the sequences used were n% alike BLOSUM62 is standard
Algorithms in Bioinformatics FOUR Pairwise Sequence Alignment. Pairwise Sequence Alignment. Convention: DNA Sequences 5. Sequence Alignment
Algorithms in Bioinformatics FOUR Sami Khuri Department of Computer Science San José State University Pairwise Sequence Alignment Homology Similarity Global string alignment Local string alignment Dot
More information3. SEQUENCE ANALYSIS BIOINFORMATICS COURSE MTAT
3. SEQUENCE ANALYSIS BIOINFORMATICS COURSE MTAT.03.239 25.09.2012 SEQUENCE ANALYSIS IS IMPORTANT FOR... Prediction of function Gene finding the process of identifying the regions of genomic DNA that encode
More informationSara C. Madeira. Universidade da Beira Interior. (Thanks to Ana Teresa Freitas, IST for useful resources on this subject)
Bioinformática Sequence Alignment Pairwise Sequence Alignment Universidade da Beira Interior (Thanks to Ana Teresa Freitas, IST for useful resources on this subject) 1 16/3/29 & 23/3/29 27/4/29 Outline
More informationSequence analysis and Genomics
Sequence analysis and Genomics October 12 th November 23 rd 2 PM 5 PM Prof. Peter Stadler Dr. Katja Nowick Katja: group leader TFome and Transcriptome Evolution Bioinformatics group Paul-Flechsig-Institute
More informationCISC 889 Bioinformatics (Spring 2004) Sequence pairwise alignment (I)
CISC 889 Bioinformatics (Spring 2004) Sequence pairwise alignment (I) Contents Alignment algorithms Needleman-Wunsch (global alignment) Smith-Waterman (local alignment) Heuristic algorithms FASTA BLAST
More informationSequence Alignments. Dynamic programming approaches, scoring, and significance. Lucy Skrabanek ICB, WMC January 31, 2013
Sequence Alignments Dynamic programming approaches, scoring, and significance Lucy Skrabanek ICB, WMC January 31, 213 Sequence alignment Compare two (or more) sequences to: Find regions of conservation
More informationAlgorithms in Bioinformatics
Algorithms in Bioinformatics Sami Khuri Department of omputer Science San José State University San José, alifornia, USA khuri@cs.sjsu.edu www.cs.sjsu.edu/faculty/khuri Pairwise Sequence Alignment Homology
More informationComputational Biology
Computational Biology Lecture 6 31 October 2004 1 Overview Scoring matrices (Thanks to Shannon McWeeney) BLAST algorithm Start sequence alignment 2 1 What is a homologous sequence? A homologous sequence,
More informationCONCEPT OF SEQUENCE COMPARISON. Natapol Pornputtapong 18 January 2018
CONCEPT OF SEQUENCE COMPARISON Natapol Pornputtapong 18 January 2018 SEQUENCE ANALYSIS - A ROSETTA STONE OF LIFE Sequence analysis is the process of subjecting a DNA, RNA or peptide sequence to any of
More informationBioinformatics (GLOBEX, Summer 2015) Pairwise sequence alignment
Bioinformatics (GLOBEX, Summer 2015) Pairwise sequence alignment Substitution score matrices, PAM, BLOSUM Needleman-Wunsch algorithm (Global) Smith-Waterman algorithm (Local) BLAST (local, heuristic) E-value
More informationPractical considerations of working with sequencing data
Practical considerations of working with sequencing data File Types Fastq ->aligner -> reference(genome) coordinates Coordinate files SAM/BAM most complete, contains all of the info in fastq and more!
More information5. MULTIPLE SEQUENCE ALIGNMENT BIOINFORMATICS COURSE MTAT
5. MULTIPLE SEQUENCE ALIGNMENT BIOINFORMATICS COURSE MTAT.03.239 03.10.2012 ALIGNMENT Alignment is the task of locating equivalent regions of two or more sequences to maximize their similarity. Homology:
More informationSequence Alignment: A General Overview. COMP Fall 2010 Luay Nakhleh, Rice University
Sequence Alignment: A General Overview COMP 571 - Fall 2010 Luay Nakhleh, Rice University Life through Evolution All living organisms are related to each other through evolution This means: any pair of
More informationTools and Algorithms in Bioinformatics
Tools and Algorithms in Bioinformatics GCBA815, Fall 2013 Week3: Blast Algorithm, theory and practice Babu Guda, Ph.D. Department of Genetics, Cell Biology & Anatomy Bioinformatics and Systems Biology
More informationBioinformatics and BLAST
Bioinformatics and BLAST Overview Recap of last time Similarity discussion Algorithms: Needleman-Wunsch Smith-Waterman BLAST Implementation issues and current research Recap from Last Time Genome consists
More informationPairwise Alignment. Guan-Shieng Huang. Dept. of CSIE, NCNU. Pairwise Alignment p.1/55
Pairwise Alignment Guan-Shieng Huang shieng@ncnu.edu.tw Dept. of CSIE, NCNU Pairwise Alignment p.1/55 Approach 1. Problem definition 2. Computational method (algorithms) 3. Complexity and performance Pairwise
More informationPairwise sequence alignments
Pairwise sequence alignments Volker Flegel VI, October 2003 Page 1 Outline Introduction Definitions Biological context of pairwise alignments Computing of pairwise alignments Some programs VI, October
More informationBiochemistry 324 Bioinformatics. Pairwise sequence alignment
Biochemistry 324 Bioinformatics Pairwise sequence alignment How do we compare genes/proteins? When we have sequenced a genome, we try and identify the function of unknown genes by finding a similar gene
More informationLecture 4: Evolutionary Models and Substitution Matrices (PAM and BLOSUM)
Bioinformatics II Probability and Statistics Universität Zürich and ETH Zürich Spring Semester 2009 Lecture 4: Evolutionary Models and Substitution Matrices (PAM and BLOSUM) Dr Fraser Daly adapted from
More informationSingle alignment: Substitution Matrix. 16 march 2017
Single alignment: Substitution Matrix 16 march 2017 BLOSUM Matrix BLOSUM Matrix [2] (Blocks Amino Acid Substitution Matrices ) It is based on the amino acids substitutions observed in ~2000 conserved block
More information20 Grundlagen der Bioinformatik, SS 08, D. Huson, May 27, Global and local alignment of two sequences using dynamic programming
20 Grundlagen der Bioinformatik, SS 08, D. Huson, May 27, 2008 4 Pairwise alignment We will discuss: 1. Strings 2. Dot matrix method for comparing sequences 3. Edit distance 4. Global and local alignment
More informationSequence Analysis 17: lecture 5. Substitution matrices Multiple sequence alignment
Sequence Analysis 17: lecture 5 Substitution matrices Multiple sequence alignment Substitution matrices Used to score aligned positions, usually of amino acids. Expressed as the log-likelihood ratio of
More informationPairwise sequence alignment
Department of Evolutionary Biology Example Alignment between very similar human alpha- and beta globins: GSAQVKGHGKKVADALTNAVAHVDDMPNALSALSDLHAHKL G+ +VK+HGKKV A+++++AH+D++ +++++LS+LH KL GNPKVKAHGKKVLGAFSDGLAHLDNLKGTFATLSELHCDKL
More informationTHEORY. Based on sequence Length According to the length of sequence being compared it is of following two types
Exp 11- THEORY Sequence Alignment is a process of aligning two sequences to achieve maximum levels of identity between them. This help to derive functional, structural and evolutionary relationships between
More informationMotivating the need for optimal sequence alignments...
1 Motivating the need for optimal sequence alignments... 2 3 Note that this actually combines two objectives of optimal sequence alignments: (i) use the score of the alignment o infer homology; (ii) use
More informationPairwise sequence alignments. Vassilios Ioannidis (From Volker Flegel )
Pairwise sequence alignments Vassilios Ioannidis (From Volker Flegel ) Outline Introduction Definitions Biological context of pairwise alignments Computing of pairwise alignments Some programs Importance
More informationModule: Sequence Alignment Theory and Applications Session: Introduction to Searching and Sequence Alignment
Module: Sequence Alignment Theory and Applications Session: Introduction to Searching and Sequence Alignment Introduction to Bioinformatics online course : IBT Jonathan Kayondo Learning Objectives Understand
More informationCollected Works of Charles Dickens
Collected Works of Charles Dickens A Random Dickens Quote If there were no bad people, there would be no good lawyers. Original Sentence It was a dark and stormy night; the night was dark except at sunny
More informationSequence Comparison. mouse human
Sequence Comparison Sequence Comparison mouse human Why Compare Sequences? The first fact of biological sequence analysis In biomolecular sequences (DNA, RNA, or amino acid sequences), high sequence similarity
More informationIntroduction to sequence alignment. Local alignment the Smith-Waterman algorithm
Lecture 2, 12/3/2003: Introduction to sequence alignment The Needleman-Wunsch algorithm for global sequence alignment: description and properties Local alignment the Smith-Waterman algorithm 1 Computational
More informationLocal Alignment: Smith-Waterman algorithm
Local Alignment: Smith-Waterman algorithm Example: a shared common domain of two protein sequences; extended sections of genomic DNA sequence. Sensitive to detect similarity in highly diverged sequences.
More informationbioinformatics 1 -- lecture 7
bioinformatics 1 -- lecture 7 Probability and conditional probability Random sequences and significance (real sequences are not random) Erdos & Renyi: theoretical basis for the significance of an alignment
More informationPairwise & Multiple sequence alignments
Pairwise & Multiple sequence alignments Urmila Kulkarni-Kale Bioinformatics Centre 411 007 urmila@bioinfo.ernet.in Basis for Sequence comparison Theory of evolution: gene sequences have evolved/derived
More informationLecture 1, 31/10/2001: Introduction to sequence alignment. The Needleman-Wunsch algorithm for global sequence alignment: description and properties
Lecture 1, 31/10/2001: Introduction to sequence alignment The Needleman-Wunsch algorithm for global sequence alignment: description and properties 1 Computational sequence-analysis The major goal of computational
More informationBLAST: Target frequencies and information content Dannie Durand
Computational Genomics and Molecular Biology, Fall 2016 1 BLAST: Target frequencies and information content Dannie Durand BLAST has two components: a fast heuristic for searching for similar sequences
More information1.5 Sequence alignment
1.5 Sequence alignment The dramatic increase in the number of sequenced genomes and proteomes has lead to development of various bioinformatic methods and algorithms for extracting information (data mining)
More informationLarge-Scale Genomic Surveys
Bioinformatics Subtopics Fold Recognition Secondary Structure Prediction Docking & Drug Design Protein Geometry Protein Flexibility Homology Modeling Sequence Alignment Structure Classification Gene Prediction
More informationSequence Analysis '17 -- lecture 7
Sequence Analysis '17 -- lecture 7 Significance E-values How significant is that? Please give me a number for......how likely the data would not have been the result of chance,......as opposed to......a
More informationQuantifying sequence similarity
Quantifying sequence similarity Bas E. Dutilh Systems Biology: Bioinformatic Data Analysis Utrecht University, February 16 th 2016 After this lecture, you can define homology, similarity, and identity
More informationPairwise alignment, Gunnar Klau, November 9, 2005, 16:
Pairwise alignment, Gunnar Klau, November 9, 2005, 16:36 2012 2.1 Growth rates For biological sequence analysis, we prefer algorithms that have time and space requirements that are linear in the length
More informationStudy and Implementation of Various Techniques Involved in DNA and Protein Sequence Analysis
Study and Implementation of Various Techniques Involved in DNA and Protein Sequence Analysis Kumud Joseph Kujur, Sumit Pal Singh, O.P. Vyas, Ruchir Bhatia, Varun Singh* Indian Institute of Information
More informationBioinformatics for Biologists
Bioinformatics for Biologists Sequence Analysis: Part I. Pairwise alignment and database searching Fran Lewitter, Ph.D. Head, Biocomputing Whitehead Institute Bioinformatics Definitions The use of computational
More informationPairwise Sequence Alignment
Introduction to Bioinformatics Pairwise Sequence Alignment Prof. Dr. Nizamettin AYDIN naydin@yildiz.edu.tr Outline Introduction to sequence alignment pair wise sequence alignment The Dot Matrix Scoring
More informationLocal Alignment Statistics
Local Alignment Statistics Stephen Altschul National Center for Biotechnology Information National Library of Medicine National Institutes of Health Bethesda, MD Central Issues in Biological Sequence Comparison
More informationTools and Algorithms in Bioinformatics
Tools and Algorithms in Bioinformatics GCBA815, Fall 2015 Week-4 BLAST Algorithm Continued Multiple Sequence Alignment Babu Guda, Ph.D. Department of Genetics, Cell Biology & Anatomy Bioinformatics and
More informationSequence comparison: Score matrices
Sequence comparison: Score matrices http://facultywashingtonedu/jht/gs559_2013/ Genome 559: Introduction to Statistical and omputational Genomics Prof James H Thomas FYI - informal inductive proof of best
More informationIn-Depth Assessment of Local Sequence Alignment
2012 International Conference on Environment Science and Engieering IPCBEE vol.3 2(2012) (2012)IACSIT Press, Singapoore In-Depth Assessment of Local Sequence Alignment Atoosa Ghahremani and Mahmood A.
More informationHidden Markov Models
Hidden Markov Models Outline 1. CG-Islands 2. The Fair Bet Casino 3. Hidden Markov Model 4. Decoding Algorithm 5. Forward-Backward Algorithm 6. Profile HMMs 7. HMM Parameter Estimation 8. Viterbi Training
More informationAn Introduction to Bioinformatics Algorithms Hidden Markov Models
Hidden Markov Models Outline 1. CG-Islands 2. The Fair Bet Casino 3. Hidden Markov Model 4. Decoding Algorithm 5. Forward-Backward Algorithm 6. Profile HMMs 7. HMM Parameter Estimation 8. Viterbi Training
More informationBioinformatics for Computer Scientists (Part 2 Sequence Alignment) Sepp Hochreiter
Bioinformatics for Computer Scientists (Part 2 Sequence Alignment) Institute of Bioinformatics Johannes Kepler University, Linz, Austria Sequence Alignment 2. Sequence Alignment Sequence Alignment 2.1
More informationSequence comparison: Score matrices. Genome 559: Introduction to Statistical and Computational Genomics Prof. James H. Thomas
Sequence comparison: Score matrices Genome 559: Introduction to Statistical and omputational Genomics Prof James H Thomas FYI - informal inductive proof of best alignment path onsider the last step in
More informationFirst generation sequencing and pairwise alignment (High-tech, not high throughput) Analysis of Biological Sequences
First generation sequencing and pairwise alignment (High-tech, not high throughput) Analysis of Biological Sequences 140.638 where do sequences come from? DNA is not hard to extract (getting DNA from a
More informationLecture 4: Evolutionary models and substitution matrices (PAM and BLOSUM).
1 Bioinformatics: In-depth PROBABILITY & STATISTICS Spring Semester 2011 University of Zürich and ETH Zürich Lecture 4: Evolutionary models and substitution matrices (PAM and BLOSUM). Dr. Stefanie Muff
More informationSequence analysis and comparison
The aim with sequence identification: Sequence analysis and comparison Marjolein Thunnissen Lund September 2012 Is there any known protein sequence that is homologous to mine? Are there any other species
More informationAn Introduction to Sequence Similarity ( Homology ) Searching
An Introduction to Sequence Similarity ( Homology ) Searching Gary D. Stormo 1 UNIT 3.1 1 Washington University, School of Medicine, St. Louis, Missouri ABSTRACT Homologous sequences usually have the same,
More informationAlignment & BLAST. By: Hadi Mozafari KUMS
Alignment & BLAST By: Hadi Mozafari KUMS SIMILARITY - ALIGNMENT Comparison of primary DNA or protein sequences to other primary or secondary sequences Expecting that the function of the similar sequence
More informationIntroduction to Bioinformatics
Introduction to Bioinformatics Jianlin Cheng, PhD Department of Computer Science Informatics Institute 2011 Topics Introduction Biological Sequence Alignment and Database Search Analysis of gene expression
More informationAlignment principles and homology searching using (PSI-)BLAST. Jaap Heringa Centre for Integrative Bioinformatics VU (IBIVU)
Alignment principles and homology searching using (PSI-)BLAST Jaap Heringa Centre for Integrative Bioinformatics VU (IBIVU) http://ibivu.cs.vu.nl Bioinformatics Nothing in Biology makes sense except in
More informationBioinformatics. Scoring Matrices. David Gilbert Bioinformatics Research Centre
Bioinformatics Scoring Matrices David Gilbert Bioinformatics Research Centre www.brc.dcs.gla.ac.uk Department of Computing Science, University of Glasgow Learning Objectives To explain the requirement
More informationSimilarity or Identity? When are molecules similar?
Similarity or Identity? When are molecules similar? Mapping Identity A -> A T -> T G -> G C -> C or Leu -> Leu Pro -> Pro Arg -> Arg Phe -> Phe etc If we map similarity using identity, how similar are
More informationHeuristic Alignment and Searching
3/28/2012 Types of alignments Global Alignment Each letter of each sequence is aligned to a letter or a gap (e.g., Needleman-Wunsch). Local Alignment An optimal pair of subsequences is taken from the two
More informationSequence comparison: Score matrices. Genome 559: Introduction to Statistical and Computational Genomics Prof. James H. Thomas
Sequence comparison: Score matrices Genome 559: Introduction to Statistical and omputational Genomics Prof James H Thomas Informal inductive proof of best alignment path onsider the last step in the best
More informationLecture 5,6 Local sequence alignment
Lecture 5,6 Local sequence alignment Chapter 6 in Jones and Pevzner Fall 2018 September 4,6, 2018 Evolution as a tool for biological insight Nothing in biology makes sense except in the light of evolution
More informationScoring Matrices. Shifra Ben-Dor Irit Orr
Scoring Matrices Shifra Ben-Dor Irit Orr Scoring matrices Sequence alignment and database searching programs compare sequences to each other as a series of characters. All algorithms (programs) for comparison
More informationProbability and random variables
Probability and random variables Events A simple event is the outcome of an experiment. For example, the experiment of tossing a coin twice has four possible outcomes: HH, HT, TH, TT. A compound event
More information8 Grundlagen der Bioinformatik, SoSe 11, D. Huson, April 18, 2011
8 Grundlagen der Bioinformatik, SoSe 11, D. Huson, April 18, 2011 2 Pairwise alignment We will discuss: 1. Strings 2. Dot matrix method for comparing sequences 3. Edit distance and alignment 4. The number
More informationMoreover, the circular logic
Moreover, the circular logic How do we know what is the right distance without a good alignment? And how do we construct a good alignment without knowing what substitutions were made previously? ATGCGT--GCAAGT
More informationLecture 2, 5/12/2001: Local alignment the Smith-Waterman algorithm. Alignment scoring schemes and theory: substitution matrices and gap models
Lecture 2, 5/12/2001: Local alignment the Smith-Waterman algorithm Alignment scoring schemes and theory: substitution matrices and gap models 1 Local sequence alignments Local sequence alignments are necessary
More informationBio nformatics. Lecture 3. Saad Mneimneh
Bio nformatics Lecture 3 Sequencing As before, DNA is cut into small ( 0.4KB) fragments and a clone library is formed. Biological experiments allow to read a certain number of these short fragments per
More informationC E N T R. Introduction to bioinformatics 2007 E B I O I N F O R M A T I C S V U F O R I N T. Lecture 5 G R A T I V. Pair-wise Sequence Alignment
C E N T R E F O R I N T E G R A T I V E B I O I N F O R M A T I C S V U Introduction to bioinformatics 2007 Lecture 5 Pair-wise Sequence Alignment Bioinformatics Nothing in Biology makes sense except in
More informationSequence Alignment Techniques and Their Uses
Sequence Alignment Techniques and Their Uses Sarah Fiorentino Since rapid sequencing technology and whole genomes sequencing, the amount of sequence information has grown exponentially. With all of this
More information8 Grundlagen der Bioinformatik, SS 09, D. Huson, April 28, 2009
8 Grundlagen der Bioinformatik, SS 09, D. Huson, April 28, 2009 2 Pairwise alignment We will discuss: 1. Strings 2. Dot matrix method for comparing sequences 3. Edit distance and alignment 4. The number
More informationCS 310 Advanced Data Structures and Algorithms
CS 310 Advanced Data Structures and Algorithms Runtime Analysis May 31, 2017 Tong Wang UMass Boston CS 310 May 31, 2017 1 / 37 Topics Weiss chapter 5 What is algorithm analysis Big O, big, big notations
More informationChapter 5. Proteomics and the analysis of protein sequence Ⅱ
Proteomics Chapter 5. Proteomics and the analysis of protein sequence Ⅱ 1 Pairwise similarity searching (1) Figure 5.5: manual alignment One of the amino acids in the top sequence has no equivalent and
More informationSequence Alignment: Scoring Schemes. COMP 571 Luay Nakhleh, Rice University
Sequence Alignment: Scoring Schemes COMP 571 Luay Nakhleh, Rice University Scoring Schemes Recall that an alignment score is aimed at providing a scale to measure the degree of similarity (or difference)
More informationSequence Alignment (chapter 6)
Sequence lignment (chapter 6) he biological problem lobal alignment Local alignment Multiple alignment Introduction to bioinformatics, utumn 6 Background: comparative genomics Basic question in biology:
More informationBIO 285/CSCI 285/MATH 285 Bioinformatics Programming Lecture 8 Pairwise Sequence Alignment 2 And Python Function Instructor: Lei Qian Fisk University
BIO 285/CSCI 285/MATH 285 Bioinformatics Programming Lecture 8 Pairwise Sequence Alignment 2 And Python Function Instructor: Lei Qian Fisk University Measures of Sequence Similarity Alignment with dot
More informationLecture 2: Pairwise Alignment. CG Ron Shamir
Lecture 2: Pairwise Alignment 1 Main source 2 Why compare sequences? Human hexosaminidase A vs Mouse hexosaminidase A 3 www.mathworks.com/.../jan04/bio_genome.html Sequence Alignment עימוד רצפים The problem:
More informationHidden Markov Models
Hidden Markov Models Outline CG-islands The Fair Bet Casino Hidden Markov Model Decoding Algorithm Forward-Backward Algorithm Profile HMMs HMM Parameter Estimation Viterbi training Baum-Welch algorithm
More informationGlobal alignments - review
Global alignments - review Take two sequences: X[j] and Y[j] M[i-1, j-1] ± 1 M[i, j] = max M[i, j-1] 2 M[i-1, j] 2 The best alignment for X[1 i] and Y[1 j] is called M[i, j] X[j] Initiation: M[,]= pply
More informationInDel 3-5. InDel 8-9. InDel 3-5. InDel 8-9. InDel InDel 8-9
Lecture 5 Alignment I. Introduction. For sequence data, the process of generating an alignment establishes positional homologies; that is, alignment provides the identification of homologous phylogenetic
More informationAdministration. ndrew Torda April /04/2008 [ 1 ]
ndrew Torda April 2008 Administration 22/04/2008 [ 1 ] Sprache? zu verhandeln (Englisch, Hochdeutsch, Bayerisch) Selection of topics Proteins / DNA / RNA Two halves to course week 1-7 Prof Torda (larger
More informationHidden Markov Models
Hidden Markov Models Slides revised and adapted to Bioinformática 55 Engª Biomédica/IST 2005 Ana Teresa Freitas CG-Islands Given 4 nucleotides: probability of occurrence is ~ 1/4. Thus, probability of
More informationAlgorithms in Bioinformatics: A Practical Introduction. Sequence Similarity
Algorithms in Bioinformatics: A Practical Introduction Sequence Similarity Earliest Researches in Sequence Comparison Doolittle et al. (Science, July 1983) searched for platelet-derived growth factor (PDGF)
More informationBINF 730. DNA Sequence Alignment Why?
BINF 730 Lecture 2 Seuence Alignment DNA Seuence Alignment Why? Recognition sites might be common restriction enzyme start seuence stop seuence other regulatory seuences Homology evolutionary common progenitor
More informationSubstitution matrices
Introduction to Bioinformatics Substitution matrices Jacques van Helden Jacques.van-Helden@univ-amu.fr Université d Aix-Marseille, France Lab. Technological Advances for Genomics and Clinics (TAGC, INSERM
More informationPairwise alignment. 2.1 Introduction GSAQVKGHGKKVADALTNAVAHVDDMPNALSALSD----LHAHKL
2 Pairwise alignment 2.1 Introduction The most basic sequence analysis task is to ask if two sequences are related. This is usually done by first aligning the sequences (or parts of them) and then deciding
More informationData structures Exercise 1 solution. Question 1. Let s start by writing all the functions in big O notation:
Data structures Exercise 1 solution Question 1 Let s start by writing all the functions in big O notation: f 1 (n) = 2017 = O(1), f 2 (n) = 2 log 2 n = O(n 2 ), f 3 (n) = 2 n = O(2 n ), f 4 (n) = 1 = O
More informationCopyright 2000 N. AYDIN. All rights reserved. 1
Introduction to Bioinformatics Prof. Dr. Nizamettin AYDIN naydin@yildiz.edu.tr Multiple Sequence Alignment Outline Multiple sequence alignment introduction to msa methods of msa progressive global alignment
More informationLecture 4: September 19
CSCI1810: Computational Molecular Biology Fall 2017 Lecture 4: September 19 Lecturer: Sorin Istrail Scribe: Cyrus Cousins Note: LaTeX template courtesy of UC Berkeley EECS dept. Disclaimer: These notes
More informationEECS730: Introduction to Bioinformatics
EECS730: Introduction to Bioinformatics Lecture 05: Index-based alignment algorithms Slides adapted from Dr. Shaojie Zhang (University of Central Florida) Real applications of alignment Database search
More informationMarkov Chains and Hidden Markov Models. = stochastic, generative models
Markov Chains and Hidden Markov Models = stochastic, generative models (Drawing heavily from Durbin et al., Biological Sequence Analysis) BCH339N Systems Biology / Bioinformatics Spring 2016 Edward Marcotte,
More informationBasic Local Alignment Search Tool
Basic Local Alignment Search Tool Alignments used to uncover homologies between sequences combined with phylogenetic studies o can determine orthologous and paralogous relationships Local Alignment uses
More informationPair Hidden Markov Models
Pair Hidden Markov Models Scribe: Rishi Bedi Lecturer: Serafim Batzoglou January 29, 2015 1 Recap of HMMs alphabet: Σ = {b 1,...b M } set of states: Q = {1,..., K} transition probabilities: A = [a ij ]
More informationString Matching Problem
String Matching Problem Pattern P Text T Set of Locations L 9/2/23 CAP/CGS 5991: Lecture 2 Computer Science Fundamentals Specify an input-output description of the problem. Design a conceptual algorithm
More informationSequence Alignment. Johannes Starlinger
Sequence Alignment Johannes Starlinger his Lecture Approximate String Matching Edit distance and alignment Computing global alignments Local alignment Johannes Starlinger: Bioinformatics, Summer Semester
More informationLecture 5: September Time Complexity Analysis of Local Alignment
CSCI1810: Computational Molecular Biology Fall 2017 Lecture 5: September 21 Lecturer: Sorin Istrail Scribe: Cyrus Cousins Note: LaTeX template courtesy of UC Berkeley EECS dept. Disclaimer: These notes
More informationLecture 2. Fundamentals of the Analysis of Algorithm Efficiency
Lecture 2 Fundamentals of the Analysis of Algorithm Efficiency 1 Lecture Contents 1. Analysis Framework 2. Asymptotic Notations and Basic Efficiency Classes 3. Mathematical Analysis of Nonrecursive Algorithms
More informationCSE : Computational Issues in Molecular Biology. Lecture 6. Spring 2004
CSE 397-497: Computational Issues in Molecular Biology Lecture 6 Spring 2004-1 - Topics for today Based on premise that algorithms we've studied are too slow: Faster method for global comparison when sequences
More informationWhole Genome Alignments and Synteny Maps
Whole Genome Alignments and Synteny Maps IINTRODUCTION It was not until closely related organism genomes have been sequenced that people start to think about aligning genomes and chromosomes instead of
More information