Thermodynamic system drift in protein evolution

Insights into the low-temperature adaptation of an enzyme as studied through ancestral sequence reconstruction

For billions of years, enzymes have evolved in response to the changing environments in which their host organisms lived. Various lines of evidence suggest the earliest primitive organisms inhabited high-temperature environments and possessed enzymes adapted to such conditions. Consequently, extant mesophilic and psychrophilic enzymes are believed to have adapted to lower temperatures during the evolutionary process. Herein, we analyzed this low-temperature adaptation using ancestral sequence reconstruction. Previously, we generated the phylogenetic tree of 3-isopropylmalate dehydrogenases (IPMDHs) and reconstructed the sequence of the last bacterial common ancestor. The corresponding ancestral enzyme displayed high thermostability and catalytic activity at elevated temperatures but moderate activity at low temperatures (Furukawa et al., Sci. Rep., 2020;10:15493). Here, to identify amino acid residues that are responsible for the low-temperature adaptation, we reconstructed and characterized all 11 evolutionary intermediates that sequentially connect the last bacterial common ancestor with extant mesophilic IPMDH from Escherichia coli. A remarkable change in catalytic properties, from those suited for high reaction temperatures to those adapted for low temperatures, occurred between two consecutive evolutionary intermediates. Using a combination of sequence comparisons between ancestral proteins and site-directed mutagenesis analyses, three key amino acid substitutions were identified that enhance low-temperature catalytic activity. Intriguingly, amino acid substitutions that had the most significant impact on activity at low temperatures displayed no discernable effect on thermostability. However, these substitutions markedly reduced the activation energy for catalysis, thereby improving low-temperature activity. The results were further investigated by molecular dynamics simulations of the predicted structures of the ancestral enzymes. Our findings exemplify how ancestral sequence reconstruction can identify residues crucial for adaptation to low temperatures.

Evolutionary rewiring of the dynamic network underpinning allosteric epistasis in NS1 of the influenza A virus

Viral proteins frequently mutate to evade host innate immune responses, yet the impact of these mutations on the molecular energy landscape remains unclear. Epistasis, the intramolecular communications between mutations, often renders the combined mutational effects unpredictable. Nonstructural protein 1 (NS1) is a major virulence factor of the influenza A virus (IAV) that activates host PI3K by binding to its p85β subunit. Here, we present a deep analysis of the impact of evolutionary mutations in NS1 that emerged between the 1918 pandemic IAV strain and its descendant PR8 strain. Our analysis reveals how the mutations rewired interresidue communications, which underlie long-range allosteric and epistatic networks in NS1. Our findings show that PR8 NS1 binds to p85β with approximately 10-fold greater affinity than 1918 NS1 due to allosteric mutational effects, which are further tuned by epistasis. NMR chemical shift perturbation and methyl-axis order parameter analyses revealed that the mutations induced long-range structural and dynamic changes in PR8 NS1, relative to 1918 NS1, enhancing its affinity to p85β. Complementary molecular dynamics simulations and graph theory-based network analysis for conformational dynamics on the submicrosecond timescales uncover how these mutations rewire the dynamic network, which underlies the allosteric epistasis. Significantly, we find that conformational dynamics of residues with high betweenness centrality play a crucial role in communications between network communities and are highly conserved across influenza A virus evolution. These findings advance our mechanistic understanding of the allosteric and epistatic communications between distant residues and provide insight into their role in the molecular evolution of NS1.

Protein stability is determined by single-site bias rather than pairwise covariance

The biases revealed in protein sequence alignments have been shown to provide information related to protein structure, stability, and function. For example, sequence biases at individual positions can be used to design consensus proteins that are often more stable than naturally occurring counterparts. Likewise, correlations between pairs of residue can be used to predict protein structures. Recent work using Potts models show that together, single-site biases and pair correlations lead to improved predictions of protein fitness, activity, and stability. Here we use a Potts model to design groups of protein sequences with different amounts of single-site biases and pair correlations, and determine the thermodynamic stabilities of a representative set of sequences from each group. Surprisingly, sequences excluding pair correlations maximize stability, whereas sequences that maximize pair correlations are less stable, suggesting that pair correlations contribute to another aspect of protein fitness. Consistent with this interpretation, we find that for adenylate kinase, enzyme activity is greatly increased by maximizing pair correlations. The finding that elimination of covariant residue pairs increases protein stability suggests a route to enhance stability of designed proteins; indeed, this strategy produces hyperstable homeodomain and adenylate kinase proteins that retain significant activity.

Deciphering Structural Traits for Thermal and Kinetic Stability across Protein Family Evolution through Ancestral Sequence Reconstruction

Natural proteins are frequently marginally stable, and an increase in environmental temperature can easily lead to unfolding. As a result, protein engineering to improve protein stability is an area of intensive research. Nonetheless, since there is usually a high degree of structural homology between proteins from thermophilic organisms and their mesophilic counterparts, the identification of structural determinants for thermoadaptation is challenging. Moreover, in many cases, it has become clear that the success of stabilization strategies is often dependent on the evolutionary history of a protein family. In the last few years, the use of ancestral sequence reconstruction (ASR) as a tool for elucidation of the evolutionary history of functional traits of a protein family has gained strength. Here, we used ASR to trace the evolutionary pathways between mesophilic and thermophilic kinases that participate in the biosynthetic pathway of vitamin B1 in bacteria. By combining biophysics approaches, X-ray crystallography, and molecular dynamics simulations, we found that the thermal stability of these enzymes correlates with their kinetic stability, where the highest thermal/kinetic stability is given by an increase in small hydrophobic amino acids that allow a higher number of interatomic hydrophobic contacts, making this type of interaction the main support for stability in this protein architecture. The results highlight the potential benefits of using ASR to explore the evolutionary history of protein sequence and structure to identify traits responsible for the kinetic and thermal stability of any protein architecture.

Ancestral Reconstruction and the Evolution of Protein Energy Landscapes

A protein's sequence determines its conformational energy landscape. This, in turn, determines the protein's function. Understanding the evolution of new protein functions therefore requires understanding how mutations alter the protein energy landscape. Ancestral sequence reconstruction (ASR) has proven a valuable tool for tackling this problem. In ASR, one phylogenetically infers the sequences of ancient proteins, allowing characterization of their properties. When coupled to biophysical, biochemical, and functional characterization, ASR can reveal how historical mutations altered the energy landscape of ancient proteins, allowing the evolution of enzyme activity, altered conformations, binding specificity, oligomerization, and many other protein features. In this article, we review how ASR studies have been used to dissect the evolution of energy landscapes. We also discuss ASR studies that reveal how energy landscapes have shaped protein evolution. Finally, we propose that thinking about evolution from the perspective of an energy landscape can improve how we approach and interpret ASR studies.

Rheostatic contributions to protein stability can obscure a position's functional role

