Adaptation of organisms to environmental niches is a hallmark of evolution. One prevalent example is that of thermal adaptation, wherein two descendants evolve at different temperature extremes1,2. Underlying the physiological differences between such organisms are changes in enzymes catalyzing essential reactions3, with orthologues from each organism undergoing adaptive mutations that preserve similar catalytic rates at their respective physiological temperatures 4,5. The sequence changes responsible for these adaptive differences, however, are often at surface exposed sites distant from the substrate binding site, leaving the active site of the enzyme structurally unperturbed6,7. How such changes are allosterically propagated to the active site, to modulate activity, is not known. Here we show that entropy-tuning changes can be engineered into distal sites of Escherichia coli adenylate kinase (AK) to quantitatively assess the role of dynamics in determining affinity, turnover, and the role in driving adaptation. The results not only reveal a dynamics-based allosteric tuning mechanism, but also uncover a spatial separation of the control of key enzymatic parameters. Fluctuations in one mobile domain (i.e. the LID) control substrate affinity, while dynamic attenuation in the other (i.e. the AMPbd) affects rate-limiting conformational changes governing enzyme turnover. Dynamics-based regulation may thus represent an elegant, widespread, and previously unrealized evolutionary adaptation mechanism that fine-tunes biological function without altering the ground state structure. Furthermore, because rigid-body conformational changes in both domains were thought to be rate limiting for turnover8,9, these adaptation studies reveal a new paradigm for understanding the relationship between dynamics and turnover in AK.

Many regulatory proteins bind peptide regions of target proteins and modulate their activity. Such regulatory proteins can often interact with highly diverse target peptides. In many instances, it is not known if the peptide-binding interface discriminates targets in a biological context, or whether biological specificity is achieved exclusively through external factors such as subcellular localization. We used an evolutionary biochemical approach to distinguish these possibilities for two such low-specificity proteins: S100A5 and S100A6. We used isothermal titration calorimetry to study the binding of peptides with diverse sequence and biochemistry to human S100A5 and S100A6. These proteins bound distinct, but overlapping, sets of peptide targets. We then studied the peptide binding properties of orthologs sampled from across five amniote species. Binding specificity was conserved along all lineages, for the last 320 million years, despite the low specificity of each protein. We used ancestral sequence reconstruction to determine the binding specificity of the last common ancestor of the paralogs. The ancestor bound the entire set of peptides bound by modern S100A5 and S100A6 proteins, suggesting that paralog specificity evolved via subfunctionalization. To rule out the possibility that specificity is conserved because it is difficult to modify, we identified a single historical mutation that, when reverted in human S100A5, gave it the ability to bind an S100A6-specific peptide. These results reveal strong evolutionary constraints on peptide binding specificity. Despite being able to bind a large number of targets, the specificity of S100 peptide interfaces is likely important for the biology of these proteins.

Copper (Cu), iron (Fe), and thyroid hormone (TH) deficiencies produce similar defects in late brain development including hypomyelination of axons and impaired synapse formation and function, suggesting that these micronutrient deficiencies share a common mechanism contributing to these derangements. We previously demonstrated that fetal/neonatal Cu and Fe deficiencies lower circulating TH concentrations in neonatal rats. Fe deficiency also reduces whole-brain T(3) content, suggesting impaired TH action in the developing Fe-deficient brain. We hypothesized that fetal/neonatal Cu and Fe deficiencies will produce mild or moderate TH deficiencies and will impair TH-responsive gene expression in the neonatal cerebral cortex and hippocampus. To test this hypothesis, we rendered pregnant Sprague Dawley rats Cu-, Fe-, or TH-deficient from early gestation through postnatal d 10 (P10). Mild and moderate TH deficiencies were induced by 1 and 3 ppm propylthiouracil treatment, respectively. Cu deficiency did not significantly alter serum or tissue TH concentrations or TH-responsive brain mRNA expression. Fe deficiency significantly lowered P10 serum total T(3) (45%), serum total T(4) (52%), whole brain T(3) (14%), and hippocampal T(3) (18%) concentrations, producing a mild TH deficiency similar to 1 ppm propylthiouracil treatment. Fe deficiency lowered Pvalb, Enpp6, and Mbp mRNA levels in the P10 hippocampus. Fe deficiency also altered Hairless, Dbm, and Dio2 mRNA levels in the P10 cerebral cortex. These results suggest that some of the brain defects associated with Fe deficiency may be mediated through altered thyroidal status and the concomitant alterations in TH-responsive gene transcription.

No abstract

The amino acid sequences of proteins have evolved over billions of years, preserving their structures and functions while responding to evolutionary forces. Are there conserved sequence and structural elements that preserve the protein folding mechanisms? The functionally diverse and ancient (βα)1–8 TIM barrel motif may answer this question. We mapped the complex six-state folding free energy surface of a ∼3.6 billion y old, bacterial indole-3-glycerol phosphate synthase (IGPS) TIM barrel enzyme by equilibrium and kinetic hydrogen–deuterium exchange mass spectrometry (HDX-MS). HDX-MS on the intact protein reported exchange in the native basin and the presence of two thermodynamically distinct on- and off-pathway intermediates in slow but dynamic equilibrium with each other. Proteolysis revealed protection in a small (α1β2) and a large cluster (β5α5β6α6β7) and that these clusters form cores of stability in Ia and Ibp. The strongest protection in both states resides in β4α4 with the highest density of branched aliphatic side chain contacts in the folded structure. Similar correlations were observed previously for an evolutionarily distinct archaeal IGPS, emphasizing a key role for hydrophobicity in stabilizing common high-energy folding intermediates. A bioinformatics analysis of IGPS sequences from the three superkingdoms revealed an exceedingly high hydrophobicity and surprising α-helix propensity for β4, preceded by a highly conserved βα-hairpin clamp that links β3 and β4. The conservation of the folding mechanisms for archaeal and bacterial IGPS proteins reflects the conservation of key elements of sequence and structure that first appeared in the last universal common ancestor of these ancient proteins.