In vivo molecular evolution reveals biophysical origins of organismal fitness

Experimental evolution, loss-of-function mutations, and "the first rule of adaptive evolution"

Adaptive evolution can cause a species to gain, lose, or modify a function; therefore, it is of basic interest to determine whether any of these modes dominates the evolutionary process under particular circumstances. Because mutation occurs at the molecular level, it is necessary to examine the molecular changes produced by the underlying mutation in order to assess whether a given adaptation is best considered as a gain, loss, or modification of function. Although that was once impossible, the advance of molecular biology in the past half century has made it feasible. In this paper, I review molecular changes underlying some adaptations, with a particular emphasis on evolutionary experiments with microbes conducted over the past four decades. I show that by far the most common adaptive changes seen in those examples are due to the loss or modification of a pre-existing molecular function, and I discuss the possible reasons for the prominence of such mutations.

Experimental evolution of adenylate kinase reveals contrasting strategies toward protein thermostability

Success in evolution depends critically upon the ability of organisms to adapt, a property that is also true for the proteins that contribute to the fitness of an organism. Successful protein evolution is enhanced by mutational pathways that generate a wide range of physicochemical mechanisms to adaptation. In an earlier study, we used a weak-link method to favor changes to an essential but maladapted protein, adenylate kinase (AK), within a microbial population. Six AK mutants (a single mutant followed by five double mutants) had success within the population, revealing a diverse range of adaptive strategies that included changes in nonpolar packing, protein folding dynamics, and formation of new hydrogen bonds and electrostatic networks. The first mutation, AK(BSUB) Q199R, was essential in defining the structural context that facilitated subsequent mutations as revealed by a considerable mutational epistasis and, in one case, a very strong dependence upon the order of mutations. Namely, whereas the single mutation AK(BSUB) G213E decreases protein stability by >25 degrees C, the same mutation in the background of AK(BSUB) Q199R increases stability by 3.4 degrees C, demonstrating that the order of mutations can play a critical role in favoring particular molecular pathways to adaptation. In turn, protein folding kinetics shows that four of the five AK(BSUB) double mutants utilize a strategy in which an increase in the folding rate accompanied by a decrease in the unfolding rate results in additional stability. However, one mutant exhibited a dramatic increase in the folding relative to a modest increase in the unfolding rate, suggesting a different adaptive strategy for thermostability. In all cases, an increase in the folding rates for the double mutants appears to be the preferred mechanism in conferring additional stability and may be an important aspect of protein evolution. The range of overlapping as well as contrasting strategies for success illustrates both the power and subtlety of adaptation at even the smallest unit of change, a single amino acid.

Evolutionary fates within a microbial population highlight an essential role for protein folding during natural selection

Systems biology can offer a great deal of insight into evolution by quantitatively linking complex properties such as protein structure, folding, and function to the fitness of an organism. Although the link between diseases such as Alzheimer's and misfolding is well appreciated, directly showing the importance of protein folding to success in evolution has been more difficult. We show here that predicting success during adaptation can depend critically on enzyme kinetic and folding models. We used a 'weak link' method to favor mutations to an essential, but maladapted, adenylate kinase gene within a microbial population that resulted in the identification of five mutants that arose nearly simultaneously and competed for success. Physicochemical characterization of these mutants showed that, although steady-state enzyme activity is important, success within the population is critically dependent on resistance to denaturation and aggregation. A fitness function based on in vitro measurements of enzyme activity, reversible and irreversible unfolding, and the physiological context reproduces in vivo evolutionary fates in the population linking organismal adaptation to its physical basis.

Thermal adaptation of viruses and bacteria

