The time-course of Ca2+ exchange with calmodulin, troponin, parvalbumin, and myosin in response to transient increases in Ca2+

A detailed mathematical model of the human atrial cardiomyocyte: integration of electrophysiology and cardiomechanics

Mechano-electric regulations (MER) play an important role in the maintenance of cardiac performance. Mechano-calcium and mechano-electric feedback (MCF and MEF) pathways adjust the cardiomyocyte contractile force according to mechanical perturbations and affects electro-mechanical coupling. MER integrates all these regulations in one unit resulting in a complex phenomenon. Computational modelling is a useful tool to accelerate the mechanistic understanding of complex experimental phenomena. We have developed a novel model that integrates the MER loop for human atrial cardiomyocytes with proper consideration of feedforward and feedback pathways. The model couples a modified version of the action potential (AP) Koivumäki model with the contraction model by Quarteroni group. The model simulates iso-sarcometric and isometric twitches and the feedback effects on AP and Ca2+-handling. The model showed a biphasic response of Ca2+ transient (CaT) peak to increasing pacing rates and highlights the possible mechanisms involved. The model has shown a shift of the threshold for AP and CaT alternans from 4.6 to 4 Hz under post-operative atrial fibrillation, induced by depressed SERCA activity. The alternans incidence was dependent on a chain of mechanisms including RyRs availability time, MCF coupling, CaMKII phosphorylation, and the stretch levels. As a result, the model predicted a 10% slowdown of conduction velocity for a 20% stretch, suggesting a role of stretch in creation of substrate formation for atrial fibrillation. Overall, we conclude that the developed model provides a physiological CaT followed by a physiological twitch. This model can open pathways for the future studies of human atrial electromechanics. KEY POINTS: With the availability of human atrial cellular data, interest in atrial-specific model integration has been enhanced. We have developed a detailed mathematical model of human atrial cardiomyocytes including the mechano-electric regulatory loop. The model has gone through calibration and evaluation phases against a wide collection of available human in-vitro data. The usefulness of the model for analysing clinical problems has been preliminaryly tested by simulating the increased incidence of Ca2+ transient and action potential alternans at high rates in post-operative atrial fibrillation condition. The model determines the possible role of mechano-electric feedback in alternans incidence, which can increase vulnerability to atrial arrhythmias by varying stretch levels. We found that our physiologically accurate description of Ca2+ handling can reproduce many experimental phenomena and can help to gain insights into the underlying pathophysiological mechanisms.

Evolution of Bioelectric Membrane Potentials: Implications in Cancer Pathogenesis and Therapeutic Strategies

Electrophysiology typically deals with the electrical properties of excitable cells like neurons and muscles. However, all other cells (non-excitable) also possess bioelectric membrane potentials for intracellular and extracellular communications. These membrane potentials are generated by different ions present in fluids available in and outside the cell, playing a vital role in communication and coordination between the cell and its organelles. Bioelectric membrane potential variations disturb cellular ionic homeostasis and are characteristic of many diseases, including cancers. A rapidly increasing interest has emerged in sorting out the electrophysiology of cancer cells. Compared to healthy cells, the distinct electrical properties exhibited by cancer cells offer a unique way of understanding cancer development, migration, and progression. Decoding the altered bioelectric signals influenced by fluctuating electric fields benefits understanding cancer more closely. While cancer research has predominantly focussed on genetic and molecular traits, the delicate area of electrophysiological characteristics has increasingly gained prominence. This review explores the historical exploration of electrophysiology in the context of cancer cells, shedding light on how alterations in bioelectric membrane potentials, mediated by ion channels and gap junctions, contribute to the pathophysiology of cancer.

Ca2+ storage function is altered in the sarcoplasmic reticulum of skeletal muscle lacking mitsugumin 23

Mitsugumin 23 (MG23) has been identified as a ball-shaped cation channel in the sarcoplasmic reticulum (SR) but its physiological role remains unclear. This study aimed to examine the contribution of MG23 to Ca2+ storage function in skeletal muscle by using Mg23-knockout (Mg23-/-) mice. There was no difference in the isometric specific force of the extensor digitorum longus (EDL) and soleus (SOL) muscles between Mg23-/- and wild-type (Wt) mice. In Mg23-/- mice, the calsequestrin 2 content in the EDL muscle and SR Ca2+-ATPase 2 content in the SOL were increased. We have examined SR and myofibril functions using mechanically skinned fibers and determined their fiber types based on the response to Sr2+, which showed that Mg23-/- mice, compared with Wt, had: 1) elevated total Ca2+ content in the membranous components including SR, mitochondria, and transverse tubular system referred to as endogenous Ca2+ content, in both type I and II fibers of the EDL and SOL; 2) increased maximal Ca2+ content in both type I and II fibers of the EDL and SOL; 3) decreased SR Ca2+ leakage in type I fibers of the SOL; and 4) enhanced SR Ca2+ uptake in type I fibers of the SOL, although myofibril function was not different in both type I and II fibers of the SOL and EDL muscles. These results suggest that MG23 decreases SR Ca2+ storage in both type I and type II fibers, likely due to increased SR Ca2+ leakage.NEW & NOTEWORTHY The function of calcium storage within sarcoplasmic reticulum (SR) plays a pivotal role in influencing the health and disease states of skeletal muscle. In the present study, we demonstrated that mitsgumin 23, a novel non-selective cation channel, modifies SR Ca2+ storage in skeletal muscle fibers. These findings provide valuable insights into the physiological regulation of Ca2+ in skeletal muscle, offering significant potential for uncovering the mechanisms underlying muscle fatigue, muscle adaptation, and muscle diseases.

Physiology of intracellular calcium buffering

Calcium signaling underlies much of physiology. Almost all the Ca2+ in the cytoplasm is bound to buffers, with typically only ∼1% being freely ionized at resting levels in most cells. Physiological Ca2+ buffers include small molecules and proteins, and experimentally Ca2+ indicators will also buffer calcium. The chemistry of interactions between Ca2+ and buffers determines the extent and speed of Ca2+ binding. The physiological effects of Ca2+ buffers are determined by the kinetics with which they bind Ca2+ and their mobility within the cell. The degree of buffering depends on factors such as the affinity for Ca2+, the Ca2+ concentration, and whether Ca2+ ions bind cooperatively. Buffering affects both the amplitude and time course of cytoplasmic Ca2+ signals as well as changes of Ca2+ concentration in organelles. It can also facilitate Ca2+ diffusion inside the cell. Ca2+ buffering affects synaptic transmission, muscle contraction, Ca2+ transport across epithelia, and the killing of bacteria. Saturation of buffers leads to synaptic facilitation and tetanic contraction in skeletal muscle and may play a role in inotropy in the heart. This review focuses on the link between buffer chemistry and function and how Ca2+ buffering affects normal physiology and the consequences of changes in disease. As well as summarizing what is known, we point out the many areas where further work is required.

Optimization and comparison of two methods for spike train estimation in an unfused tetanic contraction of low threshold motor units

BACKGROUND

Recent findings have shown that imaging voluntarily activated motor units (MUs) by decomposing ultrasound-based displacement images provides estimates of unfused tetanic signals evoked by spinal motoneurons' neural discharges (spikes). Two methods have been suggested to estimate its spike trains: band-pass filter (BPM) and Haar wavelet transform (HWM). However, the methods' optimal parameters and which method performs the best are unknown. This study will answer these questions.

METHOD

HWM and BPM were optimized using simulations. Their performance was evaluated based on simulations and 21 experimental datasets, considering their rate of agreement, spike offset, and spike offset variability to the simulated or experimental spikes.

RESULTS

A range of parameter sets that resulted in the highest possible agreement with simulated spikes was provided. Both methods highly agreed with simulated and experimental spikes, but HWM was a better spike estimation method than BPM because it had a higher agreement, less bias, and less variation (p < 0.001).

CONCLUSIONS

The optimized HWM will be an important contributor to further developing the identification and analysis of MUs using imaging, providing indirect access to the neural drive of the spinal cord to the muscle by the unfused tetanic signals.

Sarcomere dynamics revealed by a myofilament integrated FRET-based biosensor in live skeletal muscle fibers

The sarcomere is the functional unit of skeletal muscle, essential for proper contraction. Numerous acquired and inherited myopathies impact sarcomere function causing clinically significant disease. Mechanistic investigations of sarcomere activation have been challenging to undertake in the context of intact, live skeletal muscle fibers during real time physiological twitch contractions. Here, a skeletal muscle specific, intramolecular FRET-based biosensor was designed and engineered into fast skeletal muscle troponin C (TnC) to investigate the dynamics of sarcomere activation. In transgenic animals, the TnC biosensor incorporated into the skeletal muscle fiber sarcomeres by stoichiometric replacement of endogenous TnC and did not alter normal skeletal muscle contractile form or function. In intact single adult skeletal muscle fibers, real time twitch contractile data showed the TnC biosensor transient preceding the peak amplitude of contraction. Importantly, under physiological temperatures, inactivation of the TnC biosensor transient decayed significantly more slowly than the Ca²⁺ transient and contraction. The uncoupling of the TnC biosensor transient from the Ca²⁺ transient indicates the biosensor is not functioning as a Ca²⁺ transient reporter, but rather reports dynamic sarcomere activation/ inactivation that, in turn, is due to the ensemble effects of multiple activating ligands within the myofilaments. Together, these findings provide the foundation for implementing this new biosensor in future physiological studies investigating the mechanism of activation of the skeletal muscle sarcomere in health and disease.

