{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,9,23]],"date-time":"2026-09-23T15:07:54Z","timestamp":1790176074482,"version":"4.1.0"},"reference-count":39,"publisher":"Springer Science and Business Media LLC","license":[{"start":{"date-parts":[[2026,9,23]],"date-time":"2026-09-23T00:00:00Z","timestamp":1790121600000},"content-version":"tdm","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0"},{"start":{"date-parts":[[2026,9,23]],"date-time":"2026-09-23T00:00:00Z","timestamp":1790121600000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0"}],"funder":[{"name":"Centro de Biolog\u00eda Molecular Severo Ochoa","id":[{"id":"https:\/\/ror.org\/03v9e8t09","id-type":"ROR","asserted-by":"crossref"}]}],"content-domain":{"domain":["link.springer.com"],"crossmark-restriction":false},"short-container-title":["J Mol Evol"],"abstract":"<jats:title>Abstract<\/jats:title>\n                  <jats:p>Statistical couplings between protein sites have been enormously useful for protein structure prediction, but there is still debate on the selective forces that generate them. Here we explicitly predict the couplings considering selection for protein folding stability, both against unfolding and misfolding, which we model through the Stability constrained model of protein evolution (SCPE). The SCPE is based on contact interactions, and it predicts stability against misfolding through the Random Energy Model in the space of contact matrices. We adopt the minimum selection principle, assuming that the multivariate amino acid distribution has minimum Kullback-Leibler divergence with respect to a site-unspecific distribution that is not influenced by protein folding stability. We predict the couplings at first order, assuming that they are small (small coupling approximation, SCA), and obtain an explicit formula that relates them with the contact interaction matrix and with the structural properties of the pairs of sites: presence of native contact, distance along the sequence and buriedness of the sites. We test these predictions on a representative set of the Protein Data Bank. For pairs of sites in contact, both predicted and observed couplings are inversely related with the contact energies of the amino acid pairs, while for pairs not in contact this relation is positive, i.e. negative design destabilizes pairs of amino acids for short-range pairs of sites that can form wrong contacts with high probability. Our results suggest that the strongest couplings correspond to native contacts, but an even larger number of couplings are influenced by negative design. Intriguingly, the most informative pairs are formed by surface sites, suggesting a role for protein-protein interactions and functional motions. We interpret the site-unspecific couplings, which are not influenced by selection for protein folding stability, as the result of variations of global mutational or selective forces across the lineages where the proteins evolve. Consistently, the strongest global couplings are self-couplings between equal amino acids. Cys-Cys, which form disulfide bridges, show the strongest global coupling, followed by pairs of metabolically costly amino acids (Trp-Trp, Tyr-Tyr, Trp-Tyr) and positively charged amino acids that tend to interact with nucleic acids (Lys-Lys, Arg-Arg). Ser, Thr, Asn and Gln form a cluster of globally correlated polar amino acids.<\/jats:p>","DOI":"10.1007\/s00239-026-10348-x","type":"journal-article","created":{"date-parts":[[2026,9,23]],"date-time":"2026-09-23T14:30:06Z","timestamp":1790173806000},"update-policy":"https:\/\/doi.org\/10.1007\/springer_crossmark_policy","source":"Crossref","is-referenced-by-count":0,"title":["Selection for Folding Stability Predicts Observed Covariation Between Protein Positions in the PDB"],"prefix":"10.1007","author":[{"given":"Fatemeh","family":"Saebi","sequence":"first","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Jonas","family":"Minning","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Ugo","family":"Bastolla","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"297","published-online":{"date-parts":[[2026,9,23]]},"reference":[{"issue":"8","key":"10348_CR1","doi-asserted-by":"publisher","first-page":"2195","DOI":"10.1093\/molbev\/msv085","volume":"32","author":"M Arenas","year":"2015","unstructured":"Arenas M, Sanchez-Cobos A, Bastolla U (2015) Maximum-likelihood phylogenetic inference with selection on protein folding stability. Mol Biol Evol 32(8):2195\u20132207","journal-title":"Mol Biol Evol"},{"issue":"1","key":"10348_CR2","doi-asserted-by":"publisher","DOI":"10.1088\/1478-3975\/6\/1\/015001","volume":"6","author":"D Barrick","year":"2009","unstructured":"Barrick D (2009) What have we learned from the studies of two-state folders, and what are the unanswered questions about two-state protein folding? Phys Biol 6(1):015001","journal-title":"Phys Biol"},{"key":"10348_CR3","doi-asserted-by":"publisher","first-page":"59","DOI":"10.1016\/j.sbi.2016.10.020","volume":"42","author":"U Bastolla","year":"2017","unstructured":"Bastolla U, Dehouck Y, Echave J (2017) What evolution tells us about protein physics, and protein physics tells us about evolution. Curr Opin Struct Biol 42:59\u201366","journal-title":"Curr Opin Struct Biol"},{"key":"10348_CR4","doi-asserted-by":"publisher","first-page":"405","DOI":"10.1093\/protein\/gzi045","volume":"18","author":"U Bastolla","year":"2005","unstructured":"Bastolla U, Demetrius L (2005) Stability constraints and protein evolution: the role of chain length, composition and disulfide bonds. Prot Eng Des Sel 18:405\u2013415","journal-title":"Prot Eng Des Sel"},{"key":"10348_CR5","doi-asserted-by":"publisher","DOI":"10.1186\/1471-2148-6-43","volume":"6","author":"U Bastolla","year":"2006","unstructured":"Bastolla U, Porto M, Roman HE, Vendruscolo M (2006) A protein evolution model with independent sites that reproduces site-specific amino acid distributions from the Protein Data Bank. BMC Evol Biol 6:43","journal-title":"BMC Evol Biol"},{"issue":"8","key":"10348_CR6","doi-asserted-by":"publisher","first-page":"3977","DOI":"10.1073\/pnas.97.8.3977","volume":"97","author":"U Bastolla","year":"2000","unstructured":"Bastolla U, Vendruscolo M, Knapp EW (2000) A statistical mechanical method to optimize energy functions for protein folding. Proc Nat Acad Sci 97(8):3977\u20133981","journal-title":"Proc Nat Acad Sci"},{"key":"10348_CR7","doi-asserted-by":"publisher","first-page":"179","DOI":"10.2307\/2987782","volume":"24","author":"J Besag","year":"1975","unstructured":"Besag J (1975) Statistical Analysis of Non-Lattice Data. The Statistician 24:179\u2013195","journal-title":"The Statistician"},{"key":"10348_CR8","doi-asserted-by":"publisher","first-page":"249","DOI":"10.1038\/nrg3414","volume":"14","author":"D De Juan","year":"2013","unstructured":"De Juan D, Pazos F, Valencia A (2013) Emerging methods in protein co-evolution. Nat Rev Genet 14:249\u2013261","journal-title":"Nat Rev Genet"},{"key":"10348_CR9","doi-asserted-by":"publisher","first-page":"627","DOI":"10.1039\/c7ib00079k","volume":"9","author":"Y Dehouck","year":"2017","unstructured":"Dehouck Y, Bastolla U (2017) The maximum penalty criterion for ridge regression: application to the calibration of the force constant in elastic network models. Integr Biol (Camb) 9:627\u2013641","journal-title":"Integr Biol (Camb)"},{"key":"10348_CR10","doi-asserted-by":"publisher","first-page":"2613","DOI":"10.1103\/PhysRevB.24.2613","volume":"24","author":"B Derrida","year":"1981","unstructured":"Derrida B (1981) Random-energy model: an exactly solvable model of disordered systems. Phys Rev B 24:2613","journal-title":"Phys Rev B"},{"key":"10348_CR12","doi-asserted-by":"publisher","first-page":"413","DOI":"10.1016\/j.cplett.2008.04.042","volume":"457","author":"J Echave","year":"2008","unstructured":"Echave J (2008) Evolutionary divergence of protein structure: The linearly forced elastic network model. Chem Phys Lett 457:413\u2013416","journal-title":"Chem Phys Lett"},{"key":"10348_CR13","doi-asserted-by":"publisher","first-page":"109","DOI":"10.1038\/nrg.2015.18","volume":"17","author":"J Echave","year":"2016","unstructured":"Echave J, Spielman SJ, Wilke CO (2016) Causes of evolutionary rate variation among protein sites. Nat Rev Genet 17:109\u2013121","journal-title":"Nat Rev Genet"},{"key":"10348_CR14","doi-asserted-by":"publisher","DOI":"10.1103\/PhysRevE.87.012707","volume":"87","author":"M Ekeberg","year":"2013","unstructured":"Ekeberg M, L\u00f6vkvist C, Lan Y, Weigt M, Aurell E (2013) Improved contact prediction in proteins: using pseudolikelihoods to infer Potts models. Phys Rev E Stat Nonlin Soft Matter Phys 87:012707","journal-title":"Phys Rev E Stat Nonlin Soft Matter Phys"},{"key":"10348_CR16","first-page":"205","volume":"50","author":"R Fisher","year":"1930","unstructured":"Fisher R (1930) The distribution of gene ratios for rare mutations. Proc R Soc Edin 50:205\u2013220","journal-title":"Proc R Soc Edin"},{"key":"10348_CR17","doi-asserted-by":"publisher","first-page":"309","DOI":"10.1002\/prot.340180402","volume":"18","author":"U G\u00f6bel","year":"1994","unstructured":"G\u00f6bel U, Sander C, Schneider R, Valencia A (1994) Correlated mutations and residue contacts in proteins. Proteins 18:309\u2013317","journal-title":"Proteins"},{"key":"10348_CR18","doi-asserted-by":"publisher","first-page":"170","DOI":"10.1016\/j.sbi.2008.01.006","volume":"18","author":"RA Goldstein","year":"2008","unstructured":"Goldstein RA (2008) The structure of protein evolution and the evolution of protein structure. Curr op struct biol 18:170\u2013177","journal-title":"Curr Op Struct Biol"},{"key":"10348_CR19","doi-asserted-by":"publisher","DOI":"10.7554\/eLife.08932","volume":"4","author":"E Jacob","year":"2015","unstructured":"Jacob E, Unger R, Horovitz A (2015) Codon-level information improves predictions of inter-residue contacts in proteins by correlated mutation analysis. Elife 4:e08932","journal-title":"Elife"},{"key":"10348_CR21","doi-asserted-by":"publisher","first-page":"184","DOI":"10.1093\/bioinformatics\/btr638","volume":"28","author":"DT Jones","year":"2012","unstructured":"Jones DT, Buchan DW, Cozzetto D, Pontil M (2012) PSICOV: precise structural contact prediction using sparse inverse covariance estimation on large multiple sequence alignments. Bioinformatics 28:184\u2013190","journal-title":"Bioinformatics"},{"key":"10348_CR22","doi-asserted-by":"publisher","first-page":"583","DOI":"10.1038\/s41586-021-03819-2","volume":"596","author":"J Jumper","year":"2021","unstructured":"Jumper J, Evans R, Pritzel A, Green T, Figurnov M, Ronneberger O, Hassabis D (2021) Highly accurate protein structure prediction with AlphaFold. Nature 596:583\u2013589","journal-title":"Nature"},{"key":"10348_CR23","doi-asserted-by":"publisher","first-page":"15674","DOI":"10.1073\/pnas.1314045110","volume":"110","author":"H Kamisetty","year":"2013","unstructured":"Kamisetty H, Ovchinnikov S, Baker D (2013) Assessing the utility of coevolution-based residue-residue contact predictions in a sequence-and structure-rich era. Proc Nat Acad Sci 110:15674\u201315679","journal-title":"Proc Nat Acad Sci"},{"key":"10348_CR24","doi-asserted-by":"publisher","first-page":"12524","DOI":"10.1073\/pnas.96.22.12524","volume":"96","author":"P Koehl","year":"1999","unstructured":"Koehl P, Levitt M (1999) Structure-based conformational preferences of amino acids. Proc Natl Acad Sci U S A 96:12524\u20139","journal-title":"Proc Natl Acad Sci U S A"},{"key":"10348_CR44","doi-asserted-by":"crossref","unstructured":"Kullback S, Leibler RA (1951) On information and sufficiency. Ann Math Stat 22:79\u201386","DOI":"10.1214\/aoms\/1177729694"},{"key":"10348_CR26","unstructured":"Lapedes A, Giraud B, Jarzynski C (2002) Using sequence alignments to predict protein structure and stability with high accuracy. Preprint available at arXiv 1207:2484"},{"key":"10348_CR27","doi-asserted-by":"publisher","first-page":"769","DOI":"10.1002\/pro.2071","volume":"21","author":"DA Liberles","year":"2012","unstructured":"Liberles DA, Teichmann SA, Bahar I, Bastolla U, Bloom J, Bornberg-Bauer E, Colwell LJ, de Koning AP, Dokholyan NV, Echave J, Elofsson A, Gerloff DL, Goldstein RA, Grahnen JA, Holder MT, Lakner C, Lartillot N, Lovell SC, Naylor G, Perica T, Pollock DD, Pupko T, Regan L, Roger A, Rubinstein N, Shakhnovich E, Sj\u00f6lander K, Sunyaev S, Teufel AI, Thorne JL, Thornton JW, Weinreich DM, Whelan S (2012) The interface of protein structure, protein biophysics, and molecular evolution. Protein Sci 21:769\u2013785","journal-title":"Protein Sci"},{"key":"10348_CR28","doi-asserted-by":"crossref","unstructured":"Lorca-Alonso I, Otero De-Navascues F, Arenas M, Bastolla U (2025). Site-specific structure and stability constrained substitution models improve phylogenetic inference. Syst Biol: syaf007","DOI":"10.1093\/sysbio\/syaf007"},{"key":"10348_CR29","doi-asserted-by":"publisher","DOI":"10.1371\/journal.pone.0028766","volume":"6","author":"DS Marks","year":"2011","unstructured":"Marks DS, Colwell LJ, Sheridan R, Hopf TA, Pagnani A, Zecchina R, Sander C (2011) Protein 3D structure computed from evolutionary sequence variation. PLoS One 6:e28766","journal-title":"PLoS One"},{"key":"10348_CR30","doi-asserted-by":"publisher","first-page":"1072","DOI":"10.1038\/nbt.2419","volume":"30","author":"DS Marks","year":"2012","unstructured":"Marks DS, Hopf TA (2012) Sander C (2012) Protein structure prediction from sequence variation. Nat Biotechnol 30:1072\u201380","journal-title":"Nat Biotechnol"},{"key":"10348_CR31","unstructured":"Minning J (2012) Correlations in thermodynamics and evolution of proteins (Doctoral dissertation)"},{"key":"10348_CR32","doi-asserted-by":"publisher","first-page":"1102","DOI":"10.1002\/prot.24244","volume":"81","author":"J Minning","year":"2013","unstructured":"Minning J, Porto M, Bastolla U (2013) Detecting selection for negative design in proteins through an improved model of the misfolded state. Proteins Struct Func Bioinfo 81:1102\u20131112","journal-title":"Proteins Struct Func Bioinfo"},{"key":"10348_CR33","doi-asserted-by":"crossref","unstructured":"Morcos F, Pagnani A, Lunt B, Bertolino A, Marks DS, Sander C, Weigt M (2011) Direct-coupling analysis of residue coevolution captures native contacts across many protein families. Proc Nat Acad Sci 108:E1293\u2013E1301","DOI":"10.1073\/pnas.1111471108"},{"key":"10348_CR34","doi-asserted-by":"publisher","first-page":"4248","DOI":"10.1073\/pnas.0907953107","volume":"107","author":"V Mustonen","year":"2010","unstructured":"Mustonen V, L\u00e4ssig M (2010) Fitness flux and ubiquity of adaptive evolution. Proc Nat Acad Sci 107:4248\u20134253","journal-title":"Proc Nat Acad Sci"},{"issue":"1","key":"10348_CR35","doi-asserted-by":"publisher","first-page":"1046","DOI":"10.1038\/s42003-024-06616-5","volume":"7","author":"F Puente-S\u00e1nchez","year":"2024","unstructured":"Puente-S\u00e1nchez F, Pascual-Garc\u00eda A, Bastolla U, Pedr\u00f3s-Ali\u00f3 C, Tamames J (2024) Cross-biome microbial networks reveal functional redundancy and suggest genome reduction through functional complementarity. Commun Biol 7:1046","journal-title":"Commun Biol"},{"key":"10348_CR36","doi-asserted-by":"publisher","first-page":"9541","DOI":"10.1073\/pnas.0501865102","volume":"102","author":"G Sella","year":"2005","unstructured":"Sella G, Hirsh AE (2005) The application of statistical physics to evolutionary biology. Proc Nat Acad Sci 102:9541\u20139546","journal-title":"Proc Nat Acad Sci"},{"key":"10348_CR37","doi-asserted-by":"publisher","first-page":"229","DOI":"10.1016\/0959-440X(95)80081-6","volume":"5","author":"MJ Sippl","year":"1995","unstructured":"Sippl MJ (1995) Knowledge-based potentials for proteins. Curr op struct biol 5:229\u2013235","journal-title":"Curr Op Struct Biol"},{"key":"10348_CR38","doi-asserted-by":"publisher","first-page":"1815","DOI":"10.1093\/bioinformatics\/btt259","volume":"29","author":"MJ Skwark","year":"2013","unstructured":"Skwark MJ, Abdel-Rehim A, Elofsson A (2013) PconsC: combination of direct information methods and alignments improves contact prediction. Bioinformatics 29:1815\u20131816","journal-title":"Bioinformatics"},{"key":"10348_CR40","doi-asserted-by":"publisher","first-page":"67","DOI":"10.1073\/pnas.0805923106","volume":"106","author":"M Weigt","year":"2009","unstructured":"Weigt M, White RA, Szurmant H, Hoch JA, Hwa T (2009) Identification of direct residue contacts in protein-protein interaction by message passing. Proc Nat Acad Sci 106:67\u201372","journal-title":"Proc Nat Acad Sci"},{"key":"10348_CR41","doi-asserted-by":"publisher","first-page":"235","DOI":"10.1103\/RevModPhys.54.235","volume":"54","author":"FY Wu","year":"1982","unstructured":"Wu FY (1982) The Potts model. Rev Mod Phys 54:235\u2013268","journal-title":"Rev Mod Phys"},{"key":"10348_CR42","doi-asserted-by":"publisher","DOI":"10.1093\/acprof:oso\/9780198567028.001.0001","volume-title":"Computational molecular evolution","author":"Z Yang","year":"2006","unstructured":"Yang Z (2006) Computational molecular evolution. OUP"},{"key":"10348_CR43","doi-asserted-by":"publisher","first-page":"16152","DOI":"10.1073\/pnas.0705366104","volume":"104","author":"KB Zeldovich","year":"2007","unstructured":"Zeldovich KB, Chen P, Shakhnovich EI (2007) Protein stability imposes limits on organism complexity and speed of molecular evolution. Proc Nat Acad Sc 104:16152\u201316157","journal-title":"Proc Nat Acad Sc"}],"container-title":["Journal of Molecular Evolution"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/link.springer.com\/content\/pdf\/10.1007\/s00239-026-10348-x.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"text-mining"},{"URL":"https:\/\/link.springer.com\/article\/10.1007\/s00239-026-10348-x","content-type":"text\/html","content-version":"vor","intended-application":"text-mining"},{"URL":"https:\/\/link.springer.com\/content\/pdf\/10.1007\/s00239-026-10348-x.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2026,9,23]],"date-time":"2026-09-23T14:30:10Z","timestamp":1790173810000},"score":1,"resource":{"primary":{"URL":"https:\/\/link.springer.com\/10.1007\/s00239-026-10348-x"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2026,9,23]]},"references-count":39,"alternative-id":["10348"],"URL":"https:\/\/doi.org\/10.1007\/s00239-026-10348-x","relation":{},"ISSN":["0022-2844","1432-1432"],"issn-type":[{"value":"0022-2844","type":"print"},{"value":"1432-1432","type":"electronic"}],"subject":[],"published":{"date-parts":[[2026,9,23]]},"assertion":[{"value":"21 April 2026","order":1,"name":"received","label":"Received","group":{"name":"ArticleHistory","label":"Article History"}},{"value":"1 September 2026","order":2,"name":"accepted","label":"Accepted","group":{"name":"ArticleHistory","label":"Article History"}},{"value":"23 September 2026","order":3,"name":"first_online","label":"First Online","group":{"name":"ArticleHistory","label":"Article History"}},{"value":"The authors declare no competing interests.","order":1,"name":"Ethics","label":"Conflict of interest","group":{"name":"EthicsHeading","label":"Declarations"}}]}}