sequence alignment
DESCRIPTION
Sequence Alignment. Sequence s. Much of bioinformatics involves sequences DNA sequences RNA sequences Protein sequences We can think of these sequences as strings of letters DNA & RNA: alphabet of 4 letters Protein: alphabet of 20 letters. Glycine (G, GLY) Alanine (A, ALA) - PowerPoint PPT PresentationTRANSCRIPT
![Page 1: Sequence Alignment](https://reader030.vdocument.in/reader030/viewer/2022020417/56814d7d550346895dbada01/html5/thumbnails/1.jpg)
.
Sequence Alignment
![Page 2: Sequence Alignment](https://reader030.vdocument.in/reader030/viewer/2022020417/56814d7d550346895dbada01/html5/thumbnails/2.jpg)
Sequences
Much of bioinformatics involves sequences DNA sequences RNA sequences Protein sequences
We can think of these sequences as strings of letters DNA & RNA: alphabet of 4 letters Protein: alphabet of 20 letters
![Page 3: Sequence Alignment](https://reader030.vdocument.in/reader030/viewer/2022020417/56814d7d550346895dbada01/html5/thumbnails/3.jpg)
20 Amino Acids
Glycine (G, GLY) Alanine (A, ALA) Valine (V, VAL) Leucine (L, LEU) Isoleucine (I, ILE) Phenylalanine (F, PHE) Proline (P, PRO) Serine (S, SER) Threonine (T, THR) Cysteine (C, CYS) Methionine (M, MET) Tryptophan (W, TRP) Tyrosine (T, TYR) Asparagine (N, ASN) Glutamine (Q, GLN) Aspartic acid (D, ASP) Glutamic Acid (E, GLU) Lysine (K, LYS) Arginine (R, ARG) Histidine (H, HIS) START: AUG STOP: UAA, UAG, UGA
![Page 4: Sequence Alignment](https://reader030.vdocument.in/reader030/viewer/2022020417/56814d7d550346895dbada01/html5/thumbnails/4.jpg)
Sequence Comparison
Finding similarity between sequences is important for many biological questions
For example: Find genes/proteins with common origin
Allows to predict function & structure Locate common subsequences in genes/proteins
Identify common “motifs” Locate sequences that might overlap
Help in sequence assembly
![Page 5: Sequence Alignment](https://reader030.vdocument.in/reader030/viewer/2022020417/56814d7d550346895dbada01/html5/thumbnails/5.jpg)
Sequence Alignment
Input: two sequences over the same alphabet
Output: an alignment of the two sequences
Example: GCGCATGGATTGAGCGA TGCGCCATTGATGACCA
A possible alignment:
-GCGC-ATGGATTGAGCGA
TGCGCCATTGAT-GACC-A
![Page 6: Sequence Alignment](https://reader030.vdocument.in/reader030/viewer/2022020417/56814d7d550346895dbada01/html5/thumbnails/6.jpg)
Alignments
-GCGC-ATGGATTGAGCGA
TGCGCCATTGAT-GACC-A
Three elements: Perfect matches Mismatches Insertions & deletions (indel)
![Page 7: Sequence Alignment](https://reader030.vdocument.in/reader030/viewer/2022020417/56814d7d550346895dbada01/html5/thumbnails/7.jpg)
Choosing Alignments
There are many possible alignments
For example, compare:
-GCGC-ATGGATTGAGCGA
TGCGCCATTGAT-GACC-A
to
------GCGCATGGATTGAGCGA
TGCGCC----ATTGATGACCA--
Which one is better?
![Page 8: Sequence Alignment](https://reader030.vdocument.in/reader030/viewer/2022020417/56814d7d550346895dbada01/html5/thumbnails/8.jpg)
Scoring Alignments
Rough intuition: Similar sequences evolved from a common
ancestor Evolution changed the sequences from this
ancestral sequence by mutations: Replacements: one letter replaced by another Deletion: deletion of a letter Insertion: insertion of a letter
Scoring of sequence similarity should examine how many operations took place
![Page 9: Sequence Alignment](https://reader030.vdocument.in/reader030/viewer/2022020417/56814d7d550346895dbada01/html5/thumbnails/9.jpg)
Simple Scoring Rule
Score each position independently: Match: +1 Mismatch : -1 Indel -2
Score of an alignment is sum of positional scores
![Page 10: Sequence Alignment](https://reader030.vdocument.in/reader030/viewer/2022020417/56814d7d550346895dbada01/html5/thumbnails/10.jpg)
Example
Example:
-GCGC-ATGGATTGAGCGA
TGCGCCATTGAT-GACC-A
Score: (+1x13) + (-1x2) + (-2x4) = 3
------GCGCATGGATTGAGCGA
TGCGCC----ATTGATGACCA--
Score: (+1x5) + (-1x6) + (-2x11) = -23
![Page 11: Sequence Alignment](https://reader030.vdocument.in/reader030/viewer/2022020417/56814d7d550346895dbada01/html5/thumbnails/11.jpg)
More General Scores
The choice of +1,-1, and -2 scores was quite arbitrary
Depending on the context, some changes are more plausible than others Exchange of an amino-acid by one with similar
properties (size, charge, etc.)
vs. Exchange of an amino-acid by one with opposite
properties
![Page 12: Sequence Alignment](https://reader030.vdocument.in/reader030/viewer/2022020417/56814d7d550346895dbada01/html5/thumbnails/12.jpg)
For proteins
![Page 13: Sequence Alignment](https://reader030.vdocument.in/reader030/viewer/2022020417/56814d7d550346895dbada01/html5/thumbnails/13.jpg)
Additive Scoring Rules
We define a scoring function by specifying a function (x,y) is the score of replacing x by y (x,-) is the score of deleting x (-,x) is the score of inserting x
The score of an alignment is the sum of position scores
![Page 14: Sequence Alignment](https://reader030.vdocument.in/reader030/viewer/2022020417/56814d7d550346895dbada01/html5/thumbnails/14.jpg)
Edit Distance
The edit distance between two sequences is the “cost” of the “cheapest” set of edit operations needed to transform one sequence into the other
Computing edit distance between two sequences almost equivalent to finding the alignment that minimizes the distance
nment)score(aligmax),d( & of alignment 21 ss21 ss
![Page 15: Sequence Alignment](https://reader030.vdocument.in/reader030/viewer/2022020417/56814d7d550346895dbada01/html5/thumbnails/15.jpg)
Computing Edit Distance
How can we compute the edit distance?? If |s| = n and |t| = m, there are more than
alignments
The additive form of the score allows to perform dynamic programming to compute edit distance efficiently
m
nm
![Page 16: Sequence Alignment](https://reader030.vdocument.in/reader030/viewer/2022020417/56814d7d550346895dbada01/html5/thumbnails/16.jpg)
Recursive Argument
Define the notation:
Using the recursive argument, we get the following recurrence for V:
])[,(],[
)],[(],[
])[],[(],[
max],[
1jtj1iV
1is1jiV
1jt1isjiV
1j1iV
])..[],..[(],[ j1ti1sdjiV
![Page 17: Sequence Alignment](https://reader030.vdocument.in/reader030/viewer/2022020417/56814d7d550346895dbada01/html5/thumbnails/17.jpg)
Recursive Argument
Of course, we also need to handle the base cases in the recursion:
])[,(],[],[
)],[(],[],[
],[
1jtj0V1j0V
1is0iV01iV
000V
![Page 18: Sequence Alignment](https://reader030.vdocument.in/reader030/viewer/2022020417/56814d7d550346895dbada01/html5/thumbnails/18.jpg)
Dynamic Programming Algorithm
We fill the matrix using the recurrence rule
0A1
G2
C3
0
A 1
A 2
A 3
C 4
![Page 19: Sequence Alignment](https://reader030.vdocument.in/reader030/viewer/2022020417/56814d7d550346895dbada01/html5/thumbnails/19.jpg)
Dynamic Programming Algorithm
0A1
G2
C3
0 0 -2 -4 -6
A 1 -2 1 -1 -3
A 2 -4 -1 0 -2
A 3 -6 -3 -2 -1
C 4 -8 -5 -4 -1
Conclusion: d(AAAC,AGC) = -1
![Page 20: Sequence Alignment](https://reader030.vdocument.in/reader030/viewer/2022020417/56814d7d550346895dbada01/html5/thumbnails/20.jpg)
Reconstructing the Best Alignment
To reconstruct the best alignment, we record which case in the recursive rule maximized the score
0A1
G2
C3
0 0 -2 -4 -6
A 1 -2 1 -1 -3
A 2 -4 -1 0 -2
A 3 -6 -3 -2 -1
C 4 -8 -5 -4 -1
![Page 21: Sequence Alignment](https://reader030.vdocument.in/reader030/viewer/2022020417/56814d7d550346895dbada01/html5/thumbnails/21.jpg)
Reconstructing the Best Alignment
We now trace back the path the corresponds to the best alignment
0A1
G2
C3
0 0 -2 -4 -6
A 1 -2 1 -1 -3
A 2 -4 -1 0 -2
A 3 -6 -3 -2 -1
C 4 -8 -5 -4 -1
AAACAG-C
![Page 22: Sequence Alignment](https://reader030.vdocument.in/reader030/viewer/2022020417/56814d7d550346895dbada01/html5/thumbnails/22.jpg)
Reconstructing the Best Alignment
Sometimes, more than one alignment has the best score
0A1
G2
C3
0 0 -2 -4 -6
A 1 -2 1 -1 -3
A 2 -4 -1 0 -2
A 3 -6 -3 -2 -1
C 4 -8 -5 -4 -1
AAACA-GC
![Page 23: Sequence Alignment](https://reader030.vdocument.in/reader030/viewer/2022020417/56814d7d550346895dbada01/html5/thumbnails/23.jpg)
Local Alignment
Consider now a different question: Can we find similar substring of s and t Formally, given s[1..n] and t[1..m] find i,j,k, and l
such that d(s[i..j],t[k..l]) is maximal
![Page 24: Sequence Alignment](https://reader030.vdocument.in/reader030/viewer/2022020417/56814d7d550346895dbada01/html5/thumbnails/24.jpg)
Local Alignment
As before, we use dynamic programming We now want to setV[i,j] to record the best
alignment of a suffix of s[1..i] and a suffix of t[1..j]
How should we change the recurrence rule?
![Page 25: Sequence Alignment](https://reader030.vdocument.in/reader030/viewer/2022020417/56814d7d550346895dbada01/html5/thumbnails/25.jpg)
Local Alignment
New option: We can start a new match instead of extend
previous alignment
0
1jtj1iV1is1jiV
1jt1isjiV
1j1iV])[,(],[)],[(],[
])[],[(],[
max],[
Alignment of empty suffixes
]))[,(],[,max(],[
))],[(],[,max(],[
],[
1jtj0V01j0V
1is0iV001iV
000V
![Page 26: Sequence Alignment](https://reader030.vdocument.in/reader030/viewer/2022020417/56814d7d550346895dbada01/html5/thumbnails/26.jpg)
Local Alignment Example
0A1
T2
C3
T4
A5
A6
0
T 1
A 2
A 3
T 4
A 5
s = TAATAt = ATCTAA
![Page 27: Sequence Alignment](https://reader030.vdocument.in/reader030/viewer/2022020417/56814d7d550346895dbada01/html5/thumbnails/27.jpg)
Local Alignment Example
0T1
A2
C3
T4
A5
A6
0 0 0 0 0 0 0 0
T 1 0 1 0 0 1 0 0
A 2 0 0 2 0 0 2 1
A 3 0 0 1 1 0 1 3
T 4 0 0 0 0 2 0 1
A 5 0 0 1 0 0 3 1
s = TAATAt = TACTAA
![Page 28: Sequence Alignment](https://reader030.vdocument.in/reader030/viewer/2022020417/56814d7d550346895dbada01/html5/thumbnails/28.jpg)
Local Alignment Example
0T1
A2
C3
T4
A5
A6
0 0 0 0 0 0 0 0
T 1 0 1 0 0 1 0 0
A 2 0 0 2 0 0 2 1
A 3 0 0 1 1 0 1 3
T 4 0 0 0 0 2 0 1
A 5 0 0 1 0 0 3 1
s = TAATAt = TACTAA
![Page 29: Sequence Alignment](https://reader030.vdocument.in/reader030/viewer/2022020417/56814d7d550346895dbada01/html5/thumbnails/29.jpg)
Local Alignment Example
0T1
A2
C3
T4
A5
A6
0 0 0 0 0 0 0 0
T 1 0 1 0 0 1 0 0
A 2 0 0 2 0 0 2 1
A 3 0 0 1 1 0 1 3
T 4 0 0 0 0 2 0 1
A 5 0 0 1 0 0 3 1
s = TAATAt = TACTAA
![Page 30: Sequence Alignment](https://reader030.vdocument.in/reader030/viewer/2022020417/56814d7d550346895dbada01/html5/thumbnails/30.jpg)
Sequence Alignment
We seen two variants of sequence alignment: Global alignment Local alignment
Other variants: Finding best overlap (exercise)
All are based on the same basic idea of dynamic programming