Myofilament properties comprise the rate-limiting step for cardiac relaxation at body temperature in the rat

Regulating Striated Muscle Contraction: Through Thick and Thin

Force generation in striated muscle is primarily controlled by structural changes in the actin-containing thin filaments triggered by an increase in intracellular calcium concentration. However, recent studies have elucidated a new class of regulatory mechanisms, based on the myosin-containing thick filament, that control the strength and speed of contraction by modulating the availability of myosin motors for the interaction with actin. This review summarizes the mechanisms of thin and thick filament activation that regulate the contractility of skeletal and cardiac muscle. A novel dual-filament paradigm of muscle regulation is emerging, in which the dynamics of force generation depends on the coordinated activation of thin and thick filaments. We highlight the interfilament signaling pathways based on titin and myosin-binding protein-C that couple thin and thick filament regulatory mechanisms. This dual-filament regulation mediates the length-dependent activation of cardiac muscle that underlies the control of the cardiac output in each heartbeat.

Machine learning-based classification of cardiac relaxation impairment using sarcomere length and intracellular calcium transients

Impaired relaxation of cardiomyocytes leads to diastolic dysfunction in the left ventricle. Relaxation velocity is regulated in part by intracellular calcium (Ca2+) cycling, and slower outflux of Ca2+ during diastole translates to reduced relaxation velocity of sarcomeres. Sarcomere length transient and intracellular calcium kinetics are integral parts of characterizing the relaxation behavior of the myocardium. However, a classifier tool that can separate normal cells from cells with impaired relaxation using sarcomere length transient and/or calcium kinetics remains to be developed. In this work, we employed nine different classifiers to classify normal and impaired cells, using ex-vivo measurements of sarcomere kinematics and intracellular calcium kinetics data. The cells were isolated from wild-type mice (referred to as normal) and transgenic mice expressing impaired left ventricular relaxation (referred to as impaired). We utilized sarcomere length transient data with a total of n = 126 cells (n = 60 normal cells and n = 66 impaired cells) and intracellular calcium cycling measurements with a total of n = 116 cells (n = 57 normal cells and n = 59 impaired cells) from normal and impaired cardiomyocytes as inputs to machine learning (ML) models for classification. We trained all ML classifiers with cross-validation method separately using both sets of input features, and compared their performance metrics. The performance of classifiers on test data showed that our soft voting classifier outperformed all other individual classifiers on both sets of input features, with 0.94 and 0.95 area under the receiver operating characteristic curves for sarcomere length transient and calcium transient, respectively, while multilayer perceptron achieved comparable scores of 0.93 and 0.95, respectively. However, the performance of decision tree, and extreme gradient boosting was found to be dependent on the set of input features used for training. Our findings highlight the importance of selecting appropriate input features and classifiers for the accurate classification of normal and impaired cells. Layer-wise relevance propagation (LRP) analysis demonstrated that the time to 50% contraction of the sarcomere had the highest relevance score for sarcomere length transient, whereas time to 50% decay of calcium had the highest relevance score for calcium transient input features. Despite the limited dataset, our study demonstrated satisfactory accuracy, suggesting that the algorithm can be used to classify relaxation behavior in cardiomyocytes when the potential relaxation impairment of the cells is unknown.

Suppression of lusitropy as a disease mechanism in cardiomyopathies

In cardiac muscle the action of adrenaline on β1 receptors of heart muscle cells is essential to adjust cardiac output to the body's needs. Adrenergic activation leads to enhanced contractility (inotropy), faster heart rate (chronotropy) and faster relaxation (lusitropy), mainly through activation of protein kinase A (PKA). Efficient enhancement of heart output under stress requires all of these responses to work together. Lusitropy is essential for shortening the heartbeat when heart rate increases. It therefore follows that, if the lusitropic response is not present, heart function under stress will be compromised. Current literature suggests that lusitropy is primarily achieved due to PKA phosphorylation of troponin I (TnI) and phospholamban (PLB). It has been well documented that PKA-induced phosphorylation of TnI releases Ca2+ from troponin C faster and increases the rate of cardiac muscle relaxation, while phosphorylation of PLB increases SERCA activity, speeding up Ca2+ removal from the cytoplasm. In this review we consider the current scientific evidences for the connection between suppression of lusitropy and cardiac dysfunction in the context of mutations in phospholamban and thin filament proteins that are associated with cardiomyopathies. We will discuss what advances have been made into understanding the physiological mechanism of lusitropy due to TnI and PLB phosphorylation and its suppression by mutations and we will evaluate the evidence whether lack of lusitropy is sufficient to cause cardiomyopathy, and under what circumstances, and consider the range of pathologies associated with loss of lusitropy. Finally, we will discuss whether suppressed lusitropy due to mutations in thin filament proteins can be therapeutically restored.

A novel modular modeling approach for understanding different electromechanics between left and right heart in rat

While ion channels and transporters involved in excitation-contraction coupling have been linked and constructed as comprehensive computational models, validation of whether each individual component of a model can be reused has not been previously attempted. Here we address this issue while using a novel modular modeling approach to investigate the underlying mechanism for the differences between left ventricle (LV) and right ventricle (RV). Our model was developed from modules constructed using the module assembly principles of the CellML model markup language. The components of three existing separate models of cardiac function were disassembled as to create smaller modules, validated individually, and then the component parts were combined into a new integrative model of a rat ventricular myocyte. The model was implemented in OpenCOR using the CellML standard in order to ensure reproducibility. Simulated action potential (AP), Ca²⁺ transient, and tension were in close agreement with our experimental measurements: LV AP showed a prolonged duration and a more prominent plateau compared with RV AP; Ca²⁺ transient showed prolonged duration and slow decay in LV compared to RV; the peak value and relaxation of tension were larger and slower, respectively, in LV compared to RV. Our novel approach of module-based mathematical modeling has established that the ionic mechanisms underlying the APs and Ca²⁺ handling play a role in the variation in force production between ventricles. This simulation process also provides a useful way to reuse and elaborate upon existing models in order to develop a new model.

Cooling intact and demembranated trabeculae from rat heart releases myosin motors from their inhibited conformation

Myosin filament–based regulation supplements actin filament–based regulation to control the strength and speed of contraction in heart muscle. In diastole, myosin motors form a folded helical array that inhibits actin interaction; during contraction, they are released from that array. A similar structural transition has been observed in mammalian skeletal muscle, in which cooling below physiological temperature has been shown to reproduce some of the structural features of the activation of myosin filaments during active contraction. Here, we used small-angle x-ray diffraction to characterize the structural changes in the myosin filaments associated with cooling of resting and relaxed trabeculae from the right ventricle of rat hearts from 39°C to 7°C. In intact quiescent trabeculae, cooling disrupted the folded helical conformation of the myosin motors and induced extension of the filament backbone, as observed in the transition from diastole to peak systolic force at 27°C. Demembranation of trabeculae in relaxing conditions induced expansion of the filament lattice, but the structure of the myosin filaments was mostly preserved at 39°C. Cooling of relaxed demembranated trabeculae induced changes in motor conformation and filament structure similar to those observed in intact quiescent trabeculae. Osmotic compression of the filament lattice to restore its spacing to that of intact trabeculae at 39°C stabilized the helical folded state against disruption by cooling. The myosin filament structure and motor conformation of intact trabeculae at 39°C were largely preserved in demembranated trabeculae at 27°C or above in the presence of Dextran, allowing the physiological mechanisms of myosin filament–based regulation to be studied in those conditions.

Effect of Myosin Isoforms on Cardiac Muscle Twitch of Mice, Rats and Humans

To understand how pathology-induced changes in contractile protein isoforms modulate cardiac muscle function, it is necessary to quantify the temporal-mechanical properties of contractions that occur under various conditions. Pathological responses are much easier to study in animal model systems than in humans, but extrapolation between species presents numerous challenges. Employing computational approaches can help elucidate relationships that are difficult to test experimentally by translating the observations from rats and mice, as model organisms, to the human heart. Here, we use the spatially explicit MUSICO platform to model twitch contractions from rodent and human trabeculae collected in a single laboratory. This approach allowed us to identify the variations in kinetic characteristics of α- and β-myosin isoforms across species and to quantify their effect on cardiac muscle contractile responses. The simulations showed how the twitch transient varied with the ratio of the two myosin isoforms. Particularly, the rate of tension rise was proportional to the fraction of α-myosin present, while the β-isoform dominated the rate of relaxation unless α-myosin was >50%. Moreover, both the myosin isoform and the Ca2+ transient contributed to the twitch tension transient, allowing two levels of regulation of twitch contraction.

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.

