Ancestral Sequence Reconstruction: From Chemical Paleogenetics to Maximum Likelihood Algorithms and Beyond

The Characterization of Ancient Methanococcales Malate Dehydrogenases Reveals That Strong Thermal Stability Prevents Unfolding Under Intense γ-Irradiation

Malate dehydrogenases (MalDHs) (EC.1.1.1.37), which are involved in the conversion of oxaloacetate to pyruvate in the tricarboxylic acid cycle, are a relevant model for the study of enzyme evolution and adaptation. Likewise, a recent study showed that Methanococcales, a major lineage of Archaea, is a good model to study the molecular processes of proteome thermoadaptation in prokaryotes. Here, we use ancestral sequence reconstruction and paleoenzymology to characterize both ancient and extant MalDHs. We observe a good correlation between inferred optimal growth temperatures and experimental optimal temperatures for activity (A-Topt). In particular, we show that the MalDH present in the ancestor of Methanococcales was hyperthermostable and had an A-Topt of 80 °C, consistent with a hyperthermophilic lifestyle. This ancestor gave rise to two lineages with different thermal constraints: one remained hyperthermophilic, while the other underwent several independent adaptations to colder environments. Surprisingly, the enzymes of the first lineage have retained a thermoresistant behavior (i.e. strong thermostability and high A-Topt), whereas the ancestor of the second lineage shows a strong thermostability, but a reduced A-Topt. Using mutants, we mimic the adaptation trajectory toward mesophily and show that it is possible to significantly reduce the A-Topt without altering the thermostability of the enzyme by introducing a few mutations. Finally, we reveal an unexpected link between thermostability and the ability to resist γ-irradiation-induced unfolding.

Applications of ancestral sequence reconstruction for understanding the evolution of plant specialized metabolism

Studies of enzymes in modern-day plants have documented the diversity of metabolic activities retained by species today but only provide limited insight into how those properties evolved. Ancestral sequence reconstruction (ASR) is an approach that provides statistical estimates of ancient plant enzyme sequences which can then be resurrected to test hypotheses about the evolution of catalytic activities and pathway assembly. Here, I review the insights that have been obtained using ASR to study plant metabolism and highlight important methodological aspects. Overall, studies of resurrected plant enzymes show that (i) exaptation is widespread such that even low or undetectable levels of ancestral activity with a substrate can later become the apparent primary activity of descendant enzymes, (ii) intramolecular epistasis may or may not limit evolutionary paths towards catalytic or substrate preference switches, and (iii) ancient pathway flux often differs from modern-day metabolic networks. These and other insights gained from ASR would not have been possible using only modern-day sequences. Future ASR studies characterizing entire ancestral metabolic networks as well as those that link ancient structures with enzymatic properties should continue to provide novel insights into how the chemical diversity of plants evolved. This article is part of the theme issue 'The evolution of plant metabolism'.

Insertions and Deletions: Computational Methods, Evolutionary Dynamics, and Biological Applications

Insertions and deletions constitute the second most important source of natural genomic variation. Insertions and deletions make up to 25% of genomic variants in humans and are involved in complex evolutionary processes including genomic rearrangements, adaptation, and speciation. Recent advances in long-read sequencing technologies allow detailed inference of insertions and deletion variation in species and populations. Yet, despite their importance, evolutionary studies have traditionally ignored or mishandled insertions and deletions due to a lack of comprehensive methodologies and statistical models of insertions and deletion dynamics. Here, we discuss methods for describing insertions and deletion variation and modeling insertions and deletions over evolutionary time. We provide practical advice for tackling insertions and deletions in genomic sequences and illustrate our discussion with examples of insertions and deletion-induced effects in human and other natural populations and their contribution to evolutionary processes. We outline promising directions for future developments in statistical methodologies that would allow researchers to analyze insertions and deletion variation and their effects in large genomic data sets and to incorporate insertions and deletions in evolutionary inference.

Unlocking the potential of enzyme engineering via rational computational design strategies

Enzymes play a pivotal role in various industries by enabling efficient, eco-friendly, and sustainable chemical processes. However, the low turnover rates and poor substrate selectivity of enzymes limit their large-scale applications. Rational computational enzyme design, facilitated by computational algorithms, offers a more targeted and less labor-intensive approach. There has been notable advancement in employing rational computational protein engineering strategies to overcome these issues, it has not been comprehensively reviewed so far. This article reviews recent developments in rational computational enzyme design, categorizing them into three types: structure-based, sequence-based, and data-driven machine learning computational design. Case studies are presented to demonstrate successful enhancements in catalytic activity, stability, and substrate selectivity. Lastly, the article provides a thorough analysis of these approaches, highlights existing challenges and potential solutions, and offers insights into future development directions.

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.

The evolution of mammalian Rem2: unraveling the impact of purifying selection and coevolution on protein function, and implications for human disorders

