Ancestral Sequence Research Articles

Humanization and functional characterization of enhanced coagulation factor IX variants identified through ancestral sequence reconstruction

BackgroundLaboratory resurrection of ancient coagulation factor (F) IX variants generated through ancestral sequence reconstruction led to the discovery of a FIX variant, designated An96, which possesses enhanced specific activity independent of and additive to that provided by human p.Arg384Lys, referred to as FIX-Padua. ObjectivesThe goal of the current study was to identify the amino acid substitution(s) responsible for the enhanced activity of An96 and create a humanized An96 FIX transgene for gene therapy application. MethodsReductionist screening approaches, including domain swapping and scanning residue substitution, were used and guided by one-stage FIX activity assays. In vitro characterization of top candidates included recombinant high-purity preparation, specific activity determination, and enzyme kinetic analysis. Final candidates were packaged into adeno-associated viral (AAV) vectors and delivered to hemophilia B mice. ResultsFive of 42 total amino acid substitutions in An96 appear sufficient to retain the enhanced activity of An96 in an otherwise human FIX variant. Additional substitution of the Padua variant further increased the specific activity 5-fold. This candidate, designated ET9, demonstrated 51-fold greater specific activity than hFIX. AAV2/8-ET9 treated hemophilia B mice produced plasma FIX activities equivalent to those observed previously for AAV2/8-An96-Padua, which were 10-fold higher than AAV2/8-hFIX-Padua. ConclusionStarting from computationally inferred ancient FIX sequences, novel amino acid substitutions conferring activity enhancement were identified and translated into an AAV-FIX gene therapy cassette demonstrating high potency. This ancestral sequence reconstruction discovery and sequence mapping refinement approach represents a promising platform for broader protein drug and gene therapy candidate optimization.

FireProt 2.0: web-based platform for the fully automated design of thermostable proteins.

Thermostable proteins find their use in numerous biomedical and biotechnological applications. However, the computational design of stable proteins often results in single-point mutations with a limited effect on protein stability. However, the construction of stable multiple-point mutants can prove difficult due to the possibility of antagonistic effects between individual mutations. FireProt protocol enables the automated computational design of highly stable multiple-point mutants. FireProt 2.0 builds on top of the previously published FireProt web, retaining the original functionality and expanding it with several new stabilization strategies. FireProt 2.0 integrates the AlphaFold database and the homology modeling for structure prediction, enabling calculations starting from a sequence. Multiple-point designs are constructed using the Bron-Kerbosch algorithm minimizing the antagonistic effect between the individual mutations. Users can newly limit the FireProt calculation to a set of user-defined mutations, run a saturation mutagenesis of the whole protein or select rigidifying mutations based on B-factors. Evolution-based back-to-consensus strategy is complemented by ancestral sequence reconstruction. FireProt 2.0 is significantly faster and a reworked graphical user interface broadens the tool's availability even to users with older hardware. FireProt 2.0 is freely available at http://loschmidt.chemi.muni.cz/fireprotweb.

Development of a Novel L-Asparaginase with Enhanced Properties Using Ancestral Sequence Reconstruction

