The thermostability and specificity of ancient proteins

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.

A conserved π-helix plays a key role in thermoadaptation of catalysis in the glycoside hydrolase family 4

Here, we characterize the role of a π-helix in the molecular mechanisms underlying thermoadaptation in the glycoside hydrolase family 4 (GH4). The interspersed π-helix present in a subgroup is evolutionarily related to a conserved α-helix in other orthologs by a single residue insertion/deletion event. The insertional residue, Phe407, in a hyperthermophilic α-glucuronidase, makes specific interactions across the inter-subunit interface. In order to establish the sequence-structure-stability implications of the π-helix, the wild-type and the deletion variant (Δ407) were characterized. The variant showed a significant lowering of melting temperature and optimum temperature for the highest activity. Crystal structures of the proteins show a transformation of the π-helix to a continuous α-helix in the variant, identical to that in orthologs lacking this insertion. Thermodynamic parameters were determined from stability curves representing the temperature dependence of unfolding free energy. Though the proteins display maximum stabilities at similar temperatures, a higher melting temperature in the wild-type is achieved by a combination of higher enthalpy and lower heat capacity of unfolding. Comparisons of the structural changes, and the activity and thermodynamic profiles allow us to infer that specific non-covalent interactions, and the existence of residual structure in the unfolded state, are crucial determinants of its thermostability. These features permit the enzyme to balance the preservation of structure at a higher temperature with the thermodynamic stability required for optimum catalysis.

The use of consensus sequence information to engineer stability and activity in proteins

The goal of protein design is to create proteins that are stable, soluble, and active. Here we focus on one approach to protein design in which sequence information is used to create a "consensus" sequence. Such consensus sequences comprise the most common residue at each position in a multiple sequence alignment (MSA). After describing some general ideas that relate MSA and consensus sequences and presenting a statistical thermodynamic framework that relates consensus and non-consensus sequences to stability, we detail the process of designing a consensus sequence and survey reports of consensus design and characterization from the literature. Many of these consensus proteins retain native biological activities including ligand binding and enzyme activity. Remarkably, in most cases the consensus protein shows significantly higher stability than extant versions of the protein, as measured by thermal or chemical denaturation, consistent with the statistical thermodynamic model. To understand this stability increase, we compare various features of consensus sequences with the extant MSA sequences from which they were derived. Consensus sequences show enrichment in charged residues (most notably glutamate and lysine) and depletion of uncharged polar residues (glutamine, serine, and asparagine). Surprisingly, a survey of stability changes resulting from point substitutions show little correlation with residue frequencies at the corresponding positions within the MSA, suggesting that the high stability of consensus proteins may result from interactions among residue pairs or higher-order clusters. Whatever the source, the large number of reported successes demonstrates that consensus design is a viable route to generating active and in many cases highly stabilized proteins.

Resurrection of ancestral effector caspases identifies novel networks for evolution of substrate specificity

Apoptotic caspases evolved with metazoans more than 950 million years ago (MYA), and a series of gene duplications resulted in two subfamilies consisting of initiator and effector caspases. The effector caspase genes (caspases-3, -6, and -7) were subsequently fixed into the Chordata phylum more than 650 MYA when the gene for a common ancestor (CA) duplicated, and the three effector caspases have persisted throughout mammalian evolution. All caspases prefer an aspartate residue at the P1 position of substrates, so each caspase evolved discrete cellular roles through changes in substrate recognition at the P4 position combined with allosteric regulation. We examined the evolution of substrate specificity in caspase-6, which prefers valine at the P4 residue, compared with caspases-3 and -7, which prefer aspartate, by reconstructing the CA of effector caspases (AncCP-Ef1) and the CA of caspase-6 (AncCP-6An). We show that AncCP-Ef1 is a promiscuous enzyme with little distinction between Asp, Val, or Leu at P4. The specificity of caspase-6 was defined early in its evolution, where AncCP-6An demonstrates a preference for Val over Asp at P4. Structures of AncCP-Ef1 and of AncCP-6An show a network of charged amino acids near the S4 pocket that, when combined with repositioning a flexible active site loop, resulted in a more hydrophobic binding pocket in AncCP-6An. The ancestral protein reconstructions show that the caspase-hemoglobinase fold has been conserved for over 650 million years and that only three substitutions in the scaffold are necessary to shift substrate selection toward Val over Asp.

