Thermostable and promiscuous Precambrian proteins

Click, Compute, Create: A Review of Web-based Tools for Enzyme Engineering

Enzyme engineering, though pivotal across various biotechnological domains, is often plagued by its time-consuming and labor-intensive nature. This review aims to offer an overview of supportive in silico methodologies for this demanding endeavor. Starting from methods to predict protein structures, to classification of their activity and even the discovery of new enzymes we continue with describing tools used to increase thermostability and production yields of selected targets. Subsequently, we discuss computational methods to modulate both, the activity as well as selectivity of enzymes. Last, we present recent approaches based on cutting-edge machine learning methods to redesign enzymes. With exception of the last chapter, there is a strong focus on methods easily accessible via web-interfaces or simple Python-scripts, therefore readily useable for a diverse and broad community.

Evolution of CRISPR-associated endonucleases as inferred from resurrected proteins

Clustered regularly interspaced short palindromic repeats (CRISPR)-associated Cas9 is an effector protein that targets invading DNA and plays a major role in the prokaryotic adaptive immune system. Although Streptococcus pyogenes CRISPR-Cas9 has been widely studied and repurposed for applications including genome editing, its origin and evolution are poorly understood. Here, we investigate the evolution of Cas9 from resurrected ancient nucleases (anCas) in extinct firmicutes species that last lived 2.6 billion years before the present. We demonstrate that these ancient forms were much more flexible in their guide RNA and protospacer-adjacent motif requirements compared with modern-day Cas9 enzymes. Furthermore, anCas portrays a gradual palaeoenzymatic adaptation from nickase to double-strand break activity, exhibits high levels of activity with both single-stranded DNA and single-stranded RNA targets and is capable of editing activity in human cells. Prediction and characterization of anCas with a resurrected protein approach uncovers an evolutionary trajectory leading to functionally flexible ancient enzymes.

Ancestral sequence reconstruction as a tool to study the evolution of wood decaying fungi

The study of evolution is limited by the techniques available to do so. Aside from the use of the fossil record, molecular phylogenetics can provide a detailed characterization of evolutionary histories using genes, genomes and proteins. However, these tools provide scarce biochemical information of the organisms and systems of interest and are therefore very limited when they come to explain protein evolution. In the past decade, this limitation has been overcome by the development of ancestral sequence reconstruction (ASR) methods. ASR allows the subsequent resurrection in the laboratory of inferred proteins from now extinct organisms, becoming an outstanding tool to study enzyme evolution. Here we review the recent advances in ASR methods and their application to study fungal evolution, with special focus on wood-decay fungi as essential organisms in the global carbon cycling.

Ancestral Sequence Reconstruction and Alternate Amino Acid States Guide Protein Library Design for Directed Evolution

Engineered proteins possess nearly limitless possibilities in medical and industrial applications but finding a precise amino acid sequence for these applications is challenging. A robust approach for discovering protein sequences with a desired functionality uses a library design method in which combinations of mutations are applied to a robust starting point. Determining useful mutations can be tortuous, yet rewarding; in this chapter, we present a novel library design method that uses information provided by ancestral sequence reconstruction (ASR) to create a library likely to have stable proteins with diverse function. ASR computational tools use a multi-sequence alignment of homologous proteins and an evolutionary model to estimate the protein sequences of the numerous common ancestors. For all ancestors, these tools calculate the probability of every amino acid occurring at each position within the sequence alignment. The alternate amino acid states at individual positions corelate to a region of stability in sequence space around the ancestral sequence which can inform site-wise diversification within a combinatorial library. The method presented in this chapter balances the quality of results, the computational resources needed, and ease of use.

Ancestral Sequence Reconstruction of Ancient Lipase from Family I.3 Bacterial Lipolytic Enzymes

Family I.3 lipase is distinguished from other families by the amino acid sequence and secretion mechanism. Little is known about the evolutionary process driving these differences. This study attempt to understand how the diverse temperature stabilities of bacterial lipases from family I.3 evolved. To achieve that, eighty-three protein sequences sharing a minimum 30% sequence identity with Antarctic Pseudomonas sp. AMS8 lipase were used to infer phylogenetic tree. Using ancestral sequence reconstruction (ASR) technique, the last universal common ancestor (LUCA) sequence of family I.3 was reconstructed. A gene encoding LUCA was synthesized, cloned and expressed as inclusion bodies in E. coli system. Insoluble form of LUCA was refolded using urea dilution method and then purified using affinity chromatography. The purified LUCA exhibited an optimum temperature and pH at 70℃ and 10 respectively. Various metal ions increased or retained the activity of LUCA. LUCA also demonstrated tolerance towards various organic solvents in 25% v/v concentration. The finding from this study could support the understanding on temperature and environment during ancient time. In overall, reconstructed ancestral enzymes have improved physicochemical properties that make them suitable for industrial applications and ASR technique can be employed as a general technique for enzyme engineering.

