β-adrenergic effects on cardiac myofilaments and contraction in an integrated rabbit ventricular myocyte model

Phosphorylation of phospholamban promotes SERCA2a activation by dwarf open reading frame (DWORF)

In cardiac myocytes, the type 2a sarco/endoplasmic reticulum Ca-ATPase (SERCA2a) plays a key role in intracellular Ca regulation. Due to its critical role in heart function, SERCA2a activity is tightly regulated by different mechanisms, including micropeptides. While phospholamban (PLB) is a well-known SERCA2a inhibitor, dwarf open reading frame (DWORF) is a recently identified SERCA2a activator. Since PLB phosphorylation is the most recognized mechanism of SERCA2a activation during adrenergic stress, we studied whether PLB phosphorylation also affects SERCA2a regulation by DWORF. By using confocal Ca imaging in a HEK293 expressing cell system, we analyzed the effect of the co-expression of PLB and DWORF using a bicistronic construct on SERCA2a-mediated Ca uptake. Under these conditions of matched expression of PLB and DWORF, we found that SERCA2a inhibition by non-phosphorylated PLB prevails over DWORF activating effect. However, when PLB is phosphorylated at PKA and CaMKII sites, not only PLB's inhibitory effect was relieved, but SERCA2a was effectively activated by DWORF. Förster resonance energy transfer (FRET) analysis between SERCA2a and DWORF showed that DWORF has a higher relative affinity for SERCA2a when PLB is phosphorylated. Thus, SERCA2a regulation by DWORF responds to the PLB phosphorylation status, suggesting that DWORF might contribute to SERCA2a activation during conditions of adrenergic stress.

Computational Modeling of Cardiac Electrophysiology

Mathematical modeling and simulation are well-established and powerful tools to integrate experimental data of individual components of cardiac electrophysiology, excitation-contraction coupling, and regulatory signaling pathways, to gain quantitative and mechanistic insight into pathophysiological processes and guide therapeutic strategies. Here, we briefly describe the processes governing cardiac myocyte electrophysiology and Ca2+ handling and their regulation, as well as action potential propagation in tissue. We discuss the models and methods used to describe these phenomena, including procedures for model parameterization and validation, in addition to protocols for model interrogation and analysis and techniques that account for phenotypic variability and parameter uncertainty. Our objective is to provide a summary of basic concepts and approaches as a resource for scientists training in this discipline and for all researchers aiming to gain an understanding of cardiac modeling studies.

Functional reserve and contractile phenotype of atrial myocardium from patients with atrial remodeling without and with atrial fibrillation

Atrial contractility and functional reserve in atrial remodeling (AR) without (AR/-AF) or with atrial fibrillation (AR/+AF) are not well characterized. In this study, functional measurements were performed in right atrial muscle strips (n = 71) obtained from patients (N = 22) undergoing routine cardiac surgery with either no AR [left atrial (LA) diameter < 40 mm and no history of AF (hAF)], AR/-AF (LA diameter ≥ 40 mm, no hAF), or AR/+AF (hAF and LA diameter ≥ 40 mm or LAEF < 45%). AR/-AF and AR/+AF were associated with a prolongation of half-time-to-peak (HTTP, P < 0.001) and time-to-peak (TTP) contraction (P < 0.01) when compared with no AR. This effect was seen at baseline and during β-adrenergic stimulation with isoproterenol (Iso). Early relaxation assessed by half-relaxation time (HRT) was prolonged in AR/-AF (P = 0.03) but not in AR/+AF when compared with no AR at baseline, but this delay in relaxation in AR/-AF was attenuated with Iso. Late relaxation (τ) did not differ between AR/-AF and no AR but was consistently shorter in AR/+AF than no AR before (P = 0.04) and during Iso (P = 0.01), indicating accelerated late relaxation in AR/+AF. Relative force increase during Iso was higher (P = 0.01) and more dispersed (P = 0.047) in patients with AR/+AF. Relative adrenergic response was unaltered in the myocardium of patients with AR/-AF and AR/+AF. In conclusion, AR/-AF and AR/+AF are associated with changes in myocardial inotropic reserve and contractility. The changes are particularly pronounced in patients with AR/+AF, suggesting that the progression from AR/-AF to AR/+AF is associated with progressive alterations in atrial function that may contribute to arrhythmogenesis.NEW & NOTEWORTHY Mechanical alterations in atrial remodeling without (AR/-AF) and with atrial fibrillation (AR/+AF) have not been studied in detail in human atrial tissue preparations. To our knowledge, this is the first study to compare the mechanical phenotype and inotropic reserve in human atrial myocardial preparations from patients with no atrial remodeling, AR/-AF, and AR/+AF. We identify specific patterns of contractile dysfunction and heterogeneity for both, AR/-AF and AR/+AF, indicating the progression of atrial disease.

Steady-state and transient effects of SK channel block and adrenergic stimulation to counteract acetylcholine-induced arrhythmogenic effects in the human atria: A computational study

Hyperactivity of the parasympathetic nervous system has been linked to the development of paroxysmal atrial fibrillation (AF). The parasympathetic neurotransmitter acetylcholine (ACh) causes a reduction in action potential (AP) duration (APD) and an increase in resting membrane potential (RMP), both of which contribute to enhance the risk for reentry. Research suggests that small-conductance calcium activated potassium (SK) channels may be an effective target for treating AF. Therapies targeting the autonomic nervous system, either alone or in combination with other drugs, have been explored and have been shown to decrease the incidence of atrial arrhythmias. This study uses computational modeling and simulation to examine the impact of SK channel block (SKb) and β-adrenergic stimulation through Isoproterenol (Iso) on countering the negative effects of cholinergic activity in human atrial cell and 2D tissue models. The steady-state effects of Iso and/or SKb on AP shape, APD at 90% repolarization (APD₉₀) and RMP were evaluated. The ability to terminate stable rotational activity in cholinergically-stimulated 2D tissue models of AF was also investigated. A range of SKb and Iso application kinetics, which reflect varying drug binding rates, were taken into consideration. The results showed that SKb alone prolonged APD₉₀ and was able to stop sustained rotors in the presence of ACh concentrations up to 0.01 μM. Iso terminated rotors under all tested ACh concentrations, but resulted in highly-variable steady-state outcomes depending on baseline AP morphology. Importantly, the combination of SKb and Iso resulted in greater APD₉₀ prolongation and showed promising anti-arrhythmic potential by stopping stable rotors and preventing re-inducibility.

3-Iodothyronamine, a trace amine-associated receptor agonist, regulates intracellular Ca2+ increases via CaMK II through Epac2 in rat cerebral arterioles

