Ancestral protein reconstruction: techniques and 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 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.

Evolutionary divergence of TLR9 through ancestral sequence reconstruction

The transmembrane pattern recognition receptor, Toll-like receptor (TLR), are best known for their roles in innate immunity via recognition of pathogen and initiation of signaling response. Mammalian TLRs recognize molecular patterns associated with pathogens and initiate innate immune response. We have studied the evolutionary diversity of mammalian TLR genes for differences in immunological response. Reconstruction of ancestral sequences is a key aspect of the molecular evolution of TLR to track changes across the TLR genes. The comprehensive analysis of mammalian TLRs revealed a distinct pattern of evolution of TLR9. Various sequence-based features such as amino acid usage, hydrophobicity, GC content, and evolutionary constraints are found to influence the divergence of TLR9 from other TLRs. Ancestral sequence reconstruction analysis also revealed that the gradual evolution of TLR genes in several ancestral lineages leads to the distinct pattern of TLR9. It demonstrates evolutionary divergence with the progressive accumulation of mutations results in the distinct pattern of TLR9.

Computational scoring and experimental evaluation of enzymes generated by neural networks

In recent years, generative protein sequence models have been developed to sample novel sequences. However, predicting whether generated proteins will fold and function remains challenging. We evaluate a set of 20 diverse computational metrics to assess the quality of enzyme sequences produced by three contrasting generative models: ancestral sequence reconstruction, a generative adversarial network and a protein language model. Focusing on two enzyme families, we expressed and purified over 500 natural and generated sequences with 70-90% identity to the most similar natural sequences to benchmark computational metrics for predicting in vitro enzyme activity. Over three rounds of experiments, we developed a computational filter that improved the rate of experimental success by 50-150%. The proposed metrics and models will drive protein engineering research by serving as a benchmark for generative protein sequence models and helping to select active variants for experimental testing.

Click, Compute, Create: A Review of Web-based Tools for Enzyme Engineering

Enzyme engineering, though pivotal across various biotechnological domains, is often plagued by its time-consuming and labor-intensive nature. This review aims to offer an overview of supportive in silico methodologies for this demanding endeavor. Starting from methods to predict protein structures, to classification of their activity and even the discovery of new enzymes we continue with describing tools used to increase thermostability and production yields of selected targets. Subsequently, we discuss computational methods to modulate both, the activity as well as selectivity of enzymes. Last, we present recent approaches based on cutting-edge machine learning methods to redesign enzymes. With exception of the last chapter, there is a strong focus on methods easily accessible via web-interfaces or simple Python-scripts, therefore readily useable for a diverse and broad community.

A billion years of evolution manifest in nanosecond protein dynamics

Protein dynamics form a critical bridge between protein structure and function, yet the impact of evolution on ultrafast processes inside proteins remains enigmatic. This study delves deep into nanosecond-scale protein dynamics of a structurally and functionally conserved protein across species separated by almost a billion years, investigating ten homologs in complex with their ligand. By inducing a photo-triggered destabilization of the ligand inside the binding pocket, we resolved distinct kinetic footprints for each homolog via transient infrared spectroscopy. Strikingly, we found a cascade of rearrangements within the protein complex which manifest in time points of increased dynamic activity conserved over hundreds of millions of years within a narrow window. Among these processes, one displays a subtle temporal shift correlating with evolutionary divergence, suggesting reduced selective pressure in the past. Our study not only uncovers the impact of evolution on molecular processes in a specific case, but has also the potential to initiate a field of scientific inquiry within molecular paleontology, where species are compared and classified based on the rapid pace of protein dynamic processes; a field which connects the shortest conceivable time scale in living matter (10[Formula: see text] s) with the largest ones (10[Formula: see text] s).

Opportunities and Challenges for Machine Learning-Assisted Enzyme Engineering

Enzymes can be engineered at the level of their amino acid sequences to optimize key properties such as expression, stability, substrate range, and catalytic efficiency-or even to unlock new catalytic activities not found in nature. Because the search space of possible proteins is vast, enzyme engineering usually involves discovering an enzyme starting point that has some level of the desired activity followed by directed evolution to improve its "fitness" for a desired application. Recently, machine learning (ML) has emerged as a powerful tool to complement this empirical process. ML models can contribute to (1) starting point discovery by functional annotation of known protein sequences or generating novel protein sequences with desired functions and (2) navigating protein fitness landscapes for fitness optimization by learning mappings between protein sequences and their associated fitness values. In this Outlook, we explain how ML complements enzyme engineering and discuss its future potential to unlock improved engineering outcomes.

Evolution and variation in amide aminoacyl-tRNA synthesis

The amide proteogenic amino acids, asparagine and glutamine, are two of the twenty amino acids used in translation by all known life. The aminoacyl-tRNA synthetases for asparagine and glutamine, asparaginyl-tRNA synthetase and glutaminyl tRNA synthetase, evolved after the split in the last universal common ancestor of modern organisms. Before that split, life used two-step indirect pathways to synthesize asparagine and glutamine on their cognate tRNAs to form the aminoacyl-tRNA used in translation. These two-step pathways were retained throughout much of the bacterial and archaeal domains of life and eukaryotic organelles. The indirect routes use non-discriminating aminoacyl-tRNA synthetases (non-discriminating aspartyl-tRNA synthetase and non-discriminating glutamyl-tRNA synthetase) to misaminoacylate the tRNA. The misaminoacylated tRNA formed is then transamidated into the amide aminoacyl-tRNA used in protein synthesis by tRNA-dependent amidotransferases (GatCAB and GatDE). The enzymes and tRNAs involved assemble into complexes known as transamidosomes to help maintain translational fidelity. These pathways have evolved to meet the varied cellular needs across a diverse set of organisms, leading to significant variation. In certain bacteria, the indirect pathways may provide a means to adapt to cellular stress by reducing the fidelity of protein synthesis. The retention of these indirect pathways versus acquisition of asparaginyl-tRNA synthetase and glutaminyl tRNA synthetase in lineages likely involves a complex interplay of the competing uses of glutamine and asparagine beyond translation, energetic costs, co-evolution between enzymes and tRNA, and involvement in stress response that await further investigation.

Perspectives on Computational Enzyme Modeling: From Mechanisms to Design and Drug Development

