Ancestral Sequence Research Articles

Regulatory transposable elements in the encyclopedia of DNA elements

Transposable elements (TEs) comprise ~50% of our genome, but knowledge of how TEs affect genome evolution remains incomplete. Leveraging ENCODE4 data, we provide the most comprehensive study to date of TE contributions to the regulatory genome. We find 236,181 (~25%) human candidate cis-regulatory elements (cCREs) are TE-derived, with over 90% lineage-specific since the human-mouse split, accounting for 8–36% of lineage-specific cCREs. Except for SINEs, cCRE-associated transcription factor (TF) motifs in TEs are derived from ancestral TE sequence more than expected by chance. We show that TEs may adopt similar regulatory activities of elements near their integration site. Since human-mouse divergence, TEs have contributed 3–56% of TF binding site turnover events across 30 examined TFs. Finally, TE-derived cCREs are similar to non-TE cCREs in terms of MPRA activity and GWAS variant enrichment. Overall, our results substantiate the notion that TEs have played an important role in shaping the human regulatory genome.

Enhancing the Catalytic Activity of Geranylgeranyl Diphosphate Synthase through Ancestral Sequence Reconstruction and Semirational Design.

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

Engineering an Ancestral Glycosyltransferase for Biosynthesis of 2-Phenylethyl-β-d-Glucopyranoside and Salidroside.

Phenylethanoid glycosides (PhGs) are naturally occurring glycosides derived from plants with various biological activities. Glycosyltransferases catalyze the production of PhGs from phenylethanols via a transglycosylation reaction. The low activity and stability of glycosyltransferase limit its industrial application. An ancestral glycosyltransferase, UGTAn85, with heat resistance, alkali resistance, and high stability was resurrected using ancestral sequence reconstruction technology. This enzyme can efficiently convert phenylethanols to PhGs. The optimal reaction temperature and pH for UGTAn85 were found to be 70 °C and pH 10.0, respectively. This study employed a combination of structure-guided rational design and co-evolution analysis to enhance its catalytic activity. Potential mutation sites were identified through computer-aided design, including homology modeling, molecular docking, Rosetta dock design, molecular dynamics simulation, and co-evolution analysis. By targeted mutagenesis, the UGTAn85 mutant Q23E/N65D exhibited a 2.2-fold increase in enzyme activity (11.85 U/mg) and elevated affinity (Km = 0.11 mM) for 2-phenylethanol compared to UGTAn85. Following a fed-batch reaction, 36.16 g/L 2-phenylethyl-β-d-glucopyranoside and 51.49 g/L salidroside could be produced within 24 h, respectively. The findings in this study provide a new perspective on enhancing the stability and activity of glycosyltransferases, as well as a potential biocatalyst for the industrial production of PhGs.

Engineering the Substrate Specificity of UDP‐Glycosyltransferases for Synthesizing Triterpenoid Glycosides with a Linear Trisaccharide as Aided by Ancestral Sequence Reconstruction

Engineering the Substrate Specificity of UDP‐Glycosyltransferases for Synthesizing Triterpenoid Glycosides with a Linear Trisaccharide as Aided by Ancestral Sequence Reconstruction

Stable monomers in the ancestral sequence reconstruction of the last opisthokont common ancestor of dimeric triosephosphate isomerase.

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

Characterization, structural, and evolutionary analysis of an extremophilic GH5 endoglucanase from Bacillus sp. G131: Insights from ancestral sequence reconstruction