A previously established multiscale population genetics model posits that fitness can be inferred from the physical properties of proteins under the physiological assumption that a loss of stability by any protein confers the lethal phenotype to an organism. Here, we develop this model further by positing that replication rate (fitness) of a bacterial or viral strain directly depends on the copy number of folded proteins, which determine its replication rate. Using this model, and both numerical and analytical approaches, we studied the adaptation process of bacteria and viruses at varied environmental temperatures. We found that a broad distribution of protein stabilities observed in the model and in experiment is the key determinant of thermal response for viruses and bacteria. Our results explain most of the earlier experimental observations: the striking asymmetry of thermal response curves; the absence of evolutionary tradeoff, which was expected but not found in experiments; correlation between denaturation temperature for several protein families and the optimal growth temperature of their carrier organisms; and proximity of bacterial or viral optimal growth temperatures to their evolutionary temperatures. Our theory quantitatively and with high accuracy described thermal response curves for 35 bacterial species using, for each species, only two adjustable parameters-the number of rate-determining genes and the energy barrier for metabolic reactions.

Evolution of a single gene highlights the complexity underlying molecular descriptions of fitness

Evolution by natural selection is the driving force behind the endless variation we see in nature, yet our understanding of how changes at the molecular level give rise to different phenotypes and altered fitness at the population level remains inadequate. The reproductive fitness of an organism is the most basic metric that describes the chance that an organism will succeed or fail in its environment and it depends upon a complex network of inter- and intramolecular interactions. A deeper understanding of the quantitative relationships relating molecular evolution to adaptation, and consequently fitness, can guide our understanding of important issues in biomedicine such as drug resistance and the engineering of new organisms with applications to biotechnology. We have developed the "weak link" approach to determine how changes in molecular structure and function can relate to fitness and evolutionary outcomes. By replacing adenylate kinase (AK), an essential gene, in a thermophile with a homologous AK from a mesophile we have created a maladapted weak link that produces a temperature-sensitive phenotype. The recombinant strain adapts to nonpermissive temperatures through point mutations to the weak link that increase both stability and activity of the enzyme AK at higher temperatures. Here, we propose a fitness function relating enzyme activity to growth rate and use it to create a dynamic model of a population of bacterial cells. Using metabolic control analysis we show that the growth rate exhibits thresholdlike behavior, saturating at high enzyme activity as other reactions in the energy metabolism pathway become rate limiting. The dynamic model accurately recapitulates observed evolutionary outcomes. These findings suggest that in vitro enzyme kinetic data, in combination with metabolic network analysis, can be used to create fitness functions and dynamic models of evolution within simple metabolic systems.

Crystal structure of a trimeric archaeal adenylate kinase from the mesophile Methanococcus maripaludis with an unusually broad functional range and thermal stability

The structure of the trimeric adenylate kinase from the Archaebacteria Methanococcus mariplaludis (AK(MAR)) has been solved to 2.5-A resolution and the temperature dependent stability and kinetics of the enzyme measured. The K(M) and V(max) of AK(MAR) exhibit only modest temperature dependence from 30 degrees -60 degrees C. Although M. mariplaludis is a mesophile with a maximum growth temperature of 43 degrees C, AK(MAR) has a very broad functional range and stability (T(m) = 74.0 degrees C) that are more consistent with a thermophilic enzyme with high thermostability and exceptional activity over a wide range of temperatures, suggesting that this microbe may have only recently invaded a mesophilic niche and has yet to fully adapt. A comparison of the Local Structural Entropy (LSE) for AK(MAR) to the related adenylate kinases from the mesophile Methanococcus voltae and thermophile Methanococcus thermolithotrophicus show that changes in LSE are able to fully account for the intermediate stability of AK(MAR) and highlights a general mechanism for protein adaptation in this class of enzymes.

Exploring protein fitness landscapes by directed evolution

Directed evolution circumvents our profound ignorance of how a protein's sequence encodes its function by using iterative rounds of random mutation and artificial selection to discover new and useful proteins. Proteins can be tuned to adapt to new functions or environments by simple adaptive walks involving small numbers of mutations. Directed evolution studies have shown how rapidly some proteins can evolve under strong selection pressures and, because the entire 'fossil record' of evolutionary intermediates is available for detailed study, they have provided new insight into the relationship between sequence and function. Directed evolution has also shown how mutations that are functionally neutral can set the stage for further adaptation.

Zipping and unzipping of adenylate kinase: atomistic insights into the ensemble of open<-->closed transitions

