Search

Enzymes catalyze the chemical reactions of life. For nearly half of known enzymes, catalysis requires the binding of small molecules known as cofactors. Polypeptide-cofactor complexes likely formed at a primordial stage and became starting points for the evolution of many efficient enzymes. Yet, evolution has no foresight so the driver for the primordial complex formation is unknown. Here, we use a resurrected ancestral TIM-barrel protein to identify one potential driver. Heme binding at a flexible region of the ancestral structure yields a peroxidation catalyst with enhanced efficiency when compared to free heme. This enhancement, however, does not arise from protein-mediated promotion of catalysis. Rather, it reflects protection of bound heme from common degradation processes and a resulting longer life time and higher effective concentration for the catalyst. Protection of catalytic cofactors by polypeptides emerges as a general mechanism to enhance catalysis and may have plausibly benefited primordial polypeptide-cofactor associations.

The routine generation of enzymes with completely new active sites is one of the major unsolved problems in protein engineering. Advances in this field have been so far modest, perhaps due, at least in part, to the widespread use of modern natural proteins as scaffolds for de novo engineering. Most modern proteins are highly evolved and specialized, and, consequently, difficult to repurpose for completely new functionalities. Conceivably, resurrected ancestral proteins with the biophysical properties that promote evolvability, such as high stability and conformational diversity, could provide better scaffolds for de novo enzyme generation. Kemp elimination, a non-natural reaction that provides a simple model of proton abstraction from carbon, has been extensively used as a benchmark in de novo enzyme engineering. Here, we present an engineered ancestral β-lactamase with a new active site capable of efficiently catalysing the Kemp elimination. Our Kemp eliminase is the outcome of a minimalist design based on a single function-generating mutation followed by sharply-focused, low-throughput library screening. Yet, its catalytic parameters (kcat/KM=2·105 M−1s−1, kcat=635 s−1) compare favourably with the average modern natural enzyme and with the best proton-abstraction de novo Kemp eliminases reported in the literature. General implications of our results for de novo enzyme engineering are discussed.

W.D.V.H. acknowledges support from National Institutes of Health (Grant: R01GM112077). S.B.O. acknowledges support from the Gordon and Betty Moore Foundations and National Science Foundation (Awards: 1715591 and 1901709). J.M.S.R. acknowledges support from Spanish Ministry of Economy and Competitiveness/FEDER Funds (Grants BIO2015-66426-R and RTI2018-097142-B-100) and the Human Frontier Science Program (Grant RGP0041/2017). V.A.R. acknowledges support from FEDER/Junta de Andalucia-Consejeria de Economia y Conocimiento (Grant E.FQM.113.UGR18). We would like to thank the beamline staff of ID30B of the ESRF (European Synchrotron Radiation Facility, Grenoble, France) for their assistance during data collection and the ESRF for the provision of time through proposals MX-2064., TEM-1 β-lactamase degrades β-lactam antibiotics with a strong preference for penicillins. Sequence reconstruction studies indicate that it evolved from ancestral enzymes that degraded a variety of β-lactam antibiotics with moderate efficiency. This generalist to specialist conversion involved more than 100 mutational changes, but conserved fold and catalytic residues, suggesting a role for dynamics in enzyme evolution. Here, we develop a conformational dynamics computational approach to rationally mold a protein flexibility profile on the basis of a hinge-shift mechanism. By deliberately weighting and altering the conformational dynamics of a putative Precambrian β-lactamase, we engineer enzyme specificity that mimics the modern TEM-1 β-lactamase with only 21 amino acid replacements. Our conformational dynamics design thus re-enacts the evolutionary process and provides a rational allosteric approach for manipulating function while conserving the enzyme active site., United States Department of Health & Human Services National Institutes of Health (NIH) - USA R01GM112077, Gordon and Betty Moore Foundations, National Science Foundation (NSF) 1715591 1901709, Spanish Ministry of Economy and Competitiveness/FEDER Funds BIO2015-66426-R RTI2018-097142-B-100, Human Frontier Science Program RGP0041/2017, FEDER/Junta de Andalucia-Consejeria de Economia y Conocimiento E.FQM.113.UGR18

Many metabolites are generated in one step of a biochemical pathway and consumed in a subsequent step. Such metabolic intermediates are often reactive molecules which, if allowed to freely diffuse in the intracellular milieu, could lead to undesirable side reactions and even become toxic to the cell. Therefore, metabolic intermediates are often protected as protein-bound species and directly transferred between enzyme active sites in multi-function al enzymes, multi-enzyme complexes, and metabolons. Sequestration of reactive metabolic intermediates thus con tributes to metabolic efficiency. It is not known, however, whether this evolutionary adaptation can be relaxed in response to challenges to organismal survival. Here, we report evolutionary repair experiments on Escherichia coli cells in which an enzyme crucial for the biosynthesis of proline has been deleted. The deletion makes cells unable to grow in a culture medium lacking proline. Remarkably, however, cell growth is efficiently restored by many single mutations (12 at least) in the gene of glutamine synthetase. The mutations cause the leakage to the intracellular milieu of a highly reactive phosphorylated intermediate common to the biosynthetic pathways of glutamine and pro line. This intermediate is generally assumed to exist only as a protein-bound species. Nevertheless, its diffusion upon mutation-induced leakage enables a new route to proline biosynthesis. Our results support that leakage of seques tered metabolic intermediates can readily occur and contribute to organismal adaptation in some scenarios. Enhanced availability of reactive molecules may enable the generation of new biochemical pathways and the poten tial of mutation-induced leakage in metabolic engineering is noted, Human Frontier Science Program RGP0041/2017, Spanish Government RTI2018-097142-B-100, Ministry of Science and Innovation, Spain (MICINN) 80NSSC18K1277, European Commission, Junta de Andalucia, Regional Andalusian Government E-BIO-464-UGR-20 2020_DOC_00541

