Palaeotemperature trend for Precambrian life inferred from resurrected proteins

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.

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.

Primordial magnetotaxis in putative giant paleoproterozoic magnetofossils

Magnetotactic bacteria produce chains of nanoscopic iron minerals used for navigation, which can be preserved over geological timescales in the form of magnetofossils. Micrometer-sized magnetite crystals with unusual shapes suggesting a biologically controlled mineralization have been found in the geological record and termed giant magnetofossils. The biological origin and function of giant magnetofossils remains unclear, due to the lack of modern analogues to giant magnetofossils. Using distinctive Ptychographic nanotomography data of Precambrian (1.88 Ga) rocks, we recovered the morphology of micrometric cuboid grains of iron oxides embedded in an organic filamentous fossil to construct synthetic magnetosomes. Their morphology is different from that of previously found giant magnetofossils, but their occurrence in filamentous microfossils and micromagnetic simulations support the hypothesis that they could have functioned as a navigation aid, akin to modern magnetosomes.

Evolutionary analysis reveals the origin of sodium coupling in glutamate transporters

Secondary active membrane transporters harness the energy of ion gradients to concentrate their substrates. Homologous transporters evolved to couple transport to different ions in response to changing environments and needs. The bases of such diversification, and thus principles of ion coupling, are unexplored. Employing phylogenetics and ancestral protein reconstruction, we investigated sodium-coupled transport in prokaryotic glutamate transporters, a mechanism ubiquitous across life domains and critical to neurotransmitter recycling in humans. We found that the evolutionary transition from sodium-dependent to independent substrate binding to the transporter preceded changes in the coupling mechanism. Structural and functional experiments suggest that the transition entailed allosteric mutations, making sodium binding dispensable without affecting ion-binding sites. Allosteric tuning of transporters' energy landscapes might be a widespread route of their functional diversification.

Hot springs viruses at Yellowstone National Park have ancient origins and are adapted to thermophilic hosts

Geothermal springs house unicellular red algae in the class Cyanidiophyceae that dominate the microbial biomass at these sites. Little is known about host-virus interactions in these environments. We analyzed the virus community associated with red algal mats in three neighboring habitats (creek, endolithic, soil) at Lemonade Creek, Yellowstone National Park (YNP), USA. We find that despite proximity, each habitat houses a unique collection of viruses, with the giant viruses, Megaviricetes, dominant in all three. The early branching phylogenetic position of genes encoded on metagenome assembled virus genomes (vMAGs) suggests that the YNP lineages are of ancient origin and not due to multiple invasions from mesophilic habitats. The existence of genomic footprints of adaptation to thermophily in the vMAGs is consistent with this idea. The Cyanidiophyceae at geothermal sites originated ca. 1.5 Bya and are therefore relevant to understanding biotic interactions on the early Earth.

Extant Sequence Reconstruction: The Accuracy of Ancestral Sequence Reconstructions Evaluated by Extant Sequence Cross-Validation

Ancestral sequence reconstruction (ASR) is a phylogenetic method widely used to analyze the properties of ancient biomolecules and to elucidate mechanisms of molecular evolution. Despite its increasingly widespread application, the accuracy of ASR is currently unknown, as it is generally impossible to compare resurrected proteins to the true ancestors. Which evolutionary models are best for ASR? How accurate are the resulting inferences? Here we answer these questions using a cross-validation method to reconstruct each extant sequence in an alignment with ASR methodology, a method we term "extant sequence reconstruction" (ESR). We thus can evaluate the accuracy of ASR methodology by comparing ESR reconstructions to the corresponding known true sequences. We find that a common measure of the quality of a reconstructed sequence, the average probability, is indeed a good estimate of the fraction of correct amino acids when the evolutionary model is accurate or overparameterized. However, the average probability is a poor measure for comparing reconstructions from different models, because, surprisingly, a more accurate phylogenetic model often results in reconstructions with lower probability. While better (more predictive) models may produce reconstructions with lower sequence identity to the true sequences, better models nevertheless produce reconstructions that are more biophysically similar to true ancestors. In addition, we find that a large fraction of sequences sampled from the reconstruction distribution may have fewer errors than the single most probable (SMP) sequence reconstruction, despite the fact that the SMP has the lowest expected error of all possible sequences. Our results emphasize the importance of model selection for ASR and the usefulness of sampling sequence reconstructions for analyzing ancestral protein properties. ESR is a powerful method for validating the evolutionary models used for ASR and can be applied in practice to any phylogenetic analysis of real biological sequences. Most significantly, ESR uses ASR methodology to provide a general method by which the biophysical properties of resurrected proteins can be compared to the properties of the true protein.

Oxygen isotope ensemble reveals Earth's seawater, temperature, and carbon cycle history

Earth's persistent habitability since the Archean remains poorly understood. Using an oxygen isotope ensemble approach-comprising shale, iron oxide, carbonate, silica, and phosphate records-we reconcile a multibillion-year history of seawater δ18O, temperature, and marine and terrestrial clay abundance. Our results reveal a rise in seawater δ18O and a temperate Proterozoic climate distinct to interpretations of a hot early Earth, indicating a strongly buffered climate system. Precambrian sediments are enriched in marine authigenic clay, with prominent reductions occurring in concert with Paleozoic and Cenozoic cooling, the expansion of siliceous life, and the radiation of land plants. These findings support the notion that shifts in the locus and extent of clay formation contributed to seawater 18O enrichment, clement early Earth conditions, major climate transitions, and climate stability through the reverse weathering feedback.

Biotechnological innovations in nanocellulose production from waste biomass with a focus on pineapple waste

New materials' synthesis and utilization have shown many critical challenges in healthcare and other industrial sectors as most of these materials are directly or indirectly developed from fossil fuel resources. Environmental regulations and sustainability concepts have promoted the use of natural compounds with unique structures and properties that can be biodegradable, biocompatible, and eco-friendly. In this context, nanocellulose (NC) utility in different sectors and industries is reported due to their unique properties including biocompatibility and antimicrobial characteristics. The bacterial nanocellulose (BNC)-based materials have been synthesized by bacterial cells and extracted from plant waste materials including pineapple plant waste biomass. These materials have been utilized in the form of nanofibers and nanocrystals. These materials are found to have excellent surface properties, low density, and good transparency, and are rich in hydroxyl groups for their modifications to other useful products. These materials are well utilized in different sectors including biomedical or health care centres, nanocomposite materials, supercapacitors, and polymer matrix production. This review explores different approaches for NC production from pineapple waste residues using biotechnological interventions, approaches for their modification, and wider applications in different sectors. Recent technological developments in NC production by enzymatic treatment are critically discussed. The utilization of pineapple waste-derived NC from a bioeconomic perspective is summarized in the paper. The chemical composition and properties of nanocellulose extracted from pineapple waste may have unique characteristics compared to other sources. Pineapple waste for nanocellulose production aligns with the principles of sustainability, waste reduction, and innovation, making it a promising and novel approach in the field of nanocellulose materials.

