Integrating Imaging and Invasive Pressure Data into a Multiscale Whole-Heart Model

Cardiovascular diseases are the leading cause of death. Clinical data used to decide treatment are hard to integrate and interpret, making optimal treatment selection difficult. Personalized models can be used to integrate clinical data into a physics and physiology-constrained framework, but their clinical application faces limitations due to complex calibration and validation. In this study, we present a novel systematic calibration method for a whole-heart, multiscale, electromechanics model using emulators, sensitivity analysis, and history matching. Using cardiac motion derived from ECG-gated computed tomography (CT) and invasive left ventricular (LV) pressure data, we calibrated 25 model parameters to match the LV end-diastolic (ED) and peak pressure, ED and end-systolic (ES) volumes (EDV and ESV), right ventricle EDV, and the left atrium EDV, ESV, and the maximum volume during venous return. After calibration, all features were fit within [0.8, 10.8]% of the mean target value, and fell within 1.4 experimental standard deviations from the target values. We validated the model by comparing CT-derived and simulated atrioventricular plane displacement (AVPD) (8.2 versus 8.1 mm) and the ED and ES configurations against the CT images. The model replicated the measured acute hemodynamic response to biventricular (BIV) pacing (simulated: 222 mmHg/s versus clinical: 213±65 mmHg/s). This study provides a systematic method to integrate clinical data into a whole-heart, multiscale electromechanics framework. The validation shows that the model replicates local heart motion and response to therapy, demonstrating potential in assisting clinical decision-making.