Were Ancestral Proteins Less Specific?

Some have hypothesized that ancestral proteins were, on average, less specific than their descendants. If true, this would provide a universal axis along which to organize protein evolution and suggests that reconstructed ancestral proteins may be uniquely powerful tools for protein engineering. Ancestral sequence reconstruction studies are one line of evidence used to support this hypothesis. Previously, we performed such a study, investigating the evolution of peptide-binding specificity for the paralogs S100A5 and S100A6. The modern proteins appeared more specific than their last common ancestor (ancA5/A6), as each paralog bound a subset of the peptides bound by ancA5/A6. In this study, we revisit this transition, using quantitative phage display to measure the interactions of 30,533 random peptides with human S100A5, S100A6, and ancA5/A6. This unbiased screen reveals a different picture. While S100A5 and S100A6 do indeed bind to a subset of the peptides recognized by ancA5/A6, they also acquired new peptide partners outside of the set recognized by ancA5/A6. Our previous work showed that ancA5/A6 had lower specificity than its descendants when measured against biological targets; our new work shows that ancA5/A6 has similar specificity to the modern proteins when measured against a random set of peptide targets. This demonstrates that altered biological specificity does not necessarily indicate altered intrinsic specificity, and sounds a cautionary note for using ancestral reconstruction studies with biological targets as a means to infer global evolutionary trends in specificity.

Divergent Evolution of a Protein-Protein Interaction Revealed through Ancestral Sequence Reconstruction and Resurrection

The postsynaptic density extends across the postsynaptic dendritic spine with discs large (DLG) as the most abundant scaffolding protein. DLG dynamically alters the structure of the postsynaptic density, thus controlling the function and distribution of specific receptors at the synapse. DLG contains three PDZ domains and one important interaction governing postsynaptic architecture is that between the PDZ3 domain from DLG and a protein called cysteine-rich interactor of PDZ3 (CRIPT). However, little is known regarding functional evolution of the PDZ3:CRIPT interaction. Here, we subjected PDZ3 and CRIPT to ancestral sequence reconstruction, resurrection, and biophysical experiments. We show that the PDZ3:CRIPT interaction is an ancient interaction, which was likely present in the last common ancestor of Eukaryotes, and that high affinity is maintained in most extant animal phyla. However, affinity is low in nematodes and insects, raising questions about the physiological function of the interaction in species from these animal groups. Our findings demonstrate how an apparently established protein–protein interaction involved in cellular scaffolding in bilaterians can suddenly be subject to dynamic evolution including possible loss of function.

Mitochondrial copper and phosphate transporter specificity was defined early in the evolution of eukaryotes

The mitochondrial carrier family protein SLC25A3 transports both copper and phosphate in mammals, yet in Saccharomyces cerevisiae the transport of these substrates is partitioned across two paralogs: PIC2 and MIR1. To understand the ancestral state of copper and phosphate transport in mitochondria, we explored the evolutionary relationships of PIC2 and MIR1 orthologs across the eukaryotic tree of life. Phylogenetic analyses revealed that PIC2-like and MIR1-like orthologs are present in all major eukaryotic supergroups, indicating an ancient gene duplication created these paralogs. To link this phylogenetic signal to protein function, we used structural modeling and site-directed mutagenesis to identify residues involved in copper and phosphate transport. Based on these analyses, we generated an L175A variant of mouse SLC25A3 that retains the ability to transport copper but not phosphate. This work highlights the utility of using an evolutionary framework to uncover amino acids involved in substrate recognition by mitochondrial carrier family proteins.

FireProtASR: A Web Server for Fully Automated Ancestral Sequence Reconstruction