Rheostat positions, which can be substituted with various amino acids to tune protein function across a range of outcomes, are a developing area for advancing personalized medicine and bioengineering. Current methods cannot accurately predict which proteins contain rheostat positions or their substitution outcomes. To compare the prevalence of rheostat positions in homologs, we previously investigated their occurrence in two pyruvate kinase (PYK) isozymes. Human liver PYK contained numerous rheostat positions that tuned the apparent affinity for the substrate phosphoenolpyruvate (Kapp-PEP) across a wide range. In contrast, no functional rheostat positions were identified in Zymomonas mobilis PYK (ZmPYK). Further, the set of ZmPYK substitutions included an unusually large number that lacked measurable activity. We hypothesized that the inactive substitution variants had reduced protein stability, precluding detection of Kapp-PEP tuning. Using modified buffers, robust enzymatic activity was obtained for 19 previously-inactive ZmPYK substitution variants at three positions. Surprisingly, both previously-inactive and previously-active substitution variants all had Kapp-PEP values close to wild-type. Thus, none of the three positions were functional rheostat positions, and, unlike human liver PYK, ZmPYK's Kapp-PEP remained poorly tunable by single substitutions. To directly assess effects on stability, we performed thermal denaturation experiments for all ZmPYK substitution variants. Many diminished stability, two enhanced stability, and the three positions showed different thermal sensitivity to substitution, with one position acting as a "stability rheostat." The differences between the two PYK homologs raises interesting questions about the underlying mechanism(s) that permit functional tuning by single substitutions in some proteins but not in others.

The importance of input sequence set to consensus-derived proteins and their relationship to reconstructed ancestral proteins

A protein sequence encodes its energy landscape-all the accessible conformations, energetics, and dynamics. The evolutionary relationship between sequence and landscape can be probed phylogenetically by compiling a multiple sequence alignment of homologous sequences and generating common ancestors via Ancestral Sequence Reconstruction or a consensus protein containing the most common amino acid at each position. Both ancestral and consensus proteins are often more stable than their extant homologs-questioning the differences between them and suggesting that both approaches serve as general methods to engineer thermostability. We used the Ribonuclease H family to compare these approaches and evaluate how the evolutionary relationship of the input sequences affects the properties of the resulting consensus protein. While the consensus protein derived from our full Ribonuclease H sequence alignment is structured and active, it neither shows properties of a well-folded protein nor has enhanced stability. In contrast, the consensus protein derived from a phylogenetically-restricted set of sequences is significantly more stable and cooperatively folded, suggesting that cooperativity may be encoded by different mechanisms in separate clades and lost when too many diverse clades are combined to generate a consensus protein. To explore this, we compared pairwise covariance scores using a Potts formalism as well as higher-order sequence correlations using singular value decomposition (SVD). We find the SVD coordinates of a stable consensus sequence are close to coordinates of the analogous ancestor sequence and its descendants, whereas the unstable consensus sequences are outliers in SVD space.

Engineering Biocatalysts for the C-H Activation of Fatty Acids by Ancestral Sequence Reconstruction

Selective, one-step C-H activation of fatty acids from biomass is an attractive concept in sustainable chemistry. Biocatalysis has shown promise for generating high-value hydroxy acids, but to date enzyme discovery has relied on laborious screening and produced limited hits, which predominantly oxidise the subterminal positions of fatty acids. Herein we show that ancestral sequence reconstruction (ASR) is an effective tool to explore the sequence-activity landscape of a family of multidomain, self-sufficient P450 monooxygenases. We resurrected 11 catalytically active CYP116B ancestors, each with a unique regioselectivity fingerprint that varied from subterminal in the older ancestors to mid-chain in the lineage leading to the extant, P450-TT. In lineages leading to extant enzymes in thermophiles, thermostability increased from ancestral to extant forms, as expected if thermophily had arisen de novo. Our studies show that ASR can be applied to multidomain enzymes to develop active, self-sufficient monooxygenases as regioselective biocatalysts for fatty acid hydroxylation.

Empirical evidence for metabolic drift in plant and algal lipid biosynthesis pathways

Metabolic pathway drift has been formulated as a general principle to help in the interpretation of comparative analyses between biosynthesis pathways. Indeed, such analyses often indicate substantial differences, even in widespread pathways that are sometimes believed to be conserved. Here, our purpose is to check how much this interpretation fits to empirical data gathered in the field of plant and algal biosynthesis pathways. After examining several examples representative of the diversity of lipid biosynthesis pathways, we explain why it is important to compare closely related species to gain a better understanding of this phenomenon. Furthermore, this comparative approach brings us to the question of how much biotic interactions are responsible for shaping this metabolic plasticity. We end up introducing some model systems that may be promising for further exploration of this question.

Selection pressures on evolution of ribonuclease H explored with rigorous free-energy-based design

Understanding natural protein evolution and designing novel proteins are motivating interest in development of high-throughput methods to explore large sequence spaces. In this work, we demonstrate the application of multisite λ dynamics (MSλD), a rigorous free energy simulation method, and chemical denaturation experiments to quantify evolutionary selection pressure from sequence-stability relationships and to address questions of design. This study examines a mesophilic phylogenetic clade of ribonuclease H (RNase H), furthering its extensive characterization in earlier studies, focusing on E. coli RNase H (ecRNH) and a more stable consensus sequence (AncCcons) differing at 15 positions. The stabilities of 32,768 chimeras between these two sequences were computed using the MSλD framework. The most stable and least stable chimeras were predicted and tested along with several other sequences, revealing a designed chimera with approximately the same stability increase as AncCcons, but requiring only half the mutations. Comparing the computed stabilities with experiment for 12 sequences reveals a Pearson correlation of 0.86 and root mean squared error of 1.18 kcal/mol, an unprecedented level of accuracy well beyond less rigorous computational design methods. We then quantified selection pressure using a simple evolutionary model in which sequences are selected according to the Boltzmann factor of their stability. Selection temperatures from 110 to 168 K are estimated in three ways by comparing experimental and computational results to evolutionary models. These estimates indicate selection pressure is high, which has implications for evolutionary dynamics and for the accuracy required for design, and suggests accurate high-throughput computational methods like MSλD may enable more effective protein design.

Comparative analysis of reconstructed ancestral proteins with their extant counterparts suggests primitive life had an alkaline habitat

To understand the origin and early evolution of life it is crucial to establish characteristics of the primordial environment that facilitated the emergence and evolution of life. One important environmental factor is the pH of the primordial environment. Here, we assessed the pH-dependent thermal stabilities of previously reconstructed ancestral nucleoside diphosphate kinases and ribosomal protein uS8s. The selected proteins were likely to be present in ancient organisms such as the last common ancestor of bacteria and that of archaea. We also assessed the thermal stability of homologous proteins from extant acidophilic, neutralophilic, and alkaliphilic microorganisms as a function of pH. Our results indicate that the reconstructed ancestral proteins are more akin to those of extant alkaliphilic bacteria, which display greater stability under alkaline conditions. These findings suggest that the common ancestors of bacterial and archaeal species thrived in an alkaline environment. Moreover, we demonstrate the reconstruction method employed in this study is a valuable technique for generating alkali-tolerant proteins that can be used in a variety of biotechnological and environmental applications.

