Clusters of calcium release channels harness the Ising phase transition to confine their elementary intracellular signals

Elementary intracellular Ca signals approximated as a transition of release channel system from a metastable state

Cardiac muscle contraction is initiated by an elementary Ca signal (called Ca spark) which is achieved by collective action of Ca release channels in a cluster. The mechanism of this synchronization remains uncertain. We approached Ca spark activation as an emergent phenomenon of an interactive system of release channels. We constructed a weakly lumped Markov chain that applies an Ising model formalism to such release channel clusters and probable open channel configurations and demonstrated that spark activation is described as a system transition from a metastable to an absorbing state, analogous to the pressure required to overcome surface tension in bubble formation. This yielded quantitative estimates of the spark generation probability as a function of various system parameters. We performed numerical simulations to find spark probabilities as a function of sarcoplasmic reticulum Ca concentration, obtaining similar values for spark activation threshold as our analytic model, as well as those reported in experimental studies. Our parametric sensitivity analyses also showed that the spark activation threshold decreased as Ca sensitivity of RyR activation and RyR cluster size increased.

Modeling Calcium Cycling in the Heart: Progress, Pitfalls, and Challenges

Intracellular calcium (Ca) cycling in the heart plays key roles in excitation-contraction coupling and arrhythmogenesis. In cardiac myocytes, the Ca release channels, i.e., the ryanodine receptors (RyRs), are clustered in the sarcoplasmic reticulum membrane, forming Ca release units (CRUs). The RyRs in a CRU act collectively to give rise to discrete Ca release events, called Ca sparks. A cell contains hundreds to thousands of CRUs, diffusively coupled via Ca to form a CRU network. A rich spectrum of spatiotemporal Ca dynamics is observed in cardiac myocytes, including Ca sparks, spark clusters, mini-waves, persistent whole-cell waves, and oscillations. Models of different temporal and spatial scales have been developed to investigate these dynamics. Due to the complexities of the CRU network and the spatiotemporal Ca dynamics, it is challenging to model the Ca cycling dynamics in the cardiac system, particularly at the tissue sales. In this article, we review the progress of modeling of Ca cycling in cardiac systems from single RyRs to the tissue scale, the pros and cons of the current models and different modeling approaches, and the challenges to be tackled in the future.

Disorder in Ca2+ release unit locations confers robustness but cuts flexibility of heart pacemaking

Excitation-contraction coupling kinetics is dictated by the action potential rate of sinoatrial-nodal cells. These cells generate local Ca releases (LCRs) that activate Na/Ca exchanger current, which accelerates diastolic depolarization and determines the pace. LCRs are generated by clusters of ryanodine receptors, Ca release units (CRUs), residing in the sarcoplasmic reticulum. While CRU distribution exhibits substantial heterogeneity, its functional importance remains unknown. Using numerical modeling, here we show that with a square lattice distribution of CRUs, Ca-induced-Ca-release propagation during diastolic depolarization is insufficient for pacemaking within a broad range of realistic ICaL densities. Allowing each CRU to deviate randomly from its lattice position allows sparks to propagate, as observed experimentally. As disorder increases, the CRU distribution exhibits larger empty spaces and simultaneously CRU clusters, as in Poisson clumping. Propagating within the clusters, Ca release becomes synchronized, increasing action potential rate and reviving pacemaker function of dormant/nonfiring cells. However, cells with fully disordered CRU positions could not reach low firing rates and their β-adrenergic-receptor stimulation effect was substantially decreased. Inclusion of Cav1.3, a low-voltage activation L-type Ca channel isoform into ICaL, strongly increases recruitment of CRUs to fire during diastolic depolarization, increasing robustness of pacemaking and complementing effects of CRU distribution. Thus, order/disorder in CRU locations along with Cav1.3 expression regulates pacemaker function via synchronization of CRU firing. Excessive CRU disorder and/or overexpression of Cav1.3 boosts pacemaker function in the basal state, but limits the rate range, which may contribute to heart rate range decline with age and disease.

Functional Heterogeneity of Cell Populations Increases Robustness of Pacemaker Function in a Numerical Model of the Sinoatrial Node Tissue

