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.     

Accessory mutations balance the marginal stability of the HIV-1 protease in drug resistance

The HIV-1 protease is a major target of inhibitor drugs in AIDS therapies. The therapies are impaired by mutations of the HIV-1 protease that can lead to resistance to protease inhibitors. These mutations are classified into major mutations, which usually occur first and clearly reduce the susceptibility to protease inhibitors, and minor, accessory mutations that occur later and individually do not substantially affect the susceptibility to inhibitors. Major mutations are predominantly located in the active site of the HIV-1 protease and can directly interfere with inhibitor binding. Minor mutations, in contrast, are typically located distal to the active site. A central question is how these distal mutations contribute to resistance development. In this article, we present a systematic computational investigation of stability changes caused by major and minor mutations of the HIV-1 protease. As most small single-domain proteins, the HIV-1 protease is only marginally stable. Mutations that destabilize the folded, active state of the protease therefore can shift the conformational equilibrium towards the unfolded, inactive state. We find that the most frequent major mutations destabilize the HIV-1 protease, whereas roughly half of the frequent minor mutations are stabilizing. An analysis of protease sequences from patients in treatment indicates that the stabilizing minor mutations are frequently correlated with destabilizing major mutations, and that highly resistant HIV-1 proteases exhibit significant fractions of stabilizing mutations. Our results thus indicate a central role of minor mutations in balancing the marginal stability of the protease against the destabilization induced by the most frequent major mutations.     

Mutations that eliminate the requirement for the vertex protein in bacteriophage T4 capsid assembly.

The capsid of bacteriophage T4 is composed of two essential structural proteins, gp23, the major constituent of the capsid, and gp24, a less prevalent protein that is located in the pentameric vertices of the capsid. gp24 is required both to stabilize the capsid and to allow it to be further matured. This requirement can be eliminated by bypass-24 (byp24) mutations within g23. We have isolated, cloned and sequenced several new byp24 mutations. These mutations are cold-sensitive in the absence of gp24, and are located in regions of g23 not known to contain any other mutations affecting capsid assembly. The cold-sensitivity of the byp24 mutations can be reduced by further mutations within g23 (trb mutations). Cloning and sequencing of these trb mutations has revealed that they lie in regions of g23 that contain clusters of mutations that cause the production of high levels of petite and giant phage (ptg mutations). Despite the proximity of the trb mutations to the ptg mutations, none of the ptg mutations has a Trb phenotype. The mutation ptE920g, which is also located near one of the ptg clusters, and which produces only petite and wild-type phage, has been shown to confer a Trb but not a Byp24 phenotype. The relevance of these observations to our understanding of capsid assembly is discussed.     

The Let-60 Locus Controls the Switch between Vulval and Nonvulval Cell Fates in Caenorhabditis Elegans

During induction of the Caenorhabditis elegans hermaphrodite vulva by the anchor cell of the gonad, six multipotent vulval precursor cells (VPCs) have two distinct fates: three VPCs generate the vulva and the other three VPCs generate nonspecialized hypodermis. Genes that control the fates of the VPCs in response to the anchor cell signal are defined by mutations that cause all six VPCs to generate vulval tissue (Multivulva or Muv) or that cause all six VPCs to generate hypodermis (Vulvaless or Vul). Seven dominant Vul mutations were isolated as dominant suppressors of a lin-15 Muv mutation. These mutations are dominant alleles of the gene let-60, previously identified only by recessive lethal mutations. Our genetic studies of these dominant Vul recessive lethal mutations, recessive lethal mutations, intragenic revertants of the dominant Vul mutations, and the closely mapping semi-dominant multivulva lin-34 mutations suggest that: (1) loss-of-function mutations of let-60 are recessive lethal at a larval stage, but they also cause a Vul phenotype if the lethality is rescued maternally by a lin-34 gain-of-function mutation. (2) The dominant Vul alleles of let-60 are dominant negative mutations whose gene products compete with wild-type activity. (3) lin-34 semidominant Muv alleles are either gain-of-function mutations of let-60 or gain-of-function mutations of an intimately related gene that elevates let-60 activity. We propose that let-60 activity controls VPC fates. In a wild-type animal, reception by a VPC of inductive signal activates let-60, and it generates into a vulval cell type; in absence of inductive signal, let-60 activity is low and the VPC generates hypodermal cells. Our genetic interaction studies suggest that let-60 acts downstream of let-23 and lin-15 and upstream of lin-1 and lin-12 in the genetic pathway specifying the switch between vulval and nonvulval cell types.     

Genetic Constraints on Protein Evolution