The importance of input sequence set to consensus-derived proteins and their relationship to reconstructed ancestral proteins

A protein sequence encodes its energy landscape - all the accessible conformations, energetics, and dynamics. The evolutionary relationship between sequence and landscape can be probed phylogenetically by compiling a multiple sequence alignment of homologous sequences and generating common ancestors via Ancestral Sequence Reconstruction or a consensus protein containing the most common amino acid at each position. Both ancestral and consensus proteins are often more stable than their extant homologs - questioning the differences and suggesting that both approaches serve as general methods to engineer thermostability. We used the Ribonuclease H family to compare these approaches and evaluate how the evolutionary relationship of the input sequences affects the properties of the resulting consensus protein. While the overall consensus protein is structured and active, it neither shows properties of a well-folded protein nor has enhanced stability. In contrast, the consensus protein derived from a phylogenetically-restricted region is significantly more stable and cooperatively folded, suggesting that cooperativity may be encoded by different mechanisms in separate clades and lost when too many diverse clades are combined to generate a consensus protein. To explore this, we compared pairwise covariance scores using a Potts formalism as well as higher-order couplings using singular value decomposition (SVD). We find the SVD coordinates of a stable consensus sequence are close to coordinates of the analogous ancestor sequence and its descendants, whereas the unstable consensus sequences are outliers in SVD space.

Deciphering the evolution of flavin-dependent monooxygenase stereoselectivity using ancestral sequence reconstruction

Controlling the selectivity of a reaction is critical for target-oriented synthesis. Accessing complementary selectivity profiles enables divergent synthetic strategies, but is challenging to achieve in biocatalytic reactions given enzymes' innate preferences of a single selectivity. Thus, it is critical to understand the structural features that control selectivity in biocatalytic reactions to achieve tunable selectivity. Here, we investigate the structural features that control the stereoselectivity in an oxidative dearomatization reaction that is key to making azaphilone natural products. Crystal structures of enantiocomplementary biocatalysts guided the development of multiple hypotheses centered on the structural features that control the stereochemical outcome of the reaction; however, in many cases, direct substitutions of active site residues in natural proteins led to inactive enzymes. Ancestral sequence reconstruction (ASR) and resurrection were employed as an alternative strategy to probe the impact of each residue on the stereochemical outcome of the dearomatization reaction. These studies suggest that two mechanisms are active in controlling the stereochemical outcome of the oxidative dearomatization reaction: one involving multiple active site residues in AzaH and the other dominated by a single Phe to Tyr switch in TropB and AfoD. Moreover, this study suggests that the flavin-dependent monooxygenases (FDMOs) adopt simple and flexible strategies to control stereoselectivity, which has led to stereocomplementary azaphilone natural products produced by fungi. This paradigm of combining ASR and resurrection with mutational and computational studies showcases sets of tools for understanding enzyme mechanisms and provides a solid foundation for future protein engineering efforts.

Engineering functional thermostable proteins using ancestral sequence reconstruction

Natural proteins are often only slightly more stable in the native state than the denatured state, and an increase in environmental temperature can easily shift the balance toward unfolding. Therefore, the engineering of proteins to improve protein stability is an area of intensive research. Thermostable proteins are required to withstand industrial process conditions, for increased shelf-life of protein therapeutics, for developing robust 'biobricks' for synthetic biology applications, and for research purposes (e.g., structure determination). In addition, thermostability buffers the often destabilizing effects of mutations introduced to improve other properties. Rational design approaches to engineering thermostability require structural information, but even with advanced computational methods, it is challenging to predict or parameterize all the relevant structural factors with sufficient precision to anticipate the results of a given mutation. Directed evolution is an alternative when structures are unavailable but requires extensive screening of mutant libraries. Recently, however, bioinspired approaches based on phylogenetic analyses have shown great promise. Leveraging the rapid expansion in sequence data and bioinformatic tools, ancestral sequence reconstruction can generate highly stable folds for novel applications in industrial chemistry, medicine, and synthetic biology. This review provides an overview of the factors important for successful inference of thermostable proteins by ancestral sequence reconstruction and what it can reveal about the determinants of stability in proteins.

Energy landscape reshaped by strain-specific mutations underlies epistasis in NS1 evolution of influenza A virus

Elucidating how individual mutations affect the protein energy landscape is crucial for understanding how proteins evolve. However, predicting mutational effects remains challenging because of epistasis—the nonadditive interactions between mutations. Here, we investigate the biophysical mechanism of strain-specific epistasis in the nonstructural protein 1 (NS1) of influenza A viruses (IAVs). We integrate structural, kinetic, thermodynamic, and conformational dynamics analyses of four NS1s of influenza strains that emerged between 1918 and 2004. Although functionally near-neutral, strain-specific NS1 mutations exhibit long-range epistatic interactions with residues at the p85β-binding interface. We reveal that strain-specific mutations reshaped the NS1 energy landscape during evolution. Using NMR spin dynamics, we find that the strain-specific mutations altered the conformational dynamics of the hidden network of tightly packed residues, underlying the evolution of long-range epistasis. This work shows how near-neutral mutations silently alter the biophysical energy landscapes, resulting in diverse background effects during molecular evolution.

Influenza A virus (IAV) nonstructural protein 1 (NS1) is a multifunctional virulence factor that interacts with several host factors such as phosphatidylinositol-3-kinase (PI3K). NS1 binds specifically to the p85β regulatory subunit of PI3K and subsequently activates PI3K signaling. Here, Kim et al. show that functionally near-neutral, strain-specific NS1 mutations lead to variations in binding kinetics to p85β exhibit long-range epistatic interactions. Applying NMR they provide evidence that the structural dynamics of the NS1 hydrophobic core have evolved over time and contributed to epistasis.

Ancestral sequence reconstruction of a cytochrome P450 family involved in chemical defence reveals the functional evolution of a promiscuous, xenobiotic-metabolizing enzyme in vertebrates

The cytochrome P450 family 1 enzymes (CYP1s) are a diverse family of hemoprotein monooxygenases, which metabolize many xenobiotics including numerous environmental carcinogens. However, their historical function and evolution remain largely unstudied. Here we investigate CYP1 evolution via the reconstruction and characterization of the vertebrate CYP1 ancestors. Younger ancestors and extant forms generally demonstrated higher activity toward typical CYP1 xenobiotic and steroid substrates than older ancestors, suggesting significant diversification away from the original CYP1 function. Caffeine metabolism appears to be a recently evolved trait of the CYP1A subfamily, observed in the mammalian CYP1A lineage, and may parallel the recent evolution of caffeine synthesis in multiple separate plant species. Likewise, the aryl hydrocarbon receptor agonist, 6-formylindolo[3,2-b]carbazole (FICZ) was metabolized to a greater extent by certain younger ancestors and extant forms, suggesting that activity toward FICZ increased in specific CYP1 evolutionary branches, a process that may have occurred in parallel to the exploitation of land where UV-exposure was higher than in aquatic environments. As observed with previous reconstructions of P450 enzymes, thermostability correlated with evolutionary age; the oldest ancestor was up to 35 °C more thermostable than the extant forms, with a ¹⁰T₅₀ (temperature at which 50% of the hemoprotein remains intact after 10 min) of 71 °C. This robustness may have facilitated evolutionary diversification of the CYP1s by buffering the destabilizing effects of mutations that conferred novel functions, a phenomenon which may also be useful in exploiting the catalytic versatility of these ancestral enzymes for commercial application as biocatalysts.