Nature has developed extremozymes that catalyze complex reaction processes in extreme environmental conditions. Accordingly, a combined approach consisting of extremozyme screening, ancestral sequence resurrection (ASR), and molecular dynamic simulation was utilized to construct a developed endoglucanase. The primary experimental and in-silico data led to the prediction of a hypothetical sequence of endoglucanase (EG5-G131) using Bacillus sp. G131 confirmed by amplification and sequencing. EG5-G131 exhibited noticeable stability in a broad-pH range, several detergents, organic solvents, and temperatures up to 80 °C. The molecular weight, Vmax, and Km of the purified endoglucanase were estimated to be 36 kDa, 4.32 μmol/min, and 23.62 mg/ml, respectively. The calculated thermodynamic parameters for EG5-G131 confirmed its intrinsic thermostability. Computational analysis revealed Glu142 and Glu230 as active-site residues of the enzyme. Furthermore, the enzyme remained bound to cellotetraose at 298 K, 333 K, 343 K, and 353 K for 300 ns, consistent with our experimental data. ASR of EG5-G131 led to the introduction of ancestral ANC204 and ANC205, which show similar thermodynamic characteristics with the last Firmicute common ancestor. Finally, truncating loops from the N-terminal of two sequences created two variants with desirable thermal stability, suggesting the evolutionary deciphering of the functional domain of the GH5 family in Bacillus sp. G131.

Unique myoglobin adaptation to endothermy and flight since the origin of birds.

Myoglobin (Mb) mediates oxygen diffusion and storage in muscle tissue and thus is important for the energy utilization and activity of animals. Birds generally have a high body temperature, and most species also possess the capability of powered flight. Both of these require high levels of aerobic metabolism. Within endothermic mammals, bats also independently evolved flight. Although the functional evolution of myoglobins in deep-diving amniote vertebrates has been well-studied, the functional evolution of myoglobin since the origins of both birds and bats is unclear. Here, with Mb-coding sequences from >200 extant amniote species, we reconstructed ancestral sequences to estimate the functional properties of myoglobin through amniote evolution. A dramatic change in net surface charge on myoglobin occurred during the origin of Aves, which might have been driven by positively selected amino acid substitutions that occurred on the lineage leading to all birds. However, in bats, no change in net surface charge occurred and instead, the Mb genes show evidence of strong purifying selection. The increased net surface charge on bird myoglobins implies an adaptation to flight-related endothermic and higher body temperatures, possibly by reducing harmful protein aggregations. Different from the findings of net surface charge, myoglobins of extant birds show lower stability compared with other amniotes, which probably accelerates the rate of oxygen utilization in muscles. In bats and other mammals, higher stability of Mb may be an alternative pathway for adaptation to endothermy, indicating divergent evolution of myoglobin in birds and bats.

The smallest angiosperm genomes may be the price for effective traps of bladderworts.

Species of the carnivorous family Lentibulariaceae exhibit the smallest genomes in flowering plants. We explored the hypothesis that their minute genomes result from the unique mitochondrial cytochrome c oxidase (COX) mutation. The mutation may boost mitochondrial efficiency, which is especially useful for suction-bladder traps of Utricularia, but also increase DNA-damaging reactive oxygen species, leading to genome shrinkage through deletion-biased DNA repair. We aimed to explore this mutation's impact on genome size, providing insights into genetic mutation roles in plant genome evolution under environmental pressures. We compiled and measured genome and mean chromosome sizes for 127 and 67 species, respectively, representing all three genera (Genlisea, Pinguicula, and Utricularia) of Lentibulariaceae. We also isolated and analyzed COX sequences to detect the mutation. Through phylogenetic regressions and Ornstein-Uhlenbeck models of trait evolution, we assessed the impact of the COX mutation on the genome and chromosome sizes across the family. Our findings reveal significant correlations between the COX mutations and smaller genome and chromosome sizes. Specifically, species carrying the ancestral COX sequence exhibited larger genomes and chromosomes than those with the mutation. This evidence supports the notion that the COX mutation contributes to genome downsizing, with statistical analyses confirming a directional evolution towards smaller genomes in species harboring these mutations. Our study confirms that the COX mutation in Lentibulariaceae is associated with genome downsizing, likely driven by increased reactive oxygen species production and subsequent DNA damage requiring deletion-biased repair mechanisms. While boosting mitochondrial energy output, this genetic mutation compromises genome integrity and may potentially affect recombination rates, illustrating a complex trade-off between evolutionary advantages and disadvantages. Our results highlight the intricate processes by which genetic mutations and environmental pressures shape genome size evolution in carnivorous plants.

