Stochastic simulation of cardiac ventricular myocyte calcium dynamics and waves

Computational cardiology and risk stratification for sudden cardiac death: one of the grand challenges for cardiology in the 21st century

Risk stratification in the context of sudden cardiac death has been acknowledged as one of the major challenges facing cardiology for the past four decades. In recent years, the advent of high performance computing has facilitated organ-level simulation of the heart, meaning we can now examine the causes, mechanisms and impact of cardiac dysfunction in silico. As a result, computational cardiology, largely driven by the Physiome project, now stands at the threshold of clinical utility in regards to risk stratification and treatment of patients at risk of sudden cardiac death. In this white paper, we outline a roadmap of what needs to be done to make this translational step, using the relatively well-developed case of acquired or drug-induced long QT syndrome as an exemplar case.

Calcium Sparks in the Heart: Dynamics and Regulation

Calcium (Ca²⁺) plays a central role in the contraction of the heart. It is the bi-directional link between electrical excitation of the heart and contraction. Electrical excitation initiates Ca²⁺influx across the sarcolemma and T-tubular membrane that triggered calcium release from the sarcoplasmic reticulum. Ca²⁺sparks are the elementary events of calcium release from the sarcoplasmic reticulum. Therefore, understanding the dynamics of Ca²⁺sparks is essential for understanding the function of the heart. To this end, numerous experimental and computational studies have focused on this topic, exploring the mechanisms of calcium spark initiation, termination, and regulation and what role these play in normal and patho-physiology. The proper understanding of Ca²⁺ spark regulation and dynamics serves as the foundation for our insights into a multitude of pathological conditions may develop that can be the result of structural and/or functional changes at the cellular or subcellular level. Computational modeling of Ca²⁺ spark dynamics has proven to be a useful tool to understand Ca²⁺ spark dynamics. This review addresses our current understanding of Ca²⁺ sparks and how synchronized SR Ca²⁺ release, in which Ca²⁺ sparks is a major pathway, is linked to the different cardiac diseases, especially arrhythmias.

Hierarchical clustering of ryanodine receptors enables emergence of a calcium clock in sinoatrial node cells

The sinoatrial node, whose cells (sinoatrial node cells [SANCs]) generate rhythmic action potentials, is the primary pacemaker of the heart. During diastole, calcium released from the sarcoplasmic reticulum (SR) via ryanodine receptors (RyRs) interacts with membrane currents to control the rate of the heartbeat. This "calcium clock" takes the form of stochastic, partially periodic, localized calcium release (LCR) events that propagate, wave-like, for limited distances. The detailed mechanisms controlling the calcium clock are not understood. We constructed a computational model of SANCs, including three-dimensional diffusion and buffering of calcium in the cytosol and SR; explicit, stochastic gating of individual RyRs and L-type calcium channels; and a full complement of voltage- and calcium-dependent membrane currents. We did not include an anatomical submembrane space or inactivation of RyRs, the two heuristic components that have been used in prior models but are not observed experimentally. When RyRs were distributed in discrete clusters separated by >1 µm, only isolated sparks were produced in this model and LCR events did not form. However, immunofluorescent staining of SANCs for RyR revealed the presence of bridging RyR groups between large clusters, forming an irregular network. Incorporation of this architecture into the model led to the generation of propagating LCR events. Partial periodicity emerged from the interaction of LCR events, as observed experimentally. This calcium clock becomes entrained with membrane currents to accelerate the beating rate, which therefore was controlled by the activity of the SERCA pump, RyR sensitivity, and L-type current amplitude, all of which are targets of β-adrenergic-mediated phosphorylation. Unexpectedly, simulations revealed the existence of a pathological mode at high RyR sensitivity to calcium, in which the calcium clock loses synchronization with the membrane, resulting in a paradoxical decrease in beating rate in response to β-adrenergic stimulation. The model indicates that the hierarchical clustering of surface RyRs in SANCs may be a crucial adaptive mechanism. Pathological desynchronization of the clocks may explain sinus node dysfunction in heart failure and RyR mutations.