Relevant Conformational Changes Research Articles

3D variability analysis reveals a hidden conformational change controlling ammonia transport in human asparagine synthetase

Advances in X-ray crystallography and cryogenic electron microscopy (cryo-EM) offer the promise of elucidating functionally relevant conformational changes that are not easily studied by other biophysical methods. Here we show that 3D variability analysis (3DVA) of the cryo-EM map for wild-type (WT) human asparagine synthetase (ASNS) identifies a functional role for the Arg-142 side chain and test this hypothesis experimentally by characterizing the R142I variant in which Arg-142 is replaced by isoleucine. Support for Arg-142 playing a role in the intramolecular translocation of ammonia between the active site of the enzyme is provided by the glutamine-dependent synthetase activity of the R142 variant relative to WT ASNS, and MD simulations provide a possible molecular mechanism for these findings. Combining 3DVA with MD simulations is a generally applicable approach to generate testable hypotheses of how conformational changes in buried side chains might regulate function in enzymes.

Adaptive resolution force probe simulations: Coarse graining in the ideal gas approximation.

The unfolding of molecular complexes or biomolecules under the influence of external mechanical forces can routinely be simulated with atomistic resolution. To obtain a match of the characteristic time scales with those of experimental force spectroscopy, often coarse graining procedures are employed. Here, building on a previous study, we apply the adaptive resolution scheme (AdResS) to force probe molecular dynamics (FPMD) simulations using two model systems as examples: One system is the previously investigated calix[4]arene dimer that shows reversible one-step unfolding, and the other example is provided by a small peptide, a β-alanine octamer in methanol solvent. The mechanical unfolding of this peptide proceeds via a metastable intermediate and, therefore, represents a first step toward a complex unfolding pathway. We show that the average number of native contacts serves as a robust order parameter for the forced unfolding of this small peptide. In addition to increasing the complexity of the relevant conformational changes, we study the impact of the methodology used for coarse graining. Apart from the iterative Boltzmann inversion method, we apply an ideal gas approximation, and therefore, we replace the solvent by a non-interacting system of spherical particles. In all cases, we find excellent agreement between the results of FPMD simulations performed fully atomistically and those of the AdResS simulations also in the case of fast pulling. This holds for all details of the unfolding pathways, such as the distributions of the characteristic forces and also the sequence of hydrogen-bond opening in case of the β-alanine octamer. Therefore, the methodology is very well suited to simulate the mechanical unfolding of systems of experimental relevance.

3D Variability Analysis Reveals a Hidden Conformational Change Controlling Ammonia Transport in Human Asparagine Synthetase.

How motions in enzymes might be linked to catalytic function is of considerable general interest. Advances in X-ray crystallography and cryogenic electron microscopy (cryo-EM) offer the promise of elucidating functionally relevant conformational changes that are not easily studied by other biophysical methods. Here we use 3D variability analysis (3DVA) of the cryo-EM map for wild-type (WT) human asparagine synthetase (ASNS) to identify a functional role for the Arg-142 side chain as a gate that mediates ammonia access to a catalytically relevant intramolecular tunnel. Our 3DVA-derived hypothesis is assessed experimentally, using the R142I variant in which Arg-142 is replaced by isoleucine, and by molecular dynamics (MD) simulations on independent, computational models of the WT human ASNS monomer and its catalytically relevant, ternary complex with β-aspartyl-AMP and MgPPi. Residue fluctuations in the MD trajectories for the human ASNS monomer are consistent with those determined for 3DVA-derived structures. These MD simulations also indicate that the gating function of Arg-142 is separate from the molecular events that form a continuous tunnel linking the two active sites. Experimental support for Arg-142 playing a role in intramolecular ammonia translocation is provided by the glutamine-dependent synthetase activity of the R142 variant relative to WT ASNS. MD simulations of computational models for the R142I variant and the R142I/β-aspartyl-AMP/MgPPi ternary complex provide a possible molecular basis for this observation. Overall, the combination of 3DVA with MD simulations is a generally applicable approach to generate testable hypotheses of how conformational changes in buried side chains might regulate function in enzymes.

Assessing nanobody interaction with SARS-CoV-2 Nsp9.