Contraction-relaxation coupling is unaltered by exercise training and infarction in isolated canine myocardium

The two main phases of the mammalian cardiac cycle are contraction and relaxation; however, whether there is a connection between them in humans is not well understood. Routine exercise has been shown to improve cardiac function, morphology, and molecular signatures. Likewise, the acute and chronic changes that occur in the heart in response to injury, disease, and stress are well characterized, albeit not fully understood. In this study, we investigated how exercise and myocardial injury affect contraction–relaxation coupling. We retrospectively analyzed the correlation between the maximal speed of contraction and the maximal speed of relaxation of canine myocardium after receiving surgically induced myocardial infarction, followed by either sedentary recovery or exercise training for 10–12 wk. We used isolated right ventricular trabeculae, which were electrically paced at different lengths, frequencies, and with increasing β-adrenoceptor stimulation. In all conditions, contraction and relaxation were linearly correlated, irrespective of injury or training history. Based on these results and the available literature, we posit that contraction–relaxation coupling is a fundamental myocardial property that resides in the structural arrangement of proteins at the level of the sarcomere and that this may be regulated by the actions of cardiac myosin binding protein C (cMyBP-C) on actin and myosin.

Multiscale modeling of twitch contractions in cardiac trabeculae

Understanding the dynamics of a cardiac muscle twitch contraction is complex because it requires a detailed understanding of the kinetic processes of the Ca2+ transient, thin-filament activation, and the myosin-actin cross-bridge chemomechanical cycle. Each of these steps has been well defined individually, but understanding how all three of the processes operate in combination is a far more complex problem. Computational modeling has the potential to provide detailed insight into each of these processes, how the dynamics of each process affect the complexity of contractile behavior, and how perturbations such as mutations in sarcomere proteins affect the complex interactions of all of these processes. The mechanisms involved in relaxation of tension during a cardiac twitch have been particularly difficult to discern due to nonhomogeneous sarcomere lengthening during relaxation. Here we use the multiscale MUSICO platform to model rat trabecular twitches. Validation of computational models is dependent on being able to simulate different experimental datasets, but there has been a paucity of data that can provide all of the required parameters in a single experiment, such as simultaneous measurements of force, intracellular Ca2+ transients, and sarcomere length dynamics. In this study, we used data from different studies collected under similar experimental conditions to provide information for all the required parameters. Our simulations established that twitches either in an isometric sarcomere or in fixed-length, multiple-sarcomere trabeculae replicate the experimental observations if models incorporate a length-tension relationship for the nonlinear series elasticity of muscle preparations and a scheme for thick-filament regulation. The thick-filament regulation assumes an off state in which myosin heads are parked onto the thick-filament backbone and are unable to interact with actin, a state analogous to the super-relaxed state. Including these two mechanisms provided simulations that accurately predict twitch contractions over a range of different conditions.

The effect of variable troponin C mutation thin filament incorporation on cardiac muscle twitch contractions

One of the complexities of understanding the pathology of familial forms of cardiac diseases is the level of mutation incorporation in sarcomeres. Computational models of the sarcomere that are spatially explicit offer an approach to study aspects of mutational incorporation into myofilaments that are more challenging to get at experimentally. We studied two well characterized mutations of cardiac TnC, L48Q and I61Q, that decrease or increase the release rate of Ca²⁺ from cTnC, k_-Ca, resulting in HCM and DCM respectively [1]. Expression of these mutations in transgenic mice was used to provide experimental data for incorporation of 30 and 50% (respectively) into sarcomeres. Here we demonstrate that fixed length twitch contractions of trabeculae from mice containing mutant differ from WT; L48Q trabeculae have slower relaxation while I61Q trabeculae have markedly reduced peak tension. Using our multiscale modelling approach [2] we were able to describe the tension transients of WT mouse myocardium. Tension transients for the mutant cTnCs were simulated with changes in k_-Ca, measured experimentally for each cTnC mutant in whole troponin complex, a change in the affinity of cTnC for cTnI, and a reduction in the number of detached crossbridges available for binding. A major advantage of the multiscale explicit 3-D model is that it predicts the effects of variable mutation incorporation, and the effects of variations in mutation distribution within thin filaments in sarcomeres. Such effects are currently impossible to explore experimentally. We explored random and clustered distributions of mutant cTnCs in thin filaments, as well as distributions of individual thin filaments with only WT or mutant cTnCs present. The effects of variable amounts of incorporation and non-random distribution of mutant cTnCs are more marked for I61Q than L48Q cTnC. We conclude that this approach can be effective for study on mutations in multiple proteins of the sarcomere. SUMMARY: A challenge in experimental studies of diseases is accounting for the effect of variable mutation incorporation into myofilaments. Here we use a spatially explicit computational approach, informed by experimental data from transgenic mice expressing one of two mutations in cardiac Troponin C that increase or decrease calcium sensitivity. We demonstrate that the model can accurately describe twitch contractions for the data and go on to explore the effect of variable mutant incorporation and localization on simulated cardiac muscle twitches.

Biophysically detailed mathematical models of multiscale cardiac active mechanics

We propose four novel mathematical models, describing the microscopic mechanisms of force generation in the cardiac muscle tissue, which are suitable for multiscale numerical simulations of cardiac electromechanics. Such models are based on a biophysically accurate representation of the regulatory and contractile proteins in the sarcomeres. Our models, unlike most of the sarcomere dynamics models that are available in the literature and that feature a comparable richness of detail, do not require the time-consuming Monte Carlo method for their numerical approximation. Conversely, the models that we propose only require the solution of a system of PDEs and/or ODEs (the most reduced of the four only involving 20 ODEs), thus entailing a significant computational efficiency. By focusing on the two models that feature the best trade-off between detail of description and identifiability of parameters, we propose a pipeline to calibrate such parameters starting from experimental measurements available in literature. Thanks to this pipeline, we calibrate these models for room-temperature rat and for body-temperature human cells. We show, by means of numerical simulations, that the proposed models correctly predict the main features of force generation, including the steady-state force-calcium and force-length relationships, the length-dependent prolongation of twitches and increase of peak force, the force-velocity relationship. Moreover, they correctly reproduce the Frank-Starling effect, when employed in multiscale 3D numerical simulation of cardiac electromechanics.

Stretching single titin molecules from failing human hearts reveals titin's role in blunting cardiac kinetic reserve

AIMS

Heart failure (HF) patients commonly experience symptoms primarily during elevated heart rates, as a result of physical activities or stress. A main determinant of diastolic passive tension, the elastic sarcomeric protein titin, has been shown to be associated with HF, with unresolved involvement regarding its role at different heart rates. To determine whether titin is playing a role in the heart rate (frequency-) dependent acceleration of relaxation (FDAR). W, we studied the FDAR responses in live human left ventricular cardiomyocytes and the corresponding titin-based passive tension (TPT) from failing and non-failing human hearts.

METHODS AND RESULTS

Using atomic force, we developed a novel single-molecule force spectroscopy approach to detect TPT based on the frequency-modulated cardiac cycle. Mean TPT reduced upon an increased heart rate in non-failing human hearts, while this reduction was significantly blunted in failing human hearts. These mechanical changes in the titin distal Ig domain significantly correlated with the frequency-dependent relaxation kinetics of human cardiomyocytes obtained from the corresponding hearts. Furthermore, the data suggested that the higher the TPT, the faster the cardiomyocytes relaxed, but the lower the potential of myocytes to speed up relaxation at a higher heart rate. Such poorer FDAR response was also associated with a lesser reduction or a bigger increase in TPT upon elevated heart rate.

CONCLUSIONS