Understanding enzyme mechanisms is essential for unraveling the complex molecular machinery of life. In this review, we survey the field of computational enzymology, highlighting key principles governing enzyme mechanisms and discussing ongoing challenges and promising advances. Over the years, computer simulations have become indispensable in the study of enzyme mechanisms, with the integration of experimental and computational exploration now established as a holistic approach to gain deep insights into enzymatic catalysis. Numerous studies have demonstrated the power of computer simulations in characterizing reaction pathways, transition states, substrate selectivity, product distribution, and dynamic conformational changes for various enzymes. Nevertheless, significant challenges remain in investigating the mechanisms of complex multistep reactions, large-scale conformational changes, and allosteric regulation. Beyond mechanistic studies, computational enzyme modeling has emerged as an essential tool for computer-aided enzyme design and the rational discovery of covalent drugs for targeted therapies. Overall, enzyme design/engineering and covalent drug development can greatly benefit from our understanding of the detailed mechanisms of enzymes, such as protein dynamics, entropy contributions, and allostery, as revealed by computational studies. Such a convergence of different research approaches is expected to continue, creating synergies in enzyme research. This review, by outlining the ever-expanding field of enzyme research, aims to provide guidance for future research directions and facilitate new developments in this important and evolving field.

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.

Impact of ancestral sequence reconstruction on mechanistic and structural enzymology

Ancestral sequence reconstruction (ASR) provides insight into the changes within a protein sequence across evolution. More specifically, it can illustrate how specific amino acid changes give rise to different phenotypes within a protein family. Over the last few decades it has established itself as a powerful technique for revealing molecular common denominators that govern enzyme function. Here, we describe the strength of ASR in unveiling catalytic mechanisms and emerging phenotypes for a range of different proteins, also highlighting biotechnological applications the methodology can provide.

Architectural digest: Thermodynamic stability and domain structure of a consensus monomeric globin

Artificial proteins representing the consensus of a set of homologous sequences have attracted attention for their increased thermodynamic stability and conserved activity. Here, we applied the consensus approach to a b-type heme-binding protein to inspect the contribution of a dissociable cofactor to enhanced stability and the chemical consequences of creating a generic heme environment. We targeted the group 1 truncated hemoglobin (TrHb1) subfamily of proteins for their small size (∼120 residues) and ease of characterization. The primary structure, derived from a curated set of ∼300 representative sequences, yielded a highly soluble consensus globin (cGlbN) enriched in acidic residues. Optical and NMR spectroscopies revealed high-affinity heme binding in the expected site and in two orientations. At neutral pH, proximal and distal iron coordination was achieved with a pair of histidine residues, as observed in some natural TrHb1s, and with labile ligation on the distal side. As opposed to studied TrHb1s, which undergo additional folding upon heme binding, cGlbN displayed the same extent of secondary structure whether the heme was associated with the protein or not. Denaturation required guanidine hydrochloride and showed that apo- and holoprotein unfolded in two transitions-the first (occurring with a midpoint of ∼2 M) was shifted to higher denaturant concentration in the holoprotein (∼3.7 M) and reflected stabilization due to heme binding, while the second transition (∼6.2 M) was common to both forms. Thus, the consensus sequence stabilized the protein but exposed the existence of two separately cooperative subdomains within the globin architecture, masked as one single domain in TrHb1s with typical stabilities. The results suggested ways in which specific chemical or thermodynamic features may be controlled in artificial heme proteins.

Repeated evolution of similar phenotypes: Integrating comparative methods with developmental pathways

Repeated phenotypes, often referred to as 'homoplasies' in cladistic analyses, may evolve through changes in developmental processes. Genetic bases of recurrent evolution gained attention and have been studied in the past years using approaches that combine modern analytical phylogenetic tools with the stunning assemblage of new information on developmental mechanisms. In this review, we evaluated the topic under an integrated perspective, revisiting the classical definitions of convergence and parallelism and detailing comparative methods used to evaluate evolution of repeated phenotypes, which include phylogenetic inference, estimates of evolutionary rates and reconstruction of ancestral states. We provide examples to illustrate how a given methodological approach can be used to identify evolutionary patterns and evaluate developmental mechanisms associated with the intermittent expression of a given trait along the phylogeny. Finally, we address why repeated trait loss challenges strict definitions of convergence and parallelism, discussing how changes in developmental pathways might explain the high frequency of repeated trait loss in specific lineages.

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.

Ion-Powered Rotary Motors: Where Did They Come from and Where They Are Going?

Molecular motors are found in many living organisms. One such molecular machine, the ion-powered rotary motor (IRM), requires the movement of ions across a membrane against a concentration gradient to drive rotational movement. The bacterial flagellar motor (BFM) is an example of an IRM which relies on ion movement through the stator proteins to generate the rotation of the flagella. There are many ions which can be used by the BFM stators to power motility and different ions can be used by a single bacterium expressing multiple stator variants. The use of ancestral sequence reconstruction (ASR) and functional analysis of reconstructed stators shows promise for understanding how these proteins evolved and when the divergence in ion use may have occurred. In this review, we discuss extant BFM stators and the ions that power them as well as recent examples of the use of ASR to study ion-channel selectivity and how this might be applied to further study of the BFM stator complex.

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.

Immunoanalytical Approach for Detecting and Identifying Ancestral Peptide Biomarkers in Early Earth Analogue Environments

Several mass spectrometry and spectroscopic techniques have been used in the search for molecular biomarkers on Mars. A major constraint is their capability to detect and identify large and complex compounds such as peptides or other biopolymers. Multiplex immunoassays can detect these compounds, but antibodies must be produced for a large number of sequence-dependent molecular targets. Ancestral Sequence Reconstruction (ASR) followed by protein "resurrection" in the lab can help to narrow the selection of targets. Herein, we propose an immunoanalytical method to identify ancient and universally conserved protein/peptide sequences as targets for identifying ancestral biomarkers in nature. We have developed, tested, and validated this approach by producing antibodies to eight previously described ancestral resurrected proteins (three β-lactamases, three thioredoxins, one Elongation Factor Tu, and one RuBisCO, all of them theoretically dated as Precambrian), and used them as a proxy to search for any potential feature of them that could be present in current natural environments. By fluorescent sandwich microarray immunoassays (FSMI), we have detected positive immunoreactions with antibodies to the oldest β-lactamase and thioredoxin proteins (ca. 4 Ga) in samples from a hydrothermal environment. Fine epitope mapping and inhibitory immunoassays allowed the identification of well-conserved epitope peptide sequences that resulted from ASR and were present in the sample. We corroborated these results by metagenomic sequencing and found several genes encoding analogue proteins with significant matches to the peptide epitopes identified with the antibodies. The results demonstrated that peptides inferred from ASR studies have true counterpart analogues in Nature, which validates and strengthens the well-known ASR/protein resurrection technique and our immunoanalytical approach for investigating ancient environments and metabolisms on Earth and elsewhere.

