NCBI
Home
Up

The .doc format is much clearer.

BLOSUM62 Substitution Matrix

 

  C S T P A G N D E Q H R K M I L V F Y W
C

9

-1

-1

-3

0

-3

-3

-3

-4

-3

-3

-3

-3

-1

-1

-1

-1

-2

-2

-2

S

-1

4

1

-1

1

0

1

0

0

0

-1

-1

0

-1

-2

-2

-2

-2

-2

-3

T

-1

1

4

1

-1

1

0

1

0

0

0

-1

0

-1

-2

-2

-2

-2

-2

-3

P

-3

-1

1

7

-1

-2

-1

-1

-1

-1

-2

-2

-1

-2

-3

-3

-2

-4

-3

-4

A

0

1

-1

-1

4

0

-1

-2

-1

-1

-2

-1

-1

-1

-1

-1

-2

-2

-2

-3

G

-3

0

1

-2

0

6

-2

-1

-2

-2

-2

-2

-2

-3

-4

-4

0

-3

-3

-2

N

-3

1

0

-2

-2

0

6

1

0

0

-1

0

0

-2

-3

-3

-3

-3

-2

-4

D

-3

0

1

-1

-2

-1

1

6

2

0

-1

-2

-1

-3

-3

-4

-3

-3

-3

-4

E

-4

0

0

-1

-1

-2

0

2

5

2

0

0

1

-2

-3

-3

-3

-3

-2

-3

Q

-3

0

0

-1

-1

-2

0

0

2

5

0

1

1

0

-3

-2

-2

-3

-1

-2

H

-3

-1

0

-2

-2

-2

1

1

0

0

8

0

-1

-2

-3

-3

-2

-1

2

-2

R

-3

-1

-1

-2

-1

-2

0

-2

0

1

0

5

2

-1

-3

-2

-3

-3

-2

-3

K

-3

0

0

-1

-1

-2

0

-1

1

1

-1

2

5

-1

-3

-2

-3

-3

-2

-3

M

-1

-1

-1

-2

-1

-3

-2

-3

-2

0

-2

-1

-1

5

1

2

-2

0

-1

-1

I

-1

-2

-2

-3

-1

-4

-3

-3

-3

-3

-3

-3

-3

1

4

2

1

0

-1

-3

L

-1

-2

-2

-3

-1

-4

-3

-4

-3

-2

-3

-2

-2

2

2

4

3

0

-1

-2

V

-1

-2

-2

-2

0

-3

-3

-3

-2

-2

-3

-3

-2

1

3

1

4

-1

-1

-3

F

-2

-2

-2

-4

-2

-3

-3

-3

-3

-3

-1

-3

-3

0

0

0

-1

6

3

1

Y

-2

-2

-2

-3

-2

-3

-2

-3

-2

-1

2

-2

-2

-1

-1

-1

-1

3

7

2

W

-2

-3

-3

-4

-3

-2

-4

-4

-3

-2

-2

-3

-3

-1

-3

-2

-3

1

2

11

 

This table depicts the observed substitutions found in a broad sampling from the aligned segments of polypeptides. The precise method of calculation gets a bit abstruse, but an illustration may demystify some aspects.

The procedure for calculating a BLOSUM matrix is based on a likelihood method estimating the occurrence of each possible pairwise substitution. A very simple illustration of such a calculation for a very short segment is given below to illustrate the process. The polypeptides were initially aligned using an identity scoring matrix. Only aligned blocks are used to calculate the BLOSUMs.

Assume that the following 5 proteins were aligned as follows:

AVAAA

AVAAA

AVAAA

AVLAA

VVAAL

First, we have to decide whether each sequence should count equally in this process. If this database were a global representation, they should. However, databases of today typically over-represent certain classes of proteins. Therefore, the first step is to reduce this over-representation and make the dataset more representative. One method of doing this is to count all the identical blocks as if they were a single block, reducing the aligned database to:

AVAAA

AVLAA

VVAAL

Next, at each position, we must calculate the observed and expected pairwise substitutions. At position 1, we have:

A

A A-A, A-V, A-V are the observed substitutions. Without specifying the necessary calculations, one

V can see that the substitution of A for A or A for V are quite likely and that A-X is not!