Trace amines (TAs) in the nervous system bind to TA-associated receptors (TAARs) and are involved in the regulation of monoaminergic functions. Among TAAR subtypes, TAAR1 has been implicated in the development of neurological disorders, such as schizophrenia. The present study investigated the effects of the TAAR1 agonist, 3-iodothyronamine (T1AM) on cerebral arterioles using fluctuations in the intracellular concentration of Ca2+ ([Ca2+]i) as an index of contractile responses. In cerebral arterioles, most of the TAAR agonists did not increase [Ca2+]i, while only T1AM elevated [Ca2+]i in vascular smooth muscle cells. This increase involved extracellular Ca2+ influx through T-type Ca2+ channels and inositol trisphosphate- and ryanodine-receptor-mediated Ca2+ release from intracellular stores. The inhibition of the cAMP sensor, exchange protein directly activated by cAMP (Epac) 2, and calmodulin kinase (CaMK) II strongly inhibited Ca2+ elevations. The present study revealed that T1AM acted not only on the TAAR1 receptor as previously suggested, but also on other G-protein coupled receptors and/or signal transduction systems to increase intracellular Ca2+ in cerebral arteriole smooth muscle cells. These results suggest that when using T1AM in clinical practice, attention should be paid to the early rise in blood pressure.

Rad regulation of CaV1.2 channels controls cardiac fight-or-flight response

Fight-or-flight responses involve β-adrenergic-induced increases in heart rate and contractile force. In the present study, we uncover the primary mechanism underlying the heart's innate contractile reserve. We show that four protein kinase A (PKA)-phosphorylated residues in Rad, a calcium channel inhibitor, are crucial for controlling basal calcium current and essential for β-adrenergic augmentation of calcium influx in cardiomyocytes. Even with intact PKA signaling to other proteins modulating calcium handling, preventing adrenergic activation of calcium channels in Rad-phosphosite-mutant mice (4SA-Rad) has profound physiological effects: reduced heart rate with increased pauses, reduced basal contractility, near-complete attenuation of β-adrenergic contractile response and diminished exercise capacity. Conversely, expression of mutant calcium-channel β-subunits that cannot bind 4SA-Rad is sufficient to enhance basal calcium influx and contractility to adrenergically augmented levels of wild-type mice, rescuing the failing heart phenotype of 4SA-Rad mice. Hence, disruption of interactions between Rad and calcium channels constitutes the foundation toward next-generation therapeutics specifically enhancing cardiac contractility.

Effects of BNP and Sacubitrilat/Valsartan on Atrial Functional Reserve and Arrhythmogenesis in Human Myocardium

Background

Although the angiotensin receptor-neprilysin inhibitor (ARNI) sacubitril/valsartan started a new era in heart failure (HF) treatment, less is known about the tissue-level effects of the drug on the atrial myocardial functional reserve and arrhythmogenesis.

Methods and Results

Right atrial (RA) biopsies were retrieved from patients (n = 42) undergoing open-heart surgery, and functional experiments were conducted in muscle strips (n = 101). B-type natriuretic peptide (BNP) did not modulate systolic developed force in human myocardium during β-adrenergic stimulation, but it significantly reduced diastolic tension (p &lt; 0.01) and the probability of arrhythmias (p &lt; 0.01). In addition, patient's plasma NTproBNP positively correlated with isoproterenol-induced contractile reserve in atrial tissue in vitro (r = 0.65; p &lt; 0.01). Sacubitrilat+valsartan (Sac/Val) did not show positive inotropic effects on atrial trabeculae function but reduced arrhythmogeneity. Atrial and ventricular biopsies from patients with end-stage HF (n = 10) confirmed that neprilysin (NEP) is equally expressed in human atrial and ventricular myocardium. RA NEP expression correlates positively with RA ejection fraction (EF) (r = 0.806; p &lt; 0.05) and left ventricle (LV) NEP correlates inversely with left atrial (LA) volume (r = −0.691; p &lt; 0.05).

Conclusion

BNP ameliorates diastolic tension during adrenergic stress in human atrial myocardium and may have positive long-term effects on the inotropic reserve. BNP and Sac/Val reduce atrial arrhythmogeneity during adrenergic stress in vitro. Myocardial NEP expression is downregulated with declining myocardial function, suggesting a compensatory mechanism in HF.

Novel Analysis Method for Beating Cells Videomicroscopy Data: Functional Characterization of Culture Samples

Cell culture of cardiac tissue analog is becoming increasingly interesting for regenerative medicine (cell therapy and tissue engineering) and is widely used for high throughput cardiotoxicity. As a cost-effective approach to rapidly discard new compounds with high toxicity risks, cardiotoxicity evaluation is firstly done in vitro requiring cells/tissue with physiological/pathological characteristics (close to in vivo properties). Studying multicellular electrophysiological and contractile properties is needed to assess drug effects. Techniques favoring process automation which could help in simplifying screening drug candidates are thus of central importance. A lot of effort has been made to ameliorate in vitro models including several in vitro platforms for engineering neonatal rat cardiac tissues. However, most of the initial evaluation is done by studying the rate of activity. In this study, we present new approaches that use the videomicroscopy video of monolayer activity to study contractile properties of beating cells in culture. Two new variables are proposed which are linked to the contraction dynamics and are dependent on the rhythm of activity. Methods for evaluation of regional synchronicity within the image field of view are also presented that can rapidly determine regions with abnormal activity or heterogeneity in contraction dynamics.

Adrenergic Regulation of Calcium Channels in the Heart

Each heartbeat is initiated by the action potential, an electrical signal that depolarizes the plasma membrane and activates a cycle of calcium influx via voltage-gated calcium channels, calcium release via ryanodine receptors, and calcium reuptake and efflux via calcium-ATPase pumps and sodium-calcium exchangers. Agonists of the sympathetic nervous system bind to adrenergic receptors in cardiomyocytes, which, via cascading signal transduction pathways and protein kinase A (PKA), increase the heart rate (chronotropy), the strength of myocardial contraction (inotropy), and the rate of myocardial relaxation (lusitropy). These effects correlate with increased intracellular concentration of calcium, which is required for the augmentation of cardiomyocyte contraction. Despite extensive investigations, the molecular mechanisms underlying sympathetic nervous system regulation of calcium influx in cardiomyocytes have remained elusive over the last 40 years. Recent studies have uncovered the mechanisms underlying this fundamental biologic process, namely that PKA phosphorylates a calcium channel inhibitor, Rad, thereby releasing inhibition and increasing calcium influx. Here, we describe an updated model for how signals from adrenergic agonists are transduced to stimulate calcium influx and contractility in the heart.

Quantitative cross-species translators of cardiac myocyte electrophysiology: Model training, experimental validation, and applications

Animal experimentation is key in the evaluation of cardiac efficacy and safety of novel therapeutic compounds. However, interspecies differences in the mechanisms regulating excitation-contraction coupling can limit the translation of experimental findings from animal models to human physiology and undermine the assessment of drugs’ efficacy and safety. Here, we built a suite of translators for quantitatively mapping electrophysiological responses in ventricular myocytes across species. We trained these statistical operators using a broad dataset obtained by simulating populations of our biophysically detailed computational models of action potential and Ca2+ transient in mouse, rabbit, and human. We then tested our translators against experimental data describing the response to stimuli, such as ion channel block, change in beating rate, and β-adrenergic challenge. We demonstrate that this approach is well suited to predicting the effects of perturbations across different species or experimental conditions and suggest its integration into mechanistic studies and drug development pipelines.

A mathematical model of hiPSC cardiomyocytes electromechanics

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

Frank-Starling mechanism, fluid responsiveness, and length-dependent activation: Unravelling the multiscale behaviors with an in silico analysis

