Analyzing protein structure and function using ancestral gene reconstruction

Ancestral Sequence Reconstruction as a Tool for the Elucidation of a Stepwise Evolutionary Adaptation

Ancestral sequence reconstruction (ASR) is a powerful tool to infer primordial sequences from contemporary, i.e., extant ones. An essential element of ASR is the computation of a phylogenetic tree whose leaves are the chosen extant sequences. Most often, the reconstructed sequence related to the root of this tree is of greatest interest: It represents the common ancestor (CA) of the sequences under study. If this sequence encodes a protein, one can "resurrect" the CA by means of gene synthesis technology and study biochemical properties of this extinct predecessor with the help of wet-lab experiments.However, ASR deduces also sequences for all internal nodes of the tree, and the well-considered analysis of these "intermediates" can help to elucidate evolutionary processes. Moreover, one can identify key mutations that alter proteins or protein complexes and are responsible for the differing properties of extant proteins. As an illustrative example, we describe the protocol for the rapid identification of hotspots determining the binding of the two subunits within the heteromeric complex imidazole glycerol phosphate synthase.

Diversification of DNA binding specificities enabled SREBP transcription regulators to expand the repertoire of cellular functions that they govern in fungi

The Sterol Regulatory Element Binding Proteins (SREBPs) are basic-helix-loop-helix transcription regulators that control the expression of sterol biosynthesis genes in higher eukaryotes and some fungi. Surprisingly, SREBPs do not regulate sterol biosynthesis in the ascomycete yeasts (Saccharomycotina) as this role was handed off to an unrelated transcription regulator in this clade. The SREBPs, nonetheless, expanded in fungi such as the ascomycete yeasts Candida spp., raising questions about their role and evolution in these organisms. Here we report that the fungal SREBPs diversified their DNA binding preferences concomitantly with an expansion in function. We establish that several branches of fungal SREBPs preferentially bind non-palindromic DNA sequences, in contrast to the palindromic DNA motifs recognized by most basic-helix-loop-helix proteins (including SREBPs) in higher eukaryotes. Reconstruction and biochemical characterization of the likely ancestor protein suggest that an intrinsic DNA binding promiscuity in the family was resolved by alternative mechanisms in different branches of fungal SREBPs. Furthermore, we show that two SREBPs in the human commensal yeast Candida albicans drive a transcriptional cascade that inhibits a morphological switch under anaerobic conditions. Preventing this morphological transition enhances C. albicans colonization of the mammalian intestine, the fungus' natural niche. Thus, our results illustrate how diversification in DNA binding preferences enabled the functional expansion of a family of eukaryotic transcription regulators.

Mutations Utilize Dynamic Allostery to Confer Resistance in TEM-1 β-lactamase

β-lactamases are enzymes produced by bacteria to hydrolyze β-lactam antibiotics as a common mechanism of resistance. Evolution in such enzymes has been rendering a wide variety of antibiotics impotent, therefore posing a major threat. Clinical and in vitro studies of evolution in TEM-1 β-lactamase have revealed a large number of single point mutations that are responsible for driving resistance to antibiotics and/or inhibitors. The distal locations of these mutations from the active sites suggest that these allosterically modulate the antibiotic resistance. We investigated the effects of resistance driver mutations on the conformational dynamics of the enzyme to provide insights about the mechanism of their long-distance interactions. Through all-atom molecular dynamics (MD) simulations, we obtained the dynamic flexibility profiles of the variants and compared those with that of the wild type TEM-1. While the mutational sites in the variants did not have any direct van der Waals interactions with the active site position S70 and E166, we observed a change in the flexibility of these sites, which play a very critical role in hydrolysis. Such long distance dynamic interactions were further confirmed by dynamic coupling index (DCI) analysis as the sites involved in resistance driving mutations exhibited high dynamic coupling with the active sites. A more exhaustive dynamic analysis, using a selection pressure for ampicillin and cefotaxime resistance on all possible types of substitutions in the amino acid sequence of TEM-1, further demonstrated the observed mechanism. Mutational positions that play a crucial role for the emergence of resistance to new antibiotics exhibited high dynamic coupling with the active site irrespective of their locations. These dynamically coupled positions were neither particularly rigid nor particularly flexible, making them more evolvable positions. Nature utilizes these sites to modulate the dynamics of the catalytic sites instead of mutating the highly rigid positions around the catalytic site.

Ancient mechanisms for the evolution of the bicoid homeodomain's function in fly development

The ancient mechanisms that caused developmental gene regulatory networks to diversify among distantly related taxa are not well understood. Here we use ancestral protein reconstruction, biochemical experiments, and developmental assays of transgenic animals carrying reconstructed ancestral genes to investigate how the transcription factor Bicoid (Bcd) evolved its central role in anterior-posterior patterning in flies. We show that most of Bcd's derived functions are attributable to evolutionary changes within its homeodomain (HD) during a phylogenetic interval >140 million years ago. A single substitution from this period (Q50K) accounts almost entirely for the evolution of Bcd's derived DNA specificity in vitro. In transgenic embryos expressing the reconstructed ancestral HD, however, Q50K confers activation of only a few of Bcd's transcriptional targets and yields a very partial rescue of anterior development. Adding a second historical substitution (M54R) confers regulation of additional Bcd targets and further rescues anterior development. These results indicate that two epistatically interacting mutations played a major role in the evolution of Bcd's controlling regulatory role in early development. They also show how ancestral sequence reconstruction can be combined with in vivo characterization of transgenic animals to illuminate the historical mechanisms of developmental evolution.

The CaspBase: a curated database for evolutionary biochemical studies of caspase functional divergence and ancestral sequence inference