Our study established a novel approach in detecting dynamic heart rate relevant tension changes physiologically on native titin domains. Using this approach, the data suggested that the regulation of kinetic reserve in cardiac relaxation and its pathological changes were associated with the intensity and dynamic changes of passive tension by titin.

Thermal Activation of Thin Filaments in Striated Muscle

In skeletal and cardiac muscles, contraction is triggered by an increase in the intracellular Ca2+ concentration. During Ca2+ transients, Ca2+-binding to troponin C shifts the "on-off" equilibrium of the thin filament state toward the "on" sate, promoting actomyosin interaction. Likewise, recent studies have revealed that the thin filament state is under the influence of temperature; viz., an increase in temperature increases active force production. In this short review, we discuss the effects of temperature on the contractile performance of mammalian striated muscle at/around body temperature, focusing especially on the temperature-dependent shift of the "on-off" equilibrium of the thin filament state.

Length-Dependent Activation of Contractility and Ca-Transient Kinetics in Auxotonically Contracting Isolated Rat Ventricular Cardiomyocytes

Length-dependent activation (LDA) of contraction is an important mechanism of proper myocardial function that is often blunted in diseases accompanied by deficient contractility and impaired calcium homeostasis. We evaluated how the extent of LDA is related to the decreased force in healthy rat myocardium under negative inotropic conditions that affect the calcium cycle. The length-dependent effects on auxotonic twitch and Ca-transient were compared in isolated rat ventricular cardiomyocytes at room temperature (“25C”) and near-physiological temperature (“35C”) in normal Tyrode and at 25°C with thapsigargin-depleted sarcoplasmic reticulum (“25C + Thap”). At the slack length, a similar negative inotropy in “35C” and “25C + Thap” was accompanied by totally different changes in Ca-transient amplitude, time-to-peak, and time-to-decline from peak to 50% amplitude. End-systolic/end-diastolic tension-sarcomere length relationships were obtained for each individual cell, and the ratio of their slopes, the dimensionless Frank-Starling Gain index, was 2.32 ± 0.16, 1.78 ± 0.09, and 1.37 ± 0.06 in “25C,” “35C” and “25C + Thap,” respectively (mean ± S.E.M.). Ca-transient diastolic level and amplitude did not differ between “25C” and “35C” at any SL, but in “35C” it developed and declined significantly faster. In contrast, thapsigargin-induced depletion of SERCA2a significantly attenuated and retarded Ca-transient. The relative amount of Ca²⁺ utilized by troponin C, evaluated by the integral magnitude of a short-lived component of Ca-transient decline (“bump”), increased by ~25% per each 0.05 μm increase in SL in all groups. The kinetics of the Ca-TnC dissociation, evaluated by the bump time-to-peak, was significantly faster in “35C” and slower in “25C + Thap” vs. “25C” (respectively, 63.7 ± 5.3 and 253.6 ± 8.3% of the value in “25C,” mean ± S.E.M.). In conclusion, a similar inotropic effect can be observed in rat ventricular myocardium under totally different kinetics of free cytosolic calcium. The extent of LDA is not determined by actual peak systolic tension but is regulated by the level of peak systolic calcium and the kinetics of Ca-transient decline which, in turn, are governed by Ca-TnC dissociation and Ca²⁺ reuptake by the sarcoplasmic reticulum. Altogether, these findings constitute new evidence about the role of the length-dependent modulation of Ca²⁺ homeostasis in the mechanisms of calcium regulation of contraction and mechano-calcium feedback in the myocardium.

Investigating the inotropic effect of pyruvic acid on the isolated rat heart and its underlying mechanism

Pyruvic acid is important organic chemical intermediates that plays a role in cardiomyocyte pathophysiology and therapy. This study sought to explore the inotropic effects of pyruvic acid on the function of the isolated rat hearts and investigate its underlying mechanism. Pyruvic acid produced a greater negative inotropic effect compared to HCl and sodium pyruvate in a concentration-dependent pattern in the hearts. The role of low dose of pyruvic acid on heart function was regulated by pyruvic acid molecules and high dose pyruvic acid may be influenced by pyruvic acid molecules and pH. K_v channels may be involved in the pyruvic acid-induced negative inotropic effect. Finally, pyruvic acid markedly increased the level of LDH and CK and reduced the level of Ca²⁺Mg²⁺-ATPase and Na⁺K⁺-ATPase. These results suggest that pyruvic acid may modulate cardiac function at physiological or low doses but can cause damage to cardiomyocytes at high doses.

Omecamtiv Mecarbil Slows Myosin Kinetics in Skinned Rat Myocardium at Physiological Temperature

Heart failure is a life-threatening condition that occurs when the heart muscle becomes weakened and cannot adequately circulate blood and nutrients around the body. Omecamtiv mecarbil (OM) is a compound that has been developed to treat systolic heart failure via targeting the cardiac myosin heavy chain to increase myocardial contractility. Biophysical and biochemical studies have found that OM increases calcium (Ca²⁺) sensitivity of contraction by prolonging the myosin working stroke and increasing the actin-myosin cross-bridge duty ratio. Most in vitro studies probing the effects of OM on cross-bridge kinetics and muscle force production have been conducted at subphysiological temperature, even though temperature plays a critical role in enzyme activity and cross-bridge function. Herein, we used skinned, ventricular papillary muscle strips from rats to investigate the effects of [OM] on Ca²⁺-activated force production, cross-bridge kinetics, and myocardial viscoelasticity at physiological temperature (37°C). We find that OM only increases myocardial contractility at submaximal Ca²⁺ activation levels and not maximal Ca²⁺ activation levels. As [OM] increased, the kinetic rate constants for cross-bridge recruitment and detachment slowed for both submaximal and maximal Ca²⁺-activated conditions. These findings support a mechanism by which OM increases cardiac contractility at physiological temperature via increasing cross-bridge contributions to thin-filament activation as cross-bridge kinetics slow and the duration of cross-bridge attachment increases. Thus, force only increases at submaximal Ca²⁺ activation due to cooperative recruitment of neighboring cross-bridges, because thin-filament activation is not already saturated. In contrast, OM does not increase myocardial force production for maximal Ca²⁺-activated conditions at physiological temperature because cooperative activation of thin filaments may already be saturated.

Myocardial relaxation in human heart failure: Why sarcomere kinetics should be center-stage

Myocardial relaxation is critical for the heart to allow for adequate filling of the ventricles prior to the next contraction. In human heart failure, impairment of myocardial relaxation is a major problem, and impacts most patients suffering from end-stage failure. Furthering our understanding of myocardial relaxation is critical in developing future treatment strategies. This review highlights processes involved in myocardial relaxation, as well as governing processes that modulate myocardial relaxation, with a focus on impairment of myocardium-level relaxation in human end-stage heart failure.

Assessment of PKA and PKC inhibitors on force and kinetics of non-failing and failing human myocardium

AIMS

Heart failure (HF) is a prevalent disease that is considered the foremost reason for hospitalization in the United States. Most protein kinases (PK) are activated in heart disease and their inhibition has been shown to improve cardiac function in both animal and human studies. However, little is known about the direct impact of PKA and PKC inhibitors on human cardiac contractile function.

MATERIAL AND METHODS

We investigated the ex vivo effect of such inhibitors on force as well as on kinetics of left ventricular (LV) trabeculae dissected from non-failing and failing human hearts. In these experiments, we applied 0.5 μM of H-89 and GF109203X, which are PKA and PKC inhibitors, respectively, in comparison to their vehicle DMSO (0.05%).

KEY FINDINGS AND CONCLUSION

Statistical analyses revealed no significant effect for H-89 and GF109203X on either contractile force or kinetics parameters of both non-failing and failing muscles even though they were used at a concentration higher than the reported IC₅₀s and K_is. Therefore, several factors such as selectivity, concentration, and treatment time, which are related to these PK inhibitors according to previous studies require further exploration.

High tension in sarcomeres hinders myocardial relaxation: A computational study

