Crystal structure of an ancient protein: evolution by conformational epistasis

Cellular effects and epistasis among three determinants of adaptation in experimental populations of Saccharomyces cerevisiae

Epistatic interactions in which the phenotypic effect of an allele is conditional on its genetic background have been shown to play a central part in various evolutionary processes. In a previous study (J. B. Anderson et al., Curr. Biol. 20:1383-1388, 2010; J. R. Dettman, C. Sirjusingh, L. M. Kohn, and J. B. Anderson, Nature 447:585-588, 2007), beginning with a common ancestor, we identified three determinants of fitness as mutant alleles (each designated with the letter "e") that arose in replicate Saccharomyces cerevisiae populations propagated in two different environments, a low-glucose and a high-salt environment. In a low-glucose environment, MDS3e and MKT1e interacted positively to confer a fitness advantage. Also, PMA1e from a high-salt environment interacted negatively with MKT1e in a low-glucose environment, an example of a Dobzhansky-Muller incompatibility that confers reproductive isolation. Here we showed that the negative interaction between PMA1e and MKT1e is mediated by alterations in intracellular pH, while the positive interaction between MDS3e and MKT1e is mediated by changes in gene expression affecting glucose transporter genes. We specifically addressed the evolutionary significance of the positive interaction by showing that the presence of the MDS3 mutation is a necessary condition for the spread and fixation of the new mutations at the identical site in MKT1. The expected mutations in MKT1 rose to high frequencies in two of three experimental populations carrying MDS3e but not in any of three populations carrying the ancestral allele. These data show how positive and negative epistasis can contribute to adaptation and reproductive isolation.

Determinants of homodimerization specificity in histidine kinases

Two-component signal transduction pathways consisting of a histidine kinase and a response regulator are used by prokaryotes to respond to diverse environmental and intracellular stimuli. Most species encode numerous paralogous histidine kinases that exhibit significant structural similarity. Yet in almost all known examples, histidine kinases are thought to function as homodimers. We investigated the molecular basis of dimerization specificity, focusing on the model histidine kinase EnvZ and RstB, its closest paralog in Escherichia coli. Direct binding studies showed that the cytoplasmic domains of these proteins each form specific homodimers in vitro. Using a series of chimeric proteins, we identified specificity determinants at the base of the four-helix bundle in the dimerization and histidine phosphotransfer domain. Guided by molecular coevolution predictions and EnvZ structural information, we identified sets of residues in this region that are sufficient to establish homospecificity. Mutating these residues in EnvZ to the corresponding residues in RstB produced a functional kinase that preferentially homodimerized over interacting with EnvZ. EnvZ and RstB likely diverged following gene duplication to yield two homodimers that cannot heterodimerize, and the mutants we identified represent possible evolutionary intermediates in this process.

A computational-experimental approach identifies mutations that enhance surface expression of an oseltamivir-resistant influenza neuraminidase

The His274→Tyr (H274Y) oseltamivir (Tamiflu) resistance mutation causes a substantial decrease in the total levels of surface-expressed neuraminidase protein and activity in early isolates of human seasonal H1N1 influenza, and in the swine-origin pandemic H1N1. In seasonal H1N1, H274Y only became widespread after the occurrence of secondary mutations that counteracted this decrease. H274Y is currently rare in pandemic H1N1, and it remains unclear whether secondary mutations exist that might similarly counteract the decreased neuraminidase surface expression associated with this resistance mutation in pandemic H1N1. Here we investigate the possibility of predicting such secondary mutations. We first test the ability of several computational approaches to retrospectively identify the secondary mutations that enhanced levels of surface-expressed neuraminidase protein and activity in seasonal H1N1 shortly before the emergence of oseltamivir resistance. We then use the most successful computational approach to predict a set of candidate secondary mutations to the pandemic H1N1 neuraminidase. We experimentally screen these mutations, and find that several of them do indeed partially counteract the decrease in neuraminidase surface expression caused by H274Y. Two of the secondary mutations together restore surface-expressed neuraminidase activity to wildtype levels, and also eliminate the very slight decrease in viral growth in tissue-culture caused by H274Y. Our work therefore demonstrates a combined computational-experimental approach for identifying mutations that enhance neuraminidase surface expression, and describes several specific mutations with the potential to be of relevance to the spread of oseltamivir resistance in pandemic H1N1.

The reconstructed ancestral subunit a functions as both V-ATPase isoforms Vph1p and Stv1p in Saccharomyces cerevisiae

The vacuolar ATPase complex in yeast contains two isoforms of subunit a that dictate the subcellular localization of the V-ATPase enzyme. The most recent common ancestor of subunit a (Anc.a) is reconstructed, and its function and localization in modern Saccharomyces cerevisiae are characterized. Anc.a is able to replace both subunit a isoforms.

The vacuolar-type, proton-translocating ATPase (V-ATPase) is a multisubunit enzyme responsible for organelle acidification in eukaryotic cells. Many organisms have evolved V-ATPase subunit isoforms that allow for increased specialization of this critical enzyme. Differential targeting of the V-ATPase to specific subcellular organelles occurs in eukaryotes from humans to budding yeast. In Saccharomyces cerevisiae, the two subunit a isoforms are the only difference between the two V-ATPase populations. Incorporation of Vph1p or Stv1p into the V-ATPase dictates the localization of the V-ATPase to the vacuole or late Golgi/endosome, respectively. A duplication event within fungi gave rise to two subunit a genes. We used ancestral gene reconstruction to generate the most recent common ancestor of Vph1p and Stv1p (Anc.a) and tested its function in yeast. Anc.a localized to both the Golgi/endosomal network and vacuolar membrane and acidified these compartments as part of a hybrid V-ATPase complex. Trafficking of Anc.a did not require retrograde transport from the late endosome to the Golgi that has evolved for retrieval of the Stv1p isoform. Rather, Anc.a localized to both structures through slowed anterograde transport en route to the vacuole. Our results suggest an evolutionary model that describes the differential localization of the two yeast V-ATPase isoforms.

Molecular mechanisms of epistasis within and between genes

'Disease-causing' mutations do not cause disease in all individuals. One possible important reason for this is that the outcome of a mutation can depend upon other genetic variants in a genome. These epistatic interactions between mutations occur both within and between molecules, and studies in model organisms show that they are extremely prevalent. However, epistatic interactions are still poorly understood at the molecular level, and consequently difficult to predict de novo. Here I provide an overview of our current understanding of the molecular mechanisms that can cause epistasis, and areas where more research is needed. A more complete understanding of epistasis will be vital for making accurate predictions about the phenotypes of individuals.

Mechanisms for the evolution of a derived function in the ancestral glucocorticoid receptor