Structural and Functional Characterization of the Ureidoacrylate Amidohydrolase RutB from Escherichia coli

We present a detailed structure-function analysis of the ureidoacrylate amidohydrolase RutB from Eschericha coli, which is an essential enzyme of the Rut pathway for pyrimidine utilization. Crystals of selenomethionine-labeled RutB were produced, which allowed us to determine the first structure of the enzyme at a resolution of 1.9 Å and to identify it as a new member of the isochorismatase-like hydrolase family. RutB was co-crystallized with the substrate analogue ureidopropionate, revealing the mode of substrate binding. Mutation of residues constituting the catalytic triad (D24A, D24N, K133A, C166A, C166S, C166T, C166Y) resulted in complete inactivation of RutB, whereas mutation of other residues close to the active site (Y29F, Y35F, N72A, W74A, W74F, E80A, E80D, S92A, S92T, S92Y, Q105A, Y136A, Y136F) leads to distinct changes of the turnover number (k_cat) and/or the Michaelis constant (K_M). The results of our structural and mutational studies allowed us to assign specific functions to individual residues and to formulate a plausible reaction mechanism for RutB.

Ancestral Sequence Reconstruction of the Ribosomal Protein uS8 and Reduction of Amino Acid Usage to a Smaller Alphabet

Understanding the origin and early evolution of proteins is important for unveiling how the RNA world developed into an RNA-protein world. Because the composition of organic molecules in the Earth's primitive environment was plausibly not as diverse as today, the number of different amino acids used in early protein synthesis is likely to be substantially less than the current 20 proteinogenic residues. In this study, we have explored the thermal stability and RNA binding of ancestral variants of the ribosomal protein uS8 constructed from a reduced-alphabet of amino acids. First, we built a phylogenetic tree based on the amino acid sequences of uS8 from multiple extant organisms and used the tree to infer two plausible amino acid sequences corresponding to the last bacterial common ancestor of uS8. Both ancestral proteins were thermally stable and bound to an RNA fragment. By eliminating individual amino acid letters and monitoring thermal stability and RNA binding in the resulting proteins, we reduced the size of the amino acid set constituting one of the ancestral proteins, eventually finding that convergent sequences consisting of 15- or 14-amino acid alphabets still folded into stable structures that bound to the RNA fragment. Furthermore, a simplified variant reconstructed from a 13-amino-acid alphabet retained affinity for the RNA fragment, although it lost conformational stability. Collectively, RNA-binding activity may be achieved with a subset of the current 20 amino acids, raising the possibility of a simpler composition of RNA-binding proteins in the earliest stage of protein evolution.

Evolution of CRISPR-associated endonucleases as inferred from resurrected proteins

Clustered regularly interspaced short palindromic repeats (CRISPR)-associated Cas9 is an effector protein that targets invading DNA and plays a major role in the prokaryotic adaptive immune system. Although Streptococcus pyogenes CRISPR-Cas9 has been widely studied and repurposed for applications including genome editing, its origin and evolution are poorly understood. Here, we investigate the evolution of Cas9 from resurrected ancient nucleases (anCas) in extinct firmicutes species that last lived 2.6 billion years before the present. We demonstrate that these ancient forms were much more flexible in their guide RNA and protospacer-adjacent motif requirements compared with modern-day Cas9 enzymes. Furthermore, anCas portrays a gradual palaeoenzymatic adaptation from nickase to double-strand break activity, exhibits high levels of activity with both single-stranded DNA and single-stranded RNA targets and is capable of editing activity in human cells. Prediction and characterization of anCas with a resurrected protein approach uncovers an evolutionary trajectory leading to functionally flexible ancient enzymes.

Ancestral reconstruction reveals catalytic inactivation of activation-induced cytidine deaminase concomitant with cold water adaption in the Gadiformes bony fish

BACKGROUND

Antibody affinity maturation in vertebrates requires the enzyme activation-induced cytidine deaminase (AID) which initiates secondary antibody diversification by mutating the immunoglobulin loci. AID-driven antibody diversification is conserved across jawed vertebrates since bony and cartilaginous fish. Two exceptions have recently been reported, the Pipefish and Anglerfish, in which the AID-encoding aicda gene has been lost. Both cases are associated with unusual reproductive behavior, including male pregnancy and sexual parasitism. Several cold water fish in the Atlantic cod (Gadinae) family carry an aicda gene that encodes for a full-length enzyme but lack affinity-matured antibodies and rely on antibodies of broad antigenic specificity. Hence, we examined the functionality of their AID.

RESULTS

By combining genomics, transcriptomics, immune responsiveness, and functional enzymology of AID from 36 extant species, we demonstrate that AID of that Atlantic cod and related fish have extremely lethargic or no catalytic activity. Through ancestral reconstruction and functional enzymology of 71 AID enzymes, we show that this enzymatic inactivation likely took place relatively recently at the emergence of the true cod family (Gadidae) from their ancestral Gadiformes order. We show that this AID inactivation is not only concordant with the previously shown loss of key adaptive immune genes and expansion of innate and cell-based immune genes in the Gadiformes but is further reflected in the genomes of these fish in the form of loss of AID-favored sequence motifs in their immunoglobulin variable region genes.

CONCLUSIONS

Recent demonstrations of the loss of the aicda gene in two fish species challenge the paradigm that AID-driven secondary antibody diversification is absolutely conserved in jawed vertebrates. These species have unusual reproductive behaviors forming an evolutionary pressure for a certain loss of immunity to avoid tissue rejection. We report here an instance of catalytic inactivation and functional loss of AID rather than gene loss in a conventionally reproducing vertebrate. Our data suggest that an expanded innate immunity, in addition to lower pathogenic pressures in a cold environment relieved the pressure to maintain robust secondary antibody diversification. We suggest that in this unique scenario, the AID-mediated collateral genome-wide damage would form an evolutionary pressure to lose AID function.