Ancestral Sequence Reconstruction as a Tool to Detect and Study De Novo Gene Emergence.

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

Pseudo-perplexity in One Fell Swoop for Protein Fitness Estimation.

Protein language models trained on the masked language modeling objective learn to predict the identity of hidden amino acid residues within a sequence using the remaining observable sequence as context. They do so by embedding the residues into a high dimensional space that encapsulates the relevant contextual cues. These embedding vectors serve as an informative context-sensitive representation that not only aids with the defined training objective, but can also be used for other tasks by downstream models. We propose a scheme to use the embeddings of an unmasked sequence to estimate the corresponding masked probability vectors for all the positions in a single forward pass through the language model. This One Fell Swoop (OFS) approach allows us to efficiently estimate the pseudo-perplexity of the sequence, a measure of the model's uncertainty in its predictions, that can also serve as a fitness estimate. We find that ESM2 OFS pseudo-perplexity performs nearly as well as the true pseudo-perplexity at fitness estimation, and more notably it defines a new state of the art on the ProteinGym Indels benchmark. The strong performance of the fitness measure prompted us to investigate if it could be used to detect the elevated stability reported in reconstructed ancestral sequences. We find that this measure ranks ancestral reconstructions as more fit than extant sequences. Finally, we show that the computational efficiency of the technique allows for the use of Monte Carlo methods that can rapidly explore functional sequence space.

Following the Evolutionary Paths of Dscam1 Proteins toward Highly Specific Homophilic Interactions.

Many adhesion proteins, evolutionarily related through gene duplication, exhibit distinct and precise interaction preferences and affinities crucial for cell patterning. Yet, the evolutionary paths by which these proteins acquire new specificities and prevent cross-interactions within their family members remain unknown. To bridge this gap, this study focuses on Drosophila Down syndrome cell adhesion molecule-1 (Dscam1) proteins, which are cell adhesion proteins that have undergone extensive gene duplication. Dscam1 evolved under strong selective pressure to achieve strict homophilic recognition, essential for neuronal self-avoidance and patterning. Through a combination of phylogenetic analyses, ancestral sequence reconstruction, and cell aggregation assays, we studied the evolutionary trajectory of Dscam1 exon 4 across various insect lineages. We demonstrated that recent Dscam1 duplications in the mosquito lineage bind with strict homophilic specificities without any cross-interactions. We found that ancestral and intermediate Dscam1 isoforms maintained their homophilic binding capabilities, with some intermediate isoforms also engaging in promiscuous interactions with other paralogs. Our results highlight the robust selective pressure for homophilic specificity integral to the Dscam1 function within the process of neuronal self-avoidance. Importantly, our study suggests that the path to achieving such selective specificity does not introduce disruptive mutations that prevent self-binding but includes evolutionary intermediates that demonstrate promiscuous heterophilic interactions. Overall, these results offer insights into evolutionary strategies that underlie adhesion protein interaction specificities.

Enriching barcoding markers in environmental samples utilizing a phylogenetic probe design: Insights from mock communities

AbstractHybridization capture is an emerging method making use of short oligonucleotide baits to enrich DNA libraries for genomic fragments of specific organisms thus enabling detection of their presence in environmental samples. Although it offers a primer‐independent alternative to metabarcoding, little empirical work has been dedicated to characterizing the underlying biases and coupled implications for biological interpretation. Moreover, few published bioinformatic pipelines are available for designing polynucleotide capture baits from a reference sequence collection. We designed RNA‐baits specifically targeting two chloroplast barcoding genes matK and rbcL to reveal the plant taxonomic diversity present in a given environmental sample. Our approach leverages the sensitivity of hybridization capture and the capacity of high‐throughput DNA sequencing instruments. It builds on a new and universal method based on ancestral sequence reconstruction, ultimately limiting the number of bait‐probes required and reducing experimental costs, while accessing high levels of taxonomic diversity. Our bait‐set selectively targets four main plant orders (Fagales, Pinales, Asterales, and Poales), representing ~18% of all described vascular plants. This is achieved through the use of only 4084 baits, each 80 nucleotides in length (80‐mer), capturing ~1.0–1.6 k nucleotide sequences from each taxon. Tests on mock communities revealed important factors influencing capture efficiency and relative abundance estimates, including GC‐content, the overall target length per taxa, and the bait density and mean number of mismatches to the bait sequence. Our results show that hybridization capture, like metabarcoding, requires caution when interpreting results quantitatively within (paleo)‐ecological studies. Biases detected in this work have the potential to be mitigated with bait designs that avoid extreme base compositional biases and balancing bait targets across taxa. However, we strongly recommend the use of mock communities and read simulations to quantify the accuracy of taxonomic representation when using new bait designs.