Experiments have shown that the relaxation phase of cardiac sarcomeres during an isometric twitch is prolonged in muscles that reached a higher peak tension. However, the mechanism is not completely understood. We hypothesize that the binding of calcium to troponin is enhanced by the tension in the thin filament, thus contributing to the prolongation of contraction upon higher peak tension generation. To test this hypothesis, we developed a computational model of sarcomere mechanics that incorporates tension-dependence of calcium binding. The model was used to simulate isometric twitch experiments with time dependency in the form of a two-state cross-bridge cycle model and a transient intracellular calcium concentration. In the simulations, peak isometric twitch tension appeared to increase linearly by 51.1 KPa with sarcomere length from 1.9 μm to 2.2 μm. Experiments showed an increase of 47.3 KPa over the same range of sarcomere lengths. The duration of the twitch also increased with both sarcomere length and peak intracellular calcium concentration, likely to be induced by the inherently coupled increase of the peak tension in the thin filament. In the model simulations, the time to 50% relaxation (tR50) increased over the range of sarcomere lengths from 1.9 μm to 2.2 μm by 0.11s, comparable to the increased duration of 0.12s shown in experiments. Model simulated tR50 increased by 0.12s over the range of peak intracellular calcium concentrations from 0.87 μM to 1.45 μM. Our simulation results suggest that the prolongation of contraction at higher tension is a result of the tighter binding of Ca2+ to troponin in areas under higher tension, thus delaying the deactivation of the troponin.

A new myofilament contraction model with ATP consumption for ventricular cell model

A new contraction model of cardiac muscle was developed by combining previously described biochemical and biophysical models. The biochemical component of the new contraction model represents events in the presence of Ca2+-crossbridge attachment and power stroke following inorganic phosphate release, detachment evoked by the replacement of ADP by ATP, ATP hydrolysis, and recovery stroke. The biophysical component focuses on Ca2+ activation and force (F b) development assuming an equivalent crossbridge. The new model faithfully incorporates the major characteristics of the biochemical and biophysical models, such as F b activation by transient Ca2+ ([Ca2+]-F b), [Ca2+]-ATP hydrolysis relations, sarcomere length-F b, and F b recovery after jumps in length under the isometric mode and upon sarcomere shortening after a rapid release of mechanical load under the isotonic mode together with the load-velocity relationship. ATP consumption was obtained for all responses. When incorporated in a ventricular cell model, the contraction model was found to share approximately 60% of the total ATP usage in the cell model.

Human Myocardium Has a Robust α1A-Subtype Adrenergic Receptor Inotropic Response

Recent studies report that a single subtype of α1-adrenergic receptor (α1-AR), the α1A-subtype, mediates robust cardioprotective effects in multiple experimental models of heart failure, suggesting that the α1A-subtype is a potential therapeutic target for an agonist to treat heart failure. Moreover, we recently found that the α1A-subtype is present in human heart. The goal of this study was to assess the inotropic response mediated by the α1A-subtype in human myocardium, and to determine whether the response is downregulated in myocardium from failing human heart. We measured in vitro contractile responses of cardiac muscle preparations (trabeculae) isolated from the right ventricle from nonfailing and failing human hearts. Addition of the α1A-subtype agonist A61603 (100 nM) resulted in a large positive inotropic response (force increased ≈ 2-fold). This response represented ≈70% of the response mediated by the β-adrenergic receptor agonist isoproterenol (1 μM). Moreover, in myocardium from failing hearts, α1A-subtype responses remained robust, and only slightly reduced relative to nonfailing hearts. We conclude that α1A-subtype-mediated inotropy could represent a significant source of inotropic support in the human heart. Furthermore, the α1A-subtype remains functional in myocardium from failing human hearts and thus, might be a therapeutic target to support cardioprotective effects in patients with heart failure.

Synchronization of Intracellular Ca2+ Release in Multicellular Cardiac Preparations

In myocardial tissue, Ca²⁺ release from the sarcoplasmic reticulum (SR) that occurs via the ryanodine receptor (RyR2) channel complex. Ca²⁺ release through RyR2 can be either stimulated by an action potential (AP) or spontaneous. The latter is often associated with triggered afterdepolarizations, which in turn may lead to sustained arrhythmias. It is believed that some synchronization mechanism exists for afterdepolarizations and APs in neighboring myocytes, possibly a similarly timed recovery of RyR2 from refractoriness, which enables RyR2s to reach the threshold for spontaneous Ca²⁺ release simultaneously. To investigate this synchronization mechanism in absence of genetic factors that predispose arrhythmia, we examined the generation of triggered activity in multicellular cardiac preparations. In myocardial trabeculae from the rat, we demonstrated that in the presence of both isoproterenol and caffeine, neighboring myocytes within the cardiac trabeculae were able to synchronize their diastolic spontaneous SR Ca²⁺ release. Using confocal Ca²⁺ imaging, we could visualize Ca²⁺ waves in the multicellular preparation, while these waves were not always present in every myocyte within the trabeculae, we observed that, over time, the Ca²⁺ waves can synchronize in multiple myocytes. This synchronized activity was sufficiently strong that it could trigger a synchronized, propagated contraction in the whole trabecula encompassing even previously quiescent myocytes. The detection of Ca²⁺ dynamics in individual myocytes in their in situ setting at the multicellular level exposed a synchronization mechanism that could induce local triggered activity in the heart in the absence of global Ca²⁺ dysregulation.

Etiology-dependent impairment of relaxation kinetics in right ventricular end-stage failing human myocardium

BACKGROUND

In patients with end-stage heart failure, the primary etiology often originates in the left ventricle, and eventually the contractile function of the right ventricle (RV) also becomes compromised. RV tissue-level deficits in contractile force and/or kinetics need quantification to understand involvement in ischemic and non-ischemic failing human myocardium.

METHODS AND RESULTS

The human population suffering from heart failure is diverse, requiring many subjects to be studied in order to perform an adequately powered statistical analysis. From 2009-present we assessed live tissue-level contractile force and kinetics in isolated myocardial RV trabeculae from 44 non-failing and 41 failing human hearts. At 1 Hz stimulation rate (in vivo resting state) the developed active force was not different in non-failing compared to failing ischemic nor non-ischemic failing trabeculae. In sharp contrast, the kinetics of relaxation were significantly impacted by disease, with 50% relaxation time being significantly shorter in non-failing vs. non-ischemic failing, while the latter was still significantly shorter than ischemic failing. Gender did not significantly impact kinetics. Length-dependent activation was not impacted. Although baseline force was not impacted, contractile reserve was critically blunted. The force-frequency relation was positive in non-failing myocardium, but negative in both ischemic and non-ischemic myocardium, while the β-adrenergic response to isoproterenol was depressed in both pathologies.

CONCLUSIONS

Force development at resting heart rate is not impacted by cardiac pathology, but kinetics are impaired and the magnitude of the impairment depends on the underlying etiology. Focusing on restoration of myocardial kinetics will likely have greater therapeutic potential than targeting force of contraction.

Engineering Cardiac Muscle Tissue: A Maturating Field of Research

Twenty years after the initial description of a tissue engineered construct, 3-dimensional human cardiac tissues of different kinds are now generated routinely in many laboratories. Advances in stem cell biology and engineering allow for the generation of constructs that come close to recapitulating the complex structure of heart muscle and might, therefore, be amenable to industrial (eg, drug screening) and clinical (eg, cardiac repair) applications. Whether the more physiological structure of 3-dimensional constructs provides a relevant advantage over standard 2-dimensional cell culture has yet to be shown in head-to-head-comparisons. The present article gives an overview on current strategies of cardiac tissue engineering with a focus on different hydrogel methods and discusses perspectives and challenges for necessary steps toward the real-life application of cardiac tissue engineering for disease modeling, drug development, and cardiac repair.

Insights and Challenges of Multi-Scale Modeling of Sarcomere Mechanics in cTn and Tm DCM Mutants-Genotype to Cellular Phenotype