The interaction between SARS-CoV-2 non-structural protein Nsp9 and the nanobody 2NSP90 was investigated by NMR spectroscopy using the paramagnetic perturbation methodology PENELOP (Paramagnetic Equilibrium vs Nonequilibrium magnetization Enhancement or LOss Perturbation). The Nsp9 monomer is an essential component of the replication and transcription complex (RTC) that reproduces the viral gRNA for subsequent propagation. Therefore preventing Nsp9 recruitment in RTC would represent an efficient antiviral strategy that could be applied to different coronaviruses, given the Nsp9 relative invariance. The NMR results were consistent with a previous characterization suggesting a 4:4 Nsp9-to-nanobody stoichiometry with the occurrence of two epitope pairs on each of the Nsp9 units that establish the inter-dimer contacts of Nsp9 tetramer. The oligomerization state of Nsp9 was also analyzed by molecular dynamics simulations and both dimers and tetramers resulted plausible. A different distribution of the mapped epitopes on the tetramer surface with respect to the former 4:4 complex could also be possible, as well as different stoichiometries of the Nsp9-nanobody assemblies such as the 2:2 stoichiometry suggested by the recent crystal structure of the Nsp9 complex with 2NSP23 (PDB ID: 8dqu), a nanobody exhibiting essentially the same affinity as 2NSP90. The experimental NMR evidence, however, ruled out the occurrence in liquid state of the relevant Nsp9 conformational change observed in the same crystal structure.

Multilevel superposition for deciphering the conformational variability of protein ensembles.

The dynamics and variability of protein conformations are directly linked to their functions. Many comparative studies of X-ray protein structures have been conducted to elucidate the relevant conformational changes, dynamics and heterogeneity. The rapid increase in the number of experimentally determined structures has made comparison an effective tool for investigating protein structures. For example, it is now possible to compare structural ensembles formed by enzyme species, variants or the type of ligands bound to them. In this study, the author developed a multilevel model for estimating two covariance matrices that represent inter- and intra-ensemble variability in the Cartesian coordinate space. Principal component analysis using the two estimated covariance matrices identified the inter-/intra-enzyme variabilities, which seemed to be important for the enzyme functions, with the illustrative examples of cytochrome P450 family 2 enzymes and class A $\beta$-lactamases. In P450, in which each enzyme has its own active site of a distinct size, an active-site motion shared universally between the enzymes was captured as the first principal mode of the intra-enzyme covariance matrix. In this case, the method was useful for understanding the conformational variability after adjusting for the differences between enzyme sizes. The developed method is advantageous in small ensemble-size problems and hence promising for use in comparative studies on experimentally determined structures where ensemble sizes are smaller than those generated, for example, by molecular dynamics simulations.

Insights into the Structural and Energetic Descriptions of Ubiquitin Specific Protease 7 (USP7) Catalytic Mechanisms by Hybrid QM/MM Simulations.

AbstractWe have reported the first computational study about USP7 reaction mechanism with the substrate Ubiquitin‐Rhodamine 110‐G (Ub‐Rho) using a robust methodology that integrated homology modeling, classical molecular dynamics (MD) simulations, protein‐protein interaction fingerprints (IFPs) analysis, principal components analysis (PCA) and clustering, and the hybrid QM/MM simulations using the adaptive string method (ASM). The results presented in this work could be useful to understand the dynamic nature of USP7 enzyme‐substrate complexes. In this sense, we have provided a detailed structural description of the most relevant conformational changes observed in our simulations, which could be used as reference for modeling more complicated USP7’s enzyme‐substrate complexes. On the other hand, our descriptions of the protein‐protein interactions of Ub‐Rho's residues at P5 to P1’ positions with the CD's active site, and the geometries of the stationary states located along the MFEP could be considered for designing new potent and selective inhibitors for USP7’s activity.

The Complex Roles of Adenosine Triphosphate in Bioenergetics.

ATP is generally defined as the "energy currency" of the cell. Its phosphoanhydride P-O bonds are often considered to be "high energy" linkages that release free energy when broken, and its hydrolysis is described as "strongly exergonic". However, breaking bonds cannot release energy and ATP hydrolysis in motor and active transport proteins is not "strongly exergonic". So, the relevance of ATP resides elsewhere. As important as the nucleotide are the proteins that undergo functionally relevant conformational changes upon both ATP binding and release of ADP and inorganic phosphate. ATP phosphorylates proteins for signaling, active transport, and substrates in condensation reactions. The ensuing dephosphorylation has different consequences in each case. In signaling and active transport the phosphate group is hydrolyzed whereas in condensation reactions the phosphoryl fragment acts as a dehydrating agent. As it will be discussed in this article, ATP does much more than simply contribute free energy to biological processes.