Comparisons of Ribonuclease HI Homologs and Mutants Uncover a Multistate Model for Substrate Recognition

Ribonuclease HI (RNHI) nonspecifically cleaves the RNA strand of RNA:DNA hybrid duplexes in a myriad of biological processes. Several RNHI homologs contain an extended domain, termed the handle region, which is critical to substrate binding. Nuclear magnetic resonance (NMR) spectroscopy and molecular dynamics (MD) simulations have suggested a kinetic model in which the handle region can exist in open (substrate-binding competent) or closed (substrate-binding incompetent) states in homologs containing arginine or lysine at position 88 (using sequence numbering of E. coli RNHI), while the handle region populates states intermediate between the open and closed conformers in homologs with asparagine at residue 88 [Stafford, K. A., et al., PLoS Comput. Biol. 2013, 9, 1-10]. NMR parameters characterizing handle region dynamics are highly correlated with enzymatic activity for RNHI homologs with two-state (open/closed) handle regions [Martin, J. A., et al., Biochemistry 2020, 59, 3201-3205]. The work presented herein shows that homologs containing asparagine 88 display distinct structural features compared with their counterparts containing arginine or lysine 88. Comparisons of RNHI homologs and site-directed mutants with asparagine 88 support a kinetic model for handle region dynamics that includes 12 unique transitions between eight conformations. Overall, these findings present an example of the structure-function relationships of enzymes and spotlight the use of NMR spectroscopy and MD simulations in uncovering fine details of conformational preferences.

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.

System drift and speciation

Even if a species' phenotype does not change over evolutionary time, the underlying mechanism may change, as distinct molecular pathways can realize identical phenotypes. Here we use linear system theory to explore the consequences of this idea, describing how a gene network underlying a conserved phenotype evolves, as the genetic drift of small changes to these molecular pathways causes a population to explore the set of mechanisms with identical phenotypes. To do this, we model an organism's internal state as a linear system of differential equations for which the environment provides input and the phenotype is the output, in which context there exists an exact characterization of the set of all mechanisms that give the same input-output relationship. This characterization implies that selectively neutral directions in genotype space should be common and that the evolutionary exploration of these distinct but equivalent mechanisms can lead to the reproductive incompatibility of independently evolving populations. This evolutionary exploration, or system drift, is expected to proceed at a rate proportional to the amount of intrapopulation genetic variation divided by the effective population size ( Ne$N_e$ ). At biologically reasonable parameter values this could lead to substantial interpopulation incompatibility, and thus speciation, on a time scale of Ne$N_e$ generations. This model also naturally predicts Haldane's rule, thus providing a concrete explanation of why heterogametic hybrids tend to be disrupted more often than homogametes during the early stages of speciation.

Systematic investigation of the link between enzyme catalysis and cold adaptation

Cold temperature is prevalent across the biosphere and slows the rates of chemical reactions. Increased catalysis has been predicted to be a dominant adaptive trait of enzymes to reduced temperature, and this expectation has informed physical models for enzyme catalysis and influenced bioprospecting strategies. To systematically test rate enhancement as an adaptive trait to cold, we paired kinetic constants of 2223 enzyme reactions with their organism's optimal growth temperature (T_Growth) and analyzed trends of rate constants as a function of T_Growth. These data do not support a general increase in rate enhancement in cold adaptation. In the model enzyme ketosteroid isomerase (KSI), there is prior evidence for temperature adaptation from a change in an active site residue that results in a tradeoff between activity and stability. Nevertheless, we found that little of the rate constant variation for 20 KSI variants was accounted for by T_Growth. In contrast, and consistent with prior expectations, we observed a correlation between stability and T_Growth across 433 proteins. These results suggest that temperature exerts a weaker selection pressure on enzyme rate constants than stability and that evolutionary forces other than temperature are responsible for the majority of enzymatic rate constant variation.

Epistasis shapes the fitness landscape of an allosteric specificity switch

Epistasis is a major determinant in the emergence of novel protein function. In allosteric proteins, direct interactions between inducer-binding mutations propagate through the allosteric network, manifesting as epistasis at the level of biological function. Elucidating this relationship between local interactions and their global effects is essential to understanding evolution of allosteric proteins. We integrate computational design, structural and biophysical analysis to characterize the emergence of novel inducer specificity in an allosteric transcription factor. Adaptive landscapes of different inducers of the designed mutant show that a few strong epistatic interactions constrain the number of viable sequence pathways, revealing ridges in the fitness landscape leading to new specificity. The structure of the designed mutant shows that a striking change in inducer orientation still retains allosteric function. Comparing biophysical and functional properties suggests a nonlinear relationship between inducer binding affinity and allostery. Our results highlight the functional and evolutionary complexity of allosteric proteins.

Epistasis plays an important role in the evolution of novel protein functions because it determines the mutational path a protein takes. Here, the authors combine functional, structural and biophysical analyses to characterize epistasis in a computationally redesigned ligand-inducible allosteric transcription factor and found that epistasis creates distinct biophysical and biological functional landscapes.

Thermostable lipases and their dynamics of improved enzymatic properties

Thermal stability is one of the most desirable characteristics in the search for novel lipases. The search for thermophilic microorganisms for synthesising functional enzyme biocatalysts with the ability to withstand high temperature, and capacity to maintain their native state in extreme conditions opens up new opportunities for their biotechnological applications. Thermophilic organisms are one of the most favoured organisms, whose distinctive characteristics are extremely related to their cellular constituent particularly biologically active proteins. Modifications on the enzyme structure are critical in optimizing the stability of enzyme to thermophilic conditions. Thermostable lipases are one of the most favourable enzymes used in food industries, pharmaceutical field, and actively been studied as potential biocatalyst in biodiesel production and other biotechnology application. Particularly, there is a trade-off between the use of enzymes in high concentration of organic solvents and product generation. Enhancement of the enzyme stability needs to be achieved for them to maintain their enzymatic activity regardless the environment. Various approaches on protein modification applied since decades ago conveyed a better understanding on how to improve the enzymatic properties in thermophilic bacteria. In fact, preliminary approach using advanced computational analysis is practically conducted before any modification is being performed experimentally. Apart from that, isolation of novel extremozymes from various microorganisms are offering great frontier in explaining the crucial native interaction within the molecules which could help in protein engineering. In this review, the thermostability prospect of lipases and the utility of protein engineering insights into achieving functional industrial usefulness at their high temperature habitat are highlighted. Similarly, the underlying thermodynamic and structural basis that defines the forces that stabilize these thermostable lipase is discussed. KEY POINTS: • The dynamics of lipases contributes to their non-covalent interactions and structural stability. • Thermostability can be enhanced by well-established genetic tools for improved kinetic efficiency. • Molecular dynamics greatly provides structure-function insights on thermodynamics of lipase.