Protein folding/unfolding processes involve a large number of weak, non-covalent interactions and are more appropriately described in terms of the movement of a point representing protein conformation in a plot of internal free energy versus conformational degrees of freedom. While these energy landscapes have an astronomically large number of dimensions, it has been shown that many relevant aspects of protein folding can be understood in terms of their projections onto a few relevant coordinates. Remarkably, such low-dimensional free energy surfaces can be obtained from experimental DSC data using suitable analytical models. Here, we describe the experimental procedures to be used to obtain the high-quality DSC data that are required for free-energy surface analysis.

Manipulating Conformational Dynamics To Repurpose Ancient Proteins for Modern Catalytic Functions

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 evolutionarily 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. © 2019 The Author(s)., This research was supported by FEDER Funds, grant BIO2015-66426-R from the Spanish Ministry of Economy and Competitiveness ( J.M.S.-R.), grant RGP0041/2017 from the Human Frontier Science Program ( J.M.S.-R. and E.A.G.) and National Institutes of Health 1R01AR069137 (E.A.G.), Department of Defence MURI W911NF-16-1-0372 (E.A.G.)

This work was supported by Human Frontier Science Program Grant RGP0041/2017 (J.M.S.-R. and E.A.G.), National Science Foundation Award #2032315 (E.A.G.), National Institutes of Health Award #R01AR069137 (E.A.G.), Department of Defense MURI Award #W911NF-16-1-0372 (E.A.G.), Spanish Ministry of Science and Innovation/FEDER Funds Grants RTI-2018-097142-B-100 (J.M.S.-R.) and BIO2016-74875-P (J.A.G.) and the Science, Engineering and Research Board (SERB, India) Grant MTR/2019/000392 (A.N.N.). We are grateful to the European Synchrotron Radiation Facility (ESRF), Grenoble, France, for the provision of time and the staff at ID23-1 beamline for assistance during data collection. J.M.S.R. designed the research. G.G.-A. purified the modern/ancestral chimeras and the thioredoxin variants; she also performed and analysed the experiments aimed at determining their folding kinetics and biomolecular properties. V.A.R. performed experiments addressed at determining the efficiency of heterologous expression and provided essential input for the molecular interpretation of mutational effects on expression efficiency. E.A.G. provided essential input for the evolutionary interpretation of the data. J.A.G. determined the X-ray structure of the symbiont protein and provided essential input regarding its interpretation and implications. A.N.N. performed the computational simulations of the folding landscape for thioredoxins and provided essential input regarding their engineering implications. B.I.M. and J.M.S.-R. directed the project. J.M.S.-R. wrote the first draft of the manuscript to which V.A.R. J.A.G. A.N.N. and B.I.M. added crucial paragraphs and sections. All authors discussed the manuscript, suggested modifications and improvements, and contributed to the final version. The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., Obligate symbionts typically exhibit high evolutionary rates. Consequently, their proteins may differ considerably from their modern and ancestral homologs in terms of both sequence and properties, thus providing excellent models to study protein evolution. Also, obligate symbionts are challenging to culture in the lab and proteins from uncultured organisms must be produced in heterologous hosts using recombinant DNA technology. Obligate symbionts thus replicate a fundamental scenario of metagenomics studies aimed at the functional characterization and biotechnological exploitation of proteins from the bacteria in soil. Here, we use the thioredoxin from Candidatus Photodesmus katoptron, an uncultured symbiont of flashlight fish, to explore evolutionary and engineering aspects of protein folding in heterologous hosts. The symbiont protein is a standard thioredoxin in terms of 3D-structure, stability and redox activity. However, its folding outside the original host is severely impaired, as shown by a very slow refolding in vitro and an inefficient expression in E. coli that leads mostly to insoluble protein. By contrast, resurrected Precambrian thioredoxins express efficiently in E. coli, plausibly reflecting an ancient adaptation to unassisted folding. We have used a statistical-mechanical model of the folding landscape to guide back-to-ancestor engineering of the symbiont protein. Remarkably, we find that the efficiency of heterologous expression correlates with the in vitro (i.e., unassisted) folding rate and that the ancestral expression efficiency can be achieved with only 1–2 back-to-ancestor replacements. These results demonstrate a minimal-perturbation, sequence-engineering approach to rescue inefficient heterologous expression which may potentially be useful in metagenomics efforts targeting recent adaptations., National Science Foundation 2032315, National Institutes of Health 01AR069137, U.S. Department of Defense 911NF-16-1-0372, Human Frontier Science Program RGP0041/2017, European Synchrotron Radiation Facility, Science and Engineering Research Board MTR/2019/000392, Ministerio de Ciencia e Innovación RTI-2018-097142-B-100