Understanding the genetic, structural, and biophysical mechanisms that caused protein functions to evolve is a central goal of molecular evolutionary studies. Ancestral sequence reconstruction (ASR) offers an experimental approach to these questions. Here we use ASR to shed light on the earliest functions and evolution of the glucocorticoid receptor (GR), a steroid-activated transcription factor that plays a key role in the regulation of vertebrate physiology. Prior work showed that GR and its paralog, the mineralocorticoid receptor (MR), duplicated from a common ancestor roughly 450 million years ago; the ancestral functions were largely conserved in the MR lineage, but the functions of GRs-reduced sensitivity to all hormones and increased selectivity for glucocorticoids-are derived. Although the mechanisms for the evolution of glucocorticoid specificity have been identified, how reduced sensitivity evolved has not yet been studied. Here we report on the reconstruction of the deepest ancestor in the GR lineage (AncGR1) and demonstrate that GR's reduced sensitivity evolved before the acquisition of restricted hormone specificity, shortly after the GR-MR split. Using site-directed mutagenesis, X-ray crystallography, and computational analyses of protein stability to recapitulate and determine the effects of historical mutations, we show that AncGR1's reduced ligand sensitivity evolved primarily due to three key substitutions. Two large-effect mutations weakened hydrogen bonds and van der Waals interactions within the ancestral protein, reducing its stability. The degenerative effect of these two mutations is extremely strong, but a third permissive substitution, which has no apparent effect on function in the ancestral background and is likely to have occurred first, buffered the effects of the destabilizing mutations. Taken together, our results highlight the potentially creative role of substitutions that partially degrade protein structure and function and reinforce the importance of permissive mutations in protein evolution.

Reciprocal sign epistasis between frequently experimentally evolved adaptive mutations causes a rugged fitness landscape

The fitness landscape captures the relationship between genotype and evolutionary fitness and is a pervasive metaphor used to describe the possible evolutionary trajectories of adaptation. However, little is known about the actual shape of fitness landscapes, including whether valleys of low fitness create local fitness optima, acting as barriers to adaptive change. Here we provide evidence of a rugged molecular fitness landscape arising during an evolution experiment in an asexual population of Saccharomyces cerevisiae. We identify the mutations that arose during the evolution using whole-genome sequencing and use competitive fitness assays to describe the mutations individually responsible for adaptation. In addition, we find that a fitness valley between two adaptive mutations in the genes MTH1 and HXT6/HXT7 is caused by reciprocal sign epistasis, where the fitness cost of the double mutant prohibits the two mutations from being selected in the same genetic background. The constraint enforced by reciprocal sign epistasis causes the mutations to remain mutually exclusive during the experiment, even though adaptive mutations in these two genes occur several times in independent lineages during the experiment. Our results show that epistasis plays a key role during adaptation and that inter-genic interactions can act as barriers between adaptive solutions. These results also provide a new interpretation on the classic Dobzhansky-Muller model of reproductive isolation and display some surprising parallels with mutations in genes often associated with tumors.

Convergence, adaptation, and constraint

Convergent evolution of similar phenotypic features in similar environmental contexts has long been taken as evidence of adaptation. Nonetheless, recent conceptual and empirical developments in many fields have led to a proliferation of ideas about the relationship between convergence and adaptation. Despite criticism from some systematically minded biologists, I reaffirm that convergence in taxa occupying similar selective environments often is the result of natural selection. However, convergent evolution of a trait in a particular environment can occur for reasons other than selection on that trait in that environment, and species can respond to similar selective pressures by evolving nonconvergent adaptations. For these reasons, studies of convergence should be coupled with other methods-such as direct measurements of selection or investigations of the functional correlates of trait evolution-to test hypotheses of adaptation. The independent acquisition of similar phenotypes by the same genetic or developmental pathway has been suggested as evidence of constraints on adaptation, a view widely repeated as genomic studies have documented phenotypic convergence resulting from change in the same genes, sometimes even by the same mutation. Contrary to some claims, convergence by changes in the same genes is not necessarily evidence of constraint, but rather suggests hypotheses that can test the relative roles of constraint and selection in directing phenotypic evolution.

Utilizing natural diversity to evolve protein function: applications towards thermostability

Protein evolution relies on designing a library of sequences that capture meaningful functional diversity in a limited number of protein variants. Several approaches take advantage of the sequence space already explored through natural selection by incorporating sequence diversity available from modern genomes (and their ancestors) when designing these libraries. The success of these approaches is, partly, owing to the fact that modern sequence diversity has already been subjected to evolutionary selective forces and thus the diversity has already been deemed 'fit to survive'. Five of these approaches will be discussed in this review to highlight how protein engineers can use evolutionary sequence history/diversity of homologous proteins in unique ways to design protein libraries.

Fitness ranking of individual mutants drives patterns of epistatic interactions in HIV-1

Fitness interactions between mutations, referred to as epistasis, can strongly impact evolution. For RNA viruses and retroviruses with their high mutation rates, epistasis may be particularly important to overcome fitness losses due to the accumulation of deleterious mutations and thus could influence the frequency of mutants in a viral population. As human immunodeficiency virus type 1 (HIV-1) resistance to azidothymidine (AZT) requires selection of sequential mutations, it is a good system to study the impact of epistasis. Here we present a thorough analysis of a classical AZT-resistance pathway (the 41–215 cluster) of HIV-1 variants by fitness measurements in single round infection assays covering physiological drug concentrations ex vivo. The sign and value of epistasis varied and did not predict the epistatic effect on the mutant frequency. This complex behavior is explained by the fitness ranking of the variants that strongly depends on environmental factors, i.e., the presence and absence of drugs and the host cells used. Although some interactions compensate fitness losses, the observed small effect on the relative mutant frequencies suggests that epistasis might be inefficient as a buffering mechanism for fitness losses in vivo. While the use of epistasis-based hypotheses to make general assumptions on the evolutionary dynamics of viral populations is appealing, our data caution their interpretation without further knowledge on the characteristics of the viral mutant spectrum under different environmental conditions.

Evolution of promiscuous nuclear hormone receptors: LXR, FXR, VDR, PXR, and CAR

Nuclear hormone receptors (NHRs) are transcription factors that work in concert with co-activators and co-repressors to regulate gene expression. Some examples of ligands for NHRs include endogenous compounds such as bile acids, retinoids, steroid hormones, thyroid hormone, and vitamin D. This review describes the evolution of liver X receptors α and β (NR1H3 and 1H2, respectively), farnesoid X receptor (NR1H4), vitamin D receptor (NR1I1), pregnane X receptor (NR1I2), and constitutive androstane receptor (NR1I3). These NHRs participate in complex, overlapping transcriptional regulation networks involving cholesterol homeostasis and energy metabolism. Some of these receptors, particularly PXR and CAR, are promiscuous with respect to the structurally wide range of ligands that act as agonists. A combination of functional and computational analyses has shed light on the evolutionary changes of NR1H and NR1I receptors across vertebrates, and how these receptors may have diverged from ancestral receptors that first appeared in invertebrates.

Nuclear hormone receptors in nematodes: evolution and function

