{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,7,8]],"date-time":"2026-07-08T00:52:39Z","timestamp":1783471959176,"version":"3.55.0"},"reference-count":31,"publisher":"Oxford University Press (OUP)","issue":"5","license":[{"start":{"date-parts":[[2026,4,27]],"date-time":"2026-04-27T00:00:00Z","timestamp":1777248000000},"content-version":"vor","delay-in-days":1,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"funder":[{"DOI":"10.13039\/501100002347","name":"BMBF","doi-asserted-by":"publisher","id":[{"id":"10.13039\/501100002347","id-type":"DOI","asserted-by":"publisher"}]}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2026,5,3]]},"abstract":"<jats:title>Abstract<\/jats:title>\n                  <jats:sec>\n                    <jats:title>Motivation<\/jats:title>\n                    <jats:p>Substitution matrices like BLOSUM62 model the likelihood of replacement of amino acids in evolution. Substitution matrices are used in protein sequence alignment tasks. Since the introduction of BLOSUM62 over three decades ago, many matrices have been released. Yet, to date, no effort uses large amounts of 3D structures predicted by AlphaFold.<\/jats:p>\n                  <\/jats:sec>\n                  <jats:sec>\n                    <jats:title>Results<\/jats:title>\n                    <jats:p>Here, we define AFSM, the AlphaFold Substitution Matrix derived from over 20\u00a0000 predicted 3D structures following the BLOSUM methodology. We benchmark AFSM against BLOSUM62 and 16 other matrices on five tasks in multiple sequence alignment (MSA) and protein homology search. Our analysis surprisingly reveals that all matrix families perform similarly. Only when there are few sequences in an MSA do BLOSUM62 and AFSM perform better than using no matrix. This suggests that substitution matrices were most beneficial when there was little sequence data. We corroborate this argument by showing that embeddings, which are computed from billions of sequences, perform better than substitution matrices, when sequence data is sparse. Taken together, this suggests that structural data does not improve BLOSUM62. But increased sequence data makes extrapolation with substitution matrices obsolete. Nonetheless, BLOSUM62 continues to capture chemists\u2019 intuition on amino acids by providing numerical values implicitly reflecting physicochemical properties, and it remains indispensable for sparse MSAs and direct comparison of two sequences.<\/jats:p>\n                  <\/jats:sec>\n                  <jats:sec>\n                    <jats:title>Availability and implementation<\/jats:title>\n                    <jats:p>Data is available from doi.org\/10.5281\/zenodo.18777546<\/jats:p>\n                  <\/jats:sec>","DOI":"10.1093\/bioinformatics\/btag188","type":"journal-article","created":{"date-parts":[[2026,4,24]],"date-time":"2026-04-24T11:46:03Z","timestamp":1777031163000},"source":"Crossref","is-referenced-by-count":1,"title":["Much ado about nothing: modeling amino acid replacement with predicted protein structures"],"prefix":"10.1093","volume":"42","author":[{"given":"Lukas","family":"Buschmann","sequence":"first","affiliation":[{"name":"BIOTEC, CMCB, TU Dresden, Dresden 01062, Germany"},{"name":"ScaDS.AI, TU Dresden, Dresden 01062, Germany"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Sarah Naomi","family":"Bolz","sequence":"additional","affiliation":[{"name":"BIOTEC, CMCB, TU Dresden, Dresden 01062, Germany"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Ferras","family":"el-Hendi","sequence":"additional","affiliation":[{"name":"BIOTEC, CMCB, TU Dresden, Dresden 01062, Germany"},{"name":"ScaDS.AI, TU Dresden, Dresden 01062, Germany"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Negin","family":"Malekian","sequence":"additional","affiliation":[{"name":"BIOTEC, CMCB, TU Dresden, Dresden 01062, Germany"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0003-2848-6949","authenticated-orcid":false,"given":"Michael","family":"Schroeder","sequence":"additional","affiliation":[{"name":"BIOTEC, CMCB, TU Dresden, Dresden 01062, Germany"},{"name":"ScaDS.AI, TU Dresden, Dresden 01062, Germany"}],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"286","published-online":{"date-parts":[[2026,4,26]]},"reference":[{"key":"2026052023570202200_btag188-B1","doi-asserted-by":"crossref","first-page":"493","DOI":"10.1038\/s41586-024-07487-w","article-title":"Accurate structure prediction of biomolecular interactions with AlphaFold 3","volume":"630","author":"Abramson","year":"2024","journal-title":"Nature"},{"key":"2026052023570202200_btag188-B2","doi-asserted-by":"crossref","first-page":"D376","DOI":"10.1093\/nar\/gkz1064","article-title":"The SCOP database in 2020: expanded classification of representative family and superfamily domains of known protein structures","volume":"48","author":"Andreeva","year":"2020","journal-title":"Nucleic Acids Res"},{"key":"2026052023570202200_btag188-B3","doi-asserted-by":"crossref","first-page":"1323","DOI":"10.1093\/protein\/7.11.1323","article-title":"Amino acid substitution during functionally constrained divergent evolution of protein sequences","volume":"7","author":"Benner","year":"1994","journal-title":"Protein Eng"},{"key":"2026052023570202200_btag188-B4","doi-asserted-by":"crossref","first-page":"235","DOI":"10.1093\/nar\/28.1.235","article-title":"The protein data bank","volume":"28","author":"Berman","year":"2000","journal-title":"Nucleic Acids Res"},{"key":"2026052023570202200_btag188-B5","doi-asserted-by":"crossref","first-page":"321","DOI":"10.3233\/ISB-00245","article-title":"Analysis and comparison of benchmarks for multiple sequence alignment","volume":"6","author":"Blackshields","year":"2006","journal-title":"In Silico Biol"},{"key":"2026052023570202200_btag188-B6","doi-asserted-by":"crossref","first-page":"D444","DOI":"10.1093\/nar\/gkae1082","article-title":"InterPro: the protein sequence classification resource in 2025","volume":"53","author":"Blum","year":"2025","journal-title":"Nucleic Acids Res"},{"key":"2026052023570202200_btag188-B7","doi-asserted-by":"crossref","first-page":"396","DOI":"10.1186\/1471-2105-10-396","article-title":"Optimizing substitution matrix choice and gap parameters for sequence alignment","volume":"10","author":"Edgar","year":"2009","journal-title":"BMC Bioinformatics"},{"key":"2026052023570202200_btag188-B8","doi-asserted-by":"crossref","DOI":"10.1093\/bioinformatics\/btae687","article-title":"Protein structure alignment by reseek improves sensitivity to remote homologs","volume":"40","author":"Edgar","year":"2024","journal-title":"Bioinformatics"},{"key":"2026052023570202200_btag188-B9","author":"El-Hendi","year":"2026"},{"key":"2026052023570202200_btag188-B10","doi-asserted-by":"crossref","first-page":"7112","DOI":"10.1109\/TPAMI.2021.3095381","article-title":"ProtTrans: toward understanding the language of life through self-supervised learning","volume":"44","author":"Elnaggar","year":"2022","journal-title":"IEEE Trans Pattern Anal Mach Intell"},{"key":"2026052023570202200_btag188-B11","doi-asserted-by":"crossref","first-page":"862","DOI":"10.1126\/science.185.4154.862","article-title":"Amino acid difference formula to help explain protein evolution","volume":"185","author":"Grantham","year":"1974","journal-title":"Science"},{"key":"2026052023570202200_btag188-B12","doi-asserted-by":"crossref","first-page":"10915","DOI":"10.1073\/pnas.89.22.10915","article-title":"Amino acid substitution matrices from protein blocks","volume":"89","author":"Henikoff","year":"1992","journal-title":"Proc Natl Acad Sci U S A"},{"key":"2026052023570202200_btag188-B13","doi-asserted-by":"crossref","first-page":"189","DOI":"10.1186\/s12859-016-1060-3","article-title":"Addressing inaccuracies in BLOSUM computation improves homology search performance","volume":"17","author":"Hess","year":"2016","journal-title":"BMC Bioinformatics"},{"key":"2026052023570202200_btag188-B14","doi-asserted-by":"crossref","first-page":"499","DOI":"10.1002\/prot.22458","article-title":"Structure is three to ten times more conserved than sequence\u2014a study of structural response in protein cores","volume":"77","author":"Illergard","year":"2009","journal-title":"Proteins"},{"key":"2026052023570202200_btag188-B15","doi-asserted-by":"crossref","first-page":"671","DOI":"10.1002\/prot.26050","article-title":"New amino acid substitution matrix brings sequence alignments into agreement with structure matches","volume":"89","author":"Jia","year":"2021","journal-title":"Proteins"},{"key":"2026052023570202200_btag188-B16","doi-asserted-by":"crossref","first-page":"583","DOI":"10.1038\/s41586-021-03819-2","article-title":"Highly accurate protein structure prediction with AlphaFold","volume":"596","author":"Jumper","year":"2021","journal-title":"Nature"},{"key":"2026052023570202200_btag188-B17","doi-asserted-by":"crossref","first-page":"293","DOI":"10.1186\/s12859-017-1703-z","article-title":"PFASUM: a substitution matrix from Pfam structural alignments","volume":"18","author":"Keul","year":"2017","journal-title":"BMC Bioinformatics"},{"key":"2026052023570202200_btag188-B18","doi-asserted-by":"crossref","first-page":"736","DOI":"10.1186\/s40064-016-2526-5","article-title":"Determination of optimal parameters of MAFFT program based on BAliBASE3.0 database","volume":"5","author":"Long","year":"2016","journal-title":"Springerplus"},{"key":"2026052023570202200_btag188-B19","doi-asserted-by":"crossref","first-page":"D412","DOI":"10.1093\/nar\/gkaa913","article-title":"Pfam: the protein families database in 2021","volume":"49","author":"Mistry","year":"2021","journal-title":"Nucleic Acids Res"},{"key":"2026052023570202200_btag188-B20","doi-asserted-by":"crossref","first-page":"8","DOI":"10.1093\/oxfordjournals.molbev.a003985","article-title":"Estimating amino acid substitution models: a comparison of Dayhoff\u2019s estimator, the resolvent approach and a maximum likelihood method","volume":"19","author":"M\u00fcller","year":"2002","journal-title":"Mol Biol Evol"},{"key":"2026052023570202200_btag188-B21","doi-asserted-by":"crossref","first-page":"761","DOI":"10.1089\/10665270050514918","article-title":"Modeling amino acid replacement","volume":"7","author":"M\u00fcller","year":"2000","journal-title":"J Comput Biol"},{"key":"2026052023570202200_btag188-B22","doi-asserted-by":"crossref","DOI":"10.1093\/bioinformatics\/btad786","article-title":"Embedding-based alignment: combining protein language models with dynamic programming alignment to detect structural similarities in the twilight-zone","volume":"40","author":"Pantolini","year":"2024","journal-title":"Bioinformatics"},{"key":"2026052023570202200_btag188-B23","doi-asserted-by":"crossref","first-page":"452","DOI":"10.1093\/nar\/gkg062","article-title":"The CATH database: an extended protein family resource for structural and functional genomics","volume":"31","author":"Pearl","year":"2003","journal-title":"Nucleic Acids Res"},{"key":"2026052023570202200_btag188-B24","doi-asserted-by":"crossref","first-page":"3824","DOI":"10.1093\/bioinformatics\/bti627","article-title":"Statistical evaluation of pairwise protein sequence comparison with the Bayesian bootstrap","volume":"21","author":"Price","year":"2005","journal-title":"Bioinformatics"},{"key":"2026052023570202200_btag188-B25","doi-asserted-by":"crossref","first-page":"545","DOI":"10.1093\/protein\/13.8.545","article-title":"Structure-derived substitution matrices for alignment of distantly related sequences","volume":"13","author":"Prlic","year":"2000","journal-title":"Protein Eng"},{"key":"2026052023570202200_btag188-B26","doi-asserted-by":"crossref","first-page":"85","DOI":"10.1093\/protein\/12.2.85","article-title":"Twilight zone of protein sequence alignments","volume":"12","author":"Rost","year":"1999","journal-title":"Protein Eng"},{"key":"2026052023570202200_btag188-B27","doi-asserted-by":"crossref","first-page":"e79854","DOI":"10.7554\/eLife.79854","article-title":"Generative power of a protein language model trained on multiple sequence alignments","volume":"12","author":"Sgarbossa","year":"2023","journal-title":"eLife"},{"key":"2026052023570202200_btag188-B28","doi-asserted-by":"crossref","first-page":"274","DOI":"10.1038\/nbt0308-274","article-title":"BLOSUM62 miscalculations improve search performance","volume":"26","author":"Styczynski","year":"2008","journal-title":"Nat Biotechnol"},{"key":"2026052023570202200_btag188-B29","doi-asserted-by":"crossref","first-page":"e18093","DOI":"10.1371\/journal.pone.0018093","article-title":"A comprehensive benchmark study of multiple sequence alignment methods: current challenges and future perspectives","volume":"6","author":"Thompson","year":"2011","journal-title":"PLoS One"},{"key":"2026052023570202200_btag188-B30","doi-asserted-by":"crossref","first-page":"D439","DOI":"10.1093\/nar\/gkab1061","article-title":"AlphaFold protein structure database: massively expanding the structural coverage of protein- sequence space with high-accuracy models","volume":"50","author":"Varadi","year":"2022","journal-title":"Nucleic Acids Res"},{"key":"2026052023570202200_btag188-B31","doi-asserted-by":"crossref","first-page":"D368","DOI":"10.1093\/nar\/gkad1011","article-title":"AlphaFold protein structure database in 2024: providing structure coverage for over 214 million protein sequences","volume":"52","author":"Varadi","year":"2024","journal-title":"Nucleic Acids Res"}],"container-title":["Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/academic.oup.com\/bioinformatics\/advance-article-pdf\/doi\/10.1093\/bioinformatics\/btag188\/68162155\/btag188.pdf","content-type":"application\/pdf","content-version":"am","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/42\/5\/btag188\/68162155\/btag188.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/42\/5\/btag188\/68162155\/btag188.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2026,5,21]],"date-time":"2026-05-21T03:57:15Z","timestamp":1779335835000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article\/doi\/10.1093\/bioinformatics\/btag188\/8662821"}},"subtitle":[],"editor":[{"given":"Can","family":"Alkan","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"editor"}]}],"short-title":[],"issued":{"date-parts":[[2026,4,26]]},"references-count":31,"journal-issue":{"issue":"5","published-print":{"date-parts":[[2026,5,3]]}},"URL":"https:\/\/doi.org\/10.1093\/bioinformatics\/btag188","relation":{},"ISSN":["1367-4803","1367-4811"],"issn-type":[{"value":"1367-4803","type":"print"},{"value":"1367-4811","type":"electronic"}],"subject":[],"published-other":{"date-parts":[[2026,5]]},"published":{"date-parts":[[2026,4,26]]},"article-number":"btag188"}}