Gating kinetics and pharmacological properties of small-conductance Ca2+-activated potassium channels

Small-conductance (SK) calcium-activated potassium channels are a promising treatment target in atrial fibrillation. However, the functional properties that differentiate SK inhibitors remain poorly understood. The objective of this study was to determine how two unrelated SK channel inhibitors, apamin and AP14145, impact SK channel function in excised inside-out single-channel recordings. Surprisingly, both apamin and AP14145 exert much of their inhibition by inducing a class of very-long-lived channel closures (apamin: τc,vl = 11.8 ± 7.1 s, and AP14145: τc,vl = 10.3 ± 7.2 s), which were never observed under control conditions. Both inhibitors also induced changes to the three closed and two open durations typical of normal SK channel gating. AP14145 shifted the open duration distribution to favor longer open durations, whereas apamin did not alter open-state kinetics. AP14145 also prolonged the two shortest channel closed durations (AP14145: τc,s = 3.50 ± 0.81 ms, and τc,i = 32.0 ± 6.76 ms versus control: τc,s = 1.59 ± 0.19 ms, and τc,i = 13.5 ± 1.17 ms), thus slowing overall gating kinetics within bursts of channel activity. In contrast, apamin accelerated intraburst gating kinetics by shortening the two shortest closed durations (τc,s = 0.75 ± 0.10 ms and τc,i = 5.08 ± 0.49 ms) and inducing periods of flickery activity. Finally, AP14145 introduced a unique form of inhibition by decreasing unitary current amplitude. SK channels exhibited two clearly distinguishable amplitudes (control: Ahigh = 0.76 ± 0.03 pA, and Alow = 0.54 ± 0.03 pA). AP14145 both reduced the fraction of patches exhibiting the higher amplitude (AP14145: 4/9 patches versus control: 16/16 patches) and reduced the mean low amplitude (0.38 ± 0.03 pA). Here, we have demonstrated that both inhibitors introduce very long channel closures but that each also exhibits unique effects on other components of SK gating kinetics and unitary current. The combination of these effects is likely to be critical for understanding the functional differences of each inhibitor in the context of cyclical Ca2+-dependent channel activation in vivo.

Mechanisms of spontaneous Ca2+ release-mediated arrhythmia in a novel 3D human atrial myocyte model: II Ca2+ -handling protein variation

Disruption of the transverse-axial tubule system (TATS) in diseases such as heart failure and atrial fibrillation occurs in combination with changes in the expression and distribution of key Ca²⁺ -handling proteins. Together this ultrastructural and ionic remodelling is associated with aberrant Ca²⁺ cycling and electrophysiological instabilities that underlie arrhythmic activity. However, due to the concurrent changes in TATs and Ca²⁺ -handling protein expression and localization that occur in disease it is difficult to distinguish their individual contributions to the arrhythmogenic state. To investigate this, we applied our novel 3D human atrial myocyte model with spatially detailed Ca²⁺ diffusion and TATS to investigate the isolated and interactive effects of changes in expression and localization of key Ca²⁺ -handling proteins and variable TATS density on Ca²⁺ -handling abnormality driven membrane instabilities. We show that modulating the expression and distribution of the sodium-calcium exchanger, ryanodine receptors and the sarcoplasmic reticulum (SR) Ca²⁺ buffer calsequestrin have varying pro- and anti-arrhythmic effects depending on the balance of opposing influences on SR Ca²⁺ leak-load and Ca²⁺ -voltage relationships. Interestingly, the impact of protein remodelling on Ca²⁺ -driven proarrhythmic behaviour varied dramatically depending on TATS density, with intermediately tubulated cells being more severely affected compared to detubulated and densely tubulated myocytes. This work provides novel mechanistic insight into the distinct and interactive consequences of TATS and Ca²⁺ -handling protein remodelling that underlies dysfunctional Ca²⁺ cycling and electrophysiological instability in disease. KEY POINTS: In our companion paper we developed a 3D human atrial myocyte model, coupling electrophysiology and Ca²⁺ handling with subcellular spatial details governed by the transverse-axial tubule system (TATS). Here we utilize this model to mechanistically examine the impact of TATS loss and changes in the expression and distribution of key Ca²⁺ -handling proteins known to be remodelled in disease on Ca²⁺ homeostasis and electrophysiological stability. We demonstrate that varying the expression and localization of these proteins has variable pro- and anti-arrhythmic effects with outcomes displaying dependence on TATS density. Whereas detubulated myocytes typically appear unaffected and densely tubulated cells seem protected, the arrhythmogenic effects of Ca²⁺ handling protein remodelling are profound in intermediately tubulated cells. Our work shows the interaction between TATS and Ca²⁺ -handling protein remodelling that underlies the Ca²⁺ -driven proarrhythmic behaviour observed in atrial fibrillation and may help to predict the effects of antiarrhythmic strategies at varying stages of ultrastructural remodelling.

Mechanisms of spontaneous Ca2+ release-mediated arrhythmia in a novel 3D human atrial myocyte model: I. Transverse-axial tubule variation

Intracellular calcium (Ca²⁺ ) cycling is tightly regulated in the healthy heart ensuring effective contraction. This is achieved by transverse (t)-tubule membrane invaginations that facilitate close coupling of key Ca²⁺ -handling proteins such as the L-type Ca²⁺ channel and Na⁺ -Ca²⁺ exchanger (NCX) on the cell surface with ryanodine receptors (RyRs) on the intracellular Ca²⁺ store. Although less abundant and regular than in the ventricle, t-tubules also exist in atrial myocytes as a network of transverse invaginations with axial extensions known as the transverse-axial tubule system (TATS). In heart failure and atrial fibrillation, there is TATS remodelling that is associated with aberrant Ca²⁺ -handling and Ca²⁺ -induced arrhythmic activity; however, the mechanism underlying this is not fully understood. To address this, we developed a novel 3D human atrial myocyte model that couples electrophysiology and Ca²⁺ -handling with variable TATS organization and density. We extensively parameterized and validated our model against experimental data to build a robust tool examining TATS regulation of subcellular Ca²⁺ release. We found that varying TATS density and thus the localization of key Ca²⁺ -handling proteins has profound effects on Ca²⁺ handling. Following TATS loss, there is reduced NCX that results in increased cleft Ca²⁺ concentration through decreased Ca²⁺ extrusion. This elevated Ca²⁺ increases RyR open probability causing spontaneous Ca²⁺ releases and the promotion of arrhythmogenic waves (especially in the cell interior) leading to voltage instabilities through delayed afterdepolarizations. In summary, the present study demonstrates a mechanistic link between TATS remodelling and Ca²⁺ -driven proarrhythmic behaviour that probably reflects the arrhythmogenic state observed in disease. KEY POINTS: Transverse-axial tubule systems (TATS) modulate Ca²⁺ handling and excitation-contraction coupling in atrial myocytes, with TATS remodelling in heart failure and atrial fibrillation being associated with altered Ca²⁺ cycling and subsequent arrhythmogenesis. To investigate the poorly understood mechanisms linking TATS variation and spontaneous Ca²⁺ release, we built, parameterized and validated a 3D human atrial myocyte model coupling electrophysiology and spatially-detailed subcellular Ca²⁺ handling governed by the TATS. Simulated TATS loss causes diastolic Ca²⁺ and voltage instabilities through reduced Na⁺ -Ca²⁺ exchanger-mediated Ca²⁺ removal, cleft Ca²⁺ accumulation and increased ryanodine receptor open probability, resulting in spontaneous Ca²⁺ release and promotion of arrhythmogenic waves and delayed afterdepolarizations. At fast electrical rates typical of atrial tachycardia/fibrillation, spontaneous Ca²⁺ releases are larger and more frequent in the cell interior than at the periphery. Our work provides mechanistic insight into how atrial TATS remodelling can lead to Ca²⁺ -driven instabilities that may ultimately contribute to the arrhythmogenic state in disease.