The Frank-Starling mechanism is a fundamental regulatory property which underlies the cardiac output adaptation to venous filling. Length-dependent activation is generally assumed to be the cellular origin of this mechanism. At the heart scale, it is commonly admitted that an increase in preload (ventricular filling) leads to an increased cellular force and an increased volume of ejected blood. This explanation also forms the basis for vascular filling therapy. It is actually difficult to unravel the exact nature of the relationship between length-dependent activation and the Frank-Starling mechanism, as three different scales (cellular, ventricular and cardiovascular) are involved. Mathematical models are powerful tools to overcome these limitations. In this study, we use a multiscale model of the cardiovascular system to untangle the three concepts (length-dependent activation, Frank-Starling, and vascular filling). We first show that length-dependent activation is required to observe both the Frank-Starling mechanism and a positive response to high vascular fillings. Our results reveal a dynamical length dependent activation-driven response to changes in preload, which involves interactions between the cellular, ventricular and cardiovascular levels and thus highlights fundamentally multiscale behaviors. We show however that the cellular force increase is not enough to explain the cardiac response to rapid changes in preload. We also show that the absence of fluid responsiveness is not related to a saturating Frank-Starling effect. As it is challenging to study those multiscale phenomena experimentally, this computational approach contributes to a more comprehensive knowledge of the sophisticated length-dependent properties of cardiac muscle.

It is commonly admitted that the length-dependent activation is the cellular property underlying the Frank-Starling mechanism. However, it is challenging to assess the cardiac cell length in vivo, and to link it to system-level variables. In this study, we investigated the relationship between length-dependent activation (LDA), the Frank-Starling mechanism, and vascular filling therapy using a multiscale model of the cardiovascular system. We showed that the Frank-Starling mechanism is a multiscale, dynamical and LDA-driven phenomenon involving essential connections between the three cellular, ventricular and cardiovascular scales. We also demonstrated that the response to vascular filling therapy is similarly multiscale and LDA-driven, and that afterload is an important factor to account for when trying to predict the system fluid responsiveness.

Determinants of Ca2+ release restitution: Insights from genetically altered animals and mathematical modeling

Each heartbeat is followed by a refractory period. Recovery from refractoriness is known as Ca²⁺ release restitution (CRR), and its alterations are potential triggers of Ca²⁺ arrhythmias. Although the control of CRR has been associated with SR Ca²⁺ load and RYR2 Ca²⁺ sensitivity, the relative role of some of the determinants of CRR remains largely undefined. An intriguing point, difficult to dissect and previously neglected, is the possible independent effect of SR Ca²⁺ content versus the velocity of SR Ca²⁺ refilling on CRR. To assess these interrogations, we used isolated myocytes with phospholamban (PLN) ablation (PLNKO), knock-in mice with pseudoconstitutive CaMKII phosphorylation of RYR2 S2814 (S2814D), S2814D crossed with PLNKO mice (SDKO), and a previously validated human cardiac myocyte model. Restitution of cytosolic Ca²⁺ (Fura-2 AM) and L-type calcium current (ICaL; patch-clamp) was evaluated with a two-pulse (S₁/S₂) protocol. CRR and ICaL restitution increased as a function of the (S₂-S₁) coupling interval, following an exponential curve. When SR Ca²⁺ load was increased by increasing extracellular [Ca²⁺] from 2.0 to 4.0 mM, CRR and ICaL restitution were enhanced, suggesting that ICaL restitution may contribute to the faster CRR observed at 4.0 mM [Ca²⁺]. In contrast, ICaL restitution did not differ among the different mouse models. For a given SR Ca²⁺ load, CRR was accelerated in S2814D myocytes versus WT, but not in PLNKO and SDKO myocytes versus WT and S2814D, respectively. The model mimics all experimental data. Moreover, when the PLN ablation-induced decrease in RYR2 expression was corrected, the model revealed that CRR was accelerated in PLNKO and SDKO versus WT and S2814D myocytes, consistent with the enhanced velocity of refilling, SR [Ca²⁺] recovery, and CRR. We speculate that refilling rate might enhance CRR independently of SR Ca²⁺ load.

Computationally efficient model of myocardial electromechanics for multiscale simulations

A model of myocardial electromechanics is suggested. It combines modified and simplified versions of previously published models of cardiac electrophysiology, excitation-contraction coupling, and mechanics. The mechano-calcium and mechano-electrical feedbacks, including the strain-dependence of the propagation velocity of the action potential, are also accounted for. The model reproduces changes in the twitch amplitude and Ca2+-transients upon changes in muscle strain including the slow response. The model also reproduces the Bowditch effect and changes in the twitch amplitude and duration upon changes in the interstimulus interval, including accelerated relaxation at high stimulation frequency. Special efforts were taken to reduce the stiffness of the differential equations of the model. As a result, the equations can be integrated numerically with a relatively high time step making the model suitable for multiscale simulation of the human heart and allowing one to study the impact of myocardial mechanics on arrhythmias.

Evolution of mathematical models of cardiomyocyte electrophysiology

Advanced computational techniques and mathematical modeling have become more and more important to the study of cardiac electrophysiology. In this review, we provide a brief history of the evolution of cardiomyocyte electrophysiology models and highlight some of the most important ones that had a major impact on our understanding of the electrical activity of the myocardium and associated transmembrane ion fluxes in normal and pathological states. We also present the use of these models in the study of various arrhythmogenesis mechanisms, particularly the integration of experimental pharmacology data into advanced humanized models for in silico proarrhythmogenic risk prediction as an essential component of the Comprehensive in vitro Proarrhythmia Assay (CiPA) drug safety paradigm.

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.

Characterization of Spatio-Temporal Cardiac Action Potential Variability at Baseline and Under β-Adrenergic Stimulation by Combined Unscented Kalman Filter and Double Greedy Dimension Reduction

OBJECTIVE

Elevated spatio-temporal variability of human ventricular repolarization has been related to increased risk for ventricular arrhythmias and sudden cardiac death, particularly under β-adrenergic stimulation ( β-AS). This work presents a methodology for theoretical characterization of temporal and spatial repolarization variability at baseline conditions and in response to β-AS. For any measured voltage trace, the proposed methodology estimates the parameters and state variables of an underlying human ventricular action potential (AP) model by combining Double Greedy Dimension Reduction (DGDR) with automatic selection of biomarkers and the Unscented Kalman Filter (UKF). Such theoretical characterization can facilitate subsequent characterization of underlying variability mechanisms.

MATERIAL AND METHODS

Given an AP trace, initial estimates for the ionic conductances in a stochastic version of the baseline human ventricular O'Hara et al. model were obtained by DGDR. Those estimates served to initialize and update model parameter estimates by the UKF method based on formulation of an associated nonlinear state-space representation and joint estimation of model parameters and state variables. Similarly, β-AS-induced phosphorylation levels of cellular substrates were estimated by the DGDR-UKF methodology. Performance was tested by building an experimentally-calibrated population of virtual cells, from which synthetic AP traces were generated for baseline and β-AS conditions.

RESULTS

The combined DGDR-UKF methodology led to 25% reduction in the error associated with estimation of ionic current conductances at baseline conditions and phosphorylation levels under β-AS with respect to individual DGDR and UKF methods. This improvement was not at the expense of higher computational load, which was diminished by 90% with respect to the individual UKF method. Both temporal and spatial AP variability of repolarization were accurately characterized by the DGDR-UKF methodology.