Evolution of optimal growth temperature in Asgard archaea inferred from the temperature dependence of GDP binding to EF-1A

The archaeal ancestor of eukaryotes apparently belonged to the phylum Asgardarchaeota, but the ecology and evolution of Asgard archaea are poorly understood. The optimal GDP-binding temperature of a translation elongation factor (EF-1A or EF-Tu) has been previously shown to correlate with the optimal growth temperature of diverse prokaryotes. Here, we reconstruct ancestral EF-1A sequences and experimentally measure the optimal GDP-binding temperature of EF-1A from ancient and extant Asgard archaea, to infer the evolution of optimal growth temperatures in Asgardarchaeota. Our results suggest that the Asgard ancestor of eukaryotes was a moderate thermophile, with an optimal growth temperature around 53 °C. The origin of eukaryotes appears to coincide with a transition from thermophilic to mesophilic lifestyle during the evolution of Asgard archaea.

Comparative analysis of reconstructed ancestral proteins with their extant counterparts suggests primitive life had an alkaline habitat

To understand the origin and early evolution of life it is crucial to establish characteristics of the primordial environment that facilitated the emergence and evolution of life. One important environmental factor is the pH of the primordial environment. Here, we assessed the pH-dependent thermal stabilities of previously reconstructed ancestral nucleoside diphosphate kinases and ribosomal protein uS8s. The selected proteins were likely to be present in ancient organisms such as the last common ancestor of bacteria and that of archaea. We also assessed the thermal stability of homologous proteins from extant acidophilic, neutralophilic, and alkaliphilic microorganisms as a function of pH. Our results indicate that the reconstructed ancestral proteins are more akin to those of extant alkaliphilic bacteria, which display greater stability under alkaline conditions. These findings suggest that the common ancestors of bacterial and archaeal species thrived in an alkaline environment. Moreover, we demonstrate the reconstruction method employed in this study is a valuable technique for generating alkali-tolerant proteins that can be used in a variety of biotechnological and environmental applications.

Characterization of a Novel Thermostable 7α-Hydroxysteroid Dehydrogenase

BACKGROUND

7α-Hydroxysteroid dehydrogenase (7α-HSDH) plays a pivotal role in vivo in the biotransformation of secondary bile acids and has great potential in industrial biosynthesis due to its broad substrate specificity. In this study, we expressed and characterized a novel thermostable 7α-HSDH (named Sa 7α-HSDH).

METHODS

The DNA sequence was derived from the black bear gut microbiome metagenomic sequencing data, and the coding sequence of Sa 7α-HSDH was chemically synthesized. The heterologous expression of the enzyme was carried out using the pGEX-6p-1 vector. Subsequently, the activity of the purified enzyme was studied by measuring the absorbance change at 340 nm. Finally, the three-dimensional structure was predicted with AlphaFold2.

RESULTS

Coenzyme screening results confirmed it to be NAD(H) dependent. Substrate specificity test revealed that Sa 7α-HSDH could catalyze taurochenodeoxycholic acid (TCDCA) with catalytic efficiency (kcat/Km) 3.81 S-1 mM-1. The optimum temperature of Sa 7α-HSDH was measured to be 75°C, confirming that it belongs to thermophilic enzymes. Additionally, its thermostability was assessed using an accelerated stability test over 32 hours. The catalytic activity of Sa 7α-HSDH remained largely unchanged for the first 24 hours and retained over 90% of its functionality after 32 hours at 50°C. Sa 7α-HSDH exhibited maximal activity at pH 10. The effect of metal ions-K+, Na+, Mg2+ and Cu2+-on the enzymatic activity of Sa 7α-HSDH was investigated. Only Mg2+ was observed to enhance the enzyme's activity by 27% at a concentration of 300 mM. Neither K+ nor Na+ had a significant influence on activity. Only Cu2+ was found to reduce enzyme activity.

CONCLUSION

We characterized the thermostable 7α-HSDH, which provides a promising biocatalyst for bioconversion of steroids at high reaction temperatures.

Influence of the Weak Nuclear Force on Metal-Promoted Autocatalytic Strecker Synthesis of Amino Acids: Formation of a Chiral Pool of Precursors for Prebiotic Peptide and Protein Synthesis

Natural chiral amino acids typically adopt an L structural configuration. While a preference for specific molecular chiralities is observed throughout biology and cellular chemistry, the origins of this preference are unclear. In a previous report the origin of enantiomeric selectivity was analyzed in terms of an "RNA World" model, and a pathway to a chiral preference for d-ribose was proposed based on the autocatalytic transformation of glyceraldehyde as a precursor to the formation of sugars. Metal-ion-promoted catalysis allows the parity non-conserving (PNC) weak nuclear interaction to influence the chirality of a nascent chiral carbon center. Since the PNC effect is the only natural property with an inherent handedness, it is an obvious candidate to influence enantiomeric preference from a catalytic reaction performed over geologically relevant time scales. The PNC influence requires and emphasizes the important role of catalytic metal ions in primordial chemistry. In this study, the impact of geologically available divalent calcium and higher Z alkaline earth elements are examined as mediators of chiral preference. Detailed calculations of the magnitude of the effect are presented, including the influence of time, temperature, pH, and metal ion identity. It is concluded that metal ions can direct chiral preference for amino acid synthesis via a metal-promoted autocatalytic Strecker reaction within a relatively short geological timeframe, thereby providing a pool of l-amino acids for catalytic chemistry evolving either from an RNA-world model of molecular evolution or alternative pathways to protein synthesis.

Evolution towards simplicity in bacterial small heat shock protein system