Sequence databases are powerful tools for the contemporary scientists' toolkit. However, most functional annotations in public databases are determined computationally and are not verified by a human expert. While hypotheses generated from computational studies are now amenable to experimentation, the quality of the results relies on the quality of input data. We developed the CaspBase to expedite high-quality dataset compilation of annotated caspase sequences, to maximize phylogenetic signal, and to reduce the noise contributed from public databanks. We describe our methods of curation for the CaspBase and how researchers can acquire sequences from CaspBase.org. Our immediate goal for developing the CaspBase was to optimize the ancestral protein reconstruction (APR) of caspases, and we demonstrate the utility of the CaspBase in APR studies. We also developed the Common Position (CP) system for comparing human caspase family paralogs and suggest the CP system as an update to current reporting methods of caspase amino acid positions. We present a standardized multiple sequence alignment (MSA) for the CP system and show the advantage of using large databases such as the CaspBase in defining structural positions in proteins. Although the results described here pertain to caspase evolution and structure-function studies, the methods can be adapted to any gene family.

Molecular modeling of conformational dynamics and its role in enzyme evolution

With increasing computational power, biomolecular simulations have become an invaluable tool for understanding enzyme mechanisms and the origins of enzyme catalysis. More recently, computational studies have started to focus on understanding how enzyme activity itself evolves, both in terms of enhancing the native or new activities on existing enzyme scaffolds, or completely de novo on previously non-catalytic scaffolds. In this context, both experiment and molecular modeling provided strong evidence for an important role of conformational dynamics in the evolution of enzyme functions. This contribution will present a brief overview of the current state of the art for computationally exploring enzyme conformational dynamics in enzyme evolution, and, using several showcase studies, illustrate the ways molecular modeling can be used to shed light on how enzyme function evolves, at the most fundamental molecular level.

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.

Biotechnological and protein-engineering implications of ancestral protein resurrection

Approximations to the sequences of ancestral proteins can be derived from the sequences of their modern descendants. Proteins encoded by such reconstructed sequences can be prepared in the laboratory and subjected to experimental scrutiny. These 'resurrected' ancestral proteins often display remarkable properties, reflecting ancestral adaptations to intra-cellular and extra-cellular environments that differed from the environments hosting modern/extant proteins. Recent experimental and computational work has specifically discussed high stability, substrate and catalytic promiscuity, conformational flexibility/diversity and altered patterns of interaction with other sub-cellular components. In this review, we discuss these remarkable properties as well as recent attempts to explore their biotechnological and protein-engineering potential.

Molecular basis of hemoglobin adaptation in the high-flying bar-headed goose

During the adaptive evolution of a particular trait, some selectively fixed mutations may be directly causative and others may be purely compensatory. The relative contribution of these two classes of mutation to adaptive phenotypic evolution depends on the form and prevalence of mutational pleiotropy. To investigate the nature of adaptive substitutions and their pleiotropic effects, we used a protein engineering approach to characterize the molecular basis of hemoglobin (Hb) adaptation in the high-flying bar-headed goose (Anser indicus), a hypoxia-tolerant species renowned for its trans-Himalayan migratory flights. To test the effects of observed substitutions on evolutionarily relevant genetic backgrounds, we synthesized all possible genotypic intermediates in the line of descent connecting the wildtype bar-headed goose genotype with the most recent common ancestor of bar-headed goose and its lowland relatives. Site-directed mutagenesis experiments revealed one major-effect mutation that significantly increased Hb-O2 affinity on all possible genetic backgrounds. Two other mutations exhibited smaller average effect sizes and less additivity across backgrounds. One of the latter mutations produced a concomitant increase in the autoxidation rate, a deleterious side-effect that was fully compensated by a second-site mutation at a spatially proximal residue. The experiments revealed three key insights: (i) subtle, localized structural changes can produce large functional effects; (ii) relative effect sizes of function-altering mutations may depend on the sequential order in which they occur; and (iii) compensation of deleterious pleiotropic effects may play an important role in the adaptive evolution of protein function.

Unravelling the Carbohydrate-Binding Preferences of the Carbohydrate-Binding Modules of AMP-Activated Protein Kinase

The β subunit of adenosine monophosphate (AMP)-activated protein kinase (AMPK), which exists as two isoforms (β1 and β2) in humans, has a carbohydrate-binding module (CBM) that interacts with glycogen. Although the β1- and β2-CBMs are structurally similar, with strictly conserved ligand-contact residues, they show different carbohydrate affinities. β2-CBM shows the strongest affinity for both branched and unbranched oligosaccharides and it has recently been shown that a Thr insertion into β2-CBM (Thr101) forms a pocket to accommodate branches. This insertion does not explain why β2-CBM binds all carbohydrates with stronger affinity. Herein, it is shown that residue 134 (Val for β2 and Thr for β1), which does not come into contact with a carbohydrate, appears to account for the affinity difference. Characterisation by NMR spectroscopy, however, suggests that mutant β2-Thr101Δ/Val134Thr differs from that of β1-CBM, and mutant β1-Thr101ins/Thr134Val differs from that of β2-CBM. Furthermore, these mutants are less stable to chemical denaturation, relative to that of wild-type β-CBMs, which confounds the affinity analyses. To support the importance of Thr101 and Val134, the ancestral CBM has been constructed. This CBM retains Thr101 and Val134, which suggests that the extant β1-CBM has a modest loss of function in carbohydrate binding. Because the ancestor bound carbohydrate with equal affinity to that of β2-CBM, it is concluded that residue 134 plays an indirect role in carbohydrate binding.

Structural and Dynamics Comparison of Thermostability in Ancient, Modern, and Consensus Elongation Factor Tus

