Protection of Catalytic Cofactors by Polypeptides as a Driver for the Emergence of Primordial Enzymes.

Enzymes catalyze the chemical reactions of life. For nearly half of known enzymes, catalysis requires the binding of small molecules known as cofactors. Polypeptide-cofactor complexes likely formed at a primordial stage and became starting points for the evolution of many efficient enzymes. Yet, evolution has no foresight so the driver for the primordial complex formation is unknown. Here, we use a resurrected ancestral TIM-barrel protein to identify one potential driver. Heme binding at a flexible region of the ancestral structure yields a peroxidation catalyst with enhanced efficiency when compared to free heme. This enhancement, however, does not arise from protein-mediated promotion of catalysis. Rather, it reflects the protection of bound heme from common degradation processes and a resulting longer lifetime and higher effective concentration for the catalyst. Protection of catalytic cofactors by polypeptides emerges as a general mechanism to enhance catalysis and may have plausibly benefited primordial polypeptide-cofactor associations.

Evolution of bacterial steroid biosynthesis and its impact on eukaryogenesis.

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.

Phylogenetic Articulation of Uric Acid Evolution in Mammals and How It Informs a Therapeutic Uricase.

The role of uric acid during primate evolution has remained elusive ever since it was discovered over 100 years ago that humans have unusually high levels of the small molecule in our serum. It has been difficult to generate a neutral or adaptive explanation in part because the uricase enzyme evolved to become a pseudogene in apes thus masking typical signals of sequence evolution. Adding to the difficulty is a lack of clarity on the functional role of uric acid in apes. One popular hypothesis proposes that uric acid is a potent antioxidant that increased in concentration to compensate for the lack of vitamin C synthesis in primate species ∼65 Ma. Here, we have expanded on our previous work with resurrected ancient uricase proteins to better resolve the reshaping of uricase enzymatic activity prior to ape evolution. Our results suggest that the pivotal death-knell to uricase activity occurred between 20 and 30 Ma despite small sequential modifications to its catalytic efficiency for the tens of millions of years since primates lost their ability to synthesize vitamin C, and thus the two appear uncorrelated. We also use this opportunity to demonstrate how molecular evolution can contribute to biomedicine by presenting ancient uricases to human immune cells that assay for innate reactivity against foreign antigens. A highly stable and highly catalytic ancient uricase is shown to elicit a lower immune response in more human haplotypes than other uricases currently in therapeutic development.

Oral Delivery of Nanoparticles Carrying Ancestral Uricase Enzyme Protects against Hyperuricemia in Knockout Mice.

The therapeutic value of delivering recombinant uricase to human patients has been appreciated for decades. The development of therapeutic uricases has been hampered by the fact that humans do not encode an endogenous uricase and therefore most recombinant forms of the protein are recognized as foreign by the immune system and are therefore highly immunogenic. In order to both shield and stabilize the active enzyme, we encapsulated a functional ancestral uricase in recombinant, noninfectious Qß capsid nanoparticles and characterized its catalytic activity. Oral delivery of the nanoparticles moderated key symptoms of kidney dysfunction in uricase-knockout mice by lowering uric acid levels. Histological kidney samples of the treated mice suggest that delivery of recombinant uricase had a protective effect against the destructive effects of uric acid that lead to renal failure caused by hyperuricemia.

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.

Stabilization of Near-Infrared Fluorescent Proteins by Packaging in Virus-like Particles.

Near-IR fluorescent Qß virus-like particles (VLPs) were produced in a high yield by packaging highly red-shifted monomeric and dimeric versions of biliverdin-dependent fluorescent proteins within the capsid shell. The simple addition of biliverdin hydrochloride to the medium during or after Escherichia coli protein expression was enough to produce fully matured encapsidated fluorophores. The packaged near-IR proteins exhibited identical photochemical properties to their nonencapsidated analogues but were far more stable toward heat, chaotrope-induced denaturation, and proteolysis. Noninvasive in vivo imaging showed the VLPs to traffic primarily to the liver after systemic injection in mice, revealing that the particles were easily detected by a standard instrument.

Non-conservation of folding rates in the thioredoxin family reveals degradation of ancestral unassisted-folding.

