Vestigialization of an Allosteric Switch: Genetic and Structural Mechanisms for the Evolution of Constitutive Activity in a Steroid Hormone Receptor

An important goal in molecular evolution is to understand the genetic and physical mechanisms by which protein functions evolve and, in turn, to characterize how a protein's physical architecture influences its evolution. Here we dissect the mechanisms for an evolutionary shift in function in the mollusk ortholog of the steroid hormone receptors (SRs), a family of biologically essential transcription factors. In vertebrates, the activity of SRs allosterically depends on binding a hormonal ligand; in mollusks, however, the SR ortholog (called ER, because of high sequence similarity to vertebrate estrogen receptors) activates transcription in the absence of ligand and does not respond to steroid hormones. To understand how this shift in regulation evolved, we combined evolutionary, structural, and functional analyses. We first determined the X-ray crystal structure of the ER of the Pacific oyster Crassostrea gigas (CgER), and found that its ligand pocket is filled with bulky residues that prevent ligand occupancy. To understand the genetic basis for the evolution of mollusk ERs' unique functions, we resurrected an ancient SR progenitor and characterized the effect of historical amino acid replacements on its functions. We found that reintroducing just two ancient replacements from the lineage leading to mollusk ERs recapitulates the evolution of full constitutive activity and the loss of ligand activation. These substitutions stabilize interactions among key helices, causing the allosteric switch to become “stuck” in the active conformation and making activation independent of ligand binding. Subsequent changes filled the ligand pocket without further affecting activity; by degrading the allosteric switch, these substitutions vestigialized elements of the protein's architecture required for ligand regulation and made reversal to the ancestral function more complex. These findings show how the physical architecture of allostery enabled a few large-effect mutations to trigger a profound evolutionary change in the protein's function and shaped the genetics of evolutionary reversibility.</p

Pervasive contingency and entrenchment in a billion years of Hsp90 evolution [preprint]

Interactions among mutations within a protein have the potential to make molecular evolution contingent and irreversible, but the extent to which epistasis actually shaped historical evolutionary trajectories is unclear. We addressed this question by identifying all amino acid substitutions that occurred during the billion-year evolutionary history of the heat shock protein 90 (Hsp90) ATPase domain beginning from a deep eukaryotic ancestor to modern Saccharomyces cerevisiae and then precisely measuring their fitness effects when introduced into both extant and reconstructed ancestral Hsp90 proteins. We find a pervasive influence of epistasis: of 98 derived states that evolved during history, most were deleterious at times before they happened, and the vast majority also became subsequently entrenched, with the ancestral state becoming deleterious after its substitution. This epistasis was primarily caused by specific interactions among sites rather than a general permissive or restrictive effect on the protein\u27s tolerance to mutation. Our results show that epistasis continually opens and closes windows of mutational opportunity over evolutionary timescales, producing histories and biological states that reflect the transient internal constraints imposed by a protein\u27s fleeting sequence states

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

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

Family of neural wiring receptors in bilaterians defined by phylogenetic, biochemical, and structural evidence

The evolution of complex nervous systems was accompanied by the expansion of numerous protein families, including cell-adhesion molecules, surface receptors, and their ligands. These proteins mediate axonal guidance, synapse targeting, and other neuronal wiring-related functions. Recently, 32 interacting cell surface proteins belonging to two newly defined families of the Ig superfamily (IgSF) in fruit flies were discovered to label different subsets of neurons in the brain and ventral nerve cord. They have been shown to be involved in synaptic targeting and morphogenesis, retrograde signaling, and neuronal survival. Here, we show that these proteins, Dprs and DIPs, are members of a widely distributed family of two- and three-Ig domain molecules with neuronal wiring functions, which we refer to as Wirins. Beginning from a single ancestral Wirin gene in the last common ancestor of Bilateria, numerous gene duplications produced the heterophilic Dprs and DIPs in protostomes, along with two other subfamilies that diversified independently across protostome phyla. In deuterostomes, the ancestral Wirin evolved into the IgLON subfamily of neuronal receptors. We show that IgLONs interact with each other and that their complexes can be broken by mutations designed using homology models based on Dpr and DIP structures. The nematode orthologs ZIG-8 and RIG-5 also form heterophilic and homophilic complexes, and crystal structures reveal numerous apparently ancestral features shared with Dpr-DIP complexes. The evolutionary, biochemical, and structural relationships we demonstrate here provide insights into neural development and the rise of the metazoan nervous system

Assessing simulation ecosystem processes for climate variability research at Glacier National Park