CONCLUSIONS

A combined DGDR-UKF methodology is proposed for parameter and state variable estimation of human ventricular cell models from available AP traces at baseline and under β-AS. This methodology improves the estimation performance and reduces the convergence time with respect to individual DGDR and UKF methods and renders a suitable approach for computational characterization of spatio-temporal repolarization variability to be used for ascertainment of variability mechanisms and its relation to arrhythmogenesis.

Cardiomyocyte calcium handling in health and disease: Insights from in vitro and in silico studies

Calcium (Ca²⁺) plays a central role in cardiomyocyte excitation-contraction coupling. To ensure an optimal electrical impulse propagation and cardiac contraction, Ca²⁺ levels are regulated by a variety of Ca²⁺-handling proteins. In turn, Ca²⁺ modulates numerous electrophysiological processes. Accordingly, Ca²⁺-handling abnormalities can promote cardiac arrhythmias via various mechanisms, including the promotion of afterdepolarizations, ion-channel modulation and structural remodeling. In the last 30 years, significant improvements have been made in the computational modeling of cardiomyocyte Ca²⁺ handling under physiological and pathological conditions. However, numerous questions involving the Ca²⁺-dependent regulation of different macromolecular complexes, cross-talk between Ca²⁺-dependent regulatory pathways operating over a wide range of time scales, and bidirectional interactions between electrophysiology and mechanics remain to be addressed by in vitro and in silico studies. A better understanding of disease-specific Ca²⁺-dependent proarrhythmic mechanisms may facilitate the development of improved therapeutic strategies. In this review, we describe the fundamental mechanisms of cardiomyocyte Ca²⁺ handling in health and disease, and provide an overview of currently available computational models for cardiomyocyte Ca²⁺ handling. Finally, we discuss important uncertainties and open questions about cardiomyocyte Ca²⁺ handling and highlight how synergy between in vitro and in silico studies may help to answer several of these issues.

Mathematical model for β1-adrenergic regulation of the mouse ventricular myocyte contraction

The β1-adrenergic regulation of cardiac myocyte contraction plays an important role in regulating heart function. Activation of this system leads to an increased heart rate and stronger myocyte contraction. However, chronic stimulation of the β1-adrenergic signaling system can lead to cardiac hypertrophy and heart failure. To understand the mechanisms of action of β1-adrenoceptors, a mathematical model of cardiac myocyte contraction that includes the β1-adrenergic system was developed and studied. The model was able to simulate major experimental protocols for measurements of steady-state force-calcium relationships, cross-bridge release rate and force development rate, force-velocity relationship, and force redevelopment rate. It also reproduced quite well frequency and isoproterenol dependencies for intracellular Ca2+ concentration ([Ca2+]i) transients, total contraction force, and sarcomere shortening. The mathematical model suggested the mechanisms of increased contraction force and myocyte shortening on stimulation of β1-adrenergic receptors is due to phosphorylation of troponin I and myosin-binding protein C and increased [Ca2+]i transient resulting from activation of the β1-adrenergic signaling system. The model was used to simulate work-loop contractions and estimate the power during the cardiac cycle as well as the effects of 4-aminopyridine and tedisamil on the myocyte contraction. The developed mathematical model can be used further for simulations of contraction of ventricular myocytes from genetically modified mice and myocytes from mice with chronic cardiac diseases.

NEW & NOTEWORTHY A new mathematical model of mouse ventricular myocyte contraction that includes the β1-adrenergic system was developed. The model simulated major experimental protocols for myocyte contraction and predicted the effects of 4-aminopyridine and tedisamil on the myocyte contraction. The model also allowed for simulations of work-loop contractions and estimation of the power during the cardiac cycle.

Blockade of β-adrenergic signaling suppresses inflammasome and alleviates cardiac fibrosis

Background

Heart failure (HF) is an end-stage syndrome of all structural heart diseases which accompanies the loss of myocardium and cardiac fibrosis. Although the role of inflammasome in cardiac fibrosis has recently been a point of focus, the mechanism of inflammasome activation in HF has not yet been elucidated.

Methods

In this study, we investigated the expression of inflammasome proteins in a rat thoracic aorta constriction (TAC) model and cultured cardiac fibroblasts with stimulation of norepinephrine (NE).

Results

Our results showed that levels of inflammasome proteins in the myocardial of TAC rats were elevated. By blocking β-adrenergic signaling in the rats, inflammasome activation was suppressed and heart function was improved. The stimulation of cultured cardiac fibroblasts with NE activated inflammasome in vitro, which was abrogated by the inhibition of the calcium channels and reactive oxygen species (ROS). The activation of inflammasome by NE promoted cardiac fibrosis, whereas the inhibition of the calcium channels, ROS, and inflammasome reduced this effect.

Conclusions

The present study indicated that activation of inflammasome by β-adrenergic signaling promotes cardiac fibrosis. Therefore, modulation of inflammasome during HF might provide a novel strategy to treat this disease.

Whole-Cell cAMP and PKA Activity are Epiphenomena, Nanodomain Signaling Matters

Novel targeted fluorescent biosensors provide key insights into very local nanodomains of cAMP and PKA activity, and how they respond differently to β-adrenergic activation in cardiac myocytes. This unique spatiotemporal detail in living cells is not available with biochemical measurements of total cellular cAMP and PKA, and provides unique physiological insights.

A Microwell Cell Capture Device Reveals Variable Response to Dobutamine in Isolated Cardiomyocytes

Isolated ventricular cardiomyocytes exhibit substantial cell-to-cell variability, even when obtained from the same small volume of myocardium. In this study, we investigated the possibility that cardiomyocyte responses to β-adrenergic stimulus are also highly heterogeneous. To achieve the throughput and measurement duration desired for these experiments, we designed and validated a novel microwell system that immobilizes and uniformly orients isolated adult cardiomyocytes. In this configuration, detailed drug responses of dozens of cells can be followed for intervals exceeding 1 h. At the conclusion of an experiment, specific cells can also be harvested via a precision aspirator for single-cell gene expression profiling. Using this system, we followed changes in Ca²⁺ signaling and contractility of individual cells under sustained application of either dobutamine or omecamtiv mecarbil. Both compounds increased average cardiomyocyte contractility over the course of an hour, but responses of individual cells to dobutamine were significantly more variable. Surprisingly, some dobutamine-treated cardiomyocytes augmented Ca²⁺ release without increasing contractility. Other cells responded with increased contractility despite unchanged Ca²⁺ release. Single-cell gene expression analysis revealed significant co-expression of β-adrenergic pathway genes PKA regulatory subunit type I, PKA regulatory subunit type II, and Ca²⁺/calmodulin-dependent protein kinase II across cardiomyocytes. Other data supported a connection between the effects of dobutamine on relaxation rate and the expression of protein phosphatase 2. These findings suggest that variable drug responses among cells are not merely experimental artifacts. By enabling direct comparison of the functional behavior of an individual cell and the genes it expresses, this new system constitutes a unique tool for interrogating cardiomyocyte drug responses and discovering the genes that modulate them.

