Engineering ancestral protein hyperstability

Characterization, structural, and evolutionary analysis of an extremophilic GH5 endoglucanase from Bacillus sp. G131: Insights from ancestral sequence reconstruction

Nature has developed extremozymes that catalyze complex reaction processes in extreme environmental conditions. Accordingly, a combined approach consisting of extremozyme screening, ancestral sequence resurrection (ASR), and molecular dynamic simulation was utilized to construct a developed endoglucanase. The primary experimental and in-silico data led to the prediction of a hypothetical sequence of endoglucanase (EG5-G131) using Bacillus sp. G131 confirmed by amplification and sequencing. EG5-G131 exhibited noticeable stability in a broad-pH range, several detergents, organic solvents, and temperatures up to 80 °C. The molecular weight, Vmax, and Km of the purified endoglucanase were estimated to be 36 kDa, 4.32 μmol/min, and 23.62 mg/ml, respectively. The calculated thermodynamic parameters for EG5-G131 confirmed its intrinsic thermostability. Computational analysis revealed Glu142 and Glu230 as active-site residues of the enzyme. Furthermore, the enzyme remained bound to cellotetraose at 298 K, 333 K, 343 K, and 353 K for 300 ns, consistent with our experimental data. ASR of EG5-G131 led to the introduction of ancestral ANC204 and ANC205, which show similar thermodynamic characteristics with the last Firmicute common ancestor. Finally, truncating loops from the N-terminal of two sequences created two variants with desirable thermal stability, suggesting the evolutionary deciphering of the functional domain of the GH5 family in Bacillus sp. G131.

Thr-to-Ala Mutation Leads to a Larger Aromatic Pair and Reduced Packing Density in α1,α3-Helices during Thioredoxin Cold Adaptation

This study investigates the impact of aromatic-aromatic interactions on the cold adaptation of thioredoxin (Trx), a small redox protein with a conserved Trx-fold structure. Two Trx orthologs, one from the psychrophilic Arctic bacterium Sphingomonas sp. (SpTrx) and the other from the mesophilic Escherichia coli (EcTrx), display distinct aromatic interactions in their α1,α3-helices. SpTrx features a larger Trp11-Phe69 pair, while EcTrx employs a smaller Phe12-Tyr70 pair along with an additional Asp9-Thr66 hydrogen bond. Smaller aromatic residues in SpTrx (Phe-Phe or Phe-Tyr pair) lead to decreased thermal and thermodynamic stabilities, increased conformational flexibility, and reduced enzyme activity. In contrast, EcTrx's thermal stability is primarily influenced by the larger Trp residue, especially in the more hydrophobic Trp-Phe pair compared to the Trp-Tyr pair. Both SpTrx and EcTrx exhibit a strengthening of the Asp-Thr hydrogen bond by a Phe-Tyr pair and a weakening by a Trp-Phe pair. Additionally, the Asp8-Thr65 hydrogen bond in SpTrx contributes to the destabilization of the Phe-Phe pair. Molecular dynamics simulations of SpTrx indicate that a smaller aromatic pair or the Asp-Thr hydrogen bond in the α1,α3-helices further destabilizes the α2-helix across the central β-sheet. Our results suggest that the Thr-to-Ala mutation destabilizes the α1,α3-helices, resulting in a larger aromatic pair and reduced packing density in psychrophilic Trxs during cold adaptation. These findings enhance our understanding of Trx's adaptation to colder temperatures.

Active site center redesign increases protein stability preserving catalysis in thioredoxin

