Evolutionary biochemistry: revealing the historical and physical causes of protein properties

Evolutionary drivers of thermoadaptation in enzyme catalysis

With early life likely to have existed in a hot environment, enzymes had to cope with an inherent drop in catalytic speed caused by lowered temperature. Here we characterize the molecular mechanisms underlying thermoadaptation of enzyme catalysis in adenylate kinase using ancestral sequence reconstruction spanning 3 billion years of evolution. We show that evolution solved the enzyme's key kinetic obstacle-how to maintain catalytic speed on a cooler Earth-by exploiting transition-state heat capacity. Tracing the evolution of enzyme activity and stability from the hot-start toward modern hyperthermophilic, mesophilic, and psychrophilic organisms illustrates active pressure versus passive drift in evolution on a molecular level, refutes the debated activity/stability trade-off, and suggests that the catalytic speed of adenylate kinase is an evolutionary driver for organismal fitness.

Evolutionary interplay between structure, energy and epistasis in the coat protein of the ϕX174 phage family

Viral capsids are structurally constrained by interactions among the amino acids (AAs) of their constituent proteins. Therefore, epistasis is expected to evolve among physically interacting sites and to influence the rates of substitution. To study the evolution of epistasis, we focused on the major structural protein of the ϕX174 phage family by first reconstructing the ancestral protein sequences of 18 species using a Bayesian statistical framework. The inferred ancestral reconstruction differed at eight AAs, for a total of 256 possible ancestral haplotypes. For each ancestral haplotype and the extant species, we estimated, in silico, the distribution of free energies and epistasis of the capsid structure. We found that free energy has not significantly increased but epistasis has. We decomposed epistasis up to fifth order and found that higher-order epistasis sometimes compensates pairwise interactions making the free energy seem additive. The dN/dS ratio is low, suggesting strong purifying selection, and that structure is under stabilizing selection. We synthesized phages carrying ancestral haplotypes of the coat protein gene and measured their fitness experimentally. Our findings indicate that stabilizing mutations can have higher fitness, and that fitness optima do not necessarily coincide with energy minima.

Constrained evolution of a bispecific enzyme: lessons for biocatalyst design

Analysis of the natural evolution of bispecificity in triazine hydrolase highlights the importance of epistasis in protein engineering and evolution.

Improved insights into protein thermal stability: from the molecular to the structurome scale

Despite the intense efforts of the last decades to understand the thermal stability of proteins, the mechanisms responsible for its modulation still remain debated. In this investigation, we tackle this issue by showing how a multiscale perspective can yield new insights. With the help of temperature-dependent statistical potentials, we analysed some amino acid interactions at the molecular level, which are suggested to be relevant for the enhancement of thermal resistance. We then investigated the thermal stability at the protein level by quantifying its modification upon amino acid substitutions. Finally, a large scale analysis of protein stability-at the structurome level-contributed to the clarification of the relation between stability and natural evolution, thereby showing that the mutational profile of proteins differs according to their thermal properties. Some considerations on how the multiscale approach could help in unravelling the protein stability mechanisms are briefly discussed.This article is part of the themed issue 'Multiscale modelling at the physics-chemistry-biology interface'.

Resurrecting ancestral structural dynamics of an antiviral immune receptor: adaptive binding pocket reorganization repeatedly shifts RNA preference

Background

Although resurrecting ancestral proteins is a powerful tool for understanding the molecular-functional evolution of gene families, nearly all studies have examined proteins functioning in relatively stable biological processes. The extent to which more dynamic systems obey the same ‘rules’ governing stable processes is unclear. Here we present the first detailed investigation of the functional evolution of the RIG-like receptors (RLRs), a family of innate immune receptors that detect viral RNA in the cytoplasm.

Results

Using kinetic binding assays and molecular dynamics simulations of ancestral proteins, we demonstrate how a small number of adaptive protein-coding changes repeatedly shifted the RNA preference of RLRs throughout animal evolution by reorganizing the shape and electrostatic distribution across the RNA binding pocket, altering the hydrogen bond network between the RLR and its RNA target. In contrast to observations of proteins involved in metabolism and development, we find that RLR-RNA preference ‘flip flopped’ between two functional states, and shifts in RNA preference were not always coupled to gene duplications or speciation events. We demonstrate at least one reversion of RLR-RNA preference from a derived to an ancestral function through a novel structural mechanism, indicating multiple structural implementations of similar functions.

Conclusions

Our results suggest a model in which frequent shifts in selection pressures imposed by an evolutionary arms race preclude the long-term functional optimization observed in stable biological systems. As a result, the evolutionary dynamics of immune receptors may be less constrained by structural epistasis and historical contingency.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0818-6) contains supplementary material, which is available to authorized users.

Bridging the physical scales in evolutionary biology: from protein sequence space to fitness of organisms and populations

Bridging the gap between the molecular properties of proteins and organismal/population fitness is essential for understanding evolutionary processes. This task requires the integration of the several physical scales of biological organization, each defined by a distinct set of mechanisms and constraints, into a single unifying model. The molecular scale is dominated by the constraints imposed by the physico-chemical properties of proteins and their substrates, which give rise to trade-offs and epistatic (non-additive) effects of mutations. At the systems scale, biological networks modulate protein expression and can either buffer or enhance the fitness effects of mutations. The population scale is influenced by the mutational input, selection regimes, and stochastic changes affecting the size and structure of populations, which eventually determine the evolutionary fate of mutations. Here, we summarize the recent advances in theory, computer simulations, and experiments that advance our understanding of the links between various physical scales in biology.

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.

Drug targets evolve, and so should the methods

Design of specific kinase inhibitors is an appealing approach for developing new anticancer treatments. However, only a few success stories have been reported to date. Here we demonstrate how the combination of old-fashioned and new biophysical tools together with recent advances in genomics and molecular evolution can aid in overcoming existing limitations.

The Deep Origin and Recent Loss of Venom Toxin Genes in Rattlesnakes

