Data-based theoretical identification of subcellular calcium compartments and estimation of calcium dynamics in cardiac myocytes

A computational model of spatio-temporal cardiac intracellular calcium handling with realistic structure and spatial flux distribution from sarcoplasmic reticulum and t-tubule reconstructions

Intracellular calcium cycling is a vital component of cardiac excitation-contraction coupling. The key structures responsible for controlling calcium dynamics are the cell membrane (comprising the surface sarcolemma and transverse-tubules), the intracellular calcium store (the sarcoplasmic reticulum), and the co-localisation of these two structures to form dyads within which calcium-induced-calcium-release occurs. The organisation of these structures tightly controls intracellular calcium dynamics. In this study, we present a computational model of intracellular calcium cycling in three-dimensions (3-D), which incorporates high resolution reconstructions of these key regulatory structures, attained through imaging of tissue taken from the sheep left ventricle using serial block face scanning electron microscopy. An approach was developed to model the sarcoplasmic reticulum structure at the whole-cell scale, by reducing its full 3-D structure to a 3-D network of one-dimensional strands. The model reproduces intracellular calcium dynamics during control pacing and reveals the high-resolution 3-D spatial structure of calcium gradients and intracellular fluxes in both the cytoplasm and sarcoplasmic reticulum. We also demonstrated the capability of the model to reproduce potentially pro-arrhythmic dynamics under perturbed conditions, pertaining to calcium-transient alternans and spontaneous release events. Comparison with idealised cell models emphasised the importance of structure in determining calcium gradients and controlling the spatial dynamics associated with calcium-transient alternans, wherein the probabilistic nature of dyad activation and recruitment was constrained. The model was further used to highlight the criticality in calcium spark propagation in relation to inter-dyad distances. The model presented provides a powerful tool for future investigation of structure-function relationships underlying physiological and pathophysiological intracellular calcium handling phenomena at the whole-cell. The approach allows for the first time direct integration of high-resolution images of 3-D intracellular structures with models of calcium cycling, presenting the possibility to directly assess the functional impact of structural remodelling at the cellular scale.

The organisation of the membrane and sub-cellular structures of cells in the heart closely controls the coupling between its electrical and mechanical function. Computational models of the cellular calcium handling system, which is responsible for this electro-mechanical coupling, have been developed in recent years to study underlying structure-function relationships. Previous models have been largely idealised in structure; we present a new model which incorporates experimental data describing the high-resolution organisation of the primary structures involved in calcium dynamics. Significantly, the structure of the intracellular calcium store is modelled for the first time. The model is shown to reproduce calcium dynamics in control cells in both normal and abnormal conditions, demonstrating its suitability for future investigation of structure-function relationships. Thus, the model presented provides a powerful tool for the direct integration of experimentally acquired structural data in healthy and diseased cells and assessment of the role of structure in regulating normal and abnormal calcium dynamics.

The Effect of a Novel Highly Selective Inhibitor of the Sodium/Calcium Exchanger (NCX) on Cardiac Arrhythmias in In Vitro and In Vivo Experiments

Background

In this study the effects of a new, highly selective sodium-calcium exchanger (NCX) inhibitor, ORM-10962 were investigated on cardiac NCX current, Ca²⁺ transients, cell shortening and in experimental arrhythmias. The level of selectivity of the novel inhibitor on several major transmembrane ion currents (L-type Ca²⁺ current, major repolarizing K⁺ currents, late Na⁺ current, Na⁺/K⁺ pump current) was also determined.

Methods

Ion currents in single dog ventricular cells (cardiac myocytes; CM), and action potentials in dog cardiac multicellular preparations were recorded utilizing the whole-cell patch clamp and standard microelectrode techniques, respectively. Ca²⁺ transients and cell shortening were measured in fluorescent dye loaded isolated dog myocytes. Antiarrhythmic effects of ORM-10962 were studied in anesthetized ouabain (10 μg/kg/min i.v.) pretreated guinea pigs and in ischemia-reperfusion models (I/R) of anesthetized coronary artery occluded rats and Langendorff perfused guinea pigs hearts.

Results

ORM-10962 significantly reduced the inward/outward NCX currents with estimated EC50 values of 55/67 nM, respectively. The compound, even at a high concentration of 1 μM, did not modify significantly the magnitude of I_CaL in CMs, neither had any apparent influence on the inward rectifier, transient outward, the rapid and slow components of the delayed rectifier potassium currents, the late and peak sodium and Na⁺/K⁺ pump currents. NCX inhibition exerted moderate positive inotropic effect under normal condition, negative inotropy when reverse, and further positive inotropic effect when forward mode was facilitated. In dog Purkinje fibres 1 μM ORM-10962 decreased the amplitude of digoxin induced delayed afterdepolarizations (DADs). Pre-treatment with 0.3 mg/kg ORM-10962 (i.v.) 10 min before starting ouabain infusion significantly delayed the development and recurrence of ventricular extrasystoles (by about 50%) or ventricular tachycardia (by about 30%) in anesthetized guinea pigs. On the contrary, ORM-10962 pre-treatment had no apparent influence on the time of onset or the severity of I/R induced arrhythmias in anesthetized rats and in Langendorff perfused guinea-pig hearts.

