Parameter subset reduction for imaging-based digital twin generation of patients with left ventricular mechanical discoordination.

BACKGROUND: Integration of a patient's non-invasive imaging data in a digital twin (DT) of the heart can provide valuable insight into the myocardial disease substrates underlying left ventricular (LV) mechanical discoordination. However, when generating a DT, model parameters should be identifiable to obtain robust parameter estimations. In this study, we used the CircAdapt model of the human heart and circulation to find a subset of parameters which were identifiable from LV cavity volume and regional strain measurements of patients with different substrates of left bundle branch block (LBBB) and myocardial infarction (MI). To this end, we included seven patients with heart failure with reduced ejection fraction (HFrEF) and LBBB (study ID: 2018-0863, registration date: 2019-10-07), of which four were non-ischemic (LBBB-only) and three had previous MI (LBBB-MI), and six narrow QRS patients with MI (MI-only) (study ID: NL45241.041.13, registration date: 2013-11-12). Morris screening method (MSM) was applied first to find parameters which were important for LV volume, regional strain, and strain rate indices. Second, this parameter subset was iteratively reduced based on parameter identifiability and reproducibility. Parameter identifiability was based on the diaphony calculated from quasi-Monte Carlo simulations and reproducibility was based on the intraclass correlation coefficient ( ICC ) obtained from repeated parameter estimation using dynamic multi-swarm particle swarm optimization. Goodness-of-fit was defined as the mean squared error ( χ 2 ) of LV myocardial strain, strain rate, and cavity volume. RESULTS: A subset of 270 parameters remained after MSM which produced high-quality DTs of all patients ( χ 2 < 1.6), but minimum parameter reproducibility was poor ( ICC min = 0.01). Iterative reduction yielded a reproducible ( ICC min = 0.83) subset of 75 parameters, including cardiac output, global LV activation duration, regional mechanical activation delay, and regional LV myocardial constitutive properties. This reduced subset produced patient-resembling DTs ( χ 2 < 2.2), while septal-to-lateral wall workload imbalance was higher for the LBBB-only DTs than for the MI-only DTs (p < 0.05). CONCLUSIONS: By applying sensitivity and identifiability analysis, we successfully determined a parameter subset of the CircAdapt model which can be used to generate imaging-based DTs of patients with LV mechanical discoordination. Parameters were reproducibly estimated using particle swarm optimization, and derived LV myocardial work distribution was representative for the patient's underlying disease substrate. This DT technology enables patient-specific substrate characterization and can potentially be used to support clinical decision making.

Virtual pacing of a patient's digital twin to predict left ventricular reverse remodelling after cardiac resynchronization therapy.

AIMS: Identifying heart failure (HF) patients who will benefit from cardiac resynchronization therapy (CRT) remains challenging. We evaluated whether virtual pacing in a digital twin (DT) of the patient's heart could be used to predict the degree of left ventricular (LV) reverse remodelling post-CRT. METHODS AND RESULTS: Forty-five HF patients with wide QRS complex (≥130 ms) and reduced LV ejection fraction (≤35%) receiving CRT were retrospectively enrolled. Echocardiography was performed before (baseline) and 6 months after CRT implantation to obtain LV volumes and 18-segment longitudinal strain. A previously developed algorithm was used to generate 45 DTs by personalizing the CircAdapt model to each patient's baseline measurements. From each DT, baseline septal-to-lateral myocardial work difference (MWLW-S,DT) and maximum rate of LV systolic pressure rise (dP/dtmax,DT) were derived. Biventricular pacing was then simulated using patient-specific atrioventricular delay and lead location. Virtual pacing-induced changes ΔMWLW-S,DT and ΔdP/dtmax,DT were correlated with real-world LV end-systolic volume change at 6-month follow-up (ΔLVESV). The DT's baseline MWLW-S,DT and virtual pacing-induced ΔMWLW-S,DT were both significantly associated with the real patient's reverse remodelling ΔLVESV (r = -0.60, P < 0.001 and r = 0.62, P < 0.001, respectively), while correlation between ΔdP/dtmax,DT and ΔLVESV was considerably weaker (r = -0.34, P = 0.02). CONCLUSION: Our results suggest that the reduction of septal-to-lateral work imbalance by virtual pacing in the DT can predict real-world post-CRT LV reverse remodelling. This DT approach could prove to be an additional tool in selecting HF patients for CRT and has the potential to provide valuable insights in optimization of CRT delivery.

A Lumped Two-Compartment Model for Simulation of Ventricular Pump and Tissue Mechanics in Ischemic Heart Disease.