ABSTRACT

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.     

How protein stability and new functions trade off

Numerous studies have noted that the evolution of new enzymatic specificities is accompanied by loss of the protein's thermodynamic stability (DeltaDeltaG), thus suggesting a tradeoff between the acquisition of new enzymatic functions and stability. However, since most mutations are destabilizing (DeltaDeltaG>0), one should ask how destabilizing mutations that confer new or altered enzymatic functions relative to all other mutations are. We applied DeltaDeltaG computations by FoldX to analyze the effects of 548 mutations that arose from the directed evolution of 22 different enzymes. The stability effects, location, and type of function-altering mutations were compared to DeltaDeltaG changes arising from all possible point mutations in the same enzymes. We found that mutations that modulate enzymatic functions are mostly destabilizing (average DeltaDeltaG = +0.9 kcal/mol), and are almost as destabilizing as the "average" mutation in these enzymes (+1.3 kcal/mol). Although their stability effects are not as dramatic as in key catalytic residues, mutations that modify the substrate binding pockets, and thus mediate new enzymatic specificities, place a larger stability burden than surface mutations that underline neutral, non-adaptive evolutionary changes. How are the destabilizing effects of functional mutations balanced to enable adaptation? Our analysis also indicated that many mutations that appear in directed evolution variants with no obvious role in the new function exert stabilizing effects that may compensate for the destabilizing effects of the crucial function-altering mutations. Thus, the evolution of new enzymatic activities, both in nature and in the laboratory, is dependent on the compensatory, stabilizing effect of apparently "silent" mutations in regions of the protein that are irrelevant to its function.     

The Impact of Neutral Mutations on Genome Evolvability

Beneficial mutations are rare and deleterious mutations are purged by natural selection. As a result, the vast majority of mutations that accumulate in genomes belong to the class of neutral mutations. Over the last two decades, neutral mutations, despite their null effect on fitness, have been shown to affect evolvability by providing access to new phenotypes through subsequent mutations that would not have been available otherwise. Here we propose that in addition, many mutations — independent of their selective effects — can affect the mutability of neighboring DNA sequences and modulate the efficacy of homologous recombination. Such mutations do not change the spectrum of accessible phenotypes, but rather the rate at which new phenotypes will be produced. Therefore, neutral mutations that accumulate in genomes have an important long-term impact on the evolutionary fate of genomes.     

In vivo analysis of a gain-of-function mutation in the Drosophila eag-encoded K+ channel

Neuronal Na+ and K+ channels elicit currents in opposing directions and thus have opposing effects on neuronal excitability. Mutations in genes encoding Na+ or K+ channels often interact genetically, leading to either phenotypic suppression or enhancement for genes with opposing or similar effects on excitability, respectively. For example, the effects of mutations in Shaker (Sh), which encodes a K+ channel subunit, are suppressed by loss-of-function mutations in the Na+ channel structural gene para, but enhanced by loss-of-function mutations in a second K+ channel encoded by eag. Here we identify two novel mutations that suppress the effects of a Sh mutation on behavior and neuronal excitability. We used recombination mapping to localize both mutations to the eag locus, and we used sequence analysis to determine that both mutations are caused by a single amino acid substitution (G297E) in the S2-S3 linker of Eag. Because these novel eag mutations confer opposite phenotypes to eag loss-of-function mutations, we suggest that eag(G297E) causes an eag gain-of-function phenotype. We hypothesize that the G297E substitution may cause premature, prolonged, or constitutive opening of the Eag channels by favoring the "unlocked" state of the channel.     

Mutations of the CD40 ligand gene in 13 Japanese patients with X-linked hyper-IgM syndrome

X-linked hyper-IgM syndrome (XHIM) is a rare primary immunodeficiency caused by a defective CD40 ligand. We identified mutations of the CD40 ligand gene in 13 unrelated Japanese XHIM patients. Of the four patients with missense mutations, one had a mutation within the transmembrane domain, and the three others had mutations affecting the TNF homology region of the extracellular domain. Two of the missense mutations resulted in the substitution of amino acids that are highly conserved in TNF family proteins. Three patients had nonsense mutations, all of which resulted in the truncation of the TNF homology domain of the CD40 ligand. Three patients had genomic DNA deletions of 2, 3 or 4 nucleotides, respectively. All of the deletions were flanked by direct repeat sequences, suggesting that these deletions were caused by slipped mispairing. Three patients had mutations within introns resulting in altered splicing, and multiple splicing products were found in one patient. Thus, each of the 13 Japanese patients had different mutations, 9 of them being novel mutations. These results indicate that mutations in XHIM are highly heterogeneous, although codon 140 seems to be a hot spot of the CD40 ligand gene since two additional point mutations were located at Trp 140, bringing the total numbers of mutations affecting codon 140 to six. In one XHIM family with a missense mutation, prenatal diagnosis was performed by single-strand conformation polymorphism analysis of genomic DNA of a male fetus. Received: 20 August 1996     

