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.

[1]  Purificación López-García,et al.  Extending the conserved phylogenetic core of archaea disentangles the evolution of the third domain of life. , 2015, Molecular biology and evolution.

[2]  David S. Goodsell,et al.  AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility , 2009, J. Comput. Chem..

[3]  Randy J. Read,et al.  Phaser crystallographic software , 2007, Journal of applied crystallography.

[4]  V. Arcus,et al.  Rapid bursts and slow declines: on the possible evolutionary trajectories of enzymes , 2015, Journal of The Royal Society Interface.

[5]  V. Arcus,et al.  On the origin and evolution of thermophily: reconstruction of functional precambrian enzymes from ancestors of Bacillus. , 2012, Molecular biology and evolution.

[6]  Kenneth H. Wolfe,et al.  Turning a hobby into a job: How duplicated genes find new functions , 2008, Nature Reviews Genetics.

[7]  F. J. Fernández,et al.  Crystal Structure, SAXS and Kinetic Mechanism of Hyperthermophilic ADP-Dependent Glucokinase from Thermococcus litoralis Reveal a Conserved Mechanism for Catalysis , 2013, PloS one.

[8]  David R. Liu,et al.  Negative selection and stringency modulation in phage-assisted constinuous evolution , 2014, Nature chemical biology.

[9]  Arthur J. Olson,et al.  AutoDock Vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading , 2009, J. Comput. Chem..

[10]  H. Matsuzawa,et al.  Structural basis for the ADP-specificity of a novel glucokinase from a hyperthermophilic archaeon. , 2001, Structure.

[11]  V. Castro-Fernández,et al.  Bifunctional ADP‐dependent phosphofructokinase/glucokinase activity in the order Methanococcales – biochemical characterization of the mesophilic enzyme from Methanococcus maripaludis , 2014, The FEBS journal.

[12]  P. Griffin,et al.  Biophysical mechanisms for large-effect mutations in the evolution of steroid hormone receptors , 2013, Proceedings of the National Academy of Sciences.

[13]  D. Posada jModelTest: phylogenetic model averaging. , 2008, Molecular biology and evolution.

[14]  M. Sippl Recognition of errors in three‐dimensional structures of proteins , 1993, Proteins.

[15]  D. Eisenberg,et al.  A method to identify protein sequences that fold into a known three-dimensional structure. , 1991, Science.

[16]  Vincent B. Chen,et al.  Correspondence e-mail: , 2000 .

[17]  M. Pagel The Maximum Likelihood Approach to Reconstructing Ancestral Character States of Discrete Characters on Phylogenies , 1999 .

[18]  Kevin Cowtan,et al.  research papers Acta Crystallographica Section D Biological , 2005 .

[19]  Anders Gorm Pedersen,et al.  RevTrans: multiple alignment of coding DNA from aligned amino acid sequences , 2003, Nucleic Acids Res..

[20]  H. Sakuraba,et al.  ADP-dependent Glucokinase/Phosphofructokinase, a Novel Bifunctional Enzyme from the Hyperthermophilic ArchaeonMethanococcus jannaschii * , 2002, The Journal of Biological Chemistry.

[21]  H. Sakuraba,et al.  Biochemical characterization, cloning, and sequencing of ADP-dependent (AMP-forming) glucokinase from two hyperthermophilic archaea, Pyrococcus furiosus and Thermococcus litoralis. , 2000, Journal of biochemistry.

[22]  Dan S. Tawfik,et al.  The moderately efficient enzyme: evolutionary and physicochemical trends shaping enzyme parameters. , 2011, Biochemistry.

[23]  Annabel E. Todd,et al.  Evolution of function in protein superfamilies, from a structural perspective. , 2001, Journal of molecular biology.

[24]  N. Allewell,et al.  Thematic Minireview Series on Enzyme Evolution in the Post-genomic Era , 2011, The Journal of Biological Chemistry.

[25]  F. Merino,et al.  Specificity evolution of the ADP‐dependent sugar kinase family –in silico studies of the glucokinase/phosphofructokinase bifunctional enzyme from Methanocaldococcus jannaschii , 2008, The FEBS journal.

[26]  Ziheng Yang PAML 4: phylogenetic analysis by maximum likelihood. , 2007, Molecular biology and evolution.

[27]  Dan S. Tawfik,et al.  Enzyme promiscuity: a mechanistic and evolutionary perspective. , 2010, Annual review of biochemistry.

[28]  S. Fushinobu,et al.  Crystal structure of an ADP-dependent glucokinase from Pyrococcus furiosus: implications for a sugar-induced conformational change in ADP-dependent kinase. , 2003, Journal of molecular biology.

[29]  John P. Huelsenbeck,et al.  MrBayes 3: Bayesian phylogenetic inference under mixed models , 2003, Bioinform..