Cytosolic Ca2+-dependent Ca2+ release activity primarily determines the ER Ca2+ level in cells expressing the CPVT-linked mutant RYR2

Type 2 ryanodine receptor (RYR2) is a cardiac Ca2+ release channel in the ER. Mutations in RYR2 are linked to catecholaminergic polymorphic ventricular tachycardia (CPVT). CPVT is associated with enhanced spontaneous Ca2+ release, which tends to occur when [Ca2+]ER reaches a threshold. Mutations lower the threshold [Ca2+]ER by increasing luminal Ca2+ sensitivity or enhancing cytosolic [Ca2+] ([Ca2+]cyt)-dependent activity. Here, to establish the mechanism relating the change in [Ca2+]cyt-dependent activity of RYR2 and the threshold [Ca2+]ER, we carried out cell-based experiments and in silico simulations. We expressed WT and CPVT-linked mutant RYR2s in HEK293 cells and measured [Ca2+]cyt and [Ca2+]ER using fluorescent Ca2+ indicators. CPVT RYR2 cells showed higher oscillation frequency and lower threshold [Ca2+]ER than WT cells. The [Ca2+]cyt-dependent activity at resting [Ca2+]cyt, Arest, was greater in CPVT mutants than in WT, and we found an inverse correlation between threshold [Ca2+]ER and Arest. In addition, lowering RYR2 expression increased the threshold [Ca2+]ER and a product of Arest, and the relative expression level for each mutant correlated with threshold [Ca2+]ER, suggesting that the threshold [Ca2+]ER depends on the net Ca2+ release rate via RYR2. Modeling reproduced Ca2+ oscillations with [Ca2+]cyt and [Ca2+]ER changes in WT and CPVT cells. Interestingly, the [Ca2+]cyt-dependent activity of specific mutations correlated with the age of disease onset in patients carrying them. Our data suggest that the reduction in threshold [Ca2+]ER for spontaneous Ca2+ release by CPVT mutation is explained by enhanced [Ca2+]cyt-dependent activity without requiring modulation of the [Ca2+]ER sensitivity of RYR2.

Monoclonal Antibodies as Probes to Study Ligand-Induced Conformations of Troponin Subunits

Striated muscle contraction and relaxation is regulated by Ca2+ at the myofilament level via conformational modulations of the troponin complex. To understand the structure-function relationship of troponin in normal muscle and in myopathies, it is necessary to study the functional effects of troponin isoforms and mutations at the level of allosteric conformations of troponin subunits. Traditional methodologies assessing such conformational studies are laborious and require significant amounts of purified protein, while many current methodologies require non-physiological conditions or labeling of the protein, which may affect their physiological conformation and function. To address these issues, we developed a novel approach using site-specific monoclonal antibodies (mAb) as molecular probes to detect and monitor conformational changes of proteins. Here, we present examples for its application in studies of two subunits of troponin: the Ca2+-binding subunit, TnC, and the tropomyosin-binding/thin filament-anchoring subunit, TnT. Studies using a high-throughput microplate assay are compared with that using localized surface plasmon resonance (LSPR) to demonstrate the effectiveness of using mAb probes to assess ligand-induced conformations of troponin subunits in physiological conditions. The assays utilize relatively small amounts of protein and are free of protein modification, which may bias results. Detailed methodologies using various monoclonal antibodies (mAbs) are discussed with considerations for the optimization of assay conditions and the broader application in studies of other proteins as well as in screening of therapeutic reagents that bind a specific target site with conformational and functional effects.

Lower troponin expression in the right ventricle of rats explains interventricular differences in E-C coupling

Despite distinctive functional and anatomic differences, a precise understanding of the cardiac interventricular differences in excitation–contraction (E–C) coupling mechanisms is still lacking. Here, we directly compared rat right and left cardiomyocytes (RVCM and LVCM). Whole-cell patch clamp, the IonOptix system, and fura-2 fluorimetry were used to measure electrical properties (action potential and ionic currents), single-cell contractility, and cytosolic Ca²⁺ ([Ca²⁺]_i), respectively. Myofilament proteins were analyzed by immunoblotting. RVCM showed significantly shorter action potential duration (APD) and higher density of transient outward K⁺ current (I_to). However, the triggered [Ca²⁺]_i change (Ca²⁺ transient) was not different, while the decay rate of the Ca²⁺ transient was slower in RVCM. Although the relaxation speed was also slower, the sarcomere shortening amplitude (ΔSL) was smaller in RVCM. SERCA activity was ∼60% lower in RVCM, which is partly responsible for the slower decay of the Ca²⁺ transient. Immunoblot analysis revealed lower expression of the cardiac troponin complex (cTn) in RVCM, implying a smaller Ca²⁺ buffering capacity (κ_S), which was proved by in situ analysis. The introduction of these new levels of cTn, I_to, and SERCA into a mathematical model of rat LVCM reproduced the similar Ca²⁺ transient, slower Ca²⁺ decay, shorter APD, and smaller ΔSL of RVCM. Taken together, these data show reduced expression of cTn proteins in the RVCM, which provides an explanation for the interventricular difference in the E–C coupling kinetics.

Role of parvalbumin in fatigue-induced changes in force and cytosolic calcium transients in intact single mouse myofibers

One of the most important cytosolic Ca²⁺ buffers present in mouse fast-twitch myofibers, but not in human myofibers, is parvalbumin (PV). Previous work using conventional PV knockout mice suggests that lifelong PV ablation increases fatigue resistance, possibly due to compensations in mitochondrial volume. In this work, PV gene ablation was induced only in adult mice (PV-KO), and contractile and cytosolic Ca²⁺ responses during fatigue were studied in isolated muscle and intact single myofibers. Results were compared to control littermates (PV-Ctr). We hypothesized that the reduced myofiber cytosolic Ca²⁺ buffering developed only in adult PV-KO mice leads to a larger cytosolic free Ca²⁺ concentration ([Ca²⁺]_c) during repetitive contractions, increasing myofiber fatigue resistance. Extensor digitorum longus (EDL) muscles from PV-KO mice had higher force in unfused stimulations (~50%, P<0.05) and slowed relaxation (~46% higher relaxation time, P<0.05) vs PV-Ctr, but muscle fatigue resistance or fatigue-induced changes in relaxation were not different between genotypes (P>0.05). In intact single myofibers from flexor digitorum brevis (FDB) muscles, basal and tetanic [Ca²⁺]_c during fatiguing contractions were higher in PV-KO (P<0.05), accompanied by a greater slowing in estimated sarcoplasmic reticulum (SR) Ca²⁺ pumping vs PV-Ctr myofibers (~84% reduction, P<0.05), but myofiber fatigue resistance was not different between genotypes (P>0.05). Our results demonstrate that although the estimated SR Ca²⁺ uptake was accelerated in PV-KO, the total energy demand by the major energy consumers in myofibers, the cross-bridges and SR Ca²⁺ ATPase, were not altered enough to affect the energy supply for contractions, and therefore fatigue resistance remained unaffected.

A multiscale sliding filament model of lymphatic muscle pumping

The lymphatics maintain fluid balance by returning interstitial fluid to veins via contraction/compression of vessel segments with check valves. Disruption of lymphatic pumping can result in a condition called lymphedema with interstitial fluid accumulation. Lymphedema treatments are often ineffective, which is partially attributable to insufficient understanding of specialized lymphatic muscle lining the vessels. This muscle exhibits cardiac-like phasic contractions and smooth muscle-like tonic contractions to generate and regulate flow. To understand the relationship between this sub-cellular contractile machinery and organ-level pumping, we have developed a multiscale computational model of phasic and tonic contractions in lymphatic muscle and coupled it to a lymphangion pumping model. Our model uses the sliding filament model (Huxley in Prog Biophys Biophys Chem 7:255-318, 1957) and its adaptation for smooth muscle (Mijailovich in Biophys J 79(5):2667-2681, 2000). Multiple structural arrangements of contractile components and viscoelastic elements were trialed but only one provided physiologic results. We then coupled this model with our previous lumped parameter model of the lymphangion to relate results to experiments. We show that the model produces similar pressure, diameter, and flow tracings to experiments on rat mesenteric lymphatics. This model provides the first estimates of lymphatic muscle contraction energetics and the ability to assess the potential effects of sub-cellular level phenomena such as calcium oscillations on lymphangion outflow. The maximum efficiency value predicted (40%) is at the upper end of estimates for other muscle types. Spontaneous calcium oscillations during diastole were found to increase outflow up to approximately 50% in the range of frequencies and amplitudes tested.