The stabilization of natural proteins is a long-standing desired goal in protein engineering. Optimizing the hydrophobicity of the protein core often results in extensive stability enhancements. However, the presence of totally or partially buried catalytic charged residues, essential for protein function, has limited the applicability of this strategy. Here, focusing on the thioredoxin, we aimed to augment protein stability by removing buried charged residues in the active site without loss of catalytic activity. To this end, we performed a charged-to-hydrophobic substitution of a buried and functional group, resulting in a significant stability increase yet abolishing catalytic activity. Then, to simulate the catalytic role of the buried ionizable group, we designed a combinatorial library of variants targeting a set of seven surface residues adjacent to the active site. Notably, more than 50% of the library variants restored, to some extent, the catalytic activity. The combination of experimental study of 2% of the library with the prediction of the whole mutational space by partial least squares regression revealed that a single point mutation at the protein surface is sufficient to fully restore the catalytic activity without thermostability cost. As a result, we engineered one of the highest thermal stabilities reported for a protein with a natural occurring fold (137°C). Further, our hyperstable variant preserves the catalytic activity both in vitro and in vivo.

The dynamic properties of a nuclear coactivator binding domain are evolutionarily conserved

Evolution of proteins is constrained by their structure and function. While there is a consensus that the plasticity of intrinsically disordered proteins relaxes the structural constraints on evolution there is a paucity of data on the molecular details of these processes. The Nuclear Coactivator Binding Domain (NCBD) from CREB-binding protein is a protein interaction domain, which contains a hydrophobic core but is not behaving as a typical globular domain, and has been described as 'molten-globule like'. The highly dynamic properties of NCBD makes it an interesting model system for evolutionary structure-function investigation of intrinsically disordered proteins. We have here compared the structure and biophysical properties of an ancient version of NCBD present in a bilaterian animal ancestor living around 600 million years ago with extant human NCBD. Using a combination of NMR spectroscopy, circular dichroism and kinetics we show that although NCBD has increased its thermodynamic stability, it has retained its dynamic biophysical properties in the ligand-free state in the evolutionary lineage leading from the last common bilaterian ancestor to humans. Our findings suggest that the dynamic properties of NCBD have been maintained by purifying selection and thus are important for its function, which includes mediating several distinct protein-protein interactions.

Combining Ancestral Reconstruction with Folding-Landscape Simulations to Engineer Heterologous Protein Expression

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.

Exploiting enzyme evolution for computational protein design

Recent years have seen an explosion of interest in understanding the physicochemical parameters that shape enzyme evolution, as well as substantial advances in computational enzyme design. This review discusses three areas where evolutionary information can be used as part of the design process: (i) using ancestral sequence reconstruction (ASR) to generate new starting points for enzyme design efforts; (ii) learning from how nature uses conformational dynamics in enzyme evolution to mimic this process in silico; and (iii) modular design of enzymes from smaller fragments, again mimicking the process by which nature appears to create new protein folds. Using showcase examples, we highlight the importance of incorporating evolutionary information to continue to push forward the boundaries of enzyme design studies.

Non-conservation of folding rates in the thioredoxin family reveals degradation of ancestral unassisted-folding

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.

Structures of hyperstable ancestral haloalkane dehalogenases show restricted conformational dynamics

Ancestral sequence reconstruction is a powerful method for inferring ancestors of modern enzymes and for studying structure-function relationships of enzymes. We have previously applied this approach to haloalkane dehalogenases (HLDs) from the subfamily HLD-II and obtained thermodynamically highly stabilized enzymes (ΔT m up to 24 °C), showing improved catalytic properties. Here we combined crystallographic structural analysis and computational molecular dynamics simulations to gain insight into the mechanisms by which ancestral HLDs became more robust enzymes with novel catalytic properties. Reconstructed ancestors exhibited similar structure topology as their descendants with the exception of a few loop deviations. Strikingly, molecular dynamics simulations revealed restricted conformational dynamics of ancestral enzymes, which prefer a single state, in contrast to modern enzymes adopting two different conformational states. The restricted dynamics can potentially be linked to their exceptional stabilization. The study provides molecular insights into protein stabilization due to ancestral sequence reconstruction, which is becoming a widely used approach for obtaining robust protein catalysts.

Molecular origins of folding rate differences in the thioredoxin family