Nuclear hormone receptors (NHRs) are proteins that regulate gene expression in response to developmental, environmental, and nutritional signals. The activity of some NHRs is selectively and reversibly modulated by small molecular weight compounds. However, for others - termed "orphan" receptors - no such ligands have (yet) been identified, and at least some NHRs may lack natural ligands. NHRs exhibit a stereotyped architecture, with conserved N-terminal DNA-binding domains (DBDs) and more variable C-terminal ligand-binding domains (LBDs). NHRs control the transcription of remarkably diverse and specific gene networks, apparently by integrating multiple regulatory inputs that interact with distinct receptor surfaces; these inputs include small molecule ligands, transcriptional coregulators, and response elements, the genomic sites to which the receptors bind. NHRs comprise an ancient superfamily found in all metazoans, and recent findings have revealed NHR-like regulatory factors in fungi. Here, we consider NHR function and evolution in nematodes, roundworms that inhabit terrestrial, marine, and freshwater habitats; we focus in particular on the well-established experimental organism Caenorhabditis elegans. Interestingly, the C. elegans genome encodes a massively expanded NHR family; we speculate that some of the multiple physiological activities governed by individual mammalian NHRs may be distributed among multiple members of the C. elegans family, potentially focusing and simplifying functional analyses. Accordingly, investigations of relevant NHR cofactors, ligands, and response elements might also prove to be simpler; moreover, the abbreviated intergenic regions of the C. elegans genome will facilitate the assignment of response elements to target genes. Finally, the growing interest in medically relevant nematodes is providing novel insights into the function and evolution of NHRs.

What are nuclear receptor ligands?

Nuclear receptors (NRs) are a family of highly conserved transcription factors that regulate transcription in response to small lipophilic compounds. They play a role in every aspect of development, physiology and disease in humans. They are also ubiquitous in and unique to the animal kingdom suggesting that they may have played an important role in their evolution. In contrast to the classical endocrine receptors that originally defined the family, recent studies suggest that the first NRs might have been sensors of their environment, binding ligands that were external to the host organism. The purpose of this review is to provide a broad perspective on NR ligands and address the issue of exactly what constitutes a NR ligand from historical, biological and evolutionary perspectives. This discussion will lay the foundation for subsequent reviews in this issue as well as pose new questions for future investigation.

Correlated evolution of nearby residues in Drosophilid proteins

Here we investigate the correlations between coding sequence substitutions as a function of their separation along the protein sequence. We consider both substitutions between the reference genomes of several Drosophilids as well as polymorphisms in a population sample of Zimbabwean Drosophila melanogaster. We find that amino acid substitutions are “clustered” along the protein sequence, that is, the frequency of additional substitutions is strongly enhanced within ≈10 residues of a first such substitution. No such clustering is observed for synonymous substitutions, supporting a “correlation length” associated with selection on proteins as the causative mechanism. Clustering is stronger between substitutions that arose in the same lineage than it is between substitutions that arose in different lineages. We consider several possible origins of clustering, concluding that epistasis (interactions between amino acids within a protein that affect function) and positional heterogeneity in the strength of purifying selection are primarily responsible. The role of epistasis is directly supported by the tendency of nearby substitutions that arose on the same lineage to preserve the total charge of the residues within the correlation length and by the preferential cosegregation of neighboring derived alleles in our population sample. We interpret the observed length scale of clustering as a statistical reflection of the functional locality (or modularity) of proteins: amino acids that are near each other on the protein backbone are more likely to contribute to, and collaborate toward, a common subfunction.

Genes are templates for proteins, yet evolutionary studies of genes and proteins often bear little resemblance. Analyses of gene evolution typically treat each codon independently, quantifying gene evolution by summing over the constituent codons. In contrast, studies of protein evolution generally incorporate protein structure and interactions between amino acids explicitly. We investigate correlations in the evolution of codons as a function of their distance from each other along the protein coding sequence. This approach is motivated by the expectation that codons near each other in sequence often encode amino acids belonging to the same functional unit. Consequently, these amino acids are more likely to interact and/or experience similar selective regimes, introducing correlation between the evolution of the underlying codons. We find codon evolution in Drosophilids to be correlated over a characteristic length scale of ≈10 codons. Specifically, the presence of a non-synonymous substitution substantially increases the probability of further such substitutions nearby, particularly within that lineage. Further analysis suggests both functional interactions between amino acids and correlation in the strength of selection contribute to this effect. These findings are relevant for understanding the relative importance of different modes of selection, and particularly the role of epistasis, in gene and protein evolution.

Evolution of ligand specificity in vertebrate corticosteroid receptors

Background

Corticosteroid receptors include mineralocorticoid (MR) and glucocorticoid (GR) receptors. Teleost fishes have a single MR and duplicate GRs that show variable sensitivities to mineralocorticoids and glucocorticoids. How these receptors compare functionally to tetrapod MR and GR, and the evolutionary significance of maintaining two GRs, remains unclear.

Results

We used up to seven steroids (including aldosterone, cortisol and 11-deoxycorticosterone [DOC]) to compare the ligand specificity of the ligand binding domains of corticosteroid receptors between a mammal (Mus musculus) and the midshipman fish (Porichthys notatus), a teleost model for steroid regulation of neural and behavioral plasticity. Variation in mineralocorticoid sensitivity was considered in a broader phylogenetic context by examining the aldosterone sensitivity of MR and GRs from the distantly related daffodil cichlid (Neolamprologus pulcher), another teleost model for neurobehavioral plasticity. Both teleost species had a single MR and duplicate GRs. All MRs were sensitive to DOC, consistent with the hypothesis that DOC was the initial ligand of the ancestral MR. Variation in GR steroid-specificity corresponds to nine identified amino acid residue substitutions rather than phylogenetic relationships based on receptor sequences.

Conclusion

The mineralocorticoid sensitivity of duplicate GRs in teleosts is highly labile in the context of their evolutionary phylogeny, a property that likely led to neo-functionalization and maintenance of two GRs.

Tryptophanyl-tRNA synthetase Urzyme: a model to recapitulate molecular evolution and investigate intramolecular complementation

We substantiate our preliminary description of the class I tryptophanyl-tRNA synthetase minimal catalytic domain with details of its construction, structure, and steady-state kinetic parameters. Generating that active fragment involved deleting 65% of the contemporary enzyme, including the anticodon-binding domain and connecting peptide 1, CP1, a 74-residue internal segment from within the Rossmann fold. We used protein design (Rosetta), rather than phylogenetic sequence alignments, to identify mutations to compensate for the severe loss of modularity, thus restoring stability, as evidenced by renaturation described previously and by 70-ns molecular dynamics simulations. Sufficient solubility to enable biochemical studies was achieved by expressing the redesigned Urzyme as a maltose-binding protein fusion. Michaelis-Menten kinetic parameters from amino acid activation assays showed that, compared with the native full-length enzyme, TrpRS Urzyme binds ATP with similar affinity. This suggests that neither of the two deleted structural modules has a strong influence on ground-state ATP binding. However, tryptophan has 10(3) lower affinity, and the Urzyme has comparably reduced specificity relative to the related amino acid, tyrosine. Molecular dynamics simulations revealed how CP1 may contribute significantly to cognate amino acid specificity. As class Ia editing domains are nested within the CP1, this finding suggests that this module enhanced amino acid specificity continuously, throughout their evolution. We call this type of reconstructed protein catalyst an Urzyme (Ur prefix indicates original, primitive, or earliest). It establishes a model for recapitulating very early steps in molecular evolution in which fitness may have been enhanced by accumulating entire modules, rather than by discrete amino acid sequence changes.

Systematic dissection and trajectory-scanning mutagenesis of the molecular interface that ensures specificity of two-component signaling pathways