Evolution can tinker with multi-protein machines and replace them with simpler single-protein systems performing equivalent functions in an equally efficient manner. It is unclear how, on a molecular level, such simplification can arise. With ancestral reconstruction and biochemical analysis, we have traced the evolution of bacterial small heat shock proteins (sHsp), which help to refold proteins from aggregates using either two proteins with different functions (IbpA and IbpB) or a secondarily single sHsp that performs both functions in an equally efficient way. Secondarily single sHsp evolved from IbpA, an ancestor specialized in strong substrate binding. Evolution of an intermolecular binding site drove the alteration of substrate binding properties, as well as the formation of higher-order oligomers. Upon two mutations in the α-crystallin domain, secondarily single sHsp interacts with aggregated substrates less tightly. Paradoxically, less efficient binding positively influences the ability of sHsp to stimulate substrate refolding, since the dissociation of sHps from aggregates is required to initiate Hsp70-Hsp100-dependent substrate refolding. After the loss of a partner, IbpA took over its role in facilitating the sHsp dissociation from an aggregate by weakening the interaction with the substrate, which became beneficial for the refolding process. We show that the same two amino acids introduced in modern-day systems define whether the IbpA acts as a single sHsp or obligatorily cooperates with an IbpB partner. Our discoveries illuminate how one sequence has evolved to encode functions previously performed by two distinct proteins.

Ribosomal proteins can hold a more accurate record of bacterial thermal adaptation compared to rRNA

Ribosomal genes are widely used as 'molecular clocks' to infer evolutionary relationships between species. However, their utility as 'molecular thermometers' for estimating optimal growth temperature of microorganisms remains uncertain. Previously, some estimations were made using the nucleotide composition of ribosomal RNA (rRNA), but the universal application of this approach was hindered by numerous outliers. In this study, we aimed to address this problem by identifying additional indicators of thermal adaptation within the sequences of ribosomal proteins. By comparing sequences from 2021 bacteria with known optimal growth temperature, we identified novel indicators among the metal-binding residues of ribosomal proteins. We found that these residues serve as conserved adaptive features for bacteria thriving above 40°C, but not at lower temperatures. Furthermore, the presence of these metal-binding residues exhibited a stronger correlation with the optimal growth temperature of bacteria compared to the commonly used correlation with the 16S rRNA GC content. And an even more accurate correlation was observed between the optimal growth temperature and the YVIWREL amino acid content within ribosomal proteins. Overall, our work suggests that ribosomal proteins contain a more accurate record of bacterial thermal adaptation compared to rRNA. This finding may simplify the analysis of unculturable and extinct species.

In vitro allosteric transcription factor-based biosensing

A biosensor is an analytical device that converts a biological response into a measurable output signal. Bacterial allosteric transcription factors (aTFs) have been utilized as a novel class of recognition elements for in vitro biosensing, which circumvents the limitations of aTF-based whole-cell biosensors (WCBs) and helps to meet the increasing requirement of small-molecule biosensors for diverse applications. In this review, we summarize the recent advances related to the configuration of aTF-based biosensors in vitro. Particularly, we evaluate the advantages of aTFs for in vitro biosensing and highlight their great potential for the establishment of robust and easy-to-implement biosensing strategies. We argue that key technical innovations and generalizable workflows will enhance the pipeline for facile construction of diverse aTF-based small-molecule biosensors.

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.

Preparation and characterization of a laccase-like enzyme from Thermomicrobium roseum

In this study, a laccase-like gene from Thermomicrobium roseum DSM 5159 (TrLac-like) (NCBI: WP_012642205.1) was recombinantly expressed in Bacillus subtilis WB600. The optimum temperature and pH for TrLac-like were 50 °C and 6.0, respectively. TrLac-like showed high tolerance to mixed systems of water and organic solvents, indicating its potential for large-scale application in various industries. It showed 36.81 % similarity with YlmD from Geobacillus stearothermophilus (PDB:6T1B) in sequence alignment; therefore, 6T1B was employed as the template for homology modeling. To improve catalytic efficiency, amino acid substitutions within 5 Å of the inosine ligand were simulated to reduce the binding energy and promote substrate affinity. Single and double substitutions (44 and 18, respectively) were prepared, and the catalytic efficiency of the mutant A248D was increased to approximately 110-fold that of the wild type, while the thermal stability was maintained. Bioinformatics analysis revealed that the significant improvement in catalytic efficiency could be attributed to the formation of new hydrogen bonds between the enzyme and substrate. With a further decrease in the binding energy, the catalytic efficiency of the multiple mutant H129N/A248D was approximately 14-fold higher than that of the wild type but lower than that of the single mutant A248D. This is possibly because k_cat also decreased with the decrease of K_m; consequently, the substrate could not be released in time owing to the enzyme with the combination mutation not being able to release the substrate at a high rate.

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.

Efficient Exploration of Sequence Space by Sequence-Guided Protein Engineering and Design

The rapid growth of sequence databases over the past two decades means that protein engineers faced with optimizing a protein for any given task will often have immediate access to a vast number of related protein sequences. These sequences encode information about the evolutionary history of the protein and the underlying sequence requirements to produce folded, stable, and functional protein variants. Methods that can take advantage of this information are an increasingly important part of the protein engineering tool kit. In this Perspective, we discuss the utility of sequence data in protein engineering and design, focusing on recent advances in three main areas: the use of ancestral sequence reconstruction as an engineering tool to generate thermostable and multifunctional proteins, the use of sequence data to guide engineering of multipoint mutants by structure-based computational protein design, and the use of unlabeled sequence data for unsupervised and semisupervised machine learning, allowing the generation of diverse and functional protein sequences in unexplored regions of sequence space. Altogether, these methods enable the rapid exploration of sequence space within regions enriched with functional proteins and therefore have great potential for accelerating the engineering of stable, functional, and diverse proteins for industrial and biomedical applications.

Engineering the Thermostability of a d-Carbamoylase Based on Ancestral Sequence Reconstruction for the Efficient Synthesis of d-Tryptophan