Prolonged β-adrenergic stimulation disperses ryanodine receptor clusters in cardiomyocytes

Ryanodine receptors (RyRs) exhibit dynamic arrangements in cardiomyocytes, and we previously showed that ‘dispersion’ of RyR clusters disrupts Ca²⁺ homeostasis during heart failure (HF) (Kolstad et al., eLife, 2018). Here, we investigated whether prolonged β-adrenergic stimulation, a hallmark of HF, promotes RyR cluster dispersion and examined the underlying mechanisms. We observed that treatment of healthy rat cardiomyocytes with isoproterenol for 1 hr triggered progressive fragmentation of RyR clusters. Pharmacological inhibition of Ca²⁺/calmodulin-dependent protein kinase II (CaMKII) reversed these effects, while cluster dispersion was reproduced by specific activation of CaMKII, and in mice with constitutively active Ser2814-RyR. A similar role of protein kinase A (PKA) in promoting RyR cluster fragmentation was established by employing PKA activation or inhibition. Progressive cluster dispersion was linked to declining Ca²⁺ spark fidelity and magnitude, and slowed release kinetics from Ca²⁺ propagation between more numerous RyR clusters. In healthy cells, this served to dampen the stimulatory actions of β-adrenergic stimulation over the longer term and protect against pro-arrhythmic Ca²⁺ waves. However, during HF, RyR dispersion was linked to impaired Ca²⁺ release. Thus, RyR localization and function are intimately linked via channel phosphorylation by both CaMKII and PKA, which, while finely tuned in healthy cardiomyocytes, underlies impaired cardiac function during pathology.

Validating the Arrhythmogenic Potential of High-, Intermediate-, and Low-Risk Drugs in a Human-Induced Pluripotent Stem Cell-Derived Cardiac Microphysiological System

Evaluation of arrhythmogenic drugs is required by regulatory agencies before any new compound can obtain market approval. Despite rigorous review, cardiac disorders remain the second most common cause for safety-related market withdrawal. On the other hand, false-positive preclinical findings prohibit potentially beneficial candidates from moving forward in the development pipeline. Complex in vitro models using cardiomyocytes derived from human-induced pluripotent stem cells (hiPSC-CM) have been identified as a useful tool that allows for rapid and cost-efficient screening of proarrhythmic drug risk. Currently available hiPSC-CM models employ simple two-dimensional (2D) culture formats with limited structural and functional relevance to the human heart muscle. Here, we present the use of our 3D cardiac microphysiological system (MPS), composed of a hiPSC-derived heart micromuscle, as a platform for arrhythmia risk assessment. We employed two different hiPSC lines and tested seven drugs with known ion channel effects and known clinical risk: dofetilide and bepridil (high risk); amiodarone and terfenadine (intermediate risk); and nifedipine, mexiletine, and lidocaine (low risk). The cardiac MPS successfully predicted drug cardiotoxicity risks based on changes in action potential duration, beat waveform (i.e., shape), and occurrence of proarrhythmic events of healthy patient hiPSC lines in the absence of risk cofactors. We showcase examples where the cardiac MPS outperformed existing hiPSC-CM 2D models.

Arrhythmogenic influence of mutations in a myocyte-based computational model of the pulmonary vein sleeve

In the heart, electrophysiological dysregulation arises from defects at many biological levels (from point mutations in ion channel proteins to gross structural abnormalities). These defects disrupt the normal pattern of electrical activation, producing ectopic activity and reentrant arrhythmia. To interrogate mechanisms that link these primary biological defects to macroscopic electrophysiologic dysregulation most prior computational studies have utilized either (i) detailed models of myocyte ion channel dynamics at limited spatial scales, or (ii) homogenized models of action potential conduction that reproduce arrhythmic activity at tissue and organ levels. Here we apply our recent model (EMI), which integrates electrical activation and propagation across these scales, to study human atrial arrhythmias originating in the pulmonary vein (PV) sleeves. These small structures initiate most supraventricular arrhythmias and include pronounced myocyte-to-myocyte heterogeneities in ion channel expression and intercellular coupling. To test EMI's cell-based architecture in this physiological context we asked whether ion channel mutations known to underlie atrial fibrillation are capable of initiating arrhythmogenic behavior via increased excitability or reentry in a schematic PV sleeve geometry. Our results illustrate that EMI's improved spatial resolution can directly interrogate how electrophysiological changes at the individual myocyte level manifest in tissue and as arrhythmia in the PV sleeve.

Image-Driven Modeling of Nanoscopic Cardiac Function: Where Have We Come From, and Where Are We Going?

Complementary developments in microscopy and mathematical modeling have been critical to our understanding of cardiac excitation–contraction coupling. Historically, limitations imposed by the spatial or temporal resolution of imaging methods have been addressed through careful mathematical interrogation. Similarly, limitations imposed by computational power have been addressed by imaging macroscopic function in large subcellular domains or in whole myocytes. As both imaging resolution and computational tractability have improved, the two approaches have nearly merged in terms of the scales that they can each be used to interrogate. With this review we will provide an overview of these advances and their contribution to understanding ventricular myocyte function, including exciting developments over the last decade. We specifically focus on experimental methods that have pushed back limits of either spatial or temporal resolution of nanoscale imaging (e.g., DNA-PAINT), or have permitted high resolution imaging on large cellular volumes (e.g., serial scanning electron microscopy). We also review the progression of computational approaches used to integrate and interrogate these new experimental data sources, and comment on near-term advances that may unify understanding of the underlying biology. Finally, we comment on several outstanding questions in cardiac physiology that stand to benefit from a concerted and complementary application of these new experimental and computational methods.

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.

Sarcoplasmic Reticulum Calcium Release Is Required for Arrhythmogenesis in the Mouse

Dysfunctional sarcoplasmic reticulum Ca²⁺ handling is commonly observed in heart failure, and thought to contribute to arrhythmogenesis through several mechanisms. Some time ago we developed a cardiomyocyte-specific inducible SERCA2 knockout mouse, which is remarkable in the degree to which major adaptations to sarcolemmal Ca²⁺ entry and efflux overcome the deficit in SR reuptake to permit relatively normal contractile function. Conventionally, those adaptations would also be expected to dramatically increase arrhythmia susceptibility. However, that susceptibility has never been tested, and it is possible that the very rapid repolarization of the murine action potential (AP) allows for large changes in sarcolemmal Ca²⁺ transport without substantially disrupting electrophysiologic stability. We investigated this hypothesis through telemetric ECG recording in the SERCA2-KO mouse, and patch-clamp electrophysiology, Ca²⁺ imaging, and mathematical modeling of isolated SERCA2-KO myocytes. While the SERCA2-KO animals exhibit major (and unique) electrophysiologic adaptations at both the organ and cell levels, they remain resistant to arrhythmia. A marked increase in peak L-type calcium (I_CaL) current and slowed I_CaL decay elicited pronounced prolongation of initial repolarization, but faster late repolarization normalizes overall AP duration. Early afterdepolarizations were seldom observed in KO animals, and those that were observed exhibited a mechanism intermediate between murine and large mammal dynamical properties. As expected, spontaneous SR Ca²⁺ sparks and waves were virtually absent. Together these findings suggest that intact SR Ca²⁺ handling is an absolute requirement for triggered arrhythmia in the mouse, and that in its absence, dramatic changes to the major inward currents can be resisted by the substantial K⁺ current reserve, even at end-stage disease.

From Millimeters to Micrometers; Re-introducing Myocytes in Models of Cardiac Electrophysiology