Enzymatic Conjugation of Modified RNA Fragments by Ancestral RNA Ligase AncT4_2

Oligonucleotide therapeutics have great potential as a next-generation approach to treating intractable diseases. Large quantities of modified DNA/RNA containing xenobiotic nucleic acids (XNAs) must be synthesized before clinical application. In this study, the ancestral RNA ligase AncT4_2 was designed by ancestral sequence reconstruction (ASR) to perform the conjugation reaction of modified RNA fragments. AncT4_2 had superior properties to native RNA ligase 2 from T4 phage (T4Rnl2), including high productivity, a >2.5-fold-higher turnover number, and >10°C higher thermostability. One remarkable point is the broad substrate selectivity of AncT4_2; the activity of AncT4_2 toward 17 of the modified RNA fragments was higher than that of T4Rnl2. The activity was estimated by measuring the conjugation reaction of two RNA strands, 3'-OH (12 bp) and 5'-PO₄ (12 bp), in which the terminal and penultimate positions of the 3'-OH fragment and the first and second positions of the 5'-PO₄ fragment were substituted by 2'-fluoro, 2'-O-methyl, 2'-O-methoxyethyl, and 2'-H, respectively. The enzymatic properties of AncT4_2 allowed the enzyme to conjugate large quantities of double-stranded RNA coding for patisiran (>400 μM level), which was formed by four RNA fragments containing 2'-OMe-substituted nucleic acids. Structural analysis of modeled AncT4_2 suggested that protein dynamics were changed by mutation to Gly or indel during ASR and that this may positively impact the conjugation of modified RNA fragments with the enzyme. AncT4_2 is expected to be a key biocatalyst in synthesizing RNA therapeutics by an enzymatic reaction. IMPORTANCE RNA therapeutics is one of the next-generation medicines for treating various diseases. Our designed ancestral RNA ligase AncT4_2 exhibited excellent enzymatic properties, such as high thermal stability, productivity, specific activity, and broad substrate selectivity compared to native enzymes. These advantages create the potential for AncT4_2 to be applied in conjugating the modified RNA fragments containing various xenobiotic nucleic acids. In addition, patisiran, a known polyneuropathy therapeutic, could be synthesized from four fragmented oligonucleotides at a preparative scale. Taken together, these findings indicate AncT4_2 could open the door to synthesizing RNA therapeutics by enzymatic reaction at large-scale production.

Reaction Mechanism of Ancestral l-Lys α-Oxidase from Caulobacter Species Studied by Biochemical, Structural, and Computational Analysis

The flavin-dependent amine oxidase superfamily contains various l-amino acid oxidases (LAAOs) bearing different substrate specificities and enzymatic properties. LAAOs catalyze the oxidation of the α-amino group of l-amino acids (L-AAs) to produce imino acids and H₂O₂. In this study, an ancestral l-Lys α-oxidase (AncLLysO2) was designed utilizing genome-mined sequences from the Caulobacter species. The AncLLysO2 exhibited high specificity toward l-Lys; the k _cat/K _m values toward l-Lys were one and two orders larger than those of l-Arg and l-ornithine, respectively. Liquid chromatography-high resolution mass spectrometry analysis indicated that AncLLysO2 released imino acid immediately from the active site after completion of oxidation of the α-amino group. Crystal structures of the ligand-free, l-Lys- and l-Arg-bound forms of AncLLysO2 were determined at 1.4-1.6 Å resolution, indicating that the active site of AncLLysO2 kept an open state during the reaction and more likely to release products. The structures also indicated the substrate recognition mechanism of AncLLysO2; ε-amino, α-amino, and carboxyl groups of l-Lys formed interactions with Q357, A551, and R77, respectively. Biochemical and molecular dynamics simulation analysis of AncLLysO2 indicated that active site residues that indirectly interact with the substrate are also important to exhibit high activity; for example, the aromatic group of Y219 is important to ensure that the l-Lys substrate is placed in the correct position to allow the reaction to proceed efficiently. Taken together, we propose the reaction mechanism of AncLLysO2.

Probing the origin of prion protein misfolding via reconstruction of ancestral proteins

Prion diseases are fatal neurodegenerative diseases caused by pathogenic misfolding of the prion protein, PrP. They are transmissible between hosts, and sometimes between different species, as with transmission of bovine spongiform encephalopathy to humans. Although PrP is found in a wide range of vertebrates, prion diseases are seen only in certain mammals, suggesting that infectious misfolding was a recent evolutionary development. To explore when PrP acquired the ability to misfold infectiously, we reconstructed the sequences of ancestral versions of PrP from the last common primate, primate-rodent, artiodactyl, placental, bird, and amniote. Recombinant ancestral PrPs were then tested for their ability to form β-sheet aggregates, either spontaneously or when seeded with infectious prion strains from human, cervid, or rodent species. The ability to aggregate developed after the oldest ancestor (last common amniote), and aggregation capabilities diverged along evolutionary pathways consistent with modern-day susceptibilities. Ancestral bird PrP could not be seeded with modern-day prions, just as modern-day birds are resistant to prion disease. Computational modeling of structures suggested that differences in helix 2 could account for the resistance of ancestral bird PrP to seeding. Interestingly, ancestral primate PrP could be converted by all prion seeds, including both human and cervid prions, raising the possibility that species descended from an ancestral primate have retained the susceptibility to conversion by cervid prions. More generally, the results suggest that susceptibility to prion disease emerged prior to ~100 million years ago, with placental mammals possibly being generally susceptible to disease.

Substrate Affinity Versus Catalytic Efficiency: Ancestral Sequence Reconstruction of tRNA Nucleotidyltransferases Solves an Enzyme Puzzle