Optimal phylogenetic reconstruction of insertion and deletion events.

Insertions and deletions (indels) influence the genetic code in fundamentally distinct ways from substitutions, significantly impacting gene product structure and function. Despite their influence, the evolutionary history of indels is often neglected in phylogenetic tree inference and ancestral sequence reconstruction, hindering efforts to comprehend biological diversity determinants and engineer variants for medical and industrial applications. We frame determining the optimal history of indel events as a single Mixed-Integer Programming (MIP) problem, across all branch points in a phylogenetic tree adhering to topological constraints, and all sites implied by a given set of aligned, extant sequences. By disentangling the impact on ancestral sequences at each branch point, this approach identifies the minimal indel events that jointly explain the diversity in sequences mapped to the tips of that tree. MIP can recover alternate optimal indel histories, if available. We evaluated MIP for indel inference on a dataset comprising 15 real phylogenetic trees associated with protein families ranging from 165 to 2000 extant sequences, and on 60 synthetic trees at comparable scales of data and reflecting realistic rates of mutation. Across relevant metrics, MIP outperformed alternative parsimony-based approaches and reported the fewest indel events, on par or below their occurrence in synthetic datasets. MIP offers a rational justification for indel patterns in extant sequences; importantly, it uniquely identifies global optima on complex protein data sets without making unrealistic assumptions of independence or evolutionary underpinnings, promising a deeper understanding of molecular evolution and aiding novel protein design. The implementation is available via GitHub at https://github.com/santule/indelmip.

Engineering a Robust UDP-Glucose Pyrophosphorylase for Enhanced Biocatalytic Synthesis via ProteinMPNN and Ancestral Sequence Reconstruction.

UDP-glucose is a key metabolite in carbohydrate metabolism and plays a vital role in glycosyl transfer reactions. Its significance spans across the food and agricultural industries. This study focuses on UDP-glucose synthesis via multienzyme catalysis using dextrin, incorporating UTP production and ATP regeneration modules to reduce costs. To address thermal stability limitations of the key UDP-glucose pyrophosphorylase (UGP), a deep learning-based protein sequence design approach and ancestral sequence reconstruction are employed to engineer a thermally stable UGP variant. The engineered UGP variant is significantly 500-fold more thermally stable at 60 °C and has a half-life of 49.8 h compared to the wild-type enzyme. MD simulations and umbrella sampling calculations provide insights into the mechanism behind the enhanced thermal stability. Experimental validation demonstrates that the engineered UGP variant can produce 52.6 mM UDP-glucose within 6 h in an in vitro cascade reaction. This study offers practical insights for efficient UDP-glucose synthesis methods.

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

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

Evidence that the cold- and menthol-sensing functions of the human TRPM8 channel evolved separately.

Transient receptor potential melastatin 8 (TRPM8) is a temperature- and menthol-sensitive ion channel that contributes to diverse physiological roles, including cold sensing and pain perception. Clinical trials targeting TRPM8 have faced repeated setbacks predominantly due to the knowledge gap in unraveling the molecular underpinnings governing polymodal activation. A better understanding of the molecular foundations between the TRPM8 activation modes may aid the development of mode-specific, thermal-neutral therapies. Ancestral sequence reconstruction was used to explore the origins of TRPM8 activation modes. By resurrecting key TRPM8 nodes along the human evolutionary trajectory, we gained valuable insights into the trafficking, stability, and function of these ancestral forms. Notably, this approach unveiled the differential emergence of cold and menthol sensitivity over evolutionary time, providing a fresh perspective on complex polymodal behavior. These studies provide a paradigm for understanding polymodal behavior in TRPM8 and other proteins with the potential to enhance our understanding of sensory receptor biology and pave the way for innovative therapeutic interventions.