The genetic origin of novel traits is a central but challenging puzzle in evolutionary biology. Among snakes, phospholipase A2 (PLA2)-related toxins have evolved in different lineages to function as potent neurotoxins, myotoxins, or hemotoxins. Here, we traced the genomic origin and evolution of PLA2 toxins by examining PLA2 gene number, organization, and expression in both neurotoxic and non-neurotoxic rattlesnakes. We found that even though most North American rattlesnakes do not produce neurotoxins, the genes of a specialized heterodimeric neurotoxin predate the origin of rattlesnakes and were present in their last common ancestor (∼22 mya). The neurotoxin genes were then deleted independently in the lineages leading to the Western Diamondback (Crotalus atrox) and Eastern Diamondback (C. adamanteus) rattlesnakes (∼6 mya), while a PLA2 myotoxin gene retained in C. atrox was deleted from the neurotoxic Mojave rattlesnake (C. scutulatus; ∼4 mya). The rapid evolution of PLA2 gene number appears to be due to transposon invasion that provided a template for non-allelic homologous recombination.

Evolution of Ubiquinone Biosynthesis: Multiple Proteobacterial Enzymes with Various Regioselectivities To Catalyze Three Contiguous Aromatic Hydroxylation Reactions

UQ, a key molecule for cellular bioenergetics that is conserved from proteobacteria to humans, appeared in an ancestral proteobacterium more than 2 billion years ago. UQ biosynthesis has been studied only in a few model organisms, and thus, the diversity of UQ biosynthesis pathways is largely unknown. In the work reported here, we conducted a phylogenomic analysis of hydroxylases involved in UQ biosynthesis. Our results support the existence of at least two UQ hydroxylases in the proteobacterial ancestor, and yet, we show that their number varies from one to four in extant proteobacterial species. Our biochemical experiments demonstrated that bacteria containing only one or two UQ hydroxylases have developed generalist enzymes that are able to catalyze several steps of UQ biosynthesis. Our study documents a rare case where evolution favored the broadening of an enzyme’s regioselectivity, which resulted in gene loss in several proteobacterial species with small genomes.

The ubiquitous ATP synthase uses an electrochemical gradient to synthesize cellular energy in the form of ATP. The production of this electrochemical gradient relies on liposoluble proton carriers like ubiquinone (UQ), which is used in the respiratory chains of eukaryotes and proteobacteria. The biosynthesis of UQ requires three hydroxylation reactions on contiguous positions of an aromatic ring. In Escherichia coli, each of three UQ flavin monooxygenases (FMOs), called UbiF, UbiH, and UbiI, modifies a single position of the aromatic ring. This pattern of three hydroxylation reactions/three proteins has been accepted as a paradigm in UQ biology. Using a phylogenetic analysis, we found that UbiF, UbiH, and UbiI are detected only in a small fraction of proteobacteria, and we identified two new types of UQ FMOs: UbiM, which is distributed in members of the alpha, beta, and gamma classes of proteobacteria, and UbiL, which is restricted to members of the alphaproteobacteria. Remarkably, the ubiL and ubiM genes were found in genomes with fewer than three UQ hydroxylase-encoding genes. We demonstrated, using biochemical approaches, that UbiL from Rhodospirillum rubrum and UbiM from Neisseria meningitidis hydroxylate, respectively, two and three positions of the aromatic ring during UQ biosynthesis. We conclude that bacteria have evolved a large repertoire of hydroxylase combinations for UQ biosynthesis, including pathways with either three specialist enzymes or pathways with one or two generalist enzymes of broader regioselectivity. The emergence of the latter is potentially related to genome reduction events.

IMPORTANCE UQ, a key molecule for cellular bioenergetics that is conserved from proteobacteria to humans, appeared in an ancestral proteobacterium more than 2 billion years ago. UQ biosynthesis has been studied only in a few model organisms, and thus, the diversity of UQ biosynthesis pathways is largely unknown. In the work reported here, we conducted a phylogenomic analysis of hydroxylases involved in UQ biosynthesis. Our results support the existence of at least two UQ hydroxylases in the proteobacterial ancestor, and yet, we show that their number varies from one to four in extant proteobacterial species. Our biochemical experiments demonstrated that bacteria containing only one or two UQ hydroxylases have developed generalist enzymes that are able to catalyze several steps of UQ biosynthesis. Our study documents a rare case where evolution favored the broadening of an enzyme’s regioselectivity, which resulted in gene loss in several proteobacterial species with small genomes.

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.

Conformational Tinkering Drives Evolution of a Promiscuous Activity through Indirect Mutational Effects

How remote mutations can lead to changes in enzyme function at a molecular level is a central question in evolutionary biochemistry and biophysics. Here, we combine laboratory evolution with biochemical, structural, genetic, and computational analysis to dissect the molecular basis for the functional optimization of phosphotriesterase activity in a bacterial lactonase (AiiA) from the metallo-β-lactamase (MBL) superfamily. We show that a 1000-fold increase in phosphotriesterase activity is caused by a more favorable catalytic binding position of the paraoxon substrate in the evolved enzyme that resulted from conformational tinkering of the active site through peripheral mutations. A nonmutated active site residue, Phe68, was displaced by ∼3 Å through the indirect effects of two second-shell trajectory mutations, allowing molecular interactions between the residue and paraoxon. Comparative mutational scanning, i.e., examining the effects of alanine mutagenesis on different genetic backgrounds, revealed significant changes in the functional roles of Phe68 and other nonmutated active site residues caused by the indirect effects of trajectory mutations. Our work provides a quantitative measurement of the impact of second-shell mutations on the catalytic contributions of nonmutated residues and unveils the underlying intramolecular network of strong epistatic mutational relationships between active site residues and more remote residues. Defining these long-range conformational and functional epistatic relationships has allowed us to better understand the subtle, but cumulatively significant, role of second- and third-shell mutations in 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.

Venom Resistance as a Model for Understanding the Molecular Basis of Complex Coevolutionary Adaptations