In tRNA maturation, CCA-addition by tRNA nucleotidyltransferase is a unique and highly accurate reaction. While the mechanism of nucleotide selection and polymerization is well understood, it remains a mystery why bacterial and eukaryotic enzymes exhibit an unexpected and surprisingly low tRNA substrate affinity while they efficiently catalyze the CCA-addition. To get insights into the evolution of this high-fidelity RNA synthesis, the reconstruction and characterization of ancestral enzymes is a versatile tool. Here, we investigate a reconstructed candidate of a 2 billion years old CCA-adding enzyme from Gammaproteobacteria and compare it to the corresponding modern enzyme of Escherichia coli. We show that the ancestral candidate catalyzes an error-free CCA-addition, but has a much higher tRNA affinity compared with the extant enzyme. The consequence of this increased substrate binding is an enhanced reverse reaction, where the enzyme removes the CCA end from the mature tRNA. As a result, the ancestral candidate exhibits a lower catalytic efficiency in vitro as well as in vivo. Furthermore, the efficient tRNA interaction leads to a processive polymerization, while the extant enzyme catalyzes nucleotide addition in a distributive way. Thus, the modern enzymes increased their polymerization efficiency by lowering the binding affinity to tRNA, so that CCA synthesis is efficiently promoted due to a reduced reverse reaction. Hence, the puzzling and at a first glance contradicting and detrimental weak substrate interaction represents a distinct activity enhancement in the evolution of CCA-adding enzymes.

Directed reconstruction of a novel ancestral alcohol dehydrogenase featuring shifted pH-profile, enhanced thermostability and expanded substrate spectrum

Ancestral enzymes are promising for industrial biotechnology due to high stability and catalytic promiscuity. An effective protocol was developed for the directed resurrection of ancestral enzymes. Employing genome mining with diaryl alcohol dehydrogenase KpADH as the probe, descendant enzymes D10 and D11 were firstly identified. Then through ancestral sequence reconstruction, A64 was resurrected with a specific activity of 4.3 U·mg^-1. The optimum pH of A64 was 7.5, distinct from 5.5 of D10. The T15 50 and T_m values of A64 were 57.5 °C and 61.7 °C, significantly higher than those of the descendant counterpart. Substrate spectrum of A64 was quantitively characterized with a Shannon-Wiener index of 2.38, more expanded than D10, especially, towards bulky ketones in Group A and B. A64 also exhibited higher enantioselectivity. This study provides an effective protocol for constructing of ancestral enzymes and an efficient ancestral enzyme of industrial relevance for asymmetric synthesis of chiral alcohols.

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.

Preparation and Characterization of an Ancient Aminopeptidase Obtained from Ancestral Sequence Reconstruction for L-Carnosine Synthesis

As a biologically active peptide, L-carnosine has been widely used in the pharmaceutical, cosmetic and health care industries due to its various physiological properties. However, relatively little research is available regarding L-carnosine's enzymatic synthesis function. In this study, a potential enzyme sequence with the function of carnosine synthesizing was screened out using the ancestral sequence reconstruction (ASR) technique. Identified with L-carnosine synthesis activity, this enzyme was further confirmed using autoproteolytic phenomenon via Western blot and N-terminal sequencing. After purification, the enzymatic properties of LUCA-DmpA were characterized. The melting temperature (T_m) and denaturation enthalpy (ΔH) of LUCA-DmpA were 60.27 ± 1.24 °C and 1306.00 ± 26.73 kJ·mol^-1, respectively. Circular dichroism (CD) spectroscopy results showed that this ancestral enzyme was composed of α-helix (35.23 ± 0.06%), β-sheet (11.06 ± 0.06%), β-turn (23.67 ± 0.06%) and random coil (32.03 ± 0.06%). The enzyme was characterized with the optimal temperature and pH of 45 °C and 9.0, respectively. Notably, LUCA-DmpA was also characterized with remarkable pH tolerance based on the observation of more than 85% remaining enzymatic activity after incubation at different pH buffers (pH = 6-11) for 12 h. Additionally, rather than being improved or inhibited by metal ions, its enzymatic activity was found to be promoted by introducing organic solvent with a larger log P value. Based on these homology modeling results, the screened LUCA-DmpA is suggested to have further optimization potential, and thereafter to be offered as a promising candidate for real industrial applications.

Biochemical Evolution of a Potent Target of Mosquito Larvicide, 3-Hydroxykynurenine Transaminase

A specific mosquito enzyme, 3-hydroxykynurenine transaminase (HKT), is involved in the processing of toxic metabolic intermediates of the tryptophan metabolic pathway. The HKT enzymatic product, xanthurenic acid, is required for Plasmodium spp. development in the mosquito vectors. Therefore, an inhibitor of HKT may not only be a mosquitocide but also a malaria-transmission blocker. In this work, we present a study investigating the evolution of HKT, which is a lineage-specific duplication of an alanine glyoxylate aminotransferases (AGT) in mosquitoes. Synteny analyses, together with the phylogenetic history of the AGT family, suggests that HKT and the mosquito AGTs are paralogous that were formed via a duplication event in their common ancestor. Furthermore, 41 amino acid sites with significant evidence of positive selection were identified, which could be responsible for biochemical and functional evolution and the stability of conformational stabilization. To get a deeper understanding of the evolution of ligands’ capacity and the ligand-binding mechanism of HKT, the sequence and the 3D homology model of the common ancestor of HKT and AGT in mosquitoes, ancestral mosquito AGT (AncMosqAGT), were inferred and built. The homology model along with 3-hydroxykynurenine, kynurenine, and alanine were used in docking experiments to predict the binding capacity and ligand-binding mode of the new substrates related to toxic metabolites detoxification. Our study provides evidence for the dramatic biochemical evolution of the key detoxifying enzyme and provides potential sites that could hinder the detoxification function, which may be used in mosquito larvicide and design.

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.

Protposer: the web server that readily proposes protein stabilizing mutations with high PPV

Protein stability is a requisite for most biotechnological and medical applications of proteins. As natural proteins tend to suffer from a low conformational stability ex vivo, great efforts have been devoted toward increasing their stability through rational design and engineering of appropriate mutations. Unfortunately, even the best currently used predictors fail to compute the stability of protein variants with sufficient accuracy and their usefulness as tools to guide the rational stabilisation of proteins is limited. We present here Protposer, a protein stabilising tool based on a different approach. Instead of quantifying changes in stability, Protposer uses structure- and sequence-based screening modules to nominate candidate mutations for subsequent evaluation by a logistic regression model, carefully trained to avoid overfitting. Thus, Protposer analyses PDB files in search for stabilization opportunities and provides a ranked list of promising mutations with their estimated success rates (eSR), their probabilities of being stabilising by at least 0.5 kcal/mol. The agreement between eSRs and actual positive predictive values (PPV) on external datasets of mutations is excellent. When Protposer is used with its Optimal kappa selection threshold, its PPV is above 0.7. Even with less stringent thresholds, Protposer largely outperforms FoldX, Rosetta and PoPMusiC. Indicating the PDB file of the protein suffices to obtain a ranked list of mutations, their eSRs and hints on the likely source of the stabilization expected. Protposer is a distinct, straightforward and highly successful tool to design protein stabilising mutations, and it is freely available for academic use at http://webapps.bifi.es/the-protposer.

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.

