Experimental and modelling evidence of shortening heat in cardiac muscle

In silico study of the mechanisms of hypoxia and contractile dysfunction during ischemia and reperfusion of hiPSC cardiomyocytes

Interconnected mechanisms of ischemia and reperfusion (IR) has increased the interest in IR in vitro experiments using human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs). We developed a whole-cell computational model of hiPSC-CMs including the electromechanics, a metabolite-sensitive sarcoplasmic reticulum Ca2+-ATPase (SERCA) and an oxygen dynamics formulation to investigate IR mechanisms. Moreover, we simulated the effect and action mechanism of levosimendan, which recently showed promising anti-arrhythmic effects in hiPSC-CMs in hypoxia. The model was validated using hiPSC-CM and in vitro animal data. The role of SERCA in causing relaxation dysfunction in IR was anticipated to be comparable to its function in sepsis-induced heart failure. Drug simulations showed that levosimendan counteracts the relaxation dysfunction by utilizing a particular Ca2+-sensitizing mechanism involving Ca2+-bound troponin C and Ca2+ flux to the myofilament, rather than inhibiting SERCA phosphorylation. The model demonstrates extensive characterization and promise for drug development, making it suitable for evaluating IR therapy strategies based on the changing levels of cardiac metabolites, oxygen and molecular pathways.

Resolving an inconsistency in the estimation of the energy for excitation of cardiac muscle contraction

In the excitation of muscle contraction, calcium ions interact with transmembrane transporters. This process is accompanied by energy consumption and heat liberation. To quantify this activation energy or heat in the heart or cardiac muscle, two non-pharmacological approaches can be used. In one approach using the "pressure-volume area" concept, the same estimate of activation energy is obtained regardless of the mode of contraction (either isovolumic/isometric or ejecting/shortening). In the other approach, an accurate estimate of activation energy is obtained only when the muscle contracts isometrically. If the contraction involves muscle shortening, then an additional component of heat associated with shortening is liberated, over and above that of activation. The present study thus examines the reconcilability of the two approaches by performing experiments on isolated muscles measuring contractile force and heat output. A framework was devised from the experimental data to allow us to replicate several mechanoenergetics results gleaned from the literature. From these replications, we conclude that the choice of initial muscle length (or ventricular volume) underlies the divergence of the two approaches in the estimation of activation energy when the mode of contraction involves shortening (ejection). At low initial muscle lengths, the heat of shortening is relatively small, which can lead to the misconception that activation energy is contraction mode independent. In fact, because cardiac muscle liberates heat of shortening when allowed to shorten, estimation of activation heat must be performed only under isometric (isovolumic) contractions. We thus recommend caution when estimating activation energy using the "pressure-volume area" concept.

Altered contractility in mutation-specific hypertrophic cardiomyopathy: A mechano-energetic in silico study with pharmacological insights

Introduction: Mavacamten (MAVA), Blebbistatin (BLEB), and Omecamtiv mecarbil (OM) are promising drugs directly targeting sarcomere dynamics, with demonstrated efficacy against hypertrophic cardiomyopathy (HCM) in (pre)clinical trials. However, the molecular mechanism affecting cardiac contractility regulation, and the diseased cell mechano-energetics are not fully understood yet. Methods: We present a new metabolite-sensitive computational model of human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) electromechanics to investigate the pathology of R403Q HCM mutation and the effect of MAVA, BLEB, and OM on the cell mechano-energetics. Results: We offer a mechano-energetic HCM calibration of the model, capturing the prolonged contractile relaxation due to R403Q mutation (∼33%), without assuming any further modifications such as an additional Ca²⁺ flux to the thin filaments. The HCM model variant correctly predicts the negligible alteration in ATPase activity in R403Q HCM condition compared to normal hiPSC-CMs. The simulated inotropic effects of MAVA, OM, and BLEB, along with the ATPase activities in the control and HCM model variant agree with in vitro results from different labs. The proposed model recapitulates the tension-Ca²⁺ relationship and action potential duration change due to 1 µM OM and 5 µM BLEB, consistently with in vitro data. Finally, our model replicates the experimental dose-dependent effect of OM and BLEB on the normalized isometric tension. Conclusion: This work is a step toward deep-phenotyping the mutation-specific HCM pathophysiology, manifesting as altered interfilament kinetics. Accordingly, the modeling efforts lend original insights into the MAVA, BLEB, and OM contributions to a new interfilament balance resulting in a cardioprotective effect.

