Translating myosin-binding protein C and titin abnormalities to whole-heart function using a novel calcium-contraction coupling model.

Mutations in cardiac myosin-binding protein C (cMyBP-C) or titin may respectively lead to hypertrophic (HCM) or dilated (DCM) cardiomyopathies. The mechanisms leading to these phenotypes remain unclear because of the challenge of translating cellular abnormalities to whole-heart and system function. We developed and validated a novel computer model of calcium-contraction coupling incorporating the role of cMyBP-C and titin based on the key assumptions: 1) tension in the thick filament promotes cross-bridge attachment mechanochemically, 2) with increasing titin tension, more myosin heads are unlocked for attachment, and 3) cMyBP-C suppresses cross-bridge attachment. Simulated stationary calcium-tension curves, isotonic and isometric contractions, and quick release agreed with experimental data. The model predicted that a loss of cMyBP-C function decreases the steepness of the calcium-tension curve, and that more compliant titin decreases the level of passive and active tension and its dependency on sarcomere length. Integrating this cellular model in the CircAdapt model of the human heart and circulation showed that a loss of cMyBP-C function resulted in HCM-like hemodynamics with higher left ventricular end-diastolic pressures and smaller volumes. More compliant titin led to higher diastolic pressures and ventricular dilation, suggesting DCM-like hemodynamics. The novel model of calcium-contraction coupling incorporates the role of cMyBP-C and titin. Its coupling to whole-heart mechanics translates changes in cellular calcium-contraction coupling to changes in cardiac pump and circulatory function and identifies potential mechanisms by which cMyBP-C and titin abnormalities may develop into HCM and DCM phenotypes. This modeling platform may help identify distinct mechanisms underlying clinical phenotypes in cardiac diseases.

Diagnostic value of reversed differential cyanosis in (supra)cardiac total anomalous pulmonary venous return.

BACKGROUND: To investigate the occurrence of reversed differential cyanosis (RDC) in case of (supra)cardiac total anomalous pulmonary venous return (TAPVR), we explored the hemodynamic changes and oxygen saturation levels during the fetal-to-neonatal transition in (supra)cardiac TAPVR, thereby revealing determinant factors of RDC. METHODS: A computational model was used to simulate the cardiovascular fetal-to-neonatal transition up to 24 h after birth. Abnormalities associated with TAPVR, like patent ductus arteriosus (PDA) and persistent pulmonary hypertension of the neonate (PPHN), were imposed on the model. Hemodynamic impact on flow distribution and right-sided pressures as well as oxygen saturations were assessed. RESULTS: Model findings demonstrated that RDC in (supra)cardiac TAPVR was dependent on two key factors: (1) the type of pulmonary venous connection being supracardiac or cardiac, and (2) the presence of a patent ductus arteriosus exhibiting right-to-left shunting. Persistence of RDC was mainly determined by the latter; an increase in pulmonary-to-systemic pressure difference by PPHN or PDA-induced pulmonary over-circulation contributed to persistence of RDC. CONCLUSION: This study highlights the significance of RDC in (supra)cardiac TAPVR and suggests to incorporate early screening ( < 24 h after birth) and to consider RDC as an immediate fail in screening protocols to ensure prompt detection of (supra)cardiac TAPVR. IMPACT: Utilizing a validated computational model for the cardiovascular fetal-to-neonatal transition, this study sheds light on the complex hemodynamics in neonates with (supra)cardiac Total Anomalous Pulmonary Venous Return (TAPVR). Model findings suggest that the often-present pulmonary over-circulation in neonates with TAPVR might significantly contribute to the anomaly's frequent omission during pulse-oximetry screening beyond the first 24 h after birth. This study highlights the diagnostic value of reversed differential cyanosis in early screenings within the first 24 h after birth. By including RDC as an immediate fail in early pulse-oximetry screenings, the likelihood of missing (supra)cardiac TAPVR cases could be reduced.

How to assess and treat right ventricular electromechanical dyssynchrony in post-repair tetralogy of Fallot: insights from imaging, invasive studies, and computational modelling.