Dilated Cardiomyopathy (DCM) is a leading cause of sudden cardiac death characterized by impaired pump function and dilatation of cardiac ventricles. In this review we discuss various in silico approaches to elucidating the mechanisms of genetic mutations leading to DCM. The approaches covered in this review focus on bridging the spatial and temporal gaps that exist between molecular and cellular processes. Mutations in sarcomeric regulatory thin filament proteins such as the troponin complex (cTn) and Tropomyosin (Tm) have been associated with DCM. Despite the experimentally-observed myofilament measures of contractility in the case of these mutations, the mechanisms by which the underlying molecular changes and protein interactions scale up to organ failure by these mutations remains elusive. The review highlights multi-scale modeling approaches and their applicability to study the effects of sarcomeric gene mutations in-silico. We discuss some of the insights that can be gained from computational models of cardiac biomechanics when scaling from molecular states to cellular level.

Can inorganic phosphate explain sag during unfused tetanic contractions of skeletal muscle?

We test the hypothesis that cytosolic inorganic phosphate (P_i) can account for the contraction‐induced reductions in twitch duration which impair summation and cause force to decline (sag) during unfused tetanic contractions of fast‐twitch muscle. A five‐state model of crossbridge cycling was used to simulate twitch and unfused tetanic contractions. As P_i concentration ([P_i]) was increased from 0 to 30 mmol·L⁻¹, twitch duration decreased, with progressive reductions in sensitivity to P_i as [P_i] was increased. When unfused tetani were simulated with rising [P_i], sag was most pronounced when initial [P_i] was low, and when the magnitude of [P_i] increase was large. Fast‐twitch extensor digitorum longus (EDL) muscles (sag‐prone, typically low basal [P_i]) and slow‐twitch soleus muscles (sag‐resistant, typically high basal [P_i]) were isolated from 14 female C57BL/6 mice. Muscles were sequentially incubated in solutions containing either glucose or pyruvate to create typical and low P_i environments, respectively. Twitch duration was greater (P < 0.05) in pyruvate than glucose in both muscles. Stimuli applied at intervals approximately three times the time to peak twitch tension resulted in sag of 35.0 ± 3.7% in glucose and 50.5 ± 1.4% in pyruvate in the EDL (pyruvate > glucose; P < 0.05), and 3.9 ± 0.3% in glucose and 37.8 ± 2.7% in pyruvate in the soleus (pyruvate > glucose; P < 0.05). The influence of P_i on crossbridge cycling provides a tenable mechanism for sag. Moreover, the low basal [P_i] in fast‐twitch relative to slow‐twitch muscle has promise as an explanation for the fiber‐type dependency of sag.

Cardiac activation heat remains inversely dependent on temperature over the range 27-37°C

The relation between heat output and stress production (force per cross-sectional area) of isolated cardiac tissue is a key metric that provides insight into muscle energetic performance. The heat intercept of the relation, termed "activation heat," reflects the metabolic cost of restoring transmembrane gradients of Na(+) and K(+) following electrical excitation, and myoplasmic Ca(2+) concentration following its release from the sarcoplasmic reticulum. At subphysiological temperatures, activation heat is inversely dependent on temperature. Thus one may presume that activation heat would decrease even further at body temperature. However, this assumption is prima facie inconsistent with a study, using intact hearts, which revealed no apparent change in the combination of activation and basal metabolism between 27 and 37°C. It is thus desired to directly determine the change in activation heat between 27 and 37°C. In this study, we use our recently constructed high-thermal resolution muscle calorimeter to determine the first heat-stress relation of isolated cardiac muscle at 37°C. We compare the relation at 37°C to that at 27°C to examine whether the inverse temperature dependence of activation heat, observed under hypothermic conditions, prevails at body temperature. Our results show that activation heat was reduced (from 3.5 ± 0.3 to 2.3 ± 0.3 kJ/m(3)) at the higher temperature. This leads us to conclude that activation metabolism continues to decline as temperature is increased from hypothermia to normothermia and allows us to comment on results obtained from the intact heart by previous investigators.

Influence of metabolic dysfunction on cardiac mechanics in decompensated hypertrophy and heart failure

Alterations in energetic state of the myocardium are associated with decompensated heart failure in humans and in animal models. However, the functional consequences of the observed changes in energetic state on mechanical function are not known. The primary aim of the study was to quantify mechanical/energetic coupling in the heart and to determine if energetic dysfunction can contribute to mechanical failure. A secondary aim was to apply a quantitative systems pharmacology analysis to investigate the effects of drugs that target cross-bridge cycling kinetics in heart failure-associated energetic dysfunction. Herein, a model of metabolite- and calcium-dependent myocardial mechanics was developed from calcium concentration and tension time courses in rat cardiac muscle obtained at different lengths and stimulation frequencies. The muscle dynamics model accounting for the effect of metabolites was integrated into a model of the cardiac ventricles to simulate pressure-volume dynamics in the heart. This cardiac model was integrated into a simple model of the circulation to investigate the effects of metabolic state on whole-body function. Simulations predict that reductions in metabolite pools observed in canine models of heart failure can cause systolic dysfunction, blood volume expansion, venous congestion, and ventricular dilation. Simulations also predict that myosin-activating drugs may partially counteract the effects of energetic state on cross-bridge mechanics in heart failure while increasing myocardial oxygen consumption. Our model analysis demonstrates how metabolic changes observed in heart failure are alone sufficient to cause systolic dysfunction and whole-body heart failure symptoms.

Computational analysis of the regulation of Ca(2+) dynamics in rat ventricular myocytes

Force-frequency relationships of isolated cardiac myocytes show complex behaviors that are thought to be specific to both the species and the conditions associated with the experimental preparation. Ca(2+) signaling plays an important role in shaping the force-frequency relationship, and understanding the properties of the force-frequency relationship in vivo requires an understanding of Ca(2+) dynamics under physiologically relevant conditions. Ca(2+) signaling is itself a complicated process that is best understood on a quantitative level via biophysically based computational simulation. Although a large number of models are available in the literature, the models are often a conglomeration of components parameterized to data of incompatible species and/or experimental conditions. In addition, few models account for modulation of Ca(2+) dynamics via β-adrenergic and calmodulin-dependent protein kinase II (CaMKII) signaling pathways even though they are hypothesized to play an important regulatory role in vivo. Both protein-kinase-A and CaMKII are known to phosphorylate a variety of targets known to be involved in Ca(2+) signaling, but the effects of these pathways on the frequency- and inotrope-dependence of Ca(2+) dynamics are not currently well understood. In order to better understand Ca(2+) dynamics under physiological conditions relevant to rat, a previous computational model is adapted and re-parameterized to a self-consistent dataset obtained under physiological temperature and pacing frequency and updated to include β-adrenergic and CaMKII regulatory pathways. The necessity of specific effector mechanisms of these pathways in capturing inotrope- and frequency-dependence of the data is tested by attempting to fit the data while including and/or excluding those effector components. We find that: (1) β-adrenergic-mediated phosphorylation of the L-type calcium channel (LCC) (and not of phospholamban (PLB)) is sufficient to explain the inotrope-dependence; and (2) that CaMKII-mediated regulation of neither the LCC nor of PLB is required to explain the frequency-dependence of the data.

Quantifying inter-species differences in contractile function through biophysical modelling

Animal models and measurements are frequently used to guide and evaluate clinical interventions. In this context, knowledge of inter-species differences in physiology is crucial for understanding the limitations and relevance of animal experimental assays for informing clinical applications. Extensive effort has been put into studying the structure and function of cardiac contractile proteins and how differences in these translate into the functional properties of muscles. However, integrating this knowledge into a quantitative description, formalising and highlighting inter-species differences both in the kinetics and in the regulation of physiological mechanisms, remains challenging. In this study we propose and apply a novel approach for the quantification of inter-species differences between mouse, rat and human. Assuming conservation of the fundamental physiological mechanisms underpinning contraction, biophysically based computational models are fitted to simulate experimentally recorded phenotypes from multiple species. The phenotypic differences between species are then succinctly quantified as differences in the biophysical model parameter values. This provides the potential of quantitatively establishing the human relevance of both animal-based experimental and computational models for application in a clinical context. Our results indicate that the parameters related to the sensitivity and cooperativity of calcium binding to troponin C and the activation and relaxation rates of tropomyosin/crossbridge binding kinetics differ most significantly between mouse, rat and human, while for example the reference tension, as expected, shows only minor differences between the species. Hence, while confirming expected inter-species differences in calcium sensitivity due to large differences in the observed calcium transients, our results also indicate more unexpected differences in the cooperativity mechanism. Specifically, the decrease in the unbinding rate of calcium to troponin C with increasing active tension was much lower for mouse than for rat and human. Our results also predicted crossbridge binding to be slowest in human and fastest in mouse.

