Ancestral Sequence Reconstruction Research Articles

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.

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.

Adaptive evolution of scn4aa in Takifugu and Tetraodon

Tetrodotoxin (TTX) is a potent neurotoxin firstly discovered in the ovary of pufferfish. The genetic basis of voltage-gated sodium channel resistance to TTX has been widely studied, but it remains unclear in the evolution history of voltage-gated sodium channel resistance to TTX in pufferfish with different TTX concentrations. In this study, six scn4aa coding sequences of pufferfish were firstly cloned and sequenced, then used to investigate the adaptive evolution of scn4aa associated with TTX concentration and reconstruct ancestral sequences with seven scn4aa of other fishes. The result of CODEML (codon substitution model) program from the PAML (phylogenetic analysis by maximum likelihood) package shows only in the genus of Takifugu, which contains TTX highly in the liver, under positive selection. The result also indicates that three of four positively selected sites are located in the intracellular regions, which may compensate for normal function. The ancestral sequence reconstruction may suggest that the replacements providing weak toxin resistance might have appeared first in scn4aa, then the genus Takifugu evolved stronger resistance to TTX later. These results contribute to the explanation of the evolutional history of voltage-gated sodium channel resistance to TTX in pufferfish.

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.

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.

Evolution of increased complexity and specificity at the dawn of form I Rubiscos.

The evolution of ribulose-1,5-bisphosphate carboxylase/oxygenases (Rubiscos) that discriminate strongly between their substrate carbon dioxide and the undesired side substrate dioxygen was an important event for photosynthetic organisms adapting to an oxygenated environment. We use ancestral sequence reconstruction to recapitulate this event. We show that Rubisco increased its specificity and carboxylation efficiency through the gain of an accessory subunit before atmospheric oxygen was present. Using structural and biochemical approaches, we retrace how this subunit was gained and became essential. Our work illuminates the emergence of an adaptation to rising ambient oxygen levels, provides a template for investigating the function of interactions that have remained elusive because of their essentiality, and sheds light on the determinants of specificity in Rubisco.

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 (Tm) 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.

Saturation mutagenesis of a predicted ancestral Syk-family kinase.

Many tyrosine kinases cannot be expressed readily in Escherichia coli, limiting facile production of these proteins for biochemical experiments. We used ancestral sequence reconstruction to generate a spleen tyrosine kinase (Syk) variant that can be expressed in bacteria and purified in soluble form, unlike the human members of this family (Syk and zeta‐chain‐associated protein kinase of 70 kDa [ZAP‐70]). The catalytic activity, substrate specificity, and regulation by phosphorylation of this Syk variant are similar to the corresponding properties of human Syk and ZAP‐70. Taking advantage of the ability to express this novel Syk‐family kinase in bacteria, we developed a two‐hybrid assay that couples the growth of E. coli in the presence of an antibiotic to successful phosphorylation of a bait peptide by the kinase. Using this assay, we screened a site‐saturation mutagenesis library of the kinase domain of this reconstructed Syk‐family kinase. Sites of loss‐of‐function mutations identified in the screen correlate well with residues established previously as critical to function and/or structure in protein kinases. We also identified activating mutations in the regulatory hydrophobic spine and activation loop, which are within key motifs involved in kinase regulation. Strikingly, one mutation in an ancestral Syk‐family variant increases the soluble expression of the protein by 75‐fold. Thus, through ancestral sequence reconstruction followed by deep mutational scanning, we have generated Syk‐family kinase variants that can be expressed in bacteria with very high yield.

Characterization of ancestral Fe/Mn superoxide dismutases indicates their cambialistic origin.

Superoxide dismutases (SODs) are critical metalloenzymes mitigating the damages of the modern oxygenated world. However, the emergence of one family of SODs, the Fe/Mn SOD, has been recurrently proposed to predate the great oxygenation event (GOE). This ancient family lacks metal binding selectivity, but displays strong catalytic selectivity. Therefore, some homologues would only be active when bound to Fe or Mn, although others, dubbed cambialistic, would function when loaded with either ion. This posed the longstanding question about the identity of the cognate metal ion of the first SODs to emerge. In this work, we utilize ancestral sequence reconstruction techniques to infer the earliest SODs. We show that the “ancestors” are active in vivo and in vitro. Further, we test their metal specificity and demonstrate that they are cambialistic in nature. Our findings shed light on how the predicted Last Common Universal Ancestor was capable of dealing with decomposition of the superoxide anion, and the early relationship between life, oxygen, and metal ion availability.

