<p>Sean Eddy wrote an excellent explanation of the BLOSUM Matrix in his Nature Biotechnology paper <a href="https://www.nature.com/articles/nbt0804-1035" rel="nofollow noreferrer">Where did the BLOSUM62 alignment score matrix come from?</a> He also gives real examples of how one would calculate the score for different amino acid pairs and where the numbers come from.  Further references given at the end of the paper might also be useful reading, specifically: <a href="https://www.sciencedirect.com/science/article/pii/002228369190193A?via%3Dihub" rel="nofollow noreferrer">Amino acid substitution matrices from an information theoretic perspective</a>, <a href="https://www.pnas.org/content/87/6/2264.long" rel="nofollow noreferrer">Methods for assessing the statistical significance of molecular sequence features by using general scoring schemes</a>, and <a href="https://www.pnas.org/content/89/22/10915.long" rel="nofollow noreferrer">Amino acid substitution matrices from protein blocks</a>.</p>

CC BY-SA 4.0
2018-12-10T16:51:04.410Z
Could you please add a summary of what these links say to your answer? Answers [need to be able to stand on their own](https://meta.stackexchange.com/q/225370/203101) even if the inks are removed.
terdon
CC BY-SA 4.0
2018-12-11T01:54:30.137Z
Hi @Matt Bashton, I can't access the reading materials in the links as It needs subscription to read it.
Heeh wei cheng
CC BY-SA 4.0
2018-12-10T10:46:12.737Z
2018-12-11T07:30:19.643Z
2018-12-11T07:30:19.643Z
Matt Bashton
<p>BLOSUM stands for BLOcks SUbstitution Matrices (Henikoff &amp; Henikoff, 1992), and were created by observing substitution frequencies in local ungapped multiple sequence alignments. The score reflects the chance (log-odds) one amino acid is substituted for another in a set of protein multiple sequence alignments. BLOSUM scoring matrices are normally followed by a number eg BLOSUM62. The number is % sequence identity between the sequences in the multiple sequence alignment (MSA) used to create the score matrix. eg BLOSUM62 matrices were created from multiple sequence alignments with blocks that shared 62% identity. Note because these are ungapped alignments, the multiple sequence alignments were generated from conserved protein families.  With the older PAM matrices, the number is opposite, a higher number is greater distance. PAM matrices reflect evolutionarily distant in global (aka gapped) alignments of the full-length sequences that include both conserved and diverged regions.  See Wikipedia <a href="https://en.wikipedia.org/wiki/BLOSUM" rel="nofollow noreferrer">https://en.wikipedia.org/wiki/BLOSUM</a>. There are several newer substitution matrices.  The R package SubVis allows you to compare/visualize different substitution matrices. </p>

CC BY-SA 4.0
2019-01-13T07:29:40.433Z
Hi @Aedin, thanks for your input. Is it possible to provide an example to run through the formula itself and also showing how the matrix is formed as shown in the Blosum matrix in the wikipedia ?
Heeh wei cheng
CC BY-SA 4.0
2019-01-08T22:28:15.670Z
2019-01-08T22:28:15.670Z
Aedin
<p>Amino acid matrices are contructed by maximum likelihood and easily performed using PAML. You present a known tree and an amino acid alignment and ask what amino acid matrix would maximise the likelihood. Hence all amino acid matrices using this old school approach are dependent on the alignment composition. JTT was generally regarded as obtaining the highest likelihood. However, I reitereate that things have moved on a lot since this 2-D reversible amino acid substitution matrix was devised.</p>

CC BY-SA 4.0
2018-12-28T20:39:15.773Z
2018-12-28T20:39:15.773Z
M__
<p>I would like to ask how BLOSUM matrix is constructed and calculated ? I read the wikipedia and I still do not understand it. I do not understand the mathematical calculations as I have low knowledge about logarithms. </p>

<p>Maybe could someone give a simple example and then a real-scenario example of how the BLOSUM matrix can be used and calculated, maybe in relation to the Blast algorithm where it uses the BLOSUM scoring matrix to determine high scoring words for each word in the query sequence ?</p>

<p>Also, why are they different types of BLOSUM matrices like BLOSUM62 and so on ?</p>

<phylogenetics><blast>
How BLOSUM Matrix is constructed and calculated
# How BLOSUM Matrix is constructed and calculated

(Asked by: username_0 on 2018-12-10)

I would like to ask how BLOSUM matrix is constructed and calculated ? I read the wikipedia and I still do not understand it. I do not understand the mathematical calculations as I have low knowledge about logarithms.

Maybe could someone give a simple example and then a real-scenario example of how the BLOSUM matrix can be used and calculated, maybe in relation to the Blast algorithm where it uses the BLOSUM scoring matrix to determine high scoring words for each word in the query sequence ?

Also, why are they different types of BLOSUM matrices like BLOSUM62 and so on ?

**Tags:** phylogenetics, blast

---

## 3 answers

### Answer by username_1 on 2018-12-10. Score: 3

Sean Eddy wrote an excellent explanation of the BLOSUM Matrix in his Nature Biotechnology paper [Where did the BLOSUM62 alignment score matrix come from?](https://www.nature.com/articles/nbt0804-1035) He also gives real examples of how one would calculate the score for different amino acid pairs and where the numbers come from. Further references given at the end of the paper might also be useful reading, specifically: [Amino acid substitution matrices from an information theoretic perspective](https://www.sciencedirect.com/science/article/pii/002228369190193A?via%3Dihub), [Methods for assessing the statistical significance of molecular sequence features by using general scoring schemes](https://www.pnas.org/content/87/6/2264.long), and [Amino acid substitution matrices from protein blocks](https://www.pnas.org/content/89/22/10915.long).

> **Comments:**
> **username_2**: Could you please add a summary of what these links say to your answer? Answers [need to be able to stand on their own](https://meta.stackexchange.com/q/225370/203101) even if the inks are removed.  
> **username_0**: Hi @Matt Bashton, I can't access the reading materials in the links as It needs subscription to read it.  

---

### Answer by username_3 on 2019-01-08. Score: 2

BLOSUM stands for BLOcks SUbstitution Matrices (Henikoff & Henikoff, 1992), and were created by observing substitution frequencies in local ungapped multiple sequence alignments. The score reflects the chance (log-odds) one amino acid is substituted for another in a set of protein multiple sequence alignments. BLOSUM scoring matrices are normally followed by a number eg BLOSUM62. The number is % sequence identity between the sequences in the multiple sequence alignment (MSA) used to create the score matrix. eg BLOSUM62 matrices were created from multiple sequence alignments with blocks that shared 62% identity. Note because these are ungapped alignments, the multiple sequence alignments were generated from conserved protein families. With the older PAM matrices, the number is opposite, a higher number is greater distance. PAM matrices reflect evolutionarily distant in global (aka gapped) alignments of the full-length sequences that include both conserved and diverged regions. See Wikipedia <https://en.wikipedia.org/wiki/BLOSUM>. There are several newer substitution matrices. The R package SubVis allows you to compare/visualize different substitution matrices.

> **Comments:**
> **username_0**: Hi @Aedin, thanks for your input. Is it possible to provide an example to run through the formula itself and also showing how the matrix is formed as shown in the Blosum matrix in the wikipedia ?  

---

### Answer by username_4 on 2018-12-28. Score: 1

Amino acid matrices are contructed by maximum likelihood and easily performed using PAML. You present a known tree and an amino acid alignment and ask what amino acid matrix would maximise the likelihood. Hence all amino acid matrices using this old school approach are dependent on the alignment composition. JTT was generally regarded as obtaining the highest likelihood. However, I reitereate that things have moved on a lot since this 2-D reversible amino acid substitution matrix was devised.


---