Employing ancestral sequence reconstruction and consensus sequence analysis, the thermostability of a novel d-carbamoylase derived from Nitratireductor indicus (NiHyuC) was engineered through greedy-oriented iterative combinatorial mutagenesis. A mutant S202P/E208D/R277L (M4Th3) was obtained with significantly elevated thermostability. M4Th3 has a half-life of 36.5 h at 40 °C, about 28.5 times of 1.3 h of its parent M4. For the reaction at 40 °C, M4Th3 can catalyze 10 mM N-carbamoyl-d-tryptophan to produce d-tryptophan with a conversion ratio of 96.4% after 12 h, which is significantly higher than 64.1% of M4. MD simulation reveals that new hydrogen bonds emerging from E208D on the surface can increase the hydrophobicity of the protein, leading to improved stability. More importantly, R277L could contribute to enhanced interface stability of homodimeric M4. This study provides a thermostable d-carbamoylase for the "hydantoinase process", which has potential in the industrial synthesis of optically pure natural and non-natural amino acids.

Engineering functional thermostable proteins using ancestral sequence reconstruction

Natural proteins are often only slightly more stable in the native state than the denatured state, and an increase in environmental temperature can easily shift the balance toward unfolding. Therefore, the engineering of proteins to improve protein stability is an area of intensive research. Thermostable proteins are required to withstand industrial process conditions, for increased shelf-life of protein therapeutics, for developing robust 'biobricks' for synthetic biology applications, and for research purposes (e.g., structure determination). In addition, thermostability buffers the often destabilizing effects of mutations introduced to improve other properties. Rational design approaches to engineering thermostability require structural information, but even with advanced computational methods, it is challenging to predict or parameterize all the relevant structural factors with sufficient precision to anticipate the results of a given mutation. Directed evolution is an alternative when structures are unavailable but requires extensive screening of mutant libraries. Recently, however, bioinspired approaches based on phylogenetic analyses have shown great promise. Leveraging the rapid expansion in sequence data and bioinformatic tools, ancestral sequence reconstruction can generate highly stable folds for novel applications in industrial chemistry, medicine, and synthetic biology. This review provides an overview of the factors important for successful inference of thermostable proteins by ancestral sequence reconstruction and what it can reveal about the determinants of stability in proteins.

Efficient Base-Catalyzed Kemp Elimination in an Engineered Ancestral Enzyme

The routine generation of enzymes with completely new active sites is a major unsolved problem in protein engineering. Advances in this field have thus far been modest, perhaps due, at least in part, to the widespread use of modern natural proteins as scaffolds for de novo engineering. Most modern proteins are highly evolved and specialized and, consequently, difficult to repurpose for completely new functionalities. Conceivably, resurrected ancestral proteins with the biophysical properties that promote evolvability, such as high stability and conformational diversity, could provide better scaffolds for de novo enzyme generation. Kemp elimination, a non-natural reaction that provides a simple model of proton abstraction from carbon, has been extensively used as a benchmark in de novo enzyme engineering. Here, we present an engineered ancestral β-lactamase with a new active site that is capable of efficiently catalyzing Kemp elimination. The engineering of our Kemp eliminase involved minimalist design based on a single function-generating mutation, inclusion of an extra polypeptide segment at a position close to the de novo active site, and sharply focused, low-throughput library screening. Nevertheless, its catalytic parameters (kcat/KM~2·105 M-1 s-1, kcat~635 s-1) compare favorably with the average modern natural enzyme and match the best proton-abstraction de novo Kemp eliminases that are reported in the literature. The general implications of our results for de novo enzyme engineering are discussed.

Empirical comparison and analysis of machine learning-based predictors for predicting and analyzing of thermophilic proteins

Thermophilic proteins (TPPs) are critical for basic research and in the food industry due to their ability to maintain a thermodynamically stable fold at extremely high temperatures. Thus, the expeditious identification of novel TPPs through computational models from protein sequences is very desirable. Over the last few decades, a number of computational methods, especially machine learning (ML)-based methods, for in silico prediction of TPPs have been developed. Therefore, it is desirable to revisit these methods and summarize their advantages and disadvantages in order to further develop new computational approaches to achieve more accurate and improved prediction of TPPs. With this goal in mind, we comprehensively investigate a large collection of fourteen state-of-the-art TPP predictors in terms of their dataset size, feature encoding schemes, feature selection strategies, ML algorithms, evaluation strategies and web server/software usability. To the best of our knowledge, this article represents the first comprehensive review on the development of ML-based methods for in silico prediction of TPPs. Among these TPP predictors, they can be classified into two groups according to the interpretability of ML algorithms employed (i.e., computational black-box methods and computational white-box methods). In order to perform the comparative analysis, we conducted a comparative study on several currently available TPP predictors based on two benchmark datasets. Finally, we provide future perspectives for the design and development of new computational models for TPP prediction. We hope that this comprehensive review will facilitate researchers in selecting an appropriate TPP predictor that is the most suitable one to deal with their purposes and provide useful perspectives for the development of more effective and accurate TPP predictors.

Ancestral Sequence Reconstruction for Exploring Alkaloid Evolution

The complex and bioactive monoterpene indole alkaloids (MIAs) found in Catharanthus roseus and related species are the products of many millions of years of evolution through mutation and natural selection. Ancestral sequence reconstruction (ASR) is a method that combines phylogenetic analysis and experimental biochemistry to infer details about past events in protein evolution. Here, I propose that ASR could be leveraged to understand how enzymes catalyzing the formation of complex alkaloids arose over evolutionary time. I discuss the steps of ASR, including sequence selection, multiple sequence alignment, tree inference, and the generation and characterization of inferred ancestral enzymes.

Mechanochemical Evolution of Disulfide Bonds in Proteins

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

Ancestral Sequence Reconstruction and Alternate Amino Acid States Guide Protein Library Design for Directed Evolution

Engineered proteins possess nearly limitless possibilities in medical and industrial applications but finding a precise amino acid sequence for these applications is challenging. A robust approach for discovering protein sequences with a desired functionality uses a library design method in which combinations of mutations are applied to a robust starting point. Determining useful mutations can be tortuous, yet rewarding; in this chapter, we present a novel library design method that uses information provided by ancestral sequence reconstruction (ASR) to create a library likely to have stable proteins with diverse function. ASR computational tools use a multi-sequence alignment of homologous proteins and an evolutionary model to estimate the protein sequences of the numerous common ancestors. For all ancestors, these tools calculate the probability of every amino acid occurring at each position within the sequence alignment. The alternate amino acid states at individual positions corelate to a region of stability in sequence space around the ancestral sequence which can inform site-wise diversification within a combinatorial library. The method presented in this chapter balances the quality of results, the computational resources needed, and ease of use.

An Integrated Method to Reconstruct Ancient Proteins