Using the Evolutionary History of Proteins to Engineer Insertion-Deletion Mutants from Robust, Ancestral Templates Using Graphical Representation of Ancestral Sequence Predictions (GRASP)

Analyzing the natural evolution of proteins by ancestral sequence reconstruction (ASR) can provide valuable information about the changes in sequence and structure that drive the development of novel protein functions. However, ASR has also been used as a protein engineering tool, as it often generates thermostable proteins which can serve as robust and evolvable templates for enzyme engineering. Importantly, ASR has the potential to provide an insight into the history of insertions and deletions that have occurred in the evolution of a protein family. Indels are strongly associated with functional change during enzyme evolution and represent a largely unexplored source of genetic diversity for designing proteins with novel or improved properties. Current ASR methods differ in the way they handle indels; inclusion or exclusion of indels is often managed subjectively, based on assumptions the user makes about the likelihood of each recombination event, yet most currently available ASR tools provide limited, if any, opportunities for evaluating indel placement in a reconstructed sequence. Graphical Representation of Ancestral Sequence Predictions (GRASP) is an ASR tool that maps indel evolution throughout a reconstruction and enables the evaluation of indel variants. This chapter provides a general protocol for performing a reconstruction using GRASP and using the results to create indel variants. The method addresses protein template selection, sequence curation, alignment refinement, tree building, ancestor reconstruction, evaluation of indel variants and approaches to library development.

Mechanochemical Evolution of Disulfide Bonds in Proteins

Disulfide bonds play a pivotal role in the mechanical stability of proteins. Numerous proteins that are known to be exposed to mechanical forces in vivo contain disulfide bonds. The presence of cryptic disulfide bonds in a protein structure may be related to its resistance to an applied mechanical force. Disulfide bonds in proteins tend to be highly conserved but their evolution might be directly related to the evolution of the protein mechanical stability. Hence, tracking the evolution of disulfide bonds in a protein can help to derive crucial stability/function correlations in proteins that are exposed to mechanical forces. Phylogenic analysis and ancestral sequence reconstruction (ASR) allow tracking the evolution of proteins from the past ancestors to our modern days and also establish correlations between proteins from different species. In addition, ASR can be combined with single-molecule force spectroscopy (smFS) to investigate the mechanical properties of proteins including the occurrence and function of disulfide bonds. Here we present a detailed protocol to study the mechanochemical evolution of proteins using a fragment of the giant muscle protein titin as example. The protocol can be easily adapted to AFS studies of any resurrected mechanical force bearing protein of interest.

Methodologies for Microbial Ancestral Sequence Reconstruction

The reconstruction of genetic material of ancestral organisms constitutes a powerful application of evolutionary biology. A fundamental step in this inference is the ancestral sequence reconstruction (ASR), which can be performed with diverse methodologies implemented in computer frameworks. However, most of these methodologies ignore evolutionary properties frequently observed in microbes, such as genetic recombination and complex selection processes, that can bias the traditional ASR. From a practical perspective, here I review methodologies for the reconstruction of ancestral DNA and protein sequences, with particular focus on microbes, and including biases, recommendations, and software implementations. I conclude that microbial ASR is a complex analysis that should be carefully performed and that there is a need for methods to infer more realistic ancestral microbial sequences.

Ancestral Sequence Reconstruction of Ancient Lipase from Family I.3 Bacterial Lipolytic Enzymes

Family I.3 lipase is distinguished from other families by the amino acid sequence and secretion mechanism. Little is known about the evolutionary process driving these differences. This study attempt to understand how the diverse temperature stabilities of bacterial lipases from family I.3 evolved. To achieve that, eighty-three protein sequences sharing a minimum 30% sequence identity with Antarctic Pseudomonas sp. AMS8 lipase were used to infer phylogenetic tree. Using ancestral sequence reconstruction (ASR) technique, the last universal common ancestor (LUCA) sequence of family I.3 was reconstructed. A gene encoding LUCA was synthesized, cloned and expressed as inclusion bodies in E. coli system. Insoluble form of LUCA was refolded using urea dilution method and then purified using affinity chromatography. The purified LUCA exhibited an optimum temperature and pH at 70℃ and 10 respectively. Various metal ions increased or retained the activity of LUCA. LUCA also demonstrated tolerance towards various organic solvents in 25% v/v concentration. The finding from this study could support the understanding on temperature and environment during ancient time. In overall, reconstructed ancestral enzymes have improved physicochemical properties that make them suitable for industrial applications and ASR technique can be employed as a general technique for enzyme engineering.

In Silico Identification and Experimental Validation of Distal Activity-Enhancing Mutations in Tryptophan Synthase

Allostery is a central mechanism for the regulation of multi-enzyme complexes. The mechanistic basis that drives allosteric regulation is poorly understood but harbors key information for enzyme engineering. In the present study, we focus on the tryptophan synthase complex that is composed of TrpA and TrpB subunits, which allosterically activate each other. Specifically, we develop a rational approach for identifying key amino acid residues of TrpB distal from the active site. Those residues are predicted to be crucial for shifting the inefficient conformational ensemble of the isolated TrpB to a productive ensemble through intra-subunit allosteric effects. The experimental validation of the conformationally driven TrpB design demonstrates its superior stand-alone activity in the absence of TrpA, comparable to those enhancements obtained after multiple rounds of experimental laboratory evolution. Our work evidences that the current challenge of distal active site prediction for enhanced function in computational enzyme design has become within reach.

Comprehensive mutagenesis to identify amino acid residues contributing to the difference in thermostability between two originally thermostable ancestral proteins