Glacier National Park served as a test site for ecosystem analyses that involved a suite of integrated models embedded within a geographic information system. The goal of the exercise was to provide managers with maps that could illustrate probable shifts in vegetation, net primary production (NPP), and hydrologic responses associated with two selected climatic scenarios. The climatic scenarios were (a) a recent 12-yr record of weather data, and (b) a reconstituted set that sequentially introduced in repeated 3-yr intervals wetter–cooler, drier–warmer, and typical conditions. To extrapolate the implications of changes in ecosystem processes and resulting growth and distribution of vegetation and snowpack, the model incorporated geographic data. With underlying digital elevation maps, soil depth and texture, extrapolated climate, and current information on vegetation types and satellite-derived estimates of leaf area indices, simulations were extended to envision how the park might look after 120 yr. The predictions of change included underlying processes affecting the availability of water and nitrogen. Considerable field data were acquired to compare with model predictions under current climatic conditions. In general, the integrated landscape models of ecosystem processes had good agreement with measured NPP, snowpack, and streamflow, but the exercise revealed the difficulty and necessity of averaging point measurements across landscapes to achieve comparable results with modeled values. Under the extremely variable climate scenario significant changes in vegetation composition and growth as well as hydrologic responses were predicted across the park. In particular, a general rise in both the upper and lower limits of treeline was predicted. These shifts would probably occur along with a variety of disturbances (fire, insect, and disease outbreaks) as predictions of physiological stress (water, nutrients, light) altered competitive relations and hydrologic responses. The use of integrated landscape models applied in this exercise should provide managers with insights into the underlying processes important in maintaining community structure, and at the same time, locate where changes on the landscape are most likely to occur

Pyruvate Oxidase of \u3ci\u3eStreptococcus pneumoniae\u3c/i\u3e Contributes to Penumolysin Release

Background Streptococcus pneumoniae is one of the leading causes of community acquired pneumonia and acute otitis media. Certain aspects of S. pneumoniae’s virulence are dependent upon expression and release of the protein toxin pneumolysin (PLY) and upon the activity of the peroxide-producing enzyme, pyruvate oxidase (SpxB). We investigated the possible synergy of these two proteins and identified that release of PLY is enhanced by expression of SpxB prior to stationary phase growth. Results Mutants lacking the \u3c\u3espxB gene were defective in PLY release and complementation of spxB restored PLY release. This was demonstrated by cytotoxic effects of sterile filtered supernatants upon epithelial cells and red blood cells. Additionally, peroxide production appeared to contribute to the mechanism of PLY release since a significant correlation was found between peroxide production and PLY release among a panel of clinical isolates. Exogenous addition of H2O2 failed to induce PLY release and catalase supplementation prevented PLY release in some strains, indicating peroxide may exert its effect intracellularly or in a strain-dependent manner. SpxB expression did not trigger bacterial cell death or LytA-dependent autolysis, but did predispose cells to deoxycholate lysis. Conclusions Here we demonstrate a novel link between spxB expression and PLY release. These findings link liberation of PLY toxin to oxygen availability and pneumococcal metabolism

Family of neural wiring receptors in bilaterians defined by phylogenetic, biochemical, and structural evidence

The evolution of complex nervous systems was accompanied by the expansion of numerous protein families, including cell-adhesion molecules, surface receptors, and their ligands. These proteins mediate axonal guidance, synapse targeting, and other neuronal wiring-related functions. Recently, 32 interacting cell surface proteins belonging to two newly defined families of the Ig superfamily (IgSF) in fruit flies were discovered to label different subsets of neurons in the brain and ventral nerve cord. They have been shown to be involved in synaptic targeting and morphogenesis, retrograde signaling, and neuronal survival. Here, we show that these proteins, Dprs and DIPs, are members of a widely distributed family of two- and three-Ig domain molecules with neuronal wiring functions, which we refer to as Wirins. Beginning from a single ancestral Wirin gene in the last common ancestor of Bilateria, numerous gene duplications produced the heterophilic Dprs and DIPs in protostomes, along with two other subfamilies that diversified independently across protostome phyla. In deuterostomes, the ancestral Wirin evolved into the IgLON subfamily of neuronal receptors. We show that IgLONs interact with each other and that their complexes can be broken by mutations designed using homology models based on Dpr and DIP structures. The nematode orthologs ZIG-8 and RIG-5 also form heterophilic and homophilic complexes, and crystal structures reveal numerous apparently ancestral features shared with Dpr-DIP complexes. The evolutionary, biochemical, and structural relationships we demonstrate here provide insights into neural development and the rise of the metazoan nervous system

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