Recruitment of a Middling Promiscuous Enzyme Drives Adaptive Metabolic Evolution in Escherichia coli.

A key step in metabolic pathway evolution is the recruitment of promiscuous enzymes to perform new functions. Despite the recognition that promiscuity is widespread in biology, factors dictating the preferential recruitment of one promiscuous enzyme over other candidates are unknown. Escherichia coli contains four sugar kinases that are candidates for recruitment when the native glucokinase machinery is deleted-allokinase (AlsK), manno(fructo)kinase (Mak), N-acetylmannosamine kinase (NanK), and N-acetylglucosamine kinase (NagK). The catalytic efficiencies of these enzymes are 103- to 105-fold lower than native glucokinases, ranging from 2,400 M-1 s-1 for the most active candidate, NagK, to 15 M-1 s-1 for the least active candidate, AlsK. To investigate the relationship between catalytic activities of promiscuous enzymes and their recruitment, we performed adaptive evolution of a glucokinase-deficient E. coli strain to restore glycolytic metabolism. We observed preferential recruitment of NanK via a trajectory involving early mutations that facilitate glucose uptake and amplify nanK transcription, followed by nonsynonymous substitutions in NanK that enhance the enzyme's promiscuous glucokinase activity. These substitutions reduced the native activity of NanK and reduced organismal fitness during growth on an N-acetylated carbon source, indicating that enzyme recruitment comes at a cost for growth on other substrates. Notably, the two most active candidates, NagK and Mak, were not recruited, suggesting that catalytic activity alone does not dictate evolutionary outcomes. The results highlight our lack of knowledge regarding biological drivers of enzyme recruitment and emphasize the need for a systems-wide approach to identify factors facilitating or constraining this important adaptive process.

Biochemical methods to map and quantify allosteric motions in human glucokinase.

Allosteric regulation of protein function is ubiquitous in biology. Allostery originates from ligand-mediated alterations in polypeptide structure and/or dynamics, which produce a cooperative kinetic or thermodynamic response to changing ligand concentrations. Establishing a mechanistic description of individual allosteric events requires both mapping the relevant changes in protein structure and quantifying the rates of differential conformational dynamics in the absence and presence of effectors. In this chapter, we describe three biochemical approaches to understand the dynamic and structural signatures of protein allostery using the well-established cooperative enzyme glucokinase as a case study. The combined application of pulsed proteolysis, biomolecular nuclear magnetic resonance spectroscopy and hydrogen-deuterium exchange mass spectrometry offers complementary information that can used to establish molecular models for allosteric proteins, especially when differential protein dynamics are involved.

The Chemosensory Repertoire of the Eastern Diamondback Rattlesnake (Crotalus adamanteus) Reveals Complementary Genetics of Olfactory and Vomeronasal-Type Receptors.

Pitviper sensory perception incorporates diverse stimuli through the integration of trichromatic color vision, bifocal heat-sensing, and dual-system chemoperception. Chemoperception, or olfaction, is mediated by chemoreceptors in the olfactory bulb and the vomeronasal organ, but the true genomic complexity of the gene families and their relative contributions is unknown. A full genomic accounting of pitviper chemoperception directly complements our current understanding of their venoms by generating a more complete polyphenic representation of their predatory arsenal. To characterize the genetic repertoire of pitviper chemoperception, we analyzed a full-genome assembly for Crotalus adamanteus, the eastern diamondback rattlesnake. We identified hundreds of genes encoding both olfactory receptors (ORs; 362 full-length genes) and type-2 vomeronasal receptors (V2Rs; 430 full-length genes). Many chemoreceptor genes are organized into large tandem repeat arrays. Comparative analysis of V2R orthologs across squamates demonstrates how gene array expansion and contraction underlies the evolution of the chemoreceptor repertoire, which likely reflects shifts in life history traits. Chromosomal assignments of chemosensory genes identified sex chromosome specific chemoreceptor genes, providing gene candidates underlying observed sex-specific chemosensory-based behaviors. We detected widespread episodic evolution in the extracellular, ligand-binding domains of both ORs and V2Rs, suggesting the diversification of chemoreceptors is driven by transient periods of positive selection. We provide a robust genetic framework for studying pitviper chemosensory ecology and evolution.

Nanosecond-Timescale Dynamics and Conformational Heterogeneity in Human GCK Regulation and Disease.