Phylogenetic Analysis of Anti-CRISPR and Member Addition in the Families.

CRISPR-Cas is a widespread anti-viral adaptive immune system in the microorganisms. Viruses living in bacteria or some phages carry anti-CRISPR proteins to evade immunity by CRISPR-Cas. The anti-CRISPR proteins are prevalent in phages capable of lying dormant in a CRISPR-carrying host, while their orthologs frequently found in virulent phages. Here, we propose a probabilistic strategy of ancestral sequence reconstruction (ASR) and Hidden Markov Model (HMM) profile search to fish out sequences of anti-CRISPR proteins from environmental metagenomic, human microbiome metagenomic, human microbiome reference genome, and NCBI's non-redundant databases. Our results revealed that the metagenome database dark matter might contain anti-CRISPR encoding genes.

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 Tm 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 Identifies Structural Changes Underlying the Evolution of Ideonella sakaiensis PETase and Variants with Improved Stability and Activity.

The improved production, recycling, and removal of plastic waste, such as polyethylene terephthalate (PET), are pressing environmental and economic issues for society. Biocatalytic (enzymatic) PET depolymerization is potentially a sustainable, low-energy solution to PET recycling, especially when compared with current disposal methods such as landfills, incineration, or gasification. IsPETase has been extensively studied for its use in PET depolymerization; however, its evolution from cutinases is not fully understood, and most engineering studies have neglected the majority of the available sequence space remote from the active site. In this study, ancestral protein reconstruction (ASR) has been used to trace the evolutionary trajectory from ancient serine hydrolases to IsPETase, while ASR and the related design approach, protein repair one-stop shop, were used to identify enzyme variants with improved activity and stability. Kinetic and structural characterization of these variants reveals new insights into the evolution of PETase activity and the role of second-shell mutations around the active site. Among the designed and reconstructed variants, we identified several with melting points 20 °C higher than that of IsPETase and two variants with significantly higher catalytic activity.

ARPIP: Ancestral Sequence Reconstruction with Insertions and Deletions under the Poisson Indel Process.

Modern phylogenetic methods allow inference of ancestral molecular sequences given an alignment and phylogeny relating present-day sequences. This provides insight into the evolutionary history of molecules, helping to understand gene function and to study biological processes such as adaptation and convergent evolution across a variety of applications. Here, we propose a dynamic programming algorithm for fast joint likelihood-based reconstruction of ancestral sequences under the Poisson Indel Process (PIP). Unlike previous approaches, our method, named ARPIP, enables the reconstruction with insertions and deletions based on an explicit indel model. Consequently, inferred indel events have an explicit biological interpretation. Likelihood computation is achieved in linear time with respect to the number of sequences. Our method consists of two steps, namely finding the most probable indel points and reconstructing ancestral sequences. First, we find the most likely indel points and prune the phylogeny to reflect the insertion and deletion events per site. Second, we infer the ancestral states on the pruned subtree in a manner similar to FastML. We applied ARPIP (Ancestral Reconstruction under PIP) on simulated data sets and on real data from the Betacoronavirus genus. ARPIP reconstructs both the indel events and substitutions with a high degree of accuracy. Our method fares well when compared to established state-of-the-art methods such as FastML and PAML. Moreover, the method can be extended to explore both optimal and suboptimal reconstructions, include rate heterogeneity through time and more. We believe it will expand the range of novel applications of ancestral sequence reconstruction. [Ancestral sequences; dynamic programming; evolutionary stochastic process; indel; joint ancestral sequence reconstruction; maximum likelihood; Poisson Indel Process; phylogeny; SARS-CoV.].

Single mutations toggle the substrate selectivity of multifunctional Camptotheca secologanic acid synthases