Structural basis for ligand binding modes of CTP synthase

Cytidine triphosphate synthase (CTPS), which comprises an ammonia ligase domain and a glutamine amidotransferase domain, catalyzes the final step of de novo CTP biosynthesis. The activity of CTPS is regulated by the binding of four nucleotides and glutamine. While glutamine serves as an ammonia donor for the ATP-dependent conversion of UTP to CTP, the fourth nucleotide GTP acts as an allosteric activator. Models have been proposed to explain the mechanisms of action at the active site of the ammonia ligase domain and the conformational changes derived by GTP binding. However, actual GTP/ATP/UTP binding modes and relevant conformational changes have not been revealed fully. Here, we report the discovery of binding modes of four nucleotides and a glutamine analog 6-diazo-5-oxo-L-norleucine in Drosophila CTPS by cryo-electron microscopy with near-atomic resolution. Interactions between GTP and surrounding residues indicate that GTP acts to coordinate reactions at both domains by directly blocking ammonia leakage and stabilizing the ammonia tunnel. Additionally, we observe the ATP-dependent UTP phosphorylation intermediate and determine interacting residues at the ammonia ligase. A noncanonical CTP binding at the ATP binding site suggests another layer of feedback inhibition. Our findings not only delineate the structure of CTPS in the presence of all substrates but also complete our understanding of the underlying mechanisms of the allosteric regulation and CTP synthesis.

SARS-CoV-2 simulations go exascale to predict dramatic spike opening and cryptic pockets across the proteome.

SARS-CoV-2 has intricate mechanisms for initiating infection, immune evasion/suppression, and replication, which depend on the structure and dynamics of its constituent proteins. Many protein structures have been solved, but far less is known about their relevant conformational changes. To address this challenge, over a million citizen scientists banded together through the Folding@home distributed computing project to create the first exascale computer and simulate an unprecedented 0.1 seconds of the viral proteome. Our simulations capture dramatic opening of the apo Spike complex, far beyond that seen experimentally, which explains and successfully predicts the existence of ‘cryptic’ epitopes. Different Spike homologues modulate the probabilities of open versus closed structures, balancing receptor binding and immune evasion. We also observe dramatic conformational changes across the proteome, which reveal over 50 ‘cryptic’ pockets that expand targeting options for the design of antivirals. All data and models are freely available online, providing a quantitative structural atlas.

Exploring the Conformational Landscape of a Lanthipeptide Synthetase Using Native Mass Spectrometry.

Lanthipeptides are ribosomally synthesized and post-translationally modified peptide (RiPP) natural products. These genetically encoded peptides are biosynthesized by multifunctional enzymes (lanthipeptide synthetases) that possess relaxed substrate specificity and catalyze iterative rounds of post-translational modification. Recent evidence has suggested that some lanthipeptide synthetases are structurally dynamic enzymes that are allosterically activated by precursor peptide binding and that conformational sampling of the enzyme-peptide complex may play an important role in defining the efficiency and sequence of biosynthetic events. These "biophysical" processes, while critical for defining the activity and function of the synthetase, remain very challenging to study with existing methodologies. Herein, we show that native mass spectrometry coupled to ion mobility (native IM-MS) provides a powerful and sensitive means for investigating the conformational landscapes and intermolecular interactions of lanthipeptide synthetases. Namely, we demonstrate that the class II lanthipeptide synthetase (HalM2) and its noncovalent complex with the cognate HalA2 precursor peptide can be delivered into the gas phase in a manner that preserves native structures and intermolecular enzyme-peptide contacts. Moreover, gas phase ion mobility studies of the natively folded ions demonstrate that peptide binding and mutations to dynamic structural elements of HalM2 alter the conformational landscape of the enzyme. Cumulatively, these data support previous claims that lanthipeptide synthetases are structurally dynamic enzymes that undergo functionally relevant conformational changes in response to precursor peptide binding. This work establishes native IM-MS as a versatile approach for characterizing intermolecular interactions and for unraveling the relationships between protein structure and biochemical function in RiPP biosynthetic systems.

DeSiphering receptor core-induced and ligand-dependent conformational changes in arrestin via genetic encoded trimethylsilyl 1H-NMR probe

