Linking cross-bridge cycling kinetics to response to cardiac resynchronization therapy: a multiscale modelling study

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.

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.

Advancing clinical translation of cardiac biomechanics models: a comprehensive review, applications and future pathways

Cardiac mechanics models are developed to represent a high level of detail, including refined anatomies, accurate cell mechanics models, and platforms to link microscale physiology to whole-organ function. However, cardiac biomechanics models still have limited clinical translation. In this review, we provide a picture of cardiac mechanics models, focusing on their clinical translation. We review the main experimental and clinical data used in cardiac models, as well as the steps followed in the literature to generate anatomical meshes ready for simulations. We describe the main models in active and passive mechanics and the different lumped parameter models to represent the circulatory system. Lastly, we provide a summary of the state-of-the-art in terms of ventricular, atrial, and four-chamber cardiac biomechanics models. We discuss the steps that may facilitate clinical translation of the biomechanics models we describe. A well-established software to simulate cardiac biomechanics is lacking, with all available platforms involving different levels of documentation, learning curves, accessibility, and cost. Furthermore, there is no regulatory framework that clearly outlines the verification and validation requirements a model has to satisfy in order to be reliably used in applications. Finally, better integration with increasingly rich clinical and/or experimental datasets as well as machine learning techniques to reduce computational costs might increase model reliability at feasible resources. Cardiac biomechanics models provide excellent opportunities to be integrated into clinical workflows, but more refinement and careful validation against clinical data are needed to improve their credibility. In addition, in each context of use, model complexity must be balanced with the associated high computational cost of running these models.

Timing of cardiac resynchronization therapy implantation

AIMS

The optimum timing of cardiac resynchronization therapy (CRT) implantation is unknown. We explored long-term outcomes after CRT in relation to the time interval from a first heart failure hospitalization (HFH) to device implantation.

METHODS AND RESULTS

