Yoga for the sinoatrial node: sarcoplasmic reticulum calcium release confers flexibility

Electrophysiological heterogeneity of pacemaker cells in the rabbit intercaval region, including the SA node: insights from recording multiple ion currents in each cell

Cardiac pacemaker cells, including cells of the sinoatrial node, are heterogeneous in size, morphology, and electrophysiological characteristics. The exact extent to which these cells differ electrophysiologically is unclear yet is critical to understanding their functioning. We examined major ionic currents in individual intercaval pacemaker cells (IPCs) sampled from the paracristal, intercaval region (including the sinoatrial node) that were spontaneously beating after enzymatic isolation from rabbit hearts. The beating rate was measured at baseline and after inhibition of the Ca2+ pump with cyclopiazonic acid. Thereafter, in each cell, we consecutively measured the density of funny current ( If), delayed rectifier K+ current ( IK) (a surrogate of repolarization capacity), and L-type Ca2+ current ( ICa,L) using whole cell patch clamp. The ionic current densities varied to a greater extent than previously appreciated, with some IPCs demonstrating very small or zero If . The density of none of the currents was correlated with cell size, while ICa,L and If densities were related to baseline beating rates. If density was correlated with IK density but not with that of ICa,L. Inhibition of Ca2+ cycling had a greater beating rate slowing effect in IPCs with lower If densities. Our numerical model simulation indicated that 1) IPCs with small (or zero) If or small ICa,L can operate via a major contribution of Ca2+ clock, 2) If-Ca2+-clock interplay could be important for robust pacemaking function, and 3) coupled If- IK function could regulate maximum diastolic potential. Thus, we have demonstrated marked electrophysiological heterogeneity of IPCs. This heterogeneity is manifested in basal beating rate and response to interference of Ca2+ cycling, which is linked to If. NEW & NOTEWORTHY In the present study, a hitherto unrecognized range of heterogeneity of ion currents in pacemaker cells from the intercaval region is demonstrated. Relationships between basal beating rate and L-type Ca2+ current and funny current ( If) density are uncovered, along with a positive relationship between If and delayed rectifier K+ current. Links are shown between the response to Ca2+ cycling blockade and If density.

Differential roles of two delayed rectifier potassium currents in regulation of ventricular action potential duration and arrhythmia susceptibility

KEY POINTS

Arrhythmias result from disruptions to cardiac electrical activity, although the factors that control cellular action potentials are incompletely understood. We combined mathematical modelling with experiments in heart cells from guinea pigs to determine how cellular electrical activity is regulated. A mismatch between modelling predictions and the experimental results allowed us to construct an improved, more predictive mathematical model. The balance between two particular potassium currents dictates how heart cells respond to perturbations and their susceptibility to arrhythmias.

ABSTRACT

Imbalances of ionic currents can destabilize the cardiac action potential and potentially trigger lethal cardiac arrhythmias. In the present study, we combined mathematical modelling with information-rich dynamic clamp experiments to determine the regulation of action potential morphology in guinea pig ventricular myocytes. Parameter sensitivity analysis was used to predict how changes in ionic currents alter action potential duration, and these were tested experimentally using dynamic clamp, a technique that allows for multiple perturbations to be tested in each cell. Surprisingly, we found that a leading mathematical model, developed with traditional approaches, systematically underestimated experimental responses to dynamic clamp perturbations. We then re-parameterized the model using a genetic algorithm, which allowed us to estimate ionic current levels in each of the cells studied. This unbiased model adjustment consistently predicted an increase in the rapid delayed rectifier K⁺ current and a drastic decrease in the slow delayed rectifier K⁺ current, and this prediction was validated experimentally. Subsequent simulations with the adjusted model generated the clinically relevant prediction that the slow delayed rectifier is better able to stabilize the action potential and suppress pro-arrhythmic events than the rapid delayed rectifier. In summary, iterative coupling of simulations and experiments enabled novel insight into how the balance between cardiac K⁺ currents influences ventricular arrhythmia susceptibility.

There and back again: Iterating between population-based modeling and experiments reveals surprising regulation of calcium transients in rat cardiac myocytes

While many ion channels and transporters involved in cardiac cellular physiology have been identified and described, the relative importance of each in determining emergent cellular behaviors remains unclear. Here we address this issue with a novel approach that combines population-based mathematical modeling with experimental tests to systematically quantify the relative contributions of different ion channels and transporters to the amplitude of the cellular Ca(2+) transient. Sensitivity analysis of a mathematical model of the rat ventricular cardiomyocyte quantified the response of cell behaviors to changes in the level of each ion channel and transporter, and experimental tests of these predictions were performed to validate or invalidate the predictions. The model analysis found that partial inhibition of the transient outward current in rat ventricular epicardial myocytes was predicted to have a greater impact on Ca(2+) transient amplitude than either: (1) inhibition of the same current in endocardial myocytes, or (2) comparable inhibition of the sarco/endoplasmic reticulum Ca(2+) ATPase (SERCA). Experimental tests confirmed the model predictions qualitatively but showed some quantitative disagreement. This guided us to recalibrate the model by adjusting the relative importance of several Ca(2+) fluxes, thereby improving the consistency with experimental data and producing a more predictive model. Analysis of human cardiomyocyte models suggests that the relative importance of outward currents to Ca(2+) transporters is generalizable to human atrial cardiomyocytes, but not ventricular cardiomyocytes. Overall, our novel approach of combining population-based mathematical modeling with experimental tests has yielded new insight into the relative importance of different determinants of cell behavior.