Computational modeling has contributed significantly to present understanding of cardiac electrophysiology including cardiac conduction, excitation-contraction coupling, and the effects and side-effects of drugs. However, the accuracy of in silico analysis of electrochemical wave dynamics in cardiac tissue is limited by the homogenization procedure (spatial averaging) intrinsic to standard continuum models of conduction. Averaged models cannot resolve the intricate dynamics in the vicinity of individual cardiomyocytes simply because the myocytes are not present in these models. Here we demonstrate how recently developed mathematical models based on representing every myocyte can significantly increase the accuracy, and thus the utility of modeling electrophysiological function and dysfunction in collections of coupled cardiomyocytes. The present gold standard of numerical simulation for cardiac electrophysiology is based on the bidomain model. In the bidomain model, the extracellular (E) space, the cell membrane (M) and the intracellular (I) space are all assumed to be present everywhere in the tissue. Consequently, it is impossible to study biophysical processes taking place close to individual myocytes. The bidomain model represents the tissue by averaging over several hundred myocytes and this inherently limits the accuracy of the model. In our alternative approach both E, M, and I are represented in the model which is therefore referred to as the EMI model. The EMI model approach allows for detailed analysis of the biophysical processes going on in functionally important spaces very close to individual myocytes, although at the cost of significantly increased CPU-requirements.

A computational method for identifying an optimal combination of existing drugs to repair the action potentials of SQT1 ventricular myocytes

Mutations are known to cause perturbations in essential functional features of integral membrane proteins, including ion channels. Even restricted or point mutations can result in substantially changed properties of ion currents. The additive effect of these alterations for a specific ion channel can result in significantly changed properties of the action potential (AP). Both AP shortening and AP prolongation can result from known mutations, and the consequences can be life-threatening. Here, we present a computational method for identifying new drugs utilizing combinations of existing drugs. Based on the knowledge of theoretical effects of existing drugs on individual ion currents, our aim is to compute optimal combinations that can ‘repair’ the mutant AP waveforms so that the baseline AP-properties are restored. More specifically, we compute optimal, combined, drug concentrations such that the waveforms of the transmembrane potential and the cytosolic calcium concentration of the mutant cardiomyocytes (CMs) becomes as similar as possible to their wild type counterparts after the drug has been applied. In order to demonstrate the utility of this method, we address the question of computing an optimal drug for the short QT syndrome type 1 (SQT1). For the SQT1 mutation N588K, there are available data sets that describe the effect of various drugs on the mutated K⁺ channel. These published findings are the basis for our computational analysis which can identify optimal compounds in the sense that the AP of the mutant CMs resembles essential biomarkers of the wild type CMs. Using recently developed insights regarding electrophysiological properties among myocytes from different species, we compute optimal drug combinations for hiPSC-CMs, rabbit ventricular CMs and adult human ventricular CMs with the SQT1 mutation. Since the ‘composition’ of ion channels that form the AP is different for the three types of myocytes under consideration, so is the composition of the optimal drug.

Poly-pharmacology (using multiple drugs to treat disease) has been proposed for improving cardiac anti-arrhythmic therapy for at least two decades. However, the specific arrhythmia contexts in which polytherapy is likely to be both safe and effective have remained elusive. Type 1 short QT syndrome (SQT1) is a rare form of cardiac arrhythmia that results from mutations to the human Ether-á-go-go Related Gene (hERG) potassium channel. Functionally, these mutations are remarkably consistent in that they permit the channel to open earlier during each heart beat. While hundreds of compounds are known to inhibit hERG channels, the specific effect of SQT1 mutations that allows for early channel opening also limits the ability of most of those compounds to correct SQT1 dysfunction. Here, we have applied a suite of ventricular cardiomyocyte computational models to ask whether polytherapy may offer a more effective therapeutic strategy in SQT1, and if so, what the likely characteristics of that strategy are. Our analyses suggest that simultaneous induction of late sodium current and partial hERG blockade offers a promising strategy. While no activators of late sodium current have been clinically approved, several experimental compounds are available and may provide a basis for interrogating this strategy. The method presented here can be used to compute optimal drug combinations provided that the effect of each drug on every relevant ion channel is known.

Heart Muscle Microphysiological System for Cardiac Liability Prediction of Repurposed COVID-19 Therapeutics

Despite global efforts, it took 7 months between the proclamation of global SARS-CoV-2 pandemic and the first FDA-approved treatment for COVID-19. During this timeframe, clinicians focused their efforts on repurposing drugs, such as hydroxychloroquine (HCQ) or azithromycin (AZM) to treat hospitalized COVID-19 patients. While clinical trials are time-consuming, the exponential increase in hospitalizations compelled the FDA to grant an emergency use authorization for HCQ and AZM as treatment for COVID-19, although there was limited evidence of their combined efficacy and safety. The authorization was revoked 4 months later, giving rise to controversial political and scientific debates illustrating important challenges such as premature authorization of potentially ineffective or unsafe therapeutics, while diverting resources from screening of effective drugs. Here we report on a preclinical drug screening platform, a cardiac microphysiological system (MPS), to rapidly identify clinically relevant cardiac liabilities associated with HCQ and AZM. The cardiac MPS is a microfabricated fluidic system in which cardiomyocytes derived from human induced pluripotent stem cells self-arrange into a uniaxially beating tissue. The drug response was measured using outputs that correlate with clinical measurements such as action potential duration (proxy for clinical QT interval) and drug-biomarker pairing. The cardiac MPS predicted clinical arrhythmias associated with QT prolongation and rhythm instabilities in tissues treated with HCQ. We found no change in QT interval upon acute exposure to AZM, while still observing a significant increase in arrhythmic events. These results suggest that this MPS can not only predict arrhythmias, but it can also identify arrhythmias even when QT prolongation is absent. When exposed to HCQ and AZM polytherapy, this MPS faithfully reflected clinical findings, in that the combination of drugs synergistically increased QT interval when compared to single drug exposure, while not worsening the overall frequency of arrhythmic events. The high content cardiac MPS can rapidly evaluate the cardiac safety of potential therapeutics, ultimately accelerating patients' access to safe and effective treatments.

Sex-Specific Classification of Drug-Induced Torsade de Pointes Susceptibility Using Cardiac Simulations and Machine Learning

Torsade de Pointes (TdP), a rare but lethal ventricular arrhythmia, is a toxic side effect of many drugs. To assess TdP risk, safety regulatory guidelines require quantification of hERG channel block in vitro and QT interval prolongation in vivo for all new therapeutic compounds. Unfortunately, these have proven to be poor predictors of torsadogenic risk, and are likely to have prevented safe compounds from reaching clinical phases. Although this has stimulated numerous efforts to define new paradigms for cardiac safety, none of the recently developed strategies accounts for patient conditions. In particular, despite being a well-established independent risk factor for TdP, female sex is vastly under-represented in both basic research and clinical studies, and thus current TdP metrics are likely biased toward the male sex. Here, we apply statistical learning to synthetic data, generated by simulating drug effects on cardiac myocyte models capturing male and female electrophysiology, to develop new sex-specific classification frameworks for TdP risk. We show that (i) TdP classifiers require different features in females vs. males; (ii) male-based classifiers perform more poorly when applied to female data; and (iii) female-based classifier performance is largely unaffected by acute effects of hormones (i.e., during various phases of the menstrual cycle). Notably, when predicting TdP risk of intermediate drugs on female simulated data, male-biased predictive models consistently underestimate TdP risk in women. Therefore, we conclude that pipelines for preclinical cardiotoxicity risk assessment should consider sex as a key variable to avoid potentially life-threatening consequences for the female population.

In vitro safety "clinical trial" of the cardiac liability of drug polytherapy