Two-component signal transduction systems enable bacteria to sense and respond to a wide range of environmental stimuli. Sensor histidine kinases transmit signals to their cognate response regulators via phosphorylation. The faithful transmission of information through two-component pathways and the avoidance of unwanted cross-talk require exquisite specificity of histidine kinase-response regulator interactions to ensure that cells mount the appropriate response to external signals. To identify putative specificity-determining residues, we have analyzed amino acid coevolution in two-component proteins and identified a set of residues that can be used to rationally rewire a model signaling pathway, EnvZ-OmpR. To explore how a relatively small set of residues can dictate partner selectivity, we combined alanine-scanning mutagenesis with an approach we call trajectory-scanning mutagenesis, in which all mutational intermediates between the specificity residues of EnvZ and another kinase, RstB, were systematically examined for phosphotransfer specificity. The same approach was used for the response regulators OmpR and RstA. Collectively, the results begin to reveal the molecular mechanism by which a small set of amino acids enables an individual kinase to discriminate amongst a large set of highly-related response regulators and vice versa. Our results also suggest that the mutational trajectories taken by two-component signaling proteins following gene or pathway duplication may be constrained and subject to differential selective pressures. Only some trajectories allow both the maintenance of phosphotransfer and the avoidance of unwanted cross-talk.

Maintaining the specificity of signal transduction pathways is critical to the ability of cells to process information, make decisions, and regulate their behavior. Preventing cross-talk often relies predominantly on molecular recognition and a set of specificity-determining residues in cognate proteins. Identifying these residues and understanding how they dictate specificity is still a major challenge. Additionally, we have a rudimentary understanding of how specificity evolves, particularly after gene duplication events. We tackled these questions using two-component signaling proteins, the largest family of bacterial signaling proteins. Using analyses of amino acid coevolution, we pinpointed a set of specificity residues in histidine kinases and their cognate substrates. Then, using systematic mutagenesis, we characterized the complete set of intermediates between two different signaling systems, EnvZ/OmpR and RstA/RstB. The results demonstrate that specificity residues contribute unequally and, importantly, that some residues depend substantially on the identity of neighboring residues. We also demonstrate how the specificity of EnvZ/OmpR can be reprogrammed to match that of RstB/RstA through a series of individual substitutions without disrupting the kinase/regulator interaction. Notably, this property is not shared by all trajectories from EnvZ/OmpR to RstA/RstB, suggesting that the duplication/divergence process that likely produced these two pathways may have been fundamentally constrained.

Exploring the interplay of stability and function in protein evolution: new methods further elucidate why protein stability is necessarily so tenuous and stability-increasing mutations compromise biological function

A new split beta-lactamase assay promises experimental testing of the interplay of protein stability and function. Proteins are sufficiently stable to act effectively within cells. However, mutations generally destabilize structure, with effects on free energy that are comparable to the free energy of folding. Assays of protein functionality and stability in vivo enable a quick study of factors that influence these properties in response to targeted mutations. These assays can help molecular engineering but can also be used to target important questions, including why most proteins are marginally stable, how mutations alter structural makeup, and how thermodynamics, function, and environment shape molecular change. Processes of self-organization and natural selection are determinants of stability and function. Non-equilibrium thermodynamics provides crucial concepts, e.g., cells as emergent energy-dissipating entities that do work and build their own parts, and a framework to study the sculpting role of evolution at different scales.

Evolution of a signaling nexus constrained by protein interfaces and conformational States

Heterotrimeric G proteins act as the physical nexus between numerous receptors that respond to extracellular signals and proteins that drive the cytoplasmic response. The Gα subunit of the G protein, in particular, is highly constrained due to its many interactions with proteins that control or react to its conformational state. Various organisms contain differing sets of Gα-interacting proteins, clearly indicating that shifts in sequence and associated Gα functionality were acquired over time. These numerous interactions constrained much of Gα evolution; yet Gα has diversified, through poorly understood processes, into several functionally specialized classes, each with a unique set of interacting proteins. Applying a synthetic sequence-based approach to mammalian Gα subunits, we established a set of seventy-five evolutionarily important class-distinctive residues, sites where a single Gα class is differentiated from the three other classes. We tested the hypothesis that shifts at these sites are important for class-specific functionality. Importantly, we mapped known and well-studied class-specific functionalities from all four mammalian classes to sixteen of our class-distinctive sites, validating the hypothesis. Our results show how unique functionality can evolve through the recruitment of residues that were ancestrally functional. We also studied acquisition of functionalities by following these evolutionarily important sites in non-mammalian organisms. Our results suggest that many class-distinctive sites were established early on in eukaryotic diversification and were critical for the establishment of new Gα classes, whereas others arose in punctuated bursts throughout metazoan evolution. These Gα class-distinctive residues are rational targets for future structural and functional studies.

Protein evolution by molecular tinkering: diversification of the nuclear receptor superfamily from a ligand-dependent ancestor

Phylogenetic reconstruction of the structure and function of the ancestor of the nuclear receptor protein family reveals how functional diversity evolves by subtle tinkering with an ancestral template.

Understanding how protein structures and functions have diversified is a central goal in molecular evolution. Surveys of very divergent proteins from model organisms, however, are often insufficient to determine the features of ancestral proteins and to reveal the evolutionary events that yielded extant diversity. Here we combine genomic, biochemical, functional, structural, and phylogenetic analyses to reconstruct the early evolution of nuclear receptors (NRs), a diverse superfamily of transcriptional regulators that play key roles in animal development, physiology, and reproduction. By inferring the structure and functions of the ancestral NR, we show—contrary to current belief—that NRs evolved from a ligand-activated ancestral receptor that existed near the base of the Metazoa, with fatty acids as possible ancestral ligands. Evolutionary tinkering with this ancestral structure generated the extraordinary diversity of modern receptors: sensitivity to different ligands evolved because of subtle modifications of the internal cavity, and ligand-independent activation evolved repeatedly because of various mutations that stabilized the active conformation in the absence of ligand. Our findings illustrate how a mechanistic dissection of protein evolution in a phylogenetic context can reveal the deep homology that links apparently “novel” molecular functions to a common ancestral form.

Many protein families are so diverse that it is hard to determine their ancestral functions and to understand how their derived functions evolved. The existence of so many different functions within protein families often creates the impression that complex, novel functions must have evolved repeatedly and independently. Nuclear receptors (NRs) are a large family of related proteins that regulate key biological processes in animals by binding to specific DNA sequences and triggering expression of nearby target genes. Many NRs are activated by a specific hormone or other small molecule, but some do not require a ligand, and still others are incapable of activating gene expression and so act primarily as repressors of transcription. To understand how the functional diversity of NRs evolved, we reconstructed the structural and functional characteristics of the ancient protein from which the entire family evolved, using genomic, biochemical, functional, and structural analyses in a phylogenetic framework. We show, contrary to current belief, that the ancestral NR was a ligand-activated transcriptional activator that existed in the earliest period of animal evolution. Our analysis reveals how the extraordinary functional diversity of modern receptors was generated by subtle tinkering with this ancestral template—slightly reshaping the ligand cavity, stabilizing the protein's active conformation so it no longer required a ligand, or disabling the protein's capacity to activate transcription without affecting its other properties. We predict that, when sufficient data are gathered to allow detailed evolutionary reconstructions in other protein families, it will become apparent that most protein functional diversity evolved by tinkering with ancient functions; invoking the evolution of wholesale “novelty” will seldom be necessary.

Molecular and structural analysis of mosaic variants of penicillin-binding protein 2 conferring decreased susceptibility to expanded-spectrum cephalosporins in Neisseria gonorrhoeae: role of epistatic mutations