The overall survival of children with acute lymphoblastic leukemia (ALL) now exceeds 90% due to the precise disease risk stratification and intensification of multi-agent chemotherapy regimens. L-asparaginase (L-ASNase) is an indispensable biotherapeutic enzyme used in treatment of pediatric ALL. Discontinuation of L-ASNase treatment due to immunological or adverse reactions leads to significantly inferior disease-free survival (73% ±7% vs 92% ± 2%; P <0.01, Gupta et al. JCO 2020). Clinical L-ASNases are bacterial in origin, derived from either Escherichia coli or Erwinia chrysanthemi and are therefore highly immunogenic. Additionally, emerging data suggests that the concurrent glutaminase activity seen in bacterial asparaginases exacerbates its liver and pancreatic toxicity, especially in adults thereby limiting its extensive use beyond children. Therefore, finding alternative sources with preferable biochemical features remains paramount. The inferior catalytic properties of human L-ASNase restricts its functionality in ALL treatment. However, characterization has shown that guinea pig (GP) L-ASNase demonstrates significant anti-tumor properties and improved enzyme kinetics despite a 69.8% amino acid sequence identity with human L-ASNase, as compared to 30% for bacterial L-ASNases. Additionally, guinea pig L-ASNase has no glutaminase activity. Thus, to overcome the major deficiencies of bacterial L-ASNase and develop a more humanized form of L-ASNase we performed ancestral sequence reconstruction (ASR). ASR is a protein drug discovery / optimization platform that predicts ancient DNA and protein sequences from extant sequences. This approach permits higher-resolution mapping potentially highlighting functional resides. The objective therefore is to exploit functional diversity refined by natural selection using ASR to identify therapeutic variants with superior properties. Recent work on guinea pig (GP) L-ASNase (Schalk et al. J Biol Chem 2014) identified a low micromolar (μM) Michaelis constant (Km) at 57.7 ± 6.4 μM essential for hydrolyzing the extra-cellular pool of asparagine ~50 μM in human blood; thus, equipping us with the ideal mammalian L-ASNase ortholog to serve as a template to resurrect ancient ASNase sequences. Ancestral L-ASNase sequences were generated utilizing 54 extant L- asparaginase sequences, aligned using MUSCLE and an evolutionary tree inferred using MrBayes. A total of 53 ancestral nodes were identified, and the ten ancestral L-ASNase variants spanning the ancient primate and guinea pig (GP) lineage were resurrected. The sequence identity of these ancestral L-ASNase variants ranged from 81 to 98% when compared to human L-ASNase. E. coli codon optimized complementary DNA (cDNA) sequences were synthesized and subcloned intoan expression vector. Correctly cloned plasmids were transformed into E. coli BL21 (DE3) cells for protein expression. An-ASNase candidates were isolated by Ni 2+ affinity chromatography and then further purified by size exclusion chromatography. Asparaginase activity of the An-ASNase candidates were determined via a modified Nessler's reagent assay by using a continuous spectroscopic enzyme-coupled assay. Three of the ancestral L-ASNase candidates, An-88, An-104 and An-107 demonstrated excellent asparaginase activity when tested at an enzyme concentration of 0.1 mg/mL and an excess substrate concentration, with activity within a comparable range to the clinically relevant E. coli and Erwinia asparaginases. An-88 had 81% similarity whereas both An-104 and An-107 ASNases had 88% identity to human L-ASNase, which is in stark contrast to the only 30% similarity seen with bacterial L-ASNases. Detailed enzymatic characterization of these three lead drug candidates is currently being performed. Preliminary cytotoxicity data testing the ancestral L-ASNase candidates against a T-cell ALL cell line CCRF-CEM showed an expected dose response, with An-107 showing the highest cytotoxicity. Our findings demonstrate that ASR is an efficient protein drug design/optimization platform that we have now extended beyond our previous work with recombinant coagulation factors (Zakas et al. Nat Biotech 2017, Knight et al. Blood Adv 2021) to therapeutic enzymes; thereby, enabling the creation of novel recombinant L-ASNases with translational implications due to improved enzymatic specificity and reduced toxicity.

Representing and extending ensembles of parsimonious evolutionary histories with a directed acyclic graph

In many situations, it would be useful to know not just the best phylogenetic tree for a given data set, but the collection of high-quality trees. This goal is typically addressed using Bayesian techniques, however, current Bayesian methods do not scale to large data sets. Furthermore, for large data sets with relatively low signal one cannot even store every good tree individually, especially when the trees are required to be bifurcating. In this paper, we develop a novel object called the “history subpartition directed acyclic graph” (or “history sDAG” for short) that compactly represents an ensemble of trees with labels (e.g. ancestral sequences) mapped onto the internal nodes. The history sDAG can be built efficiently and can also be efficiently trimmed to only represent maximally parsimonious trees. We show that the history sDAG allows us to find many additional equally parsimonious trees, extending combinatorially beyond the ensemble used to construct it. We argue that this object could be useful as the “skeleton” of a more complete uncertainty quantification.

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

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

Dissection of the role of a Src homology 3 domain in the evolution of binding preference of paralogous proteins.