Only a handful of US Food and Drug Administration (FDA) Emergency Use Authorizations exist for drug and biologic therapeutics that treat severe acute respiratory syndrome-coronavirus 2 (SARS-CoV-2) infection. Potential therapeutics include repurposed drugs, some with cardiac liabilities. We report on a chronic preclinical drug screening platform, a cardiac microphysiological system (MPS), to assess cardiotoxicity associated with repurposed hydroxychloroquine (HCQ) and azithromycin (AZM) polytherapy in a mock phase I safety clinical trial. The MPS contained human heart muscle derived from induced pluripotent stem cells. The effect of drug response was measured using outputs that correlate with clinical measurements, such as QT interval (action potential duration) and drug-biomarker pairing. Chronic exposure (10 days) of heart muscle to HCQ alone elicited early afterdepolarizations and increased QT interval past 5 days. AZM alone elicited an increase in QT interval from day 7 onward, and arrhythmias were observed at days 8 and 10. Monotherapy results mimicked clinical trial outcomes. Upon chronic exposure to HCQ and AZM polytherapy, we observed an increase in QT interval on days 4-8. Interestingly, a decrease in arrhythmias and instabilities was observed in polytherapy relative to monotherapy, in concordance with published clinical trials. Biomarkers, most of them measurable in patients' serum, were identified for negative effects of monotherapy or polytherapy on tissue contractile function, morphology, and antioxidant protection. The cardiac MPS correctly predicted clinical arrhythmias associated with QT prolongation and rhythm instabilities. This high content system can help clinicians design their trials, rapidly project cardiac outcomes, and define new monitoring biomarkers to accelerate access of patients to safe coronavirus disease 2019 (COVID-19) therapeutics.

In Vitro Safety "Clinical Trial" of the Cardiac Liability of Hydroxychloroquine and Azithromycin as COVID19 Polytherapy

Despite global efforts, there are no effective FDA-approved medicines for the treatment of SARS-CoV-2 infection. Potential therapeutics focus on repurposed drugs, some with cardiac liabilities. Here we report on a preclinical drug screening platform, a cardiac microphysiological system (MPS), to assess cardiotoxicity associated with hydroxychloroquine (HCQ) and azithromycin (AZM) polytherapy in a mock clinical trial. The MPS contained human heart muscle derived from patient-specific induced pluripotent stem cells. The effect of drug response was measured using outputs that correlate with clinical measurements such as QT interval (action potential duration) and drug-biomarker pairing. Chronic exposure to HCQ alone elicited early afterdepolarizations (EADs) and increased QT interval from day 6 onwards. AZM alone elicited an increase in QT interval from day 7 onwards and arrhythmias were observed at days 8 and 10. Monotherapy results closely mimicked clinical trial outcomes. Upon chronic exposure to HCQ and AZM polytherapy, we observed an increase in QT interval on days 4-8.. Interestingly, a decrease in arrhythmias and instabilities was observed in polytherapy relative to monotherapy, in concordance with published clinical trials. Furthermore, biomarkers, most of them measurable in patients' serum, were identified for negative effects of single drug or polytherapy on tissue contractile function, morphology, and antioxidant protection. The cardiac MPS can predict clinical arrhythmias associated with QT prolongation and rhythm instabilities. This high content system can help clinicians design their trials, rapidly project cardiac outcomes, and define new monitoring biomarkers to accelerate access of patients to safe COVID-19 therapeutics.

Populations of in silico myocytes and tissues reveal synergy of multiatrial-predominant K+ -current block in atrial fibrillation

BACKGROUND AND PURPOSE

Pharmacotherapy of atrial fibrillation (AF), the most common cardiac arrhythmia, remains unsatisfactory due to low efficacy and safety concerns. New therapeutic strategies target atrial-predominant ion-channels and involve multichannel block (poly)therapy. As AF is characterized by rapid and irregular atrial activations, compounds displaying potent antiarrhythmic effects at fast and minimal effects at slow rates are desirable. We present a novel systems pharmacology framework to quantitatively evaluate synergistic anti-AF effects of combined block of multiple atrial-predominant K⁺ currents (ultra-rapid delayed rectifier K⁺ current, I_Kur , small conductance Ca²⁺ -activated K⁺ current, I_KCa , K_2P 3.1 2-pore-domain K⁺ current, I_K2P ) in AF.

EXPERIMENTAL APPROACH

We constructed experimentally calibrated populations of virtual atrial myocyte models in normal sinus rhythm and AF-remodelled conditions using two distinct, well-established atrial models. Sensitivity analyses on our atrial populations was used to investigate the rate dependence of action potential duration (APD) changes due to blocking I_Kur , I_K2P or I_KCa and interactions caused by blocking of these currents in modulating APD. Block was simulated in both single myocytes and one-dimensional tissue strands to confirm insights from the sensitivity analyses and examine anti-arrhythmic effects of multi-atrial-predominant K⁺ current block in single cells and coupled tissue.

KEY RESULTS

In both virtual atrial myocytes and tissues, multiple atrial-predominant K⁺ -current block promoted favourable positive rate-dependent APD prolongation and displayed positive rate-dependent synergy, that is, increasing synergistic antiarrhythmic effects at fast pacing versus slow rates.

CONCLUSION AND IMPLICATIONS

Simultaneous block of multiple atrial-predominant K⁺ currents may be a valuable antiarrhythmic pharmacotherapeutic strategy for AF.

A unified computational model for cortical post-synaptic plasticity

Signalling pathways leading to post-synaptic plasticity have been examined in many types of experimental studies, but a unified picture on how multiple biochemical pathways collectively shape neocortical plasticity is missing. We built a biochemically detailed model of post-synaptic plasticity describing CaMKII, PKA, and PKC pathways and their contribution to synaptic potentiation or depression. We developed a statistical AMPA-receptor-tetramer model, which permits the estimation of the AMPA-receptor-mediated maximal synaptic conductance based on numbers of GluR1s and GluR2s predicted by the biochemical signalling model. We show that our model reproduces neuromodulator-gated spike-timing-dependent plasticity as observed in the visual cortex and can be fit to data from many cortical areas, uncovering the biochemical contributions of the pathways pinpointed by the underlying experimental studies. Our model explains the dependence of different forms of plasticity on the availability of different proteins and can be used for the study of mental disorder-associated impairments of cortical plasticity.

Hypokalemia Promotes Arrhythmia by Distinct Mechanisms in Atrial and Ventricular Myocytes

RATIONALE

Hypokalemia occurs in up to 20% of hospitalized patients and is associated with increased incidence of ventricular and atrial fibrillation. It is unclear whether these differing types of arrhythmia result from direct and perhaps distinct effects of hypokalemia on cardiomyocytes.

OBJECTIVE

To investigate proarrhythmic mechanisms of hypokalemia in ventricular and atrial myocytes.

METHODS AND RESULTS

Experiments were performed in isolated rat myocytes exposed to simulated hypokalemia conditions (reduction of extracellular [K⁺] from 5.0 to 2.7 mmol/L) and supported by mathematical modeling studies. Ventricular cells subjected to hypokalemia exhibited Ca²⁺ overload and increased generation of both spontaneous Ca²⁺ waves and delayed afterdepolarizations. However, similar Ca²⁺-dependent spontaneous activity during hypokalemia was only observed in a minority of atrial cells that were observed to contain t-tubules. This effect was attributed to close functional pairing of the Na⁺-K⁺ ATPase and Na⁺-Ca²⁺ exchanger proteins within these structures, as reduction in Na⁺ pump activity locally inhibited Ca²⁺ extrusion. Ventricular myocytes and tubulated atrial myocytes additionally exhibited early afterdepolarizations during hypokalemia, associated with Ca²⁺ overload. However, early afterdepolarizations also occurred in untubulated atrial cells, despite Ca²⁺ quiescence. These phase-3 early afterdepolarizations were rather linked to reactivation of nonequilibrium Na⁺ current, as they were rapidly blocked by tetrodotoxin. Na⁺ current-driven early afterdepolarizations in untubulated atrial cells were enabled by membrane hyperpolarization during hypokalemia and short action potential configurations. Brief action potentials were in turn maintained by ultra-rapid K⁺ current (I_Kur); a current which was found to be absent in tubulated atrial myocytes and ventricular myocytes.

CONCLUSIONS

Distinct mechanisms underlie hypokalemia-induced arrhythmia in the ventricle and atrium but also vary between atrial myocytes depending on subcellular structure and electrophysiology.

Secretoneurin Is an Endogenous Calcium/Calmodulin-Dependent Protein Kinase II Inhibitor That Attenuates Ca2+-Dependent Arrhythmia