Rationally engineering thermostability in proteins would create enzymes and receptors that function under harsh industrial applications. Several sequence-based approaches can generate thermostable variants of mesophilic proteins. To gain insight into the mechanisms by which proteins become more stable, we use structural and dynamic analyses to compare two popular approaches, ancestral sequence reconstruction (ASR) and the consensus method, used to generate thermostable variants of Elongation Factor Thermo-unstable (EF-Tu). We present crystal structures of ancestral and consensus EF-Tus, accompanied by molecular dynamics simulations aimed at probing the strategies employed to enhance thermostability. All proteins adopt crystal structures similar to extant EF-Tus, revealing no difference in average structure between the methods. Molecular dynamics reveals that ASR-generated sequences retain dynamic properties similar to extant, thermostable EF-Tu from Thermus aquaticus, while consensus EF-Tu dynamics differ from evolution-based sequences. This work highlights the advantage of ASR for engineering thermostability while preserving natural motions in multidomain proteins.

Evolution of protein specificity: insights from ancestral protein reconstruction

Specific interactions between proteins and their molecular partners drive most biological processes, so understanding how these interactions evolve is an important question for biochemists and evolutionary biologists alike. It is often thought that ancestral proteins were systematically more promiscuous than modern proteins and that specificity usually evolves after gene duplication by partitioning and refining the activities of multifunctional ancestors. However, recent studies using ancestral protein reconstruction (APR) have found that ligand-specific functions in some modern protein families evolved de novo from ancestors that did not already have those functions. Further, the new specific interactions evolved by simple mechanisms, with just a few mutations changing classically recognized biochemical determinants of specificity, such as steric and electrostatic complementarity. Acquiring new specific interactions during evolution therefore appears to be neither difficult nor rare. Rather, it is likely that proteins continually gain and lose new activities over evolutionary time as mutations cause subtle but consequential changes in the shape and electrostatics of interaction interfaces. Only a few of these activities, however, are incorporated into the biological processes that contribute to fitness before they are lost to the ravages of further mutation.

RubisCO Early Oxygenase Activity: A Kinetic and Evolutionary Perspective

RubisCO (D-ribulose 1,5-bisphosphate carboxylase/oxygenase) is Earth's main enzyme responsible for CO₂ fixation via carboxylation of ribulose-1,5-bisphosphate (RuBP) into organic matter. Besides the carboxylation reaction, RubisCO also catalyzes the oxygenation of RuBP by O₂ , which is probably as old as its carboxylation properties. Based on molecular phylogeny, the occurrence of the reactive oxygen species (ROS)-removing system and kinetic properties of different RubisCO forms, we postulated that RubisCO oxygenase activity appeared in local microoxic areas, yet before the appearance of oxygenic photosynthesis. Here, in reviewing the literature, we present a novel hypothesis: the RubisCO early oxygenase activity hypothesis. This hypothesis may be compared with the exaptation hypothesis, according to which latent RubisCO oxygenase properties emerged later during the oxygenation of the Earth's atmosphere. The reconstruction of ancestral RubisCO forms using ancestral sequence reconstruction (ASR) techniques, as a promising way for testing of RubisCO early oxygenase activity hypothesis, is presented.

Back to the future: Rational maps for exploring acetylcholine receptor space and time

Global functions of nicotinic acetylcholine receptors, such as subunit cooperativity and compatibility, likely emerge from a network of amino acid residues distributed across the entire pentameric complex. Identification of such networks has stymied traditional approaches to acetylcholine receptor structure and function, likely due to the cryptic interdependency of their underlying amino acid residues. An emerging evolutionary biochemistry approach, which traces the evolutionary history of acetylcholine receptor subunits, allows for rational mapping of acetylcholine receptor sequence space, and offers new hope for uncovering the amino acid origins of these enigmatic properties.

Ancestral Genome Reconstruction on Whole Genome Level

Comparative genomics, evolutionary biology, and cancer researches require tools to elucidate the evolutionary trajectories and reconstruct the ancestral genomes. Various methods have been developed to infer the genome content and gene ordering of ancestral genomes by using such genomic structural variants. There are mainly two kinds of computational approaches in the ancestral genome reconstruction study. Distance/event-based approaches employ genome evolutionary models and reconstruct the ancestral genomes that minimize the total distance or events over the edges of the given phylogeny. The homology/adjacency-based approaches search for the conserved gene adjacencies and genome structures, and assemble these regions into ancestral genomes along the internal node of the given phylogeny. We review the principles and algorithms of these approaches that can reconstruct the ancestral genomes on the whole genome level. We talk about their advantages and limitations of these approaches in dealing with various genome datasets, evolutionary events, and reconstruction problems. We also talk about the improvements and developments of these approaches in the subsequent researches. We select four most famous and powerful approaches from both distance/event-based and homology/adjacency-based categories to analyze and compare their performances in dealing with different kinds of datasets and evolutionary events. Based on our experiment, GASTS has the best performance in solving the problems with equal genome contents that only have genome rearrangement events. PMAG++ achieves the best performance in solving the problems with unequal genome contents that have all possible complicated evolutionary events.

Reconstructed ancestral enzymes reveal that negative selection drove the evolution of substrate specificity in ADP-dependent kinases

One central goal in molecular evolution is to pinpoint the mechanisms and evolutionary forces that cause an enzyme to change its substrate specificity; however, these processes remain largely unexplored. Using the glycolytic ADP-dependent kinases of archaea, including the orders Thermococcales, Methanosarcinales, and Methanococcales, as a model and employing an approach involving paleoenzymology, evolutionary statistics, and protein structural analysis, we could track changes in substrate specificity during ADP-dependent kinase evolution along with the structural determinants of these changes. To do so, we studied five key resurrected ancestral enzymes as well as their extant counterparts. We found that a major shift in function from a bifunctional ancestor that could phosphorylate either glucose or fructose 6-phosphate (fructose-6-P) as a substrate to a fructose 6-P-specific enzyme was started by a single amino acid substitution resulting in negative selection with a ground-state mode against glucose and a subsequent 1,600-fold change in specificity of the ancestral protein. This change rendered the residual phosphorylation of glucose a promiscuous and physiologically irrelevant activity, highlighting how promiscuity may be an evolutionary vestige of ancestral enzyme activities, which have been eliminated over time. We also could reconstruct the evolutionary history of substrate utilization by using an evolutionary model of discrete binary characters, indicating that substrate uses can be discretely lost or acquired during enzyme evolution. These findings exemplify how negative selection and subtle enzyme changes can lead to major evolutionary shifts in function, which can subsequently generate important adaptive advantages, for example, in improving glycolytic efficiency in Thermococcales.