Automatic parameter estimation of multicompartmental neuron models via minimization of trace error with control adjustment

We describe a new technique to fit conductance-based neuron models to intracellular voltage traces from isolated biological neurons. The biological neurons are recorded in current-clamp with pink (1/f) noise injected to perturb the activity of the neuron. The new algorithm finds a set of parameters that allows a multicompartmental model neuron to match the recorded voltage trace. Attempting to match a recorded voltage trace directly has a well-known problem: mismatch in the timing of action potentials between biological and model neuron is inevitable and results in poor phenomenological match between the model and data. Our approach avoids this by applying a weak control adjustment to the model to promote alignment during the fitting procedure. This approach is closely related to the control theoretic concept of a Luenberger observer. We tested this approach on synthetic data and on data recorded from an anterior gastric receptor neuron from the stomatogastric ganglion of the crab Cancer borealis. To test the flexibility of this approach, the synthetic data were constructed with conductance models that were different from the ones used in the fitting model. For both synthetic and biological data, the resultant models had good spike-timing accuracy.

Modern perspectives on numerical modeling of cardiac pacemaker cell

Cardiac pacemaking is a complex phenomenon that is still not completely understood. Together with experimental studies, numerical modeling has been traditionally used to acquire mechanistic insights in this research area. This review summarizes the present state of numerical modeling of the cardiac pacemaker, including approaches to resolve present paradoxes and controversies. Specifically we discuss the requirement for realistic modeling to consider symmetrical importance of both intracellular and cell membrane processes (within a recent "coupled-clock" theory). Promising future developments of the complex pacemaker system models include the introduction of local calcium control, mitochondria function, and biochemical regulation of protein phosphorylation and cAMP production. Modern numerical and theoretical methods such as multi-parameter sensitivity analyses within extended populations of models and bifurcation analyses are also important for the definition of the most realistic parameters that describe a robust, yet simultaneously flexible operation of the coupled-clock pacemaker cell system. The systems approach to exploring cardiac pacemaker function will guide development of new therapies such as biological pacemakers for treating insufficient cardiac pacemaker function that becomes especially prevalent with advancing age.

Comprehensive analyses of ventricular myocyte models identify targets exhibiting favorable rate dependence

Reverse rate dependence is a problematic property of antiarrhythmic drugs that prolong the cardiac action potential (AP). The prolongation caused by reverse rate dependent agents is greater at slow heart rates, resulting in both reduced arrhythmia suppression at fast rates and increased arrhythmia risk at slow rates. The opposite property, forward rate dependence, would theoretically overcome these parallel problems, yet forward rate dependent (FRD) antiarrhythmics remain elusive. Moreover, there is evidence that reverse rate dependence is an intrinsic property of perturbations to the AP. We have addressed the possibility of forward rate dependence by performing a comprehensive analysis of 13 ventricular myocyte models. By simulating populations of myocytes with varying properties and analyzing population results statistically, we simultaneously predicted the rate-dependent effects of changes in multiple model parameters. An average of 40 parameters were tested in each model, and effects on AP duration were assessed at slow (0.2 Hz) and fast (2 Hz) rates. The analysis identified a variety of FRD ionic current perturbations and generated specific predictions regarding their mechanisms. For instance, an increase in L-type calcium current is FRD when this is accompanied by indirect, rate-dependent changes in slow delayed rectifier potassium current. A comparison of predictions across models identified inward rectifier potassium current and the sodium-potassium pump as the two targets most likely to produce FRD AP prolongation. Finally, a statistical analysis of results from the 13 models demonstrated that models displaying minimal rate-dependent changes in AP shape have little capacity for FRD perturbations, whereas models with large shape changes have considerable FRD potential. This can explain differences between species and between ventricular cell types. Overall, this study provides new insights, both specific and general, into the determinants of AP duration rate dependence, and illustrates a strategy for the design of potentially beneficial antiarrhythmic drugs.

Several drugs intended to treat cardiac arrhythmias have failed because of unfavorable rate-dependent properties. That is, the drugs fail to alter electrical activity at fast heart rates, where this would be beneficial, but they do affect electrical activity at slow rates, where this is unwanted. In targeted studies, several agents have been shown to exhibit these unfavorable properties, suggesting that these rate-dependent responses may be intrinsic to ventricular muscle. To determine whether drugs with desirable rate-dependent properties could be rationally designed, we performed comprehensive and systematic analyses of several heart cell models. These analyses calculated the rate-dependent properties of changes in any model parameter, thereby generating simultaneously a large number of model predictions. The analyses showed that targets with favorable rate-dependent properties could indeed be identified, and further simulations uncovered the mechanisms underlying these behaviors. Moreover, a quantitative comparison of results obtained in different models provided new insight in why a given drug applied to different species, or to different tissue types, might produce different rate-dependent behaviors. Overall this study shows how a comprehensive and systematic approach to heart cell models can both identify novel targets and produce more general insight into rate-dependent alterations to cardiac electrical activity.