Adenylate kinase (AdK), a phosphotransferase enzyme, plays an important role in cellular energy homeostasis. It undergoes a large conformational change between an open and a closed state, even in the absence of substrate. We investigate the apo-AdK transition at the atomic level both with free-energy calculations and with our new dynamic importance sampling (DIMS) molecular dynamics method. DIMS is shown to sample biologically relevant conformations as verified by comparing an ensemble of hundreds of DIMS transitions to AdK crystal structure intermediates. The simulations reveal in atomic detail how hinge regions partially and intermittently unfold during the transition. Conserved salt bridges are seen to have important structural and dynamic roles; in particular, four ionic bonds that open in a sequential, zipper-like fashion and, thus, dominate the free-energy landscape of the transition are identified. Transitions between the closed and open conformations only have to overcome moderate free-energy barriers. Unexpectedly, the closed state and the open state encompass broad free-energy basins that contain conformations differing in domain hinge motions by up to 40 degrees . The significance of these extended states is discussed in relation to recent experimental Förster resonance energy transfer measurements. Taken together, these results demonstrate how a small number of cooperative key interactions can shape the overall dynamics of an enzyme and suggest an "all-or-nothing" mechanism for the opening and closing of AdK. Our efficient DIMS molecular dynamics computer simulation approach can provide a detailed picture of a functionally important macromolecular transition and thus help to interpret and suggest experiments to probe the conformational landscape of dynamic proteins such as AdK.

Structure and biochemical characterization of an adenylate kinase originating from the psychrophilic organism Marinibacillus marinus

Adenylate kinases (AKs; EC 2.7.4.3) are essential members of the NMP kinase family that maintain cellular homeostasis by the interconversion of AMP, ADP and ATP. AKs play a critical role in adenylate homeostasis across all domains of life and have been used extensively as prototypes for the study of protein adaptation and the relationship of protein dynamics and stability to function. To date, kinetic studies of psychrophilic AKs have not been performed. In order to broaden understanding of extremophilic adaptation, the kinetic parameters of adenylate kinase from the psychrophile Marinibacillus marinus were examined and the crystal structure of this cold-adapted enzyme was determined at 2.0 A resolution. As expected, the overall structure and topology of the psychrophilic M. marinus AK are similar to those of mesophilic and thermophilic AKs. The thermal denaturation midpoint of M. marinus AK (321.1 K) is much closer to that of the mesophile Bacillus subtilis (320.7 K) than the more closely related psychrophile B. globisporus (316.4 K). In addition, the enzymatic properties of M. marinus AK are quite close to those of the mesophilic AK and suggests that M. marinus experiences temperature ranges in which excellent enzyme function over a broad temperature range (293-313 K) has been retained for the success of the organism. Even transient loss of AK function is lethal and as a consequence AK must be robust and be well adapted to the environment of the host organism.

Lethal mutagenesis in viruses and bacteria

In this work we study how mutations that change physical properties of cell proteins (stability) affect population survival and growth. We present a model in which the genotype is presented as a set folding free energies of cell proteins. Mutations occur upon replication, so stabilities of some proteins in daughter cells differ from those in the parent cell by amounts deduced from the distribution of mutational effects on protein stability. The genotype-phenotype relationship posits that the cell's fitness (replication rate) is proportional to the concentration of its folded proteins and that unstable essential proteins result in lethality. Simulations reveal that lethal mutagenesis occurs at a mutation rate close to seven mutations in each replication of the genome for RNA viruses and at about half that rate for DNA-based organisms, in accord with earlier predictions from analytical theory and experimental results. This number appears somewhat dependent on the number of genes in the organisms and the organism's natural death rate. Further, our model reproduces the distribution of stabilities of natural proteins, in excellent agreement with experiments. We find that species with high mutation rates tend to have less stable proteins compared to species with low mutation rates.

Direct Mg(2+) binding activates adenylate kinase from Escherichia coli