Evolution involves not only adaptation, but also the degradation of superfluous features. Many examples of degradation at the morphological level are known (vestigial organs, for instance). However, the impact of degradation on molecular evolution has been rarely addressed. Thioredoxins serve as general oxidoreductases in all cells. Here, we report extensive mutational analyses on the folding of modern and resurrected ancestral bacterial thioredoxins. Contrary to claims from recent literature, in vitro folding rates in the thioredoxin family are not evolutionarily conserved, but span at least a ∼100-fold range. Furthermore, modern thioredoxin folding is often substantially slower than ancestral thioredoxin folding. Unassisted folding, as probed in vitro, thus emerges as an ancestral vestigial feature that underwent degradation, plausibly upon the evolutionary emergence of efficient cellular folding assistance. More generally, our results provide evidence that degradation of ancestral features shapes, not only morphological evolution, but also the evolution of individual proteins.

On the Origin of Isoprenoid Biosynthesis.

Isoprenoids and their derivatives represent the largest group of organic compounds in nature and are distributed universally in the three domains of life. Isoprenoids are biosynthesized from isoprenyl diphosphate units, generated by two distinctive biosynthetic pathways: mevalonate pathway and methylerthritol 4-phosphate pathway. Archaea and eukaryotes exclusively have the former pathway, while most bacteria have the latter. Some bacteria, however, are known to possess the mevalonate pathway genes. Understanding the evolutionary history of these two isoprenoid biosynthesis pathways in each domain of life is critical since isoprenoids are so interweaved in the architecture of life that they would have had indispensable roles in the early evolution of life. Our study provides a detailed phylogenetic analysis of enzymes involved in the mevalonate pathway and sheds new light on its evolutionary history. The results suggest that a potential mevalonate pathway is present in the recently discovered superphylum Candidate Phyla Radiation (CPR), and further suggest a strong evolutionary relationship exists between archaea and CPR. Interestingly, CPR harbors the characteristics of both the bacterial-type and archaeal-type mevalonate pathways and may retain signatures regarding the ancestral isoprenoid biosynthesis pathway in the last universal common ancestor. Our study supports the ancient origin of the mevalonate pathway in the three domains of life as previously inferred, but concludes that the evolution of the mevalonate pathway was more complex.

Coevolution of URAT1 and Uricase during Primate Evolution: Implications for Serum Urate Homeostasis and Gout.

Uric acid is the highly insoluble end-product of purine metabolism in humans. Serum levels exceeding the solubility threshold can trigger formation of urate crystals resulting in gouty arthritis. Uric acid is primarily excreted through the kidneys with 90% reabsorbed back into the bloodstream through the uric acid transporter URAT1. This reabsorption process is essential for the high serum uric acid levels found in humans. We discovered that URAT1 proteins from humans and baboons have higher affinity for uric acid compared with transporters from rats and mice. This difference in transport kinetics of URAT1 orthologs, along with inability of modern apes to oxidize uric acid due to loss of the uricase enzyme, prompted us to ask whether these events occurred concomitantly during primate evolution. Ancestral URAT1 sequences were computationally inferred and ancient transporters were resurrected and assayed, revealing that affinity for uric acid was increased during the evolution of primates. This molecular fine-tuning occurred between the origins of simians and their diversification into New- and Old-World monkey and ape lineages. Remarkably, it was driven in large-part by only a few amino acid replacements within the transporter. This alteration in primate URAT1 coincided with changes in uricase that greatly diminished the enzymatic activity and took place 27-77 Ma. These results suggest that the modifications to URAT1 transporters were potentially adaptive and that maintaining more constant, high levels of serum uric acid may have provided an advantage to our primate ancestors.

Experimental Evolution of Escherichia coli Harboring an Ancient Translation Protein.

The ability to design synthetic genes and engineer biological systems at the genome scale opens new means by which to characterize phenotypic states and the responses of biological systems to perturbations. One emerging method involves inserting artificial genes into bacterial genomes and examining how the genome and its new genes adapt to each other. Here we report the development and implementation of a modified approach to this method, in which phylogenetically inferred genes are inserted into a microbial genome, and laboratory evolution is then used to examine the adaptive potential of the resulting hybrid genome. Specifically, we engineered an approximately 700-million-year-old inferred ancestral variant of tufB, an essential gene encoding elongation factor Tu, and inserted it in a modern Escherichia coli genome in place of the native tufB gene. While the ancient homolog was not lethal to the cell, it did cause a twofold decrease in organismal fitness, mainly due to reduced protein dosage. We subsequently evolved replicate hybrid bacterial populations for 2000 generations in the laboratory and examined the adaptive response via fitness assays, whole genome sequencing, proteomics, and biochemical assays. Hybrid lineages exhibit a general adaptive strategy in which the fitness cost of the ancient gene was ameliorated in part by upregulation of protein production. Our results suggest that an ancient-modern recombinant method may pave the way for the synthesis of organisms that exhibit ancient phenotypes, and that laboratory evolution of these organisms may prove useful in elucidating insights into historical adaptive processes.