Proteins have played a fundamental role throughout life's history on Earth. Despite their biological importance, ancient origin, early function, and evolution of proteins are seldom able to be directly studied because few of these attributes are preserved across geologic timescales. Ancestral sequence reconstruction (ASR) provides a method to infer ancestral amino acid sequences and determine the evolutionary predecessors of modern-day proteins using phylogenetic tools. Laboratory application of ASR allows ancient sequences to be deduced from genetic information available in extant organisms and then experimentally resurrected to elucidate ancestral characteristics. In this article, we provide a generalized, stepwise protocol that considers the major elements of a well-designed ASR study and details potential sources of reconstruction bias that can reduce the relevance of historical inferences. We underscore key stages in our approach so that it may be broadly utilized to reconstruct the evolutionary histories of proteins.

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

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

A novel sequence-based predictor for identifying and characterizing thermophilic proteins using estimated propensity scores of dipeptides

Owing to their ability to maintain a thermodynamically stable fold at extremely high temperatures, thermophilic proteins (TTPs) play a critical role in basic research and a variety of applications in the food industry. As a result, the development of computation models for rapidly and accurately identifying novel TTPs from a large number of uncharacterized protein sequences is desirable. In spite of existing computational models that have already been developed for characterizing thermophilic proteins, their performance and interpretability remain unsatisfactory. We present a novel sequence-based thermophilic protein predictor, termed SCMTPP, for improving model predictability and interpretability. First, an up-to-date and high-quality dataset consisting of 1853 TPPs and 3233 non-TPPs was compiled from published literature. Second, the SCMTPP predictor was created by combining the scoring card method (SCM) with estimated propensity scores of g-gap dipeptides. Benchmarking experiments revealed that SCMTPP had a cross-validation accuracy of 0.883, which was comparable to that of a support vector machine-based predictor (0.906-0.910) and 2-17% higher than that of commonly used machine learning models. Furthermore, SCMTPP outperformed the state-of-the-art approach (ThermoPred) on the independent test dataset, with accuracy and MCC of 0.865 and 0.731, respectively. Finally, the SCMTPP-derived propensity scores were used to elucidate the critical physicochemical properties for protein thermostability enhancement. In terms of interpretability and generalizability, comparative results showed that SCMTPP was effective for identifying and characterizing TPPs. We had implemented the proposed predictor as a user-friendly online web server at http://pmlabstack.pythonanywhere.com/SCMTPP in order to allow easy access to the model. SCMTPP is expected to be a powerful tool for facilitating community-wide efforts to identify TPPs on a large scale and guiding experimental characterization of TPPs.

Precambrian Paleobiology: Precedents, Progress, and Prospects

In 1859, C. R. Darwin highlighted the “inexplicable” absence of evidence of life prior to the beginning of the Cambrian. Given this lack of evidence and the natural rather than theological unfolding of life’s development Darwin espoused, over the following 50 years his newly minted theory was disputed. At the turn of the 19th century, beginning with the discoveries of C. D. Walcott, glimmerings of the previously “unknown and unknowable” early fossil record came to light – but Walcott’s Precambrian finds were also discounted. It was not until the breakthrough advances of the 1950’s and the identification of modern stromatolites (1956), Precambrian phytoplankton in shales (1950’s), stromatolitic microbes in cherts (1953), and terminal-Precambrian soft-bodied animal fossils (1950’s) that the field was placed on firm footing. Over the following half-century, the development and application of new analytical techniques coupled with the groundbreaking contributions of the Precambrian Paleobiology Research Group spurred the field to its international and distinctly interdisciplinary status. Significant progress has been made worldwide. Among these advances, the known fossil record has been extended sevenfold (from ∼0.5 to ∼3.5 Ga); the fossil record has been shown consistent with rRNA phylogenies (adding credence to both); and the timing and evolutionary significance of an increase of environmental oxygen (∼2.3 Ga), of eukaryotic organisms (∼2.0 Ga), and of evolution-speeding and biota-diversifying eukaryotic sexual reproduction (∼1.2 Ga) have been identified. Nevertheless, much remains to be learned. Such major unsolved problems include the absence of definitive evidence of the widely assumed life-generating “primordial soup”; the timing of the origin of oxygenic photosynthesis; the veracity of postulated changes in global photic-zone temperature from 3.5 Ga to the present; the bases of the advent of eukaryotic sexuality-requiring gametogenesis and syngamy; and the timing of origin and affinities of the small soft-bodied precursors of the Ediacaran Fauna.

Were Ancestral Proteins Less Specific?

Some have hypothesized that ancestral proteins were, on average, less specific than their descendants. If true, this would provide a universal axis along which to organize protein evolution and suggests that reconstructed ancestral proteins may be uniquely powerful tools for protein engineering. Ancestral sequence reconstruction studies are one line of evidence used to support this hypothesis. Previously, we performed such a study, investigating the evolution of peptide-binding specificity for the paralogs S100A5 and S100A6. The modern proteins appeared more specific than their last common ancestor (ancA5/A6), as each paralog bound a subset of the peptides bound by ancA5/A6. In this study, we revisit this transition, using quantitative phage display to measure the interactions of 30,533 random peptides with human S100A5, S100A6, and ancA5/A6. This unbiased screen reveals a different picture. While S100A5 and S100A6 do indeed bind to a subset of the peptides recognized by ancA5/A6, they also acquired new peptide partners outside of the set recognized by ancA5/A6. Our previous work showed that ancA5/A6 had lower specificity than its descendants when measured against biological targets; our new work shows that ancA5/A6 has similar specificity to the modern proteins when measured against a random set of peptide targets. This demonstrates that altered biological specificity does not necessarily indicate altered intrinsic specificity, and sounds a cautionary note for using ancestral reconstruction studies with biological targets as a means to infer global evolutionary trends in specificity.

The Mutational Robustness of the Genetic Code and Codon Usage in Environmental Context: A Non-Extremophilic Preference?

