{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,10,8]],"date-time":"2026-10-08T07:02:59Z","timestamp":1791442979576,"version":"4.3.3"},"reference-count":49,"publisher":"Frontiers Media SA","license":[{"start":{"date-parts":[[2022,11,9]],"date-time":"2022-11-09T00:00:00Z","timestamp":1667952000000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"content-domain":{"domain":["frontiersin.org"],"crossmark-restriction":true},"short-container-title":["Front. Physiol."],"abstract":"<jats:p>Computational physiological models continue to increase in complexity, however, the task of efficiently calibrating the model to available clinical data remains a significant challenge. One part of this challenge is associated with long calibration times, which present a barrier for the routine application of model-based prediction in clinical practice. Another aspect of this challenge is the limited available data for the unique calibration of complex models. Therefore, to calibrate a patient-specific model, it may be beneficial to verify that task-specific model predictions have acceptable uncertainty, rather than requiring all parameters to be uniquely identified. We have developed a pipeline that reduces the set of fitting parameters to make them structurally identifiable and to improve the efficiency of a subsequent Markov Chain Monte Carlo (MCMC) analysis. MCMC was used to find the optimal parameter values and to determine the confidence interval of a task-specific prediction. This approach was demonstrated on numerical experiments where a lumped parameter model of the cardiovascular system was calibrated to brachial artery cuff pressure, echocardiogram volume measurements, and synthetic cerebral blood flow data that approximates what can be obtained from 4D-flow MRI data. This pipeline provides a cerebral arterial pressure prediction that may be useful for determining the risk of hemorrhagic stroke. For a set of three patients, this pipeline successfully reduced the parameter set of a cardiovascular system model from 12 parameters to 8\u201310 structurally identifiable parameters. This enabled a significant<jats:inline-formula><mml:math xmlns:mml=\"http:\/\/www.w3.org\/1998\/Math\/MathML\" id=\"m1\"><mml:mrow><mml:mo stretchy=\"false\">(<\/mml:mo><mml:mrow><mml:mo>&gt;<\/mml:mo><mml:mn>4<\/mml:mn><mml:mo>\u00d7<\/mml:mo><\/mml:mrow><mml:mo stretchy=\"false\">)<\/mml:mo><\/mml:mrow><\/mml:math><\/jats:inline-formula>efficiency improvement in determining confidence intervals on predictions of pressure compared to performing a naive MCMC analysis with the full parameter set. This demonstrates the potential that the proposed pipeline has in helping address one of the key challenges preventing clinical application of such models. Additionally, for each patient, the MCMC approach yielded a 95% confidence interval on systolic blood pressure prediction in the middle cerebral artery smaller than \u00b110\u00a0mmHg (\u00b11.3\u00a0kPa). The proposed pipeline exploits available high-performance computing parallelism to allow straightforward automation for general models and arbitrary data sets, enabling automated calibration of a parameter set that is specific to the available clinical data with minimal user interaction.<\/jats:p>","DOI":"10.3389\/fphys.2022.1018134","type":"journal-article","created":{"date-parts":[[2022,11,9]],"date-time":"2022-11-09T08:10:12Z","timestamp":1667981412000},"update-policy":"https:\/\/doi.org\/10.3389\/crossmark-policy","source":"Crossref","is-referenced-by-count":19,"title":["Automated model calibration with parallel MCMC: Applications for a cardiovascular system model"],"prefix":"10.3389","volume":"13","author":[{"given":"Finbar","family":"Argus","sequence":"first","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Debbie","family":"Zhao","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Thiranja P.","family":"Babarenda Gamage","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Martyn P.","family":"Nash","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Gonzalo D.","family":"Maso Talou","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"1965","published-online":{"date-parts":[[2022,11,9]]},"reference":[{"key":"B1","doi-asserted-by":"publisher","first-page":"305","DOI":"10.1093\/ajh\/hpx183","article-title":"Correlation of blood pressure variability as measured by clinic, self-measurement at home, and ambulatory blood pressure monitoring","volume":"31","author":"Abell\u00e1n-Huerta","year":"2018","journal-title":"Am. J. Hypertens."},{"key":"B2","doi-asserted-by":"publisher","first-page":"H899","DOI":"10.1152\/ajpheart.00230.2014","article-title":"An integrated mathematical model of the human cardiopulmonary system: Model development","volume":"310","author":"Albanese","year":"2016","journal-title":"Am. J. Physiol. Heart Circ. Physiol."},{"key":"B3","doi-asserted-by":"publisher","DOI":"10.5281\/zenodo.7222397","article-title":"circulatory_autogen","author":"Argus","year":"2022","journal-title":"zenodo"},{"key":"B4","doi-asserted-by":"publisher","first-page":"652","DOI":"10.1016\/j.medengphy.2012.07.011","article-title":"A dimensionally-heterogeneous closed-loop model for the cardiovascular system and its applications","volume":"35","author":"Blanco","year":"2013","journal-title":"Med. Eng. Phys."},{"key":"B5","doi-asserted-by":"publisher","first-page":"1303","DOI":"10.1007\/s10237-014-0574-8","article-title":"Blood flow distribution in an anatomically detailed arterial network model: Criteria and algorithms","volume":"13","author":"Blanco","year":"2014","journal-title":"Biomech. Model. Mechanobiol."},{"key":"B6","doi-asserted-by":"publisher","first-page":"736","DOI":"10.1109\/TBME.2014.2364522","article-title":"An anatomically detailed arterial network model for one-dimensional computational hemodynamics","volume":"62","author":"Blanco","year":"2015","journal-title":"IEEE Trans. Biomed. Eng."},{"key":"B7","doi-asserted-by":"publisher","first-page":"1015","DOI":"10.1029\/2000WR900350","article-title":"Practical identifiability analysis of large environmental simulation models","volume":"37","author":"Brun","year":"2001","journal-title":"Water Resour. Res."},{"key":"B8","doi-asserted-by":"publisher","first-page":"3513","DOI":"10.1111\/j.1742-4658.2012.08725.x","article-title":"Conclusions via unique predictions obtained despite unidentifiability - new definitions and a general method","volume":"279","author":"Cedersund","year":"2012","journal-title":"FEBS J."},{"key":"B9","doi-asserted-by":"publisher","first-page":"903","DOI":"10.1111\/j.1742-4658.2008.06845.x","article-title":"Systems biology: Model based evaluation and comparison of potential explanations for given biological data","volume":"276","author":"Cedersund","year":"2009","journal-title":"FEBS J."},{"key":"B10","doi-asserted-by":"publisher","first-page":"54","DOI":"10.3182\/20100707-3-BE-2012.0045","article-title":"Global sensitivity and identifiability implications in systems biology","volume":"11","author":"Dobre","year":"2010","journal-title":"IFAC Proc. Vol."},{"key":"B11","doi-asserted-by":"publisher","first-page":"H222","DOI":"10.1152\/ajpheart.00857.2014","article-title":"Reducing the number of parameters in 1D arterial blood flow modeling: Less is more for patient-specific simulations","volume":"309","author":"Epstein","year":"2015","journal-title":"Am. J. Physiol. Heart Circ. Physiol."},{"key":"B12","doi-asserted-by":"publisher","first-page":"644","DOI":"10.1038\/msb.2012.69","article-title":"Properties of cell death models calibrated and compared using Bayesian approaches","volume":"9","author":"Eydgahi","year":"2013","journal-title":"Mol. Syst. Biol."},{"key":"B13","doi-asserted-by":"publisher","first-page":"24","DOI":"10.21105\/joss.00024","article-title":"corner. py: Scatterplot matrices in Python","volume":"1","author":"Foreman-Mackey","year":"2016","journal-title":"J. Open Source Softw."},{"key":"B14","doi-asserted-by":"publisher","first-page":"306","DOI":"10.1086\/670067","article-title":"Emcee : The MCMC Hammer","volume":"125","author":"Foreman-Mackey","year":"2013","journal-title":"Publ. Astronomical Soc. Pac."},{"key":"B15","doi-asserted-by":"publisher","first-page":"20180546","DOI":"10.1098\/rsif.2018.0546","article-title":"Optimization of topological complexity for one-dimensional arterial blood flow models","volume":"15","author":"Fossan","year":"2018","journal-title":"J. R. Soc. Interface"},{"key":"B16","doi-asserted-by":"publisher","first-page":"26","DOI":"10.3389\/fphys.2015.00026","article-title":"OpenCOR: A modular and interoperable approach to computational biology","volume":"6","author":"Garny","year":"2015","journal-title":"Front. Physiol."},{"key":"B17","doi-asserted-by":"crossref","DOI":"10.1093\/oso\/9780198522669.003.0010","article-title":"Evaluating the accuracy of sampling-based approaches to the calculation of posterior moments","volume-title":"Bayesian statistics","author":"Geweke","year":"1992"},{"key":"B18","doi-asserted-by":"publisher","first-page":"473","DOI":"10.1214\/ss\/1177011137","article-title":"Practical Markov chain Monte Carlo","volume":"7","author":"Geyer","year":"1992","journal-title":"Stat. Sci."},{"key":"B19","first-page":"107","article-title":"Strategies for improving MCMC","volume-title":"Markov chain Monte Carlo in practice","author":"Gilks","year":"2020"},{"key":"B20","doi-asserted-by":"publisher","first-page":"e1007669","DOI":"10.1371\/journal.pcbi.1007669","article-title":"Parallel Tempering with Lasso for model reduction in systems biology","volume":"16","author":"Gupta","year":"2020","journal-title":"PLoS Comput. Biol."},{"key":"B21","doi-asserted-by":"publisher","first-page":"293","DOI":"10.1016\/j.mbs.2013.04.002","article-title":"High-dimensional Bayesian parameter estimation: Case study for a model of JAK2\/STAT5 signaling","volume":"246","author":"Hug","year":"2013","journal-title":"Math. Biosci."},{"key":"B22","doi-asserted-by":"publisher","first-page":"2595","DOI":"10.1098\/rsta.2010.0048","article-title":"A vision and strategy for the virtual physiological human in 2010 and beyond","volume":"368","author":"Hunter","year":"2010","journal-title":"Philos. Trans. A Math. Phys. Eng. Sci."},{"key":"B23","doi-asserted-by":"publisher","first-page":"187","DOI":"10.1016\/j.mad.2004.09.025","article-title":"Integration from proteins to organs: The IUPS Physiome Project","volume":"126","author":"Hunter","year":"2005","journal-title":"Mech. Ageing Dev."},{"key":"B24","doi-asserted-by":"publisher","first-page":"108328","DOI":"10.1016\/j.mbs.2020.108328","article-title":"An efficient procedure to assist in the re-parametrization of structurally unidentifiable models","volume":"323","author":"Joubert","year":"2020","journal-title":"Math. Biosci."},{"key":"B25","doi-asserted-by":"publisher","first-page":"2401","DOI":"10.1002\/sim.2112","article-title":"How vague is vague? A simulation study of the impact of the use of vague prior distributions in MCMC using WinBUGS","volume":"24","author":"Lambert","year":"2005","journal-title":"Stat. Med."},{"key":"B26","doi-asserted-by":"publisher","first-page":"233","DOI":"10.1093\/ehjci\/jev014","article-title":"Recommendations for cardiac chamber quantification by echocardiography in adults: An update from the American society of echocardiography and the European association of cardiovascular imaging","volume":"16","author":"Lang","year":"2015","journal-title":"Eur. Heart J. Cardiovasc. Imaging"},{"key":"B27","doi-asserted-by":"publisher","first-page":"901","DOI":"10.1111\/echo.15073","article-title":"Variability of echocardiographic measures of left ventricular diastolic function. The HUNT study","volume":"38","author":"Letnes","year":"2021","journal-title":"Echocardiography"},{"key":"B28","doi-asserted-by":"publisher","first-page":"743","DOI":"10.1007\/s11517-009-0449-9","article-title":"Multi-scale modeling of the human cardiovascular system with applications to aortic valvular and arterial stenoses","volume":"47","author":"Liang","year":"2009","journal-title":"Med. Biol. Eng. Comput."},{"key":"B29","doi-asserted-by":"publisher","first-page":"1743","DOI":"10.1097\/HJH.0b013e3283567b94","article-title":"Evaluation of a novel sphygmomanometer, which estimates central aortic blood pressure from analysis of brachial artery suprasystolic pressure waves","volume":"30","author":"Lin","year":"2012","journal-title":"J. Hypertens."},{"key":"B30","doi-asserted-by":"publisher","first-page":"3","DOI":"10.1137\/090757009","article-title":"On identifiability of nonlinear ODE models and applications in viral dynamics","volume":"53","author":"Miao","year":"2011","journal-title":"SIAM Rev. Soc. Ind. Appl. Math."},{"key":"B31","doi-asserted-by":"publisher","first-page":"626","DOI":"10.1002\/cnm.1466","article-title":"A simple, versatile valve model for use in lumped parameter and one-dimensional cardiovascular models","volume":"28","author":"Mynard","year":"2012","journal-title":"Int. J. Numer. Method. Biomed. Eng."},{"key":"B32","doi-asserted-by":"publisher","first-page":"20200886","DOI":"10.1098\/rsif.2020.0886","article-title":"Assessing model mismatch and model selection in a Bayesian uncertainty quantification analysis of a fluid-dynamics model of pulmonary blood circulation","volume":"17","author":"Paun","year":"2020","journal-title":"J. R. Soc. Interface"},{"key":"B33","doi-asserted-by":"publisher","first-page":"1923","DOI":"10.1093\/bioinformatics\/btp358","article-title":"Structural and practical identifiability analysis of partially observed dynamical models by exploiting the profile likelihood","volume":"25","author":"Raue","year":"2009","journal-title":"Bioinformatics"},{"key":"B34","doi-asserted-by":"publisher","first-page":"148","DOI":"10.3389\/fphys.2018.00148","article-title":"Bond graph model of cerebral circulation: Toward clinically feasible systemic blood flow simulations","volume":"9","author":"Safaei","year":"2018","journal-title":"Front. Physiol."},{"key":"B35","first-page":"257","article-title":"Cvodes, the sensitivity-enabled ode solver in sundials","author":"Serban","year":"2005"},{"key":"B36","doi-asserted-by":"publisher","first-page":"20200055","DOI":"10.1098\/rsif.2020.0055","article-title":"Practical parameter identifiability for spatio-temporal models of cell invasion","volume":"17","author":"Simpson","year":"2020","journal-title":"J. R. Soc. Interface"},{"key":"B37","doi-asserted-by":"publisher","first-page":"1733","DOI":"10.1152\/ajpheart.1998.275.5.h1733","article-title":"Interaction between carotid baroregulation and the pulsating heart: A mathematical model","volume":"275","author":"Ursino","year":"1998","journal-title":"Am. J. Physiol."},{"key":"B38","doi-asserted-by":"publisher","first-page":"H149","DOI":"10.1152\/ajpheart.2000.279.1.h149","article-title":"Acute cardiovascular response to isocapnic hypoxia. I. A mathematical model","volume":"279","author":"Ursino","year":"2000","journal-title":"Am. J. Physiol. Heart Circ. Physiol."},{"key":"B39","doi-asserted-by":"publisher","first-page":"1479","DOI":"10.1152\/ajpheart.00850.2002","article-title":"Role of short-term cardiovascular regulation in heart period variability: A modeling study","volume":"284","author":"Ursino","year":"2003","journal-title":"Am. J. Physiol. Heart Circ. Physiol."},{"key":"B40","doi-asserted-by":"publisher","first-page":"656","DOI":"10.1093\/ejechocard\/jer107","article-title":"Right ventricular quantification in clinical practice: Two-dimensional vs. three-dimensional echocardiography compared with cardiac magnetic resonance imaging","volume":"12","author":"Van Der Zwaan","year":"2011","journal-title":"Eur. J. Echocardiogr."},{"key":"B41","doi-asserted-by":"publisher","first-page":"1130","DOI":"10.1093\/bioinformatics\/bts088","article-title":"An integrated strategy for prediction uncertainty analysis","volume":"28","author":"Vanlier","year":"2012","journal-title":"Bioinformatics"},{"key":"B42","doi-asserted-by":"publisher","first-page":"20130505","DOI":"10.1098\/rsif.2013.0505","article-title":"Reverse engineering and identification in systems biology: Strategies, perspectives and challenges","volume":"11","author":"Villaverde","year":"2014","journal-title":"J. R. Soc. Interface"},{"key":"B43","doi-asserted-by":"crossref","DOI":"10.1201\/b13568","volume-title":"McDonald\u2019s blood flow in arteries: Theoretical, experimental and clinical principles","author":"Vlachopoulos","year":"2011"},{"key":"B44","volume-title":"Applied biofluid mechanics","author":"Waite","year":"2007"},{"key":"B45","doi-asserted-by":"publisher","first-page":"961","DOI":"10.1051\/m2an\/2012053","article-title":"Mathematical model of blood flow in an anatomically detailed arterial network of the arm","volume":"47","author":"Watanabe","year":"2013","journal-title":"ESAIM Math. Model. Numer. Analysis"},{"key":"B46","doi-asserted-by":"publisher","first-page":"1543","DOI":"10.1002\/jmri.26564","article-title":"Test\u2013retest multisite reproducibility of neurovascular 4D flow MRI","volume":"49","author":"Wen","year":"2019","journal-title":"J. Magn. Reson. Imaging"},{"key":"B47","doi-asserted-by":"publisher","first-page":"357","DOI":"10.1016\/0021-9290(70)90036-9","article-title":"Arterial viscoelasticity: A generalized model. Effect on input impedance and wave travel in the systematic tree","volume":"3","author":"Westerhof","year":"1970","journal-title":"J. Biomech."},{"key":"B48","doi-asserted-by":"publisher","first-page":"60","DOI":"10.1016\/j.coisb.2021.03.005","article-title":"On structural and practical identifiability","volume":"25","author":"Wieland","year":"2021","journal-title":"Curr. Opin. Syst. Biol."},{"key":"B49","doi-asserted-by":"publisher","first-page":"728205","DOI":"10.3389\/fcvm.2021.728205","article-title":"Systematic comparison of left ventricular geometry between 3D-echocardiography and cardiac magnetic resonance imaging","volume":"8","author":"Zhao","year":"2021","journal-title":"Front. Cardiovasc. Med."}],"container-title":["Frontiers in Physiology"],"original-title":[],"link":[{"URL":"https:\/\/www.frontiersin.org\/articles\/10.3389\/fphys.2022.1018134\/full","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2023,11,30]],"date-time":"2023-11-30T13:37:12Z","timestamp":1701351432000},"score":1,"resource":{"primary":{"URL":"https:\/\/www.frontiersin.org\/articles\/10.3389\/fphys.2022.1018134\/full"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2022,11,9]]},"references-count":49,"alternative-id":["10.3389\/fphys.2022.1018134"],"URL":"https:\/\/doi.org\/10.3389\/fphys.2022.1018134","relation":{},"ISSN":["1664-042X"],"issn-type":[{"value":"1664-042X","type":"electronic"}],"subject":[],"published":{"date-parts":[[2022,11,9]]},"article-number":"1018134"}}