Terpene indole alkaloids (TIAs) are plant-derived specialized metabolites with widespread use in medicine. Species-specific pathways derive various TIAs from common intermediates, strictosidine or strictosidinic acid, produced by coupling tryptamine with secologanin or secologanic acid. The penultimate reaction in this pathway is catalyzed by either secologanin synthase (SLS) or secologanic acid synthase (SLAS) according to whether plants produce secologanin from loganin or secologanic acid from loganic acid. Previous work has identified SLSs and SLASs from different species, but the determinants of selectivity remain unclear. Here, combining molecular modeling, ancestral sequence reconstruction, and biochemical methodologies, we identified key residues that toggle SLS and SLAS selectivity in two CYP72A (cytochrome P450) subfamily enzymes from Camptotheca acuminata. We found that the positions of foremost importance are in substrate recognition sequence 1 (SRS1), where mutations to either of two adjacent histidine residues switched selectivity; His131Phe selects for and increases secologanin production whereas His132Asp selects for secologanic acid production. Furthermore, a change in SRS3 in the predicted substrate entry channel (Arg/Lys270Thr) and another in SRS4 at the start of the I-helix (Ser324Glu) decreased enzyme activity toward either substrate. We propose that the Camptotheca SLASs have maintained the broadened activities found in a common asterid ancestor, even as the Camptotheca lineage lost its ability to produce loganin while the campanulid and lamiid lineages specialized to produce secologanin by acquiring mutations in SRS1. The identification here of the residues essential for the broad substrate scope of SLASs presents opportunities for more tailored heterologous production of TIAs.

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.

Ancient whale rhodopsin reconstructs dim-light vision over a major evolutionary transition: Implications for ancestral diving behavior

Cetaceans are fully aquatic mammals that descended from terrestrial ancestors, an iconic evolutionary transition characterized by adaptations for underwater foraging via breath-hold diving. Although the evolutionary history of this specialized behavior is challenging to reconstruct, coevolving sensory systems may offer valuable clues. The dim-light visual pigment, rhodopsin, which initiates phototransduction in the rod photoreceptors of the eye, has provided insight into the visual ecology of depth in several aquatic vertebrate lineages. Here, we use ancestral sequence reconstruction and protein resurrection experiments to quantify light-activation metrics in rhodopsin pigments from ancestors bracketing the cetacean terrestrial-to-aquatic transition. By comparing multiple reconstruction methods on a broadly sampled cetartiodactyl species tree, we generated highly robust ancestral sequence estimates. Our experimental results provide direct support for a blue-shift in spectral sensitivity along the branch separating cetaceans from terrestrial relatives. This blue-shift was 14 nm, resulting in a deep-sea signature (λmax = 486 nm) similar to many mesopelagic-dwelling fish. We also discovered that the decay rates of light-activated rhodopsin increased in ancestral cetaceans, which may indicate an accelerated dark adaptation response typical of deeper-diving mammals. Because slow decay rates are thought to help sequester cytotoxic photoproducts, this surprising result could reflect an ecological trade-off between rod photoprotection and dark adaptation. Taken together, these ancestral shifts in rhodopsin function suggest that some of the first fully aquatic cetaceans could dive into the mesopelagic zone (>200 m). Moreover, our reconstructions indicate that this behavior arose before the divergence of toothed and baleen whales.

Ancestral sequence reconstruction of the CYP711 family reveals functional divergence in strigolactone biosynthetic enzymes associated with gene duplication events in monocot grasses.

Summary The strigolactone (SL) class of phytohormones shows broad chemical diversity, the functional importance of which remains to be fully elucidated, along with the enzymes responsible for the diversification of the SL structure. Here we explore the functional evolution of the highly conserved CYP711A P450 family, members of which catalyze several key monooxygenation reactions in the strigolactone pathway.Ancestral sequence reconstruction was utilized to infer ancestral CYP711A sequences based on a comprehensive set of extant CYP711 sequences. Eleven ancestral enzymes, corresponding to key points in the CYP711A phylogenetic tree, were resurrected and their activity was characterized towards the native substrate carlactone and the pure enantiomers of the synthetic strigolactone analogue, GR24.The ancestral and extant CYP711As tested accepted GR24 as a substrate and catalyzed several diversifying oxidation reactions on the structure. Evidence was obtained for functional divergence in the CYP711A family. The monocot group 3 ancestor, arising from gene duplication events within monocot grasses, showed both increased catalytic activity towards GR24 and high stereoselectivity towards the GR24 isomer resembling strigol‐type SLs.These results are consistent with a role for CYP711As in strigolactone diversification in early land plants, which may have extended to the diversification of strigol‐type SLs.