The Stability Landscape of de novo TIM Barrels Explored by a Modular Design Approach

•

The TIM barrel is a versatile fold to understand structure-stability relationships.

•

A collection of de novo TIM barrels with improved hydrophobic cores was designed.

•

DeNovoTIMs are reversible in chemical and thermal unfolding, which is uncommon in TIM barrels.

•

Epistatic effects play a central role in DeNovoTIMs stabilization.

•

DeNovoTIMs navigate a previously uncharted region of the stability landscape.

The ability to design stable proteins with custom-made functions is a major goal in biochemistry with practical relevance for our environment and society. Understanding and manipulating protein stability provide crucial information on the molecular determinants that modulate structure and stability, and expand the applications of de novo proteins. Since the (β/⍺)₈-barrel or TIM-barrel fold is one of the most common functional scaffolds, in this work we designed a collection of stable de novo TIM barrels (DeNovoTIMs), using a computational fixed-backbone and modular approach based on improved hydrophobic packing of sTIM11, the first validated de novo TIM barrel, and subjected them to a thorough folding analysis. DeNovoTIMs navigate a region of the stability landscape previously uncharted by natural TIM barrels, with variations spanning 60 degrees in melting temperature and 22 kcal per mol in conformational stability throughout the designs. Significant non-additive or epistatic effects were observed when stabilizing mutations from different regions of the barrel were combined. The molecular basis of epistasis in DeNovoTIMs appears to be related to the extension of the hydrophobic cores. This study is an important step towards the fine-tuned modulation of protein stability by design.

A strategy for proline and glycine mutations to proteins with alchemical free energy calculations

Computation of the thermodynamic consequences of protein mutations holds great promise in protein biophysics and design. Alchemical free energy methods can give improved estimates of mutational free energies, and are already widely used in calculations of relative and absolute binding free energies in small molecule design problems. In principle, alchemical methods can address any amino acid mutation with an appropriate alchemical pathway, but identifying a strategy that produces such a path for proline and glycine mutations is an ongoing challenge. Most current strategies perturb only side chain atoms, while proline and glycine mutations also alter the backbone parameters and backbone ring topology. Some strategies also perturb backbone parameters and enable glycine mutations. This work presents a strategy that enables both proline and glycine mutations and comprises two key elements: a dual backbone with restraints and scaling of bonded terms, facilitating backbone parameter changes, and a soft bond in the proline ring, enabling ring topology changes in proline mutations. These elements also have utility for core hopping and macrocycle studies in computer-aided drug design. This new strategy shows slight improvements over an alternative side chain perturbation strategy for a set T4 lysozyme mutations lacking proline and glycine, and yields good agreement with experiment for a set of T4 lysozyme proline and glycine mutations not previously studied. To our knowledge this is the first report comparing alchemical predictions of proline mutations with experiment. With this strategy in hand, alchemical methods now have access to the full palette of amino acid mutations.

Evolution, folding, and design of TIM barrels and related proteins

Proteins are chief actors in life that perform a myriad of exquisite functions. This diversity has been enabled through the evolution and diversification of protein folds. Analysis of sequences and structures strongly suggest that numerous protein pieces have been reused as building blocks and propagated to many modern folds. This information can be traced to understand how the protein world has diversified. In this review, we discuss the latest advances in the analysis of protein evolutionary units, and we use as a model system one of the most abundant and versatile topologies, the TIM-barrel fold, to highlight the existing common principles that interconnect protein evolution, structure, folding, function, and design.

Resurrection of Ancestral Malate Dehydrogenases Reveals the Evolutionary History of Halobacterial Proteins : Deciphering Gene Trajectories and Changes in Biochemical Properties

Extreme halophilic Archaea thrive in high salt, where, through proteomic adaptation, they cope with the strong osmolarity and extreme ionic conditions of their environment. In spite of wide fundamental interest, however, studies providing insights into this adaptation are scarce, because of practical difficulties inherent to the purification and characterization of halophilic enzymes. In this work, we describe the evolutionary history of malate dehydrogenases (MalDH) within Halobacteria (a class of the Euryarchaeota phylum). We resurrected nine ancestors along the inferred halobacterial MalDH phylogeny, including the Last Common Ancestral MalDH of Halobacteria (LCAHa) and compared their biochemical properties with those of five modern halobacterial MalDHs. We monitored the stability of these various MalDHs, their oligomeric states and enzymatic properties, as a function of concentration for different salts in the solvent. We found that a variety of evolutionary processes such as amino acid replacement, gene duplication, loss of MalDH gene and replacement owing to horizontal transfer resulted in significant differences in solubility, stability and catalytic properties between these enzymes in the three Halobacteriales, Haloferacales and Natrialbales orders since the LCAHa MalDH.We also showed how a stability trade-off might favor the emergence of new properties during adaptation to diverse environmental conditions. Altogether, our results suggest a new view of halophilic protein adaptation in Archaea.

Ancestral sequence reconstruction - An underused approach to understand the evolution of gene function in plants?

Whilst substantial research effort has been placed on understanding the interactions of plant proteins with their molecular partners, relatively few studies in plants - by contrast to work in other organisms - address how these interactions evolve. It is thought that ancestral proteins were more promiscuous than modern proteins and that specificity often evolved following gene duplication and subsequent functional refining. However, ancestral protein resurrection studies have found that some modern proteins have evolved de novo from ancestors lacking those functions. Intriguingly, the new interactions evolved as a consequence of just a few mutations and, as such, acquisition of new functions appears to be neither difficult nor rare, however, only a few of them are incorporated into biological processes before they are lost to subsequent mutations. Here, we detail the approach of ancestral sequence reconstruction (ASR), providing a primer to reconstruct the sequence of an ancestral gene. We will present case studies from a range of different eukaryotes before discussing the few instances where ancestral reconstructions have been used in plants. As ASR is used to dig into the remote evolutionary past, we will also present some alternative genetic approaches to investigate molecular evolution on shorter timescales. We argue that the study of plant secondary metabolism is particularly well suited for ancestral reconstruction studies. Indeed, its ancient evolutionary roots and highly diverse landscape provide an ideal context in which to address the focal issue around the emergence of evolutionary novelties and how this affects the chemical diversification of plant metabolism.

Evolution of frustrated and stabilising contacts in reconstructed ancient proteins

Energetic properties of a protein are a major determinant of its evolutionary fitness. Using a reconstruction algorithm, dating the reconstructed proteins and calculating the interaction network between their amino acids through a coevolutionary approach, we studied how the interactions that stabilise 890 proteins, belonging to five families, evolved for billions of years. In particular, we focused our attention on the network of most strongly attractive contacts and on that of poorly optimised, frustrated contacts. Our results support the idea that the cluster of most attractive interactions extends its size along evolutionary time, but from the data, we cannot conclude that protein stability or that the degree of frustration tends always to decrease.

Ancestral sequence reconstruction produces thermally stable enzymes with mesophilic enzyme-like catalytic properties