We report evidence that adenylate kinase (AK) from Escherichia coli can be activated by the direct binding of a magnesium ion to the enzyme, in addition to ATP-complexed Mg(2+). By systematically varying the concentrations of AMP, ATP, and magnesium in kinetic experiments, we found that the apparent substrate inhibition of AK, formerly attributed to AMP, was suppressed at low magnesium concentrations and enhanced at high magnesium concentrations. This previously unreported magnesium dependence can be accounted for by a modified random bi-bi model in which Mg(2+) can bind to AK directly prior to AMP binding. A new kinetic model is proposed to replace the conventional random bi-bi mechanism with substrate inhibition and is able to describe the kinetic data over a physiologically relevant range of magnesium concentrations. According to this model, the magnesium-activated AK exhibits a 23- +/- 3-fold increase in its forward reaction rate compared with the unactivated form. The findings imply that Mg(2+) could be an important affecter in the energy signaling network in cells.

Predictability of phenotypic differentiation across flow regimes in fishes

Fish inhabit environments greatly varying in intensity of water velocity, and these flow regimes are generally believed to be of major evolutionary significance. To what extent does water flow drive repeatable and predictable phenotypic differentiation? Although many investigators have examined phenotypic variation across flow gradients in fishes, no clear consensus regarding the nature of water velocity's effects on phenotypic diversity has yet emerged. Here, I describe a generalized model that produces testable hypotheses of morphological and locomotor differentiation between flow regimes in fishes. The model combines biomechanical information (describing how fish morphology determines locomotor abilities) with ecological information (describing how locomotor performance influences fitness) to yield predictions of divergent natural selection and phenotypic differentiation between low-flow and high-flow environments. To test the model's predictions of phenotypic differentiation, I synthesized the existing literature and conducted a meta-analysis. Based on results gathered from 80 studies, providing 115 tests of predictions, the model produced some accurate results across both intraspecific and interspecific scales, as differences in body shape, caudal fin shape, and steady-swimming performance strongly matched predictions. These results suggest that water flow drives predictable phenotypic variation in disparate groups of fish based on a common, generalized model, and that microevolutionary processes might often scale up to generate broader, interspecific patterns. However, too few studies have examined differentiation in body stiffness, muscle architecture, or unsteady-swimming performance to draw clear conclusions for those traits. The analysis revealed that, at the intraspecific scale, both genetic divergence and phenotypic plasticity play important roles in phenotypic differentiation across flow regimes, but we do not yet know the relative importance of these two sources of phenotypic variation. Moreover, while major patterns within and between species were predictable, we have little direct evidence regarding the role of water flow in driving speciation or generating broad, macroevolutionary patterns, as too few studies have addressed these topics or conducted analyses within a phylogenetic framework. Thus, flow regime does indeed drive some predictable phenotypic outcomes, but many questions remain unanswered. This study establishes a general model for predicting phenotypic differentiation across flow regimes in fishes, and should help guide future studies in fruitful directions, thereby enhancing our understanding of the predictability of phenotypic variation in nature.

Neutrality and robustness in evo-devo: emergence of lateral inhibition

Embryonic development is defined by the hierarchical dynamical process that translates genetic information (genotype) into a spatial gene expression pattern (phenotype) providing the positional information for the correct unfolding of the organism. The nature and evolutionary implications of genotype-phenotype mapping still remain key topics in evolutionary developmental biology (evo-devo). We have explored here issues of neutrality, robustness, and diversity in evo-devo by means of a simple model of gene regulatory networks. The small size of the system allowed an exhaustive analysis of the entire fitness landscape and the extent of its neutrality. This analysis shows that evolution leads to a class of robust genetic networks with an expression pattern characteristic of lateral inhibition. This class is a repertoire of distinct implementations of this key developmental process, the diversity of which provides valuable clues about its underlying causal principles.

Optimal design of thermally stable proteins

Motivation: For many biotechnological purposes, it is desirable to redesign proteins to be more structurally and functionally stable at higher temperatures. For example, chemical reactions are intrinsically faster at higher temperatures, so using enzymes that are stable at higher temperatures would lead to more efficient industrial processes. We describe an innovative and computationally efficient method called Improved Configurational Entropy (ICE), which can be used to redesign a protein to be more thermally stable (i.e. stable at high temperatures). This can be accomplished by systematically modifying the amino acid sequence via local structural entropy (LSE) minimization. The minimization problem is modeled as a shortest path problem in an acyclic graph with nonnegative weights and is solved efficiently using Dijkstra's method.