Mutations in penicillin-binding protein 2 (PBP 2) encoded by mosaic penA alleles are crucial for intermediate resistance to the expanded-spectrum cephalosporins ceftriaxone and cefixime in Neisseria gonorrhoeae. Three of the ∼60 mutations present in mosaic alleles of penA, G545S, I312M, and V316T, have been reported to be responsible for increased resistance, especially to cefixime [Takahata, S., et al. (2006) Antimicrob. Agents Chemother. 50, 3638-3645]. However, we observed that the minimum inhibitory concentrations (MICs) of penicillin, ceftriaxone, and cefixime for a wild-type strain (FA19) containing a penA gene with these three mutations increased only 1.5-, 1.5-, and 3.5-fold, respectively. In contrast, when these three mutations in a mosaic penA allele (penA35) were reverted back to the wild type and the gene was transformed into FA19, the MICs of the three antibiotics were reduced to near wild-type levels. Thus, these three mutations display epistasis, in that their capacity to increase resistance to β-lactam antibiotics is dependent on the presence of other mutations in the mosaic alleles. We also identified an additional mutation, N512Y, that contributes to the decreased susceptibility to expanded-spectrum cephalosporins. Finally, we investigated the effects of a mutation (A501V) currently found only in nonmosaic penA alleles on decreased susceptibility to ceftriaxone and cefixime, with the expectation that this mutation may arise in mosaic alleles. Transfer of the mosaic penA35 allele containing an A501V mutation to FA6140, a chromosomally mediated penicillin-resistant isolate, increased the MICs of ceftriaxone (0.4 μg/mL) and cefixime (1.2 μg/mL) to levels above their respective break points. The proposed structural mechanisms of these mutations are discussed in light of the recently published structure of PBP 2.

Permissive secondary mutations enable the evolution of influenza oseltamivir resistance

The His274-->Tyr274 (H274Y) mutation confers oseltamivir resistance on N1 influenza neuraminidase but had long been thought to compromise viral fitness. However, beginning in 2007-2008, viruses containing H274Y rapidly became predominant among human seasonal H1N1 isolates. We show that H274Y decreases the amount of neuraminidase that reaches the cell surface and that this defect can be counteracted by secondary mutations that also restore viral fitness. Two such mutations occurred in seasonal H1N1 shortly before the widespread appearance of H274Y. The evolution of oseltamivir resistance was therefore enabled by "permissive" mutations that allowed the virus to tolerate subsequent occurrences of H274Y. An understanding of this process may provide a basis for predicting the evolution of oseltamivir resistance in other influenza strains.

Hormone binding and co-regulator binding to the glucocorticoid receptor are allosterically coupled

The glucocorticoid receptor initiates the cellular response to glucocorticoid steroid hormones in vertebrates. Co-regulator proteins dock to the receptor in response to hormone binding and potentiate the transcriptional activity of the receptor by modifying DNA and recruiting essential transcription factors like RNA polymerase II. Hormones and co-regulators bind at distinct sites in the ligand binding domain yet function cooperatively to mediate transcriptional control. This study reveals and quantifies energetic coupling between two binding sites using purified components. Using a library of peptides taken from co-regulator proteins, we determine the pattern of co-regulator binding to the glucocorticoid receptor ligand binding domain. We show that peptides from co-regulators differ in their effects on hormone binding and kinetics. Peptides from DAX1 and SRC1 bind with similar affinity, but DAX1 binding is coupled to hormone binding, and SRC1 is not. Mechanistic details of co-regulator binding and coupling to the hormone binding pocket are uncovered by analysis of properties endowed by mutation of a key residue in the allosteric network connecting the sites.

Analyzing protein structure and function using ancestral gene reconstruction

Protein families with functionally diverse members can illuminate the structural determinants of protein function and the process by which protein structure and function evolve. To identify the key amino acid changes that differentiate one family member from another, most studies have taken a horizontal approach, swapping candidate residues between present-day family members. This approach has often been stymied, however, by the fact that shifts in function often require multiple interacting mutations; chimeric proteins are often nonfunctional, either because one lineage has amassed mutations that are incompatible with key residues that conferred a new function on other lineages, or because it lacks mutations required to support those key residues. These difficulties can be overcome by using a vertical strategy, which reconstructs ancestral genes and uses them as the appropriate background in which to study the effects of historical mutations on functional diversification. In this review, we discuss the advantages of the vertical strategy and highlight several exemplary studies that have used ancestral gene reconstruction to reveal the molecular underpinnings of protein structure, function, and evolution.

Robustness of ancestral sequence reconstruction to phylogenetic uncertainty

Ancestral sequence reconstruction (ASR) is widely used to formulate and test hypotheses about the sequences, functions, and structures of ancient genes. Ancestral sequences are usually inferred from an alignment of extant sequences using a maximum likelihood (ML) phylogenetic algorithm, which calculates the most likely ancestral sequence assuming a probabilistic model of sequence evolution and a specific phylogeny—typically the tree with the ML. The true phylogeny is seldom known with certainty, however. ML methods ignore this uncertainty, whereas Bayesian methods incorporate it by integrating the likelihood of each ancestral state over a distribution of possible trees. It is not known whether Bayesian approaches to phylogenetic uncertainty improve the accuracy of inferred ancestral sequences. Here, we use simulation-based experiments under both simplified and empirically derived conditions to compare the accuracy of ASR carried out using ML and Bayesian approaches. We show that incorporating phylogenetic uncertainty by integrating over topologies very rarely changes the inferred ancestral state and does not improve the accuracy of the reconstructed ancestral sequence. Ancestral state reconstructions are robust to uncertainty about the underlying tree because the conditions that produce phylogenetic uncertainty also make the ancestral state identical across plausible trees; conversely, the conditions under which different phylogenies yield different inferred ancestral states produce little or no ambiguity about the true phylogeny. Our results suggest that ML can produce accurate ASRs, even in the face of phylogenetic uncertainty. Using Bayesian integration to incorporate this uncertainty is neither necessary nor beneficial.

Cyclin-dependent kinase 5 modulates the transcriptional activity of the mineralocorticoid receptor and regulates expression of brain-derived neurotrophic factor

Glucocorticoids, major end effectors of the stress response, play an essential role in the homeostasis of the central nervous system (CNS) and contribute to memory consolidation and emotional control through their intracellular receptors, the glucocorticoid and mineralocorticoid receptors. Cyclin-dependent kinase 5 (CDK5), on the other hand, plays important roles in the morphogenesis and functions of the central nervous system, and its aberrant activation has been associated with development of neurodegenerative disorders. We previously reported that CDK5 phosphorylated the glucocorticoid receptor and modulated its transcriptional activity. Here we found that CDK5 also regulated mineralocorticoid receptor-induced transcriptional activity by phosphorylating multiple serine and threonine residues located in its N-terminal domain through physical interaction. Aldosterone and dexamethasone, respectively, increased and suppressed mRNA/protein expression of brain-derived neurotrophic factor (BDNF) in rat cortical neuronal cells, whereas the endogenous glucocorticoid corticosterone showed a biphasic effect. CDK5 enhanced the effect of aldosterone and dexamethasone on BDNF expression. Because this neurotrophic factor plays critical roles in neuronal viability, synaptic plasticity, consolidation of memory, and emotional changes, we suggest that aberrant activation of CDK5 might influence these functions through corticosteroid receptors/BDNF.

