{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2025,7,30]],"date-time":"2025-07-30T11:44:54Z","timestamp":1753875894932,"version":"3.41.2"},"reference-count":45,"publisher":"Oxford University Press (OUP)","issue":"6","license":[{"start":{"date-parts":[[2025,5,6]],"date-time":"2025-05-06T00:00:00Z","timestamp":1746489600000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"funder":[{"name":"European Union\u2019s Horizon 2020 Programme"}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2025,6,2]]},"abstract":"<jats:title>Abstract<\/jats:title>\n               <jats:sec>\n                  <jats:title>Motivation<\/jats:title>\n                  <jats:p>Sensitivity analysis is a useful tool to identify key parameters in metabolic models. It is typically only applied to the growth rate, disregarding the sensitivity of other solution variables to parameters. Further, sensitivity analysis of elementary flux modes could provide low-dimensional insights into optimal solutions, but they are not defined when a model is subject to inhomogeneous flux constraints, such as the frequently used ATP maintenance reaction.<\/jats:p>\n               <\/jats:sec>\n               <jats:sec>\n                  <jats:title>Results<\/jats:title>\n                  <jats:p>We introduce optimal flux modes (OFMs), an analogue to elementary flux modes (EFMs), but specifically applied to optimal solutions of constraint-based models. Further, we prove that implicit differentiation can always be used to efficiently calculate the sensitivities of both whole-model solutions and OFM-based solutions to model parameters. This allows for fine-grained sensitivity analysis of the optimal solution, and investigation of how these parameters exert control on the optimal composition of OFMs. This novel framework is implemented in DifferentiableMetabolism.jl, a software package designed to efficiently differentiate solutions of constraint-based models. To demonstrate scalability, we differentiate solutions of 342 yeast models; additionally we show that sensitivities of specific subsystems can guide metabolic engineering. Applying our scheme to an Escherichia coli model, we find that OFM sensitivities predict the effect of knockout experiments on waste product accumulation. Sensitivity analysis of OFMs also provides key insights into metabolic changes resulting from parameter perturbations.<\/jats:p>\n               <\/jats:sec>\n               <jats:sec>\n                  <jats:title>Availability and implementation<\/jats:title>\n                  <jats:p>Software introduced here is available as open-source Julia packages DifferentiableMetabolism.jl (https:\/\/github.com\/stelmo\/DifferentiableMetabolism.jl) and ElementaryFluxModes.jl (https:\/\/github.com\/HettieC\/ElementaryFluxModes.jl), which both work on all major operating systems and computer architectures. Code to reproduce all results is available from https:\/\/github.com\/HettieC\/DifferentiableOFMPaper, and as an archive from https:\/\/doi.org\/10.5281\/zenodo.15183208.<\/jats:p>\n               <\/jats:sec>","DOI":"10.1093\/bioinformatics\/btaf287","type":"journal-article","created":{"date-parts":[[2025,5,6]],"date-time":"2025-05-06T16:25:06Z","timestamp":1746548706000},"source":"Crossref","is-referenced-by-count":0,"title":["Algebraic differentiation for fast sensitivity analysis of optimal flux modes in metabolic models"],"prefix":"10.1093","volume":"41","author":[{"ORCID":"https:\/\/orcid.org\/0009-0004-5697-5592","authenticated-orcid":false,"given":"Hester","family":"Chapman","sequence":"first","affiliation":[{"name":"Institute of Quantitative and Theoretical Biology, Heinrich Heine University , D\u00fcsseldorf 40255,","place":["Germany"]}]},{"ORCID":"https:\/\/orcid.org\/0000-0001-7356-4075","authenticated-orcid":false,"given":"Miroslav","family":"Kratochv\u00edl","sequence":"additional","affiliation":[{"name":"Luxembourg Centre of Systems Biomedicine, University of Luxembourg , Esch-sur-Alzette L-4367,","place":["Luxembourg"]}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-7229-7398","authenticated-orcid":false,"given":"Oliver","family":"Ebenh\u00f6h","sequence":"additional","affiliation":[{"name":"Institute of Quantitative and Theoretical Biology, Heinrich Heine University , D\u00fcsseldorf 40255,","place":["Germany"]}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-4113-2590","authenticated-orcid":false,"given":"St Elmo","family":"Wilken","sequence":"additional","affiliation":[{"name":"Institute of Quantitative and Theoretical Biology, Heinrich Heine University , D\u00fcsseldorf 40255,","place":["Germany"]}]}],"member":"286","published-online":{"date-parts":[[2025,5,6]]},"reference":[{"key":"2025070408320713000_btaf287-B1","doi-asserted-by":"crossref","first-page":"27162","DOI":"10.1038\/s41598-024-77382-x","article-title":"Uncovering novel polyhydroxyalkanoate biosynthesis genes and unique pathway in yeast Hanseniaspora valbyensis for sustainable bioplastic production","volume":"14","author":"Abd-El-Haleem","year":"2024","journal-title":"Sci Rep"},{"year":"2019","author":"Agrawal","key":"2025070408320713000_btaf287-B2"},{"key":"2025070408320713000_btaf287-B3","doi-asserted-by":"crossref","first-page":"99","DOI":"10.1038\/nature15765","article-title":"Overflow metabolism in Escherichia coli results from efficient proteome allocation","volume":"528","author":"Basan","year":"2015","journal-title":"Nature"},{"key":"2025070408320713000_btaf287-B4","doi-asserted-by":"crossref","first-page":"12663","DOI":"10.1073\/pnas.0609845104","article-title":"Intracellular crowding defines the mode and sequence of substrate uptake by Escherichia coli and constrains its metabolic activity","volume":"104","author":"Beg","year":"2007","journal-title":"Proc Natl Acad Sci USA"},{"key":"2025070408320713000_btaf287-B5","doi-asserted-by":"crossref","first-page":"764","DOI":"10.3389\/fmicb.2020.00764","article-title":"Lachancea fermentati strains isolated from kombucha: fundamental insights, and practical application in low alcohol beer brewing","volume":"11","author":"Bellut","year":"2020","journal-title":"Front Microbiol"},{"key":"2025070408320713000_btaf287-B7","doi-asserted-by":"crossref","first-page":"821","DOI":"10.1093\/femsre\/fuaa034","article-title":"Searching for principles of microbial physiology","volume":"44","author":"Bruggeman","year":"2020","journal-title":"FEMS Microbiol Rev"},{"key":"2025070408320713000_btaf287-B8","doi-asserted-by":"crossref","first-page":"272","DOI":"10.1038\/nbt.4072","article-title":"Recon3d enables a three-dimensional view of gene variation in human metabolism","volume":"36","author":"Brunk","year":"2018","journal-title":"Nat Biotechnol"},{"key":"2025070408320713000_btaf287-B9","doi-asserted-by":"crossref","first-page":"786","DOI":"10.1248\/bpb.32.786","article-title":"Glycan engineering and production of \u2018humanized\u2019glycoprotein in yeast cells","volume":"32","author":"Chiba","year":"2009","journal-title":"Biol Pharm Bull"},{"key":"2025070408320713000_btaf287-B10","doi-asserted-by":"crossref","first-page":"421","DOI":"10.1038\/sj.jim.7000014","article-title":"Effects of mutations in acetate metabolism on high-cell-density growth of Escherichia coli","volume":"24","author":"Contiero","year":"2000","journal-title":"J Ind Microbiol Biotechnol"},{"key":"2025070408320713000_btaf287-B11","doi-asserted-by":"crossref","first-page":"81","DOI":"10.1016\/j.copbio.2017.02.007","article-title":"Lessons on enzyme kinetics from quantitative proteomics","volume":"46","author":"Davidi","year":"2017","journal-title":"Curr Opin Biotechnol"},{"key":"2025070408320713000_btaf287-B12","doi-asserted-by":"crossref","first-page":"e1006858","DOI":"10.1371\/journal.pcbi.1006858","article-title":"The number of active metabolic pathways is bounded by the number of cellular constraints at maximal metabolic rates","volume":"15","author":"De Groot","year":"2019","journal-title":"PLoS Comput Biol"},{"key":"2025070408320713000_btaf287-B13","doi-asserted-by":"crossref","first-page":"1945","DOI":"10.1007\/s10529-006-9182-8","article-title":"Metabolic characterisation of E. coli citrate synthase and phosphoenolpyruvate carboxylase mutants in aerobic cultures","volume":"28","author":"De Maeseneire","year":"2006","journal-title":"Biotechnol Lett"},{"key":"2025070408320713000_btaf287-B14","doi-asserted-by":"crossref","first-page":"689","DOI":"10.1007\/s10295-007-0244-2","article-title":"Minimizing acetate formation in E. coli fermentations","volume":"34","author":"De Mey","year":"2007","journal-title":"J Ind Microbiol Biotechnol"},{"key":"2025070408320713000_btaf287-B15","doi-asserted-by":"crossref","first-page":"17410","DOI":"10.1074\/jbc.274.25.17410","article-title":"Systems properties of the Haemophilus influenzae Rd metabolic genotype","volume":"274","author":"Edwards","year":"1999","journal-title":"J Biol Chem"},{"key":"2025070408320713000_btaf287-B16","doi-asserted-by":"crossref","first-page":"8024","DOI":"10.1074\/jbc.M510016200","article-title":"Latent pathway activation and increased pathway capacity enable Escherichia coli adaptation to loss of key metabolic enzymes","volume":"281","author":"Fong","year":"2006","journal-title":"J Biol Chem"},{"key":"2025070408320713000_btaf287-B17","doi-asserted-by":"crossref","first-page":"1","DOI":"10.1186\/1471-2105-5-175","article-title":"Computation of elementary modes: a unifying framework and the new binary approach","volume":"5","author":"Gagneur","year":"2004","journal-title":"BMC Bioinformatics"},{"key":"2025070408320713000_btaf287-B18","doi-asserted-by":"crossref","first-page":"23182","DOI":"10.1073\/pnas.2001562117","article-title":"Kinetic profiling of metabolic specialists demonstrates stability and consistency of in vivo enzyme turnover numbers","volume":"117","author":"Heckmann","year":"2020","journal-title":"Proc Natl Acad Sci U S A"},{"key":"2025070408320713000_btaf287-B19","first-page":"65","article-title":"Rate control of biological processes","volume":"27","author":"Kacser","year":"1973","journal-title":"Symp Soc Exp Biol"},{"key":"2025070408320713000_btaf287-B20","doi-asserted-by":"crossref","first-page":"339","DOI":"10.1039\/B916446D","article-title":"Genome-scale metabolic network analysis and drug targeting of multi-drug resistant pathogen Acinetobacter baumannii AYE","volume":"6","author":"Kim","year":"2010","journal-title":"Mol Biosyst"},{"key":"2025070408320713000_btaf287-B21","doi-asserted-by":"crossref","first-page":"e1005409","DOI":"10.1371\/journal.pcbi.1005409","article-title":"From elementary flux modes to elementary flux vectors: metabolic pathway analysis with arbitrary linear flux constraints","volume":"13","author":"Klamt","year":"2017","journal-title":"PLoS Comput Biol"},{"key":"2025070408320713000_btaf287-B22","doi-asserted-by":"publisher","first-page":"btaf056","DOI":"10.1093\/bioinformatics\/btaf056","article-title":"Cobrexa 2: tidy and scalable construction of complex metabolic models","volume":"41","author":"Kratochv\u00edl","year":"2025","journal-title":"Bioinformatics"},{"key":"2025070408320713000_btaf287-B23","doi-asserted-by":"crossref","first-page":"4139","DOI":"10.1038\/s41467-023-39840-4","article-title":"Turnover number predictions for kinetically uncharacterized enzymes using machine and deep learning","volume":"14","author":"Kroll","year":"2023","journal-title":"Nat Commun"},{"key":"2025070408320713000_btaf287-B24","doi-asserted-by":"crossref","first-page":"662","DOI":"10.1038\/s41929-022-00798-z","article-title":"Deep learning-based k cat prediction enables improved enzyme-constrained model reconstruction","volume":"5","author":"Li","year":"2022","journal-title":"Nat Catal"},{"year":"2018","author":"Liebermeister","key":"2025070408320713000_btaf287-B25","doi-asserted-by":"publisher","DOI":"10.1101\/483891"},{"key":"2025070408320713000_btaf287-B26","doi-asserted-by":"crossref","first-page":"111568","DOI":"10.1016\/j.mcat.2021.111568","article-title":"Engineering of Saccharomyces pastorianus old yellow enzyme 1 for the synthesis of pharmacologically active (s)-profen derivatives","volume":"507","author":"Liu","year":"2021","journal-title":"Mol Catal"},{"key":"2025070408320713000_btaf287-B27","doi-asserted-by":"crossref","first-page":"102","DOI":"10.1016\/j.ymben.2016.05.006","article-title":"Characterization of physiological responses to 22 gene knockouts in Escherichia coli central carbon metabolism","volume":"37","author":"Long","year":"2016","journal-title":"Metab Eng"},{"key":"2025070408320713000_btaf287-B28","doi-asserted-by":"crossref","first-page":"19","DOI":"10.1186\/s12934-019-1202-1","article-title":"Engineering protein production by rationally choosing a carbon and nitrogen source using e. coli bl21 acetate metabolism knockout strains","volume":"18","author":"Lozano Terol","year":"2019","journal-title":"Microb Cell Fact"},{"key":"2025070408320713000_btaf287-B29","doi-asserted-by":"publisher","first-page":"D479","DOI":"10.1093\/nar\/gku975","article-title":"The complex portal\u2013an encyclopaedia of macromolecular complexes","volume":"43","author":"Meldal","year":"2015","journal-title":"Nucleic Acids Res"},{"key":"2025070408320713000_btaf287-B30","doi-asserted-by":"crossref","first-page":"904","DOI":"10.1038\/nbt.3956","article-title":"iML1515, a knowledgebase that computes Escherichia coli traits","volume":"35","author":"Monk","year":"2017","journal-title":"Nat Biotechnol"},{"key":"2025070408320713000_btaf287-B31","doi-asserted-by":"crossref","first-page":"e1004913","DOI":"10.1371\/journal.pcbi.1004913","article-title":"Constrained allocation flux balance analysis","volume":"12","author":"Mori","year":"2016","journal-title":"PLoS Comput Biol"},{"key":"2025070408320713000_btaf287-B32","doi-asserted-by":"crossref","first-page":"597","DOI":"10.1007\/s002530051101","article-title":"Influence of acetic acid on the growth of escherichia coli k12 during high-cell-density cultivation in a dialysis reactor","volume":"48","author":"Nakano","year":"1997","journal-title":"Appl Microbiol Biotechnol"},{"key":"2025070408320713000_btaf287-B33","doi-asserted-by":"crossref","first-page":"245","DOI":"10.1038\/nbt.1614","article-title":"What is flux balance analysis?","volume":"28","author":"Orth","year":"2010","journal-title":"Nat Biotechnol"},{"key":"2025070408320713000_btaf287-B34","doi-asserted-by":"crossref","first-page":"89","DOI":"10.1007\/BF01032744","article-title":"Physiological constraints in increasing biomass concentration of Escherichia coli b in fed-batch culture","volume":"9","author":"Pan","year":"1987","journal-title":"Biotechnol Lett"},{"key":"2025070408320713000_btaf287-B35","doi-asserted-by":"crossref","first-page":"979","DOI":"10.1016\/j.biotechadv.2009.05.019","article-title":"Constraints-based genome-scale metabolic simulation for systems metabolic engineering","volume":"27","author":"Park","year":"2009","journal-title":"Biotechnol Adv"},{"key":"2025070408320713000_btaf287-B36","doi-asserted-by":"crossref","first-page":"301","DOI":"10.1007\/s00253-021-11727-8","article-title":"Production of galactosylated complex-type n-glycans in glycoengineered Saccharomyces cerevisiae","volume":"106","author":"Piirainen","year":"2022","journal-title":"Appl Microbiol Biotechnol"},{"key":"2025070408320713000_btaf287-B37","doi-asserted-by":"crossref","first-page":"935","DOI":"10.15252\/msb.20167411","article-title":"Improving the phenotype predictions of a yeast genome-scale metabolic model by incorporating enzymatic constraints","volume":"13","author":"S\u00e1nchez","year":"2017","journal-title":"Mol Syst Biol"},{"key":"2025070408320713000_btaf287-B38","doi-asserted-by":"crossref","first-page":"104","DOI":"10.1038\/nbt.3418","article-title":"The quantitative and condition-dependent escherichia coli proteome","volume":"34","author":"Schmidt","year":"2016","journal-title":"Nat Biotechnol"},{"key":"2025070408320713000_btaf287-B39","first-page":"79","article-title":"Refined algorithm and computer program for calculating all non\u2013negative fluxes admissible in steady states of biochemical reaction systems with or without some flux rates fixed","volume":"9","author":"Schuster","year":"1993","journal-title":"Comput Appl Biosci"},{"year":"2009","author":"Terzer","key":"2025070408320713000_btaf287-B40"},{"key":"2025070408320713000_btaf287-B6","doi-asserted-by":"publisher","first-page":"D609","DOI":"10.1093\/nar\/gkae1010","volume-title":"Nucleic Acids Res","year":"2025"},{"key":"2025070408320713000_btaf287-B41","doi-asserted-by":"crossref","first-page":"274","DOI":"10.1049\/iet-syb:20060073","article-title":"Enumerating constrained elementary flux vectors of metabolic networks","volume":"1","author":"Urbanczik","year":"2007","journal-title":"IET Syst Biol"},{"key":"2025070408320713000_btaf287-B42","doi-asserted-by":"crossref","first-page":"btae691","DOI":"10.1093\/bioinformatics\/btae691","article-title":"Sensitivities in protein allocation models reveal distribution of metabolic capacity and flux control","volume":"40","author":"van den Bogaard","year":"2024","journal-title":"Bioinformatics"},{"key":"2025070408320713000_btaf287-B43","doi-asserted-by":"crossref","first-page":"994","DOI":"10.1038\/nbt1094-994","article-title":"Metabolic flux balancing: basic concepts, scientific and practical use","volume":"12","author":"Varma","year":"1994","journal-title":"Nat Biotechnol"},{"key":"2025070408320713000_btaf287-B44","doi-asserted-by":"crossref","first-page":"72","DOI":"10.1016\/j.ymben.2022.09.002","article-title":"Interrogating the effect of enzyme kinetics on metabolism using differentiable constraint-based models","volume":"74","author":"Wilken","year":"2022","journal-title":"Metab Eng"},{"key":"2025070408320713000_btaf287-B45","doi-asserted-by":"crossref","first-page":"164","DOI":"10.1016\/j.ymben.2004.02.004","article-title":"Effect of zwf gene knockout on the metabolism of Escherichia coli grown on glucose or acetate","volume":"6","author":"Zhao","year":"2004","journal-title":"Metab Eng"}],"container-title":["Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/academic.oup.com\/bioinformatics\/advance-article-pdf\/doi\/10.1093\/bioinformatics\/btaf287\/63069969\/btaf287.pdf","content-type":"application\/pdf","content-version":"am","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/41\/6\/btaf287\/63069969\/btaf287.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/41\/6\/btaf287\/63069969\/btaf287.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2025,7,4]],"date-time":"2025-07-04T12:32:19Z","timestamp":1751632339000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article\/doi\/10.1093\/bioinformatics\/btaf287\/8125804"}},"subtitle":[],"editor":[{"given":"Jianlin","family":"Cheng","sequence":"additional","affiliation":[]}],"short-title":[],"issued":{"date-parts":[[2025,5,6]]},"references-count":45,"journal-issue":{"issue":"6","published-print":{"date-parts":[[2025,6,2]]}},"URL":"https:\/\/doi.org\/10.1093\/bioinformatics\/btaf287","relation":{},"ISSN":["1367-4811"],"issn-type":[{"type":"electronic","value":"1367-4811"}],"subject":[],"published-other":{"date-parts":[[2025,6]]},"published":{"date-parts":[[2025,5,6]]},"article-number":"btaf287"}}