Different paths, same destination: divergent action potential responses produce conserved cardiac fight-or-flight response in mouse and rabbit hearts

KEY POINTS

Cardiac electrophysiology and Ca²⁺ handling change rapidly during the fight-or-flight response to meet physiological demands. Despite dramatic differences in cardiac electrophysiology, the cardiac fight-or-flight response is highly conserved across species. In this study, we performed physiological sympathetic nerve stimulation (SNS) while optically mapping cardiac action potentials and intracellular Ca²⁺ transients in innervated mouse and rabbit hearts. Despite similar heart rate and Ca²⁺ handling responses between mouse and rabbit hearts, we found notable species differences in spatio-temporal repolarization dynamics during SNS. Species-specific computational models revealed that these electrophysiological differences allowed for enhanced Ca²⁺ handling (i.e. enhanced inotropy) in each species, suggesting that electrophysiological responses are fine-tuned across species to produce optimal cardiac fight-or-flight responses.

ABSTRACT

Sympathetic activation of the heart results in positive chronotropy and inotropy, which together rapidly increase cardiac output. The precise mechanisms that produce the electrophysiological and Ca²⁺ handling changes underlying chronotropic and inotropic responses have been studied in detail in isolated cardiac myocytes. However, few studies have examined the dynamic effects of physiological sympathetic nerve activation on cardiac action potentials (APs) and intracellular Ca²⁺ transients (CaTs) in the intact heart. Here, we performed bilateral sympathetic nerve stimulation (SNS) in fully innervated, Langendorff-perfused rabbit and mouse hearts. Dual optical mapping with voltage- and Ca²⁺ -sensitive dyes allowed for analysis of spatio-temporal AP and CaT dynamics. The rabbit heart responded to SNS with a monotonic increase in heart rate (HR), monotonic decreases in AP and CaT duration (APD, CaTD), and a monotonic increase in CaT amplitude. The mouse heart had similar HR and CaT responses; however, a pronounced biphasic APD response occurred, with initial prolongation (50.9 ± 5.1 ms at t = 0 s vs. 60.6 ± 4.1 ms at t = 15 s, P < 0.05) followed by shortening (46.5 ± 9.1 ms at t = 60 s, P = NS vs. t = 0). We determined the biphasic APD response in mouse was partly due to dynamic changes in HR during SNS and was exacerbated by β-adrenergic activation. Simulations with species-specific cardiac models revealed that transient APD prolongation in mouse allowed for greater and more rapid CaT responses, suggesting more rapid increases in contractility; conversely, the rabbit heart requires APD shortening to produce optimal inotropic responses. Thus, while the cardiac fight-or-flight response is highly conserved between species, the underlying mechanisms orchestrating these effects differ significantly.

Enhanced Depolarization Drive in Failing Rabbit Ventricular Myocytes: Calcium-Dependent and β-Adrenergic Effects on Late Sodium, L-Type Calcium, and Sodium-Calcium Exchange Currents

BACKGROUND

Heart failure (HF) is characterized by electrophysiological remodeling resulting in increased risk of cardiac arrhythmias. Previous reports suggest that elevated inward ionic currents in HF promote action potential (AP) prolongation, increased short-term variability of AP repolarization, and delayed afterdepolarizations. However, the underlying changes in late Na+ current (INaL), L-type Ca2+ current, and NCX (Na+/Ca2+ exchanger) current are often measured in nonphysiological conditions (square-pulse voltage clamp, slow pacing rates, exogenous Ca2+ buffers).

METHODS

We measured the major inward currents and their Ca2+- and β-adrenergic dependence under physiological AP clamp in rabbit ventricular myocytes in chronic pressure/volume overload-induced HF (versus age-matched control).

RESULTS

AP duration and short-term variability of AP repolarization were increased in HF, and importantly, inhibition of INaL decreased both parameters to the control level. INaL was slightly increased in HF versus control even when intracellular Ca2+ was strongly buffered. But under physiological AP clamp with normal Ca2+ cycling, INaL was markedly upregulated in HF versus control (dependent largely on CaMKII [Ca2+/calmodulin-dependent protein kinase II] activity). β-Adrenergic stimulation (often elevated in HF) further enhanced INaL. L-type Ca2+ current was decreased in HF when Ca2+ was buffered, but CaMKII-mediated Ca2+-dependent facilitation upregulated physiological L-type Ca2+ current to the control level. Furthermore, L-type Ca2+ current response to β-adrenergic stimulation was significantly attenuated in HF. Inward NCX current was upregulated at phase 3 of AP in HF when assessed by combining experimental data and computational modeling.

CONCLUSIONS

Our results suggest that CaMKII-dependent upregulation of INaL in HF significantly contributes to AP prolongation and increased short-term variability of AP repolarization, which may lead to increased arrhythmia propensity, and is further exacerbated by adrenergic stress.

Cardiac CaV1.2 channels require β subunits for β-adrenergic-mediated modulation but not trafficking

Ca2+ channel β-subunit interactions with pore-forming α-subunits are long-thought to be obligatory for channel trafficking to the cell surface and for tuning of basal biophysical properties in many tissues. Unexpectedly, we demonstrate that transgenic expression of mutant α1C subunits lacking capacity to bind CaVβ can traffic to the sarcolemma in adult cardiomyocytes in vivo and sustain normal excitation-contraction coupling. However, these β-less Ca2+ channels cannot be stimulated by β-adrenergic pathway agonists, and thus adrenergic augmentation of contractility is markedly impaired in isolated cardiomyocytes and in hearts. Similarly, viral-mediated expression of a β-subunit-sequestering peptide sharply curtailed β-adrenergic stimulation of WT Ca2+ channels, identifying an approach to specifically modulate β-adrenergic regulation of cardiac contractility. Our data demonstrate that β subunits are required for β-adrenergic regulation of CaV1.2 channels and positive inotropy in the heart, but are dispensable for CaV1.2 trafficking to the adult cardiomyocyte cell surface, and for basal function and excitation-contraction coupling.

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.

Experimental assessment of a myocyte-based multiscale model of cardiac contractile dysfunction

Cardiac contractile dysfunction (CD) is a multifactorial syndrome caused by different acute or progressive diseases which hamper assessing the role of the underlying mechanisms characterizing a defined pathological condition. Mathematical modeling can help to understand the processes involved in CD and analyze their relative impact in the overall response. The aim of this study was thus to use a myocyte-based multiscale model of the circulatory system to simulate the effects of halothane, a volatile anesthetic which at high doses elicits significant acute CD both in isolated myocytes and intact animals. Ventricular chambers built using a human myocyte model were incorporated into a whole circulatory system represented by resistances and capacitances. Halothane-induced decreased sarco(endo)plasmic reticulum Ca²⁺ (SERCA2a) reuptake pump, transient outward K⁺ (I_to), Na⁺-Ca²⁺ exchanger (I_NCX) and L-type Ca²⁺ channel (I_CaL) currents, together with ryanodine receptor (RyR2) increased open probability (P_o) and reduced myofilament Ca²⁺ sensitivity, reproduced equivalent decreased action potential duration at 90% repolarization and intracellular Ca²⁺ concentration at the myocyte level reported in the literature. In the whole circulatory system, model reduction in mean arterial pressure, cardiac output and regional wall thickening fraction was similar to experimental results in open-chest sheep subjected to acute halothane overdose. Effective model performance indicates that the model structure could be used to study other changes in myocyte targets eliciting CD.