Larger active site in an ancestral hydroxynitrile lyase increases catalytically promiscuous esterase activity

Hydroxynitrile lyases (HNL's) belonging to the α/β-hydrolase-fold superfamily evolved from esterases approximately 100 million years ago. Reconstruction of an ancestral hydroxynitrile lyase in the α/β-hydrolase fold superfamily yielded a catalytically active hydroxynitrile lyase, HNL1. Several properties of HNL1 differ from the modern HNL from rubber tree (HbHNL). HNL1 favors larger substrates as compared to HbHNL, is two-fold more catalytically promiscuous for ester hydrolysis (p-nitrophenyl acetate) as compared to mandelonitrile cleavage, and resists irreversible heat inactivation to 35 °C higher than for HbHNL. We hypothesized that the x-ray crystal structure of HNL1 may reveal the molecular basis for the differences in these properties. The x-ray crystal structure solved to 1.96-Å resolution shows the expected α/β-hydrolase fold, but a 60% larger active site as compared to HbHNL. This larger active site echoes its evolution from esterases since related esterase SABP2 from tobacco also has a 38% larger active site than HbHNL. The larger active site in HNL1 likely accounts for its ability to accept larger hydroxynitrile substrates. Site-directed mutagenesis of HbHNL to expand the active site increased its promiscuous esterase activity 50-fold, consistent with the larger active site in HNL1 being the primary cause of its promiscuous esterase activity. Urea-induced unfolding of HNL1 indicates that it unfolds less completely than HbHNL (m-value = 0.63 for HNL1 vs 0.93 kcal/mol·M for HbHNL), which may account for the ability of HNL1 to better resist irreversible inactivation upon heating. The structure of HNL1 shows changes in hydrogen bond networks that may stabilize regions of the folded structure.

Larger active site in an ancestral hydroxynitrile lyase increases catalytically promiscuous esterase activity.. [DOI: 10.1101/2020.04.06.027797] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

Hydroxynitrile lyases (HNL’s) belonging to the α/β-hydrolase-fold superfamily evolved from esterases approximately 100 million years ago. Reconstruction of an ancestral hydroxynitrile lyase in the α/β-hydrolase fold superfamily yielded a catalytically active hydroxynitrile lyase, HNL1. Several properties of HNL1 differ from the modern HNL from rubber tree (HbHNL). HNL1 favors larger substrates as compared to HbHNL, is two-fold more catalytically promiscuous for ester hydrolysis (p-nitrophenyl acetate) as compared to mandelonitrile cleavage, and resists irreversible heat inactivation to 35 °C higher than for HbHNL. We hypothesized that the x-ray crystal structure of HNL1 may reveal the molecular basis for the differences in these properties. The x-ray crystal structure solved to 1.96-Å resolution shows the expected α/β-hydrolase fold, but a 60% larger active site as compared to HbHNL. This larger active site echoes its evolution from esterases since related esterase SABP2 from tobacco also has a 38% larger active site than HbHNL. The larger active site in HNL1 likely accounts for its ability to accept larger hydroxynitrile substrates. Site-directed mutagenesis of HbHNL to expand the active site increased its promiscuous esterase activity 50-fold, consistent with the larger active site in HNL1 being the primary cause of its promiscuous esterase activity. Urea-induced unfolding of HNL1 indicates that it unfolds less completely than HbHNL (m-value = 0.63 for HNL1 vs 0.93 kcal/ mol·M for HbHNL), which may account for the ability of HNL1 to better resist irreversible inactivation upon heating. The structure of HNL1 shows changes in hydrogen bond networks that may stabilize regions of the folded structure.

Modeling Structural Constraints on Protein Evolution via Side-Chain Conformational States

Few models of sequence evolution incorporate parameters describing protein structure, despite its high conservation, essential functional role and increasing availability. We present a structurally aware empirical substitution model for amino acid sequence evolution in which proteins are expressed using an expanded alphabet that relays both amino acid identity and structural information. Each character specifies an amino acid as well as information about the rotamer configuration of its side-chain: the discrete geometric pattern of permitted side-chain atomic positions, as defined by the dihedral angles between covalently linked atoms. By assigning rotamer states in 251,194 protein structures and identifying 4,508,390 substitutions between closely related sequences, we generate a 55-state “Dayhoff-like” model that shows that the evolutionary properties of amino acids depend strongly upon side-chain geometry. The model performs as well as or better than traditional 20-state models for divergence time estimation, tree inference, and ancestral state reconstruction. We conclude that not only is rotamer configuration a valuable source of information for phylogenetic studies, but that modeling the concomitant evolution of sequence and structure may have important implications for understanding protein folding and function.

Strategies for Increasing Protein Stability

The stability of wild-type proteins is often a hurdle to their practical use in research, industry, and medicine. The route to engineering stability of a protein of interest lies largely with the available data. Where high-resolution structural data is available, rational design, based on fundamental principles of protein chemistry, can improve protein stability. Recent advances in computational biology and the use of nonnatural amino acids have also provided novel rational methods for improving protein stability. Likewise, the explosion of sequence and structural data available in public databases, in combination with improvements in freely available computational tools, has produced accessible phylogenetic approaches. Trawling modern sequence databases can identify the thermostable homologs of a target protein, and evolutionary data can be quickly generated using available phylogenetic tools. Grafting features from those thermostable homologs or ancestors provides stability improvement through a semi-rational approach. Further, molecular techniques such as directed evolution have shown great promise in delivering designer proteins. These strategies are well documented and newly accessible to the molecular biologist, allowing for rapid enhancements of protein stability.

Ancestral-sequence reconstruction unveils the structural basis of function in mammalian FMOs

Flavin-containing monooxygenases (FMOs) are ubiquitous in all domains of life and metabolize a myriad of xenobiotics, including toxins, pesticides and drugs. However, despite their pharmacological importance, structural information remains bereft. To further our understanding behind their biochemistry and diversity, we used ancestral-sequence reconstruction, kinetic and crystallographic techniques to scrutinize three ancient mammalian FMOs: AncFMO2, AncFMO3-6 and AncFMO5. Remarkably, all AncFMOs could be crystallized and were structurally resolved between 2.7- and 3.2-Å resolution. These crystal structures depict the unprecedented topology of mammalian FMOs. Each employs extensive membrane-binding features and intricate substrate-profiling tunnel networks through a conspicuous membrane-adhering insertion. Furthermore, a glutamate-histidine switch is speculated to induce the distinctive Baeyer-Villiger oxidation activity of FMO5. The AncFMOs exhibited catalysis akin to human FMOs and, with sequence identities between 82% and 92%, represent excellent models. Our study demonstrates the power of ancestral-sequence reconstruction as a strategy for the crystallization of proteins.

Deciphering protein evolution and fitness landscapes with latent space models

Protein sequences contain rich information about protein evolution, fitness landscapes, and stability. Here we investigate how latent space models trained using variational auto-encoders can infer these properties from sequences. Using both simulated and real sequences, we show that the low dimensional latent space representation of sequences, calculated using the encoder model, captures both evolutionary and ancestral relationships between sequences. Together with experimental fitness data and Gaussian process regression, the latent space representation also enables learning the protein fitness landscape in a continuous low dimensional space. Moreover, the model is also useful in predicting protein mutational stability landscapes and quantifying the importance of stability in shaping protein evolution. Overall, we illustrate that the latent space models learned using variational auto-encoders provide a mechanism for exploration of the rich data contained in protein sequences regarding evolution, fitness and stability and hence are well-suited to help guide protein engineering efforts.

Multiple sequence alignments of proteins carry information about evolution, the protein’s fitness landscape and its stability in the face of mutations. Here, the authors demonstrate the utility of latent space models learned using variational autoencoders to infer these properties from sequences.

Highly thermostable carboxylic acid reductases generated by ancestral sequence reconstruction

Carboxylic acid reductases (CARs) are biocatalysts of industrial importance. Their properties, especially their poor stability, render them sub-optimal for use in a bioindustrial pipeline. Here, we employed ancestral sequence reconstruction (ASR) - a burgeoning engineering tool that can identify stabilizing but enzymatically neutral mutations throughout a protein. We used a three-algorithm approach to reconstruct functional ancestors of the Mycobacterial and Nocardial CAR1 orthologues. Ancestral CARs (AncCARs) were confirmed to be CAR enzymes with a preference for aromatic carboxylic acids. Ancestors also showed varied tolerances to solvents, pH and in vivo-like salt concentrations. Compared to well-studied extant CARs, AncCARs had a Tm up to 35 °C higher, with half-lives up to nine times longer than the greatest previously observed. Using ancestral reconstruction we have expanded the existing CAR toolbox with three new thermostable CAR enzymes, providing access to the high temperature biosynthesis of aldehydes to drive new applications in biocatalysis.

A Second Backbone: The Contribution of a Buried Asparagine Ladder to the Global and Local Stability of a Leucine-Rich Repeat Protein

Parallel β-sheet-containing repeat proteins often display a structural motif in which conserved asparagines form a continuous ladder buried within the hydrophobic core. In such "asparagine ladders", the asparagine side-chain amides form a repetitive pattern of hydrogen bonds with neighboring main-chain NH and CO groups. Although asparagine ladders have been thought to be important for stability, there is little experimental evidence to support such speculation. Here we test the contribution of a minimal asparagine ladder from the leucine-rich repeat protein pp32 to stability and investigate lattice rigidity and hydrogen bond character using solution nuclear magnetic resonance (NMR) spectroscopy. Point substitutions of the two ladder asparagines of pp32 are strongly destabilizing and decrease the cooperativity of unfolding. The chemical shifts of the ladder side-chain HZ protons are shifted significantly downfield in the NMR spectrum and have low temperature coefficients, indicative of strong hydrogen bonding. In contrast, the HE protons are shifted upfield and have temperature coefficients close to zero, suggesting an asymmetry in hydrogen bond strength along the ladder. Ladder NH2 groups have weak 1H-15N cross-peak intensities; 1H-15N nuclear Overhauser effect and 15N CPMG experiments show this to be the result of high rigidity. Hydrogen exchange measurements demonstrate that the ladder NH2 groups exchange very slowly, with rates approaching the global exchange limit. Overall, these results show that the asparagine side chains are held in a very rigid, nondynamic structure, making a significant contribution to the overall stability. In this regard, buried asparagine ladders can be considered "second backbones" within the cores of their elongated β-sheet host proteins.

Sequence selection by FitSS4ASR alleviates ancestral sequence reconstruction as exemplified for geranylgeranylglyceryl phosphate synthase

For evolutionary studies, but also for protein engineering, ancestral sequence reconstruction (ASR) has become an indispensable tool. The first step of every ASR protocol is the preparation of a representative sequence set containing at most a few hundred recent homologs whose composition determines decisively the outcome of a reconstruction. A common approach for sequence selection consists of several rounds of manual recompilation that is driven by embedded phylogenetic analyses of the varied sequence sets. For ASR of a geranylgeranylglyceryl phosphate synthase, we additionally utilized FitSS4ASR, which replaces this time-consuming protocol with an efficient and more rational approach. FitSS4ASR applies orthogonal filters to a set of homologs to eliminate outlier sequences and those bearing only a weak phylogenetic signal. To demonstrate the usefulness of FitSS4ASR, we determined experimentally the oligomerization state of eight predecessors, which is a delicate and taxon-specific property. Corresponding ancestors deduced in a manual approach and by means of FitSS4ASR had the same dimeric or hexameric conformation; this concordance testifies to the efficiency of FitSS4ASR for sequence selection. FitSS4ASR-based results of two other ASR experiments were added to the Supporting Information. Program and documentation are available at https://gitlab.bioinf.ur.de/hek61586/FitSS4ASR.

Protein engineers turned evolutionists-the quest for the optimal starting point

The advent of laboratory directed evolution yielded a fruitful crosstalk between the disciplines of molecular evolution and bio-engineering. Here, we outline recent developments in both disciplines with respect to how one can identify the best starting points for directed evolution, such that highly efficient and robust tailor-made enzymes can be obtained with minimal optimization. Directed evolution studies have highlighted essential features of engineer-able enzymes: highly stable, mutationally robust enzymes with the capacity to accept a broad range of substrates. Robust, evolvable enzymes can be inferred from the natural sequence record. Broad substrate spectrum relates to conformational plasticity and can also be predicted by phylogenetic analyses and/or by computational design. Overall, an increasingly powerful toolkit is becoming available for identifying optimal starting points including network analyses of enzyme superfamilies and other bioinformatics methods.

Base editing: precision chemistry on the genome and transcriptome of living cells

RNA-guided programmable nucleases from CRISPR systems generate precise breaks in DNA or RNA at specified positions. In cells, this activity can lead to changes in DNA sequence or RNA transcript abundance. Base editing is a newer genome-editing approach that uses components from CRISPR systems together with other enzymes to directly install point mutations into cellular DNA or RNA without making double-stranded DNA breaks. DNA base editors comprise a catalytically disabled nuclease fused to a nucleobase deaminase enzyme and, in some cases, a DNA glycosylase inhibitor. RNA base editors achieve analogous changes using components that target RNA. Base editors directly convert one base or base pair into another, enabling the efficient installation of point mutations in non-dividing cells without generating excess undesired editing by-products. In this Review, we summarize base-editing strategies to generate specific and precise point mutations in genomic DNA and RNA, highlight recent developments that expand the scope, specificity, precision and in vivo delivery of base editors and discuss limitations and future directions of base editing for research and therapeutic applications.

Consensus designs and thermal stability determinants of a human glutamate transporter

Human excitatory amino acid transporters (EAATs) take up the neurotransmitter glutamate in the brain and are essential to maintain excitatory neurotransmission. Our understanding of the EAATs’ molecular mechanisms has been hampered by the lack of stability of purified protein samples for biophysical analyses. Here, we present approaches based on consensus mutagenesis to obtain thermostable EAAT1 variants that share up to ~95% amino acid identity with the wild type transporters, and remain natively folded and functional. Structural analyses of EAAT1 and the consensus designs using hydrogen-deuterium exchange linked to mass spectrometry show that small and highly cooperative unfolding events at the inter-subunit interface rate-limit their thermal denaturation, while the transport domain unfolds at a later stage in the unfolding pathway. Our findings provide structural insights into the kinetic stability of human glutamate transporters, and introduce general approaches to extend the lifetime of human membrane proteins for biophysical analyses.

Improving cytidine and adenine base editors by expression optimization and ancestral reconstruction

Base editors enable targeted single-nucleotide conversions in genomic DNA. Here we show that expression levels are a bottleneck in base-editing efficiency. We optimize cytidine (BE4) and adenine (ABE7.10) base editors by modification of nuclear localization signals (NLS) and codon usage, and ancestral reconstruction of the deaminase component. The resulting BE4max, AncBE4max, and ABEmax editors correct pathogenic SNPs with substantially increased efficiency in a variety of mammalian cell types.

Tracing a protein's folding pathway over evolutionary time using ancestral sequence reconstruction and hydrogen exchange

The conformations populated during protein folding have been studied for decades; yet, their evolutionary importance remains largely unexplored. Ancestral sequence reconstruction allows access to proteins across evolutionary time, and new methods such as pulsed-labeling hydrogen exchange coupled with mass spectrometry allow determination of folding intermediate structures at near amino-acid resolution. Here, we combine these techniques to monitor the folding of the ribonuclease H family along the evolutionary lineages of T. thermophilus and E. coli RNase H. All homologs and ancestral proteins studied populate a similar folding intermediate despite being separated by billions of years of evolution. Even though this conformation is conserved, the pathway leading to it has diverged over evolutionary time, and rational mutations can alter this trajectory. Our results demonstrate that evolutionary processes can affect the energy landscape to preserve or alter specific features of a protein’s folding pathway.

The burst-phase folding intermediate of ribonuclease H changes conformation over evolutionary history

The amino acid sequence encodes the energy landscape of a protein. Therefore, we expect evolutionary mutations to change features of the protein energy landscape, including the conformations adopted by a polypeptide as it folds to its native state. Ribonucleases H (RNase H) from Escherichia coli and Thermus thermophilus both fold via a partially folded intermediate in which the core region of the protein (helices A-D and strands 4-5) is structured. Strand 1, however, uniquely contributes to the T. thermophilus RNase H folding intermediate (Icore+1 ), but not the E. coli RNase H intermediate (Icore ) (Rosen & Marqusee, PLoS One 2015). We explore the origin of this difference by characterizing the folding intermediate of seven ancestral RNases H spanning the evolutionary history of these two homologs. Using fragment models with or without strand 1 and FRET probes to characterize the folding intermediate of each ancestor, we find a distinct evolutionary trend across the family-the involvement of strand 1 in the folding intermediate is an ancestral feature that is maintained in the thermophilic lineage and is gradually lost in the mesophilic lineage. Evolutionary sequence changes indeed modulate the conformations present on the folding landscape and altered the folding trajectory of RNase H.

Benchmark Analysis of Native and Artificial NAD+-Dependent Enzymes Generated by a Sequence-Based Design Method with or without Phylogenetic Data

The expansion of protein sequence databases has enabled us to design artificial proteins by sequence-based design methods, such as full-consensus design (FCD) and ancestral-sequence reconstruction (ASR). Artificial proteins with enhanced activity levels compared with native ones can potentially be generated by such methods, but successful design is rare because preparing a sequence library by curating the database and selecting a method is difficult. Utilizing a curated library prepared by reducing conservation energies, we successfully designed two artificial l-threonine 3-dehydrogenases (SDR-TDH) with higher activity levels than native SDR-TDH, FcTDH-N1, and AncTDH, using FCD and ASR, respectively. The artificial SDR-TDHs had excellent thermal stability and NAD⁺ recognition compared to native SDR-TDH from Cupriavidus necator (CnTDH); the melting temperatures of FcTDH-N1 and AncTDH were about 10 and 5 °C higher than that of CnTDH, respectively, and the dissociation constants toward NAD⁺ of FcTDH-N1 and AncTDH were 2- and 7-fold lower than that of CnTDH, respectively. Enzymatic efficiency of the artificial SDR-TDHs were comparable to that of CnTDH. Crystal structures of FcTDH-N1 and AncTDH were determined at 2.8 and 2.1 Å resolution, respectively. Structural and MD simulation analysis of the SDR-TDHs indicated that only the flexibility at specific regions was changed, suggesting that multiple mutations introduced in the artificial SDR-TDHs altered their flexibility and thereby affected their enzymatic properties. Benchmark analysis of the SDR-TDHs indicated that both FCD and ASR can generate highly functional proteins if a curated library is prepared appropriately.

Improving cytidine and adenine base editors by expression optimization and ancestral reconstruction

Base editors enable targeted single-nucleotide conversions in genomic DNA. Here we show that expression levels are a bottleneck in base-editing efficiency. We optimize cytidine (BE4) and adenine (ABE7.10) base editors by modification of nuclear localization signals (NLS) and codon usage, and ancestral reconstruction of the deaminase component. The resulting BE4max, AncBE4max, and ABEmax editors correct pathogenic SNPs with substantially increased efficiency in a variety of mammalian cell types.

Evolutionary history of carbon monoxide dehydrogenase/acetyl-CoA synthase, one of the oldest enzymatic complexes

Carbon monoxide dehydrogenase/acetyl-CoA synthase (CODH/ACS) is a five-subunit enzyme complex responsible for the carbonyl branch of the Wood-Ljungdahl (WL) pathway, considered one of the most ancient metabolisms for anaerobic carbon fixation, but its origin and evolutionary history have been unclear. While traditionally associated with methanogens and acetogens, the presence of CODH/ACS homologs has been reported in a large number of uncultured anaerobic lineages. Here, we have carried out an exhaustive phylogenomic study of CODH/ACS in over 6,400 archaeal and bacterial genomes. The identification of complete and likely functional CODH/ACS complexes in these genomes significantly expands its distribution in microbial lineages. The CODH/ACS complex displays astounding conservation and vertical inheritance over geological times. Rare intradomain and interdomain transfer events might tie into important functional transitions, including the acquisition of CODH/ACS in some archaeal methanogens not known to fix carbon, the tinkering of the complex in a clade of model bacterial acetogens, or emergence of archaeal-bacterial hybrid complexes. Once these transfers were clearly identified, our results allowed us to infer the presence of a CODH/ACS complex with at least four subunits in the last universal common ancestor (LUCA). Different scenarios on the possible role of ancestral CODH/ACS are discussed. Despite common assumptions, all are equally compatible with an autotrophic, mixotrophic, or heterotrophic LUCA. Functional characterization of CODH/ACS from a larger spectrum of bacterial and archaeal lineages and detailed evolutionary analysis of the WL methyl branch will help resolve this issue.

Conservation of Specificity in Two Low-Specificity Proteins

Many regulatory proteins bind peptide regions of target proteins and modulate their activity. Such regulatory proteins can often interact with highly diverse target peptides. In many instances, it is not known if the peptide-binding interface discriminates targets in a biological context, or whether biological specificity is achieved exclusively through external factors such as subcellular localization. We used an evolutionary biochemical approach to distinguish these possibilities for two such low-specificity proteins: S100A5 and S100A6. We used isothermal titration calorimetry to study the binding of peptides with diverse sequence and biochemistry to human S100A5 and S100A6. These proteins bound distinct, but overlapping, sets of peptide targets. We then studied the peptide binding properties of orthologs sampled from across five amniote species. Binding specificity was conserved along all lineages, for the last 320 million years, despite the low specificity of each protein. We used ancestral sequence reconstruction to determine the binding specificity of the last common ancestor of the paralogs. The ancestor bound the entire set of peptides bound by modern S100A5 and S100A6 proteins, suggesting that paralog specificity evolved via subfunctionalization. To rule out the possibility that specificity is conserved because it is difficult to modify, we identified a single historical mutation that, when reverted in human S100A5, gave it the ability to bind an S100A6-specific peptide. These results reveal strong evolutionary constraints on peptide binding specificity. Despite being able to bind a large number of targets, the specificity of S100 peptide interfaces is likely important for the biology of these proteins.

Understanding enzyme function evolution from a computational perspective

In this review, we will explore recent computational approaches to understand enzyme evolution from the perspective of protein structure, dynamics and promiscuity. We will present quantitative methods to measure the size of evolutionary steps within a structural domain, allowing the correlation between change in substrate and domain structure to be assessed, and giving insights into the evolvability of different domains in terms of the number, types and sizes of evolutionary steps observed. These approaches will help to understand the evolution of new catalytic and non-catalytic functionality in response to environmental demands, showing potential to guide de novoenzyme design and directed evolution experiments.

Resurrecting the Dead (Molecules)

Biological molecules, like organisms themselves, are subject to genetic drift and may even become "extinct". Molecules that are no longer extant in living systems are of high interest for several reasons including insight into how existing life forms evolved and the possibility that they may have new and useful properties no longer available in currently functioning molecules. Predicting the sequence/structure of such molecules and synthesizing them so that their properties can be tested is the basis of "molecular resurrection" and may lead not only to a deeper understanding of evolution, but also to the production of artificial proteins with novel properties and even to insight into how life itself began.

Reconstructing Ancient Proteins to Understand the Causes of Structure and Function

A central goal in biochemistry is to explain the causes of protein sequence, structure, and function. Mainstream approaches seek to rationalize sequence and structure in terms of their effects on function and to identify function's underlying determinants by comparing related proteins to each other. Although productive, both strategies suffer from intrinsic limitations that have left important aspects of many proteins unexplained. These limits can be overcome by reconstructing ancient proteins, experimentally characterizing their properties, and retracing their evolution through time. This approach has proven to be a powerful means for discovering how historical changes in sequence produced the functions, structures, and other physical/chemical characteristics of modern proteins. It has also illuminated whether protein features evolved because of functional optimization, historical constraint, or blind chance. Here we review recent studies employing ancestral protein reconstruction and show how they have produced new knowledge not only of molecular evolutionary processes but also of the underlying determinants of modern proteins' physical, chemical, and biological properties.

Evidence for the principle of minimal frustration in the evolution of protein folding landscapes

Theoretical and experimental studies have firmly established that protein folding can be described by a funneled energy landscape. This funneled energy landscape is the result of foldable protein sequences evolving following the principle of minimal frustration, which allows proteins to rapidly fold to their native biologically functional conformations. For a protein family with a given functional fold, the principle of minimal frustration suggests that, independent of sequence, all proteins within this family should fold with similar rates. However, depending on the optimal living temperature of the organism, proteins also need to modulate their thermodynamic stability. Consequently, the difference in thermodynamic stability should be primarily caused by differences in the unfolding rates. To test this hypothesis experimentally, we performed comprehensive thermodynamic and kinetic analyses of 15 different proteins from the thioredoxin family. Eight of these thioredoxins were extant proteins from psychrophilic, mesophilic, or thermophilic organisms. The other seven protein sequences were obtained using ancestral sequence reconstruction and can be dated back over 4 billion years. We found that all studied proteins fold with very similar rates but unfold with rates that differ up to three orders of magnitude. The unfolding rates correlate well with the thermodynamic stability of the proteins. Moreover, proteins that unfold slower are more resistant to proteolysis. These results provide direct experimental support to the principle of minimal frustration hypothesis.

Robustness of Reconstructed Ancestral Protein Functions to Statistical Uncertainty

Hypotheses about the functions of ancient proteins and the effects of historical mutations on them are often tested using ancestral protein reconstruction (APR)-phylogenetic inference of ancestral sequences followed by synthesis and experimental characterization. Usually, some sequence sites are ambiguously reconstructed, with two or more statistically plausible states. The extent to which the inferred functions and mutational effects are robust to uncertainty about the ancestral sequence has not been studied systematically. To address this issue, we reconstructed ancestral proteins in three domain families that have different functions, architectures, and degrees of uncertainty; we then experimentally characterized the functional robustness of these proteins when uncertainty was incorporated using several approaches, including sampling amino acid states from the posterior distribution at each site and incorporating the alternative amino acid state at every ambiguous site in the sequence into a single "worst plausible case" protein. In every case, qualitative conclusions about the ancestral proteins' functions and the effects of key historical mutations were robust to sequence uncertainty, with similar functions observed even when scores of alternate amino acids were incorporated. There was some variation in quantitative descriptors of function among plausible sequences, suggesting that experimentally characterizing robustness is particularly important when quantitative estimates of ancient biochemical parameters are desired. The worst plausible case method appears to provide an efficient strategy for characterizing the functional robustness of ancestral proteins to large amounts of sequence uncertainty. Sampling from the posterior distribution sometimes produced artifactually nonfunctional proteins for sequences reconstructed with substantial ambiguity.

Combining ancestral sequence reconstruction with protein design to identify an interface hotspot in a key metabolic enzyme complex

It is important to identify hotspot residues that determine protein-protein interactions in interfaces of macromolecular complexes. We have applied a combination of ancestral sequence reconstruction and protein design to identify hotspots within imidazole glycerol phosphate synthase (ImGPS). ImGPS is a key metabolic enzyme complex, which links histidine and de novo purine biosynthesis and consists of the cyclase subunit HisF and the glutaminase subunit HisH. Initial fluorescence titration experiments showed that HisH from Zymomonas mobilis (zmHisH) binds with high affinity to the reconstructed HisF from the last universal common ancestor (LUCA-HisF) but not to HisF from Pyrobaculum arsenaticum (paHisF), which differ by 103 residues. Subsequent titration experiments with a reconstructed evolutionary intermediate linking LUCA-HisF and paHisF and inspection of the subunit interface of a contemporary ImGPS allowed us to narrow down the differences crucial for zmHisH binding to nine amino acids of HisF. Homology modeling and in silico mutagenesis studies suggested that at most two of these nine HisF residues are crucial for zmHisH binding. These computational results were verified by experimental site-directed mutagenesis, which finally enabled us to pinpoint a single amino acid residue in HisF that is decisive for high-affinity binding of zmHisH. Our work shows that the identification of protein interface hotspots can be very efficient when reconstructed proteins with different binding properties are included in the analysis. Proteins 2017; 85:312-321. © 2016 Wiley Periodicals, Inc.

Evolutionary trend toward kinetic stability in the folding trajectory of RNases H

Proper folding of proteins is critical to producing the biological machinery essential for cellular function. The rates and energetics of a protein's folding process, which is described by its energy landscape, are encoded in the amino acid sequence. Over the course of evolution, this landscape must be maintained such that the protein folds and remains folded over a biologically relevant time scale. How exactly a protein's energy landscape is maintained or altered throughout evolution is unclear. To study how a protein's energy landscape changed over time, we characterized the folding trajectories of ancestral proteins of the ribonuclease H (RNase H) family using ancestral sequence reconstruction to access the evolutionary history between RNases H from mesophilic and thermophilic bacteria. We found that despite large sequence divergence, the overall folding pathway is conserved over billions of years of evolution. There are robust trends in the rates of protein folding and unfolding; both modern RNases H evolved to be more kinetically stable than their most recent common ancestor. Finally, our study demonstrates how a partially folded intermediate provides a readily adaptable folding landscape by allowing the independent tuning of kinetics and thermodynamics.

Using natural sequences and modularity to design common and novel protein topologies

Protein design is still a challenging undertaking, often requiring multiple attempts or iterations for success. Typically, the source of failure is unclear, and scoring metrics appear similar between successful and failed cases. Nevertheless, the use of sequence statistics, modularity and symmetry from natural proteins, combined with computational design both at the coarse-grained and atomistic levels is propelling a new wave of design efforts to success. Here we highlight recent examples of design, showing how the wealth of natural protein sequence and topology data may be leveraged to reduce the search space and increase the likelihood of achieving desired outcomes.