Further improvement of the thermostability of inherently thermostable proteins is an attractive challenge because more thermostable proteins are industrially more useful and serve as better scaffolds for protein engineering. To establish guidelines that can be applied for the rational design of hyperthermostable proteins, we compared the amino acid sequences of two ancestral nucleoside diphosphate kinases, Arc1 and Bac1, reconstructed in our previous study. Although Bac1 is a thermostable protein whose unfolding temperature is around 100°C, Arc1 is much more thermostable with an unfolding temperature of 114°C. However, only 12 out of 139 amino acids are different between the two sequences. In this study, one or a combination of amino acid(s) in Bac1 was/were substituted by a residue(s) found in Arc1 at the same position(s). The best mutant, which contained three amino acid substitutions (S108D, G116A and L120P substitutions), showed an unfolding temperature more than 10°C higher than that of Bac1. Furthermore, a combination of the other nine amino acid substitutions also led to improved thermostability of Bac1, although the effects of individual substitutions were small. Therefore, not only the sum of the contributions of individual amino acids, but also the synergistic effects of multiple amino acids are deeply involved in the stability of a hyperthermostable protein. Such insights will be helpful for future rational design of hyperthermostable proteins.

Catalytic mechanism of ancestral L-lysine oxidase assigned by sequence data mining

A large number of protein sequences are registered in public databases such as PubMed. Functionally uncharacterized enzymes are included in these databases, some of which likely have potential for industrial applications. However, assignment of the enzymes remained difficult tasks for now. In this study, we assigned a total of 28 original sequences to uncharacterized enzymes in the FAD-dependent oxidase family expressed in some species of bacteria including Chryseobacterium, Flavobacterium, and Pedobactor. Progenitor sequence of the assigned 28 sequences was generated by ancestral sequence reconstruction, and the generated sequence exhibited L-lysine oxidase activity; thus, we named the enzyme AncLLysO. Crystal structures of ligand-free and ligand-bound forms of AncLLysO were determined, indicating that the enzyme recognizes L-Lys by hydrogen bond formation with R76 and E383. The binding of L-Lys to AncLLysO induced dynamic structural change at a plug loop formed by residues 251 to 254. Biochemical assays of AncLLysO variants revealed the functional importance of these substrate recognition residues and the plug loop. R76A and E383D variants were also observed to lose their activity, and the kcat/Km value of G251P and Y253A mutations were approximately 800- to 1800-fold lower than that of AncLLysO, despite the indirect interaction of the substrates with the mutated residues. Taken together, our data demonstrate that combinational approaches to sequence classification from database and ancestral sequence reconstruction may be effective not only to find new enzymes using databases of unknown sequences but also to elucidate their functions.

Exogenous Enzymes as Zootechnical Additives in Animal Feed: A Review

Enzymes are widely used in the food industry. Their use as a supplement to the raw material for animal feed is a current research topic. Although there are several studies on the application of enzyme additives in the animal feed industry, it is necessary to search for new enzymes, as well as to utilize bioinformatics tools for the design of specific enzymes that work in certain environmental conditions and substrates. This will allow the improvement of the productive parameters in animals, reducing costs and making the processes more efficient. Technological needs have considered these catalysts as essential in many industrial sectors and research is constantly being carried out to optimize their use in those processes. This review describes the enzymes used in animal nutrition, their mode of action, their production and new sources of production as well as studies on different animal models to evaluate their effect on the productive performance intended for the production of animal feed.

Partial Consensus Design and Enhancement of Protein Function by Secondary-Structure-Guided Consensus Mutations

Consensus design (CD) is a representative sequence-based protein design method that enables the design of highly functional proteins by analyzing vast amounts of protein sequence data. This study proposes a partial consensus design (PCD) of a protein as a derivative approach of CD. The method replaces the target protein sequence with a consensus sequence in a secondary-structure-dependent manner (i.e., regionally dependent and divided into α-helix, β-sheet, and loop regions). In this study, we generated several artificial partial consensus l-threonine 3-dehydrogenases (PcTDHs) by PCD using the TDH from Cupriavidus necator (CnTDH) as a target protein. Structural and functional analysis of PcTDHs suggested that thermostability would be independently improved when consensus mutations are introduced into the loop region of TDHs. On the other hand, enzyme kinetic parameters (k_cat/K_m) and average productivity would be synergistically enhanced by changing the combination of the mutations-replacement of one region of CnTDH with a consensus sequence provided only negative effects, but the negative effects were nullified when the two regions were replaced simultaneously. Taken together, we propose the hypothesis that there are protein regions that encode individual protein properties, such as thermostability and activity, and that the introduction of consensus mutations into these regions could additively or synergistically modify their functions.

Revenant: a database of resurrected proteins

Revenant is a database of resurrected proteins coming from extinct organisms. Currently, it contains a manually curated collection of 84 resurrected proteins derived from bibliographic data. Each protein is extensively annotated, including structural, biochemical and biophysical information. Revenant contains a browse capability designed as a timeline from where the different proteins can be accessed. The oldest Revenant entries are between 4200 and 3500 million years ago, while the younger entries are between 8.8 and 6.3 million years ago. These proteins have been resurrected using computational tools called ancestral sequence reconstruction techniques combined with wet-laboratory synthesis and expression. Resurrected proteins are commonly used, with a noticeable increase during the past years, to explore and test different evolutionary hypotheses such as protein stability, to explore the origin of new functions, to get biochemical insights into past metabolisms and to explore specificity and promiscuous behaviour of ancient proteins.

A Reconstructed Common Ancestor of the Fatty Acid Photo-decarboxylase Clade Shows Photo-decarboxylation Activity and Increased Thermostability

Light-dependent enzymes are a rare type of biocatalyst with high potential for research and biotechnology. A recently discovered fatty acid photo-decarboxylase from Chlorella variabilis NC64A (CvFAP) converts fatty acids to the corresponding hydrocarbons only when irradiated with blue light (400 to 520 nm). To expand the available catalytic diversity for fatty acid decarboxylation, we reconstructed possible ancestral decarboxylases from a set of 12 extant sequences that were classified under the fatty acid decarboxylases clade within the glucose-methanol choline (GMC) oxidoreductase family. One of the resurrected enzymes (ANC1) showed activity in the decarboxylation of fatty acids, showing that the clade indeed contains several photo-decarboxylases. ANC1 has a 15 °C higher melting temperature (Tm ) than the extant CvFAP. Its production yielded 12-fold more protein than this wild type decarboxylase, which offers practical advantages for the biochemical investigation of this photoenzyme. Homology modelling revealed amino acid substitutions to more hydrophilic residues at the surface and shorter flexible loops compared to the wild type. Using ancestral sequence reconstruction, we have expanded the existing pool of confirmed fatty acid photo-decarboxylases, providing access to a more robust catalyst for further development via directed evolution.

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 lysosomal enzymes with increased activity harbor therapeutic potential for treatment of Hunter syndrome