SynopsisVenom and venom resistance are molecular phenotypes widely considered to have diversified through coevolution between predators and prey. However, while evolutionary and functional studies on venom have been extensive, little is known about the molecular basis, variation, and complexity of venom resistance. We review known mechanisms of venom resistance and relate these mechanisms to their predicted impact on coevolutionary dynamics with venomous enemies. We then describe two conceptual approaches which can be used to examine venom/resistance systems. At the intraspecific level, tests of local adaptation in venom and resistance phenotypes can identify the functional mechanisms governing the outcomes of coevolution. At deeper evolutionary timescales, the combination of phylogenetically informed analyses of protein evolution coupled with studies of protein function promise to elucidate the mode and tempo of evolutionary change on potentially coevolving genes. We highlight case studies that use each approach to extend our knowledge of these systems as well as address larger questions about coevolutionary dynamics. We argue that resistance and venom are phenotypic traits which hold exceptional promise for investigating the mechanisms, dynamics, and outcomes of coevolution at the molecular level. Furthermore, extending the understanding of single gene-for-gene interactions to the whole resistance and venom phenotypes may provide a model system for examining the molecular and evolutionary dynamics of complex multi-gene interactions.

Chopping and Changing: the Evolution of the Flavin-dependent Monooxygenases

Flavin-dependent monooxygenases play a variety of key physiological roles and are also very powerful biotechnological tools. These enzymes have been classified into eight different classes (A–H) based on their sequences and biochemical features. By combining structural and sequence analysis, and phylogenetic inference, we have explored the evolutionary history of classes A, B, E, F, and G and demonstrate that their multidomain architectures reflect their phylogenetic relationships, suggesting that the main evolutionary steps in their divergence are likely to have arisen from the recruitment of different domains. Additionally, the functional divergence within in each class appears to have been the result of other mechanisms such as a complex set of single-point mutations. Our results reinforce the idea that a main constraint on the evolution of cofactor-dependent enzymes is the functional binding of the cofactor. Additionally, a remarkable feature of this family is that the sequence of the key flavin adenine dinucleotide-binding domain is split into at least two parts in all classes studied here. We propose a complex set of evolutionary events that gave rise to the origin of the different classes within this family.

•

Changes in domain architectures reflect the phylogeny of flavin monooxygenases.

•

Recruitment of different domains has been the main force driving its evolution.

•

A notable feature of flavin monooxygenases is that the flavin adenine dinucleotide-binding domain is split.

•

Classes of monooxygenases emerged from an ancestral domain by structural changes.

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.

PPARGC1α gene DNA methylation variations in human placenta mediate the link between maternal hyperglycemia and leptin levels in newborns

Background

Children exposed to gestational diabetes mellitus (GDM) are at a higher risk of developing obesity and type 2 diabetes. This susceptibility might involve brown adipose tissue (BAT), which is suspected to protect against obesity. The objective of this study is to assess whether fetal exposure to maternal hyperglycemia is associated with DNA methylation variations in genes involved in BAT genesis and activation.

Methods

DNA methylation levels at the PRDM16, BMP7, CTBP2, and PPARGC1α gene loci were measured in placenta samples using bisulfite pyrosequencing in E-21 (n = 133; 33 cases of GDM) and the HumanMethylation450 array in Gen3G (n = 172, all from non-diabetic women) birth cohorts. Glucose tolerance was assessed in all women using an oral glucose tolerance test at the second trimester of pregnancy. Participating women were extensively phenotyped throughout pregnancy, and placenta and cord blood samples were collected at birth.

Results

We report that maternal glycemia at the second and third trimester of pregnancy are correlated with variations in DNA methylation levels at PRDM16, BMP7, and PPARGC1α and with cord blood leptin levels. Variations in PRDM16 and PPARGC1α DNA methylation levels were also correlated with cord blood leptin levels. Mediation analyses support that DNA methylation variations at the PPARGC1α gene locus explain 0.8 % of the cord blood leptin levels variance independently of maternal fasting glucose levels (p = 0.05).

Conclusions

These results suggest that maternal glucose in pregnancy could produce variations in DNA methylation in BAT-related genes and that some of these DNA methylation marks seem to mediate the impact of maternal glycemia on cord blood leptin levels, an adipokine regulating body weight.

Electronic supplementary material

The online version of this article (doi:10.1186/s13148-016-0239-9) contains supplementary material, which is available to authorized users.

Theoretical Insights into the Biophysics of Protein Bi-stability and Evolutionary Switches

Deciphering the effects of nonsynonymous mutations on protein structure is central to many areas of biomedical research and is of fundamental importance to the study of molecular evolution. Much of the investigation of protein evolution has focused on mutations that leave a protein’s folded structure essentially unchanged. However, to evolve novel folds of proteins, mutations that lead to large conformational modifications have to be involved. Unraveling the basic biophysics of such mutations is a challenge to theory, especially when only one or two amino acid substitutions cause a large-scale conformational switch. Among the few such mutational switches identified experimentally, the one between the G_A all-α and G_B α+β folds is extensively characterized; but all-atom simulations using fully transferrable potentials have not been able to account for this striking switching behavior. Here we introduce an explicit-chain model that combines structure-based native biases for multiple alternative structures with a general physical atomic force field, and apply this construct to twelve mutants spanning the sequence variation between G_A and G_B. In agreement with experiment, we observe conformational switching from G_A to G_B upon a single L45Y substitution in the GA98 mutant. In line with the latent evolutionary potential concept, our model shows a gradual sequence-dependent change in fold preference in the mutants before this switch. Our analysis also indicates that a sharp G_A/G_B switch may arise from the orientation dependence of aromatic π-interactions. These findings provide physical insights toward rationalizing, predicting and designing evolutionary conformational switches.

The biological functions of globular proteins are intimately related to their folded structures and their associated conformational fluctuations. Evolution of new structures is an important avenue to new functions. Although many mutations do not change the folded state, experiments indicate that a single amino acid substitution can lead to a drastic change in the folded structure. The physics of this switch-like behavior remains to be elucidated. Here we develop a computational model for the relevant physical forces, showing that mutations can lead to new folds by passing through intermediate sequences where the old and new folds occur with varying probabilities. Our approach helps provide a general physical account of conformational switching in evolution and mutational effects on conformational dynamics.