Tri-modal regulation of cardiac muscle relaxation; intracellular calcium decline, thin filament deactivation, and cross-bridge cycling kinetics

Cardiac muscle relaxation is an essential step in the cardiac cycle. Even when the contraction of the heart is normal and forceful, a relaxation phase that is too slow will limit proper filling of the ventricles. Relaxation is too often thought of as a mere passive process that follows contraction. However, many decades of advancements in our understanding of cardiac muscle relaxation have shown it is a highly complex and well-regulated process. In this review, we will discuss three distinct events that can limit the rate of cardiac muscle relaxation: the rate of intracellular calcium decline, the rate of thin-filament de-activation, and the rate of cross-bridge cycling. Each of these processes are directly impacted by a plethora of molecular events. In addition, these three processes interact with each other, further complicating our understanding of relaxation. Each of these processes is continuously modulated by the need to couple bodily oxygen demand to cardiac output by the major cardiac physiological regulators. Length-dependent activation, frequency-dependent activation, and beta-adrenergic regulation all directly and indirectly modulate calcium decline, thin-filament deactivation, and cross-bridge kinetics. We hope to convey our conclusion that cardiac muscle relaxation is a process of intricate checks and balances, and should not be thought of as a single rate-limiting step that is regulated at a single protein level. Cardiac muscle relaxation is a system level property that requires fundamental integration of three governing systems: intracellular calcium decline, thin filament deactivation, and cross-bridge cycling kinetics.

Cardiac troponin I Pro82Ser variant induces diastolic dysfunction, blunts β-adrenergic response, and impairs myofilament cooperativity

Troponin I (TnI) variant Pro82Ser (cTnIP82S) was initially considered a disease-causing mutation; however, later studies suggested the contrary. We tested the hypothesis of whether a causal link exists between cTnIP82S and cardiac structural and functional remodeling, such as during aging or chronic pressure overload. A cardiac-specific transgenic (Tg) mouse model of cTnIP82S was created to test this hypothesis. During aging, Tg cTnIP82S displayed diastolic dysfunction, characterized by longer isovolumetric relaxation time, and impaired ejection and relaxation time. In young, Tg mice in vivo pressure-volume loops and intact trabecular preparations revealed normal cardiac contractility at baseline. However, upon β-adrenergic stimulation, a blunted contractile reserve and no hastening in left ventricle relaxation were evident in vivo, whereas, in isolated muscles, Ca(2+) transient amplitude isoproterenol dose-response was blunted. In addition, when exposed to chronic pressure overload, Tg mice show exacerbated hypertrophy and decreased contractility compared with age-matched non-Tg littermates. At the molecular level, this mutation significantly impairs myofilament cooperative activation. Importantly, this occurs in the absence of alterations in TnI or myosin-binding protein C phosphorylation. The cTnIP82S variant occurs near a region of interactions with troponin T; therefore, structural changes in this region could explain its meaningful effects on myofilament cooperativity. Our data indicate that cTnIP82S mutation modifies age-dependent diastolic dysfunction and impairs overall contractility after β-adrenergic stimulation or chronic pressure overload. Thus cTnIP82S variant should be regarded as a disease-modifying factor for dysfunction and adverse remodeling with aging and chronic pressure overload.

A N-terminal truncated intracellular isoform of matrix metalloproteinase-2 impairs contractility of mouse myocardium

The full-length isoform of matrixmetalloproteinase-2 (FL-MMP-2) plays a role in turnover of the cardiac extracellular matrix. FL-MMP-2 is also present intracellularly in association with sarcomeres and, in the setting of oxidative stress, cleaves myofilament proteins with resultant impaired contractility. Recently, a novel N-terminal truncated MMP-2 isoform (NTT-MMP-2) generated during oxidative stress was identified and shown to induce severe systolic failure; however, the injury mechanisms remained unclear. In this study, cardiac-specific NTT-MMP-2 transgenic mice were used to determine the physiological effects of NTT-MMP-2 on: force development of intact myocardium; the function of cardiac myofilaments in demembranated myocardium; and on intracellular Ca²⁺ transients in isolated myocytes. We related the contractile defects arising from NTT-MMP-2 expression to the known intracellular locations of NTT-MMP-2 determined using immunohistochemistry. Comparison was made with the pathophysiology arising from cardiac-specific FL-MMP-2 transgenic mice. Consistent with previous studies, FL-MMP-2 was localized to myofilaments, while NTT-MMP-2 was concentrated within subsarcolemmal mitochondria and to sites in register with the Z-line. NTT-MMP-2 expression caused a 50% reduction of force development by intact myocardium. However, NTT-MMP-2 expression did not reduce myofilament force development, consistent with the lack of NTT-MMP-2 localization to myofilaments. NTT-MMP-2 expression caused a 50% reduction in the amplitude of Ca²⁺ transients, indicating impaired activation.

Conclusions: Unlike FL-MMP-2, NTT-MMP-2 does not mediate myofilament damage. Instead, NTT-MMP-2 causes impaired myocyte activation, which may involve effects due to localization in mitochondria and/or to transverse tubules affecting Ca²⁺ transients. Thus, FL-MMP-2 and NTT-MMP-2 have discrete intracellular locations and mediate different intracellular damage to cardiac myocytes.

Small and large animal models in cardiac contraction research: advantages and disadvantages

The mammalian heart is responsible for not only pumping blood throughout the body but also adjusting this pumping activity quickly depending upon sudden changes in the metabolic demands of the body. For the most part, the human heart is capable of performing its duties without complications; however, throughout many decades of use, at some point this system encounters problems. Research into the heart's activities during healthy states and during adverse impacts that occur in disease states is necessary in order to strategize novel treatment options to ultimately prolong and improve patients' lives. Animal models are an important aspect of cardiac research where a variety of cardiac processes and therapeutic targets can be studied. However, there are differences between the heart of a human being and an animal and depending on the specific animal, these differences can become more pronounced and in certain cases limiting. There is no ideal animal model available for cardiac research, the use of each animal model is accompanied with its own set of advantages and disadvantages. In this review, we will discuss these advantages and disadvantages of commonly used laboratory animals including mouse, rat, rabbit, canine, swine, and sheep. Since the goal of cardiac research is to enhance our understanding of human health and disease and help improve clinical outcomes, we will also discuss the role of human cardiac tissue in cardiac research. This review will focus on the cardiac ventricular contractile and relaxation kinetics of humans and animal models in order to illustrate these differences.

Rate-dependent Ca2+ signalling underlying the force-frequency response in rat ventricular myocytes: a coupled electromechanical modeling study

BACKGROUND

Rate-dependent effects on the Ca2+ sub-system in a rat ventricular myocyte are investigated. Here, we employ a deterministic mathematical model describing various Ca2+ signalling pathways under voltage clamp (VC) conditions, to better understand the important role of calmodulin (CaM) in modulating the key control variables Ca2+/calmodulin-dependent protein kinase-II (CaMKII), calcineurin (CaN), and cyclic adenosine monophosphate (cAMP) as they affect various intracellular targets. In particular, we study the frequency dependence of the peak force generated by the myofilaments, the force-frequency response (FFR).

METHODS

Our cell model incorporates frequency-dependent CaM-mediated spatially heterogenous interaction of CaMKII and CaN with their principal targets (dihydropyridine (DHPR) and ryanodine (RyR) receptors and the SERCA pump). It also accounts for the rate-dependent effects of phospholamban (PLB) on the SERCA pump; the rate-dependent role of cAMP in up-regulation of the L-type Ca2+ channel (ICa,L); and the enhancement in SERCA pump activity via phosphorylation of PLB.

RESULTS