Protein-protein interactions (PPIs) drive many cellular processes. Some interactions are directed by Src homology 3 (SH3) domains that bind proline-rich motifs on other proteins. The evolution of the binding specificity of SH3 domains is not completely understood, particularly following gene duplication. Paralogous genes accumulate mutations that can modify protein functions and, for SH3 domains, their binding preferences. Here, we examined how the binding of the SH3 domains of 2 paralogous yeast type I myosins, Myo3 and Myo5, evolved following duplication. We found that the paralogs have subtly different SH3-dependent interaction profiles. However, by swapping SH3 domains between the paralogs and characterizing the SH3 domains freed from their protein context, we find that very few of the differences in interactions, if any, depend on the SH3 domains themselves. We used ancestral sequence reconstruction to resurrect the preduplication SH3 domains and examined, moving back in time, how the binding preference changed. Although the most recent ancestor of the 2 domains had a very similar binding preference as the extant ones, older ancestral domains displayed a gradual loss of interaction with the modern interaction partners when inserted in the extant paralogs. Molecular docking and experimental characterization of the free ancestral domains showed that their affinity with the proline motifs is likely not the cause for this loss of binding. Taken together, our results suggest that a SH3 and its host protein could create intramolecular or allosteric interactions essential for the SH3-dependent PPIs, making domains not functionally equivalent even when they have the same binding specificity.

Structural and functional analysis of hyper-thermostable ancestral L-amino acid oxidase that can convert Trp derivatives to D-forms by chemoenzymatic reaction

Production of D-amino acids (D-AAs) on a large-scale enables to provide precursors of peptide therapeutics. In this study, we designed a novel L-amino acid oxidase, HTAncLAAO2, by ancestral sequence reconstruction, exhibiting high thermostability and long-term stability. The crystal structure of HTAncLAAO2 was determined at 2.2 Å by X-ray crystallography, revealing that the enzyme has an octameric form like a “ninja-star” feature. Enzymatic property analysis demonstrated that HTAncLAAO2 exhibits three-order larger kcat/Km values towards four L-AAs (L-Phe, L-Leu, L-Met, and L-Ile) than that of L-Trp. Through screening the variants, we obtained the HTAncLAAO2(W220A) variant, which shows a > 6-fold increase in kcat value toward L-Trp compared to the original enzyme. This variant applies to synthesizing enantio-pure D-Trp derivatives from L- or rac-forms at a preparative scale. Given its excellent properties, HTAncLAAO2 would be a starting point for designing novel oxidases with high activity toward various amines and AAs.

Deciphering Enzyme Mechanisms with Engineered Ancestors and Substrate Analogues

AbstractEnvironmentally friendly industrial and biotech processes greatly benefit from enzyme‐based technologies. Their use is often possible only when the enzyme‐catalytic mechanism is thoroughly known. Thus, atomic‐level knowledge of a Michaelis enzyme‐substrate complex, revealing molecular details of substrate recognition and catalytic chemistry, is crucial for understanding and then rationally extending or improving enzyme‐catalyzed reactions. However, many known enzymes sample huge protein conformational space, often preventing complete structural characterization by X‐ray crystallography. Moreover, using a cognate substrate is problematic since its conversion into a reaction product in the presence of the enzyme will prevent the capture of the enzyme‐substrate conformation in an activated state. Here, we outlined how to deal with such obstacles, focusing on the recent discovery of a Renilla‐type bioluminescence reaction mechanism facilitated by a combination of engineered ancestral enzyme and the availability of a non‐oxidizable luciferin analogue. The automated ancestral sequence reconstructions using FireProtASR provided a thermostable enzyme suited for structural studies, and a stable luciferin analogue azacoelenterazine provided a structurally cognate chemical incapable of catalyzed oxidation. We suggest that an analogous strategy can be applied to various enzymes with unknown catalytic mechanisms and poor crystallizability.

Heme Pocket Structure and Its Functional Implications in an Ancestral Globin Protein.

Proteins have undergone evolutionary processes to achieve optimal stability, increased functionality, and novel functions. Comparative analysis of existent and ancestral proteins provides insights into the factors that influence protein stability and function. Ancestral sequence reconstruction allows us to deduce the amino acid sequences of ancestral proteins. Here, we present the structural and functional characteristics of an ancestral protein, AncMH, reconstructed to be the last common ancestor of hemoglobins and myoglobins. Our findings reveal that AncMH harbors heme and that the heme binds oxygen. Furthermore, we demonstrate that the ferrous heme in AncMH is pentacoordinated, similar to that of human adult hemoglobin and horse myoglobin. A detailed comparison of the heme pocket structure indicates that the heme pocket in AncMH is more similar to that of hemoglobin than that of myoglobin. However, the autoxidation of AncMH is faster than that of both hemoglobin and myoglobin. Collectively, our results suggest that ancestral proteins of hemoglobins and myoglobins evolved in steps, including the hexa- to pentacoordination transition, followed by stabilization of the oxygen-bound form.