Contact:mitchell@biochem.wisc.edu

Bioinformatic method for protein thermal stabilization by structural entropy optimization

Engineering proteins for higher thermal stability is an important and difficult challenge. We describe a bioinformatic method incorporating sequence alignments to redesign proteins to be more stable through optimization of local structural entropy. Using this method, improved configurational entropy (ICE), we were able to design more stable variants of a mesophilic adenylate kinase with only the sequence information of one psychrophilic homologue. The redesigned proteins display considerable increases in their thermal stabilities while still retaining catalytic activity. ICE does not require a three-dimensional structure or a large number of homologous sequences, indicating a broad applicability of this method. Our results also highlight the importance of entropy in the stability of protein structures.

A study in molecular contingency: glutamine phosphoribosylpyrophosphate amidotransferase is a promiscuous and evolvable phosphoribosylanthranilate isomerase

The prevalence of paralogous enzymes implies that novel catalytic functions can evolve on preexisting protein scaffolds. The weak secondary activities of proteins, which reflect catalytic promiscuity and substrate ambiguity, are plausible starting points for this evolutionary process. In this study, we observed the emergence of a new enzyme from the ASKA (A Complete Set of E. coli K-12 ORF Archive) collection of Escherichia coli open reading frames. The overexpression of (His)(6)-tagged glutamine phosphoribosylpyrophosphate amidotransferase (PurF) unexpectedly rescued a Delta trpF E. coli strain from starvation on minimal media. The wild-type PurF and TrpF enzymes are unrelated in sequence, tertiary structure and catalytic mechanism. The promiscuous phosphoribosylanthranilate isomerase activity of the ASKA PurF variant apparently stems from a preexisting affinity for phosphoribosylated substrates. The relative fitness of the (His)(6)-PurF/Delta trpF strain was improved 4.8-fold to nearly wild-type levels by random mutagenesis of purF and genetic selection. The evolved and ancestral PurF proteins were purified and reacted with phosphoribosylanthranilate in vitro. The best evolvant (k(cat)/K(M)=0.3 s(-1) M(-1)) was approximately 25-fold more efficient than its ancestor but >10(7)-fold less efficient than the wild-type phosphoribosylanthranilate isomerase. These observations demonstrate in quantitative terms that the weak secondary activities of promiscuous enzymes can dramatically improve the fitness of contemporary organisms.

Genetic constraints on protein evolution

Evolution requires the generation and optimization of new traits ("adaptation") and involves the selection of mutations that improve cellular function. These mutations were assumed to arise by selection of neutral mutations present at all times in the population. Here we review recent evidence that indicates that deleterious mutations are more frequent in the population than previously recognized and that these mutations play a significant role in protein evolution through continuous positive selection. Positively selected mutations include adaptive mutations, i.e. mutations that directly affect enzymatic function, and compensatory mutations, which suppress the pleiotropic effects of adaptive mutations. Compensatory mutations are by far the most frequent of the two and would allow potentially adaptive but deleterious mutations to persist long enough in the population to be positively selected during episodes of adaptation. Compensatory mutations are, by definition, context-dependent and thus constrain the paths available for evolution. This provides a mechanistic basis for the examples of highly constrained evolutionary landscapes and parallel evolution reported in natural and experimental populations. The present review article describes these recent advances in the field of protein evolution and discusses their implications for understanding the genetic basis of disease and for protein engineering in vitro.

Persistence of foot-and-mouth disease virus in cell culture revisited: implications for contingency in evolution