Exploring the past and the future of protein evolution with ancestral sequence reconstruction: the 'retro' approach to protein engineering

A central goal in molecular evolution is to understand the ways in which genes and proteins evolve in response to changing environments. In the absence of intact DNA from fossils, ancestral sequence reconstruction (ASR) can be used to infer the evolutionary precursors of extant proteins. To date, ancestral proteins belonging to eubacteria, archaea, yeast and vertebrates have been inferred that have been hypothesized to date from between several million to over 3 billion years ago. ASR has yielded insights into the early history of life on Earth and the evolution of proteins and macromolecular complexes. Recently, however, ASR has developed from a tool for testing hypotheses about protein evolution to a useful means for designing novel proteins. The strength of this approach lies in the ability to infer ancestral sequences encoding proteins that have desirable properties compared with contemporary forms, particularly thermostability and broad substrate range, making them good starting points for laboratory evolution. Developments in technologies for DNA sequencing and synthesis and computational phylogenetic analysis have led to an escalation in the number of ancient proteins resurrected in the last decade and greatly facilitated the use of ASR in the burgeoning field of synthetic biology. However, the primary challenge of ASR remains in accurately inferring ancestral states, despite the uncertainty arising from evolutionary models, incomplete sequences and limited phylogenetic trees. This review will focus, firstly, on the use of ASR to uncover links between sequence and phenotype and, secondly, on the practical application of ASR in protein engineering.

Ancestral protein resurrection and engineering opportunities of the mamba aminergic toxins

Mamba venoms contain a multiplicity of three-finger fold aminergic toxins known to interact with various α-adrenergic, muscarinic and dopaminergic receptors with different pharmacological profiles. In order to generate novel functions on this structural scaffold and to avoid the daunting task of producing and screening an overwhelming number of variants generated by a classical protein engineering strategy, we accepted the challenge of resurrecting ancestral proteins, likely to have possessed functional properties. This innovative approach that exploits molecular evolution models to efficiently guide protein engineering, has allowed us to generate a small library of six ancestral toxin (AncTx) variants and associate their pharmacological profiles to key functional substitutions. Among these variants, we identified AncTx1 as the most α_1A-adrenoceptor selective peptide known to date and AncTx5 as the most potent inhibitor of the three α2 adrenoceptor subtypes. Three positions in the ρ-Da1a evolutionary pathway, positions 28, 38 and 43 have been identified as key modulators of the affinities for the α₁ and α_2C adrenoceptor subtypes. Here, we present a first attempt at rational engineering of the aminergic toxins, revealing an epistasis phenomenon.

Convergence of Domain Architecture, Structure, and Ligand Affinity in Animal and Plant RNA-Binding Proteins

Reconstruction of ancestral protein sequences using phylogenetic methods is a powerful technique for directly examining the evolution of molecular function. Although ancestral sequence reconstruction (ASR) is itself very efficient, downstream functional, and structural studies necessary to characterize when and how changes in molecular function occurred are often costly and time-consuming, currently limiting ASR studies to examining a relatively small number of discrete functional shifts. As a result, we have very little direct information about how molecular function evolves across large protein families. Here we develop an approach combining ASR with structure and function prediction to efficiently examine the evolution of ligand affinity across a large family of double-stranded RNA binding proteins (DRBs) spanning animals and plants. We find that the characteristic domain architecture of DRBs-consisting of 2-3 tandem double-stranded RNA binding motifs (dsrms)-arose independently in early animal and plant lineages. The affinity with which individual dsrms bind double-stranded RNA appears to have increased and decreased often across both animal and plant phylogenies, primarily through convergent structural mechanisms involving RNA-contact residues within the β1-β2 loop and a small region of α2. These studies provide some of the first direct information about how protein function evolves across large gene families and suggest that changes in molecular function may occur often and unassociated with major phylogenetic events, such as gene or domain duplications.

The phytochelatin synthase gene in date palm (Phoenix dactylifera L.): Phylogeny, evolution and expression

We studied date palm phytochelatin synthase type I (PdPCS1), which catalyzes the cytosolic synthesis of phytochelatins (PCs), a heavy metal binding protein, in plant cells. The gene encoding PdPCS1 (Pdpcs) consists of 8 exons and 7 introns and encodes a protein of 528 amino acids. PCs gene history was studied using Notung phylogeny. During evolution, gene loss from several lineages was predicted including Proteobacteria, Bilateria and Brassicaceae. In addition, eleven gene duplication events appeared toward interior nodes of the reconciled tree and four gene duplication events appeared toward the external nodes. These latter sequences belong to species with a second copy of PCs suggesting that this gene evolved through subfunctionalization. Pdpcs1 gene expression was measured in seedling hypocotyls exposed to Cd, Cu and Cr using quantitative real-time polymerase chain reaction (qPCR). A Pdpcs1 overexpression was evidenced in P. dactylifera seedlings exposed to metals suggesting that 1-the Pdpcs1 gene is functional, 2-there is an implication of the enzyme in metal detoxification mechanisms. Additionally, the structure of PdPCS1 was predicted using its homologue from Nostoc (cyanobacterium, NsPCS) as a template in Discovery studio and PyMol software. These analyses allowed us to identify the phytochelatin synthase type I enzyme in date palm (PdPCS1) via recognition of key consensus amino acids involved in the catalytic mechanism, and to propose a hypothetical binding and catalytic site for an additional substrate binding cavity.