Rad And Gem-Like GTP-Binding Protein 2 (Rem2), a member of the RGK family of Ras-like GTPases, is implicated in Huntington's disease and Long QT Syndrome and is highly expressed in the brain and endocrine cells. We examine the evolutionary history of Rem2 identified in various mammalian species, focusing on the role of purifying selection and coevolution in shaping its sequence and protein structural constraints. Our analysis of Rem2 sequences across 175 mammalian species found evidence for strong purifying selection in 70% of non-invariant codon sites which is characteristic of essential proteins that play critical roles in biological processes and is consistent with Rem2's role in the regulation of neuronal development and function. We inferred epistatic effects in 50 pairs of codon sites in Rem2, some of which are predicted to have deleterious effects on human health. Additionally, we reconstructed the ancestral evolutionary history of mammalian Rem2 using protein structure prediction of extinct and extant sequences which revealed the dynamics of how substitutions that change the gene sequence of Rem2 can impact protein structure in variable regions while maintaining core functional mechanisms. By understanding the selective pressures, protein- and gene - interactions that have shaped the sequence and structure of the Rem2 protein, we gain a stronger understanding of its biological and functional constraints.

Design, structure and plasma binding of ancestral β-CoV scaffold antigens

We report the application of ancestral sequence reconstruction on coronavirus spike protein, resulting in stable and highly soluble ancestral scaffold antigens (AnSAs). The AnSAs interact with plasma of patients recovered from COVID-19 but do not bind to the human angiotensin-converting enzyme 2 (ACE2) receptor. Cryo-EM analysis of the AnSAs yield high resolution structures (2.6-2.8 Å) indicating a closed pre-fusion conformation in which all three receptor-binding domains (RBDs) are facing downwards. The structures reveal an intricate hydrogen-bonding network mediated by well-resolved loops, both within and across monomers, tethering the N-terminal domain and RBD together. We show that AnSA-5 can induce and boost a broad-spectrum immune response against the wild-type RBD as well as circulating variants of concern in an immune organoid model derived from tonsils. Finally, we highlight how AnSAs are potent scaffolds by replacing the ancestral RBD with the wild-type sequence, which restores ACE2 binding and increases the interaction with convalescent plasma.

Structural and enzymatic characterization of Bacillus subtilis R,R-2,3-butanediol dehydrogenase

2,3-butanediol dehydrogenase (BDH, EC 1.1.1.76) also known as acetoin reductase (AR, EC 1.1.1.4) is the key enzyme converting acetoin (AC) into 2,3-butanediol (BD) and undertaking the irreversible conversion of diacetyl to acetoin in various microorganisms. The existence of three BDHs (R,R-, meso-, and S,S-BDH) product different BD isomers. Catalyzing mechanisms of meso- and S,S-BDH have been understood with the assistance of their X-ray crystal structures. However, the lack of structural data for R,R-BDH restricts the integral understanding of the catalytic mechanism of BDHs. In this study, we successfully crystallized and solved the X-ray crystal structure of Bacillus subtilis R,R-BDH. A zinc ion was found locating in the catalytic center and coordinated by Cys37, His70 and Glu152, helping to stabilize the chiral substrates observed in the predicted molecular docking model. The interaction patterns of different chiral substrates in the molecular docking model explained the react priority measured by the enzyme activity assay of R,R-BDH. Site-directed mutation experiments determined that the amino acids Cys37, Thr244, Ile268 and Lys340 are important in the catalytically active center. The structural information of R,R-BDH presented in this study accomplished the understanding of BDHs catalytic mechanism and more importantly provides useful guidance for the directional engineering of R,R-BDH to obtain high-purity monochiral BD and AC.

Structural evolution of an amphibian-specific globin: A computational evolutionary biochemistry approach

Studies on the globin family are continuously revealing insights into the mechanisms of gene and protein evolution. The rise of a new globin gene type in Pelobatoidea and Neobatrachia (Amphibia:Anura) from an α-globin precursor provides the opportunity to investigate the genetic and physical mechanisms underlying the origin of new protein structural and functional properties. This amphibian-specific globin (globin A/GbA) discovered in the heart of Rana catesbeiana is a monomer. As the ancestral oligomeric state of α-globins is a homodimer, we inferred that the ancestral state was lost somewhere in the GbA lineage. Here, we combined computational molecular evolution with structural bioinformatics to determine the extent to which the loss of the homodimeric state is pervasive in the GbA clade. We also characterized the loci of GbA genes in Bufo bufo. We found two GbA clades in Neobatrachia. One was deleted in Ranidae, but retained and expanded to yield a new globin cluster in Bufonidae species. Loss of the ancestral oligomeric state seems to be pervasive in the GbA clade. However, a taxonomic sampling that includes more Pelobatoidea, as well as early Neobatrachia, lineages would be necessary to determine the oligomeric state of the last common ancestor of all GbA. The evidence presented here points out a possible loss of oligomerization in Pelobatoidea GbA as a result of amino acid substitutions that weaken the homodimeric state. In contrast, the loss of oligomerization in both Neobatrachia GbA clades was linked to independent deletions that disrupted many packing contacts at the homodimer interface.

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

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