There is a great interest in increasing proteins’ stability to widen their usability in numerous biomedical and biotechnological applications. However, native proteins cannot usually withstand the harsh industrial environment, since they are evolved to function under mild conditions. Ancestral sequence reconstruction is a well-established method for deducing the evolutionary history of genes. Besides its applicability to discover the most probable evolutionary ancestors of the modern proteins, ancestral sequence reconstruction has proven to be a useful approach for the design of highly stable proteins. Recently, several computational tools were developed, which make the ancestral reconstruction algorithms accessible to the community, while leaving the most crucial steps of the preparation of the input data on users’ side. FireProt^ASR aims to overcome this obstacle by constructing a fully automated workflow, allowing even the unexperienced users to obtain ancestral sequences based on a sequence query as the only input. FireProt^ASR is complemented with an interactive, easy-to-use web interface and is freely available at https://loschmidt.chemi.muni.cz/fireprotasr/.

Resurrection of ancestral effector caspases identifies novel networks for evolution of substrate specificity

Apoptotic caspases evolved with metazoans more than 950 million years ago (MYA), and a series of gene duplications resulted in two subfamilies consisting of initiator and effector caspases. The effector caspase genes (caspases-3, -6, and -7) were subsequently fixed into the Chordata phylum more than 650 MYA when the gene for a common ancestor (CA) duplicated, and the three effector caspases have persisted throughout mammalian evolution. All caspases prefer an aspartate residue at the P1 position of substrates, so each caspase evolved discrete cellular roles through changes in substrate recognition at the P4 position combined with allosteric regulation. We examined the evolution of substrate specificity in caspase-6, which prefers valine at the P4 residue, compared with caspases-3 and -7, which prefer aspartate, by reconstructing the CA of effector caspases (AncCP-Ef1) and the CA of caspase-6 (AncCP-6An). We show that AncCP-Ef1 is a promiscuous enzyme with little distinction between Asp, Val, or Leu at P4. The specificity of caspase-6 was defined early in its evolution, where AncCP-6An demonstrates a preference for Val over Asp at P4. Structures of AncCP-Ef1 and of AncCP-6An show a network of charged amino acids near the S4 pocket that, when combined with repositioning a flexible active site loop, resulted in a more hydrophobic binding pocket in AncCP-6An. The ancestral protein reconstructions show that the caspase-hemoglobinase fold has been conserved for over 650 million years and that only three substitutions in the scaffold are necessary to shift substrate selection toward Val over Asp.

Ancestral Resurrection and Directed Evolution of Fungal Mesozoic Laccases

Ancestral sequence reconstruction and resurrection provides useful information for protein engineering, yet its alliance with directed evolution has been little explored. In this study, we have resurrected several ancestral nodes of fungal laccases dating back ∼500 to 250 million years. Unlike modern laccases, the resurrected Mesozoic laccases were readily secreted by yeast, with similar kinetic parameters, a broader stability, and distinct pH activity profiles. The resurrected Agaricomycetes laccase carried 136 ancestral mutations, a molecular testimony to its origin, and it was subjected to directed evolution in order to improve the rate of 1,3-cyclopentanedione oxidation, a β-diketone initiator commonly used in vinyl polymerization reactions.IMPORTANCE The broad variety of biotechnological uses of fungal laccases is beyond doubt (food, textiles, pulp and paper, pharma, biofuels, cosmetics, and bioremediation), and protein engineering (in particular, directed evolution) has become the key driver for adaptation of these enzymes to harsh industrial conditions. Usually, the first requirement for directed laccase evolution is heterologous expression, which presents an important hurdle and often a time-consuming process. In this work, we resurrected a fungal Mesozoic laccase node which showed strikingly high heterologous expression and pH stability. As a proof of concept that the ancestral laccase is a suitable blueprint for engineering, we performed a quick directed evolution campaign geared to the oxidation of the β-diketone 1,3-cyclopentanedione, a poor laccase substrate that is used in the polymerization of vinyl monomers.

A surrogate structural platform informed by ancestral reconstruction and resurrection of a putative carbohydrate binding module hybrid illuminates the neofunctionalization of a pectate lyase