Two frequenins in Drosophila: unveiling the evolutionary history of an unusual neuronal calcium sensor (NCS) duplication

BACKGROUND

Drosophila frequenin (Frq), the homolog of the mammalian neuronal calcium sensor-1 (NCS-1), is a high affinity calcium-binding protein with ubiquitous expression in the nervous system. This protein has an important role in the regulation of neurotransmitter release per synapse, axonal growth and bouton formation. In D. melanogaster, frequenin is encoded by two genes (frq1 and frq2), a very unexpected feature in the Frq/NCS-1 subfamily. These genes are located in tandem in the same genomic region, and their products are 95% identical in their amino acid sequence, clearly indicating their recent origin by gene duplication. Here, we have investigated the factors involved in this unusual feature by examining the molecular evolution of the two frq genes in Drosophila and the evolutionary dynamics of NCS family in a large set of bilaterian species.

RESULTS

Surprisingly, we have found no amino acid replacements fixed across the twelve Drosophila species surveyed. In contrast, synonymous substitutions have been prevalent in the evolution of the coding region of frq1 and frq2, indicating the presence of strong functional constraints following gene duplication. Despite that, we have detected that significant evolutionary rate acceleration had occurred in Frq1 in early times from the duplication, in which positive selection (likely promoting functional diversification) had probably an important role. The analysis of sequence conservation and DNA topology at the non-coding regions of both genes has allowed the identification of DNA regions candidates to be cis-regulatory elements. The results reveal a possible mechanism of regulatory diversification between frq1 and frq2.

CONCLUSIONS

The presence of two frequenins in Drosophila and the rapid accumulation of amino acid substitutions after gene duplication are very unusual features in the evolution of the Frq/NCS-1 subfamily. Here we show that the action of positive selection in concordance with some extent of regulatory diversification might explain these findings. Selected amino acid substitutions in Frq1 likely contributed to the functional divergence between the two duplicates, which, in turn, should have diverged in their regulation by ecdysone-induced early genes.

Fitness epistasis and constraints on adaptation in a human immunodeficiency virus type 1 protein region

Fitness epistasis, the interaction among alleles at different loci in their effects on fitness, has potentially important consequences for adaptive evolution. We investigated fitness epistasis among amino acids of a functionally important region of the human immunodeficiency virus type 1 (HIV-1) exterior envelope glycoprotein (gp120). Seven mutations putatively involved in the adaptation of the second conserved to third variable protein region (C2-V3) to the use of an alternative host-cell chemokine coreceptor (CXCR4) for cell entry were engineered singly and in combinations on the wild-type genetic background and their effects on viral infectivity were measured. Epistasis was found to be common and complex, involving not only pairwise interactions, but also higher-order interactions. Interactions could also be surprisingly strong, changing fitness by more than 9 orders of magnitude, which is explained by some single mutations being practically lethal. A consequence of the observed epistasis is that many of the minimum-length mutational trajectories between the wild type and the mutant with highest fitness on cells expressing the alternative coreceptor are selectively inaccessible. These results may help explain the difficulty of evolving viruses that use the alternative coreceptor in culture and the delayed evolution of this phenotype in natural infection. Knowledge of common, complex, and strong fitness interactions among amino acids is necessary for a full understanding of protein evolution.

Retracing evolution of red fluorescence in GFP-like proteins from Faviina corals

Proteins of the green fluorescent protein family represent a convenient experimental model to study evolution of novelty at the molecular level. Here, we focus on the origin of Kaede-like red fluorescent proteins characteristic of the corals of the Faviina suborder. We demonstrate, using an original approach involving resurrection and analysis of the library of possible evolutionary intermediates, that it takes on the order of 12 mutations, some of which strongly interact epistatically, to fully recapitulate the evolution of a red fluorescent phenotype from the ancestral green. Five of the identified mutations would not have been found without the help of ancestral reconstruction, because the corresponding site states are shared between extant red and green proteins due to their recent descent from a dual-function common ancestor. Seven of the 12 mutations affect residues that are not in close contact with the chromophore and thus must exert their effect indirectly through adjustments of the overall protein fold; the relevance of these mutations could not have been anticipated from the purely theoretical analysis of the protein's structure. Our results introduce a powerful experimental approach for comparative analysis of functional specificity in protein families even in the cases of pronounced epistasis, provide foundation for the detailed studies of evolutionary trajectories leading to novelty and complexity, and will help rational modification of existing fluorescent labels.

Identification of a novel set of scaffolding residues that are instrumental for the inhibitory property of Kunitz (STI) inhibitors

For canonical serine protease inhibitors (SPIs), scaffolding spacer residue Asn or Arg religates cleaved scissile peptide bond to offer efficient inhibition. However, several designed "mini-proteins," containing the inhibitory loop and the spacer(s) with trimmed scaffold behave like substrates, indicating that scaffolding region beyond the spacer is also important in the inhibitory process. To understand the loop-scaffold compatibility, we prepared three chimeric proteins ECI(L)-WCI(S), ETI(L)-WCI(S), and STI(L)-WCI(S), where the inhibitory loop of ECI, ETI, and STI is placed on the scaffold of their homolog WCI. Results show that although ECI(L)-WCI(S) and STI(L)-WCI(S) behave like good inhibitors, ETI(L)-WCI(S) behaves like a substrate. That means a set of loop residues (SRLRSAFI), offering strong trypsin inhibition in ETI, act as a substrate when they seat on the scaffold of WCI. Crystal structure of ETI(L)-WCI(S) shows that the inhibitory loop is of noncanonical conformation. We identified three novel scaffolding residues Trp88, Arg74, and Tyr113 in ETI that act as barrier to confine the inhibitory loop to canonical conformation. Absence of this barrier in the scaffold of WCI makes the inhibitory loop flexible in ETI(L)-WCI(S) leading to a loss of canonical conformation, explaining its substrate-like behavior. Incorporation of this barrier back in ETI(L)-WCI(S) through mutations increases its inhibitory power, supporting our proposition. Our study provides structural evidence for the contribution of remote scaffolding residues in the inhibitory process of canonical SPIs. Additionally, we rationalize why the loop-scaffold swapping is not permitted even among the members of highly homologous inhibitors, which might be important in the light of inhibitor design.

The human glucocorticoid receptor: molecular basis of biologic function

The characterization of the subfamily of steroid hormone receptors has enhanced our understanding of how a set of hormonally derived lipophilic ligands controls cellular and molecular functions to influence development and help achieve homeostasis. The glucocorticoid receptor (GR), the first member of this subfamily, is a ubiquitously expressed intracellular protein, which functions as a ligand-dependent transcription factor that regulates the expression of glucocorticoid-responsive genes. The effector domains of the GR mediate transcriptional activation by recruiting coregulatory multi-subunit complexes that remodel chromatin, target initiation sites, and stabilize the RNA-polymerase II machinery for repeated rounds of transcription of target genes. This review summarizes the basic aspects of the structure and actions of the human (h) GR, and the molecular basis of its biologic functions.

Evolutionary expansion and specialization of the PDZ domains