[EN] Glycosidases are phylogenetically widely distributed enzymes that are crucial for the cleavage of glycosidic bonds. Here, we present the exceptional properties of a putative ancestor of bacterial and eukaryotic family-1 glycosidases. The ancestral protein shares the TIM-barrel fold with its modern descendants but displays large regions with greatly enhanced conformational flexibility. Yet, the barrel core remains comparatively rigid and the ancestral glycosidase activity is stable, with an optimum temperature within the experimental range for thermophilic family-1 glycosidases. None of the ∼5500 reported crystallographic structures of ∼1400 modern glycosidases show a bound porphyrin. Remarkably, the ancestral glycosidase binds heme tightly and stoichiometrically at a well-defined buried site. Heme binding rigidifies this TIM-barrel and allosterically enhances catalysis. Our work demonstrates the capability of ancestral protein reconstructions to reveal valuable but unexpected biomolecular features when sampling distant sequence space. The potential of the ancestral glycosidase as a scaffold for custom catalysis and biosensor engineering is discussed., his work was supported by Human Frontier Science Program Grant RGP0041 (J.M.S.-R., E.A.G., B.S., and S.C.L.K.), NIH grant R01AR069137 (E.A.G.), Department of Defense grant MURI W911NF-16-1-0372 (E.A.G.), the Swedish Research Council (2019-03499) (S.C.L.K.), the Knut and Alice Wallenberg Foundation (2018.0140 and 2019.0431) (S.C.L.K.), Spanish Ministry of Economy and Competitiveness/FEDER Funds Grants BIO2015-66426-R (J.M.S.-R.) RTI2018-097142-B-100 (J.M.S.-R.) and BIO2016-74875-P (J.A.G.). The simulations were enabled by resources provided by the Swedish National Infrastructure for Computing (SNIC) at UPPMAX partially funded by the Swedish Research Council through grant agreement no. 2016-07213. We acknowledge the Spanish Synchrotron Radiation Facility (ALBA, Barcelona) for the provision of synchrotron radiation facilities and the staff at XALOC beamline for their invaluable support. We are also grateful to Victoria Longobardo Polanco (Proteomic Unit, Institute of Parasitology and Biomedicine “López-Neyra”) for help with mass spectrometry experiments and data analyses and to Juan Román Luque Ortega (Molecular Interactions Facility, Centro de Investigaciones Biológicas Margarita Salas) for help with ultracentrifugation experiments and data analyses.

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., This study is based upon work funded by the Spanish Government projects BIO2013-43407-R-DEWRY and BIO2016-79106-R-Lignolution and by the Bio Based Industries Joint Undertaking under the European Union’s Horizon 2020 Research and Innovation programme (grant agreement no. 886567). B.J.G.-F. was supported by FPI national fellowship BES-2014-068887.

Directed evolution has revolutionized protein engineering. Still, enzyme optimization by random library screening remains sluggish, in large part due to futile probing of mutations that are catalytically neutral and/or impair stability and folding. FuncLib is a novel approach which uses phylogenetic analysis and Rosetta design to rank enzyme variants with multiple mutations, on the basis of predicted stability. Here, we use it to target the active site region of a minimalist-designed, de novo Kemp eliminase. The similarity between the Michaelis complex and transition state for the enzymatic reaction makes this system particularly challenging to optimize. Yet, experimental screening of a small number of active-site variants at the top of the predicted stability ranking leads to catalytic efficiencies and turnover numbers ( 2 104 M 1 s 1 and 102 s 1) for this anthropogenic reaction that compare favorably to those of modern natural enzymes. This result illustrates the promise of FuncLib as a powerful tool with which to speed up directed evolution, even on scaffolds that were not originally evolved for those functions, by guiding screening to regions of the sequence space that encode stable and catalytically diverse enzymes. Empirical valence bond calculations reproduce the experimental activation energies for the optimized eliminases to within 2 kcal mol 1 and indicate that the enhanced activity is linked to better geometric preorganization of the active site. This raises the possibility of further enhancing the stabilityguidance of FuncLib by computational predictions of catalytic activity, as a generalized approach for computational enzyme design, Knut and Alice Wallenberg Foundation (Wallenberg Academy Fellowship) 2018.0140, Human Frontier Science Program RGP0041/2017, FEDER Funds/Spanish Ministry of Science, Innovation and Universities BIO2015-66426-R RTI2018-097142-B-100, FEDER/Junta de Andalucia - Consejeria de Economia y Conocimiento E.FQM.113.UGR18, Swedish National Infrastructure for computing (SNAC) 2018/2-3 2019/2-1