Functional divergence of CYP76AKs shapes the chemodiversity of abietane-type diterpenoids in genus Salvia

The genus Salvia L. (Lamiaceae) comprises myriad distinct medicinal herbs, with terpenoids as one of their major active chemical groups. Abietane-type diterpenoids (ATDs), such as tanshinones and carnosic acids, are specific to Salvia and exhibit taxonomic chemical diversity among lineages. To elucidate how ATD chemical diversity evolved, we carried out large-scale metabolic and phylogenetic analyses of 71 Salvia species, combined with enzyme function, ancestral sequence and chemical trait reconstruction, and comparative genomics experiments. This integrated approach showed that the lineage-wide ATD diversities in Salvia were induced by differences in the oxidation of the terpenoid skeleton at C-20, which was caused by the functional divergence of the cytochrome P450 subfamily CYP76AK. These findings present a unique pattern of chemical diversity in plants that was shaped by the loss of enzyme activity and associated catalytic pathways.

MiR-430 microRNA Family in Fishes: Molecular Characterization and Evolution.

The miR-430 microRNA family has been described in multiple fish species as one of the first microRNAs expressed by the zygote. It has been suggested that this family is implicated in maternal mRNA elimination, but may also play a role in steroidogenesis, sexual differentiation, and flatfish metamorphosis. The miR-430 sequences have been found in multiple-copy tandem clusters but evidence of their conservation outside of teleost fishes is scarce. In the present study, we have characterized the tandem repeats organization of these microRNAs in different fish species, both model and of interest in aquaculture. A phylogenetic analysis of this family has allowed us to identify that the miR-430 duplication, which took place before the Chondrostei and Neopterygii groups' divergence, has resulted in three variants ("a", "b", and "c"). According to our data, variant "b" is the most closely related to the ancestral sequence. Furthermore, we have detected isolated instances of the miR-430 repeat subunit in some species, which suggests that this microRNA family may be affected by DNA rearrangements. This study provides new data about the abundance, variability, and organization of the miR-430 family in fishes.

SARS-CoV-2 nonstructural protein 6 from Alpha to Omicron: evolution of a transmembrane protein.

We recently reported that mutations in both the spike glycoprotein and nonstructural protein 6 (nsp6) were associated with attenuation of the SARS-CoV-2 Omicron BA.1 variant. While mutations in spike allow evasion of neutralizing antibodies and promote specific modes of viral entry, the role of nsp6 mutations in pathogenesis is less clear. Nsp6 is essential for modifying the endoplasmic reticulum and generating double-membrane vesicles, the site of viral RNA replication. To investigate the evolution of nsp6, we evaluated 91,596 high-confidence human SARS-CoV-2 whole-genome sequences across 19 variants and lineages. While nsp6 of early variants of concern, such as Alpha, Beta, and Gamma, carried a triple amino acid deletion (106-108, termed ΔSGF), the Delta, Epsilon, and Mu lineages retained the ancestral nsp6 sequence. For nsp6 in the emerging Omicron variants, we report a transition from an amino acid 105-107 ΔLSG deletion in BA.1 to increased dominance of the ΔSGF in BA.2 and subsequent lineages. Our findings indicate that deletion within nsp6 was independently selected in multiple lineages of SARS-CoV-2, both early and late in the pandemic. Analysis of SARS-CoV-2-related coronaviruses in bats and pangolins revealed nsp6 sequences similar to the ancestral SARS-CoV-2 virus, indicating that the deletion in nsp6 may be an adaptation to replication in humans. Analysis of nsp6 sequences from multiple coronaviruses predicts a multipass transmembrane protein with a conserved C-terminal domain. Monitoring and evaluating changes in nsp6 and other nonstructural proteins will contribute to our understanding of factors associated with the attenuation of pandemic coronaviruses. IMPORTANCE There is an ongoing need to evaluate genetic changes in SARS-CoV-2 for effects on transmission and pathogenesis. We recently reported an unexpected role for replicase component nsp6, in addition to changes in spike, in the attenuation of Omicron BA.1. In this commentary, we document a triple-amino-acid deletion in a predicted lumenal domain of nsp6 that was found in multiple, independent variants of SARS-CoV-2, including all recent Omicron lineages. Furthermore, we modeled the predicted structure of nsp6, implicating a multipass transmembrane architecture as conserved in members of the Coronaviridae family. This information can guide future studies investigating the role of nsp6 in the pathogenesis of existing and emerging coronaviruses.