The genetic code was evolved, to some extent, to minimize the effects of mutations. The effects of mutations depend on the amino acid repertoire, the structure of the genetic code and frequencies of amino acids in proteomes. The amino acid compositions of proteins and corresponding codon usages are still under selection, which allows us to ask what kind of environment the standard genetic code is adapted to. Using simple computational models and comprehensive datasets comprising genomic and environmental data from all three domains of Life, we estimate the expected severity of non-synonymous genomic mutations in proteins, measured by the change in amino acid physicochemical properties. We show that the fidelity in these physicochemical properties is expected to deteriorate with extremophilic codon usages, especially in thermophiles. These findings suggest that the genetic code performs better under non-extremophilic conditions, which not only explains the low substitution rates encountered in halophiles and thermophiles but the revealed relationship between the genetic code and habitat allows us to ponder on earlier phases in the history of Life.

The Molecular Determinants of Thermoadaptation: Methanococcales as a Case Study

Previous reports have shown that environmental temperature impacts proteome evolution in Bacteria and Archaea. However, it is unknown whether thermoadaptation mainly occurs via the sequential accumulation of substitutions, massive horizontal gene transfers, or both. Measuring the real contribution of amino acid substitution to thermoadaptation is challenging, because of confounding environmental and genetic factors (e.g., pH, salinity, genomic G + C content) that also affect proteome evolution. Here, using Methanococcales, a major archaeal lineage, as a study model, we show that optimal growth temperature is the major factor affecting variations in amino acid frequencies of proteomes. By combining phylogenomic and ancestral sequence reconstruction approaches, we disclose a sequential substitutional scheme in which lysine plays a central role by fine tuning the pool of arginine, serine, threonine, glutamine, and asparagine, whose frequencies are strongly correlated with optimal growth temperature. Finally, we show that colonization to new thermal niches is not associated with high amounts of horizontal gene transfers. Altogether, although the acquisition of a few key proteins through horizontal gene transfer may have favored thermoadaptation in Methanococcales, our findings support sequential amino acid substitutions as the main factor driving thermoadaptation.

Revenant: a database of resurrected proteins

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

Evolution of bacterial steroid biosynthesis and its impact on eukaryogenesis

Steroids are one of three major lipid components of the eukaryotic cellular membrane, along with glycerophospolipids and sphingolipids. Steroids have critical roles in eukaryotic endocytosis and thus may have been structural prerequisites for the endocytic acquisition of mitochondria during eukaryogenesis. The evolutionary history of the eukaryotic cellular membrane is poorly understood and, as such, has limited our understanding of eukaryogenesis. We address the evolution of steroid biosynthesis by combining ancestral sequence reconstruction and phylogenetic analyses of steroid biosynthesis genes. Our results indicate that steroid biosynthesis evolved within bacteria in response to the rise of oxygen and was later horizontally transferred to eukaryotes. Membrane properties of early eukaryotes are inferred to have been different than that of modern eukaryotes.

Steroids are components of the eukaryotic cellular membrane and have indispensable roles in the process of eukaryotic endocytosis by regulating membrane fluidity and permeability. In particular, steroids may have been a structural prerequisite for the acquisition of mitochondria via endocytosis during eukaryogenesis. While eukaryotes are inferred to have evolved from an archaeal lineage, there is little similarity between the eukaryotic and archaeal cellular membranes. As such, the evolution of eukaryotic cellular membranes has limited our understanding of eukaryogenesis. Despite evolving from archaea, the eukaryotic cellular membrane is essentially a fatty acid bacterial-type membrane, which implies a substantial bacterial contribution to the evolution of the eukaryotic cellular membrane. Here, we address the evolution of steroid biosynthesis in eukaryotes by combining ancestral sequence reconstruction and comprehensive phylogenetic analyses of steroid biosynthesis genes. Contrary to the traditional assumption that eukaryotic steroid biosynthesis evolved within eukaryotes, most steroid biosynthesis genes are inferred to be derived from bacteria. In particular, aerobic deltaproteobacteria (myxobacteria) seem to have mediated the transfer of key genes for steroid biosynthesis to eukaryotes. Analyses of resurrected steroid biosynthesis enzymes suggest that the steroid biosynthesis pathway in early eukaryotes may have been similar to the pathway seen in modern plants and algae. These resurrected proteins also experimentally demonstrate that molecular oxygen was required to establish the modern eukaryotic cellular membrane during eukaryogenesis. Our study provides unique insight into relationships between early eukaryotes and other bacteria in addition to the well-known endosymbiosis with alphaproteobacteria.

Fully Automated Ancestral Sequence Reconstruction using FireProtASR

Protein evolution and protein engineering techniques are of great interest in basic science and industrial applications such as pharmacology, medicine, or biotechnology. Ancestral sequence reconstruction (ASR) is a powerful technique for probing evolutionary relationships and engineering robust proteins with good thermostability and broad substrate specificity. The following protocol describes the setting up and execution of an automated FireProt^ASR workflow using a dedicated web site. The service allows for inference of ancestral proteins automatically, from a single protein sequence. Once a protein sequence is submitted, the server will build a dataset of homology sequences, perform a multiple sequence alignment (MSA), build a phylogenetic tree, and reconstruct ancestral nodes. The protocol is also highly flexible and allows for multiple forms of input, advanced settings, and the ability to start jobs from: (i) a single sequence, (ii) a set of homologous sequences, (iii) an MSA, and (iv) a phylogenetic tree. This approach automates all necessary steps and offers a way for novices with limited exposure to ASR techniques to improve the properties of a protein of interest. The technique can even be used to introduce catalytic promiscuity into an enzyme. A web server for accessing the fully automated workflow is freely accessible at https://loschmidt.chemi.muni.cz/fireprotasr/. © 2021 Wiley Periodicals LLC. Basic Protocol: ASR using the Web Server FireProt^ASR.

Molecular coevolution of coagulation factor VIII and von Willebrand factor

Ancestral sequence reconstruction provides a unique platform for investigating the molecular evolution of single gene products and recently has shown success in engineering advanced biological therapeutics. To date, the coevolution of proteins within complexes and protein-protein interactions is mostly investigated in silico via proteomics and/or within single-celled systems. Herein, ancestral sequence reconstruction is used to investigate the molecular evolution of 2 proteins linked not only by stabilizing association in circulation but also by their independent roles within the primary and secondary hemostatic systems of mammals. Using sequence analysis and biochemical characterization of recombinant ancestral von Willebrand factor (VWF) and coagulation factor VIII (FVIII), we investigated the evolution of the essential macromolecular FVIII/VWF complex. Our data support the hypothesis that these coagulation proteins coevolved throughout mammalian diversification, maintaining strong binding affinities while modulating independent and distinct hemostatic activities in diverse lineages.

