Home Page Schedule Students Assignments Links

Structural Biochemistry and Bioinformatics


Dynamic programming

Now you should leave the computer and make sure you understand how dynamic programming for sequence alignments work.

Longest common substring

Simple examples

We have the following sequences: AKWIL and AWLI how do they align using:
  1. Local alignments with an identity matrix and no gap penalties
  2. Global alignments with an identity matrix and no gap penalties
  3. Local alignments with Blossum62 and no gap penalties
  4. Global alignments with Blossum62 and no gap penalties
  5. Local alignments with Blossum62 and gap penalties -10,-1
  6. Global alignments with Blossum62 and gap penalties -10,-1

Example of real sequences

Azurin (2aza) and plastocyanin (2pcy ) are related. You should use the Blossum62 matrix and calculate the alignment of the first 10 residues of these two sequences using the blossom62 substitution matrix and the gap penalties, -1 and -10. You should try both local and global alignments. In the lab report you should show the alignment and alignment matrix (including arrows).

Resources

However if you do want to use a computer there are some nice tools available:

  1. DP table

Arne Elofsson
Last modified: Tue Oct 30 18:09:09 CET 2001
Arne Elofsson
Department of Biochemistry,
Arrheniuslaboratoriet
Stockholms Universitet
10691 Stockholm, Sweden
Tel: +46-(0)8/161553
Fax: +46-(0)8/153679
Hem: +46-(0)8/6413158
Email: arne@sbc.su.se
WWW: /~arne/