A mathematical model of hiPSC cardiomyocytes electromechanics

Human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) are becoming instrumental in cardiac research, human-based cell level cardiotoxicity tests, and developing patient-specific care. As one of the principal functional readouts is contractility, we propose a novel electromechanical hiPSC-CM computational model named the hiPSC-CM-CE. This model comprises a reparametrized version of contractile element (CE) by Rice et al., 2008, with a new passive force formulation, integrated into a hiPSC-CM electrophysiology formalism by Paci et al. in 2020. Our simulated results were validated against in vitro data reported for hiPSC-CMs at matching conditions from different labs. Specifically, key action potential (AP) and calcium transient (CaT) biomarkers simulated by the hiPSC-CM-CE model were within the experimental ranges. On the mechanical side, simulated cell shortening, contraction-relaxation kinetic indices (RT50 and RT25 ), and the amplitude of tension fell within the experimental intervals. Markedly, as an inter-scale analysis, correct classification of the inotropic effects due to non-cardiomyocytes in hiPSC-CM tissues was predicted on account of the passive force expression introduced to the CE. Finally, the physiological inotropic effects caused by Verapamil and Bay-K 8644 and the aftercontractions due to the early afterdepolarizations (EADs) were simulated and validated against experimental data. In the future, the presented model can be readily expanded to take in pharmacological trials and genetic mutations, such as those involved in hypertrophic cardiomyopathy, and study arrhythmia trigger mechanisms.

Cardiac mechanical efficiency is preserved in primary cardiac hypertrophy despite impaired mechanical function

Increased heart size is a major risk factor for heart failure and premature mortality. Although abnormal heart growth subsequent to hypertension often accompanies disturbances in mechano-energetics and cardiac efficiency, it remains uncertain whether hypertrophy is their primary driver. In this study, we aimed to investigate the direct association between cardiac hypertrophy and cardiac mechano-energetics using isolated left-ventricular trabeculae from a rat model of primary cardiac hypertrophy and its control. We evaluated energy expenditure (heat output) and mechanical performance (force length work production) simultaneously at a range of preloads and afterloads in a microcalorimeter, we determined energy expenditure related to cross-bridge cycling and Ca²⁺ cycling (activation heat), and we quantified energy efficiency. Rats with cardiac hypertrophy exhibited increased cardiomyocyte length and width. Their trabeculae showed mechanical impairment, evidenced by lower force production, extent and kinetics of shortening, and work output. Lower force was associated with lower energy expenditure related to Ca²⁺ cycling and to cross-bridge cycling. However, despite these changes, both mechanical and cross-bridge energy efficiency were unchanged. Our results show that cardiac hypertrophy is associated with impaired contractile performance and with preservation of energy efficiency. These findings provide direction for future investigations targeting metabolic and Ca²⁺ disturbances underlying cardiac mechanical and energetic impairment in primary cardiac hypertrophy.

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.

Insights From Computational Modeling Into the Contribution of Mechano-Calcium Feedback on the Cardiac End-Systolic Force-Length Relationship