Introduction: Computational modeling of cardiac mechanics and hemodynamics in ischemic heart disease (IHD) is important for a better understanding of the complex relations between ischemia-induced heterogeneity of myocardial tissue properties, regional tissue mechanics, and hemodynamic pump function. We validated and applied a lumped two-compartment modeling approach for IHD integrated into the CircAdapt model of the human heart and circulation. Methods: Ischemic contractile dysfunction was simulated by subdividing a left ventricular (LV) wall segment into a hypothetical contractile and noncontractile compartment, and dysfunction severity was determined by the noncontractile volume fraction ( N C V F ). Myocardial stiffness was determined by the zero-passive stress length ( L s 0 , p a s ) and nonlinearity ( k E C M ) of the passive stress-sarcomere length relation of the noncontractile compartment. Simulated end-systolic pressure volume relations (ESPVRs) for 20% acute ischemia were qualitatively compared between a two- and one-compartment simulation, and parameters of the two-compartment model were tuned to previously published canine data of regional myocardial deformation during acute and prolonged ischemia and reperfusion. In six patients with myocardial infarction (MI), the N C V F was automatically estimated using the echocardiographic LV strain and volume measurements obtained acutely and 6 months after MI. Estimated segmental N C V F values at the baseline and 6-month follow-up were compared with percentage late gadolinium enhancement (LGE) at 6-month follow-up. Results: Simulation of 20% of N C V F shifted the ESPVR rightward while moderately reducing the slope, while a one-compartment simulation caused a leftward shift with severe reduction in the slope. Through tuning of the N C V F , L s 0 , p a s , and k E C M , it was found that manipulation of the N C V F alone reproduced the deformation during acute ischemia and reperfusion, while additional manipulations of L s 0 , p a s and k E C M were required to reproduce deformation during prolonged ischemia and reperfusion. Out of all segments with LGE>25% at the follow-up, the majority (68%) had higher estimated N C V F at the baseline than at the follow-up. Furthermore, the baseline N C V F correlated better with percentage LGE than N C V F did at the follow-up. Conclusion: We successfully used a two-compartment model for simulation of the ventricular pump and tissue mechanics in IHD. Patient-specific optimizations using regional myocardial deformation estimated the N C V F in a small cohort of MI patients in the acute and chronic phase after MI, while estimated N C V F values closely approximated the extent of the myocardial scar at the follow-up. In future studies, this approach can facilitate deformation imaging-based estimation of myocardial tissue properties in patients with cardiovascular diseases.

Uncertainty Quantification of Regional Cardiac Tissue Properties in Arrhythmogenic Cardiomyopathy Using Adaptive Multiple Importance Sampling.

Introduction: Computational models of the cardiovascular system are widely used to simulate cardiac (dys)function. Personalization of such models for patient-specific simulation of cardiac function remains challenging. Measurement uncertainty affects accuracy of parameter estimations. In this study, we present a methodology for patient-specific estimation and uncertainty quantification of parameters in the closed-loop CircAdapt model of the human heart and circulation using echocardiographic deformation imaging. Based on patient-specific estimated parameters we aim to reveal the mechanical substrate underlying deformation abnormalities in patients with arrhythmogenic cardiomyopathy (AC). Methods: We used adaptive multiple importance sampling to estimate the posterior distribution of regional myocardial tissue properties. This methodology is implemented in the CircAdapt cardiovascular modeling platform and applied to estimate active and passive tissue properties underlying regional deformation patterns, left ventricular volumes, and right ventricular diameter. First, we tested the accuracy of this method and its inter- and intraobserver variability using nine datasets obtained in AC patients. Second, we tested the trueness of the estimation using nine in silico generated virtual patient datasets representative for various stages of AC. Finally, we applied this method to two longitudinal series of echocardiograms of two pathogenic mutation carriers without established myocardial disease at baseline. Results: Tissue characteristics of virtual patients were accurately estimated with a highest density interval containing the true parameter value of 9% (95% CI [0-79]). Variances of estimated posterior distributions in patient data and virtual data were comparable, supporting the reliability of the patient estimations. Estimations were highly reproducible with an overlap in posterior distributions of 89.9% (95% CI [60.1-95.9]). Clinically measured deformation, ejection fraction, and end-diastolic volume were accurately simulated. In presence of worsening of deformation over time, estimated tissue properties also revealed functional deterioration. Conclusion: This method facilitates patient-specific simulation-based estimation of regional ventricular tissue properties from non-invasive imaging data, taking into account both measurement and model uncertainties. Two proof-of-principle case studies suggested that this cardiac digital twin technology enables quantitative monitoring of AC disease progression in early stages of disease.

Electromechanical substrate characterization in arrhythmogenic cardiomyopathy using imaging-based patient-specific computer simulations.