Conclusions

The present study provides strong evidence for a high efficacy and selectivity of the NCX-inhibitory effect of ORM-10962. Selective NCX inhibition can exert positive as well as negative inotropic effect depending on the actual operation mode of NCX. Selective NCX blockade may contribute to the prevention of DAD based arrhythmogenesis, in vivo, however, its effect on I/R induced arrhythmias is still uncertain.

Modeling Na+-Ca2+ exchange in the heart: Allosteric activation, spatial localization, sparks and excitation-contraction coupling

The cardiac sodium (Na⁺)/calcium (Ca²⁺) exchanger (NCX1) is an electrogenic membrane transporter that regulates Ca²⁺ homeostasis in cardiomyocytes, serving mainly to extrude Ca²⁺ during diastole. The direction of Ca²⁺ transport reverses at membrane potentials near that of the action potential plateau, generating an influx of Ca²⁺ into the cell. Therefore, there has been great interest in the possible roles of NCX1 in cardiac Ca²⁺-induced Ca²⁺ release (CICR). Interest has been reinvigorated by a recent super-resolution optical imaging study suggesting that ~18% of NCX1 co-localize with ryanodine receptor (RyR2) clusters, and ~30% of additional NCX1 are localized to within ~120nm of the nearest RyR2. NCX1 may therefore occupy a privileged position in which to modulate CICR. To examine this question, we have developed a mechanistic biophysically-detailed model of NCX1 that describes both NCX1 transport kinetics and Ca²⁺-dependent allosteric regulation. This NCX1 model was incorporated into a previously developed super-resolution model of the Ca²⁺ spark as well as a computational model of the cardiac ventricular myocyte that includes a detailed description of CICR with stochastic gating of L-type Ca²⁺ channels and RyR2s, and that accounts for local Ca²⁺ gradients near the dyad via inclusion of a peri-dyadic (PD) compartment. Both models predict that increasing the fraction of NCX1 in the dyad and PD decreases spark frequency, fidelity, and diastolic Ca²⁺ levels. Spark amplitude and duration are less sensitive to NCX1 spatial redistribution. On the other hand, NCX1 plays an important role in promoting Ca²⁺ entry into the dyad, and hence contributing to the trigger for RyR2 release at depolarized membrane potentials and in the presence of elevated local Na⁺ concentration. Whole-cell simulation of NCX1 tail currents are consistent with the finding that a relatively high fraction of NCX1 (~45%) resides in the dyadic and PD spaces, with a dyad-to-PD ratio of roughly 1:2. Allosteric Ca²⁺ activation of NCX1 helps to "functionally localize" exchanger activity to the dyad and PD by reducing exchanger activity in the cytosol thereby protecting the cell from excessive loss of Ca²⁺ during diastole.

Induced NCX1 overexpression attenuates pressure overload-induced pathological cardiac remodelling

AIMS

Although increased Na(+)/Ca(2+) exchanger 1 (NCX1) expression is observed during heart failure (HF), the pathological role of NCX1 during the progression of HF remains unclear. We examined alterations of NCX1 expression and activity in hearts after transverse aortic constriction (TAC) surgery and explored whether NCX1 influences pressure overload-induced pathological cardiac remodelling.

METHODS AND RESULTS

We generated novel transgenic mice in which NCX1 expression is controlled by a cardiac-specific, doxycycline (DOX)-dependent promoter. In the absence of DOX, TAC surgery caused substantial chamber dilation with a gradual decrease in contractility by 16 weeks. Cardiomyocytes showed a decline in contractility with abnormal Ca(2+) handling during excitation-contraction (E-C) coupling. Reduced NCX1 activity was observed 8 weeks after TAC and was still apparent at 17 weeks. Induced NCX1 overexpression by DOX treatment starting 8 weeks after TAC returned NCX1 activity to pre-TAC levels and prevented chamber dilation with cardiac dysfunction. DOX treatment not only upregulated NCX1 expression in TAC-operated hearts but also returned L-type Ca(2+) channel and sarcoplasmic reticulum (SR) Ca(2+) ATPase expression levels to those in sham-operated hearts. In DOX-treated myocytes, contractility, T-tubule integrity, synchrony of Ca(2+) release from the SR, and Ca(2+) handling during E-C coupling was preserved 16 weeks after TAC surgery. In addition, DOX treatment attenuated the down-regulation of survival signalling and up-regulation of apoptosis signalling 16 weeks after TAC surgery.

CONCLUSION

Induced overexpression of NCX1 attenuated pressure overload-induced pathological cardiac remodelling. Thus, maintaining NCX1 activity may be a potential therapeutic strategy for preventing the progression of HF.

A human ventricular myocyte model with a refined representation of excitation-contraction coupling