Characterization of the dynamic conformational changes in membrane protein signaling complexes by nuclear magnetic resonance (NMR) spectroscopy remains challenging. Here we report the site-specific incorporation of 4-trimethylsilyl phenylalanine (TMSiPhe) into proteins, through genetic code expansion. Crystallographic analysis revealed structural changes that reshaped the TMSiPhe-specific amino-acyl tRNA synthetase active site to selectively accommodate the trimethylsilyl (TMSi) group. The unique up-field 1H-NMR chemical shift and the highly efficient incorporation of TMSiPhe enabled the characterization of multiple conformational states of a phospho-β2 adrenergic receptor/β-arrestin-1(β-arr1) membrane protein signaling complex, using only 5 μM protein and 20 min of spectrum accumulation time. We further showed that extracellular ligands induced conformational changes located in the polar core or ERK interaction site of β-arr1 via direct receptor transmembrane core interactions. These observations provided direct delineation and key mechanism insights that multiple receptor ligands were able to induce distinct functionally relevant conformational changes of arrestin.

Significance of the Protein Interface Configuration for Allostery in Imidazole Glycerol Phosphate Synthase.

Imidazole glycerol phosphate synthase (ImGPS) from Thermotoga maritima is a model enzyme for studying allostery. The ImGPS complex consists of the cyclase subunit HisF and the glutaminase subunit HisH whose activity is stimulated by substrate binding to HisF in a V-type manner. To investigate the significance of a putative closing hinge motion at the cyclase:glutaminase interface for HisH activity, we replaced residue W123 in HisH with the light-switchable unnatural amino acid phenylalanine-4'-azobenzene (AzoF). Crystal structure analysis employing angle, buried surface area, and distance measurements showed that incorporation of AzoF at this position causes a closing of the interface by ∼18 ± 3%. This slightly different interface configuration results in a much higher catalytic efficiency in unstimulated HisH due to an elevated turnover number. Moreover, the catalytic efficiency of HisH when stimulated by binding of a substrate to HisF was also significantly increased by AzoF incorporation. This was caused by a K-type stimulation that led to a decrease in the apparent dissociation constant for its substrate, glutamine. In addition, AzoF improved the apparent binding of a substrate analogue at the HisF active site. Remarkably, light-induced isomerization of AzoF considerably enhanced these effects. In conclusion, our findings confirm that signal transduction from HisF to HisH in ImGPS involves the closing of the cyclase:glutaminase subunit interface and that incorporation of AzoF at a hinge position reinforces this catalytically relevant conformational change.

Hands on Methods for High Resolution Cryo-Electron Microscopy Structures of Heterogeneous Macromolecular Complexes.

Electron microscopy of frozen hydrated samples (cryo-EM) is a powerful structural technique that allows the direct study of functional macromolecular complexes in an almost physiological environment. Protein macromolecular complexes are dynamic structures that usually hold together by an intricate network of protein-protein interactions that can be weak and transient. Moreover, a standard feature of many of these complexes is that they behave as nanomachines able to undergo functionally relevant conformational changes in one or several complex components. Among all the other main structural biology techniques, only cryo-EM has the potential of successfully dealing at the same time with both sample heterogeneity and inherent flexibility. The cryo-EM field is currently undergoing a revolution thanks to groundbreaking technical developments that have brought within our reach the possibility of solving the structure of biological complexes at atomic resolution. These technical developments have been mostly focused on new direct electron detector technology and improved sample preparation methods leading to better image quality. This fact has in turn required the development of new and better image processing algorithms to make the most of the higher quality data. The aim of this review is to provide a brief overview of some reported examples of single particle analysis strategies designed to find different conformational and compositional states within target macromolecular complex and specifically to deal with it to reach higher resolution information. Different image processing methodologies specifically aimed to symmetric or pseudo-symmetric protein complexes will also be discussed.

Relative Principal Components Analysis: Application to Analyzing Biomolecular Conformational Changes

A new method termed “Relative Principal Components Analysis” (RPCA) is introduced that extracts optimal relevant principal components to describe the change between two data samples representing two macroscopic states. The method is widely applicable in data-driven science. Calculating the components is based on a physical framework that introduces the objective function (the Kullback–Leibler divergence) appropriate for quantifying the change of the macroscopic state affected by the changes in the microscopic features. To demonstrate the applicability of RPCA, we analyze the thermodynamically relevant conformational changes of the protein HIV-1 protease upon binding to different drug molecules. In this case, the RPCA method provides a sound thermodynamic foundation for analyzing the binding process and thus characterizing both the collective and the locally relevant conformational changes. Moreover, the relevant collective conformational changes can be reconstructed from the informative latent variables to exhibit both the enhanced and the restricted conformational fluctuations upon ligand association.