AIMS: Arrhythmogenic cardiomyopathy (AC) is an inherited cardiac disease, characterized by life-threatening ventricular arrhythmias and progressive cardiac dysfunction. The aim of this study is to use computer simulations to non-invasively estimate the individual patient's myocardial tissue substrates underlying regional right ventricular (RV) deformation abnormalities in a cohort of AC mutation carriers. METHODS AND RESULTS: In 68 AC mutation carriers and 20 control subjects, regional longitudinal deformation patterns of the RV free wall (RVfw), interventricular septum (IVS), and left ventricular free wall (LVfw) were obtained using speckle-tracking echocardiography. We developed and used a patient-specific parameter estimation protocol based on the multi-scale CircAdapt cardiovascular system model to create virtual AC subjects. Using the individual's deformation data as model input, this protocol automatically estimated regional RVfw and global IVS and LVfw tissue properties. The computational model was able to reproduce clinically measured regional deformation patterns for all subjects, with highly reproducible parameter estimations. Simulations revealed that regional RVfw heterogeneity of both contractile function and compliance were increased in subjects with clinically advanced disease compared to mutation carriers without clinically established disease (17 ± 13% vs. 8 ± 4%, P = 0.01 and 18 ± 11% vs. 10 ± 7%, P < 0.01, respectively). No significant difference in activation delay was found. CONCLUSION: Regional RV deformation abnormalities in AC mutation carriers were related to reduced regional contractile function and tissue compliance. In clinically advanced disease stages, a characteristic apex-to-base heterogeneity of tissue abnormalities was present in the majority of the subjects, with most pronounced disease in the basal region of the RVfw.

Parameter subset reduction for patient-specific modelling of arrhythmogenic cardiomyopathy-related mutation carriers in the CircAdapt model.

Arrhythmogenic cardiomyopathy (AC) is an inherited cardiac disease, clinically characterized by life-threatening ventricular arrhythmias and progressive cardiac dysfunction. Patient-specific computational models could help understand the disease progression and may help in clinical decision-making. We propose an inverse modelling approach using the CircAdapt model to estimate patient-specific regional abnormalities in tissue properties in AC subjects. However, the number of parameters (n = 110) and their complex interactions make personalized parameter estimation challenging. The goal of this study is to develop a framework for parameter reduction and estimation combining Morris screening, quasi-Monte Carlo (qMC) simulations and particle swarm optimization (PSO). This framework identifies the best subset of tissue properties based on clinical measurements allowing patient-specific identification of right ventricular tissue abnormalities. We applied this framework on 15 AC genotype-positive subjects with varying degrees of myocardial disease. Cohort studies have shown that atypical regional right ventricular (RV) deformation patterns reveal an early-stage AC disease. The CircAdapt model of cardiovascular mechanics and haemodynamics has already demonstrated its ability to capture typical deformation patterns of AC subjects. We, therefore, use clinically measured cardiac deformation patterns to estimate model parameters describing myocardial disease substrates underlying these AC-related RV deformation abnormalities. Morris screening reduced the subset to 48 parameters. qMC and PSO further reduced the subset to a final selection of 16 parameters, including regional tissue contractility, passive stiffness, activation delay and wall reference area. This article is part of the theme issue 'Uncertainty quantification in cardiac and cardiovascular modelling and simulation'.

Improving long QT syndrome diagnosis by a polynomial-based T-wave morphology characterization.

BACKGROUND: Diagnosing long QT syndrome (LQTS) remains challenging because of a considerable overlap in QT interval between patients with LQTS and healthy subjects. Characterizing T-wave morphology might improve LQTS diagnosis. OBJECTIVE: The purpose of this study was to improve LQTS diagnosis by combining new polynomial-based T-wave morphology parameters with the corrected QT interval (QTc), age, and sex in a model. METHODS: A retrospective cohort consisting of 333 patients with LQTS and 345 genotype-negative family members was used in this study. For each patient, a linear combination of the first 2 Hermite-Gauss (HG) polynomials was fitted to the STT segments of an average complex of all precordial leads and limb leads I and II. The weight coefficients as well as the error of the best fit were used to characterize T-wave morphology. Subjects were classified as patients with LQTS or controls by clinical QTc cutoffs and 3 support vector machine models fed with different features. An external cohort consisting of 72 patients and 45 controls was finally used to check the robustness of the models. RESULTS: Baseline QTc cutoffs were specific but had low sensitivity in diagnosing LQTS. The model with T-wave morphology features, QTc, age, and sex had the best overall accuracy (84%), followed by a model with QTc, age, and sex (79%). The model with T-wave morphology features especially performed better in LQTS type 3 patients (69%). CONCLUSION: T-wave morphologies can be characterized by fitting a linear combination of the first 2 Hermite-Gauss polynomials. Adding T-wave morphology characterization to age, sex, and QTc in a support vector machine model improves LQTS diagnosis.