Glycosidases are phylogenetically widely distributed enzymes that are crucial for the cleavage of glycosidic bonds. Here, we present the exceptional properties of a putative ancestor of bacterial and eukaryotic family-1 glycosidases. The ancestral protein shares the TIM-barrel fold with its modern descendants but displays large regions with greatly enhanced conformational flexibility. Yet, the barrel core remains comparatively rigid and the ancestral glycosidase activity is stable, with an optimum temperature within the experimental range for thermophilic family-1 glycosidases. None of the ~5500 reported crystallographic structures of ~1400 modern glycosidases show a bound porphyrin. Remarkably, the ancestral glycosidase binds heme tightly and stoichiometrically at a well-defined buried site. Heme binding rigidifies this TIM-barrel and allosterically enhances catalysis. Our work demonstrates the capability of ancestral protein reconstructions to reveal valuable but unexpected biomolecular features when sampling distant sequence space. The potential of the ancestral glycosidase as a scaffold for custom catalysis and biosensor engineering is discussed.

Approximations to the sequences of ancestral proteins can be derived from the sequences of their modern descendants. Proteins encoded by such reconstructed sequences can be prepared in the laboratory and subjected to experimental scrutiny. These 'resurrected' ancestral proteins often display remarkable properties, reflecting ancestral adaptations to intra-cellular and extra-cellular environments that differed from the environments hosting modern/extant proteins. Recent experimental and computational work has specifically discussed high stability, substrate and catalytic promiscuity, conformational flexibility/diversity and altered patterns of interaction with other sub-cellular components. In this review, we discuss these remarkable properties as well as recent attempts to explore their biotechnological and protein-engineering potential.

Rubisco is an ancient, catalytically conserved yet slow enzyme, which plays a central role in the biosphere’s carbon cycle. The design of Rubiscos to increase agricultural productivity has hitherto relied on the use of in vivo selection systems, precluding the exploration of biochemical traits that are not wired to cell survival. We present a directed -in vitro- evolution platform that extracts the enzyme from its biological context to provide a new avenue for Rubisco engineering. Precambrian and extant form II Rubiscos were subjected to an ensemble of directed evolution strategies aimed at improving thermostability. The most recent ancestor of proteobacteria -dating back 2.4 billion years- was uniquely tolerant to mutagenic loading. Adaptive evolution, focused evolution and genetic drift revealed a panel of thermostable mutants, some deviating from the characteristic trade-offs in CO2-fixing speed and specificity. Our findings provide a novel approach for identifying Rubisco variants with improved catalytic evolution potential., This work was supported by the REPSOL Research contracts Rubolution (RC020401120018), Rubolution 2.0 (RC 020401140042), the CSIC project PIE-201780E043 and the Australian Research Council grant CE140100015.

Many experimental analyses and proposed scenarios support that ancient life was thermophilic. In congruence with this hypothesis, proteins encoded by reconstructed sequences corresponding to ancient phylogenetic nodes often display very high stability. Here, we show that such ‘reconstructed ancestral hyperstability’ can be further engineered on the basis of a straightforward approach that uses exclusively information afforded by the ancestral reconstruction process itself. Since evolution does not imply continuous progression, screening of the mutations between two evolutionarily related resurrected ancestral proteins may identify mutations that further stabilize the most stable one. To explore this approach, we have used a resurrected thioredoxin corresponding to the last common ancestor of the cyanobacterial, Deinococcus and Thermus groups (LPBCA thioredoxin), which has a denaturation temperature of ∼123°C. This high value is within the top 0.1% of the denaturation temperatures in the ProTherm database and, therefore, achieving further stabilization appears a priori as a challenging task. Nevertheless, experimental comparison with a resurrected thioredoxin corresponding to the last common ancestor of bacteria (denaturation temperature of ∼115°C) immediately identifies three mutations that increase the denaturation temperature of LPBCA thioredoxin to ∼128°C. Comparison between evolutionarily related resurrected ancestral proteins thus emerges as a simple approach to expand the capability of ancestral reconstruction to search sequence space for extreme protein properties of biotechnological interest. The fact that ancestral sequences for many phylogenetic nodes can be derived from a single alignment of modern sequences should contribute to the general applicability of this approach.

Enzymes are flexible catalysts, and there has been substantial discussion about the extent to which this flexibility contributes to their catalytic efficiency. What has been significantly less discussed is the extent to which this flexibility contributes to their evolvability. Despite this, recent years have seen an increasing number of both experimental and computational studies that demonstrate that cooperativity and flexibility play significant roles in enzyme innovation. This review covers key developments in the field that emphasize the importance of enzyme dynamics not just to the evolution of new enzyme function(s), but also as a property that can be harnessed in the design of new artificial enzymes.

Tzul et al. (1) report different unfolding rates and similar folding rates for a number of thioredoxins. The authors interpret this result as evidence of the principle of minimal frustration. Their study includes several resurrected Precambrian thioredoxins that we have previously prepared and characterized (2⇓⇓–5). We agree that the principle of minimal frustration is essential to understand protein evolution. However, approximate folding-rate invariance is easily explained without invoking this principle. Thioredoxin kinetic stability relies on a transition state that is substantially unstructured (5, 6). Therefore, mutations that changed unfolding rates to tune kinetic stability during evolution likely had much less effect on folding rates, as implied by the well-known principles of ϕ-value analysis (7). Moreover, our experimental results are not consistent with folding-rate invariance being a general feature of thioredoxins. Fig. 1 shows folding–unfolding rates for the modern Escherichia coli thioredoxin and a resurrected Precambrian thioredoxin. The unfolding of the ancestral protein is ∼three orders-of-magnitude slower than the unfolding of the modern protein, indicating enhanced kinetic stability. However, in clear … [↵][1]2To whom correspondence should be addressed. Email: sanchezr{at}ugr.es. [1]: #xref-corresp-1-1