Diversity as Opportunity: Insights from 600 Million Years of AHR Evolution

The aryl hydrocarbon receptor (AHR) was for many years of interest only to pharmacologists and toxicologists. However, this protein has fundamental roles in biology that are being revealed through studies in diverse animal species. The AHR is an ancient protein. AHR homologs exist in most major groups of modern bilaterian animals, including deuterostomes (chordates, hemichordates, echinoderms) and the two major clades of protostome invertebrates [ecdysozoans (e.g. arthropods and nematodes) and lophotrochozoans (e.g. molluscs and annelids)]. AHR homologs also have been identified in cnidarians such as the sea anemone Nematostella and in the genome of Trichoplax, a placozoan. Bilaterians, cnidarians, and placozoans form the clade Eumetazoa, whose last common ancestor lived approximately 600 million years ago (MYA). The presence of AHR homologs in modern representatives of all these groups indicates that the original eumetazoan animal possessed an AHR homolog. Studies in invertebrates and vertebrates reveal parallel functions of AHR in the development and function of sensory neural systems, suggesting that these may be ancestral roles. Vertebrate animals are characterized by the expansion and diversification of AHRs, via gene and genome duplications, from the ancestral protoAHR into at least five classes of AHR-like proteins: AHR, AHR1, AHR2, AHR3, and AHRR. The evolution of multiple AHRs in vertebrates coincided with the acquisition of high-affinity binding of halogenated and polynuclear aromatic hydrocarbons and the emergence of adaptive functions involving regulation of xenobiotic-metabolizing enzymes and roles in adaptive immunity. The existence of multiple AHRs may have facilitated subfunction partitioning and specialization of specific AHR types in some taxa. Additional research in diverse model and non-model species will continue to enrich our understanding of AHR and its pleiotropic roles in biology and toxicology.

Experimental test and refutation of a classic case of molecular adaptation in Drosophila melanogaster

Identifying the genetic basis for adaptive differences between species requires explicit tests of historical hypotheses concerning the effects of past changes in gene sequence on molecular function, organismal phenotype and fitness. We address this challenge by combining ancestral protein reconstruction with biochemical experiments and physiological analysis of transgenic animals that carry ancestral genes. We tested a widely held hypothesis of molecular adaptation-that changes in the alcohol dehydrogenase protein (ADH) along the lineage leading to Drosophila melanogaster increased the catalytic activity of the enzyme and thereby contributed to the ethanol tolerance and adaptation of the species to its ethanol-rich ecological niche. Our experiments strongly refute the predictions of the adaptive ADH hypothesis and caution against accepting intuitively appealing accounts of historical molecular adaptation that are based on correlative evidence. The experimental strategy we employed can be used to decisively test other adaptive hypotheses and the claims they entail about past biological causality.

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.

Synthetic Protein Switches: Theoretical and Experimental Considerations

Synthetic protein switches with tailored response functions are finding increasing applications as tools in basic research and biotechnology. With a number of successful design strategies emerging, the construction of synthetic protein switches still frequently necessitates an integrated approach that combines detailed biochemical and biophysical characterization in combination with high-throughput screening to construct tailored synthetic protein switches. This is increasingly complemented by computational strategies that aim to reduce the need for costly empirical optimization and thus facilitate the protein design process. Successful computational design approaches range from analyzing phylogenetic data to infer useful structural, biophysical, and biochemical information to modeling the structure and function of proteins ab initio. The following chapter provides an overview over the theoretical considerations and experimental approaches that have been successful applied in the construction of synthetic protein switches.

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.

Experimental macroevolution

The convergence of several disparate research programmes raises the possibility that the long-term evolutionary processes of innovation and radiation may become amenable to laboratory experimentation. Ancestors might be resurrected directly from naturally stored propagules or tissues, or indirectly from the expression of ancestral genes in contemporary genomes. New kinds of organisms might be evolved through artificial selection of major developmental genes. Adaptive radiation can be studied by mimicking major ecological transitions in the laboratory. All of these possibilities are subject to severe quantitative and qualitative limitations. In some cases, however, laboratory experiments may be capable of illuminating the processes responsible for the evolution of new kinds of organisms.

Basic units of protein structure, folding, and function

Study of the hierarchy of domain structure with alternative sets of domains and analysis of discontinuous domains, consisting of remote segments of the polypeptide chain, raised a question about the minimal structural unit of the protein domain. The hypothesis on the decisive role of the polypeptide backbone in determining the elementary units of globular proteins have led to the discovery of closed loops. It is reviewed here how closed loops form the loop-n-lock structure of proteins, providing the foundation for stability and designability of protein folds/domain and underlying their co-translational folding. Simplified protein sequences are considered here with the aim to explore the basic principles that presumably dominated the folding and stability of proteins in the early stages of structural evolution. Elementary functional loops (EFLs), closed loops with one or few catalytic residues, are, in turn, units of the protein function. They are apparent descendants of the prebiotic ring-like peptides, which gave rise to the first functional folds/domains being fused in the beginning of the evolution of protein structure. It is also shown how evolutionary relations between protein functional superfamilies and folds delineated with the help of EFLs can contribute to establishing the rules for design of desired enzymatic functions. Generalized descriptors of the elementary functions are proposed to be used as basic units in the future computational design.

Protein stability: computation, sequence statistics, and new experimental methods