Evolutionary history and metabolic insights of ancient mammalian uricases.

Uricase is an enzyme involved in purine catabolism and is found in all three domains of life. Curiously, uricase is not functional in some organisms despite its role in converting highly insoluble uric acid into 5-hydroxyisourate. Of particular interest is the observation that apes, including humans, cannot oxidize uric acid, and it appears that multiple, independent evolutionary events led to the silencing or pseudogenization of the uricase gene in ancestral apes. Various arguments have been made to suggest why natural selection would allow the accumulation of uric acid despite the physiological consequences of crystallized monosodium urate acutely causing liver/kidney damage or chronically causing gout. We have applied evolutionary models to understand the history of primate uricases by resurrecting ancestral mammalian intermediates before the pseudogenization events of this gene family. Resurrected proteins reveal that ancestral uricases have steadily decreased in activity since the last common ancestor of mammals gave rise to descendent primate lineages. We were also able to determine the 3D distribution of amino acid replacements as they accumulated during evolutionary history by crystallizing a mammalian uricase protein. Further, ancient and modern uricases were stably transfected into HepG2 liver cells to test one hypothesis that uricase pseudogenization allowed ancient frugivorous apes to rapidly convert fructose into fat. Finally, pharmacokinetics of an ancient uricase injected in rodents suggest that our integrated approach provides the foundation for an evolutionarily-engineered enzyme capable of treating gout and preventing tumor lysis syndrome in human patients.

Mutational studies on resurrected ancestral proteins reveal conservation of site-specific amino acid preferences throughout evolutionary history.

Local protein interactions ("molecular context" effects) dictate amino acid replacements and can be described in terms of site-specific, energetic preferences for any different amino acid. It has been recently debated whether these preferences remain approximately constant during evolution or whether, due to coevolution of sites, they change strongly. Such research highlights an unresolved and fundamental issue with far-reaching implications for phylogenetic analysis and molecular evolution modeling. Here, we take advantage of the recent availability of phenotypically supported laboratory resurrections of Precambrian thioredoxins and ß-lactamases to experimentally address the change of site-specific amino acid preferences over long geological timescales. Extensive mutational analyses support the notion that evolutionary adjustment to a new amino acid may occur, but to a large extent this is insufficient to erase the primitive preference for amino acid replacements. Generally, site-specific amino acid preferences appear to remain conserved throughout evolutionary history despite local sequence divergence. We show such preference conservation to be readily understandable in molecular terms and we provide crystallographic evidence for an intriguing structural-switch mechanism: Energetic preference for an ancestral amino acid in a modern protein can be linked to reorganization upon mutation to the ancestral local structure around the mutated site. Finally, we point out that site-specific preference conservation naturally leads to one plausible evolutionary explanation for the existence of intragenic global suppressor mutations.

Uric acid-dependent inhibition of AMP kinase induces hepatic glucose production in diabetes and starvation: evolutionary implications of the uricase loss in hominids.

Reduced AMP kinase (AMPK) activity has been shown to play a key deleterious role in increased hepatic gluconeogenesis in diabetes, but the mechanism whereby this occurs remains unclear. In this article, we document that another AMP-dependent enzyme, AMP deaminase (AMPD) is activated in the liver of diabetic mice, which parallels with a significant reduction in AMPK activity and a significant increase in intracellular glucose accumulation in human HepG2 cells. AMPD activation is induced by a reduction in intracellular phosphate levels, which is characteristic of insulin resistance and diabetic states. Increased gluconeogenesis is mediated by reduced TORC2 phosphorylation at Ser171 by AMPK in these cells, as well as by the up-regulation of the rate-limiting enzymes PEPCK and G6Pc. The mechanism whereby AMPD controls AMPK activation depends on the production of a specific AMP downstream metabolite through AMPD, uric acid. In this regard, humans have higher uric acid levels than most mammals due to a mutation in uricase, the enzyme involved in uric acid degradation in most mammals, that developed during a period of famine in Europe 1.5 × 10(7) yr ago. Here, working with resurrected ancestral uricases obtained from early hominids, we show that their expression on HepG2 cells is enough to blunt gluconeogenesis in parallel with an up-regulation of AMPK activity. These studies identify a key role AMPD and uric acid in mediating hepatic gluconeogenesis in the diabetic state, via a mechanism involving AMPK down-regulation and overexpression of PEPCK and G6Pc. The uricase mutation in the Miocene likely provided a survival advantage to help maintain glucose levels under conditions of near starvation, but today likely has a role in the pathogenesis of diabetes.

