Understanding how the catalytic mechanisms of enzymes are optimized through evolution remains a major challenge in molecular biology. The concept of co-evolution implicates that compensatory mutations occur to preserve the structure and function of proteins. We have combined statistical analysis of protein sequences with the sensitivity of single molecule forceclamp spectroscopy to probe how catalysis is affected by structurally distant correlated mutations in Escherichia coli thioredoxin. Our findings show that evolutionary anti-correlated mutations have an inhibitory effect on enzyme catalysis, whereas positively correlated mutations rescue the catalytic activity. We interpret these results in terms of an evolutionary tuning of both the enzyme-substrate binding process and the chemistry of the active site. Our results constitute a direct observation of distant residue co-evolution in enzyme catalysis.The acquisition of adequate activity by an enzyme is the result of natural selection through mutagenesis (1). Mutations that represent a functional advantage may prevail, reflecting not only the optimization of function but also the emergence of a new enzymatic mechanism (1-4). Disadvantageous mutations can be compensated by mutations in other positions (5-8), as proposed by the concept of co-evolution. This hypothesis suggests that pairs of positions in the sequence evolve together in such a way that mutations in one of the residues induce a coherent response in the other to compensate for a certain functional or structural effect (5, 7). Over the last decades, numerous statistical methods have been applied to detect residue co-evolution in protein sequences (9, 10). Analysis of co-evolving residues has been used to explore functional coupling in processes like protein-protein interaction (11), gating of potassium channels (12), and allosteric communications in proteins (13). Experimental validation has been provided for the identification of networks of energetically coupled residues mediating allostery (14 -16). However, in the case of enzyme catalysis, only theoretical studies have been reported (17). The limitation of the available techniques in providing insights into the mechanism of catalysis has made it difficult to experimentally validate co-evolution as a natural means of catalysis optimization. Here, we have combined statistical sequence analysis and single molecule force-clamp spectroscopy to probe the effect of evolutionary correlated mutations in the Escherichia coli Trx catalytic mechanism, a ubiquitous disulfide bond reductase.In a previous study, we developed single molecule forceclamp spectroscopy as a novel sensitive tool to study the mechanism of the catalytic reaction of E. coli Trx (18). We showed that Trx reduces disulfide bonds using two distinct force-dependent catalytic mechanisms. The first mechanism is favored at low forces (up to 200 pN), 3 whereas the second mechanism is predominant at high forces (over 200 pN). These two chemical reduction mechanisms are likely to have evolved to al...