Engineering indel and substitution variants of diverse and ancient enzymes using Graphical Representation of Ancestral Sequence Predictions (GRASP)

Ancestral sequence reconstruction is a technique that is gaining widespread use in molecular evolution studies and protein engineering. Accurate reconstruction requires the ability to handle appropriately large numbers of sequences, as well as insertion and deletion (indel) events, but available approaches exhibit limitations. To address these limitations, we developed Graphical Representation of Ancestral Sequence Predictions (GRASP), which efficiently implements maximum likelihood methods to enable the inference of ancestors of families with more than 10,000 members. GRASP implements partial order graphs (POGs) to represent and infer insertion and deletion events across ancestors, enabling the identification of building blocks for protein engineering. To validate the capacity to engineer novel proteins from realistic data, we predicted ancestor sequences across three distinct enzyme families: glucose-methanol-choline (GMC) oxidoreductases, cytochromes P450, and dihydroxy/sugar acid dehydratases (DHAD). All tested ancestors demonstrated enzymatic activity. Our study demonstrates the ability of GRASP (1) to support large data sets over 10,000 sequences and (2) to employ insertions and deletions to identify building blocks for engineering biologically active ancestors, by exploring variation over evolutionary time.

Hunt for α-amylase from metagenome and strategies to improve its thermostability: a systematic review

With the advent of green chemistry, the use of enzymes in industrial processes serves as an alternative to the conventional chemical catalysts. A high demand for sustainable processes for catalysis has brought a significant attention to hunt for novel enzymes. Among various hydrolases, the α-amylase has a gamut of biotechnological applications owing to its pivotal role in starch-hydrolysis. Industrial demand requires enzymes with thermostability and to ameliorate this crucial property, various methods such as protein engineering, directed evolution and enzyme immobilisation strategies are devised. Besides the traditional culture-dependent approach, metagenome from uncultured bacteria serves as a bountiful resource for novel genes/biocatalysts. Exploring the extreme-niches metagenome, advancements in protein engineering and biotechnology tools encourage the mining of novel α-amylase and its stable variants to tap its robust biotechnological and industrial potential. This review outlines α-amylase and its genetics, its catalytic domain architecture and mechanism of action, and various molecular methods to ameliorate its production. It aims to impart understanding on mechanisms involved in thermostability of α-amylase, cover strategies to screen novel genes from futile habitats and some molecular methods to ameliorate its properties.

Consequences of Substitution Model Selection on Protein Ancestral Sequence Reconstruction

The selection of the best-fitting substitution model of molecular evolution is a traditional step for phylogenetic inferences, including ancestral sequence reconstruction (ASR). However, a few recent studies suggested that applying this procedure does not affect the accuracy of phylogenetic tree reconstruction. Here, we revisited this debate topic by analyzing the influence of selection among substitution models of protein evolution, with focus on exchangeability matrices, on the accuracy of ASR using simulated and real data. We found that the selected best-fitting substitution model produces the most accurate ancestral sequences, especially if the data present large genetic diversity. Indeed, ancestral sequences reconstructed under substitution models with similar exchangeability matrices were similar, suggesting that if the selected best-fitting model cannot be used for the reconstruction, applying a model similar to the selected one is preferred. We conclude that selecting among substitution models of protein evolution is recommended for reconstructing accurate ancestral sequences.

Evolving Perspective on the Origin and Diversification of Cellular Life and the Virosphere

The tree of life (TOL) is a powerful framework to depict the evolutionary history of cellular organisms through time, from our microbial origins to the diversification of multicellular eukaryotes that shape the visible biosphere today. During the past decades, our perception of the TOL has fundamentally changed, in part, due to profound methodological advances, which allowed a more objective approach to studying organismal and viral diversity and led to the discovery of major new branches in the TOL as well as viral lineages. Phylogenetic and comparative genomics analyses of these data have, among others, revolutionized our understanding of the deep roots and diversity of microbial life, the origin of the eukaryotic cell, eukaryotic diversity, as well as the origin, and diversification of viruses. In this review, we provide an overview of some of the recent discoveries on the evolutionary history of cellular organisms and their viruses and discuss a variety of complementary techniques that we consider crucial for making further progress in our understanding of the TOL and its interconnection with the virosphere.

A cocktail of protein engineering strategies: breaking the enzyme bottleneck one by one for high UTP production in vitro