In experimental studies on cardiac tissue, the end-systolic force-length relation (ESFLR) has been shown to depend on the mode of contraction: isometric or isotonic. The isometric ESFLR is derived from isometric contractions spanning a range of muscle lengths while the isotonic ESFLR is derived from shortening contractions across a range of afterloads. The ESFLR of isotonic contractions consistently lies below its isometric counterpart. Despite the passing of over a hundred years since the first insight by Otto Frank, the mechanism(s) underlying this protocol-dependent difference in the ESFLR remain incompletely explained. Here, we investigate the role of mechano-calcium feedback in accounting for the difference between these two ESFLRs. Previous studies have compared the dynamics of isotonic contractions to those of a single isometric contraction at a length that produces maximum force, without considering isometric contractions at shorter muscle lengths. We used a mathematical model of cardiac excitation-contraction to simulate isometric and force-length work-loop contractions (the latter being the 1D equivalent of the whole-heart pressure-volume loop), and compared Ca²⁺ transients produced under equivalent force conditions. We found that the duration of the simulated Ca²⁺ transient increases with decreasing sarcomere length for isometric contractions, and increases with decreasing afterload for work-loop contractions. At any given force, the Ca²⁺ transient for an isometric contraction is wider than that during a work-loop contraction. By driving simulated work-loops with wider Ca²⁺ transients generated from isometric contractions, we show that the duration of muscle shortening was prolonged, thereby shifting the work-loop ESFLR toward the isometric ESFLR. These observations are explained by an increase in the rate of binding of Ca²⁺ to troponin-C with increasing force. However, the leftward shift of the work-loop ESFLR does not superimpose on the isometric ESFLR, leading us to conclude that while mechano-calcium feedback does indeed contribute to the difference between the two ESFLRs, it does not completely account for it.

An Equivocal Final Link - Quantitative Determination of the Thermodynamic Efficiency of ATP Hydrolysis - Sullies the Chain of Electric, Ionic, Mechanical and Metabolic Steps Underlying Cardiac Contraction

Each beat of the heart completes the final step in a sequence of events commencing with electrical excitation-triggered release of Ca2+ from the sarcoplasmic reticulum which, in turn, triggers ATP-hydrolysis-dependent mechanical contraction. Given that Thermodynamics is inherently detail-independent, the heart can be thus be viewed as a mechanical pump - the generator of pressure that drives blood through the systemic and pulmonary circulations. The beat-to-beat pressure-volume work (W) of the heart is relatively straightforward to measure experimentally. Given an ability to measure, simultaneously, the accompanying heat production or oxygen consumption, it is trivial to calculate the mechanical efficiency: ε = W/ΔH where ΔH is the change of enthalpy: (W + Q), Q representing the accompanying production of heat. But it is much less straightforward to measure the thermodynamic efficiency: η = W/ΔG ATP , where ΔG ATP signifies the Gibbs Free Energy of ATP hydrolysis. The difficulty arises because of uncertain quantification of the substrate-dependent yield of ATP - conveniently expressed as the P/O2 ratio. P/O2 ratios, originally ("classically") inferred from thermal studies, have been considerably reduced over the past several decades by re-analysis of the stoichiometric coefficients separating sequential steps in the electron transport system - in particular, dropping the requirement that the coefficients have integer values. Since the early classical values are incompatible with the more recent estimates, we aim to probe this discrepancy with a view to its reconciliation. Our probe consists of a simple, thermodynamically constrained, algebraic model of cardiac mechano-energetics. Our analysis fails to reconcile recent and classical estimates of PO2 ratios; hence, we are left with a conundrum.

Energetics Equivalent of the Cardiac Force-Length End-Systolic Zone: Implications for Contractility and Economy of Contraction

We have recently demonstrated the existence of a region on the cardiac mechanics stress-length plane, which we have designated "The cardiac end-systolic zone." The zone is defined as the area on the pressure-volume (or stress-length) plane within which all stress-length contraction profiles reach their end-systolic points. It is enclosed by three boundaries: the isometric end-systolic relation, the work-loop (shortening) end-systolic relation, and the zero-active stress isotonic end-systolic relation. The existence of this zone reflects the contraction-mode dependence of the cardiac end-systolic force-length relations, and has been confirmed in a range of cardiac preparations at the whole-heart, tissue and myocyte levels. This finding has led us to speculate that a comparable zone prevails for cardiac metabolism. Specifically, we hypothesize the existence of an equivalent zone on the energetics plane (heat vs. stress), and that it can be attributed to the recently-revealed heat of shortening in cardiac muscle. To test these hypotheses, we subjected trabeculae to both isometric contractions and work-loop contractions over wide ranges of preloads and afterloads. We found that the heat-stress relations for work-loop contractions were distinct from those of isometric contractions, mirroring the contraction mode-dependence of the stress-length relation. The zone bounded by these contraction-mode dependent heat-stress relations reflects the heat of shortening. Isoproterenol-induced enhancement of contractility led to proportional increases in the zones on both the mechanics and energetics planes, thereby supporting our hypothesis.