Cardiac Ca(2+)-induced Ca(2+) release (CICR) occurs by a regenerative activation of ryanodine receptors (RyRs) within each Ca(2+)-releasing unit, triggered by the activation of L-type Ca(2+) channels (LCCs). CICR is then terminated, most probably by depletion of Ca(2+) in the junctional sarcoplasmic reticulum (SR). Hinch et al. previously developed a tightly coupled LCC-RyR mathematical model, known as the Hinch model, that enables simulations to deal with a variety of functional states of whole-cell populations of a Ca(2+)-releasing unit using a personal computer. In this study, we developed a membrane excitation-contraction model of the human ventricular myocyte, which we call the human ventricular cell (HuVEC) model. This model is a hybrid of the most recent HuVEC models and the Hinch model. We modified the Hinch model to reproduce the regenerative activation and termination of CICR. In particular, we removed the inactivated RyR state and separated the single step of RyR activation by LCCs into triggering and regenerative steps. More importantly, we included the experimental measurement of a transient rise in Ca(2+) concentrations ([Ca(2+)], 10-15 μM) during CICR in the vicinity of Ca(2+)-releasing sites, and thereby calculated the effects of the local Ca(2+) gradient on CICR as well as membrane excitation. This HuVEC model successfully reconstructed both membrane excitation and key properties of CICR. The time course of CICR evoked by an action potential was accounted for by autonomous changes in an instantaneous equilibrium open probability of couplons. This autonomous time course was driven by a core feedback loop including the pivotal local [Ca(2+)], influenced by a time-dependent decay in the SR Ca(2+) content during CICR.

Computational modeling of anoctamin 1 calcium-activated chloride channels as pacemaker channels in interstitial cells of Cajal

Interstitial cells of Cajal (ICC) act as pacemaker cells in the gastrointestinal tract by generating electrical slow waves to regulate rhythmic smooth muscle contractions. Intrinsic Ca(2+) oscillations in ICC appear to produce the slow waves by activating pacemaker currents, currently thought to be carried by the Ca(2+)-activated Cl(-) channel anoctamin 1 (Ano1). In this article we present a novel model of small intestinal ICC pacemaker activity that incorporates store-operated Ca(2+) entry and a new model of Ano1 current. A series of simulations were carried out with the ICC model to investigate current controversies about the reversal potential of the Ano1 Cl(-) current in ICC and to predict the characteristics of the other ion channels that are necessary to generate slow waves. The model results show that Ano1 is a plausible pacemaker channel when coupled to a store-operated Ca(2+) channel but suggest that small cyclical depolarizations may still occur in ICC in Ano1 knockout mice. The results predict that voltage-dependent Ca(2+) current is likely to be negligible during the slow wave plateau phase. The model shows that the Cl(-) equilibrium potential is an important modulator of slow wave morphology, highlighting the need for a better understanding of Cl(-) dynamics in ICC.

Cardiac Na+-Ca2+ exchanger: dynamics of Ca2+-dependent activation and deactivation in intact myocytes

Cardiac Na(+)-Ca(2+) exchange (NCX) activity is regulated by [Ca(2+)]i. The physiological role and dynamics of this process in intact cardiomyocytes are largely unknown. We examined NCX Ca(2+) activation in intact rabbit and mouse cardiomyocytes at 37°C. Sarcoplasmic reticulum (SR) function was blocked, and cells were bathed in 2 mm Ca(2+). We probed Ca(2+) activation without voltage clamp by applying Na(+)-free (0 Na(+)) solution for 5 s bouts, repeated each 10 s, which should evoke [Ca(2+)]i transients due to Ca(2+) influx via NCX. In rested rabbit myocytes, Ca(2+) influx was undetectable even after 0 Na(+) applications were repeated for 2-5 min or more, suggesting that NCX was inactive. After external electric field stimulation pulses were applied, to admit Ca(2+) via L-type Ca(2+) channels, 0 Na(+) bouts activated Ca(2+) influx efficaciously, indicating that NCX had become active. Calcium activation increased with more field pulses, reaching a maximum typically after 15-20 pulses (1 Hz). At rest, NCX deactivated with a time constant typically of 20-40 s. An increase in [Na(+)]i, either in rabbit cardiomyocytes as a result of inhibition of Na(+)-K(+) pumping, or in mouse cardiomyocytes where normal [Na(+)]i is higher vs. rabbit, sensitized NCX to self-activation by 0 Na(+) bouts. In experiments with the SR functional but initially empty, the activation time course was slowed. It is possible that the SR initially accumulated Ca(2+) that would otherwise cause activation. We modelled Ca(2+) activation as a fourth-order highly co-operative process ([Ca]i required for half-activation K0.5act = 375 nm), with dynamics severalfold slower than the cardiac cycle. We incorporated this NCX model into an established ventricular myocyte model, which allowed us to predict responses to twitch stimulation in physiological conditions with the SR intact. Model NCX fractional activation increased from 0.1 to 1.0 as the frequency was increased from 0.2 to 2 Hz. By adjusting Ca(2+) activation on a multibeat time scale, NCX might better maintain a stable long-term Ca(2+) balance while contributing to the ability of myocytes to produce Ca(2+) transients over a wide range of intensity.