BACKGROUND

Circulating SN (secretoneurin) concentrations are increased in patients with myocardial dysfunction and predict poor outcome. Because SN inhibits CaMKIIδ (Ca²⁺/calmodulin-dependent protein kinase IIδ) activity, we hypothesized that upregulation of SN in patients protects against cardiomyocyte mechanisms of arrhythmia.

METHODS

Circulating levels of SN and other biomarkers were assessed in patients with catecholaminergic polymorphic ventricular tachycardia (CPVT; n=8) and in resuscitated patients after ventricular arrhythmia-induced cardiac arrest (n=155). In vivo effects of SN were investigated in CPVT mice (RyR2 [ryanodine receptor 2]-R2474S) using adeno-associated virus-9-induced overexpression. Interactions between SN and CaMKIIδ were mapped using pull-down experiments, mutagenesis, ELISA, and structural homology modeling. Ex vivo actions were tested in Langendorff hearts and effects on Ca²⁺ homeostasis examined by fluorescence (fluo-4) and patch-clamp recordings in isolated cardiomyocytes.

RESULTS

SN levels were elevated in patients with CPVT and following ventricular arrhythmia-induced cardiac arrest. In contrast to NT-proBNP (N-terminal pro-B-type natriuretic peptide) and hs-TnT (high-sensitivity troponin T), circulating SN levels declined after resuscitation, as the risk of a new arrhythmia waned. Myocardial pro-SN expression was also increased in CPVT mice, and further adeno-associated virus-9-induced overexpression of SN attenuated arrhythmic induction during stress testing with isoproterenol. Mechanistic studies mapped SN binding to the substrate binding site in the catalytic region of CaMKIIδ. Accordingly, SN attenuated isoproterenol induced autophosphorylation of Thr287-CaMKIIδ in Langendorff hearts and inhibited CaMKIIδ-dependent RyR phosphorylation. In line with CaMKIIδ and RyR inhibition, SN treatment decreased Ca²⁺ spark frequency and dimensions in cardiomyocytes during isoproterenol challenge, and reduced the incidence of Ca²⁺ waves, delayed afterdepolarizations, and spontaneous action potentials. SN treatment also lowered the incidence of early afterdepolarizations during isoproterenol; an effect paralleled by reduced magnitude of L-type Ca²⁺ current.

CONCLUSIONS

SN production is upregulated in conditions with cardiomyocyte Ca²⁺ dysregulation and offers compensatory protection against cardiomyocyte mechanisms of arrhythmia, which may underlie its putative use as a biomarker in at-risk patients.

General Principles for the Validation of Proarrhythmia Risk Prediction Models: An Extension of the CiPA In Silico Strategy

This white paper presents principles for validating proarrhythmia risk prediction models for regulatory use as discussed at the In Silico Breakout Session of a Cardiac Safety Research Consortium/Health and Environmental Sciences Institute/US Food and Drug Administration-sponsored Think Tank Meeting on May 22, 2018. The meeting was convened to evaluate the progress in the development of a new cardiac safety paradigm, the Comprehensive in Vitro Proarrhythmia Assay (CiPA). The opinions regarding these principles reflect the collective views of those who participated in the discussion of this topic both at and after the breakout session. Although primarily discussed in the context of in silico models, these principles describe the interface between experimental input and model-based interpretation and are intended to be general enough to be applied to other types of nonclinical models for proarrhythmia assessment. This document was developed with the intention of providing a foundation for more consistency and harmonization in developing and validating different models for proarrhythmia risk prediction using the example of the CiPA paradigm.

Arrhythmogenic Current Generation by Myofilament-Triggered Ca2+ Release and Sarcomere Heterogeneity

Heterogeneous mechanical dyskinesis has been implicated in many arrhythmogenic phenotypes. Strain-dependent perturbations to cardiomyocyte electrophysiology may contribute to this arrhythmogenesis through processes referred to as mechanoelectric feedback. Although the role of stretch-activated ion currents has been investigated using computational models, experimental studies suggest that mechanical strain may also promote arrhythmia by facilitating calcium wave propagation. To investigate whether strain-dependent changes in calcium affinity to the myofilament may promote arrhythmogenic intracellular calcium waves, we modified a mathematical model of rabbit excitation-contraction coupling coupled to a model of myofilament activation and force development. In a one-dimensional compartmental analysis, we bidirectionally coupled 50 sarcomere models in series to model calcium diffusion and stress transfer between adjacent sarcomeres. These considerations enabled the model to capture 1) the effects of mechanical feedback on calcium homeostasis at the sarcomeric level and 2) the combined effects of mechanical and calcium heterogeneities at the cellular level. The results suggest that in conditions of calcium overload, the vulnerable window of stretch-release to trigger suprathreshold delayed afterdepolarizations can be affected by heterogeneity in sarcomere length. Furthermore, stretch and sarcomere heterogeneity may modulate the susceptibility threshold for delayed afterdepolarizations and the aftercontraction wave propagation velocity.

Properties of cardiac conduction in a cell-based computational model

The conduction of electrical signals through cardiac tissue is essential for maintaining the function of the heart, and conduction abnormalities are known to potentially lead to life-threatening arrhythmias. The properties of cardiac conduction have therefore been the topic of intense study for decades, but a number of questions related to the mechanisms of conduction still remain unresolved. In this paper, we demonstrate how the so-called EMI model may be used to study some of these open questions. In the EMI model, the extracellular space, the cell membrane, the intracellular space and the cell connections are all represented as separate parts of the computational domain, and the model therefore allows for study of local properties that are hard to represent in the classical homogenized bidomain or monodomain models commonly used to study cardiac conduction. We conclude that a non-uniform sodium channel distribution increases the conduction velocity and decreases the time delays over gap junctions of reduced coupling in the EMI model simulations. We also present a theoretical optimal cell length with respect to conduction velocity and consider the possibility of ephaptic coupling (i.e. cell-to-cell coupling through the extracellular potential) acting as an alternative or supporting mechanism to gap junction coupling. We conclude that for a non-uniform distribution of sodium channels and a sufficiently small intercellular distance, ephaptic coupling can influence the dynamics of the sodium channels and potentially provide cell-to-cell coupling when the gap junction connection is absent.

The electrochemical wave traversing the heart during every beat is essential for cardiac pumping function and supply of blood to the body. Understanding the stability of this wave is crucial to understanding how lethal arrhythmias are generated. Despite this importance, our knowledge of the physical determinants of wave propagation are still evolving. One particular challenge has been the lack of accurate mathematical models of conduction at the cellular level. Because cardiac muscle is an electrical syncytium, in which direct charge transfer between cells drives wave propagation, classical bidomain and monodomain tissue models employ a homogenized approximation of this process. This approximation is not valid at the length scale of single cells, and prevents any analysis of how cellular structures impact cardiac conduction. Instead, so-called microdomain models must be used for these questions. Here we utilize a recently developed modelling framework that is well suited to represent small collections of cells. By applying this framework, we show that concentration of sodium channels at the longitudinal borders of myocytes accelerates cardiac conduction. We also demonstrate that when juxtaposed cells are sufficiently close, this non-uniform distribution induces large ephaptic currents, which contribute to intercellular coupling.

Alterations in Schizophrenia-Associated Genes Can Lead to Increased Power in Delta Oscillations

Genome-wide association studies have implicated many ion channels in schizophrenia pathophysiology. Although the functions of these channels are relatively well characterized by single-cell studies, the contributions of common variation in these channels to neurophysiological biomarkers and symptoms of schizophrenia remain elusive. Here, using computational modeling, we show that a common biomarker of schizophrenia, namely, an increase in delta-oscillation power, may be a direct consequence of altered expression or kinetics of voltage-gated ion channels or calcium transporters. Our model of a circuit of layer V pyramidal cells highlights multiple types of schizophrenia-related variants that contribute to altered dynamics in the delta-frequency band. Moreover, our model predicts that the same membrane mechanisms that increase the layer V pyramidal cell network gain and response to delta-frequency oscillations may also cause a deficit in a single-cell correlate of the prepulse inhibition, which is a behavioral biomarker highly associated with schizophrenia.