The pyrimidine metabolic pathway is tightly regulated in microorganisms, allowing limited success in metabolic engineering for the production of pathway-related substances. Here, we constructed a four-enzyme coupled system for the in vitro production of uridine triphosphate (UTP). The enzymes used include nucleoside kinase, uridylate kinase, nucleoside diphosphate kinase, and polyphosphate kinase for energy regeneration. All these enzymes are derived from extremophiles. To increase the total and unit time yield of the product, three enzymes other than polyphosphate kinase were modified separately by multiple protein engineering strategies. A nucleoside kinase variant with increased specific activity from 2.7 U/mg to 36.5 U/mg, a uridylate kinase variant (specific activity of 37.1 U/mg) with a 5.2-fold increase in thermostability, and a nucleoside diphosphate kinase variant with a 2-fold increase in specific activity to over 900 U/mg were obtained, respectively. The reaction conditions of the coupled system were further optimized, and a two-stage method was taken to avoid the problem of enzymatic pH adaptation mismatch. Under optimal conditions, this system can produce more than 65 mM UTP (31.5 g/L) in 3.0 h. The substrate conversion rate exceeded 98% and the maximum UTP productivity reached 40 mM/h. This article is protected by copyright. All rights reserved.

Ancestral Sequence Reconstruction for Exploring Alkaloid Evolution

The complex and bioactive monoterpene indole alkaloids (MIAs) found in Catharanthus roseus and related species are the products of many millions of years of evolution through mutation and natural selection. Ancestral sequence reconstruction (ASR) is a method that combines phylogenetic analysis and experimental biochemistry to infer details about past events in protein evolution. Here, I propose that ASR could be leveraged to understand how enzymes catalyzing the formation of complex alkaloids arose over evolutionary time. I discuss the steps of ASR, including sequence selection, multiple sequence alignment, tree inference, and the generation and characterization of inferred ancestral enzymes.

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.

Phylogenetic articulation of uric acid evolution in mammals and how it informs a therapeutic uricase

The role of uric acid during primate evolution has remained elusive ever since it was discovered over 100 years ago that humans have unusually high levels of the small molecule in our serum. It has been difficult to generate a neutral or adaptive explanation in part because the uricase enzyme evolved to become a pseudogene in apes thus masking typical signals of sequence evolution. Adding to the difficulty is a lack of clarity on the functional role of uric acid in apes. One popular hypothesis proposes that uric acid is a potent antioxidant that increased in concentration to compensate for the lack of vitamin C synthesis in primate species ∼65 million years ago (Mya). Here, we have expanded on our previous work with resurrected ancient uricase proteins to better resolve the reshaping of uricase enzymatic activity prior to ape evolution. Our results suggest that the pivotal death-knell to uricase activity occurred between 20-30 Mya despite small sequential modifications to its catalytic efficiency for the tens of millions of years since primates lost their ability to synthesize vitamin C, and thus the two appear uncorrelated. We also use this opportunity to demonstrate how molecular evolution can contribute to biomedicine by presenting ancient uricases to human immune cells that assay for innate reactivity against foreign antigens. A highly stable and highly catalytic ancient uricase is shown to elicit a lower immune response in more human haplotypes than other uricases currently in therapeutic development.

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.

CRISPR-Cas9-mediated reactivation of the uricase pseudogene in human cells prevents acute hyperuricemia

The utility of CRISPR-Cas9 to repair or reverse diseased states that arise from recent genetic mutations in the human genome is now widely appreciated. The use of CRISPR to "design" the outcomes of biology is challenged by both specialized ethicists and the general public. Less of a focus, however, is the ability of CRISPR to provide metabolic supplements or prophylactic molecules that improve long-term human health by overwriting ancient evolutionary events. Here, we use CRISPR to genomically integrate a functional uricase gene that encodes an enzymatically active protein into the human genome. These uricase-producing cells are able to reduce or even eliminate high concentrations of exogenous uric acid despite the enzyme being localized to peroxisomes. Our evolutionary engineered cells represent the first instance of the primate ape lineage expressing a functional uricase encoded in the genome within the last 20 million years. We anticipate that human cells expressing uricase will help prevent hyperuricemia (including gout) as well as hypertension and will help protect against fatty liver disease in the future.

Setting the stage for evolution of a new enzyme

The evolution of novel enzymes has fueled the diversification of life on earth for billions of years. Insights into events that set the stage for the evolution of a new enzyme can be obtained from ancestral reconstruction and laboratory evolution. Ancestral reconstruction can reveal the emergence of a promiscuous activity in a pre-existing protein and the impact of subsequent mutations that enhance a new activity. Laboratory evolution provides a more holistic view by revealing mutations elsewhere in the genome that indirectly enhance the level of a newly important enzymatic activity. This review will highlight recent studies that probe the early stages of the evolution of a new enzyme from these complementary points of view.

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.