Our model reproduces positive peak FFR observed in rat ventricular myocytes during voltage-clamp studies both in the presence/absence of cAMP mediated β-adrenergic stimulation. This study provides quantitative insight into the rate-dependence of Ca2+-induced Ca2+-release (CICR) by investigating the frequency-dependence of the trigger current (ICa,L) and RyR-release. It also highlights the relative role of the sodium-calcium exchanger (NCX) and the SERCA pump at higher frequencies, as well as the rate-dependence of sarcoplasmic reticulum (SR) Ca2+ content. A rigorous Ca2+ balance imposed on our investigation of these Ca2+ signalling pathways clarifies their individual roles. Here, we present a coupled electromechanical study emphasizing the rate-dependence of isometric force developed and also investigate the temperature-dependence of FFR.

CONCLUSIONS

Our model provides mechanistic biophysically based explanations for the rate-dependence of CICR, generating useful and testable hypotheses. Although rat ventricular myocytes exhibit a positive peak FFR in the presence/absence of beta-adrenergic stimulation, they show a characteristic increase in the positive slope in FFR due to the presence of Norepinephrine or Isoproterenol. Our study identifies cAMP-mediated stimulation, and rate-dependent CaMKII-mediated up-regulation of ICa,L as the key mechanisms underlying the aforementioned positive FFR.

Intracellular dyssynchrony of diastolic cytosolic [Ca²⁺] decay in ventricular cardiomyocytes in cardiac remodeling and human heart failure

RATIONALE

Synchronized release of Ca²⁺ into the cytosol during each cardiac cycle determines cardiomyocyte contraction.

OBJECTIVE

We investigated synchrony of cytosolic [Ca²⁺] decay during diastole and the impact of cardiac remodeling.

METHODS AND RESULTS

Local cytosolic [Ca²⁺] transients (1-µm intervals) were recorded in murine, porcine, and human ventricular single cardiomyocytes. We identified intracellular regions of slow (slowCaR) and fast (fastCaR) [Ca²⁺] decay based on the local time constants of decay (TAUlocal). The SD of TAUlocal as a measure of dyssynchrony was not related to the amplitude or the timing of local Ca²⁺ release. Stimulation of sarcoplasmic reticulum Ca²⁺ ATPase with forskolin or istaroxime accelerated and its inhibition with cyclopiazonic acid slowed TAUlocal significantly more in slowCaR, thus altering the relationship between SD of TAUlocal and global [Ca²⁺] decay (TAUglobal). Na⁺/Ca²⁺ exchanger inhibitor SEA0400 prolonged TAUlocal similarly in slowCaR and fastCaR. FastCaR were associated with increased mitochondrial density and were more sensitive to the mitochondrial Ca²⁺ uniporter blocker Ru360. Variation in TAUlocal was higher in pig and human cardiomyocytes and higher with increased stimulation frequency (2 Hz). TAUlocal correlated with local sarcomere relengthening. In mice with myocardial hypertrophy after transverse aortic constriction, in pigs with chronic myocardial ischemia, and in end-stage human heart failure, variation in TAUlocal was increased and related to cardiomyocyte hypertrophy and increased mitochondrial density.

CONCLUSIONS

In cardiomyocytes, cytosolic [Ca²⁺] decay is regulated locally and related to local sarcomere relengthening. Dyssynchronous intracellular [Ca²⁺] decay in cardiac remodeling and end-stage heart failure suggests a novel mechanism of cellular contractile dysfunction.

Effect of muscle length on cross-bridge kinetics in intact cardiac trabeculae at body temperature

Dynamic force generation in cardiac muscle, which determines cardiac pumping activity, depends on both the number of sarcomeric cross-bridges and on their cycling kinetics. The Frank–Starling mechanism dictates that cardiac force development increases with increasing cardiac muscle length (corresponding to increased ventricular volume). It is, however, unclear to what extent this increase in cardiac muscle length affects the rate of cross-bridge cycling. Previous studies using permeabilized cardiac preparations, sub-physiological temperatures, or both have obtained conflicting results. Here, we developed a protocol that allowed us to reliably and reproducibly measure the rate of tension redevelopment (k_tr; which depends on the rate of cross-bridge cycling) in intact trabeculae at body temperature. Using K⁺ contractures to induce a tonic level of force, we showed the k_tr was slower in rabbit muscle (which contains predominantly β myosin) than in rat muscle (which contains predominantly α myosin). Analyses of k_tr in rat muscle at optimal length (L_opt) and 90% of optimal length (L₉₀) revealed that k_tr was significantly slower at L_opt (27.7 ± 3.3 and 27.8 ± 3.0 s⁻¹ in duplicate analyses) than at L₉₀ (45.1 ± 7.6 and 47.5 ± 9.2 s⁻¹). We therefore show that k_tr can be measured in intact rat and rabbit cardiac trabeculae, and that the k_tr decreases when muscles are stretched to their optimal length under near-physiological conditions, indicating that the Frank–Starling mechanism not only increases force but also affects cross-bridge cycling kinetics.

Endocardial endothelium is a key determinant of force-frequency relationship in rat ventricular myocardium

We tested the hypothesis that removing endocardial endothelium (EE) negatively impacts the force-frequency relationship (FFR) of ventricular myocardium and dissected the signaling that underlies this phenomenon. EE of rat trabeculae was selectively damaged by brief (<1 s) exposure to 0.1% Triton X-100. Force, intracellular Ca(2+) transient (iCa(2+)), and activity of protein kinase A (PKA) and protein kinase C (PKC) were determined. In control muscles, force and iCa(2+) increased as the stimulation frequency increased in steps of 0.5 Hz up to 3.0 Hz. However, EE-denuded (EED) muscles exhibited a markedly blunted FFR. Neither isoproterenol (ISO; 0.1-5 nmol/l) nor endothelin-1 (ET-1; 10-100 nmol/l) alone restored the slope of FFR in EED muscles. Intriguingly, however, a positive FFR was restored in EED preparations by combining low concentrations of ISO (0.1 nmol/l) and ET-1 (20 nmol/l). In intact muscles, PKA and PKC activity increased proportionally with the increase in frequency. This effect was completely lost in EED muscles. Again, combining ISO and ET-1 fully restored the frequency-dependent rise in PKA and PKC activity in EED muscles. In conclusion, selective damage of EE leads to significantly blunted FFR. A combination of low concentrations of ISO and ET-1 successfully restores FFR in EED muscles. The interdependence of ISO and ET-1 in this process indicates cross-talk between the β1-PKA and ET-1-PKC pathways for a normal (positive) FFR. The results also imply that dysfunction of EE and/or EE-myocyte coupling may contribute to flat (or even negative) FFR in heart failure.

The positive inotropic effect of pyruvate involves an increase in myofilament calcium sensitivity

Pyruvate is a metabolic fuel that is a potent inotropic agent. Despite its unique inotropic and antioxidant properties, the molecular mechanism of its inotropic mechanism is still largely unknown. To examine the inotropic effect of pyruvate in parallel with intracellular calcium handling under near physiological conditions, we measured pH, myofilament calcium sensitivity, developed force, and calcium transients in ultra thin rabbit heart trabeculae at 37 °C loaded iontophoretically with the calcium indicator bis-fura-2. By contrasting conditions of control versus sarcoplasmic reticulum block (with either cyclopiazonic acid and ryanodine or with thapsigargin) we were able to characterize and isolate the effects of pyruvate on sarcoplasmic reticulum calcium handling and developed force. A potassium contracture technique was subsequently utilized to assess the force-calcium relationship and thus the myofilament calcium sensitivity. Pyruvate consistently increased developed force whether or not the sarcoplasmic reticulum was blocked (16.8±3.5 to 24.5±5.1 vs. 6.9±2.6 to 12.5±4.4 mN/mm(2), non-blocked vs. blocked sarcoplasmic reticulum respectively, p<0.001, n = 9). Furthermore, the sensitizing effect of pyruvate on the myofilaments was demonstrated by potassium contractures (EC50 at baseline versus 20 minutes of pyruvate infusion (peak force development) was 701±94 vs. 445±65 nM, p<0.01, n = 6). This study is the first to demonstrate that a leftward shift in myofilament calcium sensitivity is an important mediator of the inotropic effect of pyruvate. This finding can have important implications for future development of therapeutic strategies in the management of heart failure.