Biophysical Psychiatry-How Computational Neuroscience Can Help to Understand the Complex Mechanisms of Mental Disorders

The brain is the most complex of human organs, and the pathophysiology underlying abnormal brain function in psychiatric disorders is largely unknown. Despite the rapid development of diagnostic tools and treatments in most areas of medicine, our understanding of mental disorders and their treatment has made limited progress during the last decades. While recent advances in genetics and neuroscience have a large potential, the complexity and multidimensionality of the brain processes hinder the discovery of disease mechanisms that would link genetic findings to clinical symptoms and behavior. This applies also to schizophrenia, for which genome-wide association studies have identified a large number of genetic risk loci, spanning hundreds of genes with diverse functionalities. Importantly, the multitude of the associated variants and their prevalence in the healthy population limit the potential of a reductionist functional genetics approach as a stand-alone solution to discover the disease pathology. In this review, we outline the key concepts of a "biophysical psychiatry," an approach that employs large-scale mechanistic, biophysics-founded computational modelling to increase transdisciplinary understanding of the pathophysiology and strive toward robust predictions. We discuss recent scientific advances that allow a synthesis of previously disparate fields of psychiatry, neurophysiology, functional genomics, and computational modelling to tackle open questions regarding the pathophysiology of heritable mental disorders. We argue that the complexity of the increasing amount of genetic data exceeds the capabilities of classical experimental assays and requires computational approaches. Biophysical psychiatry, based on modelling diseased brain networks using existing and future knowledge of basic genetic, biochemical, and functional properties on a single neuron to a microcircuit level, may allow a leap forward in deriving interpretable biomarkers and move the field toward novel treatment options.

3D dSTORM imaging reveals novel detail of ryanodine receptor localization in rat cardiac myocytes

KEY POINTS

Using 3D direct stochastic optical reconstruction microscopy (dSTORM), we developed novel approaches to quantitatively describe the nanoscale, 3D organization of ryanodine receptors (RyRs) in cardiomyocytes. Complex arrangements of RyR clusters were observed in 3D space, both at the cell surface and within the cell interior, with allocation to dyadic and non-dyadic pools. 3D imaging importantly allowed discernment of clusters overlapping in the z-axis, for which detection was obscured by conventional 2D imaging techniques. Thus, RyR clusters were found to be significantly smaller than previous 2D estimates. Ca2+ release units (CRUs), i.e. functional groupings of neighbouring RyR clusters, were similarly observed to be smaller than earlier reports. Internal CRUs contained more RyRs in more clusters than CRUs on the cell surface, and yielded longer duration Ca2+ sparks.

ABSTRACT

Cardiomyocyte contraction is dependent on Ca2+ release from ryanodine receptors (RyRs). However, the precise localization of RyRs remains unknown, due to shortcomings of imaging techniques which are diffraction limited or restricted to 2D. We aimed to determine the 3D nanoscale organization of RyRs in rat cardiomyocytes by employing direct stochastic optical reconstruction microscopy (dSTORM) with phase ramp technology. Initial observations at the cell surface showed an undulating organization of RyR clusters, resulting in their frequent overlap in the z-axis and obscured detection by 2D techniques. Non-overlapping clusters were imaged to create a calibration curve for estimating RyR number based on recorded fluorescence blinks. Employing this method at the cell surface and interior revealed smaller RyR clusters than 2D estimates, as erroneous merging of axially aligned RyRs was circumvented. Functional groupings of RyR clusters (Ca2+ release units, CRUs), contained an average of 18 and 23 RyRs at the surface and interior, respectively, although half of all CRUs contained only a single 'rogue' RyR. Internal CRUs were more tightly packed along z-lines than surface CRUs, contained larger and more numerous RyR clusters, and constituted ∼75% of the roughly 1 million RyRs present in an average cardiomyocyte. This complex internal 3D geometry was underscored by correlative imaging of RyRs and t-tubules, which enabled quantification of dyadic and non-dyadic RyR populations. Mirroring differences in CRU size and complexity, Ca2+ sparks originating from internal CRUs were of longer duration than those at the surface. These data provide novel, nanoscale insight into RyR organization and function across cardiomyocytes.

Arrhythmia initiation in catecholaminergic polymorphic ventricular tachycardia type 1 depends on both heart rate and sympathetic stimulation

Aims

Catecholaminergic polymorphic ventricular tachycardia type 1 (CPVT1) predisposes to ventricular tachyarrhythmias (VTs) during high heart rates due to physical or psychological stress. The essential role of catecholaminergic effects on ventricular cardiomyocytes in this situation is well documented, but the importance of heart rate per se for arrhythmia initiation in CPVT1 is largely unexplored.

Methods and results

Sixteen CPVT1 patients performed a bicycle stress-test. Occurrence of VT triggers, i.e. premature ventricular complexes (PVC), depended on high heart rate, with individual thresholds. Atrial pacing above the individual PVC threshold in three patients did not induce PVCs.

The underlying mechanism for the clinical observation was explored using cardiomyocytes from mice with the RyR2-R2474S (RyR2-RS) mutation, which exhibit exercise-induced VTs. While rapid pacing increased the number of Ca²⁺ waves in both RyR2-RS and wild-type (p<0.05), β-adrenoceptor (βAR) stimulation induced more Ca²⁺ waves in RyR2-RS (p<0.05). Notably, Ca²⁺ waves occurred despite decreased sarcoplasmic reticulum (SR) Ca²⁺ content in RyR2-RS (p<0.05), suggesting increased cytosolic RyR2 Ca²⁺ sensitivity.

A computational model of mouse ventricular cardiomyocyte electrophysiology reproduced the cellular CPVT1 phenotype when RyR2 Ca²⁺ sensitivity was increased. Importantly, diastolic fluctuations in phosphorylation of RyR2 and SR Ca²⁺ content determined Ca²⁺ wave initiation. These factors were modulated towards increased propensity for arrhythmia initiation by increased pacing rates, but even more by βAR stimulation.

Conclusion

In CPVT1, VT propensity depends on individual heart rate thresholds for PVCs. Through converging data from clinical exercise stress-testing, cellular studies and computational modelling, we confirm the heart rate-independent pro-arrhythmic effects of βAR stimulation in CPVT1, but also identify an independent and synergistic contribution from effects of high heart rate.

Ryanodine receptor dispersion disrupts Ca2+ release in failing cardiac myocytes

Reduced cardiac contractility during heart failure (HF) is linked to impaired Ca²⁺ release from Ryanodine Receptors (RyRs). We investigated whether this deficit can be traced to nanoscale RyR reorganization. Using super-resolution imaging, we observed dispersion of RyR clusters in cardiomyocytes from post-infarction HF rats, resulting in more numerous, smaller clusters. Functional groupings of RyR clusters which produce Ca²⁺ sparks (Ca²⁺ release units, CRUs) also became less solid. An increased fraction of small CRUs in HF was linked to augmented ‘silent’ Ca²⁺ leak, not visible as sparks. Larger multi-cluster CRUs common in HF also exhibited low fidelity spark generation. When successfully triggered, sparks in failing cells displayed slow kinetics as Ca²⁺ spread across dispersed CRUs. During the action potential, these slow sparks protracted and desynchronized the overall Ca²⁺ transient. Thus, nanoscale RyR reorganization during HF augments Ca²⁺ leak and slows Ca²⁺ release kinetics, leading to weakened contraction in this disease.

The muscle cells of the heart coordinate how they contract and relax in order to produce the heartbeat. During heart failure, these cells become less able to contract. As a result the heart becomes inefficient, pumping less blood around the body.

For the cardiac muscle cells to contract, the levels of calcium ions in the cells needs to rapidly increase. In failing hearts, these increases in calcium ion levels are smaller, slower and less well coordinated. It was not known what causes these changes, making it difficult to treat heart failure.