Chimeric HCN Channels for Studying Camp-Induced Conformational Changes in the C-Linker

Hyperpolarization-activated cyclic nucleotide-gated (HCN1-4) channels are activated by the hyperpolarization of the membrane potential and further regulated by the direct binding of cyclic AMP (cAMP). The binding of cAMP to the cytoplasmic Cyclic Nucleotide Binding Domain (CNBD) increases channel open probability and speeds activation kinetics. cAMP modulation of HCN channels controls heart rate and working memory networks in prefrontal cortex, but also chronic pain at the level of the peripheral nervous system. Thus, a detailed description of the cAMP conformational changes in HCN channels is crucial for understanding several physiological functions, as well as diseases affecting cardiac and neuronal functions. The cAMP-induced conformational changes propagate from the CNBD to the pore through the C-linker region. We and others have previously described, at the atomic level, the conformational changes at the level of the CNBD. The goal of the present work is to describe the propagation of the movements through the C-linker to the pore. To this end, I have constructed a chimeric channel by fusing the C-linker/CNBD of HCN4 to the prokaryotic pore domain of KcsA (hereafter K-H4). Isothermal Titration Calorimetry (ITC) measurements demonstrated that K-H4 is able to bind to cAMP with a KD value of 2μM, which agrees with the one previously published for the isolated C-linker/CNBD fragment of HCN4 and Differential Scanning Calorimetry (DSC) indicated that cAMP causes conformational changes in the chimeric protein. Recently, Double Electron Electron Resonance (DEER) experiments showed that the binding of cAMP at the level of the CNBD domain of K-H4 causes relevant conformational changes in the C-linker region which transits from a dimer of dimers conformation to a 4-fold symmetrical gating ring as suggested by previous functional and biochemical data.

Cooperative Changes in Solvent Exposure Identify Cryptic Pockets, Switches, and Allosteric Coupling

Proteins are dynamic molecules that undergo conformational changes to a broad spectrum of different excited states. Unfortunately, the small populations of these states make it difficult to determine their structures or functional implications. Computer simulations are an increasingly powerful means to identify and characterize functionally relevant excited states. However, this advance has uncovered a further challenge: it can be extremely difficult to identify the most salient features of large simulation data sets. We reasoned that many functionally relevant conformational changes are likely to involve large, cooperative changes to the surfaces that are available to interact with potential binding partners. To examine this hypothesis, we introduce a method that returns a prioritized list of potentially functional conformational changes by segmenting protein structures into clusters of residues that undergo cooperative changes in their solvent exposure, along with the hierarchy of interactions between these groups. We term these groups exposons to distinguish them from other types of clusters that arise in this analysis and others. We demonstrate, using three different model systems, that this method identifies experimentally validated and functionally relevant conformational changes, including conformational switches, allosteric coupling, and cryptic pockets. Our results suggest that key functional sites are hubs in the network of exposons. As a further test of the predictive power of this approach, we apply it to discover cryptic allosteric sites in two different β-lactamase enzymes that are widespread sources of antibiotic resistance. Experimental tests confirm our predictions for both systems. Importantly, we provide the first evidence, to our knowledge, for a cryptic allosteric site in CTX-M-9 β-lactamase. Experimentally testing this prediction did not require any mutations and revealed that this site exerts the most potent allosteric control over activity of any pockets found in β-lactamases to date. Discovery of a similar pocket that was previously overlooked in the well-studied TEM-1 β-lactamase demonstrates the utility of exposons.

Equilibrium free energies from non-equilibrium trajectories with relaxation fluctuation spectroscopy

Recent advances in non-equilibrium statistical mechanics and single-molecule measurements have enabled the determination of equilibrium free energies from non-equilibrium work measurements for fluctuating systems ranging from biological molecules to quantum oscillators. However, for many important non-equilibrium processes, it is difficult or impossible to apply and measure the work required to drive the system through the relevant conformational changes. Here, we show that it is possible, with an appropriate extrapolation to infinite temporal scale and zero spatial scale, to determine equilibrium free energies, without work measurement, by analysing the stochastic trajectories of single biomolecules or other nanoscale, fluctuating systems as they spontaneously relax from a non-equilibrium initial state. We validate the method with simulations and demonstrate its application by determining the free-energy profile for DNA molecules in a structured nanofluidic environment with an experimental protocol that mimics many natural processes with energy injection followed by thermal relaxation.