We show the successful application of ancestral sequence reconstruction to enhance the activity of iduronate-2-sulfatase (IDS), thereby increasing its therapeutic potential for the treatment of Hunter syndrome—a lysosomal storage disease caused by impaired function of IDS. Current treatment, enzyme replacement therapy with recombinant human IDS, does not alleviate all symptoms, and an unmet medical need remains. We reconstructed putative ancestral sequences of mammalian IDS and compared them with extant IDS. Some ancestral variants displayed up to 2-fold higher activity than human IDS in in vitro assays and cleared more substrate in ex vivo experiments in patient fibroblasts. This could potentially allow for lower dosage or enhanced therapeutic effect in enzyme replacement therapy, thereby improving treatment outcomes and cost efficiency, as well as reducing treatment burden. In summary, we showed that ancestral sequence reconstruction can be applied to lysosomal enzymes that function in concert with modern enzymes and receptors in cells.

•

Reconstruction of ancestral lysosomal enzymes that function in complex cellular context

•

Ancestral iduronate-2-sulfatases with increased activity compared with the human enzyme

•

Increased clearance of substrate in patient fibroblasts indicates therapeutic potential

Engineering of Ancestors as a Tool to Elucidate Structure, Mechanism, and Specificity of Extant Terpene Cyclase

Structural information is crucial for understanding catalytic mechanisms and to guide enzyme engineering efforts of biocatalysts, such as terpene cyclases. However, low sequence similarity can impede homology modeling, and inherent protein instability presents challenges for structural studies. We hypothesized that X-ray crystallography of engineered thermostable ancestral enzymes can enable access to reliable homology models of extant biocatalysts. We have applied this concept in concert with molecular modeling and enzymatic assays to understand the structure activity relationship of spiroviolene synthase, a class I terpene cyclase, aiming to engineer its specificity. Engineering a surface patch in the reconstructed ancestor afforded a template structure for generation of a high-confidence homology model of the extant enzyme. On the basis of structural considerations, we designed and crystallized ancestral variants with single residue exchanges that exhibited tailored substrate specificity and preserved thermostability. We show how the two single amino acid alterations identified in the ancestral scaffold can be transferred to the extant enzyme, conferring a specificity switch that impacts the extant enzyme's specificity for formation of the diterpene spiroviolene over formation of sesquiterpenes hedycaryol and farnesol by up to 25-fold. This study emphasizes the value of ancestral sequence reconstruction combined with enzyme engineering as a versatile tool in chemical biology.

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.

Experimental Resurrection of Ancestral Mammalian CPEB3 Ribozymes Reveals Deep Functional Conservation

Self-cleaving ribozymes are genetic elements found in all domains of life, but their evolution remains poorly understood. A ribozyme located in the second intron of the cytoplasmic polyadenylation binding protein 3 gene (CPEB3) shows high sequence conservation in mammals, but little is known about the functional conservation of self-cleaving ribozyme activity across the mammalian tree of life or during the course of mammalian evolution. Here, we use a phylogenetic approach to design a mutational library and a deep sequencing assay to evaluate the in vitro self-cleavage activity of numerous extant and resurrected CPEB3 ribozymes that span over 100 My of mammalian evolution. We found that the predicted sequence at the divergence of placentals and marsupials is highly active, and this activity has been conserved in most lineages. A reduction in ribozyme activity appears to have occurred multiple different times throughout the mammalian tree of life. The in vitro activity data allow an evaluation of the predicted mutational pathways leading to extant ribozyme as well as the mutational landscape surrounding these ribozymes. The results demonstrate that in addition to sequence conservation, the self-cleavage activity of the CPEB3 ribozyme has persisted over millions of years of mammalian evolution.

De novo proteins from random sequences through in vitro evolution

Natural proteins are the result of billions of years of evolution. The earliest predecessors of today's proteins are believed to have emerged from random polypeptides. While we have no means to determine how this process exactly happened, there is great interest in understanding how it reasonably could have happened. We are reviewing how researchers have utilized in vitro selection and molecular evolution methods to investigate plausible scenarios for the emergence of early functional proteins. The studies range from analyzing general properties and structural features of unevolved random polypeptides to isolating de novo proteins with specific functions from synthetic randomized sequence libraries or generating novel proteins by combining evolution with rational design. While the results are exciting, more work is needed to fully unravel the mechanisms that seeded protein-dominated biology.

Reconstructed evolutionary history of the yeast septins Cdc11 and Shs1

Septins are GTP-binding proteins conserved across metazoans. They can polymerize into extended filaments and, hence, are considered a component of the cytoskeleton. The number of individual septins varies across the tree of life—yeast (Saccharomyces cerevisiae) has seven distinct subunits, a nematode (Caenorhabditis elegans) has two, and humans have 13. However, the overall geometric unit (an apolar hetero-octameric protomer and filaments assembled there from) has been conserved. To understand septin evolutionary variation, we focused on a related pair of yeast subunits (Cdc11 and Shs1) that appear to have arisen from gene duplication within the fungal clade. Either Cdc11 or Shs1 occupies the terminal position within a hetero-octamer, yet Cdc11 is essential for septin function and cell viability, whereas Shs1 is not. To discern the molecular basis of this divergence, we utilized ancestral gene reconstruction to predict, synthesize, and experimentally examine the most recent common ancestor (“Anc.11-S”) of Cdc11 and Shs1. Anc.11-S was able to occupy the terminal position within an octamer, just like the modern subunits. Although Anc.11-S supplied many of the known functions of Cdc11, it was unable to replace the distinct function(s) of Shs1. To further evaluate the history of Shs1, additional intermediates along a proposed trajectory from Anc.11-S to yeast Shs1 were generated and tested. We demonstrate that multiple events contributed to the current properties of Shs1: (1) loss of Shs1–Shs1 self-association early after duplication, (2) co-evolution of heterotypic Cdc11–Shs1 interaction between neighboring hetero-octamers, and (3) eventual repurposing and acquisition of novel function(s) for its C-terminal extension domain. Thus, a pair of duplicated proteins, despite constraints imposed by assembly into a highly conserved multi-subunit structure, could evolve new functionality via a complex evolutionary pathway.