Calcium ions are released in cardiac muscle cells through protein channels called ryanodine receptors. These receptors form clusters that allow them to synchronize when they open and close. Could the reorganization of ryanodine receptors account for the problems seen in failing hearts? To investigate, Kolstad et al. examined rat hearts using a technique called super-resolution microscopy. This showed that the clusters of ryanodine receptors break apart during heart failure to form smaller clusters. Further experiments showed that calcium ions ‘leak’ from these smaller clusters, reducing the amount of calcium that can be released into cardiac muscle cells during each heartbeat. Released calcium also spreads between the dispersed clusters, resulting in a slower rise of the calcium levels in the cells. Both changes contribute to weakened contractions of cells in failing hearts.

Therefore, heart failure can be traced back to very small rearrangements of the ryanodine receptors. This understanding will help researchers as they investigate new ways to treat heart failure.

Computational Modeling of Electrophysiology and Pharmacotherapy of Atrial Fibrillation: Recent Advances and Future Challenges

The pathophysiology of atrial fibrillation (AF) is broad, with components related to the unique and diverse cellular electrophysiology of atrial myocytes, structural complexity, and heterogeneity of atrial tissue, and pronounced disease-associated remodeling of both cells and tissue. A major challenge for rational design of AF therapy, particularly pharmacotherapy, is integrating these multiscale characteristics to identify approaches that are both efficacious and independent of ventricular contraindications. Computational modeling has long been touted as a basis for achieving such integration in a rapid, economical, and scalable manner. However, computational pipelines for AF-specific drug screening are in their infancy, and while the field is progressing quite rapidly, major challenges remain before computational approaches can fill the role of workhorse in rational design of AF pharmacotherapies. In this review, we briefly detail the unique aspects of AF pathophysiology that determine requirements for compounds targeting AF rhythm control, with emphasis on delimiting mechanisms that promote AF triggers from those providing substrate or supporting reentry. We then describe modeling approaches that have been used to assess the outcomes of drugs acting on established AF targets, as well as on novel promising targets including the ultra-rapidly activating delayed rectifier potassium current, the acetylcholine-activated potassium current and the small conductance calcium-activated potassium channel. Finally, we describe how heterogeneity and variability are being incorporated into AF-specific models, and how these approaches are yielding novel insights into the basic physiology of disease, as well as aiding identification of the important molecular players in the complex AF etiology.

Pleiotropic effects of schizophrenia-associated genetic variants in neuron firing and cardiac pacemaking revealed by computational modeling

Schizophrenia patients have an increased risk of cardiac dysfunction. A possible factor underlying this comorbidity are the common variants in the large set of genes that have recently been discovered in genome-wide association studies (GWASs) as risk genes of schizophrenia. Many of these genes control the cell electrogenesis and calcium homeostasis. We applied biophysically detailed models of layer V pyramidal cells and sinoatrial node cells to study the contribution of schizophrenia-associated genes on cellular excitability. By including data from functional genomics literature to simulate the effects of common variants of these genes, we showed that variants of voltage-gated Na⁺ channel or hyperpolarization-activated cation channel-encoding genes cause qualitatively similar effects on layer V pyramidal cell and sinoatrial node cell excitability. By contrast, variants of Ca²⁺ channel or transporter-encoding genes mostly have opposite effects on cellular excitability in the two cell types. We also show that the variants may crucially affect the propagation of the cardiac action potential in the sinus node. These results may help explain some of the cardiac comorbidity in schizophrenia, and may facilitate generation of effective antipsychotic medications without cardiac side-effects such as arrhythmia.

An integrative appraisal of mechano-electric feedback mechanisms in the heart

Mechanically-induced alterations in cardiac electrophysiology are referred to as mechano-electric feedback (MEF), and play an important role in electrical regulation of cardiac performance. The influence of mechanical stress and strain on electrophysiology has been investigated at all levels, however the role of MEF in arrhythmia remains poorly understood. During the normal contraction of the heart, mechano-sensitive processes are an implicit component of cardiac activity. Under abnormal mechanical events, stretch-activated mechanisms may contribute to local or global changes in electrophysiology (EP). While such mechanisms have been hypothesised to be involved in mechanically-initiated arrhythmias, the details of these mechanisms and their importance remain elusive. We assess the theoretical role of stretch mechanisms using coupled models of cellular electrophysiology and sarcomere contraction dynamics. Using models of single ventricular myocytes, we first investigated the potential MEF contributions of stretch-activated currents (SAC), and stretch-induced myofilament calcium release, to test how strain and fibrosis may alter cellular electrophysiology. For all models investigated, SACs were alone not sufficient to create a pro-arrhythmic perturbation of the action potential with stretch. However, when combined with stretch-induced myofilament calcium release, the action potential could be shortened depending on the timing of the strain. This effect was highly model dependent, with a canine epicardial EP model being the most sensitive. These model results suggest that known mechanisms of mechano-electric coupling in cardiac myocyte may be sufficient to be pro-arrhythmic, but only in combination and under specific strain patterns.

An Evaluation of the Accuracy of Classical Models for Computing the Membrane Potential and Extracellular Potential for Neurons

Two mathematical models are part of the foundation of Computational neurophysiology; (a) the Cable equation is used to compute the membrane potential of neurons, and, (b) volume-conductor theory describes the extracellular potential around neurons. In the standard procedure for computing extracellular potentials, the transmembrane currents are computed by means of (a) and the extracellular potentials are computed using an explicit sum over analytical point-current source solutions as prescribed by volume conductor theory. Both models are extremely useful as they allow huge simplifications of the computational efforts involved in computing extracellular potentials. However, there are more accurate, though computationally very expensive, models available where the potentials inside and outside the neurons are computed simultaneously in a self-consistent scheme. In the present work we explore the accuracy of the classical models (a) and (b) by comparing them to these more accurate schemes. The main assumption of (a) is that the ephaptic current can be ignored in the derivation of the Cable equation. We find, however, for our examples with stylized neurons, that the ephaptic current is comparable in magnitude to other currents involved in the computations, suggesting that it may be significant-at least in parts of the simulation. The magnitude of the error introduced in the membrane potential is several millivolts, and this error also translates into errors in the predicted extracellular potentials. While the error becomes negligible if we assume the extracellular conductivity to be very large, this assumption is, unfortunately, not easy to justify a priori for all situations of interest.

Studying dyadic structure-function relationships: a review of current modeling approaches and new insights into Ca2+ (mis)handling

Excitation–contraction coupling in cardiac myocytes requires calcium influx through L-type calcium channels in the sarcolemma, which gates calcium release through sarcoplasmic reticulum ryanodine receptors in a process known as calcium-induced calcium release, producing a myoplasmic calcium transient and enabling cardiomyocyte contraction. The spatio-temporal dynamics of calcium release, buffering, and reuptake into the sarcoplasmic reticulum play a central role in excitation–contraction coupling in both normal and diseased cardiac myocytes. However, further quantitative understanding of these cells’ calcium machinery and the study of mechanisms that underlie both normal cardiac function and calcium-dependent etiologies in heart disease requires accurate knowledge of cardiac ultrastructure, protein distribution and subcellular function. As current imaging techniques are limited in spatial resolution, limiting insight into changes in calcium handling, computational models of excitation–contraction coupling have been increasingly employed to probe these structure–function relationships. This review will focus on the development of structural models of cardiac calcium dynamics at the subcellular level, orienting the reader broadly towards the development of models of subcellular calcium handling in cardiomyocytes. Specific focus will be given to progress in recent years in terms of multi-scale modeling employing resolved spatial models of subcellular calcium machinery. A review of the state-of-the-art will be followed by a review of emergent insights into calcium-dependent etiologies in heart disease and, finally, we will offer a perspective on future directions for related computational modeling and simulation efforts.

Species-Dependent Mechanisms of Cardiac Arrhythmia: A Cellular Focus