Proviral factors are host proteins hijacked by viruses for processes essential for virus propagation such as cellular entry and replication. Pathogens and their hosts co-evolve. It follows that replacing a proviral factor with a functional ancestral form of the same protein could prevent viral propagation without fatally compromising organismal fitness. Here, we provide proof of concept of this notion. Thioredoxins serve as general oxidoreductases in all known cells. We report that several laboratory resurrections of Precambrian thioredoxins display substantial levels of functionality within Escherichia coli. Unlike E. coli thioredoxin, however, these ancestral thioredoxins are not efficiently recruited by the bacteriophage T7 for its replisome and therefore prevent phage propagation in E. coli. These results suggest an approach to the engineering of virus resistance. Diseases caused by viruses may have a devastating effect in agriculture. We discuss how the suggested approach could be applied to the engineering of plant virus resistance., This work was supported by grants BIO2012-34937, CSD2009-00088, and BIO2015-66426-R (J.M.S.-R.) from MINECO/FEFER and grant P09-CVI-5073 (B.I.-M.) from the “Junta de Andalucía” and Feder Funds.

Consensus-sequence engineering has generated protein variants with enhanced stability, and sometimes, with modulated biological function. Consensus mutations are often interpreted as the introduction of ancestral amino acid residues. However, the precise relationship between consensus engineering and ancestral protein resurrection is not fully understood. Here, we report the properties of proteins encoded by consensus sequences derived from a multiple sequence alignment of extant, class A β-lactamases, as compared with the properties of ancient Precambrian β-lactamases resurrected in the laboratory. These comparisons considered primary sequence, secondary, and tertiary structure, as well as stability and catalysis against different antibiotics. Out of the three consensus variants generated, one could not be expressed and purified (likely due to misfolding and/or low stability) and only one displayed substantial stability having substrate promiscuity, although to a lower extent than ancient β-lactamases. These results: (i) highlight the phenotypic differences between consensus variants and laboratory resurrections of ancestral proteins; (ii) question interpretations of consensus proteins as phenotypic proxies of ancestral proteins; and (iii) support the notion that ancient proteins provide a robust approach toward the preparation of protein variants having large numbers of mutational changes while possessing unique biomolecular properties.

It is generally assumed that the amino acids that exist in all homologous enzymes correspond to residues that participate in catalysis, or that are essential for folding and stability. Although this holds for catalytic residues, the function of conserved noncatalytic residues is not clear. It is not known if such residues are of equal importance and have the same role in different homologous enzymes. In humans, the E104D mutation in triosephosphate isomerase (TIM) is the most frequent mutation in the autosomal diseases named "TPI deficiencies." We explored if the E104D mutation has the same impact in TIMs from four different organisms (Homo sapiens, Giardia lamblia, Trypanosoma cruzi, and T. brucei). The catalytic properties were not significantly affected by the mutation, but it affected the rate and extent of formation of active dimers from unfolded monomers differently. Scanning calorimetry experiments indicated that the mutation was in all cases destabilizing, but the mutation effect on rates of irreversible denaturation and transition-state energetics were drastically dependent on the TIM background. For instance, the E104D mutation produce changes in activation energy ranging from 430 kJ mol(-1) in HsTIM to -78 kJ mol(-1) in TcTIM. Thus, in TIM the role of a conserved noncatalytic residue is drastically dependent on its molecular background. Accordingly, it would seem that because each protein has a particular sequence, and a distinctive set of amino acid interactions, it should be regarded as a unique entity that has evolved for function and stability in the organisms to which it belongs.

Little is known with certainty about the evolution of protein structures in general and the degree of protein structure conservation over planetary time scales in particular. Here we report the X-ray crystal structures of seven laboratory resurrections of Precambrian thioredoxins dating back up to ~4 billion years before present. Despite considerable sequence differences compared with extant enzymes, the ancestral proteins display the canonical thioredoxin fold while only small structural changes have occurred over 4 billion years. This remarkable degree of structure conservation since a time near the last common ancestor of life supports a punctuated-equilibrium model of structure evolution in which the generation of new folds occurs over comparatively short periods of time and is followed by long periods of structural stasis.

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., This work was supported by Feder Funds, Grants from the Spanish Ministry of Economy and Competitiveness BIO2015-66426-R (J.M.S.-R.), CSD2009-00088 (J.M.S.-R.), CTQ2011-29299-C02-01 (F.S.-G.), CTQ2011-22514 (M.B.), BIO2016-74875-P (J.A.G.), ‘Factoría Española de Cristalización˜’, Consolider-Ingenio 2010 (J.A.G.) and CEI BioTic V19-2015 (V.A.R.), a Wallenberg Academy Fellowship (S.C.L.K.) and DuPont Young Professor Award (E.A.G.) and Grants NNX13AI08G and NNX13AI10G (E.A.G.) from NASA Exobiology. The European Research Council has provided financial support under the European Community’s Seventh Framework Programme (FP7/2007–2013)/ERC Grant Agreement No. 306474.

