Mathematical models that integrate multi-scale physiological data can offer insight into physiological and pathophysiological function, and may eventually assist in individualized predictive medicine. We present a methodology for performing systematic analyses of multi-parameter interactions in such complex, multi-scale models. Human physiology models are often based on or inspired by Arthur Guyton's whole-body circulatory regulation model. Despite the significance of this model, it has not been the subject of a systematic and comprehensive sensitivity study. Therefore, we use this model as a case study for our methodology. Our analysis of the Guyton model reveals how the multitude of model parameters combine to affect the model dynamics, and how interesting combinations of parameters may be identified. It also includes a “virtual population” from which “virtual individuals” can be chosen, on the basis of exhibiting conditions similar to those of a real-world patient. This lays the groundwork for using the Guyton model for in silico exploration of pathophysiological states and treatment strategies. The results presented here illustrate several potential uses for the entire dataset of sensitivity results and the “virtual individuals” that we have generated, which are included in the supplementary material. More generally, the presented methodology is applicable to modern, more complex multi-scale physiological models.

This paper presents a contribution to the definition of the interfaces required to perform heterogeneous model integration in the context of integrative physiology. A formalization of the model integration problem is proposed and a coupling method is presented. The extension of the classic Guyton model, a multi-organ, integrated systems model of blood pressure regulation, is used as an example of the application of the proposed method. To this end, the Guyton model has been restructured, extensive sensitivity analyses have been performed, and appropriate transformations have been applied to replace a subset of its constituting modules by integrating a pulsatile heart and an updated representation of the renin-angiotensin system. Simulation results of the extended integrated model are presented and the impacts of their integration within the original model are evaluated.