PDZ domains are protein-protein interaction modules widely used to assemble membranous signaling complexes including those found in the neuronal synapse. PDZ-containing genes encoded in metazoan genomes vastly outnumber those in prokaryotes, plants, and fungi. By comparing 40 proteomes to track the evolutionary history of the PDZ domain, we observed that the variety of associations between PDZ and other domains expands greatly along the stem leading to metazoans and choanoflagellates. We asked whether the expansion of PDZ domains was due to random or specific sequence changes. Studying the sequence signatures of 58 PDZ lineages that are common to bilaterian animals, we showed that six common amino acid residues are able to classify 96% of PDZ domains to their correct evolutionary lineage. In PDZ domain-ligand cocrystals, four of these "classifying positions" lie in direct contact with the -1 and -3 residues of the ligand. This suggests coevolution of the more flexible regions of the binding interaction as a central mechanism of specialization inherent within the PDZ domain. To identify these positions, we devised two independent algorithms--a metric termed within-clade entropy (WCE) and an average mutual information (AvgMI) score--that both reached similar results. Extending these tools to the choanoflagellate, Monosiga brevicollis, we compared its PDZ domains with their putative metazoan orthologs. Interestingly, the M. brevicollis genes lack conservation at the classifying positions suggesting dissociation between domain organization in multidomain proteins and specific changes within the PDZ domain.

Stepwise modification of a modular enhancer underlies adaptation in a Drosophila population

The evolution of cis regulatory elements (enhancers) of developmentally regulated genes plays a large role in the evolution of animal morphology. However, the mutational path of enhancer evolution--the number, origin, effect, and order of mutations that alter enhancer function--has not been elucidated. Here, we localized a suite of substitutions in a modular enhancer of the ebony locus responsible for adaptive melanism in a Ugandan Drosophila population. We show that at least five mutations with varied effects arose recently from a combination of standing variation and new mutations and combined to create an allele of large phenotypic effect. We underscore how enhancers are distinct macromolecular entities, subject to fundamentally different, and generally more relaxed, functional constraints relative to protein sequences.

An epistatic ratchet constrains the direction of glucocorticoid receptor evolution

The extent to which evolution is reversible has long fascinated biologists. Most previous work on the reversibility of morphological and life-history evolution has been indecisive, because of uncertainty and bias in the methods used to infer ancestral states for such characters. Further, despite theoretical work on the factors that could contribute to irreversibility, there is little empirical evidence on its causes, because sufficient understanding of the mechanistic basis for the evolution of new or ancestral phenotypes is seldom available. By studying the reversibility of evolutionary changes in protein structure and function, these limitations can be overcome. Here we show, using the evolution of hormone specificity in the vertebrate glucocorticoid receptor as a case-study, that the evolutionary path by which this protein acquired its new function soon became inaccessible to reverse exploration. Using ancestral gene reconstruction, protein engineering and X-ray crystallography, we demonstrate that five subsequent 'restrictive' mutations, which optimized the new specificity of the glucocorticoid receptor, also destabilized elements of the protein structure that were required to support the ancestral conformation. Unless these ratchet-like epistatic substitutions are restored to their ancestral states, reversing the key function-switching mutations yields a non-functional protein. Reversing the restrictive substitutions first, however, does nothing to enhance the ancestral function. Our findings indicate that even if selection for the ancestral function were imposed, direct reversal would be extremely unlikely, suggesting an important role for historical contingency in protein evolution.

Following evolutionary paths to protein-protein interactions with high affinity and selectivity

How do intricate multi-residue features such as protein-protein interfaces evolve? To address this question, we evolved a new colicin-immunity binding interaction. We started with Im9, which inhibits its cognate DNase ColE9 at 10(-14) M affinity, and evolved it toward ColE7, which it inhibits promiscuously (Kd > 10(-8) M). Iterative rounds of random mutagenesis and selection toward higher affinity for ColE7, and selectivity (against ColE9 inhibition), led to an approximately 10(5)-fold increase in affinity and a 10(8)-fold increase in selectivity. Analysis of intermediates along the evolved variants revealed that changes in the binding configuration of the Im protein uncovered a latent set of interactions, thus providing the key to the rapid divergence of new Im7 variants. Overall, protein-protein interfaces seem to share the evolvability features of enzymes, that is, the exploitation of promiscuous interactions and alternative binding configurations via 'generalist' intermediates, and the key role of compensatory stabilizing mutations in facilitating the divergence of new functions.

Pharmacological characterization of intracellular, membrane, and plasma binding sites for corticosterone in house sparrows

The diversity and specificity of glucocorticoid effects are dependent on cell-specific receptor mechanisms. Three known corticosteroid receptors mediate tissue effects of glucocorticoids in vertebrates: two intracellular receptors that act primarily as ligand-activated transcription factors, and a membrane-associated receptor. The intracellular receptor sub-types have been well characterized in mammals, however relatively little is known about them across non-mammalian vertebrates. The membrane-associated receptors are poorly characterized in most vertebrate taxa. To explore the basis for glucocorticoid action in birds, we pharmacologically characterized the three putative corticosteroid receptors in the brain, as well as a plasma corticosterone binding globulin, in the house sparrow (Passer domesticus). We found that house sparrow brain cytosol contained two distinguishable binding sites for corticosterone. A high affinity, mineralocorticoid-like receptor had subnanomolar affinity for corticosterone (K(d) approximately 0.2 nM). However, this 'MR-like' high-affinity receptor did not bind RU28318 or canrenoic acid, two compounds that bind mammalian MR with high affinity. A lower-affinity, glucocorticoid-like receptor in brain cytosol bound corticosterone with an average K(d)=5.61 nM. This GR-like receptor showed subnanomolar affinity for RU 486. MR- and GR-like receptors were found in equal numbers in whole brain assays (average B(max)=69 and 62 fmol/mg protein, respectively). House sparrow brain membranes contain a single binding site specific for glucocorticoids, with characteristics consistent with a steroid/receptor interaction. Corticosterone affinity for this putative membrane receptor was approximately 24 nM, with apparent B(max)=177 fmol/mg protein. House sparrow plasma contained a single binding site for [(3)H]corticosterone. Specific binding to plasma sites was inhibited by glucocorticoids, progesterone, and testosterone. Testosterone binding to this corticosteroid binding globulin is noteworthy as sex steroid-specific binding globulins have not been identified in birds. Taken together, these data extend our ability to evaluate the comparative actions of glucocorticoids, increase our understanding of mechanisms behind the tissue specificity of glucocorticoid action, and offer insight into the evolution of glucocorticoid action in vertebrates.

Evolutionary triplet models of structured RNA