Phylogenetic Reconstruction and Functional Characterization of the Ancestral Nef Protein of Primate Lentiviruses.

Nef is an accessory protein unique to the primate HIV-1, HIV-2, and SIV lentiviruses. During infection, Nef functions by interacting with multiple host proteins within infected cells to evade the immune response and enhance virion infectivity. Notably, Nef can counter immune regulators such as CD4 and MHC-I, as well as the SERINC5 restriction factor in infected cells. In this study, we generated a posterior sample of time-scaled phylogenies relating SIV and HIV Nef sequences, followed by reconstruction of ancestral sequences at the root and internal nodes of the sampled trees up to the HIV-1 Group M ancestor. Upon expression of the ancestral primate lentivirus Nef protein within CD4+ HeLa cells, flow cytometry analysis revealed that the primate lentivirus Nef ancestor robustly downregulated cell-surface SERINC5, yet only partially downregulated CD4 from the cell surface. Further analysis revealed that the Nef-mediated CD4 downregulation ability evolved gradually, while Nef-mediated SERINC5 downregulation was recovered abruptly in the HIV-1/M ancestor. Overall, this study provides a framework to reconstruct ancestral viral proteins and enable the functional characterization of these proteins to delineate how functions could have changed throughout evolutionary history.

An Inferred Ancestral CotA Laccase with Improved Expression and Kinetic Efficiency.

Laccases are widely used in industrial production due to their broad substrate availability and environmentally friendly nature. However, the pursuit of laccases with superior stability and increased heterogeneous expression to meet industry demands appears to be an ongoing challenge. To address this challenge, we resurrected five ancestral sequences of laccase BsCotA and their homologues. All five variants were successfully expressed in soluble and functional forms with improved expression levels in Escherichia coli. Among the five variants, three exhibited higher catalytic rates, thermal stabilities, and acidic stabilities. Notably, AncCotA2, the best-performing variant, displayed a kcat/KM of 7.5 × 105 M-1·s-1, 5.2-fold higher than that of the wild-type BsCotA, an improved thermo- and acidic stability, and better dye decolorization ability. This study provides a laccase variant with high application potential and presents a new starting point for future enzyme engineering.

Evolutionary origin of B family DNA-dependent DNA polymerases from retrotranscriptases

The evolution of DNA and DNA polymerases was a crucial step in the evolution of life on Earth. In the present work, we reconstruct the ancestral sequence and structure for the B family polymerases. Using comparative analyses, we infer the transient state between the ancestor retrotranscriptase and the contemporary B family DNA polymerases. Exonuclease motif was detected in the primary ancestral sequence, as well as an elongation-functioning motif. It is remarkable that the ancestral molecule is more comparable to the retrotranscriptases in terms of structural domains, even though we discovered similarities in the primary sequence with proteins from the B family of DNA polymerases. The present B family proteins differ structurally from retrotranscriptases the most, although the reconstruction of the ancestor protein was able to capture the transitional steps between these two families of polymerases.

Extensive Diversity in Escherichia coli Group 3 Capsules Is Driven by Recombination and Plasmid Transfer from Multiple Species.