Enzymes have high catalytic efficiency and low environmental impact, and are therefore potentially useful tools for various industrial processes. Crucially, however, natural enzymes do not always have the properties required for specific processes. It may be necessary, therefore, to design, engineer, and evolve enzymes with properties that are not found in natural enzymes. In particular, the creation of enzymes that are thermally stable and catalytically active at low temperature is desirable for processes involving both high and low temperatures. In the current study, we designed two ancestral sequences of 3-isopropylmalate dehydrogenase by an ancestral sequence reconstruction technique based on a phylogenetic analysis of extant homologous amino acid sequences. Genes encoding the designed sequences were artificially synthesized and expressed in Escherichia coli. The reconstructed enzymes were found to be slightly more thermally stable than the extant thermophilic homologue from Thermus thermophilus. Moreover, they had considerably higher low-temperature catalytic activity as compared with the T. thermophilus enzyme. Detailed analyses of their temperature-dependent specific activities and kinetic properties showed that the reconstructed enzymes have catalytic properties similar to those of mesophilic homologues. Collectively, our study demonstrates that ancestral sequence reconstruction can produce a thermally stable enzyme with catalytic properties adapted to low-temperature reactions.

Larger active site in an ancestral hydroxynitrile lyase increases catalytically promiscuous esterase activity