Ryanodine receptor leak triggers fiber Ca2+ redistribution to preserve force and elevate basal metabolism in skeletal muscle

Muscle contraction depends on tightly regulated Ca²⁺ release. Aberrant Ca²⁺ leak through ryanodine receptor 1 (RyR1) on the sarcoplasmic reticulum (SR) membrane can lead to heatstroke and malignant hyperthermia (MH) susceptibility, as well as severe myopathy. However, the mechanism by which Ca²⁺ leak drives these pathologies is unknown. Here, we investigate the effects of four mouse genotypes with increasingly severe RyR1 leak in skeletal muscle fibers. We find that RyR1 Ca²⁺ leak initiates a cascade of events that cause precise redistribution of Ca²⁺ among the SR, cytoplasm, and mitochondria through altering the Ca²⁺ permeability of the transverse tubular system membrane. This redistribution of Ca²⁺ allows mice with moderate RyR1 leak to maintain normal function; however, severe RyR1 leak with RYR1 mutations reduces the capacity to generate force. Our results reveal the mechanism underlying force preservation, increased ATP metabolism, and susceptibility to MH in individuals with gain-of-function RYR1 mutations.

Inhibition of Peripheral ERK Signaling Ameliorates Persistent Muscle Pain Around Trigger Points in Rats

The purpose of this study was to investigate whether the ERK signaling pathway was involved in ameliorating chronic myofascial hyperalgesia from contused gastrocnemius muscle in rats. We established an animal model associated with myofascial pain syndrome and described the mechanism of muscle pain in an animal model. Changes in the mechanical pain threshold were observed 0.5, 1, 2, 3, 4, 5, 8, 12, 18, and 24 h after ERK inhibitor injection around myofascial trigger points (MTrPs) of the gastrocnemius muscle in rats. Morphological changes in gastrocnemius muscle cells were observed by hematoxylin and eosin (H&E) staining. ERK signaling pathway activation was detected through immunohistochemistry and Western blotting. The main morphological characteristics of injured muscle fibers around MTrPs include gathered circular or elliptical shapes of different sizes in the cross-section and continuous inflated and tapering fibers in the longitudinal section. After intramuscular injection of U0126 (ERK inhibitor), the mechanical pain threshold significantly increased. The reduction in mechanical hyperalgesia was accompanied by reduced ERK protein phosphorylation, myosin light chain kinase (MLCK) protein, p-MLC protein expression, and the cross-sectional area of skeletal muscle cells around MTrPs. An ERK inhibitor contributed to the attenuation of mechanical hyperalgesia in the rat myofascial pain model, and the increase in pain threshold may be related to MLCK downregulation and other related contraction-associated proteins by ERK.

SpiCee: A Genetic Tool for Subcellular and Cell-Specific Calcium Manipulation

Calcium is a second messenger crucial to a myriad of cellular processes ranging from regulation of metabolism and cell survival to vesicle release and motility. Current strategies to directly manipulate endogenous calcium signals lack cellular and subcellular specificity. We introduce SpiCee, a versatile and genetically encoded chelator combining low- and high-affinity sites for calcium. This scavenger enables altering endogenous calcium signaling and functions in single cells in vitro and in vivo with biochemically controlled subcellular resolution. SpiCee paves the way to investigate local calcium signaling in vivo and directly manipulate this second messenger for therapeutic use.

Binding of calcium and magnesium to human cardiac troponin C

Cardiac muscle thin filaments are composed of actin, tropomyosin, and troponin that change conformation in response to Ca²⁺ binding, triggering muscle contraction. Human cardiac troponin C (cTnC) is the Ca²⁺-sensing component of the thin filament. It contains structural sites (III/IV) that bind both Ca²⁺ and Mg²⁺ and a regulatory site (II) that has been thought to bind only Ca²⁺. Binding of Ca²⁺ at this site initiates a series of conformational changes that culminate in force production. However, the mechanisms that underpin the regulation of binding at site II remain unclear. Here, we have quantified the interaction between site II and Ca²⁺/Mg²⁺ through isothermal titration calorimetry and thermodynamic integration simulations. Direct and competitive binding titrations with WT N-terminal cTnC and full-length cTnC indicate that physiologically relevant concentrations of both Ca²⁺/Mg²⁺ interacted with the same locus. Moreover, the D67A/D73A N-terminal cTnC construct in which two coordinating residues within site II were removed was found to have significantly reduced affinity for both cations. In addition, 1 mM Mg²⁺ caused a 1.4-fold lower affinity for Ca²⁺. These experiments strongly suggest that cytosolic-free Mg²⁺ occupies a significant population of the available site II. Interaction of Mg²⁺ with site II of cTnC likely has important functional consequences for the heart both at baseline as well as in diseased states that decrease or increase the availability of Mg²⁺, such as secondary hyperparathyroidism or ischemia, respectively.

The ForceLAB simulator: Application to the comparison of current models of cardiomyocyte contraction

Mathematical models are useful tools in the study of physiological phenomena. However, due to differences in assumptions and formulations, discrepancy in simulations may occur. Among the models for cardiomyocyte contraction based on Huxley's cross-bridge cycling, those proposed by Negroni and Lascano (NL) and Rice et al. (RWH) are the most frequently used. This study was aimed at developing a computational tool, ForceLAB, which allows implementing different contraction models and modifying several functional parameters. As an application, electrically-stimulated twitches triggered by an equal Ca²⁺ input and steady-state force x pCa relationship (pCa = -log of the molar free Ca²⁺ concentration) simulated with the NL and RWH models were compared. The equilibrium Ca²⁺-troponin C (TnC) dissociation constant (K_d) was modified by changing either the association (k_on) or the dissociation (k_off) rate constant. With the NL model, raising K_d by either maneuver decreased monotonically twitch amplitude and duration, as expected. With the RWH model, in contrast, the same K_d variation caused increase or decrease of peak force depending on which rate constant was modified. Additionally, force x pCa curves simulated using Ca²⁺ binding constants estimated in cardiomyocytes bearing wild-type and mutated TnC were compared to curves previously determined in permeabilized fibers. Mutations increased k_on and k_off, and decreased K_d. Both models produced curves fairly comparable to the experimental ones, although sensitivity to Ca²⁺ was greater, especially with RWH model. The NL model reproduced slightly better the qualitative changes associated with the mutations. It is expected that this tool can be useful for teaching and investigation.

Buffering and total calcium levels determine the presence of oscillatory regimes in cardiac cells

Calcium oscillations and waves induce depolarization in cardiac cells which are believed to cause life-threathening arrhythimas. In this work, we study the conditions for the appearance of calcium oscillations in both a detailed subcellular model of calcium dynamics and a minimal model that takes into account just the minimal ingredients of the calcium toolkit. To avoid the effects of homeostatic changes and the interaction with the action potential we consider the somewhat artificial condition of a cell without pacing and with no calcium exchange with the extracellular medium. Both the full subcellular model and the minimal model present the same scenarios depending on the calcium load: two stationary states, one with closed ryanodine receptors (RyR) and most calcium in the cell stored in the sarcoplasmic reticulum (SR), and another, with open RyRs and a depleted SR. In between, calcium oscillations may appear. The robustness of these oscillations is determined by the amount of calsequestrin (CSQ). The lack of this buffer in the SR enhances the appearance of oscillations. The minimal model allows us to relate the stability of the oscillating state to the nullcline structure of the system, and find that its range of existence is bounded by a homoclinic and a Hopf bifurcation, resulting in a sudden transition to the oscillatory regime as the cell calcium load is increased. Adding a small amount of noise to the RyR behavior increases the parameter region where oscillations appear and provides a gradual transition from the resting state to the oscillatory regime, as observed in the subcellular model and experimentally.

In cardiac cells, calcium plays a very important role. An increase in calcium levels is the trigger used by the cell to initiate contraction. Besides, calcium modulates several transmembrane currents, affecting the cell transmembrane potential. Thus, dysregulations in calcium handling have been associated with the appearance of arrhythmias. Often, this dysregulation results in the appearance of periodic calcium waves or global oscillations, providing a pro-arrhythmic substrate. In this paper, we study the onset of calcium oscillations in cardiac cells using both a detailed subcellular model of calcium dynamics and a minimal model that takes into account the essential ingredients of the calcium toolkit. Both reproduce the main experimental results and link this behavior with the presence of different steady-state solutions and bifurcations that depend on the total amount of calcium in the cell and in the level of buffering present. We expect that this work will help to clarify the conditions under which calcium oscillations appear in cardiac myocytes and, therefore, will represent a step further in the understanding of the origin of cardiac arrhythmias.

Sarcomere integrated biosensor detects myofilament-activating ligands in real time during twitch contractions in live cardiac muscle

