{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,1,7]],"date-time":"2026-01-07T23:00:24Z","timestamp":1767826824494,"version":"3.49.0"},"reference-count":61,"publisher":"Oxford University Press (OUP)","issue":"14","license":[{"start":{"date-parts":[[2018,3,7]],"date-time":"2018-03-07T00:00:00Z","timestamp":1520380800000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/academic.oup.com\/journals\/pages\/about_us\/legal\/notices"}],"funder":[{"name":"European Union\u2019s Horizon 2020","award":["675585"],"award-info":[{"award-number":["675585"]}]},{"name":"MSCA ITN \u2018SyMBioSys\u2019"},{"name":"Spanish MINECO\/FEDER"},{"name":"SYNBIOFACTORY","award":["DPI2014-55276-C5-2-R"],"award-info":[{"award-number":["DPI2014-55276-C5-2-R"]}]},{"name":"SYNBIOCONTROL","award":["DPI2017-82896-C2-2-R"],"award-info":[{"award-number":["DPI2017-82896-C2-2-R"]}]},{"name":"IMPROWINE","award":["AGL2015-67504-C3-2-R"],"award-info":[{"award-number":["AGL2015-67504-C3-2-R"]}]},{"DOI":"10.13039\/100003973","name":"MSCA","doi-asserted-by":"publisher","id":[{"id":"10.13039\/100003973","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/501100003961","name":"ESR","doi-asserted-by":"publisher","id":[{"id":"10.13039\/501100003961","id-type":"DOI","asserted-by":"publisher"}]}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2018,7,15]]},"abstract":"<jats:title>Abstract<\/jats:title>\n               <jats:sec>\n                  <jats:title>Motivation<\/jats:title>\n                  <jats:p>Optimality principles have been used to explain many biological processes and systems. However, the functions being optimized are in general unknown a priori. Here we present an inverse optimal control framework for modeling dynamics in systems biology. The objective is to identify the underlying optimality principle from observed time-series data and simultaneously estimate unmeasured time-dependent inputs and time-invariant model parameters. As a special case, we also consider the problem of optimal simultaneous estimation of inputs and parameters from noisy data. After presenting a general statement of the inverse optimal control problem, and discussing special cases of interest, we outline numerical strategies which are scalable and robust.<\/jats:p>\n               <\/jats:sec>\n               <jats:sec>\n                  <jats:title>Results<\/jats:title>\n                  <jats:p>We discuss the existence, relevance and implications of identifiability issues in the above problems. We present a robust computational approach based on regularized cost functions and the use of suitable direct numerical methods based on the control-vector parameterization approach. To avoid convergence to local solutions, we make use of hybrid global-local methods. We illustrate the performance and capabilities of this approach with several challenging case studies, including simulated and real data. We pay particular attention to the computational scalability of our approach (with the objective of considering large numbers of inputs and states). We provide a software implementation of both the methods and the case studies.<\/jats:p>\n               <\/jats:sec>\n               <jats:sec>\n                  <jats:title>Availability and implementation<\/jats:title>\n                  <jats:p>The code used to obtain the results reported here is available at https:\/\/zenodo.org\/record\/1009541.<\/jats:p>\n               <\/jats:sec>\n               <jats:sec>\n                  <jats:title>Supplementary information<\/jats:title>\n                  <jats:p>Supplementary data are available at Bioinformatics online.<\/jats:p>\n               <\/jats:sec>","DOI":"10.1093\/bioinformatics\/bty139","type":"journal-article","created":{"date-parts":[[2018,3,6]],"date-time":"2018-03-06T20:11:00Z","timestamp":1520367060000},"page":"2433-2440","source":"Crossref","is-referenced-by-count":25,"title":["Optimality and identification of dynamic models in systems biology: an inverse optimal control framework"],"prefix":"10.1093","volume":"34","author":[{"ORCID":"https:\/\/orcid.org\/0000-0003-4199-6126","authenticated-orcid":false,"given":"Nikolaos","family":"Tsiantis","sequence":"first","affiliation":[{"name":"Bioprocess Engineering Group, Spanish National Research Council, IIM-CSIC Vigo, Spain"},{"name":"Department of Chemical Engineering, University of Vigo Vigo, Spain"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-1978-2626","authenticated-orcid":false,"given":"Eva","family":"Balsa-Canto","sequence":"additional","affiliation":[{"name":"Bioprocess Engineering Group, Spanish National Research Council, IIM-CSIC Vigo, Spain"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-4245-0320","authenticated-orcid":false,"given":"Julio R","family":"Banga","sequence":"additional","affiliation":[{"name":"Bioprocess Engineering Group, Spanish National Research Council, IIM-CSIC Vigo, Spain"}]}],"member":"286","published-online":{"date-parts":[[2018,3,7]]},"reference":[{"key":"2023012713012032700_bty139-B1","volume-title":"Optima for Animals","author":"Alexander","year":"1996"},{"key":"2023012713012032700_bty139-B2","doi-asserted-by":"crossref","first-page":"3357","DOI":"10.1093\/bioinformatics\/btw411","article-title":"Amigo2, a toolbox for dynamic modeling, optimization and control in systems biology","volume":"32","author":"Balsa-Canto","year":"2016","journal-title":"Bioinformatics"},{"key":"2023012713012032700_bty139-B3","doi-asserted-by":"crossref","first-page":"47.","DOI":"10.1186\/1752-0509-2-47","article-title":"Optimization in computational systems biology","volume":"2","author":"Banga","year":"2008","journal-title":"BMC Syst. Biol"},{"key":"2023012713012032700_bty139-B4","doi-asserted-by":"crossref","first-page":"407","DOI":"10.1016\/j.jbiotec.2005.02.013","article-title":"Dynamic optimization of bioprocesses: efficient and robust numerical strategies","volume":"117","author":"Banga","year":"2005","journal-title":"J. Biotechnol"},{"key":"2023012713012032700_bty139-B5","doi-asserted-by":"crossref","first-page":"67","DOI":"10.1016\/j.biosystems.2010.04.007","article-title":"Modelling the optimal timing in metabolic pathway activation\u2014use of pontryagin\u2019s maximum principle and role of the golden section","volume":"101","author":"Bartl","year":"2010","journal-title":"Biosystems"},{"key":"2023012713012032700_bty139-B6","doi-asserted-by":"crossref","DOI":"10.1038\/ncomms3243","article-title":"Dynamic optimization identifies optimal programmes for pathway regulation in prokaryotes","volume":"4","author":"Bartl","year":"2013","journal-title":"Nat. Commun"},{"key":"2023012713012032700_bty139-B7","article-title":"Perspectives on theory at the interface of physics and biology","volume":"81","author":"Bialek","year":"2017","journal-title":"Reports on Progress in Physics"},{"key":"2023012713012032700_bty139-B8","doi-asserted-by":"crossref","first-page":"2610","DOI":"10.1093\/bioinformatics\/btr431","article-title":"GenSSI: a software toolbox for structural identifiability analysis of biological models","volume":"27","author":"Chis","year":"2011","journal-title":"Bioinformatics"},{"key":"2023012713012032700_bty139-B9","doi-asserted-by":"crossref","first-page":"e27755","DOI":"10.1371\/journal.pone.0027755","article-title":"Structural identifiability of systems biology models: a critical comparison of methods","volume":"6","author":"Chis","year":"2011","journal-title":"PloS One"},{"key":"2023012713012032700_bty139-B10","first-page":"1192","article-title":"Inverse optimal control based identification of optimality criteria in whole-body human walking on level ground","author":"Clever","year":"2016","journal-title":"IEEE EMBS Conf. Biomed. Robot. Biomech"},{"key":"2023012713012032700_bty139-B11","doi-asserted-by":"crossref","first-page":"1.","DOI":"10.1186\/1752-0509-8-1","article-title":"Global dynamic optimization approach to predict activation in metabolic pathways","volume":"8","author":"de Hijas-Liste","year":"2014","journal-title":"BMC Syst. Biol"},{"key":"2023012713012032700_bty139-B12","doi-asserted-by":"crossref","first-page":"163.","DOI":"10.1186\/s12859-015-0587-z","article-title":"Optimal programs of pathway control: dissecting the influence of pathway topology and feedback inhibition on pathway regulation","volume":"16","author":"de Hijas-Liste","year":"2015","journal-title":"BMC Bioinformatics"},{"key":"2023012713012032700_bty139-B13","doi-asserted-by":"crossref","first-page":"797","DOI":"10.1016\/j.copbio.2013.04.008","article-title":"Optimality in evolution: new insights from synthetic biology","volume":"24","author":"de Vos","year":"2013","journal-title":"Curr. Opin. Biotechnol"},{"key":"2023012713012032700_bty139-B14","doi-asserted-by":"crossref","first-page":"588.","DOI":"10.1038\/nature03842","article-title":"Optimality and evolutionary tuning of the expression level of a protein","volume":"436","author":"Dekel","year":"2005","journal-title":"Nature"},{"key":"2023012713012032700_bty139-B15","doi-asserted-by":"crossref","first-page":"603","DOI":"10.1098\/rsif.2006.0143","article-title":"Systems interface biology","volume":"3","author":"Doyle","year":"2006","journal-title":"J. R. Soc. Interface"},{"key":"2023012713012032700_bty139-B16","doi-asserted-by":"crossref","first-page":"4388","DOI":"10.1021\/ie801717t","article-title":"Dynamic optimization of nonlinear processes with an enhanced scatter search method","volume":"48","author":"Egea","year":"2009","journal-title":"Ind. Eng. Chem. Res"},{"key":"2023012713012032700_bty139-B17","doi-asserted-by":"crossref","DOI":"10.1038\/srep20772","article-title":"Learning (from) the errors of a systems biology model","volume":"6","author":"Engelhardt","year":"2016","journal-title":"Sci. Rep"},{"key":"2023012713012032700_bty139-B18","doi-asserted-by":"crossref","first-page":"20170332.","DOI":"10.1098\/rsif.2017.0332","article-title":"A bayesian approach to estimating hidden variables as well as missing and wrong molecular interactions in ordinary differential equation-based mathematical models","volume":"14","author":"Engelhardt","year":"2017","journal-title":"J. R. Soc. Interface"},{"key":"2023012713012032700_bty139-B19","doi-asserted-by":"crossref","first-page":"1035","DOI":"10.1042\/BST20170137","article-title":"Deciphering the regulation of metabolism with dynamic optimization: an overview of recent advances","volume":"45","author":"Ewald","year":"2017","journal-title":"Biochem. Soc. Trans"},{"key":"2023012713012032700_bty139-B20","doi-asserted-by":"crossref","first-page":"127","DOI":"10.1038\/nrn2787","article-title":"The free-energy principle: a unified brain theory?","volume":"11","author":"Friston","year":"2010","journal-title":"Nat. Rev. Neurosci"},{"key":"2023012713012032700_bty139-B21","doi-asserted-by":"crossref","first-page":"74.","DOI":"10.1186\/s12918-015-0219-2","article-title":"Robust and efficient parameter estimation in dynamic models of biological systems","volume":"9","author":"G\u00e1bor","year":"2015","journal-title":"BMC Syst. Biol"},{"key":"2023012713012032700_bty139-B22","doi-asserted-by":"crossref","first-page":"e1004802.","DOI":"10.1371\/journal.pcbi.1004802","article-title":"Dynamical allocation of cellular resources as an optimal control problem: novel insights into microbial growth strategies","volume":"12","author":"Giordano","year":"2016","journal-title":"PLoS Comput. Biol"},{"key":"2023012713012032700_bty139-B23","doi-asserted-by":"crossref","first-page":"279.","DOI":"10.1109\/TCBB.2007.070203","article-title":"Multiobjective optimization in bioinformatics and computational biology","volume":"4","author":"Handl","year":"2007","journal-title":"IEEE\/ACM Trans. Comput. Biol. Bioinf"},{"key":"2023012713012032700_bty139-B24","author":"Hatz","year":"2014"},{"key":"2023012713012032700_bty139-B25","doi-asserted-by":"crossref","first-page":"A1707","DOI":"10.1137\/110823390","article-title":"Estimating parameters in optimal control problems","volume":"34","author":"Hatz","year":"2012","journal-title":"SIAM J. Sci. Comput"},{"key":"2023012713012032700_bty139-B26","doi-asserted-by":"crossref","first-page":"1277","DOI":"10.1002\/aic.690420509","article-title":"Analysis and design of metabolic reaction networks via mixed-integer linear optimization","volume":"42","author":"Hatzimanikatis","year":"1996","journal-title":"AIChE J"},{"key":"2023012713012032700_bty139-B27","doi-asserted-by":"crossref","first-page":"20151046.","DOI":"10.1098\/rsif.2015.1046","article-title":"Synthetic biology and regulatory networks: where metabolic systems biology meets control engineering","volume":"13","author":"He","year":"2016","journal-title":"J. R. Soc. Interface"},{"key":"2023012713012032700_bty139-B28","doi-asserted-by":"crossref","DOI":"10.1007\/978-1-4613-1161-4","volume-title":"The Regulation of Cellular Systems","author":"Heinrich","year":"1996"},{"key":"2023012713012032700_bty139-B29","first-page":"1","article-title":"Mathematical analysis of enzymic reaction systems using optimization principles","volume":"201","author":"Heinrich","year":"1991","journal-title":"febs J"},{"key":"2023012713012032700_bty139-B30","doi-asserted-by":"crossref","first-page":"191","DOI":"10.1111\/j.1432-1033.1997.0191a.x","article-title":"Theoretical approaches to the evolutionary optimization of glycolysis","volume":"243","author":"Heinrich","year":"1997","journal-title":"Eur. J. Biochem"},{"key":"2023012713012032700_bty139-B31","doi-asserted-by":"crossref","first-page":"e1002548.","DOI":"10.1371\/journal.pcbi.1002548","article-title":"Potassium starvation in yeast: mechanisms of homeostasis revealed by mathematical modeling","volume":"8","author":"Kahm","year":"2012","journal-title":"PLoS Comput. Biol"},{"key":"2023012713012032700_bty139-B32","doi-asserted-by":"crossref","first-page":"13","DOI":"10.1109\/LLS.2016.2615081","article-title":"Testing the pattern of akt activation by variational parameter estimation","volume":"2","author":"Kaschek","year":"2016","journal-title":"IEEE Life Sci. Lett"},{"key":"2023012713012032700_bty139-B33","doi-asserted-by":"crossref","first-page":"5406","DOI":"10.1046\/j.1432-1033.2002.03223.x","article-title":"Prediction of temporal gene expression","volume":"269","author":"Klipp","year":"2002","journal-title":"Eur. J. Biochem"},{"key":"2023012713012032700_bty139-B34","doi-asserted-by":"crossref","first-page":"B988","DOI":"10.1137\/15M103306X","article-title":"Modular parameter identification of biomolecular networks","volume":"38","author":"Lang","year":"2016","journal-title":"SIAM J. Sci. Comput"},{"key":"2023012713012032700_bty139-B35","doi-asserted-by":"crossref","first-page":"208301\/1","DOI":"10.1103\/PhysRevLett.91.208301","article-title":"Manipulation of self-aggregation patterns and waves in a reaction-diffusion system by optimal boundary control strategies","volume":"91","author":"Lebiedz","year":"2003","journal-title":"Phys. Rev. Lett"},{"key":"2023012713012032700_bty139-B36","doi-asserted-by":"crossref","first-page":"15","DOI":"10.1038\/269015a0","article-title":"Decision making in animals","volume":"269","author":"McFarland","year":"1977","journal-title":"Nature"},{"key":"2023012713012032700_bty139-B37","doi-asserted-by":"crossref","first-page":"869","DOI":"10.1093\/bioinformatics\/14.10.869","article-title":"Non-linear optimization of biochemical pathways: applications to metabolic engineering and parameter estimation","volume":"14","author":"Mendes","year":"1998","journal-title":"Bioinformatics"},{"key":"2023012713012032700_bty139-B38","first-page":"103","volume-title":"Optimal Control for Applications in Medical and Rehabilitation Technology: Challenges and Solutions","author":"Mombaur","year":"2016"},{"key":"2023012713012032700_bty139-B39","doi-asserted-by":"crossref","first-page":"86.","DOI":"10.1186\/s12918-016-0328-6","article-title":"Dynamic optimization of biological networks under parametric uncertainty","volume":"10","author":"Nimmegeers","year":"2016","journal-title":"BMC Syst. Biol"},{"key":"2023012713012032700_bty139-B40","doi-asserted-by":"crossref","first-page":"110","DOI":"10.1049\/iet-syb.2010.0044","article-title":"Optimal control of metabolic networks with saturable enzyme kinetics","volume":"5","author":"Oyarz\u00fan","year":"2011","journal-title":"IET Syst. Biol"},{"key":"2023012713012032700_bty139-B41","doi-asserted-by":"crossref","first-page":"1851","DOI":"10.1007\/s11538-009-9427-5","article-title":"Sequential activation of metabolic pathways: a dynamic optimization approach","volume":"71","author":"Oyarz\u00fan","year":"2009","journal-title":"Bull. Math. Biol"},{"key":"2023012713012032700_bty139-B42","doi-asserted-by":"crossref","first-page":"4910","DOI":"10.1109\/ACC.2015.7172103","volume-title":"2015 American Control Conference (ACC)","author":"Panchea","year":"2015"},{"key":"2023012713012032700_bty139-B43","doi-asserted-by":"crossref","first-page":"27","DOI":"10.1038\/348027a0","article-title":"Optimality theory in evolutionary biology","volume":"348","author":"Parker","year":"1990","journal-title":"Nature"},{"key":"2023012713012032700_bty139-B44","doi-asserted-by":"crossref","first-page":"299","DOI":"10.1023\/A:1001841615640","article-title":"Bionics, biological systems and the principle of optimal design","volume":"46","author":"Popescu","year":"1998","journal-title":"Acta Biotheor"},{"key":"2023012713012032700_bty139-B45","volume-title":"Mathematical Principles in Biology and Their Applications","author":"Rashevsky","year":"1961"},{"key":"2023012713012032700_bty139-B46","doi-asserted-by":"crossref","DOI":"10.1007\/978-1-4899-6419-9","volume-title":"Optimality Principles in Biology","author":"Rosen","year":"1967"},{"key":"2023012713012032700_bty139-B47","doi-asserted-by":"crossref","first-page":"418","DOI":"10.1086\/284070","article-title":"The application of optimal control theory to the general life history problem","volume":"121","author":"Schaffer","year":"1983","journal-title":"Am. Nat"},{"key":"2023012713012032700_bty139-B48","doi-asserted-by":"crossref","first-page":"i529","DOI":"10.1093\/bioinformatics\/bts393","article-title":"Comprehensive estimation of input signals and dynamics in biochemical reaction networks","volume":"28","author":"Schelker","year":"2012","journal-title":"Bioinformatics"},{"key":"2023012713012032700_bty139-B49","doi-asserted-by":"crossref","first-page":"795","DOI":"10.1016\/j.jprocont.2003.12.005","article-title":"Optimal adaptive control of (bio) chemical reactors: past, present and future","volume":"14","author":"Smets","year":"2004","journal-title":"J. Process Control"},{"key":"2023012713012032700_bty139-B50","doi-asserted-by":"crossref","first-page":"31","DOI":"10.1146\/annurev.es.09.110178.000335","article-title":"Optimization theory in evolution","volume":"9","author":"Smith","year":"1978","journal-title":"Annu. Rev. Ecol. Syst"},{"key":"2023012713012032700_bty139-B51","doi-asserted-by":"crossref","first-page":"254","DOI":"10.1016\/j.automatica.2015.10.037","article-title":"Shaping pulses to control bistable systems: analysis, computation and counterexamples","volume":"63","author":"Sootla","year":"2016","journal-title":"Automatica"},{"key":"2023012713012032700_bty139-B52","doi-asserted-by":"crossref","first-page":"569.","DOI":"10.1038\/435569a","article-title":"The best solution","volume":"435","author":"Sutherland","year":"2005","journal-title":"Nature"},{"key":"2023012713012032700_bty139-B53","doi-asserted-by":"crossref","first-page":"907.","DOI":"10.1038\/nn1309","article-title":"Optimality principles in sensorimotor control","volume":"7","author":"Todorov","year":"2004","journal-title":"Nat. Neurosci"},{"key":"2023012713012032700_bty139-B54","doi-asserted-by":"crossref","first-page":"207","DOI":"10.1007\/s10928-016-9467-z","article-title":"Input estimation for drug discovery using optimal control and markov chain monte carlo approaches","volume":"43","author":"Tr\u00e4g\u00e5rdh","year":"2016","journal-title":"J. Pharmacokinet. Pharm"},{"key":"2023012713012032700_bty139-B55","doi-asserted-by":"crossref","DOI":"10.3389\/fbioe.2017.00024","article-title":"Input estimation for extended-release formulations exemplified with exenatide","volume":"5","author":"Tr\u00e4g\u00e5rdh","year":"2017","journal-title":"Front. Bioeng. Biotechnol"},{"key":"2023012713012032700_bty139-B56","doi-asserted-by":"crossref","first-page":"2123","DOI":"10.1021\/ie00033a015","article-title":"Solution of a class of multistage dynamic optimization problems. 2. problems with path constraints","volume":"33","author":"Vassiliadis","year":"1994","journal-title":"Ind. Eng. Chem. Res"},{"key":"2023012713012032700_bty139-B57","doi-asserted-by":"crossref","first-page":"469","DOI":"10.1016\/j.jtbi.2014.10.035","article-title":"Dynamic optimization of metabolic networks coupled with gene expression","volume":"365","author":"Waldherr","year":"2015","journal-title":"J. Theor. Biol"},{"key":"2023012713012032700_bty139-B58","volume-title":"Identification of Parametric Models from Experimental Data","author":"Walter","year":"1997"},{"key":"2023012713012032700_bty139-B59","doi-asserted-by":"crossref","first-page":"33","DOI":"10.1016\/j.arcontrol.2008.02.001","article-title":"The role of control and system theory in systems biology","volume":"32","author":"Wellstead","year":"2008","journal-title":"Annu. Rev. Control"},{"key":"2023012713012032700_bty139-B60","doi-asserted-by":"crossref","first-page":"515.","DOI":"10.1038\/msb.2011.46","article-title":"Optimal regulatory strategies for metabolic pathways in Escherichia coli depending on protein costs","volume":"7","author":"Wessely","year":"2011","journal-title":"Mol. Syst. Biol"},{"key":"2023012713012032700_bty139-B61","doi-asserted-by":"crossref","first-page":"486","DOI":"10.1038\/ng1348","article-title":"Just-in-time transcription program in metabolic pathways","volume":"36","author":"Zaslaver","year":"2004","journal-title":"Nat. Genet"}],"container-title":["Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/34\/14\/2433\/48917317\/bioinformatics_34_14_2433.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/34\/14\/2433\/48917317\/bioinformatics_34_14_2433.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2023,1,27]],"date-time":"2023-01-27T13:52:55Z","timestamp":1674827575000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article\/34\/14\/2433\/4924215"}},"subtitle":[],"editor":[{"given":"Jonathan","family":"Wren","sequence":"additional","affiliation":[]}],"short-title":[],"issued":{"date-parts":[[2018,3,7]]},"references-count":61,"journal-issue":{"issue":"14","published-print":{"date-parts":[[2018,7,15]]}},"URL":"https:\/\/doi.org\/10.1093\/bioinformatics\/bty139","relation":{},"ISSN":["1367-4803","1367-4811"],"issn-type":[{"value":"1367-4803","type":"print"},{"value":"1367-4811","type":"electronic"}],"subject":[],"published-other":{"date-parts":[[2018,7,15]]},"published":{"date-parts":[[2018,3,7]]}}}