Multiphysics model of a rat ventricular myocyte: a voltage-clamp study

BACKGROUND

The objective of this study is to develop a comprehensive model of the electromechanical behavior of the rat ventricular myocyte to investigate the various factors influencing its contractile response.

METHODS

Here, we couple a model of Ca2 + dynamics described in our previous work, with a well-known model of contractile mechanics developed by Rice, Wang, Bers and de Tombe to develop a composite multiphysics model of excitation-contraction coupling. This comprehensive cell model is studied under voltage clamp (VC) conditions, since it allows to focus our study on the elaborate Ca2 + signaling system that controls the contractile mechanism.

RESULTS

We examine the role of various factors influencing cellular contractile response. In particular, direct factors such as the amount of activator Ca2 + available to trigger contraction and the type of mechanical load applied (resulting in isosarcometric, isometric or unloaded contraction) are investigated. We also study the impact of temperature (22 to 38°C) on myofilament contractile response. The critical role of myofilament Ca2 + sensitivity in modulating developed force is likewise studied, as is the indirect coupling of intracellular contractile mechanism with the plasma membrane via the Na + /Ca2 + exchanger (NCX). Finally, we demonstrate a key linear relationship between the rate of contraction and relaxation, which is shown here to be intrinsically coupled over the full range of physiological perturbations.

CONCLUSIONS

Extensive testing of the composite model elucidates the importance of various direct and indirect modulatory influences on cellular twitch response with wide agreement with measured data on all accounts. Thus, the model provides mechanistic insights into whole-cell responses to a wide variety of testing approaches used in studies of cardiac myofilament contractility that have appeared in the literature over the past several decades.

The force-temperature relationship in healthy and dystrophic mouse diaphragm; implications for translational study design

In the field of muscular dystrophy, striated muscle function is often assessed in vitro in dystrophin-deficient mdx mice in order to test the impact of a potential treatment strategy. Although many past studies have assessed diaphragm contractile function at or near room temperature, the diaphragm performs in vivo at 37°C. To improve translation of bench-top results to possible clinical application, we studied temperature-dependence of contractile performance in wild-type (C57BL/10) and mdx muscle strips at temperatures from 25°C to 37°C. Maximal tetanic force in wild-type muscles was higher at 37°C (198 ± 11 vs. 155 ± 9 mN/mm² at 25°C), while the difference between wild-type and mdx was extremely similar: wild-type muscles produced 45.9% and 45.1% more force at 25°C and 37°C respectively. At 37°C twitch contraction kinetics and 50% rise time to tetanic plateau were slower in mdx diaphragm. A fatigue/injury protocol indicated 2-fold fatigue/contraction-induced force deficit in mdx muscles. We conclude that assessment of diaphragm muscle strips can be reliably and reproducibly performed at 37°C.

Impact of hydroxyl radical-induced injury on calcium handling and myofilament sensitivity in isolated myocardium

Hydroxyl radicals (OH) are involved in the pathogenesis of reperfusion injury and are observed in acute heart failure, stroke, and myocardial infarction. Two different subcellular defects are involved in the pathogenesis of OH injury, deranged calcium handling, and alterations of myofilament responsiveness, but their temporal impact on contractile function is not resolved. Initially, after brief OH exposure, there is a corresponding marked increase in diastolic calcium and diastolic force. We followed these parameters until a new steady-state level was reached at ~45 min post-OH exposure. At this new baseline, diastolic calcium had returned to near-normal, pre-OH levels, whereas diastolic force remained markedly elevated. An increased calcium sensitivity was observed at the new baseline after OH-induced injury compared with the pre-OH state. The acute injury that occurs after OH exposure is mainly due to calcium overload, while the later sustained myocardial dysfunction is mainly due to the altered/increased myofilament responsiveness.

The rates of Ca2+ dissociation and cross-bridge detachment from ventricular myofibrils as reported by a fluorescent cardiac troponin C

The rate-limiting step of cardiac muscle relaxation has been proposed to reside in the myofilament. Both the rates of cross-bridge detachment and Ca(2+) dissociation from troponin C (TnC) have been hypothesized to rate-limit myofilament inactivation. In this study we used a fluorescent TnC to measure both the rate of Ca(2+) dissociation from TnC and the rate of cross-bridge detachment from several different species of ventricular myofibrils. The fluorescently labeled TnC was sensitive to both Ca(2+) dissociation and cross-bridge detachment at low Ca(2+) (presence of EGTA), allowing for a direct comparison between the two proposed rates of myofilament inactivation. Unlike Ca(2+) dissociation from TnC, cross-bridge detachment varied in myofibrils from different species and was rate-limited by ADP release. At subphysiological temperatures (<20 °C), the rate of Ca(2+) dissociation from TnC was faster than the rate of cross-bridge detachment in the presence of ADP. These results support the hypothesis that cross-bridge detachment rate-limits relaxation. However, Ca(2+) dissociation from TnC was not as temperature-sensitive as cross-bridge detachment. At a near physiological temperature (35 °C) and ADP, the rate of cross-bridge detachment may actually be faster than the rate of Ca(2+) dissociation. This provides evidence that there may not be a simple, single rate-limiting step of myofilament inactivation.

Contractile parameters and occurrence of alternans in isolated rat myocardium at supra-physiological stimulation frequency

The cardiac refractory period prevents the heart from tetanic activation that is typically used in noncardiac striated muscle tissue. To what extent the refractory period prevents successive action potentials to activate the excitation-contraction coupling process and contractile machinery at supra-physiological rates, such as those present during ventricular fibrillation, is unknown. Using multicellular trabeculae isolated from rat hearts, we studied amplitude and kinetics of contraction at rates well above the normal in vivo rat heart range. We show that even at twice the maximal heart rate of the rat, little or no mechanical instability is observed; twitch contractions are at steady state, albeit with an elevated active diastolic force. Although the amplitude of contraction increased within in vivo heart rates (positive force-frequency response), at frequencies beyond the maximal heart rate (10-30 Hz) a steady decline of contractile amplitude is observed. Not until 30 Hz do the majority of the isolated muscle preparations show mechanical alternans, where strong and weak beats alternate. Interestingly, unlike striated limb skeletal muscle, fusing of twitch contractions did not cause a continuous increase in peak force: at frequencies of 10 Hz and above, systolic force declines with relatively little elevation in diastolic force. Contractile kinetics continued to accelerate, from 1 Hz up to 30 Hz, whereas the relative speed of contraction and relaxation remained closely coupled, reflected by a singular linear relationship between the maximal and minimal derivative of force (dF/dt). We conclude that cardiac muscle can produce mechanically stable steady-state contractions at supra-physiological pacing rates, while these contractions continue to decline in amplitude and increase in diastolic force past maximal heart rate.

Micromechanical thermal assays of Ca2+-regulated thin-filament function and modulation by hypertrophic cardiomyopathy mutants of human cardiac troponin

Microfabricated thermoelectric controllers can be employed to investigate mechanisms underlying myosin-driven sliding of Ca(2+)-regulated actin and disease-associated mutations in myofilament proteins. Specifically, we examined actin filament sliding-with or without human cardiac troponin (Tn) and α-tropomyosin (Tm)-propelled by rabbit skeletal heavy meromyosin, when temperature was varied continuously over a wide range (~20-63°C). At the upper end of this temperature range, reversible dysregulation of thin filaments occurred at pCa 9 and 5; actomyosin function was unaffected. Tn-Tm enhanced sliding speed at pCa 5 and increased a transition temperature (T(t)) between a high activation energy (E(a)) but low temperature regime and a low E(a) but high temperature regime. This was modulated by factors that alter cross-bridge number and kinetics. Three familial hypertrophic cardiomyopathy (FHC) mutations, cTnI R145G, cTnI K206Q, and cTnT R278C, cause dysregulation at temperatures ~5-8°C lower; the latter two increased speed at pCa 5 at all temperatures.