The sarcomere is the functional unit of cardiac muscle, essential for normal heart function. To date, it has not been possible to study, in real time, thin filament-based activation dynamics in live cardiac muscle. We report here results from a cardiac troponin C (TnC) FRET-based biosensor integrated into the cardiac sarcomere via stoichiometric replacement of endogenous TnC. The TnC biosensor provides, for the first time, evidence of multiple thin filament activating ligands, including troponin I interfacing with TnC and cycling myosin, during a cardiac twitch. Results show that the TnC FRET biosensor transient significantly precedes that of peak twitch force. Using small molecules and genetic modifiers known to alter sarcomere activation, independently of the intracellular Ca²⁺ transient, the data show that the TnC biosensor detects significant effects of the troponin I switch domain as a sarcomere-activating ligand. Interestingly, the TnC biosensor also detected the effects of load-dependent altered myosin cycling, as shown by a significant delay in TnC biosensor transient inactivation during the isometric twitch. In addition, the TnC biosensor detected the effects of myosin as an activating ligand during the twitch by using a small molecule that directly alters cross-bridge cycling, independently of the intracellular Ca²⁺ transient. Collectively, these results aid in illuminating the basis of cardiac muscle contractile activation with implications for gene, protein, and small molecule-based strategies designed to target the sarcomere in regulating beat-to-beat heart performance in health and disease.

Influence of the tubular network on the characteristics of calcium transients in cardiac myocytes

Transverse and axial tubules (TATS) are an essential ingredient of the excitation-contraction machinery that allow the effective coupling of L-type Calcium Channels (LCC) and ryanodine receptors (RyR2). They form a regular network in ventricular cells, while their presence in atrial myocytes is variable regionally and among animal species We have studied the effect of variations in the TAT network using a bidomain computational model of an atrial myocyte with variable density of tubules. At each z-line the t-tubule length is obtained from an exponential distribution, with a given mean penetration length. This gives rise to a distribution of t-tubules in the cell that is characterized by the fractional area (F.A.) occupied by the t-tubules. To obtain consistent results, we average over different realizations of the same mean penetration length. To this, in some simulations we add the effect of a network of axial tubules. Then we study global properties of calcium signaling, as well as regional heterogeneities and local properties of sparks and RyR2 openings. In agreement with recent experiments in detubulated ventricular and atrial cells, we find that detubulation reduces the calcium transient and synchronization in release. However, it does not affect sarcoplasmic reticulum (SR) load, so the decrease in SR calcium release is due to regional differences in Ca²⁺ release, that is restricted to the cell periphery in detubulated cells. Despite the decrease in release, the release gain is larger in detubulated cells, due to recruitment of orphaned RyR2s, i.e, those that are not confronting a cluster of LCCs. This probably provides a safeguard mechanism, allowing physiological values to be maintained upon small changes in the t-tubule density. Finally, we do not find any relevant change in spark properties between tubulated and detubulated cells, suggesting that the differences found in experiments could be due to differential properties of the RyR2s in the membrane and in the t-tubules, not incorporated in the present model. This work will help understand the effect of detubulation, that has been shown to occur in disease conditions such as heart failure (HF) in ventricular cells, or atrial fibrillation (AF) in atrial cells.

Calcium Buffering in the Heart in Health and Disease

Changes of intracellular Ca²⁺ concentration regulate many aspects of cardiac myocyte function. About 99% of the cytoplasmic calcium in cardiac myocytes is bound to buffers, and their properties will therefore have a major influence on Ca²⁺ signaling. This article considers the fundamental properties and identities of the buffers and how to measure them. It reviews the effects of buffering on the systolic Ca²⁺ transient and how this may change physiologically, and in heart failure and both atrial and ventricular arrhythmias, as well. It is concluded that the consequences of this strong buffering may be more significant than currently appreciated, and a fuller understanding is needed for proper understanding of cardiac calcium cycling and contractility.

Components of activation heat in skeletal muscle

Activation heat (q_A) production by muscle is the thermal accompaniment of the release of Ca²⁺ from the sarcoplasmic reticulum (SR) into the cytoplasm, its interactions with regulatory proteins and other cytoplasmic Ca²⁺ buffers and its return to the SR. The contribution of different Ca²⁺-related reactions to q_A is difficult to determine empirically and therefore, for this study, a mathematical model was developed to describe Ca²⁺ movements and accompanying thermal changes in muscle fibres in response to stimulation. The major sources of heat within a few milliseconds of the initiation of Ca²⁺ release are Ca²⁺ binding to Tn and Pv. Ca²⁺ binding to ATP produces a relatively small amount of heat. Ca²⁺ dissociation from ATP and Tn, with heat absorption, are of similar time course to the decline of force. In muscle lacking Pv (e.g. mouse soleus), Ca²⁺ is then rapidly pumped into the SR. In muscles with Pv, Ca²⁺ that dissociates from Tn and ATP binds to Pv and then dissociates slowly (over 10 s of seconds) and is then pumped into the SR; the net effect of these two processes is heat absorption. It is proposed that this underlies Hill's "negative delayed heat". After all the Ca²⁺ is returned to the SR, q_A is proportional to the amount of Ca²⁺ released into the cytoplasm. In muscles with Pv this is 20-60 s after Ca²⁺ release; in muscles without Pv, all Ca²⁺ is returned to the SR soon after the end of force relaxation.

Y-27632 preconditioning enhances transplantation of human-induced pluripotent stem cell-derived cardiomyocytes in myocardial infarction mice

Aims

The effectiveness of cell-based treatments for regenerative myocardial therapy is limited by low rates of cell engraftment. Y-27632 inhibits Rho-associated protein kinase (ROCK), which regulates the cytoskeletal changes associated with cell adhesion, and has been used to protect cultured cells during their passaging. Here, we investigated whether preconditioning of cardiomyocytes, derived from human-induced pluripotent stem cells (hiPSC-CM), with Y-27632 improves their survival and engraftment in a murine model of acute myocardial infarction (MI).

Methods and results

After MI induction, mice were subjected to intramyocardial injections of phosphate-buffered saline, hiPSC-CM cultured under standard conditions (hiPSC-CM-RI), or Y-27632-preconditioned hiPSC-CM (hiPSC-CM+RI). The resulting engraftment rate calculated 4 weeks after implantation was significantly higher and the abundance of apoptotic transplanted cells was significantly lower in hiPSC-CM+RI recipients than in hiPSC-CM-RI animals. In cultured hiPSC-CM, Y-27632-preconditioning reversibly reduced contractile activity and the expression of troponin genes, while increasing their attachment to an underlying mouse cardiomyocyte (HL1) monolayer. Y-27632 preconditioning also increased the expression of N-cadherin and integrin ß1, the two cell junction proteins. hiPSC-CM+RI were also larger in cell area with greater cytoskeletal alignment and a more rod-like shape than hiPSC-CM-RI, both after transplantation (in vivo) and in culture. The effects of Y-27632 preconditioning on contractile activity and morphology of hiPSC-CMs in culture, as well as on their engraftment rate and apoptotic death in MI mouse grafts, could be recapitulated by hiPSC-CM treatment with the L-type calcium-channel blocker verapamil.

Conclusion

Preconditioning with the ROCK inhibitor Y-27632 increased the engraftment of transplanted hiPSC-CM in a murine MI model, while reversibly impairing hiPSC-CM contractility and promoting adhesion.

A short history of the development of mathematical models of cardiac mechanics

Cardiac mechanics plays a crucial role in atrial and ventricular function, in the regulation of growth and remodelling, in the progression of disease, and the response to treatment. The spatial scale of the critical mechanisms ranges from nm (molecules) to cm (hearts) with the fastest events occurring in milliseconds (molecular events) and the slowest requiring months (growth and remodelling). Due to its complexity and importance, cardiac mechanics has been studied extensively both experimentally and through mathematical models and simulation. Models of cardiac mechanics evolved from seminal studies in skeletal muscle, and developed into cardiac specific, species specific, human specific and finally patient specific calculations. These models provide a formal framework to link multiple experimental assays recorded over nearly 100 years into a single unified representation of cardiac function. This review first provides a summary of the proteins, physiology and anatomy involved in the generation of cardiac pump function. We then describe the evolution of models of cardiac mechanics starting with the early theoretical frameworks describing the link between sarcomeres and muscle contraction, transitioning through myosin-level models to calcium-driven systems, and ending with whole heart patient-specific models.

Hereditary heart disease: pathophysiology, clinical presentation, and animal models of HCM, RCM, and DCM associated with mutations in cardiac myosin light chains

Genetic cardiomyopathies, a group of cardiovascular disorders based on ventricular morphology and function, are among the leading causes of morbidity and mortality worldwide. Such genetically driven forms of hypertrophic (HCM), dilated (DCM), and restrictive (RCM) cardiomyopathies are chronic, debilitating diseases that result from biomechanical defects in cardiac muscle contraction and frequently progress to heart failure (HF). Locus and allelic heterogeneity, as well as clinical variability combined with genetic and phenotypic overlap between different cardiomyopathies, have challenged proper clinical prognosis and provided an incentive for identification of pathogenic variants. This review attempts to provide an overview of inherited cardiomyopathies with a focus on their genetic etiology in myosin regulatory (RLC) and essential (ELC) light chains, which are EF-hand protein family members with important structural and regulatory roles. From the clinical discovery of cardiomyopathy-linked light chain mutations in patients to an array of exploratory studies in animals, and reconstituted and recombinant systems, we have summarized the current state of knowledge on light chain mutations and how they induce physiological disease states via biochemical and biomechanical alterations at the molecular, tissue, and organ levels. Cardiac myosin RLC phosphorylation and the N-terminus ELC have been discussed as two important emerging modalities with important implications in the regulation of myosin motor function, and thus cardiac performance. A comprehensive understanding of such triggers is absolutely necessary for the development of target-specific rescue strategies to ameliorate or reverse the effects of myosin light chain-related inherited cardiomyopathies.