Adaptive evolution of complex innovations through stepwise metabolic niche expansion

A central challenge in evolutionary biology concerns the mechanisms by which complex metabolic innovations requiring multiple mutations arise. Here, we propose that metabolic innovations accessible through the addition of a single reaction serve as stepping stones towards the later establishment of complex metabolic features in another environment. We demonstrate the feasibility of this hypothesis through three complementary analyses. First, using genome-scale metabolic modelling, we show that complex metabolic innovations in Escherichia coli can arise via changing nutrient conditions. Second, using phylogenetic approaches, we demonstrate that the acquisition patterns of complex metabolic pathways during the evolutionary history of bacterial genomes support the hypothesis. Third, we show how adaptation of laboratory populations of E. coli to one carbon source facilitates the later adaptation to another carbon source. Our work demonstrates how complex innovations can evolve through series of adaptive steps without the need to invoke non-adaptive processes.

Amino Acid Change in an Orchid Desaturase Enables Mimicry of the Pollinator's Sex Pheromone

Mimicry illustrates the power of selection to produce phenotypic convergence in biology [1]. A striking example is the imitation of female insects by plants that are pollinated by sexual deception of males of the same insect species [2-4]. This involves mimicry of visual, tactile, and chemical signals of females [2-7], especially their sex pheromones [8-11]. The Mediterranean orchid Ophrys exaltata employs chemical mimicry of cuticular hydrocarbons, particularly the 7-alkenes, in an insect sex pheromone to attract and elicit mating behavior in its pollinators, males of the cellophane bee Colletes cunicularius [11-13]. A difference in alkene double-bond positions is responsible for reproductive isolation between O. exaltata and closely related species, such as O. sphegodes [13-16]. We show that these 7-alkenes are likely determined by the action of the stearoyl-acyl-carrier-protein desaturase (SAD) homolog SAD5. After gene duplication, changes in subcellular localization relative to the ancestral housekeeping desaturase may have allowed proto-SAD5's reaction products to undergo further biosynthesis to both 7- and 9-alkenes. Such ancestral coproduction of two alkene classes may have led to pollinator-mediated deleterious pleiotropy. Despite possible evolutionary intermediates with reduced activity, amino acid changes at the bottom of the substrate-binding cavity have conferred enzyme specificity for 7-alkene biosynthesis by preventing the binding of longer-chained fatty acid (FA) precursors by the enzyme. This change in desaturase function enabled the orchid to perfect its chemical mimicry of pollinator sex pheromones by escape from deleterious pleiotropy, supporting a role of pleiotropy in determining the possible trajectories of adaptive evolution.

Functional evolution of IGF2:IGF2R domain 11 binding generates novel structural interactions and a specific IGF2 antagonist

Among the 15 extracellular domains of the mannose 6-phosphate/insulin-like growth factor-2 receptor (M6P/IGF2R), domain 11 has evolved a binding site for IGF2 to negatively regulate ligand bioavailability and mammalian growth. Despite the highly evolved structural loops of the IGF2:domain 11 binding site, affinity-enhancing AB loop mutations suggest that binding is modifiable. Here we examine the extent to which IGF2:domain 11 affinity, and its specificity over IGF1, can be enhanced, and we examine the structural basis of the mechanistic and functional consequences. Domain 11 binding loop mutants were selected by yeast surface display combined with high-resolution structure-based predictions, and validated by surface plasmon resonance. We discovered previously unidentified mutations in the ligand-interacting surface binding loops (AB, CD, FG, and HI). Five combined mutations increased rigidity of the AB loop, as confirmed by NMR. When added to three independently identified CD and FG loop mutations that reduced the koff value by twofold, these mutations resulted in an overall selective 100-fold improvement in affinity. The structural basis of the evolved affinity was improved shape complementarity established by interloop (AB-CD) and intraloop (FG-FG) side chain interactions. The high affinity of the combinatorial domain 11 Fc fusion proteins functioned as ligand-soluble antagonists or traps that depleted pathological IGF2 isoforms from serum and abrogated IGF2-dependent signaling in vivo. An evolved and reengineered high-specificity M6P/IGF2R domain 11 binding site for IGF2 may improve therapeutic targeting of the frequent IGF2 gain of function observed in human cancer.

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.

Genetically Determined Variation in Lysis Time Variance in the Bacteriophage φX174

Researchers in evolutionary genetics recently have recognized an exciting opportunity in decomposing beneficial mutations into their proximal, mechanistic determinants. The application of methods and concepts from molecular biology and life history theory to studies of lytic bacteriophages (phages) has allowed them to understand how natural selection sees mutations influencing life history. This work motivated the research presented here, in which we explored whether, under consistent experimental conditions, small differences in the genome of bacteriophage φX174 could lead to altered life history phenotypes among a panel of eight genetically distinct clones. We assessed the clones' phenotypes by applying a novel statistical framework to the results of a serially sampled parallel infection assay, in which we simultaneously inoculated each of a large number of replicate host volumes with ∼1 phage particle. We sequentially plated the volumes over the course of infection and counted the plaques that formed after incubation. These counts served as a proxy for the number of phage particles in a single volume as a function of time. From repeated assays, we inferred significant, genetically determined heterogeneity in lysis time and burst size, including lysis time variance. These findings are interesting in light of the genetic and phenotypic constraints on the single-protein lysis mechanism of φX174. We speculate briefly on the mechanisms underlying our results, and we discuss the potential importance of lysis time variance in viral evolution.

Causes of molecular convergence and parallelism in protein evolution

To what extent is the convergent evolution of protein function attributable to convergent or parallel changes at the amino acid level? The mutations that contribute to adaptive protein evolution may represent a biased subset of all possible beneficial mutations owing to mutation bias and/or variation in the magnitude of deleterious pleiotropy. A key finding is that the fitness effects of amino acid mutations are often conditional on genetic background. This context dependence (epistasis) can reduce the probability of convergence and parallelism because it reduces the number of possible mutations that are unconditionally acceptable in divergent genetic backgrounds. Here, I review factors that influence the probability of replicated evolution at the molecular level.

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.