COATi: Statistical Pairwise Alignment of Protein-Coding Sequences.

Sequence alignment is an essential method in bioinformatics and the basis of many analyses, including phylogenetic inference, ancestral sequence reconstruction, and gene annotation. Sequencing artifacts and errors made during genome assembly, such as abiological frameshifts and incorrect early stop codons, can impact downstream analyses leading to erroneous conclusions in comparative and functional genomic studies. More significantly, while indels can occur both within and between codons in natural sequences, most amino-acid- and codon-based aligners assume that indels only occur between codons. This mismatch between biology and alignment algorithms produces suboptimal alignments and errors in downstream analyses. To address these issues, we present COATi, a statistical, codon-aware pairwise aligner that supports complex insertion-deletion models and can handle artifacts present in genomic data. COATi allows users to reduce the amount of discarded data while generating more accurate sequence alignments. COATi can infer indels both within and between codons, leading to improved sequence alignments. We applied COATi to a dataset containing orthologous protein-coding sequences from humans and gorillas and conclude that 41% of indels occurred between codons, agreeing with previous work in other species. We also applied COATi to semiempirical benchmark alignments and find that it outperforms several popular alignment programs on several measures of alignment quality and accuracy.

Please Mind the Gap: Indel-Aware Parsimony for Fast and Accurate Ancestral Sequence Reconstruction and Multiple Sequence Alignment Including Long Indels.

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

Structure and Activity of Reconstructed Pseudo-Ancestral Cyclotides.

Cyclotides are cyclic peptides that are promising scaffolds for the design of drug candidates and chemical tools. However, despite there being hundreds of reported cyclotides, drug design studies have commonly focussed on a select few prototypic examples. Here, we explored whether ancestral sequence reconstruction could be used to generate new cyclotides for further optimization. We show that the reconstructed 'pseudo-ancestral' sequences, named Ancy-m (for the ancestral cyclotide of the Möbius sub-family) and Ancy-b (for the bracelet sub-family), have well-defined structures like their extant members, comprising the core structural feature of a cyclic cystine knot. This motif underpins efforts to re-engineer cyclotides for agrochemical and therapeutic applications. We further show that the reconstructed sequences are resistant to temperatures approaching boiling, bind to phosphatidyl-ethanolamine lipid bilayers at micromolar affinity, and inhibit the growth of insect cells at inhibitory concentrations in the micromolar range. Interestingly, the Ancy-b cyclotide had a higher oxidative folding yield than its comparator cyclotide cyO2, which belongs to the bracelet cyclotide subfamily known to be notoriously difficult to fold. Overall, this study provides new cyclotide sequences not yet found naturally that could be valuable starting points for the understanding of cyclotide evolution and for further optimization as drug leads.

Structure, mechanism, and evolution of the last step in vitamin C biosynthesis

Photosynthetic organisms, fungi, and animals comprise distinct pathways for vitamin C biosynthesis. Besides this diversity, the final biosynthetic step consistently involves an oxidation reaction carried out by the aldonolactone oxidoreductases. Here, we study the origin and evolution of the diversified activities and substrate preferences featured by these flavoenzymes using molecular phylogeny, kinetics, mutagenesis, and crystallographic experiments. We find clear evidence that they share a common ancestor. A flavin-interacting amino acid modulates the reactivity with the electron acceptors, including oxygen, and determines whether an enzyme functions as an oxidase or a dehydrogenase. We show that a few side chains in the catalytic cavity impart the reaction stereoselectivity. Ancestral sequence reconstruction outlines how these critical positions were affixed to specific amino acids along the evolution of the major eukaryotic clades. During Eukarya evolution, the aldonolactone oxidoreductases adapted to the varying metabolic demands while retaining their overarching vitamin C-generating function.