BACKGROUND AND AIMS: Right bundle branch block (RBBB) and resulting right ventricular (RV) electromechanical discoordination are thought to play a role in the disease process of subpulmonary RV dysfunction that frequently occur post-repair tetralogy of Fallot (ToF). We sought to describe this disease entity, the role of pulmonary re-valvulation, and the potential added value of RV cardiac resynchronization therapy (RV-CRT). METHODS: Two patients with repaired ToF, complete RBBB, pulmonary regurgitation, and significantly decreased RV function underwent echocardiography, cardiac magnetic resonance, and an invasive study to evaluate the potential for RV-CRT as part of the management strategy. The data were used to personalize the CircAdapt model of the human heart and circulation. Resulting Digital Twins were analysed to quantify the relative effects of RV pressure and volume overload and to predict the effect of RV-CRT. RESULTS: Echocardiography showed components of a classic RV dyssynchrony pattern which could be reversed by RV-CRT during invasive study and resulted in acute improvement in RV systolic function. The Digital Twins confirmed a contribution of electromechanical RV dyssynchrony to RV dysfunction and suggested improvement of RV contraction efficiency after RV-CRT. The one patient who underwent successful permanent RV-CRT as part of the pulmonary re-valvulation procedure carried improvements that were in line with the predictions based on his Digital Twin. CONCLUSION: An integrative diagnostic approach to RV dysfunction, including the construction of Digital Twins may help to identify candidates for RV-CRT as part of the lifetime management of ToF and similar congenital heart lesions.

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.

Exploring QRS Area beyond Patient Selection in CRT-Can It Guide Left Ventricular Lead Placement?

Vectorcardiographic QRS area is a promising tool for patient selection and implantation guidance in cardiac resynchronization therapy (CRT). Research has mainly focused on the role of QRS area in patient selection for CRT. Recently, QRS area has been proposed as a tool to guide left ventricular lead placement in CRT. Theoretically, vector-based electrical information of ventricular fusion pacing, calculated from the basic 12-lead ECG, can give real-time insight into the extent of resynchronization at any LV lead position, as well as any selected electrode on the LV lead. The objective of this review is to provide an overview of the background of vectorcardiographic QRS area and its potential in optimizing LV lead location in order to optimize the benefits of CRT.

Mitral regurgitation in heart failure with preserved ejection fraction: The interplay of valve, ventricle, and atrium.

Mitral regurgitation (MR) is highly prevalent among patients with heart failure and preserved ejection fraction (HFpEF). Despite this combination being closely associated with unfavourable outcomes, it remains relatively understudied. This is partly due to the inherent heterogeneity of patients with HFpEF. To address this gap, dissecting HFpEF into mechanism-based phenotypes may offer a promising avenue for advancing our comprehension of these complex intertwined conditions. This review employs the validated CircAdapt model to explore the haemodynamic implications of moderate to severe MR across a well-defined spectrum of myocardial disease, characterized by impaired relaxation and reduced myocardial compliance. Both heart failure and mitral valve disease share overlapping symptomatology, primarily attributed to elevated pulmonary pressures. The intricate mechanisms contributing to these elevated pressures are multifaceted, potentially influenced by diastolic dysfunction, left atrial myopathy, and MR. Accurate evaluation of the haemodynamic and clinical impact of MR necessitates a comprehensive approach, taking into account the characteristics of both the left atrium and left ventricle, as well as their intricate interactions, which may currently be underemphasized in diagnostic practice. This holistic assessment is imperative for enhancing our understanding and refining therapeutic strategies within this patient cohort.

Increased Chamber Resting Tone Is a Key Determinant of Left Ventricular Diastolic Dysfunction.

BACKGROUND: Twitch-independent tension has been demonstrated in cardiomyocytes, but its role in heart failure (HF) is unclear. We aimed to address twitch-independent tension as a source of diastolic dysfunction by isolating the effects of chamber resting tone (RT) from impaired relaxation and stiffness. METHODS: We invasively monitored pressure-volume data during cardiopulmonary exercise in 20 patients with hypertrophic cardiomyopathy, 17 control subjects, and 35 patients with HF with preserved ejection fraction. To measure RT, we developed a new method to fit continuous pressure-volume measurements, and first validated it in a computational model of loss of cMyBP-C (myosin binding protein-C). RESULTS: In hypertrophic cardiomyopathy, RT (estimated marginal mean [95% CI]) was 3.4 (0.4-6.4) mm Hg, increasing to 18.5 (15.5-21.5) mm Hg with exercise (P<0.001). At peak exercise, RT was responsible for 64% (53%-76%) of end-diastolic pressure, whereas incomplete relaxation and stiffness accounted for the rest. RT correlated with the levels of NT-proBNP (N-terminal pro-B-type natriuretic peptide; R=0.57; P=0.02) and with pulmonary wedge pressure but following different slopes at rest and during exercise (R2=0.49; P<0.001). In controls, RT was 0.0 mm Hg and 1.2 (0.3-2.8) mm Hg in HF with preserved ejection fraction patients and was also exacerbated by exercise. In silico, RT increased in parallel to the loss of cMyBP-C function and correlated with twitch-independent myofilament tension (R=0.997). CONCLUSIONS: Augmented RT is the major cause of LV diastolic chamber dysfunction in hypertrophic cardiomyopathy and HF with preserved ejection fraction. RT transients determine diastolic pressures, pulmonary pressures, and functional capacity to a greater extent than relaxation and stiffness abnormalities. These findings support antimyosin agents for treating HF.