High stability and enhanced promiscuity (likely linked to conformational flexibility/diversity) contribute to evolvability and are advantageous features in protein scaffolds for laboratory-directed evolution and molecular design. Furthermore, the two features are not necessarily incompatible, and proteins may simultaneously be promiscuous/flexible and highly stable. In fact, it appears plausible that the combination of the two features was not uncommon among the most ancient proteins because (i) ancient life was likely thermophilic and (ii) ancient proteins were likely promiscuous generalists with broad functionalities. Phylogenetic analyses allow the reconstruction of ancestral sequences and provide an approach to explore the properties of ancient proteins. High stability and promiscuity have been often found for proteins encoded by reconstructed ancestral sequences, i.e., for “resurrected” ancestral proteins. The combination of the two features, i.e., the ancestral hyperstable generalist phenotype, has actually been obtained in recent studies. Ancestral protein resurrection thus emerges as a useful source of scaffolds for protein engineering.

The relationship between the denaturation temperatures of proteins (Tm values) and the living temperatures of their host organisms (environmental temperatures: TENV values) is poorly understood. Since different proteins in the same organism may show widely different Tm's, no simple universal relationship between Tm and TENV should hold, other than Tm≥TENV. Yet, when analyzing a set of homologous proteins from different hosts, Tm's are oftentimes found to correlate with TENV's but this correlation is shifted upward on the Tm axis. Supporting this trend, we recently reported Tm's for resurrected Precambrian thioredoxins that mirror a proposed environmental cooling over long geological time, while remaining a shocking ~50°C above the proposed ancestral ocean temperatures. Here, we show that natural selection for protein kinetic stability (denaturation rate) can produce a Tm↔TENV correlation with a large upward shift in Tm. A model for protein stability evolution suggests a link between the Tm shift and the in vivo lifetime of a protein and, more specifically, allows us to estimate ancestral environmental temperatures from experimental denaturation rates for resurrected Precambrian thioredoxins. The TENV values thus obtained match the proposed ancestral ocean cooling, support comparatively high Archaean temperatures, and are consistent with a recent proposal for the environmental temperature (above 75°C) that hosted the last universal common ancestor. More generally, this work provides a framework for understanding how features of protein stability reflect the environmental temperatures of the host organisms.

CBS (cystathionine β-synthase) is a multidomain tetrameric enzyme essential in the regulation of homocysteine metabolism, whose activity is enhanced by the allosteric regulator SAM (S-adenosylmethionine). Missense mutations in CBS are the major cause of inherited HCU (homocystinuria). In the present study we apply a novel approach based on a combination of calorimetric methods, functional assays and kinetic modelling to provide structural and energetic insight into the effects of SAM on the stability and activity of WT (wild-type) CBS and seven HCU-causing mutants. We found two sets of SAM-binding sites in the C-terminal regulatory domain with different structural and energetic features: a high affinity set of two sites, probably involved in kinetic stabilization of the regulatory domain, and a low affinity set of four sites, which are involved in the enzyme activation. We show that the regulatory domain displays a low kinetic stability in WT CBS, which is further decreased in many HCU-causing mutants. We propose that the SAM-induced stabilization may play a key role in modulating steady-state levels of WT and mutant CBS in vivo. Our strategy may be valuable for understanding ligand effects on proteins with a complex architecture and their role in human genetic diseases and for the development of novel pharmacological strategies.

Full or partial burial of ionizable groups in the hydrophobic interior of proteins underlies the large modulation in group properties (modified pK value, high nucleophilicity, enhanced capability of interaction with chemical moieties of the substrate, etc.) linked to biological function. Indeed, the few internal ionizable residues found in proteins are known to play important functional roles in catalysis and, in general, in energy transduction processes. However, ionizable-group burial is expected to be seriously disruptive and, it is important to note, most functional sites contain not just one, but several ionizable residues. Hence, the adaptations involved in the development of function in proteins (through in vitro engineering or during the course of natural evolution) are not fully understood. Here, we explore experimentally how proteins respond to the accumulation of hydrophobic-to-ionizable residue substitutions. For this purpose, we have constructed a combinatorial library targeting a hydrophobic cluster in a consensus-engineered, stabilized form of a small model protein. Contrary to naive expectation, half of the variants randomly selected from the library are soluble, folded, and active, despite including up to four mutations. Furthermore, for these variants, the dependence of stability with the number of mutations is not synergistic and catastrophic, but smooth and approximately linear. Clearly, stabilized protein scaffolds may be robust enough to withstand many disruptive hydrophobic-to-ionizable residue mutations, even when they are introduced in the same region of the structure. These results should be relevant for protein engineering and may have implications for the understanding of the early evolution of enzymes. Proteins 2012; © 2011 Wiley Periodicals, Inc.

