{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,2,27]],"date-time":"2026-02-27T04:11:10Z","timestamp":1772165470312,"version":"3.50.1"},"reference-count":44,"publisher":"Springer Science and Business Media LLC","issue":"1","license":[{"start":{"date-parts":[[2024,7,18]],"date-time":"2024-07-18T00:00:00Z","timestamp":1721260800000},"content-version":"tdm","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0"},{"start":{"date-parts":[[2024,7,18]],"date-time":"2024-07-18T00:00:00Z","timestamp":1721260800000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0"}],"content-domain":{"domain":["link.springer.com"],"crossmark-restriction":false},"short-container-title":["BMC Bioinformatics"],"abstract":"<jats:title>Abstract<\/jats:title>\n                  <jats:sec>\n                    <jats:title>Background<\/jats:title>\n                    <jats:p>Metabolic pathways support the enzyme flux that converts input chemicals into energy and cellular building blocks. With a constant rate of input, steady-state flux is achieved when metabolite concentrations and reaction rates remain constant over time. Individual genes undergo mutation, while selection acts on higher level functions of the pathway, such as steady-state flux where applicable. Modeling the evolution of metabolic pathways through mechanistic sets of ordinary differential equations is a piece of the genotype\u2013phenotype map model for interpreting genetic variation and inter-specific differences. Such models can generate distinct compensatory changes and adaptive changes from directional selection, indicating single nucleotide polymorphisms and fixed differences that could affect phenotype. If used for inference, this would ultimately enable detection of selection on metabolic pathways as well as inference of ancestral states for metabolic pathway function.<\/jats:p>\n                  <\/jats:sec>\n                  <jats:sec>\n                    <jats:title>Results<\/jats:title>\n                    <jats:p>A software tool for simulating the evolution of metabolic pathways based upon underlying biochemistry, phylogenetics, and evolutionary considerations is presented. The Python program, Phylogenetic Evolution of Metabolic Pathway Simulator (PEMPS), implements a mutation-selection framework to simulate the evolution of the pathway over a phylogeny by interfacing with COPASI to calculate the steady-state flux of the metabolic network, introducing mutations as alterations in parameter values according to a model, and calculating a fitness score and corresponding probability of fixation based on the change in steady-state flux value(s). Results from simulations are consistent with a priori expectations of fixation probabilities and systematic change in model parameters.<\/jats:p>\n                  <\/jats:sec>\n                  <jats:sec>\n                    <jats:title>Conclusions<\/jats:title>\n                    <jats:p>\n                      The PEMPS program simulates the evolution of a metabolic pathway with a mutation-selection modeling framework based on criteria like steady-state flux that is designed to work with SBML-formatted kinetic models, and Newick-formatted phylogenetic trees. The Python software is run on the Linux command line and is available at\n                      <jats:ext-link xmlns:xlink=\"http:\/\/www.w3.org\/1999\/xlink\" ext-link-type=\"uri\" xlink:href=\"https:\/\/github.com\/nmccloskey\/PEMPS\">https:\/\/github.com\/nmccloskey\/PEMPS<\/jats:ext-link>\n                      .\n                    <\/jats:p>\n                  <\/jats:sec>","DOI":"10.1186\/s12859-024-05867-7","type":"journal-article","created":{"date-parts":[[2024,7,18]],"date-time":"2024-07-18T12:03:00Z","timestamp":1721304180000},"update-policy":"https:\/\/doi.org\/10.1007\/springer_crossmark_policy","source":"Crossref","is-referenced-by-count":0,"title":["PEMPS: a phylogenetic software tool to model the evolution of metabolic pathways"],"prefix":"10.1186","volume":"25","author":[{"given":"Nicholas S.","family":"McCloskey","sequence":"first","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Ayna","family":"Mammedova","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"David A.","family":"Liberles","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"297","published-online":{"date-parts":[[2024,7,18]]},"reference":[{"key":"5867_CR1","doi-asserted-by":"publisher","first-page":"697","DOI":"10.1016\/j.biochi.2009.01.004","volume":"91","author":"NJ Kruger","year":"2009","unstructured":"Kruger NJ, Ratcliffe RG. Insights into plant metabolic networks from steady-state metabolic flux analysis. Biochimie. 2009;91:697\u2013702.","journal-title":"Biochimie"},{"key":"5867_CR2","doi-asserted-by":"publisher","first-page":"31","DOI":"10.1186\/s13062-016-0133-6","volume":"11","author":"A Orlenko","year":"2016","unstructured":"Orlenko A, Teufel AI, Chi PB, Liberles DA. Selection on metabolic pathway function in the presence of mutation-selection-drift balance leads to rate-limiting steps that are not evolutionarily stable. Biol Direct. 2016;11:31.","journal-title":"Biol Direct"},{"key":"5867_CR3","doi-asserted-by":"publisher","first-page":"1763","DOI":"10.1104\/pp.109.149237","volume":"152","author":"K Schallau","year":"2010","unstructured":"Schallau K, Junker BH. Simulating plant metabolic pathways with enzyme-kinetic models. Plant Physiol. 2010;152:1763\u201371.","journal-title":"Plant Physiol"},{"key":"5867_CR4","doi-asserted-by":"publisher","first-page":"146","DOI":"10.1007\/s00239-016-9731-2","volume":"82","author":"A Orlenko","year":"2016","unstructured":"Orlenko A, Hermansen RA, Liberles DA. Flux control in glycolysis varies across the tree of life. J Mol Evol. 2016;82:146\u201361.","journal-title":"J Mol Evol"},{"key":"5867_CR5","doi-asserted-by":"publisher","first-page":"410","DOI":"10.1016\/j.tree.2010.04.001","volume":"25","author":"K Csill\u00e9ry","year":"2010","unstructured":"Csill\u00e9ry K, Blum MGB, Gaggiotti OE, Fran\u00e7ois O. Approximate Bayesian computation (ABC) in practice. Trends Ecol Evol. 2010;25:410\u20138.","journal-title":"Trends Ecol Evol"},{"key":"5867_CR6","doi-asserted-by":"publisher","first-page":"D587","DOI":"10.1093\/nar\/gkac963","volume":"51","author":"M Kanehisa","year":"2023","unstructured":"Kanehisa M, Furumichi M, Sato Y, Kawashima M, Ishiguro-Watanabe M. KEGG for taxonomy-based analysis of pathways and genomes. Nucleic Acids Res. 2023;51:D587\u201392.","journal-title":"Nucleic Acids Res"},{"key":"5867_CR7","doi-asserted-by":"publisher","first-page":"1528","DOI":"10.1093\/bioinformatics\/btp262","volume":"25","author":"A Mithani","year":"2009","unstructured":"Mithani A, Preston GM, Hein J. A stochastic model for the evolution of metabolic networks with neighbor dependence. Bioinformatics. 2009;25:1528\u201335.","journal-title":"Bioinformatics"},{"key":"5867_CR8","doi-asserted-by":"publisher","first-page":"479","DOI":"10.1038\/msb.2011.11","volume":"7","author":"K Yizhak","year":"2011","unstructured":"Yizhak K, Tuller T, Papp B, Ruppin E. Metabolic modeling of endosymbiont genome reduction on a temporal scale. Mol Syst Biol. 2011;7:479.","journal-title":"Mol Syst Biol"},{"key":"5867_CR9","first-page":"mgen000378","volume":"6","author":"RJ Hall","year":"2020","unstructured":"Hall RJ, Thorpe S, Thomas GH, Wood AJ. Simulating the evolutionary trajectories of metabolic pathways for insect symbionts in the genus Sodalis. Microb Genom. 2020;6:mgen000378.","journal-title":"Microb Genom"},{"key":"5867_CR10","doi-asserted-by":"publisher","first-page":"e1000868","DOI":"10.1371\/journal.pcbi.1000868","volume":"6","author":"A Mithani","year":"2010","unstructured":"Mithani A, Preston GM, Hein J. A Bayesian approach to the evolution of metabolic networks on a phylogeny. PLoS Comput Biol. 2010;6:e1000868.","journal-title":"PLoS Comput Biol"},{"key":"5867_CR11","doi-asserted-by":"publisher","first-page":"483","DOI":"10.1534\/genetics.109.110411","volume":"184","author":"KM Wright","year":"2010","unstructured":"Wright KM, Rausher MD. The evolution of control and distribution of adaptive mutations in a metabolic pathway. Genetics. 2010;184:483\u2013502.","journal-title":"Genetics"},{"key":"5867_CR12","doi-asserted-by":"publisher","first-page":"117","DOI":"10.1186\/s12862-017-0962-7","volume":"17","author":"A Orlenko","year":"2017","unstructured":"Orlenko A, Chi PB, Liberles DA. Characterizing the roles of changing population size and selection on the evolution of flux control in metabolic pathways. BMC Evol Biol. 2017;17:117.","journal-title":"BMC Evol Biol"},{"key":"5867_CR13","doi-asserted-by":"publisher","first-page":"232","DOI":"10.1186\/s12862-015-0515-x","volume":"15","author":"RA Hermansen","year":"2015","unstructured":"Hermansen RA, Mannakee BK, Knecht W, Liberles DA, Gutenkunst RN. Characterizing selective pressures on the pathway for de novo biosynthesis of pyrimidines in yeast. BMC Evol Biol. 2015;15:232.","journal-title":"BMC Evol Biol"},{"key":"5867_CR14","doi-asserted-by":"publisher","first-page":"707","DOI":"10.1093\/molbev\/msv263","volume":"33","author":"X Zhao","year":"2016","unstructured":"Zhao X, Bergland AO, Behrman EL, Gregory BD, Petrov DA, Schmidt PS. Global transcriptional profiling of diapause and climatic adaptation in Drosophila melanogaster. Mol Biol Evol. 2016;33:707\u201320.","journal-title":"Mol Biol Evol"},{"key":"5867_CR15","doi-asserted-by":"publisher","first-page":"2","DOI":"10.1007\/s00239-022-10084-y","volume":"91","author":"DA Liberles","year":"2023","unstructured":"Liberles DA. The memory problem for neutral mutational models of evolution. J Mol Evol. 2023;91:2\u20135.","journal-title":"J Mol Evol"},{"key":"5867_CR16","doi-asserted-by":"publisher","first-page":"3","DOI":"10.1016\/j.pisc.2014.12.002","volume":"4","author":"A Cornish-Bowden","year":"2015","unstructured":"Cornish-Bowden A. One hundred years of Michaelis\u2013Menten kinetics. Perspect Sci (Neth). 2015;4:3\u20139.","journal-title":"Perspect Sci (Neth)"},{"key":"5867_CR17","doi-asserted-by":"publisher","first-page":"1168","DOI":"10.1002\/pro.2886","volume":"25","author":"PB Chi","year":"2016","unstructured":"Chi PB, Liberles DA. Selection on protein structure, interaction, and sequence. Protein Sci. 2016;25:1168\u201378.","journal-title":"Protein Sci"},{"key":"5867_CR18","doi-asserted-by":"publisher","first-page":"1569","DOI":"10.1093\/bioinformatics\/btq228","volume":"26","author":"J Sukumaran","year":"2010","unstructured":"Sukumaran J, Holder MT. DendroPy: a Python library for phylogenetic computing. Bioinformatics. 2010;26:1569\u201371.","journal-title":"Bioinformatics"},{"key":"5867_CR19","doi-asserted-by":"publisher","first-page":"baaa062","DOI":"10.1093\/database\/baaa062","volume":"2020","author":"CL Schoch","year":"2020","unstructured":"Schoch CL, Ciufo S, Domrachev M, Hotton CL, Kannan S, Khovanskaya R, et al. NCBI taxonomy: a comprehensive update on curation, resources and tools. Database. 2020;2020:baaa062.","journal-title":"Database"},{"key":"5867_CR20","doi-asserted-by":"publisher","first-page":"msac174","DOI":"10.1093\/molbev\/msac174","volume":"39","author":"S Kumar","year":"2022","unstructured":"Kumar S, Suleski M, Craig JM, Kasprowicz AE, Sanderford M, Li M, et al. TimeTree 5: an expanded resource for species divergence times. Mol Biol Evol. 2022;39:msac174.","journal-title":"Mol Biol Evol."},{"key":"5867_CR21","doi-asserted-by":"publisher","first-page":"240","DOI":"10.1007\/s00239-005-0096-1","volume":"63","author":"A-C Berglund-Sonnhammer","year":"2006","unstructured":"Berglund-Sonnhammer A-C, Steffansson P, Betts MJ, Liberles DA. Optimal gene trees from sequences and species trees using a soft interpretation of parsimony. J Mol Evol. 2006;63:240\u201350.","journal-title":"J Mol Evol"},{"key":"5867_CR22","doi-asserted-by":"publisher","first-page":"429","DOI":"10.1089\/106652700750050871","volume":"7","author":"K Chen","year":"2000","unstructured":"Chen K, Durand D, Farach-Colton M. NOTUNG: a program for dating gene duplications and optimizing gene family trees. J Comput Biol. 2000;7:429\u201347.","journal-title":"J Comput Biol"},{"key":"5867_CR23","doi-asserted-by":"publisher","first-page":"215","DOI":"10.1016\/j.jbiotec.2017.06.1200","volume":"261","author":"FT Bergmann","year":"2017","unstructured":"Bergmann FT, Hoops S, Klahn B, Kummer U, Mendes P, Pahle J, et al. COPASI and its applications in biotechnology. J Biotechnol. 2017;261:215\u201320.","journal-title":"J Biotechnol"},{"key":"5867_CR24","doi-asserted-by":"publisher","first-page":"3067","DOI":"10.1093\/bioinformatics\/btl485","volume":"22","author":"S Hoops","year":"2006","unstructured":"Hoops S, Sahle S, Gauges R, Lee C, Pahle J, Simus N, et al. COPASI\u2014a complex pathway simulator. Bioinformatics. 2006;22:3067\u201374.","journal-title":"Bioinformatics"},{"key":"5867_CR25","doi-asserted-by":"publisher","first-page":"5553","DOI":"10.21105\/joss.05553","volume":"8","author":"FT Bergmann","year":"2023","unstructured":"Bergmann FT. BASICO: a simplified python interface to COPASI. J Open Source Softw. 2023;8:5553.","journal-title":"J Open Source Softw"},{"key":"5867_CR26","doi-asserted-by":"publisher","first-page":"409","DOI":"10.3390\/genes9080409","volume":"9","author":"A Teufel","year":"2018","unstructured":"Teufel A, Ritchie A, Wilke C, Liberles D. Using the mutation-selection framework to characterize selection on protein sequences. Genes (Basel). 2018;9:409.","journal-title":"Genes (Basel)"},{"key":"5867_CR27","doi-asserted-by":"publisher","first-page":"1101","DOI":"10.1534\/genetics.111.136432","volume":"190","author":"AU Tamuri","year":"2012","unstructured":"Tamuri AU, dos Reis M, Goldstein RA. Estimating the distribution of selection coefficients from phylogenetic data using sitewise mutation-selection models. Genetics. 2012;190:1101\u201315.","journal-title":"Genetics"},{"key":"5867_CR28","doi-asserted-by":"publisher","first-page":"713","DOI":"10.1093\/genetics\/47.6.713","volume":"47","author":"M Kimura","year":"1962","unstructured":"Kimura M. On the probability of fixation of mutant genes in a population. Genetics. 1962;47:713\u20139.","journal-title":"Genetics"},{"key":"5867_CR29","doi-asserted-by":"publisher","unstructured":"Otto SP, Whitlock MC. Fixation probabilities and times. In: eLS. John Wiley & Sons, Ltd: Chichester; 2013. https:\/\/doi.org\/10.1002\/9780470015902.a0005464.pub3.","DOI":"10.1002\/9780470015902.a0005464.pub3"},{"key":"5867_CR30","doi-asserted-by":"publisher","first-page":"20160906","DOI":"10.1098\/rsif.2016.0906","volume":"14","author":"AI Teufel","year":"2017","unstructured":"Teufel AI, Wilke CO. Accelerated simulation of evolutionary trajectories in origin-fixation models. J R Soc Interface. 2017;14:20160906.","journal-title":"J R Soc Interface"},{"key":"5867_CR31","first-page":"1","volume":"121","author":"IJ Lovette","year":"2004","unstructured":"Lovette IJ. Mitochondrial dating and mixed support for the \u201c2% rule\u201d in birds. Auk. 2004;121:1\u20136.","journal-title":"Auk"},{"key":"5867_CR32","doi-asserted-by":"publisher","first-page":"68","DOI":"10.1186\/s12885-022-09178-z","volume":"22","author":"T Huang","year":"2022","unstructured":"Huang T, Li J, Wang SM. Core promoter mutation contributes to abnormal gene expression in bladder cancer. BMC Cancer. 2022;22:68.","journal-title":"BMC Cancer"},{"issue":"Pt 4","key":"5867_CR33","doi-asserted-by":"publisher","first-page":"1003","DOI":"10.1099\/00221287-148-4-1003","volume":"148","author":"MHN Hoefnagel","year":"2002","unstructured":"Hoefnagel MHN, et al. Metabolic engineering of lactic acid bacteria, the combined approach: kinetic modelling, metabolic control and experimental analysis. Microbiology. 2002;148(Pt 4):1003\u201313.","journal-title":"Microbiology"},{"key":"5867_CR34","doi-asserted-by":"publisher","first-page":"e1003371","DOI":"10.1371\/journal.pcbi.1003371","volume":"9","author":"EJ Kerkhoven","year":"2013","unstructured":"Kerkhoven EJ, Achcar F, Alibu VP, Burchmore RJ, Gilbert IH, Trybi\u0142o M, et al. Handling uncertainty in dynamic models: the pentose phosphate pathway in Trypanosoma brucei. PLoS Comput Biol. 2013;9:e1003371.","journal-title":"PLoS Comput Biol"},{"key":"5867_CR35","doi-asserted-by":"publisher","first-page":"e1002483","DOI":"10.1371\/journal.pcbi.1002483","volume":"8","author":"K van Eunen","year":"2012","unstructured":"van Eunen K, Kiewiet JAL, Westerhoff HV, Bakker BM. Testing biochemistry revisited: how in vivo metabolism can be understood from in vitro enzyme kinetics. PLoS Comput Biol. 2012;8:e1002483.","journal-title":"PLoS Comput Biol"},{"key":"5867_CR36","doi-asserted-by":"publisher","first-page":"e1006132","DOI":"10.1371\/journal.pgen.1006132","volume":"12","author":"BK Mannakee","year":"2016","unstructured":"Mannakee BK, Gutenkunst RN. Selection on network dynamics drives differential rates of protein domain evolution. PLoS Genet. 2016;12:e1006132.","journal-title":"PLoS Genet"},{"key":"5867_CR37","doi-asserted-by":"publisher","first-page":"E2821","DOI":"10.1073\/pnas.1310980110","volume":"110","author":"M Lynch","year":"2013","unstructured":"Lynch M. Evolutionary diversification of the multimeric states of proteins. Proc Natl Acad Sci. 2013;110:E2821\u20138.","journal-title":"Proc Natl Acad Sci"},{"key":"5867_CR38","doi-asserted-by":"publisher","first-page":"436","DOI":"10.1038\/msb.2010.88","volume":"6","author":"J Kim","year":"2010","unstructured":"Kim J, Kershner JP, Novikov Y, Shoemaker RK, Copley SD. Three serendipitous pathways in E. coli can bypass a block in pyridoxal-5\u2032-phosphate synthesis. Mol Syst Biol. 2010;6:436.","journal-title":"Mol Syst Biol"},{"key":"5867_CR39","doi-asserted-by":"publisher","first-page":"e1838","DOI":"10.1371\/journal.pone.0001838","volume":"3","author":"DJ Kliebenstein","year":"2008","unstructured":"Kliebenstein DJ. A role for gene duplication and natural variation of gene expression in the evolution of metabolism. PLoS ONE. 2008;3:e1838.","journal-title":"PLoS ONE"},{"key":"5867_CR40","doi-asserted-by":"publisher","first-page":"28","DOI":"10.1016\/j.cell.2020.09.014","volume":"183","author":"T Misteli","year":"2020","unstructured":"Misteli T. The self-organizing genome: principles of genome architecture and function. Cell. 2020;183:28\u201345.","journal-title":"Cell"},{"key":"5867_CR41","doi-asserted-by":"publisher","first-page":"687710","DOI":"10.1007\/s00239-023-10127-y","volume":"91","author":"MA Savageau","year":"2023","unstructured":"Savageau MA. Phenotype design space provides a mechanistic framework relating molecular parameters to phenotype diversity available for selection. J Mol Evol. 2023;91:687710.","journal-title":"J Mol Evol"},{"key":"5867_CR42","doi-asserted-by":"publisher","first-page":"534","DOI":"10.1002\/bit.20558","volume":"91","author":"JJ Heijnen","year":"2005","unstructured":"Heijnen JJ. Approximative kinetic formats used in metabolic network modeling. Biotechnol Bioeng. 2005;91:534\u201345.","journal-title":"Biotechnol Bioeng"},{"key":"5867_CR43","doi-asserted-by":"crossref","unstructured":"Gillespie DT. Stochastic chemical kinetics. In: Handbook of materials modeling. Springer, Dordrecht; 2005. p. 1735\u201352.","DOI":"10.1007\/1-4020-3286-2_87"},{"key":"5867_CR44","doi-asserted-by":"publisher","first-page":"2340","DOI":"10.1021\/j100540a008","volume":"81","author":"DT Gillespie","year":"1977","unstructured":"Gillespie DT. Exact stochastic simulation of coupled chemical reactions. J Phys Chem. 1977;81:2340\u201361.","journal-title":"J Phys Chem"}],"container-title":["BMC Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/link.springer.com\/content\/pdf\/10.1186\/s12859-024-05867-7.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"text-mining"},{"URL":"https:\/\/link.springer.com\/article\/10.1186\/s12859-024-05867-7\/fulltext.html","content-type":"text\/html","content-version":"vor","intended-application":"text-mining"},{"URL":"https:\/\/link.springer.com\/content\/pdf\/10.1186\/s12859-024-05867-7.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2024,7,18]],"date-time":"2024-07-18T13:01:56Z","timestamp":1721307716000},"score":1,"resource":{"primary":{"URL":"https:\/\/bmcbioinformatics.biomedcentral.com\/articles\/10.1186\/s12859-024-05867-7"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2024,7,18]]},"references-count":44,"journal-issue":{"issue":"1","published-online":{"date-parts":[[2024,12]]}},"alternative-id":["5867"],"URL":"https:\/\/doi.org\/10.1186\/s12859-024-05867-7","relation":{"has-preprint":[{"id-type":"doi","id":"10.1101\/2024.01.04.574206","asserted-by":"object"},{"id-type":"doi","id":"10.21203\/rs.3.rs-3793038\/v1","asserted-by":"object"}]},"ISSN":["1471-2105"],"issn-type":[{"value":"1471-2105","type":"electronic"}],"subject":[],"published":{"date-parts":[[2024,7,18]]},"assertion":[{"value":"22 December 2023","order":1,"name":"received","label":"Received","group":{"name":"ArticleHistory","label":"Article History"}},{"value":"11 July 2024","order":2,"name":"accepted","label":"Accepted","group":{"name":"ArticleHistory","label":"Article History"}},{"value":"18 July 2024","order":3,"name":"first_online","label":"First Online","group":{"name":"ArticleHistory","label":"Article History"}},{"order":1,"name":"Ethics","group":{"name":"EthicsHeading","label":"Declarations"}},{"value":"Not applicable.","order":2,"name":"Ethics","group":{"name":"EthicsHeading","label":"Ethics approval and consent to participate"}},{"value":"Not applicable.","order":3,"name":"Ethics","group":{"name":"EthicsHeading","label":"Consent for publication"}},{"value":"The authors declare no competing interests.","order":4,"name":"Ethics","group":{"name":"EthicsHeading","label":"Competing interests"}}],"article-number":"244"}}