Genetically Encoded Fluorescent Biosensors Illuminate the Spatiotemporal Regulation of Signaling Networks

Cellular signaling networks are the foundation which determines the fate and function of cells as they respond to various cues and stimuli. The discovery of fluorescent proteins over 25 years ago enabled the development of a diverse array of genetically encodable fluorescent biosensors that are capable of measuring the spatiotemporal dynamics of signal transduction pathways in live cells. In an effort to encapsulate the breadth over which fluorescent biosensors have expanded, we endeavored to assemble a comprehensive list of published engineered biosensors, and we discuss many of the molecular designs utilized in their development. Then, we review how the high temporal and spatial resolution afforded by fluorescent biosensors has aided our understanding of the spatiotemporal regulation of signaling networks at the cellular and subcellular level. Finally, we highlight some emerging areas of research in both biosensor design and applications that are on the forefront of biosensor development.

Loss of S100A1 expression leads to Ca2+ release potentiation in mutant mice with disrupted CaM and S100A1 binding to CaMBD2 of RyR1

Calmodulin (CaM) and S100A1 fine-tune skeletal muscle Ca2+ release via opposite modulation of the ryanodine receptor type 1 (RyR1). Binding to and modulation of RyR1 by CaM and S100A1 occurs predominantly at the region ranging from amino acid residue 3614-3640 of RyR1 (here referred to as CaMBD2). Using synthetic peptides, it has been shown that CaM binds to two additional regions within the RyR1, specifically residues 1975-1999 and 4295-4325 (CaMBD1 and CaMBD3, respectively). Because S100A1 typically binds to similar motifs as CaM, we hypothesized that S100A1 could also bind to CaMBD1 and CaMBD3. Our goals were: (1) to establish whether S100A1 binds to synthetic peptides containing CaMBD1 and CaMBD3 using isothermal calorimetry (ITC), and (2) to identify whether S100A1 and CaM modulate RyR1 Ca2+ release activation via sites other than CaMBD2 in RyR1 in its native cellular context. We developed the mouse model (RyR1D-S100A1KO), which expresses point mutation RyR1-L3625D (RyR1D) that disrupts the modulation of RyR1 by CaM and S100A1 at CaMBD2 and also lacks S100A1 (S100A1KO). ITC assays revealed that S100A1 binds with different affinities to CaMBD1 and CaMBD3. Using high-speed Ca2+ imaging and a model for Ca2+ binding and transport, we show that the RyR1D-S100A1KO muscle fibers exhibit a modest but significant increase in myoplasmic Ca2+ transients and enhanced Ca2+ release flux following field stimulation when compared to fibers from RyR1D mice, which were used as controls to eliminate any effect of binding at CaMBD2, but with preserved S100A1 expression. Our results suggest that S100A1, similar to CaM, binds to CaMBD1 and CaMBD3 within the RyR1, but that CaMBD2 appears to be the primary site of RyR1 regulation by CaM and S100A1.

Voltage modulates halothane-triggered Ca2+ release in malignant hyperthermia-susceptible muscle

Malignant hyperthermia can result from mutations in the ryanodine receptor that favor anesthetic-induced Ca²⁺ release. Zullo et al. find that membrane potential modulates the effect of the volatile anesthetic halothane on skeletal muscle ryanodine receptors possessing the Y524S mutation.

Malignant hyperthermia (MH) is a fatal hypermetabolic state that may occur during general anesthesia in susceptible individuals. It is often caused by mutations in the ryanodine receptor RyR1 that favor drug-induced release of Ca²⁺ from the sarcoplasmic reticulum. Here, knowing that membrane depolarization triggers Ca²⁺ release in normal muscle function, we study the cross-influence of membrane potential and anesthetic drugs on Ca²⁺ release. We used short single muscle fibers of knock-in mice heterozygous for the RyR1 mutation Y524S combined with microfluorimetry to measure intracellular Ca²⁺ signals. Halothane, a volatile anesthetic used in contracture testing for MH susceptibility, was equilibrated with the solution superfusing the cells by means of a vaporizer system. In the range 0.2 to 3%, the drug causes significantly larger elevations of free myoplasmic [Ca²⁺] in mutant (YS) compared with wild-type (WT) fibers. Action potential–induced Ca²⁺ signals exhibit a slowing of their time course of relaxation that can be attributed to a component of delayed Ca²⁺ release turnoff. In further experiments, we applied halothane to single fibers that were voltage-clamped using two intracellular microelectrodes and studied the effect of small (10-mV) deviations from the holding potential (−80 mV). Untreated WT fibers show essentially no changes in [Ca²⁺], whereas the Ca²⁺ level of YS fibers increases and decreases on depolarization and hyperpolarization, respectively. The drug causes a significant enhancement of this response. Depolarizing pulses reveal a substantial negative shift in the voltage dependence of activation of Ca²⁺ release. This behavior likely results from the allosteric coupling between RyR1 and its transverse tubular voltage sensor. We conclude that the binding of halothane to RyR1 alters the voltage dependence of Ca²⁺ release in MH-susceptible muscle fibers such that the resting membrane potential becomes a decisive factor for the efficiency of the drug to trigger Ca²⁺ release.

Doublet stimulation increases Ca2+ binding to troponin C to ensure rapid force development in skeletal muscle

Fast-twitch skeletal muscle fibers are often exposed to motor neuron double discharges (≥200 Hz), which markedly increase both the rate of contraction and the magnitude of the resulting force responses. However, the mechanism responsible for these effects is poorly understood, likely because of technical limitations in previous studies. In this study, we measured cytosolic Ca2+ during doublet activation using the low-affinity indicator Mag-Fluo-4 at high temporal resolution and modeled the effects of doublet stimulation on sarcoplasmic reticulum (SR) Ca2+ release, binding of Ca2+ to cytosolic buffers, and force enhancement in fast-twitch fibers. Single isolated fibers respond to doublet pulses with two clear Ca2+ spikes, at doublet frequencies up to 1 KHz. A 200-Hz doublet at the start of a tetanic stimulation train (70 Hz) decreases the drop in free Ca2+ between the first three Ca2+ spikes of the transient, maintaining a higher overall free Ca2+ level during first 20-30 ms of the response. Doublet stimulation also increased the rate of force development in isolated fast-twitch muscles. We also modeled SR Ca2+ release rates during doublet stimulation and showed that Ca2+-dependent inactivation of ryanodine receptor activity is rapid, occurring ≤1ms after initial release. Furthermore, we modeled Ca2+ binding to the main intracellular Ca2+ buffers of troponin C (TnC), parvalbumin, and the SR Ca2+ pump during Ca2+ release and found that the main effect of the second response in the doublet is to more rapidly increase the occupation of the second Ca2+-binding site on TnC (TnC2), resulting in earlier activation of force. We conclude that doublet stimulation maintains high cytosolic Ca2+ levels for longer in the early phase of the Ca2+ response, resulting in faster saturation of TnC2 with Ca2+, faster initiation of cross-bridge cycling, and more rapid force development.

Energetics of muscle contraction: further trials

Knowledge accumulated in the field of energetics of muscle contraction has been reviewed in this article. Active muscle converts chemical energy into heat and work. Therefore, measurements of heat production and mechanical work provide the framework for understanding the process of energy conversion in contraction. In the 1970s, precise comparison between energy output and the associated chemical reactions was performed. It has been found that the two do not match in several situations, resulting in an energy balance discrepancy. More recently, efforts in resolving these discrepancies in the energy balance have been made involving chemical analysis, phosphorus nuclear magnetic resonance spectroscopy, and microcalorimetry. Through reviewing the evidence from these studies, the energy balance discrepancy developed early during isometric contraction has become well understood on a quantitative basis. In this situation energy balance is established when we take into account the binding of Ca to sarcoplasmic proteins such as troponin and parvalbumin, and also the shift of cross-bridge states. On the other hand, the energy balance discrepancy observed during rapid shortening still remains to be clarified. The problem may be related to the essential mechanism of cross-bridge action.

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.

Modeling Ca2+-Bound Troponin in Excitation Contraction Coupling