Why is BLOSUM62 called BLOSUM62? Basically, this is because all blocks whose members shared at least 62% identity with ANY other member of that block were averaged and represented as 1 sequence.

How would this work with our example?

AVAAA

AVAAA

AVAAA

AVLAA

VVAAL

1-4 are all 80% identical to each other. Sequence 5 is less than 62% (it is 60% identical to 1-3 and t0 #4).

This means that the BLOCK used to make a BLOSUM62 would be: (note the averaging!)

AVA(3/4)AA

 

L(1/4)

VVA*****AL

How does the matrix relate to structural similarity? Most biochemists group the amino acids as follows:

G,A,V,L,I, M aliphatic (though some would not include G)

S,T,C hydroxyl, sulfhydryl, polar

N,Q amide side chains

F,W,Y aromatic

H,K,R basic

D,E acidic

Some rather anomalous substitutions relative to these groupings are highlighted below. For example, it seems VERY surprising to me that K-E substitution is not unusual. That changes charge! This suggests that what evolution thinks is "similar" is not necessarily similar to the molecular biologist.

(Of course, it is also possible that K-E substitution in one position is also often correlated with a E-K substitution elsewhere, and that what one retains is the electrostatic interaction between charged residues. Much like G-C or C-G compensating mutations in stem-loops of RNA secondary structures.)

  C S T P A G N D E Q H R K M I L V F Y W
C

9

-1

-1

-3

0

-3

-3

-3

-4

-3

-3

-3

-3

-1

-1

-1

-1

-2

-2

-2

S

-1

4

1

-1

1

0

1

0

0

0

-1

-1

0

-1

-2

-2

-2

-2

-2

-3

T

-1

1

4

1

-1

1

0

1

0

0

0

-1

0

-1

-2

-2

-2

-2

-2

-3

P

-3

-1

1

7

-1

-2

-1

-1

-1

-1

-2

-2

-1

-2

-3

-3

-2

-4

-3

-4

A

0

1

-1

-1

4

0

-1

-2

-1

-1

-2

-1

-1

-1

-1

-1

-2

-2

-2

-3

G

-3

0

1

-2

0

6

-2

-1

-2

-2

-2

-2

-2

-3

-4

-4

0

-3

-3

-2

N

-3

1

0

-2

-2

0

6

1

0

0

1

0

0

-2

-3

-3

-3

-3

-2

-4

D

-3

0

1

-1

-2

-1

1

6

2

0

-1

-2

-1

-3

-3

-4

-3

-3

-3

-4

E

-4

0

0

-1

-1

-2

0

2

5

2

0

0

1

-2

-3

-3

-3

-3

-2

-3

Q

-3

0

0

-1

-1

-2

0

0

2

5

0

1

1

0

-3

-2

-2

-3

-1

-2

H

-3

-1

0

-2

-2

-2

1

-1

0

0

8

0

-1

-2

-3

-3

-2

-1

2

-2

R

-3

-1

-1

-2

-1

-2

0

-2

0

1

0

5

2

-1

-3

-2

-3

-3

-2

-3

K

-3

0

0

-1

-1

-2

0

-1

1

1

-1

2

5

-1

-3

-2

-3

-3

-2

-3

M

-1

-1

-1

-2

-1

-3

-2

-3

-2

0

-2

-1

-1

5

1

2

-2

0

-1

-1

I

-1

-2

-2

-3

-1

-4

-3

-3

-3

-3

-3

-3

-3

1

4

2

1

0

-1

-3

L

-1

-2

-2

-3

-1

-4

-3

-4

-3

-2

-3

-2

-2

2

2

4

3

0

-1

-2

V

-1

-2

-2

-2

0

-3

-3

-3

-2

-2

-3

-3

-2

1

3

1

4

-1

-1

-3

F

-2

-2

-2

-4

-2

-3

-3

-3

-3

-3

-1

-3

-3

0

0

0

-1

6

3

1

Y

-2

-2

-2

-3

-2

-3

-2

-3

-2

-1

2

-2

-2

-1

-1

-1

-1

3

7

2

W

-2

-3

-3

-4

-3

-2

-4

-4

-3

-2

-2

-3

-3

-1

-3

-2

-3

1

2

11

 

University of KentuckyMorgan School of Biological SciencesNSF-CCD Support wpe1.jpg (5798 bytes)Chuck Staben, copyright reserved || 09/25/98