Human glucokinase (GCK) is the prototypic example of an emerging class of proteins with allosteric-like behavior that originates from intrinsic polypeptide dynamics. High-resolution NMR investigations of GCK have elucidated millisecond-timescale dynamics underlying allostery. In contrast, faster motions have remained underexplored, hindering the development of a comprehensive model of cooperativity. Here, we map nanosecond-timescale dynamics and structural heterogeneity in GCK using a combination of unnatural amino acid incorporation, time-resolved fluorescence, and 19F nuclear magnetic resonance spectroscopy. We find that a probe inserted within the enzyme's intrinsically disordered loop samples multiple conformations in the unliganded state. Glucose binding and disease-associated mutations that suppress cooperativity alter the number and/or relative population of these states. Together, the nanosecond kinetics characterized here and the millisecond motions known to be essential for cooperativity provide a dynamical framework with which we address the origins of cooperativity and the mechanism of activated, hyperinsulinemia-associated, noncooperative variants.

Biophysical Spandrels form a Hot-Spot for Kosmotropic Mutations in Bacteriophage Thermal Adaptation.

Temperature plays a dominating role in protein structure and function, and life has evolved myriad strategies to adapt proteins to environmental thermal stress. Cellular systems can utilize kosmotropic osmolytes, the products of complex biochemical pathways, to act as chemical chaperones. These extrinsic molecules, e.g., trehalose, alter local water structure to modulate the strength of the hydrophobic effect and increase protein stability. In contrast, simpler genetic systems must rely on intrinsic mutation to affect protein stability. In naturally occurring microvirid bacteriophages of the subfamily Bullavirinae, capsid stability is randomly distributed across the phylogeny, suggesting it is not phylogenetically linked and could be altered through adaptive mutation. We hypothesized that these phages could utilize an adaptive mechanism that mimics the stabilizing effects of the kosmotrope trehalose through mutation. Kinetic stability of wild-type ID8, a relative of ΦX174, displays a saturable response to trehalose. Thermal adaptation mutations in ID8 improve capsid stability and reduce responsiveness to trehalose suggesting the mutations move stability closer to the kosmotropic saturation point, mimicking the kosmotropic effect of trehalose. These mutations localize to and modulate the hydrophobicity of a cavern formation at the interface of phage coat and spike proteins-an evolutionary spandrel. Across a series of genetically distinct phages, responsiveness to trehalose correlates positively with cavern hydrophobicity suggesting that the level of hydrophobicity of the cavern may provide a biophysical gating mechanism constraining or permitting adaptation in a lineage-specific manner. Our results demonstrate that a single mutation can exploit pre-existing, non-adaptive structural features to mimic the adaptive effects of complex biochemical pathways.

A Single Mutation Unlocks Cascading Exaptations in the Origin of a Potent Pitviper Neurotoxin.

Evolutionary innovations and complex phenotypes seemingly require an improbable amount of genetic change to evolve. Rattlesnakes display two dramatically different venom phenotypes. Type I venoms are hemorrhagic with low systemic toxicity and high expression of tissue-destroying snake venom metalloproteinases. Type II venoms are highly neurotoxic and lack snake venom metalloproteinase expression and associated hemorrhagic activity. This dichotomy hinges on Mojave toxin (MTx), a phospholipase A2 (PLA2) based ß-neurotoxin expressed in Type II venoms. MTx is comprised of a nontoxic acidic subunit that undergoes extensive proteolytic processing and allosterically regulates activity of a neurotoxic basic subunit. Evolution of the acidic subunit presents an evolutionary challenge because the need for high expression of a nontoxic venom component and the proteolytic machinery required for processing suggests genetic changes of seemingly little immediate benefit to fitness. We showed that MTx evolved through a cascading series of exaptations unlocked by a single nucleotide change. The evolution of one new cleavage site in the acidic subunit unmasked buried cleavage sites already present in ancestral PLA2s, enabling proteolytic processing. Snake venom serine proteases, already present in the venom to disrupt prey hemostasis, possess the requisite specificities for MTx acidic subunit proteolysis. The dimerization interface between MTx subunits evolved by exploiting a latent, but masked, hydrophobic interaction between ancestral PLA2s. The evolution of MTx through exaptation of existing functional and structural features suggests complex phenotypes that depend on evolutionary innovations can arise from minimal genetic change enabled by prior evolution.

Dual allosteric activation mechanisms in monomeric human glucokinase.