A database covering the population of England (56.3 million in 2019) was used to quantify clinical outcomes after CRT implantation in relation to first HFHs. From 2010 to 2019, 64 968 patients [age: 71.4 ± 11.7 years; 48 606 (74.8%) male] underwent CRT implantation, 57% in the absence of a previous HFH, 12.9% during the first HFH, and 30.1% after ≥1 HFH. Over 4.54 (2.80-6.71) years [median (interquartile range); 272 989 person-years], the time in years from the first HFH to CRT implantation was associated with a higher risk of total mortality [hazard ratio (HR); 95% confidence intervals (95% CI)] (1.15; 95% CI 1.14-1.16, HFH (HR: 1.26; 95% CI 1.24-1.28), and the combined endpoint of total mortality or HFH (HR: 1.19; 95% CI 1.27-1.20) than CRT in patients with no previous HFHs, after co-variate adjustment. Total mortality (HR: 1.67), HFH (HR: 2.63), and total mortality or HFH (HR: 1.92) (all P < 0.001) were highest in patients undergoing CRT ≥2 years after the first HFH.

CONCLUSION

In this study of a healthcare system covering an entire nation, delays from a first HFH to CRT implantation were associated with progressively worse long-term clinical outcomes. The best clinical outcomes were observed in patients with no previous HFH and in those undergoing CRT implantation during the first HFH.

CONDENSED ABSTRACT

The optimum timing of CRT implantation is unknown. In this study of 64 968 consecutive patients, delays from a first heart failure hospitalization (HFH) to CRT implantation were associated with progressively worse long-term clinical outcomes. Each year from a first HFH to CRT implantation was associated with a 21% higher risk of total mortality and a 34% higher risk of HFH. The best outcomes after CRT were observed in patients with no previous HFHs and in those undergoing implantation during their first HFH.

A multiscale model of the cardiovascular system that regulates arterial pressure via closed loop baroreflex control of chronotropism, cell-level contractility, and vascular tone

Multiscale models of the cardiovascular system can provide new insights into physiological and pathological processes. PyMyoVent is a computer model that bridges from molecular- to organ-level function and which simulates a left ventricle pumping blood through the systemic circulation. Initial work with PyMyoVent focused on the end-systolic pressure volume relationship and ranked potential therapeutic strategies by their impact on contractility. This manuscript extends the PyMyoVent framework by adding closed-loop feedback control of arterial pressure. The control algorithm mimics important features of the physiological baroreflex and was developed as part of a long-term program that focuses on growth and biological remodeling. Inspired by the underlying biology, the reflex algorithm uses an afferent signal derived from arterial pressure to drive a kinetic model that mimics the net result of neural processing in the medulla and cell-level responses to autonomic drive. The kinetic model outputs control signals that are constrained between limits that represent maximum parasympathetic and maximum sympathetic drive and which modulate heart rate, intracellular Ca2+ dynamics, the molecular-level function of both the thick and the thin myofilaments, and vascular tone. Simulations show that the algorithm can regulate mean arterial pressure at user-defined setpoints as well as maintaining arterial pressure when challenged by changes in blood volume and/or valve resistance. The reflex also regulates arterial pressure when cell-level contractility is modulated to mimic the idealized impact of myotropes. These capabilities will be important for future work that uses computer modeling to investigate clinical conditions and treatments.

Differentiating the effects of β-adrenergic stimulation and stretch on calcium and force dynamics using a novel electromechanical cardiomyocyte model

Cardiac electrophysiology and mechanics are strongly interconnected. Calcium is crucial in this complex interplay through its role in cellular electrophysiology and sarcomere contraction. We aim to differentiate the effects of acute β-adrenergic stimulation (β-ARS) and cardiomyocyte stretch (increased sarcomere length) on calcium-transient dynamics and force generation, using a novel computational model of cardiac electromechanics. We implemented a bidirectional coupling between the O'Hara-Rudy model of human ventricular electrophysiology and the MechChem model of sarcomere mechanics through the buffering of calcium by troponin. The coupled model was validated using experimental data from large mammals or human samples. Calcium transient and force were simulated for various degrees of β-ARS and initial sarcomere lengths. The model reproduced force-frequency, quick-release, and isotonic contraction experiments, validating the bidirectional electromechanical interactions. An increase in β-ARS increased the amplitudes of force (augmented inotropy) and calcium transient, and shortened both force and calcium-transient duration (lusitropy). An increase in sarcomere length increased force amplitude even more, but decreased calcium-transient amplitude and increased both force and calcium-transient duration. Finally, a gradient in relaxation along the thin filament may explain the nonmonotonic decay in cytosolic calcium observed with high tension. Using a novel coupled human electromechanical model, we identified differential effects of β-ARS and stretch on calcium and force. Stretch mostly contributed to increased force amplitude and β-ARS to the reduction of calcium and force duration. We showed that their combination, rather than individual contributions, is key to ensure force generation, rapid relaxation, and low diastolic calcium levels.NEW & NOTEWORTHY This work identifies the contribution of electrical and mechanical alterations to regulation of calcium and force under exercise-like conditions using a novel human electromechanical model integrating ventricular electrophysiology and sarcomere mechanics. By better understanding their individual and combined effects, this can uncover arrhythmogenic mechanisms in exercise-like situations. This publicly available model is a crucial step toward understanding the complex interplay between cardiac electrophysiology and mechanics to improve arrhythmia risk prediction and treatment.

Mechano-calcium and mechano-electric feedbacks in the human cardiomyocyte analyzed in a mathematical model

Experiments on animal hearts (rat, rabbit, guinea pig, etc.) have demonstrated that mechano-calcium feedback (MCF) and mechano-electric feedback (MEF) are very important for myocardial self-regulation because they adjust the cardiomyocyte contractile function to various mechanical loads and to mechanical interactions between heterogeneous myocardial segments in the ventricle walls. In in vitro experiments on these animals, MCF and MEF manifested themselves in several basic classical phenomena (e.g., load dependence, length dependence of isometric twitches, etc.), and in the respective responses of calcium transients and action potentials. However, it is extremely difficult to study simultaneously the electrical, calcium, and mechanical activities of the human heart muscle in vitro. Mathematical modeling is a useful tool for exploring these phenomena. We have developed a novel model to describe electromechanical coupling and mechano-electric feedbacks in the human cardiomyocyte. It combines the ‘ten Tusscher–Panfilov’ electrophysiological model of the human cardiomyocyte with our module of myocardium mechanical activity taken from the ‘Ekaterinburg–Oxford’ model and adjusted to human data. Using it, we simulated isometric and afterloaded twitches and effects of MCF and MEF on excitation–contraction coupling. MCF and MEF were found to affect significantly the duration of the calcium transient and action potential in the human cardiomyocyte model in response to both smaller afterloads as compared to bigger ones and various mechanical interventions applied during isometric and afterloaded twitches.