Ancestral Protein Reconstruction Yields Insights into Adaptive Evolution of Binding Specificity in Solute-Binding Proteins

The promiscuous functions of proteins are an important reservoir of functional novelty in protein evolution, but the molecular basis for binding promiscuity remains elusive. We used ancestral protein reconstruction to experimentally characterize evolutionary intermediates in the functional expansion of the polar amino acid-binding protein family, which has evolved to bind a variety of amino acids with high affinity and specificity. High-resolution crystal structures of an ancestral arginine-binding protein in complex with l-arginine and l-glutamine show that the promiscuous binding of l-glutamine is enabled by multi-scale conformational plasticity, water-mediated interactions, and selection of an alternative conformational substate productive for l-glutamine binding. Evolution of specialized glutamine-binding proteins from this ancestral protein was achieved by displacement of water molecules from the protein-ligand interface, reducing the entropic penalty associated with the promiscuous interaction. These results provide a structural and thermodynamic basis for the co-option of a promiscuous interaction in the evolution of binding specificity.

Catalytic Promiscuity of Ancestral Esterases and Hydroxynitrile Lyases

Catalytic promiscuity is a useful, but accidental, enzyme property, so finding catalytically promiscuous enzymes in nature is inefficient. Some ancestral enzymes were branch points in the evolution of new enzymes and are hypothesized to have been promiscuous. To test the hypothesis that ancestral enzymes were more promiscuous than their modern descendants, we reconstructed ancestral enzymes at four branch points in the divergence hydroxynitrile lyases (HNL's) from esterases ∼ 100 million years ago. Both enzyme types are α/β-hydrolase-fold enzymes and have the same catalytic triad, but differ in reaction type and mechanism. Esterases catalyze hydrolysis via an acyl enzyme intermediate, while lyases catalyze an elimination without an intermediate. Screening ancestral enzymes and their modern descendants with six esterase substrates and six lyase substrates found higher catalytic promiscuity among the ancestral enzymes (P < 0.01). Ancestral esterases were more likely to catalyze a lyase reaction than modern esterases, and the ancestral HNL was more likely to catalyze ester hydrolysis than modern HNL's. One ancestral enzyme (HNL1) along the path from esterase to hydroxynitrile lyases was especially promiscuous and catalyzed both hydrolysis and lyase reactions with many substrates. A broader screen tested mechanistically related reactions that were not selected for by evolution: decarboxylation, Michael addition, γ-lactam hydrolysis and 1,5-diketone hydrolysis. The ancestral enzymes were more promiscuous than their modern descendants (P = 0.04). Thus, these reconstructed ancestral enzymes are catalytically promiscuous, but HNL1 is especially so.

How mutational epistasis impairs predictability in protein evolution and design

There has been much debate about the extent to which mutational epistasis, that is, the dependence of the outcome of a mutation on the genetic background, constrains evolutionary trajectories. The degree of unpredictability introduced by epistasis, due to the non-additivity of functional effects, strongly hinders the strategies developed in protein design and engineering. While many studies have addressed this issue through systematic characterization of evolutionary trajectories within individual enzymes, the field lacks a consensus view on this matter. In this work, we performed a comprehensive analysis of epistasis by analyzing the mutational effects from nine adaptive trajectories toward new enzymatic functions. We quantified epistasis by comparing the effect of mutations occurring between two genetic backgrounds: the starting enzyme (for example, wild type) and the intermediate variant on which the mutation occurred during the trajectory. We found that most trajectories exhibit positive epistasis, in which the mutational effect is more beneficial when it occurs later in the evolutionary trajectory. Approximately half (49%) of functional mutations were neutral or negative on the wild-type background, but became beneficial at a later stage in the trajectory, indicating that these functional mutations were not predictable from the initial starting point. While some cases of strong epistasis were associated with direct interaction between residues, many others were caused by long-range indirect interactions between mutations. Our work highlights the prevalence of epistasis in enzyme adaptive evolution, in particular positive epistasis, and suggests the necessity of incorporating mutational epistasis in protein engineering and design to create highly efficient catalysts.

Causes of evolutionary rate variation among protein sites

It has long been recognized that certain sites within a protein, such as sites in the protein core or catalytic residues in enzymes, are evolutionarily more conserved than other sites. However, our understanding of rate variation among sites remains surprisingly limited. Recent progress to address this includes the development of a wide array of reliable methods to estimate site-specific substitution rates from sequence alignments. In addition, several molecular traits have been identified that correlate with site-specific mutation rates, and novel mechanistic biophysical models have been proposed to explain the observed correlations. Nonetheless, current models explain, at best, approximately 60% of the observed variance, highlighting the limitations of current methods and models and the need for new research 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.

Evolution of a Catalytic Mechanism

The means by which superfamilies of specialized enzymes arise by gene duplication and functional divergence are poorly understood. The escape from adaptive conflict hypothesis, which posits multiple copies of a gene encoding a primitive inefficient and highly promiscuous generalist ancestor, receives support from experiments showing that resurrected ancestral enzymes are indeed more substrate-promiscuous than their modern descendants. Here, we provide evidence in support of an alternative model, the innovation-amplification-divergence hypothesis, which posits a single-copied ancestor as efficient and specific as any modern enzyme. We argue that the catalytic mechanisms of plant esterases and descendent acetone cyanohydrin lyases are incompatible with each other (e.g., the reactive substrate carbonyl must bind in opposite orientations in the active site). We then show that resurrected ancestral plant esterases are as catalytically specific as modern esterases, that the ancestor of modern acetone cyanohydrin lyases was itself only very weakly promiscuous, and that improvements in lyase activity came at the expense of esterase activity. These observations support the innovation-amplification-divergence hypothesis, in which an ancestor gains a weak promiscuous activity that is improved by selection at the expense of the ancestral activity, and not the escape from adaptive conflict in which an inefficient generalist ancestral enzyme steadily loses promiscuity throughout the transition to a highly active specialized modern enzyme.

Convergent Evolution of Hemoglobin Function in High-Altitude Andean Waterfowl Involves Limited Parallelism at the Molecular Sequence Level

A fundamental question in evolutionary genetics concerns the extent to which adaptive phenotypic convergence is attributable to convergent or parallel changes at the molecular sequence level. Here we report a comparative analysis of hemoglobin (Hb) function in eight phylogenetically replicated pairs of high- and low-altitude waterfowl taxa to test for convergence in the oxygenation properties of Hb, and to assess the extent to which convergence in biochemical phenotype is attributable to repeated amino acid replacements. Functional experiments on native Hb variants and protein engineering experiments based on site-directed mutagenesis revealed the phenotypic effects of specific amino acid replacements that were responsible for convergent increases in Hb-O₂ affinity in multiple high-altitude taxa. In six of the eight taxon pairs, high-altitude taxa evolved derived increases in Hb-O₂ affinity that were caused by a combination of unique replacements, parallel replacements (involving identical-by-state variants with independent mutational origins in different lineages), and collateral replacements (involving shared, identical-by-descent variants derived via introgressive hybridization). In genome scans of nucleotide differentiation involving high- and low-altitude populations of three separate species, function-altering amino acid polymorphisms in the globin genes emerged as highly significant outliers, providing independent evidence for adaptive divergence in Hb function. The experimental results demonstrate that convergent changes in protein function can occur through multiple historical paths, and can involve multiple possible mutations. Most cases of convergence in Hb function did not involve parallel substitutions and most parallel substitutions did not affect Hb-O₂ affinity, indicating that the repeatability of phenotypic evolution does not require parallelism at the molecular level.

The convergent evolution of similar traits in different species could be due to repeated changes at the genetic level or different changes that produce the same phenotypic effect. To investigate the extent to which convergence in phenotype is caused by repeated mutations, we investigated the molecular basis of convergent changes in the oxygenation properties of hemoglobin (Hb) in eight pairs of high- and low-altitude waterfowl taxa from the Andes. The results revealed that convergent increases in Hb-O₂ affinity in highland taxa involved a combination of unique and repeated amino acid replacements. However, convergent changes in Hb function generally did not involve parallel substitutions, indicating that repeatability in the evolution of protein function does not require repeatability at the sequence level.

Large-Scale Analysis Exploring Evolution of Catalytic Machineries and Mechanisms in Enzyme Superfamilies

Enzymes, as biological catalysts, form the basis of all forms of life. How these proteins have evolved their functions remains a fundamental question in biology. Over 100 years of detailed biochemistry studies, combined with the large volumes of sequence and protein structural data now available, means that we are able to perform large-scale analyses to address this question. Using a range of computational tools and resources, we have compiled information on all experimentally annotated changes in enzyme function within 379 structurally defined protein domain superfamilies, linking the changes observed in functions during evolution to changes in reaction chemistry. Many superfamilies show changes in function at some level, although one function often dominates one superfamily. We use quantitative measures of changes in reaction chemistry to reveal the various types of chemical changes occurring during evolution and to exemplify these by detailed examples. Additionally, we use structural information of the enzymes active site to examine how different superfamilies have changed their catalytic machinery during evolution. Some superfamilies have changed the reactions they perform without changing catalytic machinery. In others, large changes of enzyme function, in terms of both overall chemistry and substrate specificity, have been brought about by significant changes in catalytic machinery. Interestingly, in some superfamilies, relatives perform similar functions but with different catalytic machineries. This analysis highlights characteristics of functional evolution across a wide range of superfamilies, providing insights that will be useful in predicting the function of uncharacterised sequences and the design of new synthetic enzymes.

•

Examining how enzyme function evolves using sequence, structure, and reaction mechanism data.

•

Quantifying changes in reaction mechanisms reveals how function has diverged in many superfamilies.

•

Homologous domains frequently use different catalytic residues, which sometimes perform the same enzyme chemistry.

•

This large-scale analysis has significance in protein function prediction and enzyme design.

Coevolutionary Landscape Inference and the Context-Dependence of Mutations in Beta-Lactamase TEM-1

The quantitative characterization of mutational landscapes is a task of outstanding importance in evolutionary and medical biology: It is, for example, of central importance for our understanding of the phenotypic effect of mutations related to disease and antibiotic drug resistance. Here we develop a novel inference scheme for mutational landscapes, which is based on the statistical analysis of large alignments of homologs of the protein of interest. Our method is able to capture epistatic couplings between residues, and therefore to assess the dependence of mutational effects on the sequence context where they appear. Compared with recent large-scale mutagenesis data of the beta-lactamase TEM-1, a protein providing resistance against beta-lactam antibiotics, our method leads to an increase of about 40% in explicative power as compared with approaches neglecting epistasis. We find that the informative sequence context extends to residues at native distances of about 20 Å from the mutated site, reaching thus far beyond residues in direct physical contact.

Lessons in Protein Design from Combined Evolution and Conformational Dynamics

Protein-protein interactions play important roles in the control of every cellular process. How natural selection has optimized protein design to produce molecules capable of binding to many partner proteins is a fascinating problem but not well understood. Here, we performed a combinatorial analysis of protein sequence evolution and conformational dynamics to study how calmodulin (CaM), which plays essential roles in calcium signaling pathways, has adapted to bind to a large number of partner proteins. We discovered that amino acid residues in CaM can be partitioned into unique classes according to their degree of evolutionary conservation and local stability. Holistically, categorization of CaM residues into these classes reveals enriched physico-chemical interactions required for binding to diverse targets, balanced against the need to maintain the folding and structural modularity of CaM to achieve its overall function. The sequence-structure-function relationship of CaM provides a concrete example of the general principle of protein design. We have demonstrated the synergy between the fields of molecular evolution and protein biophysics and created a generalizable framework broadly applicable to the study of protein-protein interactions.

The Role of Conformational Dynamics and Allostery in the Disease Development of Human Ferritin

Determining the three-dimensional structure of myoglobin, the first solved structure of a protein, fundamentally changed the way protein function was understood. Even more revolutionary was the information that came afterward: protein dynamics play a critical role in biological functions. Therefore, understanding conformational dynamics is crucial to obtaining a more complete picture of protein evolution. We recently analyzed the evolution of different protein families including green fluorescent proteins (GFPs), β-lactamase inhibitors, and nuclear receptors, and we observed that the alteration of conformational dynamics through allosteric regulation leads to functional changes. Moreover, proteome-wide conformational dynamics analysis of more than 100 human proteins showed that mutations occurring at rigid residue positions are more susceptible to disease than flexible residue positions. These studies suggest that disease-associated mutations may impair dynamic allosteric regulations, leading to loss of function. Thus, in this study, we analyzed the conformational dynamics of the wild-type light chain subunit of human ferritin protein along with the neutral and disease forms. We first performed replica exchange molecular dynamics simulations of wild-type and mutants to obtain equilibrated dynamics and then used perturbation response scanning (PRS), where we introduced a random Brownian kick to a position and computed the fluctuation response of the chain using linear response theory. Using this approach, we computed the dynamic flexibility index (DFI) for each position in the chain for the wild-type and the mutants. DFI quantifies the resilience of a position to a perturbation and provides a flexibility/rigidity measurement for a given position in the chain. The DFI analysis reveals that neutral variants and the wild-type exhibit similar flexibility profiles in which experimentally determined functionally critical sites act as hinges in controlling the overall motion. However, disease mutations alter the conformational dynamic profile, making hinges more loose (i.e., softening the hinges), thus impairing the allosterically regulated dynamics.

Reconstructed Ancestral Enzymes Impose a Fitness Cost upon Modern Bacteria Despite Exhibiting Favourable Biochemical Properties

Ancestral sequence reconstruction has been widely used to study historical enzyme evolution, both from biochemical and cellular perspectives. Two properties of reconstructed ancestral proteins/enzymes are commonly reported--high thermostability and high catalytic activity--compared with their contemporaries. Increased protein stability is associated with lower aggregation rates, higher soluble protein abundance and a greater capacity to evolve, and therefore, these proteins could be considered "superior" to their contemporary counterparts. In this study, we investigate the relationship between the favourable in vitro biochemical properties of reconstructed ancestral enzymes and the organismal fitness they confer in vivo. We have previously reconstructed several ancestors of the enzyme LeuB, which is essential for leucine biosynthesis. Our initial fitness experiments revealed that overexpression of ANC4, a reconstructed LeuB that exhibits high stability and activity, was only able to partially rescue the growth of a ΔleuB strain, and that a strain complemented with this enzyme was outcompeted by strains carrying one of its descendants. When we expanded our study to include five reconstructed LeuBs and one contemporary, we found that neither in vitro protein stability nor the catalytic rate was correlated with fitness. Instead, fitness showed a strong, negative correlation with estimated evolutionary age (based on phylogenetic relationships). Our findings suggest that, for reconstructed ancestral enzymes, superior in vitro properties do not translate into organismal fitness in vivo. The molecular basis of the relationship between fitness and the inferred age of ancestral LeuB enzymes is unknown, but may be related to the reconstruction process. We also hypothesise that the ancestral enzymes may be incompatible with the other, contemporary enzymes of the metabolic network.

Adaptive evolutionary paths from UV reception to sensing violet light by epistatic interactions

Ultraviolet (UV) reception is useful for such basic behaviors as mate choice, foraging, predator avoidance, communication, and navigation, whereas violet reception improves visual resolution and subtle contrast detection. UV and violet reception are mediated by the short wavelength-sensitive (SWS1) pigments that absorb light maximally (λmax) at ~360 nm and ~395 to 440 nm, respectively. Because of strong nonadditive (epistatic) interactions among amino acid changes in the pigments, the adaptive evolutionary mechanisms of these phenotypes are not well understood. Evolution of the violet pigment of the African clawed frog (Xenopus laevis, λmax = 423 nm) from the UV pigment in the amphibian ancestor (λmax = 359 nm) can be fully explained by eight mutations in transmembrane (TM) I-III segments. We show that epistatic interactions involving the remaining TM IV-VII segments provided evolutionary potential for the frog pigment to gradually achieve its violet-light reception by tuning its color sensitivity in small steps. Mutants in these segments also impair pigments that would cause drastic spectral shifts and thus eliminate them from viable evolutionary pathways. The overall effects of epistatic interactions involving TM IV-VII segments have disappeared at the last evolutionary step and thus are not detectable by studying present-day pigments. Therefore, characterizing the genotype-phenotype relationship during each evolutionary step is the key to uncover the true nature of epistasis.

Recapitulating the Structural Evolution of Redox Regulation in Adenosine 5'-Phosphosulfate Kinase from Cyanobacteria to Plants

In plants, adenosine 5'-phosphosulfate (APS) kinase (APSK) is required for reproductive viability and the production of 3'-phosphoadenosine 5'-phosphosulfate (PAPS) as a sulfur donor in specialized metabolism. Previous studies of the APSK from Arabidopsis thaliana (AtAPSK) identified a regulatory disulfide bond formed between the N-terminal domain (NTD) and a cysteine on the core scaffold. This thiol switch is unique to mosses, gymnosperms, and angiosperms. To understand the structural evolution of redox control of APSK, we investigated the redox-insensitive APSK from the cyanobacterium Synechocystis sp. PCC 6803 (SynAPSK). Crystallographic analysis of SynAPSK in complex with either APS and a non-hydrolyzable ATP analog or APS and sulfate revealed the overall structure of the enzyme, which lacks the NTD found in homologs from mosses and plants. A series of engineered SynAPSK variants reconstructed the structural evolution of the plant APSK. Biochemical analyses of SynAPSK, SynAPSK H23C mutant, SynAPSK fused to the AtAPSK NTD, and the fusion protein with the H23C mutation showed that the addition of the NTD and cysteines recapitulated thiol-based regulation. These results reveal the molecular basis for structural changes leading to the evolution of redox control of APSK in the green lineage from cyanobacteria to plants.

Reverse evolution leads to genotypic incompatibility despite functional and active site convergence

Understanding the extent to which enzyme evolution is reversible can shed light on the fundamental relationship between protein sequence, structure, and function. Here, we perform an experimental test of evolutionary reversibility using directed evolution from a phosphotriesterase to an arylesterase, and back, and examine the underlying molecular basis. We find that wild-type phosphotriesterase function could be restored (>10⁴-fold activity increase), but via an alternative set of mutations. The enzyme active site converged towards its original state, indicating evolutionary constraints imposed by catalytic requirements. We reveal that extensive epistasis prevents reversions and necessitates fixation of new mutations, leading to a functionally identical sequence. Many amino acid exchanges between the new and original enzyme are not tolerated, implying sequence incompatibility. Therefore, the evolution was phenotypically reversible but genotypically irreversible. Our study illustrates that the enzyme's adaptive landscape is highly rugged, and different functional sequences may constitute separate fitness peaks.

DOI:http://dx.doi.org/10.7554/eLife.06492.001

Enzymes in bacteria and other organisms are built following instructions contained within each cell's DNA. Changes in the DNA, that is to say, mutations, can alter the shape and activity of the enzymes that are produced, which can ultimately affect the ability of the organism to survive and reproduce. Mutations that are beneficial to the organism are more likely to be passed on to future generations, which can lead to populations changing over time.

The DNA sequences that an organism carries are referred to as its ‘genotype’ and the resulting physical characteristics of the organism are known as its ‘phenotype’. Studies of evolution tend to focus on how particular species or molecules become more different over time. However, one area that remains controversial is whether it is possible for evolution to be reversed so that an organism or molecule returns to a previous form.

An enzyme called PTE is said to have phosphotriesterase activity because it catalyzes this particular type of chemical reaction. Recently, a group of researchers used a method called ‘directed evolution’ to demonstrate that it is possible for PTE to evolve in a way that means it loses its phosphotriesterase activity and becomes able to catalyze a different type of chemical reaction. Here, Kaltenbach et al.—including some of the researchers from the previous work—investigated whether it was possible to use the same method to reverse this evolution and restore the enzyme's original activity.

The experiments show that reverse evolution is possible as phosphotriesterase activity was restored to the PTE enzyme from the previous study. However, although the phenotype of the final enzyme matched that of the original PTE enzyme, the genotypes did not match as the DNA sequences of the genes that encode these enzymes differ. The DNA does not revert to its original sequence because the effect of individual mutations on the phenotype depends on what other mutations are present. For example, as the enzyme evolved its new activity, additional mutations accumulated that did not alter enzyme activity. During the reverse evolution experiment, some of these mutations could have started to exert influence on the phenotype so that different mutations were required to restore the phosphotriesterase activity.

In the future, Kaltenbach et al.'s findings may aid efforts to engineer artificial enzymes for use in medicine or industry.

DOI:http://dx.doi.org/10.7554/eLife.06492.002

Genetic approaches in comparative and evolutionary physiology

Whole animal physiological performance is highly polygenic and highly plastic, and the same is generally true for the many subordinate traits that underlie performance capacities. Quantitative genetics, therefore, provides an appropriate framework for the analysis of physiological phenotypes and can be used to infer the microevolutionary processes that have shaped patterns of trait variation within and among species. In cases where specific genes are known to contribute to variation in physiological traits, analyses of intraspecific polymorphism and interspecific divergence can reveal molecular mechanisms of functional evolution and can provide insights into the possible adaptive significance of observed sequence changes. In this review, we explain how the tools and theory of quantitative genetics, population genetics, and molecular evolution can inform our understanding of mechanism and process in physiological evolution. For example, lab-based studies of polygenic inheritance can be integrated with field-based studies of trait variation and survivorship to measure selection in the wild, thereby providing direct insights into the adaptive significance of physiological variation. Analyses of quantitative genetic variation in selection experiments can be used to probe interrelationships among traits and the genetic basis of physiological trade-offs and constraints. We review approaches for characterizing the genetic architecture of physiological traits, including linkage mapping and association mapping, and systems approaches for dissecting intermediary steps in the chain of causation between genotype and phenotype. We also discuss the promise and limitations of population genomic approaches for inferring adaptation at specific loci. We end by highlighting the role of organismal physiology in the functional synthesis of evolutionary biology.

Site-Specific Amino Acid Preferences Are Mostly Conserved in Two Closely Related Protein Homologs

Evolution drives changes in a protein’s sequence over time. The extent to which these changes in sequence lead to shifts in the underlying preference for each amino acid at each site is an important question with implications for comparative sequence-analysis methods, such as molecular phylogenetics. To quantify the extent that site-specific amino acid preferences shift during evolution, we performed deep mutational scanning on two homologs of human influenza nucleoprotein with 94% amino acid identity. We found that only a modest fraction of sites exhibited shifts in amino acid preferences that exceeded the noise in our experiments. Furthermore, even among sites that did exhibit detectable shifts, the magnitude tended to be small relative to differences between nonhomologous proteins. Given the limited change in amino acid preferences between these close homologs, we tested whether our measurements could inform site-specific substitution models that describe the evolution of nucleoproteins from more diverse influenza viruses. We found that site-specific evolutionary models informed by our experiments greatly outperformed nonsite-specific alternatives in fitting phylogenies of nucleoproteins from human, swine, equine, and avian influenza. Combining the experimental data from both homologs improved phylogenetic fit, partly because measurements in multiple genetic contexts better captured the evolutionary average of the amino acid preferences for sites with shifting preferences. Our results show that site-specific amino acid preferences are sufficiently conserved that measuring mutational effects in one protein provides information that can improve quantitative evolutionary modeling of nearby homologs.