If we could rewind the tape of evolution and play it again, would it turn out to be similar to or different from what we know? Obviously, this key question can only be addressed by fragmentary experimental approaches. Twenty-two years ago, we described the establishment of BHK-21 cells persistently infected with foot-and-mouth disease virus (FMDV), a system that displayed as its major biological feature a coevolution of the cells and the resident virus in the course of persistence. Now we report the establishment of two persistently infected cell lines in parallel, starting with the same clones of FMDV and BHK-21 cells used 22 years ago. We have asked whether the evolution of the two newly established cell lines and of the earlier cell line would be similar or different. The main conclusions of the study are: (i) the basic behaviour characterized by virus–cell coevolution is similar in the three carrier cell lines, despite differences in some genetic alterations of FMDV; (ii) a strikingly parallel behaviour has been observed with the two newly established cell lines passaged in parallel, unveiling a deterministic virus behaviour during persistence; and (iii) selective RT-PCR amplifications have detected imbalances in the proportion of positive- versus negative-strand viral RNA, mediated by both viral and cellular factors. The results confirm coevolution of cells and virus as a major and reproducible feature of FMDV persistence in cell culture, and suggest that rapidly evolving viruses may constitute adequate test systems to probe the influence of historical contingency on evolutionary events.

Empirical fitness landscapes reveal accessible evolutionary paths

When attempting to understand evolution, we traditionally rely on analysing evolutionary outcomes, despite the fact that unseen intermediates determine its course. A handful of recent studies has begun to explore these intermediate evolutionary forms, which can be reconstructed in the laboratory. With this first view on empirical evolutionary landscapes, we can now finally start asking why particular evolutionary paths are taken.

Thermostability/infectivity defect caused by deletion of the core protein V gene in human adenovirus type 5 is rescued by thermo-selectable mutations in the core protein X precursor

Mastadenoviruses represent one of the four major genera of the Adenoviridae family comprising a variety of mammalian pathogens including human adenovirus (Ad), whose genomes encode a gene for minor core protein V (pV), not found in other genera of Adenoviridae. Deletion of other genus-specific genes (gene IX and E3 genes) from the Ad type 5 (Ad5) genome has been studied experimentally in vitro and the results on biological characterization of the mutants support the phylogenetic evidence of those genes being non-essential for Ad viability. On this basis it seemed logical to suggest that a deletion of gene V from the Ad5 genome could also be tolerated. To test this hypothesis we constructed and rescued the first pV-deletion mutant of human Ad5. As compared to Ad5, this mutant formed small plaques, had dramatically reduced thermostability and lower infectivity. A subsequent thermoselection screen of the pV-deleted Ad5 allowed isolation of a suppressor mutant Ad5-dV/TSB with restored biological characteristics. Since replication and viral assembly of Ad5-dV/TSB could still occur in the absence of pV, we conclude that pV is a non-essential component of the virion. The observed rescue of the biological defects appears to be associated with a cluster of point mutations in the gene encoding the precursor for the other core protein, X/Mu. This finding, thus, suggests possible roles of pV and protein X/Mu precursor in viral assembly. It also provides an interesting insight into genetic events that mediate molecular adaptation of viruses to possible changes in the genetic background in the course of their evolutionary divergence. The possible mechanism of the observed genetic suppression is discussed.

Towards genetic markers in animal populations as biomonitors for human-induced environmental change

Genetic markers provide potentially sensitive indicators of changes in environmental conditions because the genetic constitution of populations is normally altered well before populations become extinct. Genetic indicators in populations include overall genetic diversity, genetic changes in traits measured at the phenotypic level, and evolution at specific loci under selection. While overall genetic diversity has rarely been successfully related to environmental conditions, genetically based changes in traits have now been linked to the presence of toxins and both local and global temperature shifts. Candidate loci for monitoring stressors are emerging from information on how specific genes influence traits, and from screens of random loci across environmental gradients. Drosophila research suggests that chromosomal regions under recent intense selection can be identified from patterns of molecular variation and a high frequency of transposable element insertions. Allele frequency changes at candidate loci have been linked to pesticides, pollutants and climate change. Nevertheless, there are challenges in interpreting allele frequencies in populations, particularly when a large number of loci control a trait and when interactions between alleles influence trait expression. To meet these challenges, population samples should be collected for longitudinal studies, and experimental programmes should be undertaken to link variation at candidate genes to ecological processes.