Cooperativity in human glucokinase (GCK), the body's primary glucose sensor and a major determinant of glucose homeostatic diseases, is fundamentally different from textbook models of allostery because GCK is monomeric and contains only one glucose-binding site. Prior work has demonstrated that millisecond timescale order-disorder transitions within the enzyme's small domain govern cooperativity. Here, using limited proteolysis, we map the site of disorder in unliganded GCK to a 30-residue active-site loop that closes upon glucose binding. Positional randomization of the loop, coupled with genetic selection in a glucokinase-deficient bacterium, uncovers a hyperactive GCK variant with substantially reduced cooperativity. Biochemical and structural analysis of this loop variant and GCK variants associated with hyperinsulinemic hypoglycemia reveal two distinct mechanisms of enzyme activation. In α-type activation, glucose affinity is increased, the proteolytic susceptibility of the active site loop is suppressed and the (1)H-(13)C heteronuclear multiple quantum coherence (HMQC) spectrum of (13)C-Ile-labeled enzyme resembles the glucose-bound state. In ß-type activation, glucose affinity is largely unchanged, proteolytic susceptibility of the loop is enhanced, and the (1)H-(13)C HMQC spectrum reveals no perturbation in ensemble structure. Leveraging both activation mechanisms, we engineer a fully noncooperative GCK variant, whose functional properties are indistinguishable from other hexokinase isozymes, and which displays a 100-fold increase in catalytic efficiency over wild-type GCK. This work elucidates specific structural features responsible for generating allostery in a monomeric enzyme and suggests a general strategy for engineering cooperativity into proteins that lack the structural framework typical of traditional allosteric systems.

Functional and structural characterization of a eurytolerant calsequestrin from the intertidal teleost Fundulus heteroclitus.

Calsequestrins (CSQ) are high capacity, medium affinity, calcium-binding proteins present in the sarcoplasmic reticulum (SR) of cardiac and skeletal muscles. CSQ sequesters Ca²âº during muscle relaxation and increases the Ca²âº-storage capacity of the SR. Mammalian CSQ has been well studied as a model of human disease, but little is known about the environmental adaptation of CSQ isoforms from poikilothermic organisms. The mummichog, Fundulus heteroclitus, is an intertidal fish that experiences significant daily and seasonal environmental fluctuations and is an interesting study system for investigations of adaptation at the protein level. We determined the full-length coding sequence of a CSQ isoform from skeletal muscle of F. heteroclitus (FCSQ) and characterized the function and structure of this CSQ. The dissociation constant (K(d)) of FCSQ is relatively insensitive to changes in temperature and pH, thus indicating that FCSQ is a eurytolerant protein. We identified and characterized a highly conserved salt bridge network in FCSQ that stabilizes the formation of front-to-front dimers, a process critical to CSQ function. The functional profile of FCSQ correlates with the natural history of F. heteroclitus suggesting that the eurytolerant function of FCSQ may be adaptive.

Resurrecting prehistoric parvalbumins to explore the evolution of thermal compensation in extant Antarctic fish parvalbumins.

Parvalbumins (PVs) from Antarctic notothenioid fishes display a pattern of thermal adaptation that likely reflects evolutionary changes in protein conformational flexibility. We have used ancestral sequence reconstruction and homology modeling to identify two amino acid changes that could potentially account for the present thermal sensitivity pattern of Antarctic fish PVs compared with a PV from a theoretical warm-adapted ancestral fish. To test this hypothesis, ancient PVs were resurrected in the lab using PV from the notothenioid Gobionotothen gibberifrons as a platform for introducing mutations comparable to the reconstructed ancestral PV sequences. The wild-type PV (WT) as well as three mutant expression constructs were engineered: lysine 8 to asparagine (K8N), lysine 26 to asparagine (K26N) and a double mutant (DM). Calcium equilibrium dissociation constants (K(d)) versus temperature curves for all mutants were right-shifted, as predicted, relative to that of WT PV. The K(d) values for the K8N and K26N single mutants were virtually identical at all temperatures and showed an intermediate level of thermal sensitivity. The DM construct displayed a full conversion of thermal sensitivity pattern to that of a PV from a warm/temperate-adapted fish. Additionally, the K(d) versus temperature curve for the WT construct revealed greater thermal sensitivity compared with the mutant constructs. Measurements of the rates of Ca(2+) dissociation (k(off)) showed that all mutants generally had slower k(off) values than WT at all temperatures. Calculated rates of Ca(2+) binding (k(on)) for the K8N and K26N mutants were similar to values for the WT PV at all temperatures. In contrast, the calculated k(on) values for the DM PV were faster, providing mechanistic insights into the nature of potentially adaptive changes in Ca(2+) binding in this PV. The overall results suggest that the current thermal phenotype of Antarctic PVs can be recapitulated by just two amino acid substitutions.