Thioredoxins are a family of conserved oxidoreductases responsible for maintaining redox balance within cells. They have also served as excellent model systems for protein design and engineering studies particularly through ancestral sequence reconstruction methods. The recent work by Gamiz-Arco et al. [Biochem J (2019) 476, 3631–3647] answers fundamental questions on how specific sequence differences can contribute to differences in folding rates between modern and ancient thioredoxins but also among a selected subset of modern thioredoxins. They surprisingly find that rapid unassisted folding, a feature of ancient thioredoxins, is not conserved in the modern descendants suggestive of co-evolution of better folding machinery that likely enabled the accumulation of mutations that slow-down folding. The work thus provides an interesting take on the expected folding-stability-function constraint while arguing for additional factors that contribute to sequence evolution and hence impact folding efficiency.

The Role of Conformational Dynamics and Allostery in Modulating Protein Evolution

Advances in sequencing techniques and statistical methods have made it possible not only to predict sequences of ancestral proteins but also to identify thousands of mutations in the human exome, some of which are disease associated. These developments have motivated numerous theories and raised many questions regarding the fundamental principles behind protein evolution, which have been traditionally investigated horizontally using the tip of the phylogenetic tree through comparative studies of extant proteins within a family. In this article, we review a vertical comparison of the modern and resurrected ancestral proteins. We focus mainly on the dynamical properties responsible for a protein's ability to adapt new functions in response to environmental changes. Using the Dynamic Flexibility Index and the Dynamic Coupling Index to quantify the relative flexibility and dynamic coupling at a site-specific, single-amino-acid level, we provide evidence that the migration of hinges, which are often functionally critical rigid sites, is a mechanism through which proteins can rapidly evolve. Additionally, we show that disease-associated mutations in proteins often result in flexibility changes even at positions distal from mutational sites, particularly in the modulation of active site dynamics.

Functional analyses of ancestral thioredoxins provide insights into their evolutionary history

Thioredoxin (Trx) is a conserved, cytosolic reductase in all known organisms. The enzyme receives two electrons from NADPH via thioredoxin reductase (TrxR) and passes them on to multiple cellular reductases via disulfide exchange. Despite the ubiquity of thioredoxins in all taxa, little is known about the functions of resurrected ancestral thioredoxins in the context of a modern mesophilic organism. Here, we report on functional in vitro and in vivo analyses of seven resurrected Precambrian thioredoxins, dating back 1–4 billion years, in the Escherichia coli cytoplasm. Using synthetic gene constructs for recombinant expression of the ancestral enzymes, along with thermodynamic and kinetic assays, we show that all ancestral thioredoxins, as today's thioredoxins, exhibit strongly reducing redox potentials, suggesting that thioredoxins served as catalysts of cellular reduction reactions from the beginning of evolution, even before the oxygen catastrophe. A detailed, quantitative characterization of their interactions with the electron donor TrxR from Escherichia coli and the electron acceptor methionine sulfoxide reductase, also from E. coli, strongly hinted that thioredoxins and thioredoxin reductases co-evolved and that the promiscuity of thioredoxins toward downstream electron acceptors was maintained during evolution. In summary, our findings suggest that thioredoxins evolved high specificity for their sole electron donor TrxR while maintaining promiscuity to their multiple electron acceptors.

Selection and analyses of variants of a designed protein suggest importance of hydrophobicity of partially buried sidechains for protein stability at high temperatures

Computationally designed proteins of high stability provide specimen in addition to natural proteins for the study of sequence-structure stability relationships at the very high end of protein stability spectrum. The melting temperature of E_1r26, a protein we previously designed using the A Backbone-based Amino aCid Usage Survey (ABACUS) sequence design program, is above 110 °C, more than 50 °C higher than that of the natural thioredoxin protein whose backbone (PDB ID 1R26) has been used as the design target. Using an experimental selection approach, we obtained variants of E_1r26 that remain folded but are of reduced stability, including one whose unfolding temperature and denaturing guanidine concentration are similar to those of 1r26. The mutant unfolds with a certain degree of cooperativity. Its structure solved by X-ray crystallography agrees with that of 1r26 by a root mean square deviation of 1.3 Å, adding supports to the accuracy of the ABACUS method. Analyses of intermediate mutants indicate that the substitution of two partially buried hydrophobic residues (isoleucine and leucine) by polar residues (threonine and serine, respectively) are responsible for the dramatic change in the unfolding temperature. It is suggested that the effects of mutations located in rigid secondary structure regions, but not those in loops, may be well predicted through ABACUS mutation energy analysis. The results also suggest that hydrophobic effects involving intermediately buried sidechains can be critically important for protein stability at high temperatures.

