High tension in sarcomeres hinders myocardial relaxation: A computational study

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.

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.

A Novel In Silico Electromechanical Model of Human Ventricular Cardiomyocyte

Contractility has become one of the main readouts in computational and experimental studies on cardiomyocytes. Following this trend, we propose a novel mathematical model of human ventricular cardiomyocytes electromechanics, BPSLand, by coupling a recent human contractile element to the BPS2020 model of electrophysiology. BPSLand is the result of a hybrid optimization process and it reproduces all the electrophysiology experimental indices captured by its predecessor BPS2020, simultaneously enabling the simulation of realistic human active tension and its potential abnormalities. The transmural heterogeneity in both electrophysiology and contractility departments was simulated consistent with previous computational and in vitro studies. Furthermore, our model could capture delayed afterdepolarizations (DADs), early afterdepolarizations (EADs), and contraction abnormalities in terms of aftercontractions triggered by either drug action or special pacing modes. Finally, we further validated the mechanical results of the model against previous experimental and in silico studies, e.g., the contractility dependence on pacing rate. Adding a new level of applicability to the normative models of human cardiomyocytes, BPSLand represents a robust, fully-human in silico model with promising capabilities for translational cardiology.

Reduced preload increases Mechanical Control (strain-rate dependence) of Relaxation by modifying myosin kinetics

Rapid myocardial relaxation is essential in maintaining cardiac output, and impaired relaxation is an early indicator of diastolic dysfunction. While the biochemical modifiers of relaxation are well known to include calcium handling, thin filament activation, and myosin kinetics, biophysical and biomechanical modifiers can also alter relaxation. We have previously shown that the relaxation rate is increased by an increasing strain rate, not a reduction in afterload. The slope of the relaxation rate to strain rate relationship defines Mechanical Control of Relaxation (MCR). To investigate MCR further, we performed in vitro experiments and computational modeling of preload-adjustment using intact rat cardiac trabeculae. Trabeculae studies are often performed using isometric (fixed-end) muscles at optimal length (Lo, length producing maximal developed force). We determined that reducing muscle length from Lo increased MCR by 20%, meaning that reducing preload could substantially increase the sensitivity of the relaxation rate to the strain rate. We subsequently used computational modeling to predict mechanisms that might underlie this preload-dependence. Computational modeling was not able to fully replicate experimental data, but suggested that thin-filament properties are not sufficient to explain preload-dependence of MCR because the model required the thin-filament to become more activated at reduced preloads. The models suggested that myosin kinetics may underlie the increase in MCR at reduced preload, an effect that can be enhanced by force-dependence. Relaxation can be modified and enhanced by reduced preload. Computational modeling implicates myosin-based targets for treatment of diastolic dysfunction, but further model refinements are needed to fully replicate experimental data.

Disruption of transverse-tubular network reduces energy efficiency in cardiac muscle contraction

AIM

Altered organization of the transverse-tubular network is an early pathological event occurring even prior to the onset of heart failure. Such t-tubular remodelling disturbs the synchrony and signalling between membranous and intracellular ion channels, exchangers, receptors and ATPases essential in the dynamics of excitation-contraction coupling, leading to ionic abnormality and mechanical dysfunction in heart disease progression. In this study, we investigated whether a disrupted t-tubular network has a direct effect on cardiac mechano-energetics. Our aim was to understand the fundamental link between t-tubular remodelling and impaired energy metabolism, both of which are characteristics of heart failure. We thus studied healthy tissue preparations in which cellular processes are not altered by any disease event.

METHODS

We exploited the "formamide-detubulation" technique to acutely disrupt the t-tubular network in rat left-ventricular trabeculae. We assessed the energy utilization by cellular Ca²⁺ cycling and by crossbridge cycling, and quantified the change of energy efficiency following detubulation. For these measurements, trabeculae were mounted in a microcalorimeter where force and heat output were simultaneously measured.

RESULTS

Following structural disorganization from detubulation, muscle heat output associated with Ca²⁺ cycling was reduced, indicating impaired intracellular Ca²⁺ homeostasis. This led to reduced force production and heat output by crossbridge cycling. The reduction in force-length work was not paralleled by proportionate reduction in the heat output and, as such, energy efficiency was reduced.

CONCLUSIONS

These results reveal the direct energetic consequences of disrupted t-tubular network, linking the energy disturbance and the t-tubular remodelling typically observed in heart failure.

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.

The role of pacing rate in the modulation of mechano-induced immediate and delayed changes in the force and Ca-transient of cardiac muscle

Myocardial function is tuned by dynamic changes in length and load via mechano-calcium feedback. This regulation may be significantly affected by heart rhythm. We evaluated the mechano-induced modulation of contractility and Ca-transient (CaT) in the rat myocardium subjected to twitch-by-twitch shortening-re-lengthening (↓-↑) trains of different lengths (N = 1 … 720 cycles) at low (1 Hz) and near-physiological (3.5 Hz) pacing rates. Force/CaT characteristics were evaluated in the first post-train isometric twitch (immediate effect) and during slow changes (delayed maximal elevation/decrease) and compared with those of the pre-train twitch. The immediate inotropic effect was positive for N = 30 … 720 and negative for N = 1 … 20, while the delayed effect was always positive. The immediate and delayed inotropic effects were significantly higher at 3.5-Hz vs 1-Hz (P < 0.05). The prominent inotropism was accompanied by much smaller changes in the CaT diastolic level/amplitude. The shortening-re-lengthening train induced oscillations of the slow change in force at 3.5-Hz (always) and at 1-Hz (∼50% of muscles), which were dependent of the train length and independent of the pacing rate. We suggest that twitch-by-twitch shortening-re-lengthening of cardiac muscle decreases Ca²⁺ buffering by troponin C and elevates Ca²⁺ loading of the sarcoplasmic reticulum (SR); the latter cumulatively depends on the train length. A high pacing rate intensifies the cumulative transient shift in the SR Ca²⁺ loading, augmenting the post-train inotropic response and prolonging its recovery to the pre-train level. The pacing-dependent mechano-induced inotropic effects remain to be elucidated in the myocardium with impaired Ca handling.