Many aspects of protein folding can be understood in terms of projections of the highly dimensional energy landscape onto a few (or even only one) particularly relevant coordinates. These free-energy surfaces can be probed conveniently from experimental differential scanning calorimetry (DSC) thermograms, as DSC provides a direct relation with the protein partition function. Free-energy surfaces thus obtained are consistent with two fundamental scenarios predicted by the energy-landscape perspective: (a) well-defined macrostates separated by significant free-energy barriers, in some cases, and, in many other cases, (b) marginal or even vanishingly small barriers, which furthermore show a good correlation with kinetics for fast- and ultrafast-folding proteins. Overall, the potential of DSC to assess free-energy surfaces for a wide variety of proteins makes it possible to address fundamental issues, such as the molecular basis of the barrier modulations produced by natural selection in response to functional requirements or to ensure kinetic stability.

It is possible to travel back in time at the molecular level by reconstructing proteins from extinct organisms. Here we report the reconstruction, based on sequence predicted by phylogenetic analysis, of seven Precambrian thioredoxin enzymes (Trx) dating back between ~1.4 and ~4 billion years (Gyr). The reconstructed enzymes are up to 32 °C more stable than modern enzymes, and the oldest show markedly higher activity than extant ones at pH 5. We probed the mechanisms of reduction of these enzymes using single-molecule force spectroscopy. From the force dependency of the rate of reduction of an engineered substrate, we conclude that ancient Trxs use chemical mechanisms of reduction similar to those of modern enzymes. Although Trx enzymes have maintained their reductase chemistry unchanged, they have adapted over 4 Gyr to the changes in temperature and ocean acidity that characterize the evolution of the global environment from ancient to modern Earth.

Enzymes are dynamic entities, and their dynamic properties are clearly linked to their biological function. It follows that dynamics ought to play an essential role in enzyme evolution. Indeed, a link between conformational diversity and the emergence of new enzyme functionalities has been recognized for many years. However, it is only recently that state-of-the-art computational and experimental approaches are revealing the crucial molecular details of this link. Specifically, evolutionary trajectories leading to functional optimization for a given host environment or to the emergence of a new function typically involve enriching catalytically competent conformations and/or the freezing out of non-competent conformations of an enzyme. In some cases, these evolutionary changes are achieved through distant mutations that shift the protein ensemble towards productive conformations. Multifunctional intermediates in evolutionary trajectories are probably multi-conformational, i.e. able to switch between different overall conformations, each competent for a given function. Conformational diversity can assist the emergence of a completely new active site through a single mutation by facilitating transition-state binding. We propose that this mechanism may have played a role in the emergence of enzymes at the primordial, progenote stage, where it was plausibly promoted by high environmental temperatures and the possibility of additional phenotypic mutations.

Protein physicochemical properties must undergo complex changes during evolution, as a response to modifications in the organism environment, the result of the proteins taking up new roles or because of the need to cope with the evolution of molecular interacting partners. Recent work has emphasized the role of stability and stability–function trade-offs in these protein adaptation processes. In the present study, on the other hand, we report that combinations of a few conservative, high-frequency-of-fixation mutations in the thioredoxin molecule lead to largely independent changes in both stability and the diversity of catalytic mechanisms, as revealed by single-molecule atomic force spectroscopy. Furthermore, the changes found are evolutionarily significant, as they combine typically hyperthermophilic stability enhancements with modulations in function that span the ranges defined by the quite different catalytic patterns of thioredoxins from bacterial and eukaryotic origin. These results suggest that evolutionary protein adaptation may use, in some cases at least, the potential of conservative mutations to originate a multiplicity of evolutionarily allowed mutational paths leading to a variety of protein modulation patterns. In addition the results support the feasibility of using evolutionary information to achieve protein multi-feature optimization, an important biotechnological goal.

A large collection of structural snapshots along a full catalytic cycle of Escherichia coli thioredoxin reductase (TrxR) has been generated and characterized using a combination of theoretical methods. Molecular models were built starting from the available X-ray crystallographic structures of dimeric wild-type TrxR in the flavin-oxidizing conformation and a C135S TrxR mutant enzyme in a flavin-reducing conformation “trapped” by a cross-link between Cys138 of TrxR and Cys32 of C35S mutant thioredoxin (Trx). The transition between these two extreme states, which is shown to be reproduced in a normal mode analysis, as well as natural cofactor binding and dissociation, were simulated for the wild-type species using unrestrained and targeted molecular dynamics following docking of oxidized Trx to reduced TrxR. The whole set of simulations provides a comprehensive structural framework for understanding the mechanism of disulfide reduction in atomic detail and identifying the most likely intermediates that facilitate entry of NADPH and exit of NADP+. The crucial role assigned to Arg73 and Lys36 of Trx in substrate binding and complex stabilization was ascertained when R73G, R73D, and K36A site-directed mutants of Trx were shown to be impaired to different extents in their ability to be reduced by TrxR. On the basis of previous findings and the results reported herein, E. coli TrxR appears as a beautifully engineered molecular machine that is capable of synchronizing cofactor capture and ejection with substrate binding and redox activity through an interdomain twisting motion. Proteins 2010. © 2009 Wiley-Liss, Inc.