The reconstruction and synthesis of ancestral RNAs is a feasible goal for paleogenetics. This will require new bioinformatics methods, including a robust statistical framework for reconstructing histories of substitutions, indels and structural changes. We describe a "transducer composition" algorithm for extending pairwise probabilistic models of RNA structural evolution to models of multiple sequences related by a phylogenetic tree. This algorithm draws on formal models of computational linguistics as well as the 1985 protosequence algorithm of David Sankoff. The output of the composition algorithm is a multiple-sequence stochastic context-free grammar. We describe dynamic programming algorithms, which are robust to null cycles and empty bifurcations, for parsing this grammar. Example applications include structural alignment of non-coding RNAs, propagation of structural information from an experimentally-characterized sequence to its homologs, and inference of the ancestral structure of a set of diverged RNAs. We implemented the above algorithms for a simple model of pairwise RNA structural evolution; in particular, the algorithms for maximum likelihood (ML) alignment of three known RNA structures and a known phylogeny and inference of the common ancestral structure. We compared this ML algorithm to a variety of related, but simpler, techniques, including ML alignment algorithms for simpler models that omitted various aspects of the full model and also a posterior-decoding alignment algorithm for one of the simpler models. In our tests, incorporation of basepair structure was the most important factor for accurate alignment inference; appropriate use of posterior-decoding was next; and fine details of the model were least important. Posterior-decoding heuristics can be substantially faster than exact phylogenetic inference, so this motivates the use of sum-over-pairs heuristics where possible (and approximate sum-over-pairs). For more exact probabilistic inference, we discuss the use of transducer composition for ML (or MCMC) inference on phylogenies, including possible ways to make the core operations tractable.

A synthetic protein selected for ligand binding affinity mediates ATP hydrolysis

How primitive enzymes emerged from a primordial pool remains a fundamental unanswered question with important practical implications in synthetic biology. Here we show that a de novo evolved ATP binding protein, selected solely on the basis of its ability to bind ATP, mediates the regiospecific hydrolysis of ATP to ADP when crystallized with 1 equiv of ATP. Structural insights into this reaction were obtained by growing protein crystals under saturating ATP conditions. The resulting crystal structure refined to 1.8 A resolution reveals that this man-made protein binds ATP in an unusual bent conformation that is metal-independent and held in place by a key bridging water molecule. Removal of this interaction using a null mutant results in a variant that binds ATP in a normal linear geometry and is incapable of ATP hydrolysis. Biochemical analysis, including high-resolution mass spectrometry performed on dissolved protein crystals, confirms that the reaction is accelerated in the crystalline environment. This observation suggests that proteins with weak chemical reactivity can emerge from high affinity ligand binding sites and that constrained ligand-binding geometries could have helped to facilitate the emergence of early protein enzymes.

Metabolic evolution suggests an explanation for the weakness of antioxidant defences in beta-cells

The lack of an effective antioxidant system in beta-cells, which renders them susceptible to oxidative stress, is to date without explanation. The particular weakness of beta-cells in females, in both humans and mice, is another unexplained observation. We hypothesise that reactive oxygen species (ROS) in beta-cells, by their negative effect on insulin synthesis/secretion, play a fitness-enhancing role for the whole organism. Under stress conditions, the release of stress hormones produces insulin resistance and, owing to ROS preventing beta-cells from secreting insulin at the level required to maintain homeostasis, diverts glucose to insulin-independent tissues such as the brain and the foetus. We suggest that pancreatic beta-cells lost part of their antioxidant defence in association with brain evolution, and lost even more in females when placental mammals evolved. The unusual antioxidant status of beta-cells may thus be explained as an instance of co-evolution of the brain, cortisol and corticosteroid receptors, and beta-cells in the endocrine pancreas.

Rapid sequence evolution of transcription factors controlling neuron differentiation in Caenorhabditis

Whether phenotypic evolution proceeds predominantly through changes in regulatory sequences is a controversial issue in evolutionary genetics. Ample evidence indicates that the evolution of gene regulatory networks via changes in cis-regulatory sequences is an important determinant of phenotypic diversity. However, recent experimental work suggests that the role of transcription factor (TF) divergence in developmental evolution may be underestimated. In order to help understand what levels of constraints are acting on the coding sequence of developmental regulatory genes, evolutionary rates were investigated among 48 TFs required for neuronal development in Caenorhabditis elegans. Allelic variation was then sampled for 28 of these genes within a population of the related species Caenorhabditis remanei. Neuronal TFs are more divergent, both within and between species, than structural genes. TFs affecting different neuronal classes are under different levels of selective constraints. The regulatory genes controlling the differentiation of chemosensory neurons evolve particularly fast and exhibit higher levels of within- and between-species nucleotide variation than TFs required for the development of several neuronal classes and TFs required for motorneuron differentiation. The TFs affecting chemosensory neuron development are also more divergent than chemosensory genes expressed in the neurons they differentiate. These results illustrate that TFs are not as highly constrained as commonly thought and suggest that the role of divergence in developmental regulatory genes during the evolution of gene regulatory networks requires further attention.

In the light of directed evolution: pathways of adaptive protein evolution

Directed evolution is a widely-used engineering strategy for improving the stabilities or biochemical functions of proteins by repeated rounds of mutation and selection. These experiments offer empirical lessons about how proteins evolve in the face of clearly-defined laboratory selection pressures. Directed evolution has revealed that single amino acid mutations can enhance properties such as catalytic activity or stability and that adaptation can often occur through pathways consisting of sequential beneficial mutations. When there are no single mutations that improve a particular protein property experiments always find a wealth of mutations that are neutral with respect to the laboratory-defined measure of fitness. These neutral mutations can open new adaptive pathways by at least 2 different mechanisms. Functionally-neutral mutations can enhance a protein's stability, thereby increasing its tolerance for subsequent functionally beneficial but destabilizing mutations. They can also lead to changes in "promiscuous" functions that are not currently under selective pressure, but can subsequently become the starting points for the adaptive evolution of new functions. These lessons about the coupling between adaptive and neutral protein evolution in the laboratory offer insight into the evolution of proteins in nature.

Hormone-activated estrogen receptors in annelid invertebrates: implications for evolution and endocrine disruption

As the primary mediators of estrogen signaling in vertebrates, estrogen receptors (ERs) play crucial roles in reproduction, development, and behavior. They are also the major mediators of endocrine disruption by xenobiotic pollutants that mimic or block estrogen action. ERs that are sensitive to estrogen and endocrine disrupters have long been thought to be restricted to vertebrates: although there is evidence for estrogen signaling in invertebrates, the only ERs studied to date, from mollusks and cephalochordates, have been insensitive to estrogen and therefore incapable of mediating estrogen signaling or disruption. To determine whether estrogen sensitivity is ancestral or a unique characteristic of vertebrate ERs, we isolated and characterized ERs from two annelids, Platynereis dumerilii and Capitella capitata, because annelids are the sister phylum to mollusks and have been shown to produce and respond to estrogens. Functional assays show that annelid ERs specifically activate transcription in response to low estrogen concentrations and bind estrogen with high affinity. Furthermore, numerous known endocrine-disrupting chemicals activate or antagonize the annelid ER. This is the first report of a hormone-activated invertebrate ER. Our results indicate that estrogen signaling via the ER is as ancient as the ancestral bilaterian animal and corroborate the estrogen sensitivity of the ancestral steroid receptor. They suggest that the taxonomic scope of endocrine disruption by xenoestrogens may be very broad and reveal how functional diversity evolved in a gene family central to animal endocrinology.

Enzyme (re)design: lessons from natural evolution and computation

The (re)design of enzymes to catalyze 'new' reactions is a topic of considerable practical and intellectual interest. Directed evolution (random mutagenesis followed by screening/selection) has been used widely to identify novel biocatalysts. However, 'rational' approaches using either natural divergent evolution or computational predictions based on chemical principles have been less successful. This review summarizes recent progress in evolution-based and computation-based (re)design.