Quantitative systems models illuminate arrhythmia mechanisms in heart failure: Role of the Na+ -Ca2+ -Ca2+ /calmodulin-dependent protein kinase II-reactive oxygen species feedback

Quantitative systems modeling aims to integrate knowledge in different research areas with models describing biological mechanisms and dynamics to gain a better understanding of complex clinical syndromes. Heart failure (HF) is a chronic complex cardiac disease that results from structural or functional disorders impairing the ability of the ventricle to fill with or eject blood. Highly interactive and dynamic changes in mechanical, structural, neurohumoral, metabolic, and electrophysiological properties collectively predispose the failing heart to cardiac arrhythmias, which are responsible for about a half of HF deaths. Multiscale cardiac modeling and simulation integrate structural and functional data from HF experimental models and patients to improve our mechanistic understanding of this complex arrhythmia syndrome. In particular, they allow investigating how disease-induced remodeling alters the coupling of electrophysiology, Ca²⁺ and Na⁺ handling, contraction, and energetics that lead to rhythm derangements. The Ca²⁺ /calmodulin-dependent protein kinase II, which expression and activity are enhanced in HF, emerges as a critical hub that modulates the feedbacks between these various subsystems and promotes arrhythmogenesis. This article is categorized under: Physiology > Mammalian Physiology in Health and Disease Models of Systems Properties and Processes > Mechanistic Models Models of Systems Properties and Processes > Cellular Models Models of Systems Properties and Processes > Organ, Tissue, and Physiological Models.

Automatic Optimization of an in Silico Model of Human iPSC Derived Cardiomyocytes Recapitulating Calcium Handling Abnormalities

The growing importance of human induced pluripotent stem cell-derived cardiomyoyctes (hiPSC-CMs), as patient-specific and disease-specific models for studying cellular cardiac electrophysiology or for preliminary cardiotoxicity tests, generated better understanding of hiPSC-CM biophysical mechanisms and great amount of action potential and calcium transient data. In this paper, we propose a new hiPSC-CM in silico model, with particular attention to Ca²⁺ handling. We used (i) the hiPSC-CM Paci2013 model as starting point, (ii) a new dataset of Ca²⁺ transient measurements to tune the parameters of the inward and outward Ca²⁺ fluxes of sarcoplasmic reticulum, and (iii) an automatic parameter optimization to fit action potentials and Ca²⁺ transients. The Paci2018 model simulates, together with the typical hiPSC-CM spontaneous action potentials, more refined Ca²⁺ transients and delayed afterdepolarizations-like abnormalities, which the old Paci2013 was not able to predict due to its mathematical formulation. The Paci2018 model was validated against (i) the same current blocking experiments used to validate the Paci2013 model, and (ii) recently published data about effects of different extracellular ionic concentrations. In conclusion, we present a new and more versatile in silico model, which will provide a platform for modeling the effects of drugs or mutations that affect Ca²⁺ handling in hiPSC-CMs.

Cardiophysiology Illustrated by Comparing Ventricular Volumes in Healthy Adult Males and Females

Recent advances in cardiac imaging techniques have substantially contributed to a growing interest in the analysis of global cardiac chamber dimensions and regional myocardial deformation. During the cardiac cycle, ventricular luminal volume varies due to the contraction process, which also confers a shape change including substantial alteration of long axis length, as well as rotation of the base compared to the apex. Local deformation can be assessed by strain (rate) analysis. Reviewing the present literature, it must be concluded that there is no single metric available to comprehensively characterize ventricular function. Every candidate advanced thus far has been found to incompletely reflect ventricular performance. This observation is not surprising in view of the complexity of the cardiac pump system. Additionally, sex-specific modifiers may play a role. More than three decades ago, it was shown that on average the ventricular volume is smaller in healthy women compared to matched males. Therefore, the present contribution concerns the interpretation of data derived from the healthy heart in both men and women. Starting from the classical Starling concept, we apply a simple mathematical transformation which permits an insightful representation of ventricular mechanics. Relating end-systolic volume (ESV) to end-diastolic volume creates the ventricular volume regulation graph which features the pertinent working point of an individual heart. This fundamental approach illustrates why certain proposed performance indexes cannot individually reveal the essence of ventricular systolic function. We demonstrate that particular metrics are highly interconnected and just tell us the same story in a different disguise. It is imperative to understand which associations exist and if they expectedly are (nearly) linear or frankly nonlinear. Notably, ejection fraction (EF) is primarily determined by ESV, while in turn EF is not much different from ventriculo-arterial coupling (VAC). Insight into cardiac function is promoted by identification of the paramount/essential components involved. The smaller ESV (p < 0.0001) implies that EF is higher in women and may also have consequences for VAC.

Stretch-Activated Current Can Promote or Suppress Cardiac Alternans Depending on Voltage-Calcium Interaction

Cardiac alternans has been linked to the onset of ventricular fibrillation and ventricular tachycardia, leading to life-threatening arrhythmias. Here, we investigated the effects of stretch-activated currents (ISAC) on alternans using a physiologically detailed model of the ventricular myocyte. We found that increasing ISAC suppresses alternans if the voltage-Ca coupling is positive or the alternans is voltage driven. However, for electromechanically discordant alternans, which occurs when the alternans is Ca driven with negative voltage-Ca coupling, increasing ISAC promotes Ca alternans. In addition, if action potential duration-Ca transients show quasiperiodicity, we observe a biphasic effect of ISAC, i.e., suppressing quasiperiodic oscillation at small stretch but promoting electromechanically discordant alternans at larger stretch. Our results demonstrate how ISAC interacts with coupled voltage-Ca dynamical systems with respect to alternans.

β-Adrenergic induced SR Ca2+ leak is mediated by an Epac-NOS pathway

