Protein engineering studies often suggest the emergence of completely new enzyme functionalities to be highly improbable. However, enzymes likely catalysed many different reactions already in the last universal common ancestor. Mechanisms for the emergence of completely new active sites must therefore either plausibly exist or at least have existed at the primordial protein stage. Here, we use resurrected Precambrian proteins as scaffolds for protein engineering and demonstrate that a new active site can be generated through a single hydrophobic-to-ionizable amino acid replacement that generates a partially buried group with perturbed physico-chemical properties. We provide experimental and computational evidence that conformational flexibility can assist the emergence and subsequent evolution of new active sites by improving substrate and transition-state binding, through the sampling of many potentially productive conformations. Our results suggest a mechanism for the emergence of primordial enzymes and highlight the potential of ancestral reconstruction as a tool for protein engineering.

SummaryZipcode binding protein 1 (ZBP1) is an oncofetal RNA-binding protein that mediates the transport and local translation of β-actin mRNA by the KH3-KH4 di-domain, which is essential for neuronal development. The high-resolution structures of KH3-KH4 with their respective target sequences show that KH4 recognizes a non-canonical GGA sequence via an enlarged and dynamic hydrophobic groove, whereas KH3 binding to a core CA sequence occurs with low specificity. A data-informed kinetic simulation of the two-step binding reaction reveals that the overall reaction is driven by the second binding event and that the moderate affinities of the individual interactions favor RNA looping. Furthermore, the concentration of ZBP1, but not of the target RNA, modulates the interaction, which explains the functional significance of enhanced ZBP1 expression during embryonic development.

In eukaryotes, RNA-binding proteins that contain multiple K homology (KH) domains play a key role in coordinating the different steps of RNA synthesis, metabolism and localization. Understanding how the different KH modules participate in the recognition of the RNA targets is necessary to dissect the way these proteins operate. We have designed a KH mutant with impaired RNA-binding capability for general use in exploring the role of individual KH domains in the combinatorial functional recognition of RNA targets. A double mutation in the hallmark GxxG loop (GxxG-to-GDDG) impairs nucleic acid binding without compromising the stability of the domain. We analysed the impact of the GDDG mutations in individual KH domains on the functional properties of KSRP as a prototype of multiple KH domain-containing proteins. We show how the GDDG mutant can be used to directly link biophysical information on the sequence specificity of the different KH domains of KSRP and their role in mRNA recognition and decay. This work defines a general molecular biology tool for the investigation of the function of individual KH domains in nucleic acid binding proteins.

The relative promiscuity of hub proteins such as postsynaptic density protein-95 (PSD-95) can be achieved by alternative splicing, allosteric regulation, and post-translational modifications, the latter of which is the most efficient method of accelerating cellular responses to environmental changes in vivo. Here, a mutational approach was used to determine the impact of phosphorylation and succinimidation post-translational modifications on the binding affinity of the postsynaptic density protein-95/discs large/zonula occludens-1 (PDZ3) domain of PSD-95. Molecular dynamics simulations revealed that the binding affinity of this domain is influenced by an interplay between salt-bridges linking the α3 helix, the β2–β3 loop and the positively charged Lys residues in its high-affinity hexapeptide ligand KKETAV. The α3 helix is an extra structural element that is not present in other PDZ domains, which links PDZ3 with the following SH3 domain in the PSD-95 protein. This regulatory mechanism was confirmed experimentally via thermodynamic and NMR chemical shift perturbation analyses, discarding intra-domain long-range effects. Taken together, the results presented here reveal the molecular basis of the regulatory role of the α3 extra-element and the effects of post-translational modifications of PDZ3 on its binding affinity, both energetically and dynamically.

Evolution involves not only adaptation, but also the degradation of superfluous features. Many examples of degradation at the morphological level are known (vestigial organs, for instance). However, the impact of degradation on molecular evolution has been rarely addressed. Thioredoxins serve as general oxidoreductases in all cells. Here, we report extensive mutational analyses on the folding of modern and resurrected ancestral bacterial thioredoxins. Contrary to claims from recent literature, in vitro folding rates in the thioredoxin family are not evolutionary conserved, but span at least a ∼100-fold range. Furthermore, modern thioredoxin folding is often substantially slower than ancestral thioredoxin folding. Unassisted folding, as probed in vitro, thus emerges as an ancestral vestigial feature that underwent degradation, plausibly upon the evolutionary emergence of efficient cellular folding-assistance. More generally, our results provide evidence that degradation of ancestral features shapes, not only morphological evolution, but also the evolution of individual proteins.