Function and structure are strongly coupled in obligated oligomers such as Triosephosphate isomerase (TIM). In animals and fungi, TIM monomers are inactive and unstable. Previously, we used ancestral sequence reconstruction to study TIM evolution and found that before these lineages diverged, the last opisthokonta common ancestor of TIM (LOCATIM) was an obligated oligomer that resembles those of extant TIMs. Notably, calorimetric evidence indicated that ancestral TIM monomers are more structured than extant ones. To further increase confidence about the function, structure, and stability of the LOCATIM, in this work, we applied two different inference methodologies and the worst plausible case scenario for both of them, to infer four sequences of this ancestor and test the robustness of their physicochemical properties. The extensive biophysical characterization of the four reconstructed sequences of LOCATIM showed very similar hydrodynamic and spectroscopic properties, as well as ligand-binding energetics and catalytic parameters. Their 3D structures were also conserved. Although differences were observed in melting temperature, all LOCATIMs showed reversible urea-induced unfolding transitions, and for those that reached equilibrium, high conformational stability was estimated (ΔGTot = 40.6-46.2 kcal/mol). The stability of the inactive monomeric intermediates was also high (ΔGunf = 12.6-18.4 kcal/mol), resembling some protozoan TIMs rather than the unstable monomer observed in extant opisthokonts. A comparative analysis of the 3D structure of ancestral and extant TIMs shows a correlation between the higher stability of the ancestral monomers with the presence of several hydrogen bonds located in the \"bottom\" part of the barrel.

Geranylgeranyl diphosphate synthase (GGPPS) is the crucial bottleneck in carotenoid biosynthesis. However, low activity limits the broad application of GGPPS. In this study, OsGGPPS1 in rice was engineered based on ancestral sequence reconstruction (ASR) and semirational design to improve the catalytic performances of existing GGPPS. The better mutant of A22R/A26P with improved enzyme activity was generated based on ASR. Additionally, the improved enzyme activity of mutants as V162A/M218S/F227Y was designed using a semirational design. The combinatorial assembly of the d-OsGGPPS1 mutant (A22R/A26P/V162A/M218S/F227Y) exhibited higher conversion of IPP and each cosubstrate of DMAPP for 9.8-fold in GPP production, GPP for 6.4-fold in FPP production, and FPP for 1.4-fold in GGPP production relative to wild-type OsGGPPS1 at 25 °C, which showed higher conversion than wild-type OsGGPPS1 at temperatures as high as 50 °C. The successful design of OsGGPPS1 was representative of protein engineering, which will shed new light on GGPPS engineering and active plant pigment resource utilization.

New protein-coding genes can evolve from previously noncoding genomic regions through a process known as de novo gene emergence. Evidence suggests that this process has likely occurred throughout evolution and across the tree of life. Yet, confidently identifying de novo emerged genes remains challenging. Ancestral sequence reconstruction is a promising approach for inferring whether a gene has emerged de novo or not, as it allows us to inspect whether a given genomic locus ancestrally harbored protein-coding capacity. However, the use of ancestral sequence reconstruction in the context of de novo emergence is still in its infancy and its capabilities, limitations, and overall potential are largely unknown. Notably, it is difficult to formally evaluate the protein-coding capacity of ancestral sequences, particularly when new gene candidates are short. How well-suited is ancestral sequence reconstruction as a tool for the detection and study of de novo genes? Here, we address this question by designing an ancestral sequence reconstruction workflow incorporating different tools and sets of parameters and by introducing a formal criterion that allows to estimate, within a desired level of confidence, when protein-coding capacity originated at a particular locus. Applying this workflow on ∼2,600 short, annotated budding yeast genes (<1,000 nucleotides), we found that ancestral sequence reconstruction robustly predicts an ancient origin for the most widely conserved genes, which constitute \"easy\" cases. For less robust cases, we calculated a randomization-based empirical P-value estimating whether the observed conservation between the extant and ancestral reading frame could be attributed to chance. This formal criterion allowed us to pinpoint a branch of origin for most of the less robust cases, identifying 49 genes that can unequivocally be considered de novo originated since the split of the Saccharomyces genus, including 37 Saccharomyces cerevisiae-specific genes. We find that for the remaining equivocal cases we cannot rule out different evolutionary scenarios including rapid evolution, multiple gene losses, or a recent de novo origin. Overall, our findings suggest that ancestral sequence reconstruction is a valuable tool to study de novo gene emergence but should be applied with caution and awareness of its limitations.

Terpenes constitute one of the largest family of natural products with potent applications as renewable platform chemicals and medicines. The low activity, selectivity and stability displayed by terpene biosynthetic machineries can constitute an obstacle towards achieving expedient biosynthesis of terpenoids in processes that adhere to the 12 principles of green chemistry. Accordingly, engineering of terpene synthase enzymes is a prerequisite for industrial biotechnology applications, but obstructed by their complex catalysis that depend on reactive carbocationic intermediates that are prone to undergo bifurcation mechanisms. Rational redesign of terpene synthases can be tedious and requires high-resolution structural information, which is not always available. Furthermore, it has proven difficult to link sequence space of terpene synthase enzymes to specific product profiles. Herein, the author shows how ancestral sequence reconstruction (ASR) can favorably be used as a protein engineering tool in the redesign of terpene synthases without the need of a structure, and without excessive screening. A detailed workflow of ASR is presented along with associated limitations, with a focus on applying this methodology on terpene synthases. From selected examples of both class I and II enzymes, the author advocates that ancestral terpene cyclases constitute valuable assets to shed light on terpene-synthase catalysis and in enabling accelerated biosynthesis.

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.