Each heartbeat is initiated by specialized pacemaker cells operating within the sinoatrial node (SAN). While individual cells within SAN tissue exhibit substantial heterogeneity of their electrophysiological parameters and Ca cycling, the role of this heterogeneity for cardiac pacemaker function remains mainly unknown. Here we investigated the problem numerically in a 25 × 25 square grid of connected coupled-clock Maltsev-Lakatta cell models. The tissue models were populated by cells with different degree of heterogeneity of the two key model parameters, maximum L-type Ca current conductance (g_CaL) and sarcoplasmic reticulum Ca pumping rate (P_up). Our simulations showed that in the areas of P_up-g_CaL parametric space at the edge of the system stability, where action potential (AP) firing is absent or dysrhythmic in SAN tissue models populated with identical cells, rhythmic AP firing can be rescued by populating the tissues with heterogeneous cells. This robust SAN function is synergistic with respect to heterogeneity in g_CaL and P_up and can be further strengthened by clustering of cells with similar properties. The effect of cell heterogeneity is not due to a simple summation of activity of intrinsically firing cells naturally present in heterogeneous SAN; rather AP firing cells locally and critically interact with non-firing/dormant cells. When firing cells prevail, they recruit many dormant cells to fire, strongly enhancing overall SAN function; and vice versa, prevailing dormant cells suppress AP firing in cells with intrinsic automaticity and halt SAN function. The transitions between firing and non-firing states of the system are sharp, resembling phase transitions in statistical physics. Furthermore, robust function of heterogeneous SAN tissue requires weak cell coupling, a known property of the central area of SAN where cardiac impulse emerges; stronger cell coupling reduces AP firing rate and ultimately halts SAN automaticity at the edge of stability.

Collective Dynamics of Ion Channels on Bilayer Lipid Membranes

Ion channels self-organize on cellular and organelle membranes as clusters and mutually modulate their gating behavior. It has been reported that the efficient information transfer is achieved by cooperative clustering of ion channels. To address the origin and nature of collective dynamics in ion channel clusters, a statistical mechanical model, namely, the Zimm-Bragg-type model in two dimensions with unequal weight distribution in channel-channel interactions, has been proposed. Nearest neighbor interaction along with next-nearest neighbor interaction has been considered, assuming symmetric spatial organization. The multichannel bilayer electrophysiology recordings of the voltage-dependent anion channel (VDAC) from rat brain mitochondria have been analyzed in order to test and further extend the model. The model successfully describes the multichannel gating behavior and self-organization of the VDAC cluster.

Simulating cardiac Ca2+ release units: effects of RyR cluster size and Ca2+ buffers on diastolic Ca2+ leak

Leak of Ca²⁺ out of the cardiac sarcoplasmic reticulum (SR) via ryanodine receptors (RyRs) during diastole is vital to regulate SR Ca²⁺ levels. This leak can become deleterious when large spontaneous RyR-mediated Ca²⁺ release events evoke proarrhythmic Ca²⁺ waves that can lead to delayed after-depolarizations. Here, we model diastolic SR Ca²⁺ leak at individual SR Ca²⁺ release sites using computer simulations of RyR arrays like those in the dyadic cleft. The results show that RyR arrays size has a significant effect on SR Ca²⁺ leak, with bigger arrays producing larger and more frequent Ca²⁺ release events. Moreover, big RyR arrays are more susceptible to small changes in the levels of dyadic Ca²⁺ buffers. Such changes in buffering shift Ca²⁺ leak from small Ca²⁺ release events (involving few open RyRs) to larger events (with many open RyRs). Moreover, by analyzing a large parameter space of possible buffering and SR Ca²⁺ loads, we find further evidence for the hypothesis that SR Ca²⁺ leak by RyR arrays can undergo a sudden phase transition.

Social heterogeneity drives complex patterns of the COVID-19 pandemic: insights from a novel Stochastic Heterogeneous Epidemic Model (SHEM)

In todays absence of a vaccine and impactful treatments, the most effective way to combat the virus is to find and implement mitigation strategies. An invaluable resource in this task is numerical modeling that can reveal key factors in COVID-19 pandemic development. On the other hand, it has become evident that regional infection curves of COVID-19 exhibit complex patterns which often differ from curves predicted by forecasting models. The wide variations in attack rate observed among different social strata suggest that this may be due to social heterogeneity not accounted for by regional models. We investigated this hypothesis by developing and using a new Stochastic Heterogeneous Epidemic Model (SHEM) that focuses on subpopulations that are vulnerable in the sense of having an increased likelihood of spreading infection among themselves. We found that the isolation or embedding of vulnerable sub-clusters in a major population hub generated complex stochastic infection patterns which included multiple peaks and growth periods, an extended plateau, a prolonged tail, or a delayed second wave of infection. Embedded vulnerable groups became hotspots that drove infection despite efforts of the main population to socially distance, while isolated groups suffered delayed but intense infection. Amplification of infection by these hotspots facilitated transmission from one urban area to another, causing the epidemic to hopscotch in a stochastic manner to places it would not otherwise reach, resembling a microcosm of the situation worldwide as of September 2020. Our results suggest that social heterogeneity is a key factor in the formation of complex infection propagation patterns. Thus, the mitigation of vulnerable groups is essential to control the COVID-19 pandemic worldwide. The design of our new model allows it to be applied in future studies of real-world scenarios on any scale, limited only by computing memory and the ability to determine the underlying topology and parameters.

Recruiting RyRs to Open in a Ca2+ Release Unit: Single-RyR Gating Properties Make RyR Group Dynamics