To explain disparate decay rates of cytosolic Ca²⁺ and structural changes in the thin filaments during a twitch, we model the time course of Ca²⁺-bound troponin (Tn) resulting from the free Ca²⁺ transient of fast skeletal muscle. In fibers stretched beyond overlap, the decay of Ca²⁺ as measured by a change in fluo-3 fluorescence is significantly slower than the intensity decay of the meridional 1/38.5 nm⁻¹ reflection of Tn; this is not simply explained by considering only the Ca²⁺ binding properties of Tn alone (Matsuo et al., 2010). We apply a comprehensive model that includes the known Ca²⁺ binding properties of Tn in the context of the thin filament with and without cycling crossbridges. Calculations based on the model predict that the transient of Ca²⁺-bound Tn correlates with either the fluo-3 time course in muscle with overlapping thin and thick filaments or the intensity of the meridional 1/38.5 nm⁻¹ reflection in overstretched muscle. Hence, cycling crossbridges delay the dissociation of Ca²⁺ from Tn. Correlation with the fluo-3 fluorescence change is not causal given that the transient of Ca²⁺-bound Tn depends on sarcomere length, whereas the fluo-3 fluorescence change does not. Transient positions of tropomyosin calculated from the time course of Ca²⁺-bound Tn are in reasonable agreement with the transient of measured perturbations of the Tn repeat in overlap and non-overlap muscle preparations.

Assessment of Myofilament Ca2+ Sensitivity Underlying Cardiac Excitation-contraction Coupling

Heart failure and cardiac arrhythmias are the leading causes of mortality and morbidity worldwide. However, the mechanism of pathogenesis and myocardial malfunction in the diseased heart remains to be fully clarified. Recent compelling evidence demonstrates that changes in the myofilament Ca(2+) sensitivity affect intracellular Ca(2+) homeostasis and ion channel activities in cardiac myocytes, the essential mechanisms responsible for the cardiac action potential and contraction in healthy and diseased hearts. Indeed, activities of ion channels and transporters underlying cardiac action potentials (e.g., Na(+), Ca(2+) and K(+) channels and the Na(+)-Ca(2+) exchanger) and intracellular Ca(2+) handling proteins (e.g., ryanodine receptors and Ca(2+)-ATPase in sarcoplasmic reticulum (SERCA2a) or phospholamban and its phosphorylation) are conventionally measured to evaluate the fundamental mechanisms of cardiac excitation-contraction (E-C) coupling. Both electrical activities in the membrane and intracellular Ca(2+) changes are the trigger signals of E-C coupling, whereas myofilament is the functional unit of contraction and relaxation, and myofilament Ca(2+) sensitivity is imperative in the implementation of myofibril performance. Nevertheless, few studies incorporate myofilament Ca(2+) sensitivity into the functional analysis of the myocardium unless it is the focus of the study. Here, we describe a protocol that measures sarcomere shortening/re-lengthening and the intracellular Ca(2+) level using Fura-2 AM (ratiometric detection) and evaluate the changes of myofilament Ca(2+) sensitivity in cardiac myocytes from rat hearts. The main aim is to emphasize that myofilament Ca(2+) sensitivity should be taken into consideration in E-C coupling for mechanistic analysis. Comprehensive investigation of ion channels, ion transporters, intracellular Ca(2+) handling, and myofilament Ca(2+) sensitivity that underlie myocyte contractility in healthy and diseased hearts will provide valuable information for designing more effective strategies of translational and therapeutic value.

Effect of 23-day muscle disuse on sarcoplasmic reticulum Ca2+ properties and contractility in human type I and type II skeletal muscle fibers

Inactivity negatively impacts on skeletal muscle function mainly through muscle atrophy. However, recent evidence suggests that the quality of individual muscle fibers is also altered. This study examined the effects of 23 days of unilateral lower limb suspension (ULLS) on specific force and sarcoplasmic reticulum (SR) Ca2+ content in individual skinned muscle fibers. Muscle biopsies of the vastus lateralis were taken from six young healthy adults prior to and following ULLS. After disuse, the endogenous SR Ca2+ content was ∼8% lower in type I fibers and maximal SR Ca2+ capacity was lower in both type I and type II fibers (−11 and −5%, respectively). The specific force, measured in single skinned fibers from three subjects, decreased significantly after ULLS in type II fibers (−23%) but not in type I fibers (−9%). Western blot analyses showed no significant change in the amounts of myosin heavy chain (MHC) I and MHC IIa following the disuse, whereas the amounts of sarco(endo)plasmic reticulum Ca2+-ATPase 1 (SERCA1) and calsequestrin increased by ∼120 and ∼20%, respectively, and the amount of troponin I decreased by ∼21%. These findings suggest that the decline in force and power occurring with muscle disuse is likely to be exacerbated in part by reductions in maximum specific force in type II fibers, and in the amount of releasable SR Ca2+ in both fiber types, the latter not being attributable to a reduced calsequestrin level. Furthermore, the ∼3-wk disuse in human elicits change in SR properties, in particular a more than twofold upregulation in SERCA1 density, before any fiber-type shift.

Energetics of contraction

Muscles convert energy from ATP into useful work, which can be used to move limbs and to transport ions across membranes. The energy not converted into work appears as heat. At the start of contraction heat is also produced when Ca(2+) binds to troponin-C and to parvalbumin. Muscles use ATP throughout an isometric contraction at a rate that depends on duration of stimulation, muscle type, temperature and muscle length. Between 30% and 40% of the ATP used during isometric contraction fuels the pumping Ca(2+) and Na(+) out of the myoplasm. When shortening, muscles produce less force than in an isometric contraction but use ATP at a higher rate and when lengthening force output is higher than the isometric force but rate of ATP splitting is lower. Efficiency quantifies the fraction of the energy provided by ATP that is converted into external work. Each ATP molecule provides 100 zJ of energy that can potentially be converted into work. The mechanics of the myosin cross-bridge are such that at most 50 zJ of work can be done in one ATP consuming cycle; that is, the maximum efficiency of a cross-bridge is ∼50%. Cross-bridges in tortoise muscle approach this limit, producing over 90% of the possible work per cycle. Other muscles are less efficient but contract more rapidly and produce more power.

New method for determining total calcium content in tissue applied to skeletal muscle with and without calsequestrin

The concentration of total calcium in a skeletal muscle appears to be correlated with the muscle’s likely force requirements given by the ratio of body weight to muscle weight.

We describe a new method for determining the concentration of total Ca in whole skeletal muscle samples ([Ca_T]_WM in units of mmoles/kg wet weight) using the Ca-dependent UV absorbance spectra of the Ca chelator BAPTA (1,2-bis(2-aminophenoxy)ethane-N,N,N′,N′-tetraacetic acid). Muscle tissue was homogenized in a solution containing 0.15 mM BAPTA and 0.5% sodium dodecyl sulfate (to permeabilize membranes and denature proteins) and then centrifuged. The solution volume was adjusted so that BAPTA captured essentially all of the Ca. [Ca_T]_WM was obtained with Beer’s law from the absorbance change produced by adding 1 mM EGTA to capture Ca from BAPTA. Results from mouse, rat, and frog muscles were reasonably consistent with results obtained using other methods for estimating total [Ca] in whole muscles and in single muscle fibers. Results with external Ca removed before determining [Ca_T]_WM indicate that most of the Ca was intracellular, indicative of a lack of bound Ca in the extracellular space. In both fast-twitch (extensor digitorum longus, EDL) and slow-twitch (soleus) muscles from mice, [Ca_T]_WM increased approximately linearly with decreasing muscle weight, increasing approximately twofold with a twofold decrease in muscle weight. This suggests that the Ca concentration of smaller muscles might be increased relative to that in larger muscles, thereby increasing the specific force to compensate for the smaller mass. Knocking out the high capacity Ca-binding protein calsequestrin (CSQ) did not significantly reduce [Ca_T]_WM in mouse EDL or soleus muscle. However, in EDL muscles lacking CSQ, muscle weights were significantly lower than in wild-type (WT) muscles and the values of [Ca_T]_WM were, on average, about half the expected WT values, taking into account the above [Ca_T]_WM versus muscle weight relationship. Because greater reductions in [Ca_T]_WM would be predicted in both muscle types, we hypothesize that there is a substantial increase in Ca bound to other sites in the CSQ knockout muscles.

Altered Ca(2+) signaling in skeletal muscle fibers of the R6/2 mouse, a model of Huntington's disease

Huntington's disease (HD) is caused by an expanded CAG trinucleotide repeat within the gene encoding the protein huntingtin. The resulting elongated glutamine (poly-Q) sequence of mutant huntingtin (mhtt) affects both central neurons and skeletal muscle. Recent reports suggest that ryanodine receptor-based Ca(2+) signaling, which is crucial for skeletal muscle excitation-contraction coupling (ECC), is changed by mhtt in HD neurons. Consequently, we searched for alterations of ECC in muscle fibers of the R6/2 mouse, a mouse model of HD. We performed fluorometric recordings of action potentials (APs) and cellular Ca(2+) transients on intact isolated toe muscle fibers (musculi interossei), and measured L-type Ca(2+) inward currents on internally dialyzed fibers under voltage-clamp conditions. Both APs and AP-triggered Ca(2+) transients showed slower kinetics in R6/2 fibers than in fibers from wild-type mice. Ca(2+) removal from the myoplasm and Ca(2+) release flux from the sarcoplasmic reticulum were characterized using a Ca(2+) binding and transport model, which indicated a significant reduction in slow Ca(2+) removal activity and Ca(2+) release flux both after APs and under voltage-clamp conditions. In addition, the voltage-clamp experiments showed a highly significant decrease in L-type Ca(2+) channel conductance. These results indicate profound changes of Ca(2+) turnover in skeletal muscle of R6/2 mice and suggest that these changes may be associated with muscle pathology in HD.