Improving the Thermostability of Xylanase A from Bacillus subtilis by Combining Bioinformatics and Electrostatic Interactions Optimization

The rational improvement of the enzyme catalytic activity is one of the most significant challenges in biotechnology. Most conventional strategies used to engineer enzymes involve selecting mutations to increase their thermostability. Determining good criteria for choosing these substitutions continues to be a challenge. In this work, we combine bioinformatics, electrostatic analysis, and molecular dynamics to predict beneficial mutations that may improve the thermostability of XynA from Bacillus subtilis. First, the Tanford-Kirkwood surface accessibility method is used to characterize each ionizable residue contribution to the protein native state stability. Residues identified to be destabilizing were mutated with the corresponding residues determined by the consensus or ancestral sequences at the same locations. Five mutants (K99T/N151D, K99T, S31R, N151D, and K154A) were investigated and compared with 12 control mutants derived from experimental approaches from the literature. Molecular dynamics results show that the mutants exhibited folding temperatures in the order K99T > K99T/N151D > S31R > N151D > WT > K154A. The combined approaches employed provide an effective strategy for low-cost enzyme optimization needed for large-scale biotechnological and medical applications.

A meta-analysis of the activity, stability, and mutational characteristics of temperature-adapted enzymes

Understanding the characteristics that define temperature-adapted enzymes has been a major goal of extremophile enzymology in recent decades. In the present study, we explore these characteristics by comparing psychrophilic, mesophilic, and thermophilic enzymes. Through a meta-analysis of existing data, we show that psychrophilic enzymes exhibit a significantly larger gap (T_g) between their optimum and melting temperatures compared with mesophilic and thermophilic enzymes. These results suggest that T_g may be a useful indicator as to whether an enzyme is psychrophilic or not and that models of psychrophilic enzyme catalysis need to account for this gap. Additionally, by using predictive protein stability software, HoTMuSiC and PoPMuSiC, we show that the deleterious nature of amino acid substitutions to protein stability increases from psychrophiles to thermophiles. How this ultimately affects the mutational tolerance and evolutionary rate of temperature adapted organisms is currently unknown.

Kinetic Analysis Suggests Evolution of Ribosome Specificity in Modern Elongation Factor-Tus from "Generalist" Ancestors

It has been hypothesized that early enzymes are more promiscuous than their extant orthologs. Whether or not this hypothesis applies to the translation machinery, the oldest molecular machine of life, is not known. Efficient protein synthesis relies on a cascade of specific interactions between the ribosome and the translation factors. Here, using elongation factor-Tu (EF-Tu) as a model system, we have explored the evolution of ribosome specificity in translation factors. Employing presteady state fast kinetics using quench flow, we have quantitatively characterized the specificity of two sequence-reconstructed 1.3- to 3.3-Gy-old ancestral EF-Tus toward two unrelated bacterial ribosomes, mesophilic Escherichia coli and thermophilic Thermus thermophilus. Although the modern EF-Tus show clear preference for their respective ribosomes, the ancestral EF-Tus show similar specificity for diverse ribosomes. In addition, despite increase in the catalytic activity with temperature, the ribosome specificity of the thermophilic EF-Tus remains virtually unchanged. Our kinetic analysis thus suggests that EF-Tu proteins likely evolved from the catalytically promiscuous, “generalist” ancestors. Furthermore, compatibility of diverse ribosomes with the modern and ancestral EF-Tus suggests that the ribosomal core probably evolved before the diversification of the EF-Tus. This study thus provides important insights regarding the evolution of modern translation machinery.

Ancestral sequence reconstruction - An underused approach to understand the evolution of gene function in plants?

Whilst substantial research effort has been placed on understanding the interactions of plant proteins with their molecular partners, relatively few studies in plants - by contrast to work in other organisms - address how these interactions evolve. It is thought that ancestral proteins were more promiscuous than modern proteins and that specificity often evolved following gene duplication and subsequent functional refining. However, ancestral protein resurrection studies have found that some modern proteins have evolved de novo from ancestors lacking those functions. Intriguingly, the new interactions evolved as a consequence of just a few mutations and, as such, acquisition of new functions appears to be neither difficult nor rare, however, only a few of them are incorporated into biological processes before they are lost to subsequent mutations. Here, we detail the approach of ancestral sequence reconstruction (ASR), providing a primer to reconstruct the sequence of an ancestral gene. We will present case studies from a range of different eukaryotes before discussing the few instances where ancestral reconstructions have been used in plants. As ASR is used to dig into the remote evolutionary past, we will also present some alternative genetic approaches to investigate molecular evolution on shorter timescales. We argue that the study of plant secondary metabolism is particularly well suited for ancestral reconstruction studies. Indeed, its ancient evolutionary roots and highly diverse landscape provide an ideal context in which to address the focal issue around the emergence of evolutionary novelties and how this affects the chemical diversification of plant metabolism.

Convergent Biochemical Pathways for Xanthine Alkaloid Production in Plants Evolved from Ancestral Enzymes with Different Catalytic Properties

Convergent evolution is widespread but the extent to which common ancestral conditions are necessary to facilitate the independent acquisition of similar traits remains unclear. In order to better understand how ancestral biosynthetic catalytic capabilities might lead to convergent evolution of similar modern-day biochemical pathways, we resurrected ancient enzymes of the caffeine synthase (CS) methyltransferases that are responsible for theobromine and caffeine production in flowering plants. Ancestral CS enzymes of Theobroma, Paullinia, and Camellia exhibited similar substrate preferences but these resulted in the formation of different sets of products. From these ancestral enzymes, descendants with similar substrate preference and product formation independently evolved after gene duplication events in Theobroma and Paullinia. Thus, it appears that the convergent modern-day pathways likely originated from ancestral pathways with different inferred flux. Subsequently, the modern-day enzymes originated independently via gene duplication and their convergent catalytic characteristics evolved to partition the multiple ancestral activities by different mutations that occurred in homologous regions of the ancestral proteins. These results show that even when modern-day pathways and recruited genes are similar, the antecedent conditions may be distinctive such that different evolutionary steps are required to generate convergence.

Evolution of frustrated and stabilising contacts in reconstructed ancient proteins