Yersinia enterocolitica is a pectinolytic zoonotic foodborne pathogen, the genome of which contains pectin-binding proteins and several different classes of pectinases, including polysaccharide lyases (PLs) and an exopolygalacturonase. These proteins operate within a coordinated pathway to completely saccharify homogalacturonan (HG). Polysaccharide lyase family 2 (PL2) is divided into two major subfamilies that are broadly-associated with contrasting 'endolytic' (PL2A) or 'exolytic' (PL2B) activities on HG. In the Y. enterocolitica genome, the PL2A gene is adjacent to an independent carbohydrate binding module from family 32 (YeCBM32), which possesses a N-terminal secretion tag and is known to specifically bind HG. Independent CBMs are rare in nature and, most commonly, are fused to enzymes in order to potentiate catalysis. The unconventional gene architecture of YePL2A and YeCBM32, therefore, may represent an ancestral relic of a fission event that decoupled PL2A from its cognate CBM. To provide further insight into the evolution of this pectinolytic locus and the molecular basis of HG depolymerisation within Y. enterocolitica, we have resurrected a YePL2A-YeCBM32 chimera and demonstrated that the extant PL2A digests HG more efficiently. In addition, we have engineered a tryptophan from the active site of the exolytic YePL2B into YePL2A (YePL2A-K291W) and demonstrated, using X-ray crystallography of substrate complexes, that it is a structural determinant of exo-activity within the PL2 family. In this manner, surrogate structural platforms may assist in the study of phylogenetic relationships informed by extant and resurrected sequences, and can be used to overcome challenging structural problems within carbohydrate active enzyme families.

Conservation of Specificity in Two Low-Specificity Proteins

Many regulatory proteins bind peptide regions of target proteins and modulate their activity. Such regulatory proteins can often interact with highly diverse target peptides. In many instances, it is not known if the peptide-binding interface discriminates targets in a biological context, or whether biological specificity is achieved exclusively through external factors such as subcellular localization. We used an evolutionary biochemical approach to distinguish these possibilities for two such low-specificity proteins: S100A5 and S100A6. We used isothermal titration calorimetry to study the binding of peptides with diverse sequence and biochemistry to human S100A5 and S100A6. These proteins bound distinct, but overlapping, sets of peptide targets. We then studied the peptide binding properties of orthologs sampled from across five amniote species. Binding specificity was conserved along all lineages, for the last 320 million years, despite the low specificity of each protein. We used ancestral sequence reconstruction to determine the binding specificity of the last common ancestor of the paralogs. The ancestor bound the entire set of peptides bound by modern S100A5 and S100A6 proteins, suggesting that paralog specificity evolved via subfunctionalization. To rule out the possibility that specificity is conserved because it is difficult to modify, we identified a single historical mutation that, when reverted in human S100A5, gave it the ability to bind an S100A6-specific peptide. These results reveal strong evolutionary constraints on peptide binding specificity. Despite being able to bind a large number of targets, the specificity of S100 peptide interfaces is likely important for the biology of these proteins.

Computational tools help improve protein stability but with a solubility tradeoff

Accurately predicting changes in protein stability upon amino acid substitution is a much sought after goal. Destabilizing mutations are often implicated in disease, whereas stabilizing mutations are of great value for industrial and therapeutic biotechnology. Increasing protein stability is an especially challenging task, with random substitution yielding stabilizing mutations in only ∼2% of cases. To overcome this bottleneck, computational tools that aim to predict the effect of mutations have been developed; however, achieving accuracy and consistency remains challenging. Here, we combined 11 freely available tools into a meta-predictor (meieringlab.uwaterloo.ca/stabilitypredict/). Validation against ∼600 experimental mutations indicated that our meta-predictor has improved performance over any of the individual tools. The meta-predictor was then used to recommend 10 mutations in a previously designed protein of moderate thermodynamic stability, ThreeFoil. Experimental characterization showed that four mutations increased protein stability and could be amplified through ThreeFoil's structural symmetry to yield several multiple mutants with >2-kcal/mol stabilization. By avoiding residues within functional ties, we could maintain ThreeFoil's glycan-binding capacity. Despite successfully achieving substantial stabilization, however, almost all mutations decreased protein solubility, the most common cause of protein design failure. Examination of the 600-mutation data set revealed that stabilizing mutations on the protein surface tend to increase hydrophobicity and that the individual tools favor this approach to gain stability. Thus, whereas currently available tools can increase protein stability and combining them into a meta-predictor yields enhanced reliability, improvements to the potentials/force fields underlying these tools are needed to avoid gaining protein stability at the cost of solubility.

Combining ancestral sequence reconstruction with protein design to identify an interface hotspot in a key metabolic enzyme complex