Calculating protein stability and predicting stabilizing mutations remain exceedingly difficult tasks, largely due to the inadequacy of potential functions, the difficulty of modeling entropy and the unfolded state, and challenges of sampling, particularly of backbone conformations. Yet, computational design has produced some remarkably stable proteins in recent years, apparently owing to near ideality in structure and sequence features. With caveats, computational prediction of stability can be used to guide mutation, and mutations derived from consensus sequence analysis, especially improved by recent co-variation filters, are very likely to stabilize without sacrificing function. The combination of computational and statistical approaches with library approaches, including new technologies such as deep sequencing and high throughput stability measurements, point to a very exciting near term future for stability engineering, even with difficult computational issues remaining.

Evolution of a Protein Interaction Domain Family by Tuning Conformational Flexibility

Conformational flexibility allows proteins to adopt multiple functionally important conformations but can also lead to nonfunctional structures. We analyzed the dynamic behavior of the enzyme guanylate kinase as it evolved into the GK protein interaction domain (GK_PID) to investigate the role of flexibility in the evolution of new protein functions. We found that the ancestral enzyme is very flexible, allowing it to adopt open conformations that can bind nucleotide and closed ones that enable catalysis of phosphotransfer from ATP to GMP. Historical mutations that converted the GK from an enzyme to a protein interaction domain dramatically reduce flexibility, predominantly by inhibiting rotations of the protein backbone that are coupled to the global closing motion. Removing flexibility prevents adoption of conformations that cannot fit the protein partner in the binding site. Our results highlight the importance of mutations that optimize protein conformational flexibility with function during evolution.

Emergence of de novo proteins from 'dark genomic matter' by 'grow slow and moult'

Proteins are the workhorses of the cell and, over billions of years, they have evolved an amazing plethora of extremely diverse and versatile structures with equally diverse functions. Evolutionary emergence of new proteins and transitions between existing ones are believed to be rare or even impossible. However, recent advances in comparative genomics have repeatedly called some 10%-30% of all genes without any detectable similarity to existing proteins. Even after careful scrutiny, some of those orphan genes contain protein coding reading frames with detectable transcription and translation. Thus some proteins seem to have emerged from previously non-coding 'dark genomic matter'. These 'de novo' proteins tend to be disordered, fast evolving, weakly expressed but also rapidly assuming novel and physiologically important functions. Here we review mechanisms by which 'de novo' proteins might be created, under which circumstances they may become fixed and why they are elusive. We propose a 'grow slow and moult' model in which first a reading frame is extended, coding for an initially disordered and non-globular appendage which, over time, becomes more structured and may also become associated with other proteins.

Analyzing endocrine system conservation and evolution

Analyzing variation in rates of evolution can provide important insights into the factors that constrain trait evolution, as well as those that promote diversification. Metazoan endocrine systems exhibit apparent variation in evolutionary rates of their constituent components at multiple levels, yet relatively few studies have quantified these patterns and analyzed them in a phylogenetic context. This may be in part due to historical and current data limitations for many endocrine components and taxonomic groups. However, recent technological advancements such as high-throughput sequencing provide the opportunity to collect large-scale comparative data sets for even non-model species. Such ventures will produce a fertile data landscape for evolutionary analyses of nucleic acid and amino acid based endocrine components. Here I summarize evolutionary rate analyses that can be applied to categorical and continuous endocrine traits, and also those for nucleic acid and protein-based components. I emphasize analyses that could be used to test whether other variables (e.g., ecology, ontogenetic timing of expression, etc.) are related to patterns of rate variation and endocrine component diversification. The application of phylogenetic-based rate analyses to comparative endocrine data will greatly enhance our understanding of the factors that have shaped endocrine system evolution.

PhyloBot: A Web Portal for Automated Phylogenetics, Ancestral Sequence Reconstruction, and Exploration of Mutational Trajectories

The method of phylogenetic ancestral sequence reconstruction is a powerful approach for studying evolutionary relationships among protein sequence, structure, and function. In particular, this approach allows investigators to (1) reconstruct and “resurrect” (that is, synthesize in vivo or in vitro) extinct proteins to study how they differ from modern proteins, (2) identify key amino acid changes that, over evolutionary timescales, have altered the function of the protein, and (3) order historical events in the evolution of protein function. Widespread use of this approach has been slow among molecular biologists, in part because the methods require significant computational expertise. Here we present PhyloBot, a web-based software tool that makes ancestral sequence reconstruction easy. Designed for non-experts, it integrates all the necessary software into a single user interface. Additionally, PhyloBot provides interactive tools to explore evolutionary trajectories between ancestors, enabling the rapid generation of hypotheses that can be tested using genetic or biochemical approaches. Early versions of this software were used in previous studies to discover genetic mechanisms underlying the functions of diverse protein families, including V-ATPase ion pumps, DNA-binding transcription regulators, and serine/threonine protein kinases. PhyloBot runs in a web browser, and is available at the following URL: http://www.phylobot.com. The software is implemented in Python using the Django web framework, and runs on elastic cloud computing resources from Amazon Web Services. Users can create and submit jobs on our free server (at the URL listed above), or use our open-source code to launch their own PhyloBot server.

Epistasis in protein evolution

The structure, function, and evolution of proteins depend on physical and genetic interactions among amino acids. Recent studies have used new strategies to explore the prevalence, biochemical mechanisms, and evolutionary implications of these interactions-called epistasis-within proteins. Here we describe an emerging picture of pervasive epistasis in which the physical and biological effects of mutations change over the course of evolution in a lineage-specific fashion. Epistasis can restrict the trajectories available to an evolving protein or open new paths to sequences and functions that would otherwise have been inaccessible. We describe two broad classes of epistatic interactions, which arise from different physical mechanisms and have different effects on evolutionary processes. Specific epistasis-in which one mutation influences the phenotypic effect of few other mutations-is caused by direct and indirect physical interactions between mutations, which nonadditively change the protein's physical properties, such as conformation, stability, or affinity for ligands. In contrast, nonspecific epistasis describes mutations that modify the effect of many others; these typically behave additively with respect to the physical properties of a protein but exhibit epistasis because of a nonlinear relationship between the physical properties and their biological effects, such as function or fitness. Both types of interaction are rampant, but specific epistasis has stronger effects on the rate and outcomes of evolution, because it imposes stricter constraints and modulates evolutionary potential more dramatically; it therefore makes evolution more contingent on low-probability historical events and leaves stronger marks on the sequences, structures, and functions of protein families.