Molecular Evolution of Far-Red Light-Acclimated Photosystem II.

Cyanobacteria are major contributors to global carbon fixation and primarily use visible light (400−700 nm) to drive oxygenic photosynthesis. When shifted into environments where visible light is attenuated, a small, but highly diverse and widespread number of cyanobacteria can express modified pigments and paralogous versions of photosystem subunits and phycobiliproteins that confer far-red light (FRL) absorbance (700−800 nm), a process termed far-red light photoacclimation, or FaRLiP. During FaRLiP, alternate photosystem II (PSII) subunits enable the complex to bind chlorophylls d and f, which absorb at lower energy than chlorophyll a but still support water oxidation. How the FaRLiP response arose remains poorly studied. Here, we report ancestral sequence reconstruction and structure-based molecular evolutionary studies of the FRL-specific subunits of FRL-PSII. We show that the duplications leading to the origin of two PsbA (D1) paralogs required to make chlorophyll f and to bind chlorophyll d in water-splitting FRL-PSII are likely the first to have occurred prior to the diversification of extant cyanobacteria. These duplications were followed by those leading to alternative PsbC (CP43) and PsbD (D2) subunits, occurring early during the diversification of cyanobacteria, and culminating with those leading to PsbB (CP47) and PsbH paralogs coincident with the radiation of the major groups. We show that the origin of FRL-PSII required the accumulation of a relatively small number of amino acid changes and that the ancestral FRL-PSII likely contained a chlorophyll d molecule in the electron transfer chain, two chlorophyll f molecules in the antenna subunits at equivalent positions, and three chlorophyll a molecules whose site energies were altered. The results suggest a minimal model for engineering far-red light absorbance into plant PSII for biotechnological applications.

Ancestral reconstruction of duplicated signaling proteins reveals the evolution of signaling specificity.

Gene duplication is crucial to generating novel signaling pathways during evolution. However, it remains unclear how the redundant proteins produced by gene duplication ultimately acquire new interaction specificities to establish insulated paralogous signaling pathways. Here, we used ancestral sequence reconstruction to resurrect and characterize a bacterial two-component signaling system that duplicated in α-proteobacteria. We determined the interaction specificities of the signaling proteins that existed before and immediately after this duplication event and then identified key mutations responsible for establishing specificity in the two systems. Just three mutations, in only two of the four interacting proteins, were sufficient to establish specificity of the extant systems. Some of these mutations weakened interactions between paralogous systems to limit crosstalk. However, others strengthened interactions within a system, indicating that the ancestral interaction, although functional, had the potential to be strengthened. Our work suggests that protein-protein interactions with such latent potential may be highly amenable to duplication and divergence.

Comparative Analysis and Ancestral Sequence Reconstruction of Bacterial Sortase Family Proteins Generates Functional Ancestral Mutants with Different Sequence Specificities

Gram-positive bacteria are some of the earliest known life forms, diverging from gram-negative bacteria 2 billion years ago. These organisms utilize sortase enzymes to attach proteins to their peptidoglycan cell wall, a structural feature that distinguishes the two types of bacteria. The transpeptidase activity of sortases make them an important tool in protein engineering applications, e.g., in sortase-mediated ligations or sortagging. However, due to relatively low catalytic efficiency, there are ongoing efforts to create better sortase variants for these uses. Here, we use bioinformatics tools, principal component analysis and ancestral sequence reconstruction, in combination with protein biochemistry, to analyze natural sequence variation in these enzymes. Principal component analysis on the sortase superfamily distinguishes previously described classes and identifies regions of relatively high sequence variation in structurally-conserved loops within each sortase family, including those near the active site. Using ancestral sequence reconstruction, we determined sequences of ancestral Staphylococcus and Streptococcus Class A sortase proteins. Enzyme assays revealed that the ancestral Streptococcus enzyme is relatively active and shares similar sequence variation with other Class A Streptococcus sortases. Taken together, we highlight how natural sequence variation can be utilized to investigate this important protein family, arguing that these and similar techniques may be used to discover or design sortases with increased catalytic efficiency and/or selectivity for sortase-mediated ligation experiments.