Energetic properties of a protein are a major determinant of its evolutionary fitness. Using a reconstruction algorithm, dating the reconstructed proteins and calculating the interaction network between their amino acids through a coevolutionary approach, we studied how the interactions that stabilise 890 proteins, belonging to five families, evolved for billions of years. In particular, we focused our attention on the network of most strongly attractive contacts and on that of poorly optimised, frustrated contacts. Our results support the idea that the cluster of most attractive interactions extends its size along evolutionary time, but from the data, we cannot conclude that protein stability or that the degree of frustration tends always to decrease.

Ancestral Sequence Reconstruction: From Chemical Paleogenetics to Maximum Likelihood Algorithms and Beyond

As both a computational and an experimental endeavor, ancestral sequence reconstruction remains a timely and important technique. Modern approaches to conduct ancestral sequence reconstruction for proteins are built upon a conceptual framework from journal founder Emile Zuckerkandl. On top of this, work on maximum likelihood phylogenetics published in Journal of Molecular Evolution in 1996 was one of the first approaches for generating maximum likelihood ancestral sequences of proteins. From its computational history, future model development needs as well as potential applications in areas as diverse as computational systems biology, molecular community ecology, infectious disease therapeutics and other biomedical applications, and biotechnology are discussed. From its past in this journal, there is a bright future for ancestral sequence reconstruction in the field of evolutionary biology.

Divergent Evolution of a Protein-Protein Interaction Revealed through Ancestral Sequence Reconstruction and Resurrection

The postsynaptic density extends across the postsynaptic dendritic spine with discs large (DLG) as the most abundant scaffolding protein. DLG dynamically alters the structure of the postsynaptic density, thus controlling the function and distribution of specific receptors at the synapse. DLG contains three PDZ domains and one important interaction governing postsynaptic architecture is that between the PDZ3 domain from DLG and a protein called cysteine-rich interactor of PDZ3 (CRIPT). However, little is known regarding functional evolution of the PDZ3:CRIPT interaction. Here, we subjected PDZ3 and CRIPT to ancestral sequence reconstruction, resurrection, and biophysical experiments. We show that the PDZ3:CRIPT interaction is an ancient interaction, which was likely present in the last common ancestor of Eukaryotes, and that high affinity is maintained in most extant animal phyla. However, affinity is low in nematodes and insects, raising questions about the physiological function of the interaction in species from these animal groups. Our findings demonstrate how an apparently established protein-protein interaction involved in cellular scaffolding in bilaterians can suddenly be subject to dynamic evolution including possible loss of function.

Evolution of Protein Structure and Stability in Global Warming

This review focuses on the molecular signatures of protein structures in relation to evolution and survival in global warming. It is based on the premise that the power of evolutionary selection may lead to thermotolerant organisms that will repopulate the planet and continue life in general, but perhaps with different kinds of flora and fauna. Our focus is on molecular mechanisms, whereby known examples of thermoresistance and their physicochemical characteristics were noted. A comparison of interactions of diverse residues in proteins from thermophilic and mesophilic organisms, as well as reverse genetic studies, revealed a set of imprecise molecular signatures that pointed to major roles of hydrophobicity, solvent accessibility, disulfide bonds, hydrogen bonds, ionic and π-electron interactions, and an overall condensed packing of the higher-order structure, especially in the hydrophobic regions. Regardless of mutations, specialized protein chaperones may play a cardinal role. In evolutionary terms, thermoresistance to global warming will likely occur in stepwise mutational changes, conforming to the molecular signatures, such that each "intermediate" fits a temporary niche through punctuated equilibrium, while maintaining protein functionality. Finally, the population response of different species to global warming may vary substantially, and, as such, some may evolve while others will undergo catastrophic mass extinction.

The Inflection Point Hypothesis: The Relationship between the Temperature Dependence of Enzyme-Catalyzed Reaction Rates and Microbial Growth Rates

The temperature dependence of biological rates at different scales (from individual enzymes to isolated organisms to ecosystem processes such as soil respiration and photosynthesis) is the subject of much historical and contemporary research. The precise relationship between the temperature dependence of enzyme rates and those at larger scales is not well understood. We have developed macromolecular rate theory (MMRT) to describe the temperature dependence of biological processes at all scales. Here we formalize the scaling relationship by investigating MMRT both at the molecular scale (constituent enzymes) and for growth of the parent organism. We demonstrate that the inflection point (T_inf) for the temperature dependence of individual metabolic enzymes coincides with the optimal growth temperature for the parent organism, and we rationalize this concordance in terms of the necessity for linearly correlated rates for metabolic enzymes over fluctuating environmental temperatures to maintain homeostasis. Indeed, T_inf is likely to be under strong selection pressure to maintain coordinated rates across environmental temperature ranges. At temperatures at which rates become uncorrelated, we postulate a regulatory catastrophe and organism growth rates precipitously decline at temperatures where this occurs. We show that the curvature in the plots of the natural log of the rate versus temperature for individual enzymes determines the curvature for the metabolic process overall and the curvature for the temperature dependence of the growth of the organism. We have called this "the inflection point hypothesis", and this hypothesis suggests many avenues for future investigation, including avenues for engineering the thermal tolerance of organisms.

Ancestral sequence reconstruction produces thermally stable enzymes with mesophilic enzyme-like catalytic properties

Enzymes have high catalytic efficiency and low environmental impact, and are therefore potentially useful tools for various industrial processes. Crucially, however, natural enzymes do not always have the properties required for specific processes. It may be necessary, therefore, to design, engineer, and evolve enzymes with properties that are not found in natural enzymes. In particular, the creation of enzymes that are thermally stable and catalytically active at low temperature is desirable for processes involving both high and low temperatures. In the current study, we designed two ancestral sequences of 3-isopropylmalate dehydrogenase by an ancestral sequence reconstruction technique based on a phylogenetic analysis of extant homologous amino acid sequences. Genes encoding the designed sequences were artificially synthesized and expressed in Escherichia coli. The reconstructed enzymes were found to be slightly more thermally stable than the extant thermophilic homologue from Thermus thermophilus. Moreover, they had considerably higher low-temperature catalytic activity as compared with the T. thermophilus enzyme. Detailed analyses of their temperature-dependent specific activities and kinetic properties showed that the reconstructed enzymes have catalytic properties similar to those of mesophilic homologues. Collectively, our study demonstrates that ancestral sequence reconstruction can produce a thermally stable enzyme with catalytic properties adapted to low-temperature reactions.