Adaptive evolution of asexual populations under Muller's ratchet

We study the population genetics of adaptation in nonequilibrium haploid asexual populations. We find that the accumulation of deleterious mutations, due to the operation of Muller's ratchet, can considerably reduce the rate of fixation of advantageous alleles. Such reduction can be approximated reasonably well by a reduction in the effective population size. In the absence of Muller's ratchet, a beneficial mutation can only become fixed if it creates the best possible genotype; if Muller's ratchet operates, however, mutations initially arising in a nonoptimal genotype can also become fixed in the population, since the loss of the least-loaded class implies that an initially nonoptimal background can become optimal. We show that, while the rate at which adaptive mutations become fixed is reduced, the rate of fixation of deleterious mutations due to the ratchet is not changed by the presence of beneficial mutations as long as the rate of their occurrence is low and the deleterious effects of mutations (s(d)) are higher than the beneficial effects (s(a)). When s(a) > s(d), the advantage of a beneficial mutation can outweigh the deleterious effects of associated mutations. Under these conditions, a beneficial allele can drag to fixation deleterious mutations initially associated with it at a higher rate than in the absence of advantageous alleles. We propose analytical approximations for the rates of accumulation of deleterious and beneficial mutations. Furthermore, when allowing for the possible occurrence of interference between beneficial alleles, we find that the presence of deleterious mutations of either very weak or very strong effect can marginally increase the rate of accumulation of beneficial mutations over that observed in the absence of such deleterious mutations.     

Epistasis between beneficial mutations and the phenotype-to-fitness Map for a ssDNA virus

Epistatic interactions between genes and individual mutations are major determinants of the evolutionary properties of genetic systems and have therefore been well documented, but few quantitative data exist on epistatic interactions between beneficial mutations, presumably because such mutations are so much rarer than deleterious ones. We explored epistasis for beneficial mutations by constructing genotypes with pairs of mutations that had been previously identified as beneficial to the ssDNA bacteriophage ID11 and by measuring the effects of these mutations alone and in combination. We constructed 18 of the 36 possible double mutants for the nine available beneficial mutations. We found that epistatic interactions between beneficial mutations were all antagonistic-the effects of the double mutations were less than the sums of the effects of their component single mutations. We found a number of cases of decompensatory interactions, an extreme form of antagonistic epistasis in which the second mutation is actually deleterious in the presence of the first. In the vast majority of cases, recombination uniting two beneficial mutations into the same genome would not be favored by selection, as the recombinant could not outcompete its constituent single mutations. In an attempt to understand these results, we developed a simple model in which the phenotypic effects of mutations are completely additive and epistatic interactions arise as a result of the form of the phenotype-to-fitness mapping. We found that a model with an intermediate phenotypic optimum and additive phenotypic effects provided a good explanation for our data and the observed patterns of epistatic interactions.     

The Environment Affects Epistatic Interactions to Alter the Topology of an Empirical Fitness Landscape

The fitness effect of mutations can be influenced by their interactions with the environment, other mutations, or both. Previously, we constructed 32 ( = 2⁵) genotypes that comprise all possible combinations of the first five beneficial mutations to fix in a laboratory-evolved population of Escherichia coli. We found that (i) all five mutations were beneficial for the background on which they occurred; (ii) interactions between mutations drove a diminishing returns type epistasis, whereby epistasis became increasingly antagonistic as the expected fitness of a genotype increased; and (iii) the adaptive landscape revealed by the mutation combinations was smooth, having a single global fitness peak. Here we examine how the environment influences epistasis by determining the interactions between the same mutations in two alternative environments, selected from among 1,920 screened environments, that produced the largest increase or decrease in fitness of the most derived genotype. Some general features of the interactions were consistent: mutations tended to remain beneficial and the overall pattern of epistasis was of diminishing returns. Other features depended on the environment; in particular, several mutations were deleterious when added to specific genotypes, indicating the presence of antagonistic interactions that were absent in the original selection environment. Antagonism was not caused by consistent pleiotropic effects of individual mutations but rather by changing interactions between mutations. Our results demonstrate that understanding adaptation in changing environments will require consideration of the combined effect of epistasis and pleiotropy across environments.     