Although ventricular arrhythmia remains a leading cause of morbidity and mortality, available antiarrhythmic drugs have limited efficacy. Disappointing progress in the development of novel, clinically relevant antiarrhythmic agents may partly be attributed to discrepancies between humans and animal models used in preclinical testing. However, such differences are at present difficult to predict, requiring improved understanding of arrhythmia mechanisms across species. To this end, we presently review interspecies similarities and differences in fundamental cardiomyocyte electrophysiology and current understanding of the mechanisms underlying the generation of afterdepolarizations and reentry. We specifically highlight patent shortcomings in small rodents to reproduce cellular and tissue-level arrhythmia substrate believed to be critical in human ventricle. Despite greater ease of translation from larger animal models, discrepancies remain and interpretation can be complicated by incomplete knowledge of human ventricular physiology due to low availability of explanted tissue. We therefore point to the benefits of mathematical modeling as a translational bridge to understanding and treating human arrhythmia.

Computing rates of Markov models of voltage-gated ion channels by inverting partial differential equations governing the probability density functions of the conducting and non-conducting states

Markov models are ubiquitously used to represent the function of single ion channels. However, solving the inverse problem to construct a Markov model of single channel dynamics from bilayer or patch-clamp recordings remains challenging, particularly for channels involving complex gating processes. Methods for solving the inverse problem are generally based on data from voltage clamp measurements. Here, we describe an alternative approach to this problem based on measurements of voltage traces. The voltage traces define probability density functions of the functional states of an ion channel. These probability density functions can also be computed by solving a deterministic system of partial differential equations. The inversion is based on tuning the rates of the Markov models used in the deterministic system of partial differential equations such that the solution mimics the properties of the probability density function gathered from (pseudo) experimental data as well as possible. The optimization is done by defining a cost function to measure the difference between the deterministic solution and the solution based on experimental data. By evoking the properties of this function, it is possible to infer whether the rates of the Markov model are identifiable by our method. We present applications to Markov model well-known from the literature.

The outcome of recurrent varicose vein surgery: the patients' perspective

Objectives: This study was performed to assess the patient's perception of the outcome of recurrent varicose vein surgery.

Methods: Patients who underwent surgery for recurrent varicose veins between 1 January 1995 and 31 December 1998 were sent a questionnaire and Aberdeen Varicose Vein Severity Score (AVVSS) in September 2003.

Results: A total of 48 questionnaires were distributed and 40 responses (83%) were received. The mean length of follow-up was 81.78 months. Outcome was 'excellent' (no varicose veins) for four patients, 'satisfactory' for 22 (small, asymptomatic recurrence), 'unsatisfactory' for 10 (symptomatic recurrence) and 'worse' for four patients; 65% thought that the outcome was successful (excellent or satisfactory) and 35% a failure. The median AVVSS in the success group was 11.50 (IQR 7.94–19.69) compared with 29.63 (IQR 18.46–40.31) in the failure group ( P = 0.0009, Wilcoxon).

Conclusions: Two-thirds of patients undergoing surgery for recurrent varicose veins are satisfied with the outcome, nearly seven years after surgery. One-third will regard surgery as having failed. Patients should be advised of these potential outcomes prior to surgery.

Field-measured drag area is a key correlate of level cycling time trial performance

Drag area (Ad ) is a primary factor determining aerodynamic resistance during level cycling and is therefore a key determinant of level time trial performance. However, Ad has traditionally been difficult to measure. Our purpose was to determine the value of adding field-measured Ad as a correlate of level cycling time trial performance. In the field, 19 male cyclists performed a level (22.1 km) time trial. Separately, field-determined Ad and rolling resistance were calculated for subjects along with projected frontal area assessed directly (AP ) and indirectly (Est AP ). Also, a graded exercise test was performed to determine [Formula: see text] peak, lactate threshold (LT), and economy. [Formula: see text] peak ([Formula: see text]) and power at LT were significantly correlated to power measured during the time trial (r = 0.83 and 0.69, respectively) but were not significantly correlated to performance time (r = - 0.42 and -0.45). The correlation with performance time improved significantly (p < 0.05) when these variables were normalized to Ad . Of note, Ad alone was better correlated to performance time (r = 0.85, p < 0.001) than any combination of non-normalized physiological measure. The best correlate with performance time was field-measured power output during the time trial normalized to Ad (r = - 0.92). AP only accounted for 54% of the variability in Ad . Accordingly, the correlation to performance time was significantly lower using power normalized to AP (r = - 0.75) or Est AP (r = - 0.71). In conclusion, unless normalized to Ad , level time trial performance in the field was not highly correlated to common laboratory measures. Furthermore, our field-measured Ad is easy to determine and was the single best predictor of level time trial performance.

Atrial-selective targeting of arrhythmogenic phase-3 early afterdepolarizations in human myocytes

BACKGROUND

We have previously shown that non-equilibrium Na(+) current (INa) reactivation drives isoproterenol-induced phase-3 early afterdepolarizations (EADs) in mouse ventricular myocytes. In these cells, EAD initiation occurs secondary to potentiated sarcoplasmic reticulum Ca(2+) release and enhanced Na(+)/Ca(2+) exchange (NCX). This can be abolished by tetrodotoxin-blockade of INa, but not ranolazine, which selectively inhibits ventricular late INa.

AIM

Since repolarization of human atrial myocytes is similar to mouse ventricular myocytes in that it is relatively rapid and potently modulated by Ca(2+), we investigated whether similar mechanisms can evoke EADs in human atrium. Indeed, phase-3 EADs have been shown to re-initiate atrial fibrillation (AF) during autonomic stimulation, which is a well-recognized initiator of AF.

METHODS

We integrated a Markov model of INa gating in our human atrial myocyte model. To simulate experimental results, we rapidly paced this cell model at 10Hz in the presence of 0.1μM acetylcholine and 1μM isoproterenol, and assessed EAD occurrence upon return to sinus rhythm (1Hz).

RESULTS

Cellular Ca(2+) loading during fast pacing results in a transient period of hypercontractility after return to sinus rhythm. Here, fast repolarization and enhanced NCX facilitate INa reactivation via the canonical gating mode (i.e., not late INa burst mode), which drives EAD initiation. Simulating ranolazine administration reduces atrial peak INa and leads to faster repolarization, during which INa fails to reactivate and EADs are prevented.

CONCLUSIONS

Non-equilibrium INa reactivation can critically contribute to arrhythmias, specifically in human atrial myocytes. Ranolazine might be beneficial in this context by blocking peak (not late) atrial INa.

Caveolae in ventricular myocytes are required for stretch-dependent conduction slowing

Mechanical stretch of cardiac muscle modulates action potential propagation velocity, causing potentially arrhythmogenic conduction slowing. The mechanisms by which stretch alters cardiac conduction remain unknown, but previous studies suggest that stretch can affect the conformation of caveolae in myocytes and other cell types. We tested the hypothesis that slowing of action potential conduction due to cardiac myocyte stretch is dependent on caveolae. Cardiac action potential propagation velocities, measured by optical mapping in isolated mouse hearts and in micropatterned mouse cardiomyocyte cultures, decreased reversibly with volume loading or stretch, respectively (by 19±5% and 26±4%). Stretch-dependent conduction slowing was not altered by stretch-activated channel blockade with gadolinium or by GsMTx-4 peptide, but was inhibited when caveolae were disrupted via genetic deletion of caveolin-3 (Cav3 KO) or membrane cholesterol depletion by methyl-β-cyclodextrin. In wild-type mouse hearts, stretch coincided with recruitment of caveolae to the sarcolemma, as observed by electron microscopy. In myocytes from wild-type but not Cav3 KO mice, stretch significantly increased cell membrane capacitance (by 98±64%), electrical time constant (by 285±149%), and lipid recruitment to the bilayer (by 84±39%). Recruitment of caveolae to the sarcolemma during physiologic cardiomyocyte stretch slows ventricular action potential propagation by increasing cell membrane capacitance.