Cardiac β-adrenergic receptors (β-AR) and Ca²⁺-Calmodulin dependent protein kinase (CaMKII) regulate both physiological and pathophysiological Ca²⁺ signaling. Elevated diastolic Ca²⁺ leak from the sarcoplasmic reticulum (SR) contributes to contractile dysfunction in heart failure and to arrhythmogenesis. β-AR activation is known to increase SR Ca²⁺ leak via CaMKII-dependent phosphorylation of the ryanodine receptor. Two independent and reportedly parallel pathways have been implicated in this β-AR-CaMKII cascade, one involving exchange protein directly activated by cAMP (Epac2) and another involving nitric oxide synthase 1 (NOS1). Here we tested whether Epac and NOS function in a single series pathway to increase β-AR induced and CaMKII-dependent SR Ca²⁺ leak. Leak was measured as both Ca²⁺ spark frequency and tetracaine-induced shifts in SR Ca²⁺, in mouse and rabbit ventricular myocytes. Direct Epac activation by 8-CPT (8-(4-chlorophenylthio)-2'-O-methyl-cAMP) mimicked β-AR-induced SR Ca²⁺ leak, and both were blocked by NOS inhibition. The same was true for myocyte CaMKII activation (assessed via a FRET-based reporter) and ryanodine receptor phosphorylation. Inhibitor and phosphorylation studies also implicated phosphoinositide 3-kinase (PI3K) and protein kinase B (Akt) downstream of Epac and above NOS activation in this pathway. We conclude that these two independently characterized parallel pathways function mainly via a single series arrangement (β-AR-cAMP-Epac-PI3K-Akt-NOS1-CaMKII) to mediate increased SR Ca²⁺ leak. Thus, for β-AR activation the cAMP-PKA branch effects inotropy and lusitropy (by effects on Ca²⁺ current and SR Ca²⁺-ATPase), this cAMP-Epac-NOS pathway increases pathological diastolic SR Ca²⁺leak. This pathway distinction may allow novel SR Ca²⁺ leak therapeutic targeting in treatment of arrhythmias in heart failure that spare the inotropic and lusitropic effects of the PKA branch.

FRET biosensor uncovers cAMP nano-domains at β-adrenergic targets that dictate precise tuning of cardiac contractility

Compartmentalized cAMP/PKA signalling is now recognized as important for physiology and pathophysiology, yet a detailed understanding of the properties, regulation and function of local cAMP/PKA signals is lacking. Here we present a fluorescence resonance energy transfer (FRET)-based sensor, CUTie, which detects compartmentalized cAMP with unprecedented accuracy. CUTie, targeted to specific multiprotein complexes at discrete plasmalemmal, sarcoplasmic reticular and myofilament sites, reveals differential kinetics and amplitudes of localized cAMP signals. This nanoscopic heterogeneity of cAMP signals is necessary to optimize cardiac contractility upon adrenergic activation. At low adrenergic levels, and those mimicking heart failure, differential local cAMP responses are exacerbated, with near abolition of cAMP signalling at certain locations. This work provides tools and fundamental mechanistic insights into subcellular adrenergic signalling in normal and pathological cardiac function.

Impact of RyR2 potentiation on myocardial function

This perspective attempts to shed light on an old and not yet solved controversy in cardiac physiology, i.e., the impact of increasing ryanodine receptor (RyR)2 open probability on myocardial function. Based on an already proven myocyte model, it was shown that increasing RyR2 open probability results in a purely short-lived increase in Ca²⁺ transient amplitude, and, therefore, it does not increase cardiac contractility. However, potentiation of RyR2 activity permanently enhances fractional Ca²⁺ release, shifting the intracellular Ca²⁺ transient versus sarcoplasmic reticulum (SR) Ca²⁺ content curve to a new state of higher efficiency. This would allow the heart to maintain a given contractility despite a decrease in SR Ca²⁺ content, to enhance contractility if SR Ca²⁺ content is simultaneously preserved or to successfully counteract the effects of a negative inotropic intervention.NEW & NOTEWORTHY Increasing ryanodine receptor (RyR)2 open probability does not increase cardiac contractility. However, RyR2 potentiation shifts the intracellular Ca²⁺ transient-sarcoplasmic reticulum (SR) Ca²⁺ content relationship toward an enhanced efficiency state, which may contribute to a positive inotropic effect, preserve contractility despite decreased SR Ca²⁺ content, or successfully counteract the effects of a negative inotropic action.

Quantitative analysis of the Ca2+ -dependent regulation of delayed rectifier K+ current IKs in rabbit ventricular myocytes

KEY POINTS

[Ca²⁺ ]_i enhanced rabbit ventricular slowly activating delayed rectifier K⁺ current (I_Ks ) by negatively shifting the voltage dependence of activation and slowing deactivation, similar to perfusion of isoproterenol. Rabbit ventricular rapidly activating delayed rectifier K⁺ current (I_Kr ) amplitude and voltage dependence were unaffected by high [Ca²⁺ ]_i . When measuring or simulating I_Ks during an action potential, I_Ks was not different during a physiological Ca²⁺ transient or when [Ca²⁺ ]_i was buffered to 500 nm.

ABSTRACT

The slowly activating delayed rectifier K⁺ current (I_Ks ) contributes to repolarization of the cardiac action potential (AP). Intracellular Ca²⁺ ([Ca²⁺ ]_i ) and β-adrenergic receptor (β-AR) stimulation modulate I_Ks amplitude and kinetics, but details of these important I_Ks regulators and their interaction are limited. We assessed the [Ca²⁺ ]_i dependence of I_Ks in steady-state conditions and with dynamically changing membrane potential and [Ca²⁺ ]_i during an AP. I_Ks was recorded from freshly isolated rabbit ventricular myocytes using whole-cell patch clamp. With intracellular pipette solutions that controlled free [Ca²⁺ ]_i , we found that raising [Ca²⁺ ]_i from 100 to 600 nm produced similar increases in I_Ks as did β-AR activation, and the effects appeared additive. Both β-AR activation and high [Ca²⁺ ]_i increased maximally activated tail I_Ks , negatively shifted the voltage dependence of activation, and slowed deactivation kinetics. These data informed changes in our well-established mathematical model of the rabbit myocyte. In both AP-clamp experiments and simulations, I_Ks recorded during a normal physiological Ca²⁺ transient was similar to I_Ks measured with [Ca²⁺ ]_i clamped at 500-600 nm. Thus, our study provides novel quantitative data as to how physiological [Ca²⁺ ]_i regulates I_Ks amplitude and kinetics during the normal rabbit AP. Our results suggest that micromolar [Ca²⁺ ]_i , in the submembrane or junctional cleft space, is not required to maximize [Ca²⁺ ]_i -dependent I_Ks activation during normal Ca²⁺ transients.

Multiscale model of the human cardiovascular system: Description of heart failure and comparison of contractility indices

A multiscale model of the cardiovascular system is presented. Hemodynamics is described by a lumped parameter model, while heart contraction is described at the cellular scale. An electrophysiological model and a mechanical model were coupled and adjusted so that the pressure and volume of both ventricles are linked to the force and length of a half-sarcomere. Particular attention was paid to the extreme values of the sarcomere length, which must keep physiological values. This model is able to reproduce healthy behavior, preload variations experiments, and ventricular failure. It also allows to compare the relevance of standard cardiac contractility indices. This study shows that the theoretical gold standard for assessing cardiac contractility, namely the end-systolic elastance, is actually load-dependent and therefore not a reliable index of cardiac contractility.

Phospholamban ablation rescues the enhanced propensity to arrhythmias of mice with CaMKII-constitutive phosphorylation of RyR2 at site S2814

KEY POINTS

Mice with Ca(2+) -calmodulin-dependent protein kinase (CaMKII) constitutive pseudo-phosphorylation of the ryanodine receptor RyR2 at Ser2814 (S2814D(+/+) mice) exhibit a higher open probability of RyR2, higher sarcoplasmic reticulum (SR) Ca(2+) leak in diastole and increased propensity to arrhythmias under stress conditions. We generated phospholamban (PLN)-deficient S2814D(+/+) knock-in mice by crossing two colonies, S2814D(+/+) and PLNKO mice, to test the hypothesis that PLN ablation can prevent the propensity to arrhythmias of S2814D(+/+) mice. PLN ablation partially rescues the altered intracellular Ca(2+) dynamics of S2814D(+/+) hearts and myocytes, but enhances SR Ca(2+) sparks and leak on confocal microscopy. PLN ablation diminishes ventricular arrhythmias promoted by CaMKII phosphorylation of S2814 on RyR2. PLN ablation aborts the arrhythmogenic SR Ca(2+) waves of S2814D(+/+) and transforms them into non-propagating events. A mathematical human myocyte model replicates these results and predicts the increase in SR Ca(2+) uptake required to prevent the arrhythmias induced by a CaMKII-dependent leaky RyR2.

ABSTRACT

Mice with constitutive pseudo-phosphorylation at Ser2814-RyR2 (S2814D(+/+) ) have increased propensity to arrhythmias under β-adrenergic stress conditions. Although abnormal Ca(2+) release from the sarcoplasmic reticulum (SR) has been linked to arrhythmogenesis, the role played by SR Ca(2+) uptake remains controversial. We tested the hypothesis that an increase in SR Ca(2+) uptake is able to rescue the increased arrhythmia propensity of S2814D(+/+) mice. We generated phospholamban (PLN)-deficient/S2814D(+/+) knock-in mice by crossing two colonies, S2814D(+/+) and PLNKO mice (SD(+/+) /KO). SD(+/+) /KO myocytes exhibited both increased SR Ca(2+) uptake seen in PLN knock-out (PLNKO) myocytes and diminished SR Ca(2+) load (relative to PLNKO), a characteristic of S2814D(+/+) myocytes. Ventricular arrhythmias evoked by catecholaminergic challenge (caffeine/adrenaline) in S2814D(+/+) mice in vivo or programmed electric stimulation and high extracellular Ca(2+) in S2814D(+) /(-) hearts ex vivo were significantly diminished by PLN ablation. At the myocyte level, PLN ablation converted the arrhythmogenic Ca(2+) waves evoked by high extracellular Ca(2+) provocation in S2814D(+/+) mice into non-propagated Ca(2+) mini-waves on confocal microscopy. Myocyte Ca(2+) waves, typical of S2814D(+/+) mice, could be evoked in SD(+/+) /KO cells by partially inhibiting SERCA2a. A mathematical human myocyte model replicated these results and allowed for predicting the increase in SR Ca(2+) uptake required to prevent the arrhythmias induced by a Ca(2+) -calmodulin-dependent protein kinase (CaMKII)-dependent leaky RyR2. Our results demonstrate that increasing SR Ca(2+) uptake by PLN ablation can prevent the arrhythmic events triggered by SR Ca(2+) leak due to CaMKII-dependent phosphorylation of the RyR2-S2814 site and underscore the benefits of increasing SERCA2a activity on SR Ca(2+) -triggered arrhythmias.

Modeling calcium regulation of contraction, energetics, signaling, and transcription in the cardiac myocyte

Calcium (Ca(2+)) plays many important regulatory roles in cardiac muscle cells. In the initial phase of the action potential, influx of Ca(2+) through sarcolemmal voltage-gated L-type Ca(2+) channels (LCCs) acts as a feed-forward signal that triggers a large release of Ca(2+) from the junctional sarcoplasmic reticulum (SR). This Ca(2+) drives heart muscle contraction and pumping of blood in a process known as excitation-contraction coupling (ECC). Triggered and released Ca(2+) also feed back to inactivate LCCs, attenuating the triggered Ca(2+) signal once release has been achieved. The process of ECC consumes large amounts of ATP. It is now clear that in a process known as excitation-energetics coupling, Ca(2+) signals exert beat-to-beat regulation of mitochondrial ATP production that closely couples energy production with demand. This occurs through transport of Ca(2+) into mitochondria, where it regulates enzymes of the tricarboxylic acid cycle. In excitation-signaling coupling, Ca(2+) activates a number of signaling pathways in a feed-forward manner. Through effects on their target proteins, these interconnected pathways regulate Ca(2+) signals in complex ways to control electrical excitability and contractility of heart muscle. In a process known as excitation-transcription coupling, Ca(2+) acting primarily through signal transduction pathways also regulates the process of gene transcription. Because of these diverse and complex roles, experimentally based mechanistic computational models are proving to be very useful for understanding Ca(2+) signaling in the cardiac myocyte.

S-Nitrosylation of Calcium-Handling Proteins in Cardiac Adrenergic Signaling and Hypertrophy

RATIONALE

The regulation of calcium (Ca(2+)) homeostasis by β-adrenergic receptor (βAR) activation provides the essential underpinnings of sympathetic regulation of myocardial function, as well as a basis for understanding molecular events that result in hypertrophic signaling and heart failure. Sympathetic stimulation of the βAR not only induces protein phosphorylation but also activates nitric oxide-dependent signaling, which modulates cardiac contractility. Nonetheless, the role of nitric oxide in βAR-dependent regulation of Ca(2+) handling has not yet been explicated fully.

OBJECTIVE

To elucidate the role of protein S-nitrosylation, a major transducer of nitric oxide bioactivity, on βAR-dependent alterations in cardiomyocyte Ca(2+) handling and hypertrophy.

METHODS AND RESULTS

Using transgenic mice to titrate the levels of protein S-nitrosylation, we uncovered major roles for protein S-nitrosylation, in general, and for phospholamban and cardiac troponin C S-nitrosylation, in particular, in βAR-dependent regulation of Ca(2+) homeostasis. Notably, S-nitrosylation of phospholamban consequent upon βAR stimulation is necessary for the inhibitory pentamerization of phospholamban, which activates sarcoplasmic reticulum Ca(2+)-ATPase and increases cytosolic Ca(2+) transients. Coincident S-nitrosylation of cardiac troponin C decreases myocardial sensitivity to Ca(2+). During chronic adrenergic stimulation, global reductions in cellular S-nitrosylation mitigate hypertrophic signaling resulting from Ca(2+) overload.

CONCLUSIONS

S-Nitrosylation operates in concert with phosphorylation to regulate many cardiac Ca(2+)-handling proteins, including phospholamban and cardiac troponin C, thereby playing an essential and previously unrecognized role in cardiac Ca(2+) homeostasis. Manipulation of the S-nitrosylation level may prove therapeutic in heart failure.

Pompe disease: Shared and unshared features of lysosomal storage disorders

Pompe disease, an inherited deficiency of lysosomal acid α-glucosidase (GAA), is a severe metabolic myopathy with a wide range of clinical manifestations. It is the first recognized lysosomal storage disorder and the first neuromuscular disorder for which a therapy (enzyme replacement) has been approved. As GAA is the only enzyme that hydrolyses glycogen to glucose in the acidic environment of the lysosome, its deficiency leads to glycogen accumulation within and concomitant enlargement of this organelle. Since the introduction of the therapy, the overall understanding of the disease has progressed significantly, but the pathophysiology of muscle damage is still not fully understood. The emerging complex picture of the pathological cascade involves disturbance of calcium homeostasis, mitochondrial abnormalities, dysfunctional autophagy, accumulation of toxic undegradable materials, and accelerated production of lipofuscin deposits that are unrelated to aging. The relationship of Pompe disease to other lysosomal storage disorders and potential therapeutic interventions for Pompe disease are discussed.