Contractile properties and sarcoplasmic reticulum calcium content in type I and type II skeletal muscle fibres in active aged humans

KEY POINTS

Muscle weakness in old age is due in large part to an overall loss of skeletal muscle tissue, but it remains uncertain how much also stems from alterations in the properties of the individual muscle fibres. This study examined the contractile properties and amount of stored intracellular calcium in single muscle fibres of Old (70 ± 4 years) and Young (22 ± 3 years) adults. The maximum level of force production (per unit cross-sectional area) in fast twitch fibres in Old subjects was lower than in Young subjects, and the fibres were also less sensitive to activation by calcium. The amount of calcium stored inside muscle fibres and available to trigger contraction was also lower in both fast- and slow-twitch muscle fibres in the Old subjects. These findings indicate that muscle weakness in old age stems in part from an impaired capacity for force production in the individual muscle fibres.

ABSTRACT

This study examined the contractile properties and sarcoplasmic reticulum (SR) Ca(2+) content in mechanically skinned vastus lateralis muscle fibres of Old (70 ± 4 years) and Young (22 ± 3 years) humans to investigate whether changes in muscle fibre properties contribute to muscle weakness in old age. In type II fibres of Old subjects, specific force was reduced by ∼17% and Ca(2+) sensitivity was also reduced (pCa50 decreased ∼0.05 pCa units) relative to that in Young. S-Glutathionylation of fast troponin I (TnIf ) markedly increased Ca(2+) sensitivity in type II fibres, but the increase was significantly smaller in Old versus Young (+0.136 and +0.164 pCa unit increases, respectively). Endogenous and maximal SR Ca(2+) content were significantly smaller in both type I and type II fibres in Old subjects. In fibres of Young, the SR could be nearly fully depleted of Ca(2+) by a combined caffeine and low Mg(2+) stimulus, whereas in fibres of Old the amount of non-releasable Ca(2+) was significantly increased (by > 12% of endogenous Ca(2+) content). Western blotting showed an increased proportion of type I fibres in Old subjects, and increased amounts of calsequestrin-2 and calsequestrin-like protein. The findings suggest that muscle weakness in old age is probably attributable in part to (i) an increased proportion of type I fibres, (ii) a reduction in both maximum specific force and Ca(2+) sensitivity in type II fibres, and also a decreased ability of S-glutathionylation of TnIf to counter the fatiguing effects of metabolites on Ca(2+) sensitivity, and (iii) a reduction in the amount of releasable SR Ca(2+) in both fibre types.

Balanced changes in Ca buffering by SERCA and troponin contribute to Ca handling during β-adrenergic stimulation in cardiac myocytes

AIMS

During activation of cardiac myocytes, less than 1% of cytosolic Ca is free; the rest is bound to buffers, largely SERCA, and troponin C. Signalling by phosphorylation, as occurs during β-adrenergic stimulation, changes the Ca-binding affinity of these proteins and may affect the systolic Ca transient. Our aim was to determine the effects of β-adrenergic stimulation on Ca buffering and to differentiate between the roles of SERCA and troponin.

METHODS AND RESULTS

Ca buffering was studied in cardiac myocytes from mice: wild-type (WT), phospholamban-knockout (PLN-KO), and mice expressing slow skeletal troponin I (ssTnI) that is not protein kinase A phosphorylatable. WT cells showed no change in Ca buffering in response to the β-adrenoceptor agonist isoproterenol (ISO). However, ISO decreased Ca buffering in PLN-KO myocytes, presumably unmasking the role of troponin. This effect was confirmed in WT cells in which SERCA activity was blocked with the application of thapsigargin. In contrast, ISO increased Ca buffering in ssTnI cells, presumably revealing the effect of an increase in Ca binding to SERCA.

CONCLUSIONS

These data indicate the individual roles played by SERCA and troponin in Ca buffering during β-adrenergic stimulation and that these two buffers effectively counterbalance each other so that Ca buffering remains constant during β-adrenergic stimulation, a factor which may be physiologically important. This study also emphasizes the importance of taking into account Ca buffering, particularly in disease states where Ca binding to myofilaments or SERCA may be altered.

In vitro rescue study of a malignant familial hypertrophic cardiomyopathy phenotype by pseudo-phosphorylation of myosin regulatory light chain

Pseudo-phosphorylation of cardiac myosin regulatory light chain (RLC) has never been examined as a rescue method to alleviate a cardiomyopathy phenotype brought about by a disease causing mutation in the myosin RLC. This study focuses on the aspartic acid to valine substitution (D166V) in the myosin RLC shown to be associated with a malignant phenotype of familial hypertrophic cardiomyopathy (FHC). The mutation has also been demonstrated to cause severe functional abnormalities in transgenic mice expressing D166V in the heart. To explore this novel rescue strategy, pseudo-phosphorylation of D166V was used to determine whether the D166V-induced detrimental phenotype could be brought back to the level of wild-type (WT) RLC. The S15D substitution at the phosphorylation site of RLC was inserted into the recombinant WT and D166V mutant to mimic constitutively phosphorylated RLC proteins. Non-phosphorylatable (S15A) constructs were used as controls. A multi-faceted approach was taken to determine the effect of pseudo-phosphorylation on the ability of myosin to generate force and motion. Using mutant reconstituted porcine cardiac muscle preparations, we showed an S15D-induced rescue of both the enzymatic and binding properties of D166V-myosin to actin. A significant increase in force production capacity was noted in the in vitro motility assays for S15D-D166V vs. D166V reconstituted myosin. A similar pseudo-phosphorylation induced effect was observed on the D166V-elicited abnormal Ca(2+) sensitivity of force in porcine papillary muscle strips reconstituted with phosphomimic recombinant RLCs. Results from this study demonstrate a novel in vitro rescue strategy that could be utilized in vivo to ameliorate a malignant cardiomyopathic phenotype. We show for the first time that pseudo-RLC phosphorylation can reverse the majority of the mutation-induced phenotypes highlighting the importance of RLC phosphorylation in combating cardiac disease.

Gene regulation, alternative splicing, and posttranslational modification of troponin subunits in cardiac development and adaptation: a focused review

Troponin plays a central role in regulating the contraction and relaxation of vertebrate striated muscles. This review focuses on the isoform gene regulation, alternative RNA splicing, and posttranslational modifications of troponin subunits in cardiac development and adaptation. Transcriptional and posttranscriptional regulations such as phosphorylation and proteolysis modifications, and structure-function relationships of troponin subunit proteins are summarized. The physiological and pathophysiological significances are discussed for impacts on cardiac muscle contractility, heart function, and adaptations in health and diseases.

Molecular dynamics simulations of the cardiac troponin complex performed with FRET distances as restraints

Cardiac troponin (cTn) is the Ca²⁺-sensitive molecular switch that controls cardiac muscle activation and relaxation. However, the molecular detail of the switching mechanism and how the Ca²⁺ signal received at cardiac troponin C (cTnC) is communicated to cardiac troponin I (cTnI) are still elusive. To unravel the structural details of troponin switching, we performed ensemble Förster resonance energy transfer (FRET) measurements and molecular dynamic (MD) simulations of the cardiac troponin core domain complex. The distance distributions of forty five inter-residue pairs were obtained under Ca²⁺-free and saturating Ca²⁺ conditions from time-resolved FRET measurements. These distances were incorporated as restraints during the MD simulations of the cardiac troponin core domain. Compared to the Ca²⁺-saturated structure, the absence of regulatory Ca²⁺ perturbed the cTnC N-domain hydrophobic pocket which assumed a closed conformation. This event partially unfolded the cTnI regulatory region/switch. The absence of Ca²⁺, induced flexibility to the D/E linker and the cTnI inhibitory region, and rotated the cTnC N-domain with respect to rest of the troponin core domain. In the presence of saturating Ca²⁺ the above said phenomenon were absent. We postulate that the secondary structure perturbations experienced by the cTnI regulatory region held within the cTnC N-domain hydrophobic pocket, coupled with the rotation of the cTnC N-domain would control the cTnI mobile domain interaction with actin. Concomitantly the rotation of the cTnC N-domain and perturbation of the D/E linker rigidity would control the cTnI inhibitory region interaction with actin to effect muscle relaxation.

Endogenous and maximal sarcoplasmic reticulum calcium content and calsequestrin expression in type I and type II human skeletal muscle fibres