Secondary structure and domain architecture of the 23S and 5S rRNAs.

We present a de novo re-determination of the secondary (2°) structure and domain architecture of the 23S and 5S rRNAs, using 3D structures, determined by X-ray diffraction, as input. In the traditional 2° structure, the center of the 23S rRNA is an extended single strand, which in 3D is seen to be compact and double helical. Accurately assigning nucleotides to helices compels a revision of the 23S rRNA 2° structure. Unlike the traditional 2° structure, the revised 2° structure of the 23S rRNA shows architectural similarity with the 16S rRNA. The revised 2° structure also reveals a clear relationship with the 3D structure and is generalizable to rRNAs of other species from all three domains of life. The 2° structure revision required us to reconsider the domain architecture. We partitioned the 23S rRNA into domains through analysis of molecular interactions, calculations of 2D folding propensities and compactness. The best domain model for the 23S rRNA contains seven domains, not six as previously ascribed. Domain 0 forms the core of the 23S rRNA, to which the other six domains are rooted. Editable 2° structures mapped with various data are provided (http://apollo.chemistry.gatech.edu/RibosomeGallery).

Phenotypic comparisons of consensus variants versus laboratory resurrections of Precambrian proteins.

Consensus-sequence engineering has generated protein variants with enhanced stability, and sometimes, with modulated biological function. Consensus mutations are often interpreted as the introduction of ancestral amino acid residues. However, the precise relationship between consensus engineering and ancestral protein resurrection is not fully understood. Here, we report the properties of proteins encoded by consensus sequences derived from a multiple sequence alignment of extant, class A ß-lactamases, as compared with the properties of ancient Precambrian ß-lactamases resurrected in the laboratory. These comparisons considered primary sequence, secondary, and tertiary structure, as well as stability and catalysis against different antibiotics. Out of the three consensus variants generated, one could not be expressed and purified (likely due to misfolding and/or low stability) and only one displayed substantial stability having substrate promiscuity, although to a lower extent than ancient ß-lactamases. These results: (i) highlight the phenotypic differences between consensus variants and laboratory resurrections of ancestral proteins; (ii) question interpretations of consensus proteins as phenotypic proxies of ancestral proteins; and (iii) support the notion that ancient proteins provide a robust approach toward the preparation of protein variants having large numbers of mutational changes while possessing unique biomolecular properties.

Palaeotemperature trend for Precambrian life inferred from resurrected proteins.

Biosignatures and structures in the geological record indicate that microbial life has inhabited Earth for the past 3.5 billion years or so. Research in the physical sciences has been able to generate statements about the ancient environment that hosted this life. These include the chemical compositions and temperatures of the early ocean and atmosphere. Only recently have the natural sciences been able to provide experimental results describing the environments of ancient life. Our previous work with resurrected proteins indicated that ancient life lived in a hot environment. Here we expand the timescale of resurrected proteins to provide a palaeotemperature trend of the environments that hosted life from 3.5 to 0.5 billion years ago. The thermostability of more than 25 phylogenetically dispersed ancestral elongation factors suggest that the environment supporting ancient life cooled progressively by 30 degrees C during that period. Here we show that our results are robust to potential statistical bias associated with the posterior distribution of inferred character states, phylogenetic ambiguity, and uncertainties in the amino-acid equilibrium frequencies used by evolutionary models. Our results are further supported by a nearly identical cooling trend for the ancient ocean as inferred from the deposition of oxygen isotopes. The convergence of results from natural and physical sciences suggest that ancient life has continually adapted to changes in environmental temperatures throughout its evolutionary history.

Experimental evolution of protein-protein interaction networks.

The modern synthesis of evolutionary theory and genetics has enabled us to discover underlying molecular mechanisms of organismal evolution. We know that in order to maximize an organism's fitness in a particular environment, individual interactions among components of protein and nucleic acid networks need to be optimized by natural selection, or sometimes through random processes, as the organism responds to changes and/or challenges in the environment. Despite the significant role of molecular networks in determining an organism's adaptation to its environment, we still do not know how such inter- and intra-molecular interactions within networks change over time and contribute to an organism's evolvability while maintaining overall network functions. One way to address this challenge is to identify connections between molecular networks and their host organisms, to manipulate these connections, and then attempt to understand how such perturbations influence molecular dynamics of the network and thus influence evolutionary paths and organismal fitness. In the present review, we discuss how integrating evolutionary history with experimental systems that combine tools drawn from molecular evolution, synthetic biology and biochemistry allow us to identify the underlying mechanisms of organismal evolution, particularly from the perspective of protein interaction networks.

Reconstitution of translation from Thermus thermophilus reveals a minimal set of components sufficient for protein synthesis at high temperatures and functional conservation of modern and ancient translation components.

Thermus thermophilus is a thermophilic model organism distantly related to the mesophilic model organism E. coli. We reconstituted protein translation of Thermus thermophilus in vitro from purified ribosomes, transfer ribonucleic acids (tRNAs) and 33 recombinant proteins. This reconstituted system was fully functional, capable of translating natural messenger RNA (mRNA) into active full-length proteins at temperatures up to 65°C and with yields up to 60 µg/ml. Surprisingly, the synthesis of active proteins also occurred at 37°C, a temperature well below the minimal growth temperature for T. thermophilus. A polyamine was required, with tetraamine being most effective, for translation at both high and low temperatures. Using such a defined in vitro system, we demonstrated a minimal set of components that are sufficient for synthesizing active proteins at high temperatures, the functional compatibility of key translation components between T. thermophilus and E. coli, and the functional conservation of a number of resurrected ancient elongation factors. This work sets the stage for future experiments that apply abundant structural information to biochemical characterization of protein translation and folding in T. thermophilus. Because it contains significantly reduced nucleases and proteases, this reconstituted thermostable cell-free protein synthesis system may enable in vitro engineering of proteins with improved thermostability.

Resurrecting ancestral alcohol dehydrogenases from yeast.

Modern yeast living in fleshy fruits rapidly convert sugars into bulk ethanol through pyruvate. Pyruvate loses carbon dioxide to produce acetaldehyde, which is reduced by alcohol dehydrogenase 1 (Adh1) to ethanol, which accumulates. Yeast later consumes the accumulated ethanol, exploiting Adh2, an Adh1 homolog differing by 24 (of 348) amino acids. As many microorganisms cannot grow in ethanol, accumulated ethanol may help yeast defend resources in the fruit. We report here the resurrection of the last common ancestor of Adh1 and Adh2, called Adh(A). The kinetic behavior of Adh(A) suggests that the ancestor was optimized to make (not consume) ethanol. This is consistent with the hypothesis that before the Adh1-Adh2 duplication, yeast did not accumulate ethanol for later consumption but rather used Adh(A) to recycle NADH generated in the glycolytic pathway. Silent nucleotide dating suggests that the Adh1-Adh2 duplication occurred near the time of duplication of several other proteins involved in the accumulation of ethanol, possibly in the Cretaceous age when fleshy fruits arose. These results help to connect the chemical behavior of these enzymes through systems analysis to a time of global ecosystem change, a small but useful step towards a planetary systems biology.

Impact of steroid biosynthesis on the aerobic adaptation of eukaryotes.

Steroids are indispensable components of the eukaryotic cellular membrane and the acquisition of steroid biosynthesis was a key factor that enabled the evolution of eukaryotes. The polycyclic carbon structures of steroids can be preserved in sedimentary rocks as chemical fossils for billions of years and thus provide invaluable clues to trace eukaryotic evolution from the distant past. Steroid biosynthesis consists of (1) the production of protosteroids and (2) the subsequent modifications toward "modern-type" steroids such as cholesterol and stigmasterol. While protosteroid biosynthesis requires only two genes for the cyclization of squalene, complete modification of protosteroids involves ~10 additional genes. Eukaryotes universally possess at least some of those additional genes and thus produce modern-type steroids as major final products. The geological biomarker records suggest a prolonged period of solely protosteroid production in the mid-Proterozoic before the advent of modern-type steroids in the Neoproterozoic. It has been proposed that mid-Proterozoic protosteroids were produced by hypothetical stem-group eukaryotes that presumably possessed genes only for protosteroid production, even though in modern environments protosteroid production as a final product is found exclusively in bacteria. The host identity of mid-Proterozoic steroid producers is crucial for understanding the early evolution of eukaryotes. In this perspective, we discuss how geological biomarker data and genetic data complement each other and potentially provide a more coherent scenario for the evolution of steroids and associated early eukaryotes. We further discuss the potential impacts that steroids had on the evolution of aerobic metabolism in eukaryotes, which may have been an important factor for the eventual ecological dominance of eukaryotes in many modern environments.