It is important to identify hotspot residues that determine protein-protein interactions in interfaces of macromolecular complexes. We have applied a combination of ancestral sequence reconstruction and protein design to identify hotspots within imidazole glycerol phosphate synthase (ImGPS). ImGPS is a key metabolic enzyme complex, which links histidine and de novo purine biosynthesis and consists of the cyclase subunit HisF and the glutaminase subunit HisH. Initial fluorescence titration experiments showed that HisH from Zymomonas mobilis (zmHisH) binds with high affinity to the reconstructed HisF from the last universal common ancestor (LUCA-HisF) but not to HisF from Pyrobaculum arsenaticum (paHisF), which differ by 103 residues. Subsequent titration experiments with a reconstructed evolutionary intermediate linking LUCA-HisF and paHisF and inspection of the subunit interface of a contemporary ImGPS allowed us to narrow down the differences crucial for zmHisH binding to nine amino acids of HisF. Homology modeling and in silico mutagenesis studies suggested that at most two of these nine HisF residues are crucial for zmHisH binding. These computational results were verified by experimental site-directed mutagenesis, which finally enabled us to pinpoint a single amino acid residue in HisF that is decisive for high-affinity binding of zmHisH. Our work shows that the identification of protein interface hotspots can be very efficient when reconstructed proteins with different binding properties are included in the analysis. Proteins 2017; 85:312-321. © 2016 Wiley Periodicals, Inc.

Experimental recreation of the evolution of lignin-degrading enzymes from the Jurassic to date

BACKGROUND

Floudas et al. (Science 336: 1715) established that lignin-degrading fungi appeared at the end of Carboniferous period associated with the production of the first ligninolytic peroxidases. Here, the subsequent evolution of these enzymes in Polyporales, where most wood-rotting fungi are included, is experimentally recreated using genomic information.

RESULTS

With this purpose, we analyzed the evolutionary pathway leading to the most efficient lignin-degrading peroxidases characterizing Polyporales species. After sequence reconstruction from 113 genes of ten sequenced genomes, the main enzyme intermediates were resurrected and characterized. Biochemical changes were analyzed together with predicted sequences and structures, to understand how these enzymes acquired the ability to degrade lignin and how this ability changed with time. The most probable first peroxidase in Polyporales would be a manganese peroxidase (Mn³⁺ oxidizing phenolic lignin) that did not change substantially until the appearance of an exposed tryptophan (oxidizing nonphenolic lignin) originating an ancestral versatile peroxidase. Later, a quick evolution, with loss of the Mn²⁺-binding site, generated the first lignin peroxidase that evolved to the extant form by improving the catalytic efficiency. Increased stability at acidic pH, which strongly increases the oxidizing power of these enzymes, was observed paralleling the appearance of the exposed catalytic tryptophan.

CONCLUSIONS

We show how the change in peroxidase catalytic activities meant an evolutionary exploration for more efficient ways of lignin degradation by fungi, a key step for carbon recycling in land ecosystems. The study provides ancestral enzymes with a potential biotechnological interest for the sustainable production of fuels and chemicals in a biomass-based economy.

Ancestral Protein Reconstruction and Circular Permutation for Improving the Stability and Dynamic Range of FRET Sensors

Small molecule biosensors based on Förster resonance energy transfer (FRET) enable small molecule signaling to be monitored with high spatial and temporal resolution in complex cellular environments. FRET sensors can be constructed by fusing a pair of fluorescent proteins to a suitable recognition domain, such as a member of the solute-binding protein (SBP) superfamily. However, naturally occurring SBPs may be unsuitable for incorporation into FRET sensors due to their low thermostability, which may preclude imaging under physiological conditions, or because the positions of their N- and C-termini may be suboptimal for fusion of fluorescent proteins, which may limit the dynamic range of the resulting sensors. Here, we show how these problems can be overcome using ancestral protein reconstruction and circular permutation. Ancestral protein reconstruction, used as a protein engineering strategy, leverages phylogenetic information to improve the thermostability of proteins, while circular permutation enables the termini of an SBP to be repositioned to maximize the dynamic range of the resulting FRET sensor. We also provide a protocol for cloning the engineered SBPs into FRET sensor constructs using Golden Gate assembly and discuss considerations for in situ characterization of the FRET sensors.

Ancestral protein reconstruction: techniques and applications