New insights into the membrane association mechanism of the glycosyltransferase WaaG from Escherichia coli

Monotopic glycosyltransferases (GTs) interact with membranes via electrostatic interactions. The N-terminal domain is permanently anchored to the membrane while the membrane interaction of the C-terminal domain is believed to be weaker so that it undergoes a functionally relevant conformational change upon donor or acceptor binding. Here, we studied the applicability of this model to the glycosyltransferase WaaG. WaaG is involved in the synthesis of lipopolysaccharides (LPS) in Gram-negative bacteria and was previously categorized as a monotopic GT. We analyzed the binding of WaaG to membranes by stopped-flow fluorescence and NMR diffusion experiments. We find that electrostatic interactions are required to bind WaaG to membranes while mere hydrophobic interactions are not sufficient. WaaG senses the membrane's surface charge density but there is no preferential binding to specific anionic lipids. However, the binding is weaker than expected for monotopic GTs but similar to peripheral GTs. Therefore, WaaG may be a peripheral GT and this could be of functional relevance in vivo since LPS synthesis occurs only when WaaG is membrane-bound. We could not observe a C-terminal domain movement under our experimental conditions.

Antibacterial Electrospun Poly(vinyl alcohol)/Enzymatic Synthesized Poly(catechol) Nanofibrous Midlayer Membrane for Ultrafiltration.

Two different nanofibrous antibacterial membranes containing enzymatically synthesized poly(catechol) (PC) or silver nitrate (AgNO3, positive control) blended with poly(vinyl alcohol) (PVA) and electrospun onto a poly(vinylidene fluoride) (PVDF) basal disc to generate thin-film composite midlayers were produced for water ultrafiltration applications. The developed membranes were thoroughly characterized in terms of morphology, chemical composition, and general mechanical and thermal features, antimicrobial activity, and ultrafiltration capabilities. The electrospun blends were recognized as homogeneous. Data revealed relevant conformational changes in the PVA side groups, attributed to hydrogen bonding, high thermal stability, and residual mass. PVDF+PVA/AgNO3 membrane displayed 100% growth inhibition of both Gram-positive and Gram-negative bacteria strains, despite the wide range of fiber diameters generated, from 24 to 125 nm, formation of numerous beads, and irregular morphology. The PVDF+PVA/PC membrane showed a good growth inhibition of Staphylococcus aureus (92%) and revealed a smooth morphology with no relevant bead formations and diameters ranging from 68 to 131 nm. The ultrafiltration abilities of the membrane containing PVA/PC were tested in a dead-end high-pressure cell (4 bar) using a reactive dye in distilled water and seawater. After 5 cycles, a maximum rejection of ≈85% with an average flux rate of 70 L m-2 h-1 for distilled water and ≈64% with an average flux rate of 62 L m-2 h-1 for seawater were determined with an overall salt rejection of ≈5%.

Reconstructed ancestral enzymes reveal that negative selection drove the evolution of substrate specificity in ADP-dependent kinases

One central goal in molecular evolution is to pinpoint the mechanisms and evolutionary forces that cause an enzyme to change its substrate specificity; however, these processes remain largely unexplored. Using the glycolytic ADP-dependent kinases of archaea, including the orders Thermococcales, Methanosarcinales, and Methanococcales, as a model and employing an approach involving paleoenzymology, evolutionary statistics, and protein structural analysis, we could track changes in substrate specificity during ADP-dependent kinase evolution along with the structural determinants of these changes. To do so, we studied five key resurrected ancestral enzymes as well as their extant counterparts. We found that a major shift in function from a bifunctional ancestor that could phosphorylate either glucose or fructose 6-phosphate (fructose-6-P) as a substrate to a fructose 6-P-specific enzyme was started by a single amino acid substitution resulting in negative selection with a ground-state mode against glucose and a subsequent 1,600-fold change in specificity of the ancestral protein. This change rendered the residual phosphorylation of glucose a promiscuous and physiologically irrelevant activity, highlighting how promiscuity may be an evolutionary vestige of ancestral enzyme activities, which have been eliminated over time. We also could reconstruct the evolutionary history of substrate utilization by using an evolutionary model of discrete binary characters, indicating that substrate uses can be discretely lost or acquired during enzyme evolution. These findings exemplify how negative selection and subtle enzyme changes can lead to major evolutionary shifts in function, which can subsequently generate important adaptive advantages, for example, in improving glycolytic efficiency in Thermococcales.