It is widely recognized that enhancement of protein stability is an important biotechnological goal. However, some applications at least, could actually benefit from stability being strongly dependent on a suitable environment variable, in such a way that enhanced stability or decreased stability could be realized as required. In therapeutic applications, for instance, a long shelf-life under storage conditions may be convenient, but a sufficiently fast degradation of the protein after it has performed the planned molecular task in vivo may avoid side effects and toxicity. Undesirable effects associated to high stability are also likely to occur in food-industry applications. Clearly, one fundamental factor involved here is the kinetic stability of the protein, which relates to the time-scale of the irreversible denaturation processes and which is determined to some significant extent by the free-energy barrier for unfolding (the barrier that “separates” the native state from the highly-susceptible-to-irreversible-alterations nonnative states). With an appropriate experimental model, we show that strong environment-dependencies of the thermodynamic and kinetic stabilities can be achieved using robust protein engineering. We use sequence-alignment analysis and simple computational electrostatics to design stabilizing and destabilizing mutations, the latter introducing interactions between like charges which are screened out at high salt. Our design procedures lead naturally to mutating regions which are mostly unstructured in the transition state for unfolding. As a result, the large salt effect on the thermodynamic stability of our consensus plus charge-reversal variant translates into dramatic changes in the time-scale associated to the unfolding barrier: from the order of years at high salt to the order of days at low salt. Certainly, large changes in salt concentration are not expected to occur in biological systems in vivo. Hence, proteins with strong salt-dependencies of the thermodynamic and kinetic stabilities are more likely to be of use in those cases in which high-stability is required only under storage conditions. A plausible scenario is that inclusion of high salt in liquid formulations will contribute to a long protein shelf-life, while the lower salt concentration under the conditions of the application will help prevent the side effects associated with high-stability which may potentially arise in some therapeutic and food-industry applications. From a more general viewpoint, this work shows that consensus engineering and electrostatic engineering can be readily combined and clarifies relevant aspects of the relation between thermodynamic stability and kinetic stability in proteins. Proteins 2008. © 2007 Wiley-Liss, Inc.

Differential scanning calorimetry (DSC) is a very powerful tool for investigating protein folding and stability because its experimental output reflects the energetics of all conformations that become minimally populated during thermal unfolding. Accordingly, analysis of DSC experiments with simple thermodynamic models has been key for developing our understanding of protein stability during the past five decades. The discovery of ultrafast folding proteins, which have naturally broad conformational ensembles and minimally cooperative unfolding, opens the possibility of probing the complete folding free energy landscape, including those conformations at the top of the barrier to folding, via DSC. Exploiting this opportunity requires high-quality experiments and the implementation of novel analytical methods based on statistical mechanics. Here, we cover the recent exciting developments in this front, describing the new analytical procedures in detail as well as providing experimental guidelines for performing such analysis.

We have found that, contrary to naïve intuition, the degree of operational reversibility in the thermal denaturation of lipase from Thermomyces lanuginosa (an important industrial enzyme) in urea solutions is maximum when the protein is heated several degrees above the end of the temperature-induced denaturation transition. Upon cooling to room temperature, the protein seems to reach a state with enzymatic activity similar to that of the initial native state, but with higher denaturation temperature and radically different behavior in terms of susceptibility to irreversible denaturation. These results show that patterns of operational reversibility/irreversibility in protein denaturation may be more complex than the often-taken-for-granted, two-situation classification (reversible vs. irreversible). Furthermore, they are consistent with the possibility of existence of different native or native-like states separated by high kinetic barriers under native conditions and they suggest experimental procedures to reach and study such "alternative" native states.

It appears plausible that natural selection constrains, to some extent at least, the stability in many natural proteins. If, during protein evolution, stability fluctuates within a comparatively narrow range, then mutations are expected to be fixed with frequencies that reflect mutational effects on stability. Indeed, we recently reported a robust correlation between the effect of 27 conservative mutations on the thermodynamic stability (unfolding free energy) of Escherichia coli thioredoxin and the frequencies of residues occurrences in sequence alignments. We show here that this correlation likely implies a lower limit to thermodynamic stability of only a few kJ/mol below the unfolding free energy of the wild-type (WT) protein. We suggest, therefore, that the correlation does not reflect natural selection of thermodynamic stability by itself, but of some other factor which is linked to thermodynamic stability for the mutations under study. We propose that this other factor is the kinetic stability of thioredoxin in vivo, since( i) kinetic stability relates to irreversible denaturation, (ii) the rate of irreversible denaturation in a crowded cellular environment (or in a harsh extracellular environment) is probably determined by the rate of unfolding, and (iii) the half-life for unfolding changes in an exponential manner with activation free energy and, consequently, comparatively small free energy effects can have deleterious consequences for kinetic stability. This proposal is supported by the results of a kinetic study of the WT form and the 27 single-mutant variants of E. coli thioredoxin based on the global analyses of chevron plots and equilibrium unfolding profiles determined from double-jump unfolding assays. This kinetic study suggests, furthermore, one of the factors that may contribute to the high activation free energy for unfolding in thioredoxin (required for kinetic stability), namely the energetic optimization of native-state residue environments in regions, which become disrupted in the transition state for unfolding.