[30]  N. Grishin,et al.  PROMALS3D: a tool for multiple protein sequence and structure alignments , 2008, Nucleic acids research.

[31]  M. Pagel,et al.  Bayesian Analysis of Correlated Evolution of Discrete Characters by Reversible‐Jump Markov Chain Monte Carlo , 2006, The American Naturalist.

[32]  B. Hall,et al.  Simple and accurate estimation of ancestral protein sequences. , 2006, Proceedings of the National Academy of Sciences of the United States of America.

[33]  F. Merino,et al.  Catalytic and regulatory roles of divalent metal cations on the phosphoryl-transfer mechanism of ADP-dependent sugar kinases from hyperthermophilic archaea. , 2012, Biochimie.

[34]  Jonathan P. Bollback,et al.  Empirical and hierarchical Bayesian estimation of ancestral states. , 2001, Systematic biology.

[35]  Expression and role in glycolysis of human ADP-dependent glucokinase , 2012, Molecular and Cellular Biochemistry.

[36]  Philip R. Evans,et al.  How good are my data and what is the resolution? , 2013, Acta crystallographica. Section D, Biological crystallography.

[37]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[38]  M. Pagel,et al.  Bayesian estimation of ancestral character states on phylogenies. , 2004, Systematic biology.

[39]  Laxmikant V. Kalé,et al.  Scalable molecular dynamics with NAMD , 2005, J. Comput. Chem..

[40]  Adi Stern,et al.  Assessing the prediction fidelity of ancestral reconstruction by a library approach. , 2015, Protein engineering, design & selection : PEDS.

[41]  Dan S. Tawfik,et al.  Accuracy-rate tradeoffs: how do enzymes meet demands of selectivity and catalytic efficiency? , 2014, Current opinion in chemical biology.

[42]  Randy J. Read,et al.  Acta Crystallographica Section D Biological , 2003 .

[43]  Felipe Merino,et al.  The ADP‐dependent sugar kinase family: Kinetic and evolutionary aspects , 2009, IUBMB life.

[44]  Z. Jia,et al.  ADP-dependent 6-Phosphofructokinase from Pyrococcus horikoshii OT3 , 2009, The Journal of Biological Chemistry.

[45]  Zaida Luthey-Schulten,et al.  MultiSeq: unifying sequence and structure data for evolutionary analysis , 2006, BMC Bioinformatics.

[46]  H. A. Orr,et al.  The genetic theory of adaptation: a brief history , 2005, Nature Reviews Genetics.

[47]  J. Thornton,et al.  PROCHECK: a program to check the stereochemical quality of protein structures , 1993 .

[48]  Bryan D. Kolaczkowski,et al.  Robustness of Ancestral Sequence Reconstruction to Phylogenetic Uncertainty , 2010, Molecular biology and evolution.

[49]  Lenwood S. Heath,et al.  H++: a server for estimating pKas and adding missing hydrogens to macromolecules , 2005, Nucleic Acids Res..

[50]  R. MacLean,et al.  Fitness Is Strongly Influenced by Rare Mutations of Large Effect in a Microbial Mutation Accumulation Experiment , 2014, Genetics.

[51]  K. Verstrepen,et al.  Reconstruction of Ancestral Metabolic Enzymes Reveals Molecular Mechanisms Underlying Evolutionary Innovation through Gene Duplication , 2012, PLoS biology.

[52]  E. Rocha,et al.  Horizontal Transfer, Not Duplication, Drives the Expansion of Protein Families in Prokaryotes , 2011, PLoS genetics.

[53]  M. Nei,et al.  A new method of inference of ancestral nucleotide and amino acid sequences. , 1995, Genetics.

[54]  T. Blundell,et al.  Comparative protein modelling by satisfaction of spatial restraints. , 1993, Journal of molecular biology.

[55]  M. Harms,et al.  Analyzing protein structure and function using ancestral gene reconstruction. , 2010, Current opinion in structural biology.

[56]  S. Benner,et al.  Molecular paleoscience: systems biology from the past. , 2007, Advances in enzymology and related areas of molecular biology.

[57]  W. D. de Vos,et al.  Molecular and Biochemical Characterization of the ADP-dependent Phosphofructokinase from the Hyperthermophilic Archaeon Pyrococcus furiosus * , 1999, The Journal of Biological Chemistry.

[58]  R. Moreno-Sánchez,et al.  The nutritional status of Methanosarcina acetivorans regulates glycogen metabolism and gluconeogenesis and glycolysis fluxes , 2016, The FEBS journal.

[59]  W. D. de Vos,et al.  Evidence for the operation of a novel Embden-Meyerhof pathway that involves ADP-dependent kinases during sugar fermentation by Pyrococcus furiosus. , 1994, The Journal of biological chemistry.