Pulmonary arterial hypertension reduces energy efficiency of right, but not left, rat ventricular trabeculae

KEY POINTS

Pulmonary arterial hypertension (PAH) triggers right ventricle (RV) hypertrophy and left ventricle (LV) atrophy, which progressively leads to heart failure. We designed experiments under conditions mimicking those encountered by the heart in vivo that allowed us to investigate whether consequent structural and functional remodelling of the ventricles affects their respective energy efficiencies. We found that peak work output was lower in RV trabeculae from PAH rats due to reduced extent and velocity of shortening. However, their suprabasal enthalpy was unaffected due to increased activation heat, resulting in reduced suprabasal efficiency. There was no effect of PAH on LV suprabasal efficiency. We conclude that the mechanism underlying the reduced energy efficiency of hypertrophied RV tissues is attributable to the increased energy cost of Ca²⁺ cycling, whereas atrophied LV tissues still maintain normal mechano-energetic performance.

ABSTRACT

Pulmonary arterial hypertension (PAH) greatly increases the afterload on the right ventricle (RV), triggering RV hypertrophy, which progressively leads to RV failure. In contrast, the disease reduces the passive filling pressure of the left ventricle (LV), resulting in LV atrophy. We investigated whether these distinct structural and functional consequences to the ventricles affect their respective energy efficiencies. We studied trabeculae isolated from both ventricles of Wistar rats with monocrotaline-induced PAH and their respective Control groups. Trabeculae were mounted in a calorimeter at 37°C. While contracting at 5 Hz, they were subjected to stress-length work-loops over a wide range of afterloads. They were subsequently required to undergo a series of isometric contractions at various muscle lengths. In both protocols, stress production, length change and suprabasal heat output were simultaneously measured. We found that RV trabeculae from PAH rats generated higher activation heat, but developed normal active stress. Their peak external work output was lower due to reduced extent and velocity of shortening. Despite lower peak work output, suprabasal enthalpy was unaffected, thereby rendering suprabasal efficiency lower. Crossbridge efficiency, however, was unaffected. In contrast, LV trabeculae from PAH rats maintained normal mechano-energetic performance. Pulmonary arterial hypertension reduces the suprabasal energy efficiency of hypertrophied right ventricular tissues as a consequence of the increased energy cost of Ca²⁺ cycling.

Left-Ventricular Energetics in Pulmonary Arterial Hypertension-Induced Right-Ventricular Hypertrophic Failure

Pulmonary arterial hypertension (PAH) alters the geometries of both ventricles of the heart. While the right ventricle (RV) hypertrophies, the left ventricle (LV) atrophies. Multiple lines of clinical and experimental evidence lead us to hypothesize that the impaired stroke volume and systolic pressure of the LV are a direct consequence of the effect of pressure overload in the RV, and that atrophy in the LV plays only a minor role. In this study, we tested this hypothesis by examining the mechanoenergetic response of the atrophied LV to RV hypertrophy in rats treated with monocrotaline. Experiments were performed across multiple-scales: the whole-heart in vivo and ex vivo, and its trabeculae in vitro. Under the in vivo state where the RV was pressure-overloaded, we measured reduced systemic blood pressure and LV ventricular pressure. In contrast, under both ex vivo and in vitro conditions, where the effect of RV pressure overload was circumvented, we found that LV was capable of developing normal systolic pressure and stress. Nevertheless, LV atrophy played a minor role in that LV stroke volume remained lower, thereby contributing to lower LV mechanical work output. Concomitantly lower oxygen consumption and change of enthalpy were observed, and hence LV energy efficiency was unchanged. Our internally consistent findings between working-heart and trabecula experiments explain the rapid improvement of LV systolic function observed in patients with chronic pulmonary hypertension following surgical relief of RV pressure overload.