The thermostability and specificity of ancient proteins

Were ancient proteins systematically different than modern proteins? The answer to this question is profoundly important, shaping how we understand the origins of protein biochemical, biophysical, and functional properties. Ancestral sequence reconstruction (ASR), a phylogenetic approach to infer the sequences of ancestral proteins, may reveal such trends. We discuss two proposed trends: a transition from higher to lower thermostability and a tendency for proteins to acquire higher specificity over time. We review the evidence for elevated ancestral thermostability and discuss its possible origins in a changing environmental temperature and/or reconstruction bias. We also conclude that there is, as yet, insufficient data to support a trend from promiscuity to specificity. Finally, we propose future work to understand these proposed evolutionary trends.

Functional Divergence of the Nuclear Receptor NR2C1 as a Modulator of Pluripotentiality During Hominid Evolution

Genes encoding nuclear receptors (NRs) are attractive as candidates for investigating the evolution of gene regulation because they (1) have a direct effect on gene expression and (2) modulate many cellular processes that underlie development. We employed a three-phase investigation linking NR molecular evolution among primates with direct experimental assessment of NR function. Phase 1 was an analysis of NR domain evolution and the results were used to guide the design of phase 2, a codon-model-based survey for alterations of natural selection within the hominids. By using a series of reliability and robustness analyses we selected a single gene, NR2C1, as the best candidate for experimental assessment. We carried out assays to determine whether changes between the ancestral and extant NR2C1s could have impacted stem cell pluripotency (phase 3). We evaluated human, chimpanzee, and ancestral NR2C1 for transcriptional modulation of Oct4 and Nanog (key regulators of pluripotency and cell lineage commitment), promoter activity for Pepck (a proxy for differentiation in numerous cell types), and average size of embryological stem cell colonies (a proxy for the self-renewal capacity of pluripotent cells). Results supported the signal for alteration of natural selection identified in phase 2. We suggest that adaptive evolution of gene regulation has impacted several aspects of pluripotentiality within primates. Our study illustrates that the combination of targeted evolutionary surveys and experimental analysis is an effective strategy for investigating the evolution of gene regulation with respect to developmental phenotypes.

Evolution of estrogen receptors in ray-finned fish and their comparative responses to estrogenic substances

In vertebrates, estrogens play fundamental roles in regulating reproductive activities through estrogen receptors (ESRs), and disruption of estrogen signaling is now of global concern for both wildlife and human health. To date, ESRs of only a limited number of species have been characterized. We investigated the functional diversity and molecular basis or ligand sensitivity of ESRs among ray-finned fish species (Actinopterygii), the most variable group within vertebrates. We cloned and characterized ESRs from several key species in the evolution of ray-finned fish including bichir (Polypteriformes, ESR1 and ESR2) at the basal lineage of ray-finned fish, and arowana (Osteoglossiformes, ESR1 and ESR2b) and eel (Anguilliformes, ESR1, ESR2a and ESR2b) both belonging to ancient early-branching lineages of teleosts, and suggest that ESR2a and ESR2b emerged through teleost-specific whole genome duplication, but an ESR1 paralogue has been lost in the early lineage of euteleost fish species. All cloned ESR isoforms showed similar responses to endogenous and synthetic steroidal estrogens, but they responded differently to non-steroidal estrogenic endocrine disrupting chemicals (EDCs) (e.g., ESR2a exhibits a weaker reporter activity compared with ESR2b). We show that variation in ligand sensitivity of ESRs can be attributed to phylogeny among species of different taxonomic groups in ray-finned fish. The molecular information provided contributes both to understanding of the comparative role of ESRs in the reproductive biology of fish and their comparative responses to EDCs.

A simple method for studying the molecular mechanisms of ultraviolet and violet reception in vertebrates

BACKGROUND

Many vertebrate species use ultraviolet (UV) reception for such basic behaviors as foraging and mating, but many others switched to violet reception and improved their visual resolution. The respective phenotypes are regulated by the short wavelength-sensitive (SWS1) pigments that absorb light maximally (λmax) at ~360 and 395-440 nm. Because of strong epistatic interactions, the biological significance of the extensive mutagenesis results on the molecular basis of spectral tuning in SWS1 pigments and the mechanisms of their phenotypic adaptations remains uncertain.

RESULTS

The magnitudes of the λmax-shifts caused by mutations in a present-day SWS1 pigment and by the corresponding forward mutations in its ancestral pigment are often dramatically different. To resolve these mutagenesis results, the A/B ratio, in which A and B are the areas formed by amino acids at sites 90, 113 and 118 and by those at sites 86, 90 and 118 and 295, respectively, becomes indispensable. Then, all critical mutations that generated the λmax of a SWS1 pigment can be identified by establishing that 1) the difference between the λmax of the ancestral pigment with these mutations and that of the present-day pigment is small (3 ~ 5 nm, depending on the entire λmax-shift) and 2) the difference between the corresponding A/B ratios is < 0.002.

CONCLUSION

Molecular adaptation has been studied mostly by using comparative sequence analyses. These statistical results provide biological hypotheses and need to be tested using experimental means. This is an opportune time to explore the currently available and new genetic systems and test these statistical hypotheses. Evaluating the λmaxs and A/B ratios of mutagenized present-day and their ancestral pigments, we now have a method to identify all critical mutations that are responsible for phenotypic adaptation of SWS1 pigments. The result also explains spectral tuning of the same pigments, a central unanswered question in phototransduction.