Ancestral sequence reconstruction (ASR) is the calculation of ancient protein sequences on the basis of extant ones. It is most powerful in combination with the experimental characterization of the corresponding proteins. Such analyses allow for the study of problems that are otherwise intractable. For example, ASR has been used to characterize ancestral enzymes dating back to the Paleoarchean era and to deduce properties of the corresponding habitats. In addition, the historical approach underlying ASR enables the identification of amino acid residues key to protein function, which is often not possible by only comparing extant proteins. Along these lines, residues responsible for the spectroscopic properties of protein pigments were identified as well as residues determining the binding specificity of steroid receptors. Further applications are studies related to the longevity of mutations, the contribution of gene duplications to enzyme functionalization, and the evolution of protein complexes. For these applications of ASR, we discuss recent examples; moreover, we introduce the basic principles of the underlying algorithms and present state-of-the-art protocols.

The thermostability and specificity of ancient proteins

Were ancient proteins systematically different than modern proteins? The answer to this question is profoundly important, shaping how we understand the origins of protein biochemical, biophysical, and functional properties. Ancestral sequence reconstruction (ASR), a phylogenetic approach to infer the sequences of ancestral proteins, may reveal such trends. We discuss two proposed trends: a transition from higher to lower thermostability and a tendency for proteins to acquire higher specificity over time. We review the evidence for elevated ancestral thermostability and discuss its possible origins in a changing environmental temperature and/or reconstruction bias. We also conclude that there is, as yet, insufficient data to support a trend from promiscuity to specificity. Finally, we propose future work to understand these proposed evolutionary trends.

Mutational studies on resurrected ancestral proteins reveal conservation of site-specific amino acid preferences throughout evolutionary history

Local protein interactions (“molecular context” effects) dictate amino acid replacements and can be described in terms of site-specific, energetic preferences for any different amino acid. It has been recently debated whether these preferences remain approximately constant during evolution or whether, due to coevolution of sites, they change strongly. Such research highlights an unresolved and fundamental issue with far-reaching implications for phylogenetic analysis and molecular evolution modeling. Here, we take advantage of the recent availability of phenotypically supported laboratory resurrections of Precambrian thioredoxins and β-lactamases to experimentally address the change of site-specific amino acid preferences over long geological timescales. Extensive mutational analyses support the notion that evolutionary adjustment to a new amino acid may occur, but to a large extent this is insufficient to erase the primitive preference for amino acid replacements. Generally, site-specific amino acid preferences appear to remain conserved throughout evolutionary history despite local sequence divergence. We show such preference conservation to be readily understandable in molecular terms and we provide crystallographic evidence for an intriguing structural-switch mechanism: Energetic preference for an ancestral amino acid in a modern protein can be linked to reorganization upon mutation to the ancestral local structure around the mutated site. Finally, we point out that site-specific preference conservation naturally leads to one plausible evolutionary explanation for the existence of intragenic global suppressor mutations.

Evolution of conformational dynamics determines the conversion of a promiscuous generalist into a specialist enzyme

β-Lactamases are produced by many modern bacteria as a mechanism of resistance toward β-lactam antibiotics, the most common antibiotics in use. β-Lactamases, however, are ancient enzymes that originated billions of years ago. Recently, proteins corresponding to 2- to 3-Gy-old Precambrian nodes in the evolution of Class A β-lactamases have been prepared and shown to be moderately efficient promiscuous catalysts, able to degrade a variety of antibiotics with catalytic efficiency levels similar to those of an average modern enzyme. Remarkably, there are few structural differences (in particular at the active-site regions) between the resurrected enzymes and a penicillin-specialist modern β-lactamase. Here, we propose that the ancestral promiscuity originates from conformational dynamics. We investigate the differences in conformational dynamics of the ancient and extant β-lactamases through MD simulations and quantify the contribution of each position to functionally related dynamics through Dynamic Flexibility Index. The modern TEM-1 lactamase shows a comparatively rigid active-site region, likely reflecting adaptation for efficient degradation of a specific substrate (penicillin), whereas enhanced deformability at the active-site neighborhood in the ancestral resurrected proteins likely accounts for the binding and subsequent degradation of antibiotic molecules of different size and shape. Clustering of the conformational dynamics on the basis of Principal Component Analysis is in agreement with the functional divergence, as the ancient β-lactamases cluster together, separated from their modern descendant. Finally, our analysis leads to testable predictions, as sites of potential relevance for the evolution of dynamics are identified and mutations at those sites are expected to alter substrate-specificity.