The relationship between sarcoplasmic reticulum (SR) Ca(2+) content and calsequestrin (CSQ) isoforms was investigated in human skeletal muscle. A fibre-lysing assay was used to quantify the endogenous Ca(2+) content and maximal Ca(2+) capacity of the SR in skinned segments of type I and type II fibres from vastus lateralis muscles of young healthy adults. Western blotting of individual fibres showed the great majority contained either all fast or all slow isoforms of myosin heavy chain (MHC), troponins C and I, tropomyosin and SERCA, and that the strontium sensitivity of the force response was closely indicative of the troponin C isoform present. The endogenous SR Ca(2+) content was slightly lower in type I compared to type II fibres (0.76 ± 0.03 and 0.85 ± 0.02 mmol Ca(2+) per litre of fibre, respectively), with virtually all of this Ca(2+) evidently being in the SR, as it could be rapidly released with a caffeine-low [Mg(2+)] solution (only 0.08 ± 0.01 and <0.07 mmol l(-1), respectively, remaining). The maximal Ca(2+) content that could be reached with SR Ca(2+) loading was 1.45 ± 0.04 and 1.79 ± 0.03 mmol l(-1) in type I and type II fibres, respectively (P < 0.05). In non-lysed skinned fibres, where the SR remained functional, repeated cycles of caffeine-induced Ca(2+) release and subsequent Ca(2+) reloading similarly indicated that (i) maximal SR Ca(2+) content was lower in type I fibres than in type II fibres (P < 0.05), and (ii) the endogenous Ca(2+) content represented a greater percentage of maximal content in type I fibres compared to type II fibres (∼59% and 41%, respectively, P < 0.05). Type II fibres were found on average to contain ∼3-fold more CSQ1 and ∼5-fold less CSQ2 than type I fibres (P < 0.001). The findings are consistent with the SR Ca(2+) content characteristics in human type II fibres being primarily determined by the CSQ1 abundance, and in type I fibres by the combined amounts of both CSQ1 and CSQ2.

Evoked centripetal Ca(2+) mobilization in cardiac Purkinje cells: insight from a model of three Ca(2+) release regions

Despite strong suspicion that abnormal Ca(2+) handling in Purkinje cells (P-cells) is implicated in life-threatening forms of ventricular tachycardias, the mechanism underlying the Ca(2+) cycling of these cells under normal conditions is still unclear. There is mounting evidence that P-cells have a unique Ca(2+) handling system. Notably complex spontaneous Ca(2+) activity was previously recorded in canine P-cells and was explained by a mechanistic hypothesis involving a triple layered system of Ca(2+) release channels. Here we examined the validity of this hypothesis for the electrically evoked Ca(2+) transient which was shown, in the dog and rabbit, to occur progressively from the periphery to the interior of the cell. To do so, the hypothesis was incorporated in a model of intracellular Ca(2+) dynamics which was then used to reproduce numerically the Ca(2+) activity of P-cells under stimulated conditions. The modelling was thus performed through a 2D computational array that encompassed three distinct Ca(2+) release nodes arranged, respectively, into three consecutive adjacent regions. A system of partial differential equations (PDEs) expressed numerically the principal cellular functions that modulate the local cytosolic Ca(2+) concentration (Cai). The apparent node-to-node progression of elevated Cai was obtained by combining Ca(2+) diffusion and 'Ca(2+)-induced Ca(2+) release'. To provide the modelling with a reliable experimental reference, we first re-examined the Ca(2+) mobilization in swine stimulated P-cells by 2D confocal microscopy. As reported earlier for the dog and rabbit, a centripetal Ca(2+) transient was readily visible in 22 stimulated P-cells from six adult Yucatan swine hearts (pacing rate: 0.1 Hz; pulse duration: 25 ms, pulse amplitude: 10% above threshold; 1 mm Ca(2+); 35°C; pH 7.3). An accurate replication of the observed centripetal Ca(2+) propagation was generated by the model for four representative cell examples and confirmed by statistical comparisons of simulations against cell data. Selective inactivation of Ca(2+) release regions of the computational array showed that an intermediate layer of Ca(2+) release nodes with an ~30-40% lower Ca(2+) activation threshold was required to reproduce the phenomenon. Our computational analysis was therefore fully consistent with the activation of a triple layered system of Ca(2+) release channels as a mechanism of centripetal Ca(2+) signalling in P-cells. Moreover, the model clearly indicated that the intermediate Ca(2+) release layer with increased sensitivity for Ca(2+) plays an important role in the specific intracellular Ca(2+) mobilization of Purkinje fibres and could therefore be a relevant determinant of cardiac conduction.

Comparison of myoplasmic calcium movements during excitation-contraction coupling in frog twitch and mouse fast-twitch muscle fibers

Single twitch fibers from frog leg muscles were isolated by dissection and micro-injected with furaptra, a rapidly responding fluorescent Ca(2+) indicator. Indicator resting fluorescence (FR) and the change evoked by an action potential (ΔF) were measured at long sarcomere length (16°C); ΔF/FR was scaled to units of ΔfCaD, the change in fraction of the indicator in the Ca(2+)-bound form. ΔfCaD was simulated with a multicompartment model of the underlying myoplasmic Ca(2+) movements, and the results were compared with previous measurements and analyses in mouse fast-twitch fibers. In frog fibers, sarcoplasmic reticulum (SR) Ca(2+) release evoked by an action potential appears to be the sum of two components. The time course of the first component is similar to that of the entire Ca(2+) release waveform in mouse fibers, whereas that of the second component is severalfold slower; the fractional release amounts are ~0.8 (first component) and ~0.2 (second component). Similar results were obtained in frog simulations with a modified model that permitted competition between Mg(2+) and Ca(2+) for occupancy of the regulatory sites on troponin. An anatomical basis for two release components in frog fibers is the presence of both junctional and parajunctional SR Ca(2+) release channels (ryanodine receptors [RyRs]), whereas mouse fibers (usually) have only junctional RyRs. Also, frog fibers have two RyR isoforms, RyRα and RyRβ, whereas the mouse fibers (usually) have only one, RyR1. Our simulations suggest that the second release component in frog fibers functions to supply extra Ca(2+) to activate troponin, which, in mouse fibers, is not needed because of the more favorable location of their triadic junctions (near the middle of the thin filament). We speculate that, in general, parajunctional RyRs permit increased myofilament activation in fibers whose triadic junctions are located at the z-line.

Regulation of ion gradients across myocardial ischemic border zones: a biophysical modelling analysis

The myocardial ischemic border zone is associated with the initiation and sustenance of arrhythmias. The profile of ionic concentrations across the border zone play a significant role in determining cellular electrophysiology and conductivity, yet their spatial-temporal evolution and regulation are not well understood. To investigate the changes in ion concentrations that regulate cellular electrophysiology, a mathematical model of ion movement in the intra and extracellular space in the presence of ionic, potential and material property heterogeneities was developed. The model simulates the spatial and temporal evolution of concentrations of potassium, sodium, chloride, calcium, hydrogen and bicarbonate ions and carbon dioxide across an ischemic border zone. Ischemia was simulated by sodium-potassium pump inhibition, potassium channel activation and respiratory and metabolic acidosis. The model predicted significant disparities in the width of the border zone for each ionic species, with intracellular sodium and extracellular potassium having discordant gradients, facilitating multiple gradients in cellular properties across the border zone. Extracellular potassium was found to have the largest border zone and this was attributed to the voltage dependence of the potassium channels. The model also predicted the efflux of [Formula: see text] from the ischemic region due to electrogenic drift and diffusion within the intra and extracellular space, respectively, which contributed to [Formula: see text] depletion in the ischemic region.

The green tea polyphenol (-)-epigallocatechin-3-gallate inhibits magnesium binding to the C-domain of cardiac troponin C

Cardiac muscle contraction is activated via the single Ca(2+)-binding site (site II) in the N-domain of troponin C (cTnC). The two Ca(2+)/Mg(2+) binding sites in the C-domain of cTnC (sites III and IV) have been considered to play a purely structural role in anchoring cTnC to the thin filament. However, several recent discoveries suggest a possible role of this domain in contractile regulation. The green tea polyphenol (-)-epigallocatechin 3-gallate (EGCg), which binds specifically to the C-domain of cTnC, reduces cardiac myofilament Ca(2+) sensitivity along with maximum force and acto-myosin ATPase activity. We have determined the effect of EGCg on Ca(2+) and Mg(2+) binding to the C-domain of cTnC. In the absence of Mg(2+) there was no significant effect of EGCg on the Ca(2+)-cTnC affinity. Surprisingly, in the presence of Mg(2+) EGCg caused an increase in Ca(2+) affinity for sites III and IV of cTnC. However, in the absence of Ca(2+) the addition of EGCg caused a significant reduction in Mg(2+)-cTnC affinity. This reduction is presumably responsible for the increase in Ca(2+)-cTnC affinity produced by EGCg in the presence of Mg(2+). We propose that the inhibitory effect of EGCg on myofilament Ca(2+) activation may be related to an enhanced Ca(2+)-Mg(2+)exchange at sites III and IV of cTnC, which might reduce the myosin crossbridge dependent component of thin filament activation.