Data-driven engineering of protein therapeutics

Protein therapeutics requires a series of properties beyond biochemical activity, including serum stability, low immunogenicity, and manufacturability. Mutations that improve one property often decrease one or more of the other essential requirements for therapeutic efficacy, making the protein engineering challenge difficult. The past decade has seen an explosion of new techniques centered around cheaply reading and writing DNA. This review highlights the recent use of such high throughput technologies for engineering protein therapeutics. Examples include the use of human antibody repertoire sequence data to pair antibody heavy and light chains, comprehensive mutational analysis for engineering antibody specificity, and the use of ancestral and inter-species sequence data to engineer simultaneous improvements in enzyme catalytic efficiency and stability. We conclude with a perspective on further ways to integrate mature protein engineering pipelines with the exponential increases in the volume of sequencing data expected in the forthcoming decade.

Protein engineers turned evolutionists-the quest for the optimal starting point

The advent of laboratory directed evolution yielded a fruitful crosstalk between the disciplines of molecular evolution and bio-engineering. Here, we outline recent developments in both disciplines with respect to how one can identify the best starting points for directed evolution, such that highly efficient and robust tailor-made enzymes can be obtained with minimal optimization. Directed evolution studies have highlighted essential features of engineer-able enzymes: highly stable, mutationally robust enzymes with the capacity to accept a broad range of substrates. Robust, evolvable enzymes can be inferred from the natural sequence record. Broad substrate spectrum relates to conformational plasticity and can also be predicted by phylogenetic analyses and/or by computational design. Overall, an increasingly powerful toolkit is becoming available for identifying optimal starting points including network analyses of enzyme superfamilies and other bioinformatics methods.

Ancestral Sequence Reconstruction as a Tool for the Elucidation of a Stepwise Evolutionary Adaptation

Ancestral sequence reconstruction (ASR) is a powerful tool to infer primordial sequences from contemporary, i.e., extant ones. An essential element of ASR is the computation of a phylogenetic tree whose leaves are the chosen extant sequences. Most often, the reconstructed sequence related to the root of this tree is of greatest interest: It represents the common ancestor (CA) of the sequences under study. If this sequence encodes a protein, one can "resurrect" the CA by means of gene synthesis technology and study biochemical properties of this extinct predecessor with the help of wet-lab experiments.However, ASR deduces also sequences for all internal nodes of the tree, and the well-considered analysis of these "intermediates" can help to elucidate evolutionary processes. Moreover, one can identify key mutations that alter proteins or protein complexes and are responsible for the differing properties of extant proteins. As an illustrative example, we describe the protocol for the rapid identification of hotspots determining the binding of the two subunits within the heteromeric complex imidazole glycerol phosphate synthase.

Ancestral diterpene cyclases show increased thermostability and substrate acceptance