Bacterial capsules provide protection against environmental challenges and host immunity. Historically, Escherichia coli K serotyping scheme, which relies on the hypervariable capsules, has identified around 80 K forms that fall into four distinct groups. Based on recent work by us and others, we predicted that E. coli capsular diversity is grossly underestimated. We exploited group 3 capsule gene clusters, the best genetically defined capsule group in E. coli, to analyze publicly available E. coli sequences for overlooked capsular diversity within the species. We report the discovery of seven novel group 3 clusters that fall into two distinct subgroups (3A and 3B). The majority of the 3B capsule clusters were found on plasmids, contrary to the defining feature of group 3 capsule genes localizing at the serA locus on the E. coli chromosome. Other new group 3 capsule clusters were derived from ancestral sequences through recombination events between shared genes found within the serotype variable central region 2. Intriguingly, flanking regions 1 and 3, known to be conserved areas among capsule clusters, showed considerable intra-subgroup variation in clusters from the 3B subgroup, containing genes of shared ancestry with other Enterobacteriaceae species. Variation of group 3 kps clusters within dominant E. coli lineages, including multidrug-resistant pathogenic lineages, further supports that E. coli capsules are undergoing rigorous change. Given the pivotal role of capsular polysaccharides in phage predation, our findings raise attention to the need of monitoring kps evolutionary dynamics in pathogenic E. coli in supporting phage therapy. IMPORTANCE Capsular polysaccharides protect pathogenic bacteria against environmental challenges, host immunity, and phage predations. The historical Escherichia coli K typing scheme, which relies on the hypervariable capsular polysaccharide, has identified around 80 different K forms that fall into four distinct groups. Taking advantage of the supposedly compact and genetically well-defined group 3 gene clusters, we analyzed published E. coli sequences to identify seven new gene clusters and revealed an unexpected capsular diversity. Genetic analysis revealed that group 3 gene clusters shared closely related serotype-specific region 2 and were diversified through recombination events and plasmid transfer between multiple Enterobacteriaceae species. Overall, capsular polysaccharides in E. coli are undergoing rigorous change. Given the pivotal role capsules play in phage interactions, this work highlighted the need to monitor the evolutionary dynamics of capsules in pathogenic E. coli for effective phage therapy.

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.

Inference of ancestral alleles in the pig reference genome.

The authors declare no conflict of interest. The final ancestral sequences in FASTA (https://doi.org/10.57760/sciencedb.07853), the multi-genome alignment in HAL (https://doi.org/10.57760/sciencedb.07855) and the ancestral states of all Sscrofa11.1 positions validated by re-sequencing data (https://doi.org/10.57760/sciencedb.07858) are accessible at ScienceDB. Appendix S1 Appendix S2 Please note: The publisher is not responsible for the content or functionality of any supporting information supplied by the authors. Any queries (other than missing content) should be directed to the corresponding author for the article.

Identification of a covert evolutionary pathway between two protein folds

Although homologous protein sequences are expected to adopt similar structures, some amino acid substitutions can interconvert α-helices and β-sheets. Such fold switching may have occurred over evolutionary history, but supporting evidence has been limited by the: (1) abundance and diversity of sequenced genes, (2) quantity of experimentally determined protein structures, and (3) assumptions underlying the statistical methods used to infer homology. Here, we overcome these barriers by applying multiple statistical methods to a family of ~600,000 bacterial response regulator proteins. We find that their homologous DNA-binding subunits assume divergent structures: helix-turn-helix versus α-helix + β-sheet (winged helix). Phylogenetic analyses, ancestral sequence reconstruction, and AlphaFold2 models indicate that amino acid substitutions facilitated a switch from helix-turn-helix into winged helix. This structural transformation likely expanded DNA-binding specificity. Our approach uncovers an evolutionary pathway between two protein folds and provides a methodology to identify secondary structure switching in other protein families.

Virus Pop-Expanding Viral Databases by Protein Sequence Simulation.

The improvement of our knowledge of the virosphere, which includes unknown viruses, is a key area in virology. Metagenomics tools, which perform taxonomic assignation from high throughput sequencing datasets, are generally evaluated with datasets derived from biological samples or in silico spiked samples containing known viral sequences present in public databases, resulting in the inability to evaluate the capacity of these tools to detect novel or distant viruses. Simulating realistic evolutionary directions is therefore key to benchmark and improve these tools. Additionally, expanding current databases with realistic simulated sequences can improve the capacity of alignment-based searching strategies for finding distant viruses, which could lead to a better characterization of the "dark matter" of metagenomics data. Here, we present Virus Pop, a novel pipeline for simulating realistic protein sequences and adding new branches to a protein phylogenetic tree. The tool generates simulated sequences with substitution rate variations that are dependent on protein domains and inferred from the input dataset, allowing for a realistic representation of protein evolution. The pipeline also infers ancestral sequences corresponding to multiple internal nodes of the input data phylogenetic tree, enabling new sequences to be inserted at various points of interest in the group studied. We demonstrated that Virus Pop produces simulated sequences that closely match the structural and functional characteristics of real protein sequences, taking as an example the spike protein of sarbecoviruses. Virus Pop also succeeded at creating sequences that resemble real sequences not included in the databases, which facilitated the identification of a novel pathogenic human circovirus not included in the input database. In conclusion, Virus Pop is helpful for challenging taxonomic assignation tools and could help improve databases to better detect distant viruses.