Mutational effects on the clonal interference phenomenon

We study the process of fixation of beneficial mutations in an asexual population by means of a theoretical model. Particularly, we wish to investigate how the supply of deleterious and beneficial mutations influences the dynamics of the adaptive process of an evolving population. It is well known that the deleterious mutations drastically affect the fate of beneficial mutations. In addition, an increasing supply of favorable mutations, to compensate the decay of the fitness due to the accumulation of deleterious mutations, produces the clonal interference phenomenon where advantageous mutations in distinct lineages compete to reach fixation. This competition imposes a limit to the speed of adaptation of the population. Intuitively, we would expect that the interplay of the two mechanisms would conspire to ensure fixation of only large-effect beneficial mutations. Our results, however, show that beneficial mutations of small effect have an increased probability of fixation when both beneficial and deleterious mutations rates are increased.     

Dynamic mutation-selection balance as an evolutionary attractor

The vast majority of mutations are deleterious and are eliminated by purifying selection. Yet in finite asexual populations, purifying selection cannot completely prevent the accumulation of deleterious mutations due to Muller's ratchet: once lost by stochastic drift, the most-fit class of genotypes is lost forever. If deleterious mutations are weakly selected, Muller's ratchet can lead to a rapid degradation of population fitness. Evidently, the long-term stability of an asexual population requires an influx of beneficial mutations that continuously compensate for the accumulation of the weakly deleterious ones. Hence any stable evolutionary state of a population in a static environment must involve a dynamic mutation-selection balance, where accumulation of deleterious mutations is on average offset by the influx of beneficial mutations. We argue that such a state can exist for any population size N and mutation rate U and calculate the fraction of beneficial mutations, ε, that maintains the balanced state. We find that a surprisingly low ε suffices to achieve stability, even in small populations in the face of high mutation rates and weak selection, maintaining a well-adapted population in spite of Muller's ratchet. This may explain the maintenance of mitochondria and other asexual genomes.     

Spectrum of small mutations in the dystrophin coding region.

Duchenne and Becker muscular dystrophies (DMD and BMD) are caused by defects in the dystrophin gene. About two-thirds of the affected patients have large deletions or duplications, which occur in the 5' and central portion of the gene. The nondeletion/duplication cases are most likely the result of smaller mutations that cannot be identified by current diagnostic screening strategies. We screened approximately 80% of the dystrophin coding sequence for small mutations in 158 patients without deletions or duplications and identified 29 mutations. The study indicates that many of the DMD and the majority of the BMD small mutations lie in noncoding regions of the gene. All of the mutations identified were unique to single patients, and most of the mutations resulted in protein truncation. We did not find a clustering of small mutations similar to the deletion distribution but found > 40% of the small mutations 3' of exon 55. The extent of protein truncation caused by the 3' mutations did not determine the phenotype, since even the exon 76 nonsense mutation resulted in the severe DMD phenotype. Our study confirms that the dystrophin gene is subject to a high rate of mutation in CpG sequences. As a consequence of not finding any hotspots or prevalent small mutations, we conclude that it is presently not possible to perform direct carrier and prenatal diagnostics for many families without deletions or duplications.     

Fitness Is Strongly Influenced by Rare Mutations of Large Effect in a Microbial Mutation Accumulation Experiment

Our understanding of the evolutionary consequences of mutation relies heavily on estimates of the rate and fitness effect of spontaneous mutations generated by mutation accumulation (MA) experiments. We performed a classic MA experiment in which frequent sampling of MA lines was combined with whole genome resequencing to develop a high-resolution picture of the effect of spontaneous mutations in a hypermutator (ΔmutS) strain of the bacterium Pseudomonas aeruginosa. After ∼644 generations of mutation accumulation, MA lines had accumulated an average of 118 mutations, and we found that average fitness across all lines decayed linearly over time. Detailed analyses of the dynamics of fitness change in individual lines revealed that a large fraction of the total decay in fitness (42.3%) was attributable to the fixation of rare, highly deleterious mutations (comprising only 0.5% of fixed mutations). Furthermore, we found that at least 0.64% of mutations were beneficial and probably fixed due to positive selection. The majority of mutations that fixed (82.4%) were base substitutions and we failed to find any signatures of selection on nonsynonymous or intergenic mutations. Short indels made up a much smaller fraction of the mutations that were fixed (17.4%), but we found evidence of strong selection against indels that caused frameshift mutations in coding regions. These results help to quantify the amount of natural selection present in microbial MA experiments and demonstrate that changes in fitness are strongly influenced by rare mutations of large effect.