The dynamics behavior of a protein is essential for its functionality. Here, Doucet et al. demonstrate how the evolutionary analysis of conformational pathways within a protein family serves to identify common core scaffolds that accommodate branch-specific functional regions controlled by flexibility switches, offering a model for evolutionary-dynamics based protein design.

Life depends on a conserved set of chemical energy currencies that are relics of early biochemistry. One of these is ATP, a molecule that, when paired with a divalent metal ion such as Mg2+, can be hydrolyzed to support numerous cellular and molecular processes. Despite its centrality to extant biochemistry, it is unclear whether ATP supported the function of ancient enzymes. We investigate the evolutionary necessity of ATP by experimentally reconstructing an ancestral variant of the N2-reducing enzyme nitrogenase. The Proterozoic ancestor is predicted to be ~540-2,300 million years old, post-dating the Great Oxidation Event. Growth rates under nitrogen-fixing conditions are ~80% of those of wild type in Azotobacter vinelandii. In the extant enzyme, the hydrolysis of two MgATP is coupled to electron transfer to support substrate reduction. The ancestor has a strict requirement for ATP with no other nucleotide triphosphate analogs (GTP, ITP, and UTP) supporting activity. Alternative divalent metal ions (Fe2+, Co2+, and Mn2+) support activity with ATP but with diminished activities compared to Mg2+, similar to the extant enzyme. Additionally, it is shown that the ancestor has an identical efficiency in ATP hydrolyzed per electron transferred to the extant of two. Our results provide direct laboratory evidence of ATP usage by an ancient enzyme.IMPORTANCELife depends on energy-carrying molecules to power many sustaining processes. There is evidence that these molecules may predate the rise of life on Earth, but how and when these dependencies formed is unknown. The resurrection of ancient enzymes provides a unique tool to probe the enzyme\'s function and usage of energy-carrying molecules, shedding light on their biochemical origins. Through experimental reconstruction, this research investigates the ancestral dependence of a nitrogen-fixing enzyme on the energy carrier ATP, a requirement for function in the modern enzyme. We show that the resurrected ancestor does not have generalist nucleotide specificity. Rather, the ancestor has a strict requirement for ATP, like the modern enzyme, with similar function and efficiency. The findings elucidate the early-evolved necessity of energy-yielding molecules, delineating their role in ancient biochemical processes. Ultimately, these insights contribute to unraveling the intricate tapestry of evolutionary biology and the origins of life-sustaining dependencies.

Despite having important biological implications, insertion, and deletion (indel) events are often disregarded or mishandled during phylogenetic inference. In multiple sequence alignment, indels are represented as gaps and are estimated without considering the distinct evolutionary history of insertions and deletions. Consequently, indels are usually excluded from subsequent inference steps, such as ancestral sequence reconstruction and phylogenetic tree search. Here, we introduce indel-aware parsimony (indelMaP), a novel way to treat gaps under the parsimony criterion by considering insertions and deletions as separate evolutionary events and accounting for long indels. By identifying the precise location of an evolutionary event on the tree, we can separate overlapping indel events and use affine gap penalties for long indel modeling. Our indel-aware approach harnesses the phylogenetic signal from indels, including them into all inference stages. Validation and comparison to state-of-the-art inference tools on simulated data show that indelMaP is most suitable for densely sampled datasets with closely to moderately related sequences, where it can reach alignment quality comparable to probabilistic methods and accurately infer ancestral sequences, including indel patterns. Due to its remarkable speed, our method is well suited for epidemiological datasets, eliminating the need for downsampling and enabling the exploitation of the additional information provided by dense taxonomic sampling. Moreover, indelMaP offers new insights into the indel patterns of biologically significant sequences and advances our understanding of genetic variability by considering gaps as crucial evolutionary signals rather than mere artefacts.

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.

Xylanase is an essential component used to hydrolyze the xylan in wheat flour to enhance the quality of bread. Presently, cold-activated xylanase is popularly utilized to aid in the development of dough. In this study, ancestral sequence reconstruction and molecular docking of xylanase and wheat xylan were used to enhance the activity and stability of a thermophilic xylanase. The results indicated that the ancestral enzyme TmxN3 exhibited significantly improved activity and thermal stability. The Vmax increased by 2.7 times, and the catalytic efficiency (Kcat/Km) increased by 1.7 times in comparison to TmxB. After being incubated at 100 °C for 120 min, it still retained 87.3% of its activity, and the half-life in 100 °C was 330 min, while the wild type xylanase was only 55 min. This resulted in an improved shelf life of bread, while adding TmxN3 considerably enhanced its quality with excellent volume and reduced hardness, chewiness, and gumminess. The results showed that the hardness was reduced by 55.2%, the chewiness was reduced by 40.11%, and the gumminess was reduced by 53.52%. To facilitate its industrial application, we further optimized the production conditions in a 5L bioreactor, and the xylanase activity reached 1.52 × 106 U/mL culture.