Rapid bursts and slow declines: on the possible evolutionary trajectories of enzymes

The evolution of enzymes is often viewed as following a smooth and steady trajectory, from barely functional primordial catalysts to the highly active and specific enzymes that we observe today. In this review, we summarize experimental data that suggest a different reality. Modern examples, such as the emergence of enzymes that hydrolyse human-made pesticides, demonstrate that evolution can be extraordinarily rapid. Experiments to infer and resurrect ancient sequences suggest that some of the first organisms present on the Earth are likely to have possessed highly active enzymes. Reconciling these observations, we argue that rapid bursts of strong selection for increased catalytic efficiency are interspersed with much longer periods in which the catalytic power of an enzyme erodes, through neutral drift and selection for other properties such as cellular energy efficiency or regulation. Thus, many enzymes may have already passed their catalytic peaks.

Plant synthetic biology for molecular engineering of signalling and development

Molecular genetic studies of model plants in the past few decades have identified many key genes and pathways controlling development, metabolism and environmental responses. Recent technological and informatics advances have led to unprecedented volumes of data that may uncover underlying principles of plants as biological systems. The newly emerged discipline of synthetic biology and related molecular engineering approaches is built on this strong foundation. Today, plant regulatory pathways can be reconstituted in heterologous organisms to identify and manipulate parameters influencing signalling outputs. Moreover, regulatory circuits that include receptors, ligands, signal transduction components, epigenetic machinery and molecular motors can be engineered and introduced into plants to create novel traits in a predictive manner. Here, we provide a brief history of plant synthetic biology and significant recent examples of this approach, focusing on how knowledge generated by the reference plant Arabidopsis thaliana has contributed to the rapid rise of this new discipline, and discuss potential future directions.

Protein sequence design and its applications

Design of proteins has far-reaching potentials in diverse areas that span repurposing of the protein scaffold for reactions and substrates that they were not naturally meant for, to catching a glimpse of the ephemeral proteins that nature might have sampled during evolution. These non-natural proteins, either in synthesized or virtual form have opened the scope for the design of entities that not only rival their natural counterparts but also offer a chance to visualize the protein space continuum that might help to relate proteins and understand their associations. Here, we review the recent advances in protein engineering and design, in multiple areas, with a view to drawing attention to their future potential.

Evolution of an ancient protein function involved in organized multicellularity in animals

To form and maintain organized tissues, multicellular organisms orient their mitotic spindles relative to neighboring cells. A molecular complex scaffolded by the GK protein-interaction domain (GKPID) mediates spindle orientation in diverse animal taxa by linking microtubule motor proteins to a marker protein on the cell cortex localized by external cues. Here we illuminate how this complex evolved and commandeered control of spindle orientation from a more ancient mechanism. The complex was assembled through a series of molecular exploitation events, one of which - the evolution of GKPID's capacity to bind the cortical marker protein - can be recapitulated by reintroducing a single historical substitution into the reconstructed ancestral GKPID. This change revealed and repurposed an ancient molecular surface that previously had a radically different function. We show how the physical simplicity of this binding interface enabled the evolution of a new protein function now essential to the biological complexity of many animals.

Antiplasmodial Activity Is an Ancient and Conserved Feature of Tick Defensins

Ancestral sequence reconstruction has been widely used to test evolution-based hypotheses. The genome of the European tick vector, Ixodes ricinus, encodes for defensin peptides with diverse antimicrobial activities against distantly related pathogens. These pathogens include fungi, Gram-negative, and Gram-positive bacteria, i.e., a wide antimicrobial spectrum. Ticks do not transmit these pathogens, suggesting that these defensins may act against a wide range of microbes encountered by ticks during blood feeding or off-host periods. As demonstrated here, these I. ricinus defensins are also effective against the apicomplexan parasite Plasmodium falciparum. To study the general evolution of antimicrobial activity in tick defensins, the ancestral amino acid sequence of chelicerate defensins, which existed approximately 444 million years ago, was reconstructed using publicly available scorpion and tick defensin sequences (named Scorpions-Ticks Defensins Ancestor, STiDA). The activity of STiDA was tested against P. falciparum and the same Gram-negative and Gram-positive bacteria that were used for the I. ricinus defensins. While some extant tick defensins exhibit a wide antimicrobial spectrum, the ancestral defensin showed moderate activity against one of the tested microbes, P. falciparum. This study suggests that amino acid variability and defensin family expansion increased the antimicrobial spectrum of ancestral tick defensins.

Subdomain interactions foster the design of two protein pairs with ∼80% sequence identity but different folds

Metamorphic proteins, including proteins with high levels of sequence identity but different folds, are exceptions to the long-standing rule-of-thumb that proteins with as little as 30% sequence identity adopt the same fold. Which topologies can be bridged by these highly identical sequences remains an open question. Here we bridge two 3-α-helix bundle proteins with two radically different folds. Using a straightforward approach, we engineered the sequences of one subdomain within maltose binding protein (MBP, α/β/α-sandwich) and another within outer surface protein A (OspA, β-sheet) to have high sequence identity (80 and 77%, respectively) with engineered variants of protein G (GA, 3-α-helix bundle). Circular dichroism and nuclear magnetic resonance spectra of all engineered variants demonstrate that they maintain their native conformations despite substantial sequence modification. Furthermore, the MBP variant (80% identical to GA) remained active. Thermodynamic analysis of numerous GA and MBP variants suggests that the key to our approach involved stabilizing the modified MBP and OspA subdomains via external interactions with neighboring substructures, indicating that subdomain interactions can stabilize alternative folds over a broad range of sequence variation. These findings suggest that it is possible to bridge one fold with many other topologies, which has implications for protein folding, evolution, and misfolding diseases.