Bacterial diterpene cyclases are receiving increasing attention in biocatalysis and synthetic biology for the sustainable generation of complex multicyclic building blocks. Herein, we explore the potential of ancestral sequence reconstruction (ASR) to generate remodeled cyclases with enhanced stability, activity, and promiscuity. Putative ancestors of spiroviolene synthase, a bacterial class I diterpene cyclase, display an increased yield of soluble protein of up to fourfold upon expression in the model organism Escherichia coli. Two of the resurrected enzymes, with an estimated age of approximately 1.7 million years, display an upward shift in thermostability of 7-13 °C. Ancestral spiroviolene synthases catalyze cyclization of the natural C₂₀ -substrate geranylgeranyl diphosphate (GGPP) and also accept C₁₅ farnesyl diphosphate (FPP), which is not converted by the extant enzyme. In contrast, the consensus sequence generated from the corresponding multiple sequence alignment was found to be inactive toward both substrates. Mutation of a nonconserved position within the aspartate-rich motif of the reconstructed ancestral cyclases was associated with modest effects on activity and relative substrate specificity (i.e., k_cat /K_M for GGPP over k_cat /K_M for FPP). Kinetic analyses performed at different temperatures reveal a loss of substrate saturation, when going from the ancestor with highest thermostability to the modern enzyme. The kinetics data also illustrate how an increase in temperature optimum of biocatalysis is reflected in altered entropy and enthalpy of activation. Our findings further highlight the potential and limitations of applying ASR to biosynthetic machineries in secondary metabolism.

Ancient thioredoxins evolved to modern-day stability-function requirement by altering native state ensemble

Thioredoxins (THRXs)-small globular proteins that reduce other proteins-are ubiquitous in all forms of life, from Archaea to mammals. Although ancestral thioredoxins share sequential and structural similarity with the modern-day (extant) homologues, they exhibit significantly different functional activity and stability. We investigate this puzzle by comparative studies of their (ancient and modern-day THRXs') native state ensemble, as quantified by the dynamic flexibility index (DFI), a metric for the relative resilience of an amino acid to perturbations in the rest of the protein. Clustering proteins using DFI profiles strongly resemble an alternative classification scheme based on their activity and stability. The DFI profiles of the extant proteins are substantially different around the α3, α4 helices and catalytic regions. Likewise, allosteric coupling of the active site with the rest of the protein is different between ancient and extant THRXs, possibly explaining the decreased catalytic activity at low pH with evolution. At a global level, we note that the population of low-flexibility (called hinges) and high-flexibility sites increases with evolution. The heterogeneity (quantified by the variance) in DFI distribution increases with the decrease in the melting temperature typically associated with the evolution of ancient proteins to their modern-day counterparts.This article is part of a discussion meeting issue 'Allostery and molecular machines'.

Biotechnological and protein-engineering implications of ancestral protein resurrection

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.

Making sense of the past: hyperstability of ancestral thioredoxins explained by free energy simulations

Thioredoxin (Trx), a small and globular protein, is present in all kinds of organisms, from Archea to higher mammals. Throughout evolution, the Trx sequence has undergone subtle modifications to adapt to varying environmental conditions. The high degree of sequence conservation makes Trx very amenable to ancestral protein reconstruction techniques. In this work, we address the study of the structural and energetic determinants of thermostability in E. coli Trx using a dataset of mutations inspired by ancestral reconstruction. We compute, from first principles, the expected contribution of 19 different amino acid substitutions to the stability (ΔΔG) and the melting temperature (ΔT_m) of the protein. We also describe the specific changes in structure and protein dynamics responsible for the stabilizing or destabilizing effects of these mutations. Our results point to local and independent changes for most of the variants. Our predictions are accurate enough to substantiate the proposal of new hypotheses regarding evolutionary relationships between mutations, as in the case of T89R, P68A and G74S or K90L and F102A, and reach beyond the initial set to suggest improved variants, such as K90I or K90Y.

Enzyme engineering: reaching the maximal catalytic efficiency peak

The practical need for highly efficient enzymes presents new challenges in enzyme engineering, in particular, the need to improve catalytic turnover (k_cat) or efficiency (k_cat/K_M) by several orders of magnitude. However, optimizing catalysis demands navigation through complex and rugged fitness landscapes, with optimization trajectories often leading to strong diminishing returns and dead-ends. When no further improvements are observed in library screens or selections, it remains unclear whether the maximal catalytic efficiency of the enzyme (the catalytic 'fitness peak') has been reached; or perhaps, an alternative combination of mutations exists that could yield additional improvements. Here, we discuss fundamental aspects of the process of catalytic optimization, and offer practical solutions with respect to overcoming optimization plateaus.