Hydroxynitrile lyases (HNL's) belonging to the α/β-hydrolase-fold superfamily evolved from esterases approximately 100 million years ago. Reconstruction of an ancestral hydroxynitrile lyase in the α/β-hydrolase fold superfamily yielded a catalytically active hydroxynitrile lyase, HNL1. Several properties of HNL1 differ from the modern HNL from rubber tree (HbHNL). HNL1 favors larger substrates as compared to HbHNL, is two-fold more catalytically promiscuous for ester hydrolysis (p-nitrophenyl acetate) as compared to mandelonitrile cleavage, and resists irreversible heat inactivation to 35 °C higher than for HbHNL. We hypothesized that the x-ray crystal structure of HNL1 may reveal the molecular basis for the differences in these properties. The x-ray crystal structure solved to 1.96-Å resolution shows the expected α/β-hydrolase fold, but a 60% larger active site as compared to HbHNL. This larger active site echoes its evolution from esterases since related esterase SABP2 from tobacco also has a 38% larger active site than HbHNL. The larger active site in HNL1 likely accounts for its ability to accept larger hydroxynitrile substrates. Site-directed mutagenesis of HbHNL to expand the active site increased its promiscuous esterase activity 50-fold, consistent with the larger active site in HNL1 being the primary cause of its promiscuous esterase activity. Urea-induced unfolding of HNL1 indicates that it unfolds less completely than HbHNL (m-value = 0.63 for HNL1 vs 0.93 kcal/mol·M for HbHNL), which may account for the ability of HNL1 to better resist irreversible inactivation upon heating. The structure of HNL1 shows changes in hydrogen bond networks that may stabilize regions of the folded structure.

Interplay between Conformational Strain and Intramolecular Interaction in Protein Structures: Which of Them Is Evolutionarily Conserved?

By computing strain energies of peptide fragments within protein structures and their intramolecular interaction energies, we attempt to reveal general biophysical trends behind the secondary structure formation in the context of protein evolution. Our "protein basis set" consisted of 1143 representatives of different folds obtained from curated SCOPe database, and for each member of the set, the strain and intramolecular energy was calculated on the "rolling tripeptide" basis, employing the DFT-D3/COSMO-RS method for the former and the QM-calibrated force field method (MM) for the latter. The calculated data, strain and interactions, were correlated with the conservation of amino acid residues in secondary structure elements and also with the level of the residue burial within the protein three-dimensional structure. It allowed us to formulate several observations concerning fundamental differences between two main secondary structure motifs: α-helices and β-strands. We have shown that a strong interaction is one of the determining characteristics of the β-sheet formation, at least at the level of tripeptides (and likely penta- or heptapeptides, too), and that the β-strand is a prevailing secondary structure in the strongly-interacting regions of the protein folds conserved by evolution. On the other hand, low strain was neither proven to be an important physicochemical property conserved by evolution nor does it correlate with the propensity for the α-helix and β-strand. Finally, it has been demonstrated that the strong interaction has a certain level of connection with residue burial; however, we demonstrate that these two characteristics should be rather regarded as two complementary factors. These findings represent an important contribution to understanding protein folding from first principles, which is a complementary approach to ongoing efforts to solve the protein folding problem by knowledge-based approaches and machine-learning.

Larger active site in an ancestral hydroxynitrile lyase increases catalytically promiscuous esterase activity.. [DOI: 10.1101/2020.04.06.027797] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

Hydroxynitrile lyases (HNL’s) belonging to the α/β-hydrolase-fold superfamily evolved from esterases approximately 100 million years ago. Reconstruction of an ancestral hydroxynitrile lyase in the α/β-hydrolase fold superfamily yielded a catalytically active hydroxynitrile lyase, HNL1. Several properties of HNL1 differ from the modern HNL from rubber tree (HbHNL). HNL1 favors larger substrates as compared to HbHNL, is two-fold more catalytically promiscuous for ester hydrolysis (p-nitrophenyl acetate) as compared to mandelonitrile cleavage, and resists irreversible heat inactivation to 35 °C higher than for HbHNL. We hypothesized that the x-ray crystal structure of HNL1 may reveal the molecular basis for the differences in these properties. The x-ray crystal structure solved to 1.96-Å resolution shows the expected α/β-hydrolase fold, but a 60% larger active site as compared to HbHNL. This larger active site echoes its evolution from esterases since related esterase SABP2 from tobacco also has a 38% larger active site than HbHNL. The larger active site in HNL1 likely accounts for its ability to accept larger hydroxynitrile substrates. Site-directed mutagenesis of HbHNL to expand the active site increased its promiscuous esterase activity 50-fold, consistent with the larger active site in HNL1 being the primary cause of its promiscuous esterase activity. Urea-induced unfolding of HNL1 indicates that it unfolds less completely than HbHNL (m-value = 0.63 for HNL1 vs 0.93 kcal/ mol·M for HbHNL), which may account for the ability of HNL1 to better resist irreversible inactivation upon heating. The structure of HNL1 shows changes in hydrogen bond networks that may stabilize regions of the folded structure.

Inferring Biochemical Reactions and Metabolite Structures to Understand Metabolic Pathway Drift

Inferring genome-scale metabolic networks in emerging model organisms is challenged by incomplete biochemical knowledge and partial conservation of biochemical pathways during evolution. Therefore, specific bioinformatic tools are necessary to infer biochemical reactions and metabolic structures that can be checked experimentally. Using an integrative approach combining genomic and metabolomic data in the red algal model Chondrus crispus, we show that, even metabolic pathways considered as conserved, like sterols or mycosporine-like amino acid synthesis pathways, undergo substantial turnover. This phenomenon, here formally defined as "metabolic pathway drift," is consistent with findings from other areas of evolutionary biology, indicating that a given phenotype can be conserved even if the underlying molecular mechanisms are changing. We present a proof of concept with a methodological approach to formalize the logical reasoning necessary to infer reactions and molecular structures, abstracting molecular transformations based on previous biochemical knowledge.

Conformational Flexibility Drives Cold Adaptation in Pseudoalteromonas haloplanktis TAC125 Globins

Significance: Temperature is one of the most important drivers in shaping protein adaptations. Many biochemical and physiological processes are influenced by temperature. Proteins and enzymes from organisms living at low temperature are less stable in comparison to high-temperature adapted proteins. The lower stability is generally due to greater conformational flexibility. Recent Advances: Adaptive changes in the structure of cold-adapted proteins may occur at subunit interfaces, distant from the active site, thus producing energy changes associated with conformational transitions transmitted to the active site by allosteric modulation, valid also for monomeric proteins in which tertiary structural changes may play an essential role. Critical Issues: Despite efforts, the current experimental and computational methods still fail to produce general principles on protein evolution, since many changes are protein and species dependent. Environmental constraints or other biological cellular signals may override the ancestral information included in the structure of the protein, thus introducing inaccuracy in estimates and predictions on the evolutionary adaptations of proteins in response to cold adaptation. Future Directions: In this review, we describe the studies and approaches used to investigate stability and flexibility in the cold-adapted globins of the Antarctic marine bacterium Pseudoalteromonas haloplanktis TAC125. In fact, future research directions will be prescient on more detailed investigation of cold-adapted proteins and the role of fluctuations between different conformational states.

Modeling Structural Constraints on Protein Evolution via Side-Chain Conformational States

Few models of sequence evolution incorporate parameters describing protein structure, despite its high conservation, essential functional role and increasing availability. We present a structurally aware empirical substitution model for amino acid sequence evolution in which proteins are expressed using an expanded alphabet that relays both amino acid identity and structural information. Each character specifies an amino acid as well as information about the rotamer configuration of its side-chain: the discrete geometric pattern of permitted side-chain atomic positions, as defined by the dihedral angles between covalently linked atoms. By assigning rotamer states in 251,194 protein structures and identifying 4,508,390 substitutions between closely related sequences, we generate a 55-state “Dayhoff-like” model that shows that the evolutionary properties of amino acids depend strongly upon side-chain geometry. The model performs as well as or better than traditional 20-state models for divergence time estimation, tree inference, and ancestral state reconstruction. We conclude that not only is rotamer configuration a valuable source of information for phylogenetic studies, but that modeling the concomitant evolution of sequence and structure may have important implications for understanding protein folding and function.

Deciphering protein evolution and fitness landscapes with latent space models

Protein sequences contain rich information about protein evolution, fitness landscapes, and stability. Here we investigate how latent space models trained using variational auto-encoders can infer these properties from sequences. Using both simulated and real sequences, we show that the low dimensional latent space representation of sequences, calculated using the encoder model, captures both evolutionary and ancestral relationships between sequences. Together with experimental fitness data and Gaussian process regression, the latent space representation also enables learning the protein fitness landscape in a continuous low dimensional space. Moreover, the model is also useful in predicting protein mutational stability landscapes and quantifying the importance of stability in shaping protein evolution. Overall, we illustrate that the latent space models learned using variational auto-encoders provide a mechanism for exploration of the rich data contained in protein sequences regarding evolution, fitness and stability and hence are well-suited to help guide protein engineering efforts.

How to resurrect ancestral proteins as proxies for ancient biogeochemistry

Throughout the history of life, enzymes have served as the primary molecular mediators of biogeochemical cycles by catalyzing the metabolic pathways that interact with geochemical substrates. The byproducts of enzymatic activities have been preserved as chemical and isotopic signatures in the geologic record. However, interpretations of these signatures are limited by the assumption that such enzymes have remained functionally conserved over billions of years of molecular evolution. By reconstructing ancient genetic sequences in conjunction with laboratory enzyme resurrection, preserved biogeochemical signatures can instead be related to experimentally constrained, ancestral enzymatic properties. We may thereby investigate instances within molecular evolutionary trajectories potentially tied to significant biogeochemical transitions evidenced in the geologic record. Here, we survey recent enzyme resurrection studies to provide a reasoned assessment of areas of success and common pitfalls relevant to ancient biogeochemical applications. We conclude by considering the Great Oxidation Event, which provides a constructive example of a significant biogeochemical transition that warrants investigation with ancestral enzyme resurrection. This event also serves to highlight the pitfalls of facile interpretation of paleophenotype models and data, as applied to two examples of enzymes that likely both influenced and were influenced by the rise of atmospheric oxygen - RuBisCO and nitrogenase.

Evolutionary Modes in Protein Observable Space: The Case of Thioredoxins

In this article, we investigated the structural and dynamical evolutionary behaviour of a set of ten thioredoxin proteins as formed by three extant forms and seven resurrected ones in laboratory. Starting from the crystallographic structures, we performed all-atom molecular dynamics simulations and compare the trajectories in terms of structural and dynamical properties. Interestingly, the structural properties related to the protein density (i.e. the number of residues divided by the excluded molecular volume) well describe the protein evolutionary behaviour. Our results also suggest that the changes in sequence as occurred during the evolution have affected the protein essential motions, allowing us to discriminate between ancient and extant proteins in terms of their dynamical behaviour. Such results are yet more evident when the bacterial, archaeal and eukaryotic thioredoxins are separately analysed.

Ancient mechanisms for the evolution of the bicoid homeodomain's function in fly development

The ancient mechanisms that caused developmental gene regulatory networks to diversify among distantly related taxa are not well understood. Here we use ancestral protein reconstruction, biochemical experiments, and developmental assays of transgenic animals carrying reconstructed ancestral genes to investigate how the transcription factor Bicoid (Bcd) evolved its central role in anterior-posterior patterning in flies. We show that most of Bcd's derived functions are attributable to evolutionary changes within its homeodomain (HD) during a phylogenetic interval >140 million years ago. A single substitution from this period (Q50K) accounts almost entirely for the evolution of Bcd's derived DNA specificity in vitro. In transgenic embryos expressing the reconstructed ancestral HD, however, Q50K confers activation of only a few of Bcd's transcriptional targets and yields a very partial rescue of anterior development. Adding a second historical substitution (M54R) confers regulation of additional Bcd targets and further rescues anterior development. These results indicate that two epistatically interacting mutations played a major role in the evolution of Bcd's controlling regulatory role in early development. They also show how ancestral sequence reconstruction can be combined with in vivo characterization of transgenic animals to illuminate the historical mechanisms of developmental evolution.

Tracing a protein's folding pathway over evolutionary time using ancestral sequence reconstruction and hydrogen exchange

The conformations populated during protein folding have been studied for decades; yet, their evolutionary importance remains largely unexplored. Ancestral sequence reconstruction allows access to proteins across evolutionary time, and new methods such as pulsed-labeling hydrogen exchange coupled with mass spectrometry allow determination of folding intermediate structures at near amino-acid resolution. Here, we combine these techniques to monitor the folding of the ribonuclease H family along the evolutionary lineages of T. thermophilus and E. coli RNase H. All homologs and ancestral proteins studied populate a similar folding intermediate despite being separated by billions of years of evolution. Even though this conformation is conserved, the pathway leading to it has diverged over evolutionary time, and rational mutations can alter this trajectory. Our results demonstrate that evolutionary processes can affect the energy landscape to preserve or alter specific features of a protein’s folding pathway.

The burst-phase folding intermediate of ribonuclease H changes conformation over evolutionary history

The amino acid sequence encodes the energy landscape of a protein. Therefore, we expect evolutionary mutations to change features of the protein energy landscape, including the conformations adopted by a polypeptide as it folds to its native state. Ribonucleases H (RNase H) from Escherichia coli and Thermus thermophilus both fold via a partially folded intermediate in which the core region of the protein (helices A-D and strands 4-5) is structured. Strand 1, however, uniquely contributes to the T. thermophilus RNase H folding intermediate (Icore+1 ), but not the E. coli RNase H intermediate (Icore ) (Rosen & Marqusee, PLoS One 2015). We explore the origin of this difference by characterizing the folding intermediate of seven ancestral RNases H spanning the evolutionary history of these two homologs. Using fragment models with or without strand 1 and FRET probes to characterize the folding intermediate of each ancestor, we find a distinct evolutionary trend across the family-the involvement of strand 1 in the folding intermediate is an ancestral feature that is maintained in the thermophilic lineage and is gradually lost in the mesophilic lineage. Evolutionary sequence changes indeed modulate the conformations present on the folding landscape and altered the folding trajectory of RNase H.

Evolutionary Divergent Suppressor Mutations in Conformational Diseases

Neutral and adaptive mutations are key players in the evolutionary dynamics of proteins at molecular, cellular and organismal levels. Conversely, largely destabilizing mutations are rarely tolerated by evolution, although their occurrence in diverse human populations has important roles in the pathogenesis of conformational diseases. We have recently proposed that divergence at certain sites from the consensus (amino acid) state during mammalian evolution may have rendered some human proteins more vulnerable towards disease-associated mutations, primarily by decreasing their conformational stability. We herein extend and refine this hypothesis discussing results from phylogenetic and structural analyses, structure-based energy calculations and structure-function studies at molecular and cellular levels. As proof-of-principle, we focus on different mammalian orthologues of the NQO1 (NAD(P)H:quinone oxidoreductase 1) and AGT (alanine:glyoxylate aminotransferase) proteins. We discuss the different loss-of-function pathogenic mechanisms associated with diseases involving the two enzymes, including enzyme inactivation, accelerated degradation, intracellular mistargeting, and aggregation. Last, we take into account the potentially higher robustness of mammalian orthologues containing certain consensus amino acids as suppressors of human disease, and their relation with different intracellular post-translational modifications and protein quality control capacities, to be discussed as sources of phenotypic variability between human and mammalian models of disease and as tools for improving current therapeutic approaches.

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.

Proteome Evolution of Deep-Sea Hydrothermal Vent Alvinellid Polychaetes Supports the Ancestry of Thermophily and Subsequent Adaptation to Cold in Some Lineages

Temperature, perhaps more than any other environmental factor, is likely to influence the evolution of all organisms. It is also a very interesting factor to understand how genomes are shaped by selection over evolutionary timescales, as it potentially affects the whole genome. Among thermophilic prokaryotes, temperature affects both codon usage and protein composition to increase the stability of the transcriptional/translational machinery, and the resulting proteins need to be functional at high temperatures. Among eukaryotes less is known about genome evolution, and the tube-dwelling worms of the family Alvinellidae represent an excellent opportunity to test hypotheses about the emergence of thermophily in ectothermic metazoans. The Alvinellidae are a group of worms that experience varying thermal regimes, presumably having evolved into these niches over evolutionary times. Here we analyzed 423 putative orthologous loci derived from 6 alvinellid species including the thermophilic Alvinella pompejana and Paralvinella sulfincola. This comparative approach allowed us to assess amino acid composition, codon usage, divergence, direction of residue changes and the strength of selection along the alvinellid phylogeny, and to design a new eukaryotic thermophilic criterion based on significant differences in the residue composition of proteins. Contrary to expectations, the alvinellid ancestor of all present-day species seems to have been thermophilic, a trait subsequently maintained by purifying selection in lineages that still inhabit higher temperature environments. In contrast, lineages currently living in colder habitats likely evolved under selective relaxation, with some degree of positive selection for low-temperature adaptation at the protein level.

Characterization of Reconstructed Ancestral Proteins Suggests a Change in Temperature of the Ancient Biosphere

Understanding the evolution of ancestral life, and especially the ability of some organisms to flourish in the variable environments experienced in Earth’s early biosphere, requires knowledge of the characteristics and the environment of these ancestral organisms. Information about early life and environmental conditions has been obtained from fossil records and geological surveys. Recent advances in phylogenetic analysis, and an increasing number of protein sequences available in public databases, have made it possible to infer ancestral protein sequences possessed by ancient organisms. However, the in silico studies that assess the ancestral base content of ribosomal RNAs, the frequency of each amino acid in ancestral proteins, and estimate the environmental temperatures of ancient organisms, show conflicting results. The characterization of ancestral proteins reconstructed in vitro suggests that ancient organisms had very thermally stable proteins, and therefore were thermophilic or hyperthermophilic. Experimental data supports the idea that only thermophilic ancestors survived the catastrophic increase in temperature of the biosphere that was likely associated with meteorite impacts during the early history of Earth. In addition, by expanding the timescale and including more ancestral proteins for reconstruction, it appears as though the Earth’s surface temperature gradually decreased over time, from Archean to present.

Exploring the past and the future of protein evolution with ancestral sequence reconstruction: the 'retro' approach to protein engineering

A central goal in molecular evolution is to understand the ways in which genes and proteins evolve in response to changing environments. In the absence of intact DNA from fossils, ancestral sequence reconstruction (ASR) can be used to infer the evolutionary precursors of extant proteins. To date, ancestral proteins belonging to eubacteria, archaea, yeast and vertebrates have been inferred that have been hypothesized to date from between several million to over 3 billion years ago. ASR has yielded insights into the early history of life on Earth and the evolution of proteins and macromolecular complexes. Recently, however, ASR has developed from a tool for testing hypotheses about protein evolution to a useful means for designing novel proteins. The strength of this approach lies in the ability to infer ancestral sequences encoding proteins that have desirable properties compared with contemporary forms, particularly thermostability and broad substrate range, making them good starting points for laboratory evolution. Developments in technologies for DNA sequencing and synthesis and computational phylogenetic analysis have led to an escalation in the number of ancient proteins resurrected in the last decade and greatly facilitated the use of ASR in the burgeoning field of synthetic biology. However, the primary challenge of ASR remains in accurately inferring ancestral states, despite the uncertainty arising from evolutionary models, incomplete sequences and limited phylogenetic trees. This review will focus, firstly, on the use of ASR to uncover links between sequence and phenotype and, secondly, on the practical application of ASR in protein engineering.