In cardiac myocytes, clusters of type-2 ryanodine receptors (RyR2s) release Ca²⁺ from the sarcoplasmic reticulum (SR) via a positive feedback mechanism in which fluxed Ca²⁺ activates nearby RyRs. Although the general principles of this are understood, less is known about how single-RyR gating properties define the RyR group dynamics in an array of many channels. Here, we examine this using simulations with three models of RyR gating that have identical open probabilities: the commonly used two-state Markov gating model, one that utilizes multiple exponentials to fit single-channel open time (OT) and closed time (CT) distributions, and an extension of this multiexponential model that also includes experimentally measured correlations between single-channel OTs and CTs. The simulations of RyR clusters that utilize the multiexponential gating model produce infrequent Ca²⁺ release events with relatively few open RyRs. Ca²⁺ release events become even smaller when OT/CT correlations are included. This occurs because the correlations produce a small but consistent bias against recruiting more RyRs to open during the middle of a Ca²⁺ release event, between the initiation and termination phases (which are unaltered compared to the uncorrelated simulations). In comparison, the two-state model produces frequent, large, and long Ca²⁺ release events because it had a recruitment bias in favor of opening more RyRs. This difference stems from the two-state model's single-RyR OT and CT distributions being qualitatively different from the experimental ones. Thus, the details of single-RyR gating can profoundly affect SR Ca²⁺ release even if open probability and mean OTs and CTs are identical. We also show that Ca²⁺ release events can terminate spontaneously without any reduction in SR [Ca²⁺], luminal regulation, Ca²⁺-dependent inactivation, or physical coupling between RyRs when Ca²⁺ flux is below a threshold value. This supports and extends the pernicious attrition/induction decay hypothesis that SR Ca²⁺ release events terminate below a threshold Ca²⁺ flux.

Mechanisms of Calcium Leak from Cardiac Sarcoplasmic Reticulum Revealed by Statistical Mechanics

Heart muscle contraction is normally activated by a synchronized Ca release from sarcoplasmic reticulum (SR), a major intracellular Ca store. However, under abnormal conditions, Ca leaks from the SR, decreasing heart contraction amplitude and increasing risk of life-threatening arrhythmia. The mechanisms and regimes of SR operation generating the abnormal Ca leak remain unclear. Here, we employed both numerical and analytical modeling to get mechanistic insights into the emergent Ca leak phenomenon. Our numerical simulations using a detailed realistic model of the Ca release unit reveal sharp transitions resulting in Ca leak. The emergence of leak is closely mapped mathematically to the Ising model from statistical mechanics. The system steady-state behavior is determined by two aggregate parameters: the analogs of magnetic field (h) and the inverse temperature (β) in the Ising model, for which we have explicit formulas in terms of SR [Ca] and release channel opening and closing rates. The classification of leak regimes takes the shape of a phase β-h diagram, with the regime boundaries occurring at h = 0 and a critical value of β (β^∗) that we estimate using a classical Ising model and mean field theory. Our theory predicts that a synchronized Ca leak will occur when h > 0 and β >β^∗, and a disordered leak occurs when β <β^∗ and h is not too negative. The disorder leak is distinguished from synchronized leak (in long-lasting sparks) by larger Peierls contour lengths, an output parameter reflecting degree of disorder. Thus, in addition to our detailed numerical model approach, we also offer an instantaneous computational tool using analytical formulas of the Ising model for respective ryanodine receptor parameters and SR Ca load that describe and classify phase transitions and leak emergence.

Progress in Mathematical Modeling of Gastrointestinal Slow Wave Abnormalities

Gastrointestinal (GI) motility is regulated in part by electrophysiological events called slow waves, which are generated by the interstitial cells of Cajal (ICC). Slow waves propagate by a process of "entrainment," which occurs over a decreasing gradient of intrinsic frequencies in the antegrade direction across much of the GI tract. Abnormal initiation and conduction of slow waves have been demonstrated in, and linked to, a number of GI motility disorders. A range of mathematical models have been developed to study abnormal slow waves and applied to propose novel methods for non-invasive detection and therapy. This review provides a general outline of GI slow wave abnormalities and their recent classification using multi-electrode (high-resolution) mapping methods, with a particular emphasis on the spatial patterns of these abnormal activities. The recently-developed mathematical models are introduced in order of their biophysical scale from cellular to whole-organ levels. The modeling techniques, main findings from the simulations, and potential future directions arising from notable studies are discussed.

New Insights in Cardiac Calcium Handling and Excitation-Contraction Coupling

Excitation-contraction (EC) coupling denotes the conversion of electric stimulus in mechanic output in contractile cells. Several studies have demonstrated that calcium (Ca²⁺) plays a pivotal role in this process. Here we present a comprehensive and updated description of the main systems involved in cardiac Ca²⁺ handling that ensure a functional EC coupling and their pathological alterations, mainly related to heart failure.