Properties of cardiac conduction in a cell-based computational model

Up-regulated novel-miR-17 promotes hypothermic reperfusion arrhythmias by negatively targeting Gja1 and mediating activation of the PKC/c-Jun signaling pathway

BACKGROUND

Hypothermic ischemia-reperfusion arrhythmia is a common complication of cardiothoracic surgery under cardiopulmonary bypass, but few studies have focused on this type of arrhythmia. Our prior study discovered reduced myocardial Cx43 protein levels may be linked to hypothermic reperfusion arrhythmias. However, more detailed molecular mechanism research is required.

METHOD

The microRNA and mRNA expression levels in myocardial tissues were detected by real-time quantitative PCR (RT-qPCR). Besides, the occurrence of hypothermic reperfusion arrhythmias and changes in myocardial electrical conduction were assessed by electrocardiography and ventricular epicardial activation mapping. Furthermore, bioinformatics analysis, applying antagonists of miRNA, western blotting, immunohistochemistry, a dual luciferase assay, and pearson correlation analysis were performed to investigate the underlying molecular mechanisms.

RESULTS

The expression level of novel-miR-17 was up-regulated in hypothermic ischemia-reperfusion myocardial tissues. Inhibition of novel-miR-17 upregulation ameliorated cardiomyocyte edema, reduced apoptosis, increased myocardial electrical conduction velocity, and shortened the duration of reperfusion arrhythmias. Mechanistic studies showed that novel-miR-17 reduced the expression of Cx43 by directly targeting Gja1 while mediating the activation of the PKC/c-Jun signaling pathway.

CONCLUSION

Up-regulated novel-miR-17 is a newly discovered pro-arrhythmic microRNA that may serve as a potential therapeutic target and biomarker for hypothermic reperfusion arrhythmias.

Evaluating computational efforts and physiological resolution of mathematical models of cardiac tissue

Computational techniques have significantly advanced our understanding of cardiac electrophysiology, yet they have predominantly concentrated on averaged models that do not represent the intricate dynamics near individual cardiomyocytes. Recently, accurate models representing individual cells have gained popularity, enabling analysis of the electrophysiology at the micrometer level. Here, we evaluate five mathematical models to determine their computational efficiency and physiological fidelity. Our findings reveal that cell-based models introduced in recent literature offer both efficiency and precision for simulating small tissue samples (comprising thousands of cardiomyocytes). Conversely, the traditional bidomain model and its simplified counterpart, the monodomain model, are more appropriate for larger tissue masses (encompassing millions to billions of cardiomyocytes). For simulations requiring detailed parameter variations along individual cell membranes, the EMI model emerges as the only viable choice. This model distinctively accounts for the extracellular (E), membrane (M), and intracellular (I) spaces, providing a comprehensive framework for detailed studies. Nonetheless, the EMI model's applicability to large-scale tissues is limited by its substantial computational demands for subcellular resolution.

A possible path to persistent re-entry waves at the outlet of the left pulmonary vein

Atrial fibrillation (AF) is the most common form of cardiac arrhythmia, often evolving from paroxysmal episodes to persistent stages over an extended timeframe. While various factors contribute to this progression, the precise biophysical mechanisms driving it remain unclear. Here we explore how rapid firing of cardiomyocytes at the outlet of the pulmonary vein of the left atria can create a substrate for a persistent re-entry wave. This is grounded in a recently formulated mathematical model of the regulation of calcium ion channel density by intracellular calcium concentration. According to the model, the number of calcium channels is controlled by the intracellular calcium concentration. In particular, if the concentration increases above a certain target level, the calcium current is weakened to restore the target level of calcium. During rapid pacing, the intracellular calcium concentration of the cardiomyocytes increases leading to a substantial reduction of the calcium current across the membrane of the myocytes, which again reduces the action potential duration. In a spatially resolved cell-based model of the outlet of the pulmonary vein of the left atria, we show that the reduced action potential duration can lead to re-entry. Initiated by rapid pacing, often stemming from paroxysmal AF episodes lasting several days, the reduction in calcium current is a critical factor. Our findings illustrate how such episodes can foster a conducive environment for persistent AF through electrical remodeling, characterized by diminished calcium currents. This underscores the importance of promptly addressing early AF episodes to prevent their progression to chronic stages.

New focus on cardiac voltage-gated sodium channel β1 and β1B: Novel targets for treating and understanding arrhythmias?

Voltage-gated sodium channels (VGSCs) are transmembrane protein complexes that are vital to the generation and propagation of action potentials in nerve and muscle fibers. The canonical VGSC is generally conceived as a heterotrimeric complex formed by 2 classes of membrane-spanning subunit: an α-subunit (pore forming) and 2 β-subunits (non-pore forming). NaV1.5 is the main sodium channel α-subunit of mammalian ventricle, with lower amounts of other α-subunits, including NaV1.6, being present. There are 4 β-subunits (β1-β4) encoded by 4 genes (SCN1B-SCN4B), each of which is expressed in cardiac tissues. Recent studies suggest that in addition to assignments in channel gating and trafficking, products of Scn1b may have novel roles in conduction of action potential in the heart and intracellular signaling. This includes evidence that the β-subunit extracellular amino-terminal domain facilitates adhesive interactions in intercalated discs and that its carboxyl-terminal region is a substrate for a regulated intramembrane proteolysis (RIP) signaling pathway, with a carboxyl-terminal peptide generated by β1 RIP trafficked to the nucleus and altering transcription of various genes, including NaV1.5. In addition to β1, the Scn1b gene encodes for an alternative splice variant, β1B, which contains an identical extracellular adhesion domain to β1 but has a unique carboxyl-terminus. Although β1B is generally understood to be a secreted variant, evidence indicates that when co-expressed with NaV1.5, it is maintained at the cell membrane, suggesting potential unique roles for this understudied protein. In this review, we focus on what is known of the 2 β-subunit variants encoded by Scn1b in heart, with particular focus on recent findings and the questions raised by this new information. We also explore data that indicate β1 and β1B may be attractive targets for novel antiarrhythmic therapeutics.

Adenosine monophosphate-regulated protein kinase inhibition modulates electrophysiological characteristics and calcium homeostasis of rabbit right ventricular outflow tract

BACKGROUND

Metabolic stress predisposes to ventricular arrhythmias and sudden cardiac death. Right ventricular outflow tract (RVOT) is the common origin of ventricular arrhythmias. Adenosine monophosphate-regulated protein kinase (AMPK) activation is an important compensatory mechanism for cardiac remodeling during metabolic stress.

OBJECTIVES

The purpose of this study was to access whether AMPK inhibition would modulate RVOT electrophysiology, calcium (Ca2+ ) regulation, and RVOT arrhythmogenesis or not.

METHODS

Conventional microelectrodes were used to record electrical activity before and after compound C (10 µM, an AMPK inhibitor) in isoproterenol (1 µM)-treated rabbit RVOT tissue preparations under electrical pacing. Whole-cell patch-clamp and confocal microscopic examinations were performed in baseline and compound C-treated rabbit RVOT cardiomyocytes to investigate ionic currents and intracellular Ca2+ transients in isolated rabbit RVOT cardiomyocytes.

RESULTS

Compound C decreased RVOT contractility, and reversed isoproterenol increased RVOT contractility. Compound C decreased the incidence, rate, and duration of isoproterenol-induced RVOT burst firing under rapid pacing. Compared to baseline, compound C-treated RVOT cardiomyocytes had a longer action potential duration, smaller intracellular Ca2+ transients, late sodium (Na+ ), peak L-type Ca2+ current density, Na+ -Ca2+ exchanger, transient outward potassium (K+ ) current, and rapid and slow delayed rectifier K+ currents.

CONCLUSION

AMPK inhibition modulates RVOT electrophysiological characteristics and Ca2+ homeostasis, contributing to lower RVOT arrhythmogenic activity. Accordingly, AMPK inhibition might potentially reduce ventricular tachyarrhythmias.

Myocardial ischemia simulation based on a multi-scale heart electrophysiology model

BACKGROUND

Myocardial ischemia, caused by insufficient myocardial blood supply, is a leading cause of human death worldwide. Therefore, it is crucial to prioritize the prevention and treatment of this condition. Mathematical modeling is a powerful technique for studying heart diseases.

OBJECTIVE

The aim of this study was to discuss the quantitative relationship between extracellular potassium concentration and the degree of myocardial ischemia directly related to it.

METHODS

A human cardiac electrophysiological multiscale model was developed to calculate action potentials of all cells simultaneously, enhancing efficiency over traditional reaction-diffusion models.

RESULTS

Contrary to the commonly held view that myocardial ischemia is caused by an increase in extracellular potassium concentration, our simulation results indicate that level 1 ischemia is associated with a decrease in extracellular potassium concentration.

CONCLUSION

This unusual finding provides a new perspective on the mechanisms underlying myocardial ischemia and has the potential to lead to the development of new diagnostic and treatment strategies.

The simplified Kirchhoff network model (SKNM): a cell-based reaction-diffusion model of excitable tissue

Cell-based models of excitable tissues offer the advantage of cell-level precision, which cannot be achieved using traditional homogenized electrophysiological models. However, this enhanced accuracy comes at the cost of increased computational demands, necessitating the development of efficient cell-based models. The widely-accepted bidomain model serves as the standard in computational cardiac electrophysiology, and under certain anisotropy ratio conditions, it is well known that it can be reduced to the simpler monodomain model. Recently, the Kirchhoff Network Model (KNM) was developed as a cell-based counterpart to the bidomain model. In this paper, we aim to demonstrate that KNM can be simplified using the same steps employed to derive the monodomain model from the bidomain model. We present the cell-based Simplified Kirchhoff Network Model (SKNM), which produces results closely aligned with those of KNM while requiring significantly less computational resources.

Efficient, cell-based simulations of cardiac electrophysiology; The Kirchhoff Network Model (KNM)

Mathematical models based on homogenized representation of cardiac tissue have greatly improved our understanding of cardiac electrophysiology. However, these models are too coarse to investigate the dynamics at the level of the myocytes since the cells are not present in homogenized models. Recently, fine scale models have been proposed to allow for cell-level resolution of the dynamics, but these models are too computationally expensive to be used in applications like whole heart simulations of large animals. To address this issue, we propose a model that balances computational demands and physiological accuracy. The model is founded on Kirchhoff's current law, and represents every myocyte in the tissue. This allows specific properties to be assigned to individual cardiomyocytes, and other cell types like fibroblasts can be added to the model in an accurate manner while keeping the computing efforts reasonable.

Nano-scale solution of the Poisson-Nernst-Planck (PNP) equations in a fraction of two neighboring cells reveals the magnitude of intercellular electrochemical waves

The basic building blocks of the electrophysiology of cardiomyocytes are ion channels integrated in the cell membranes. Close to the ion channels there are very strong electrical and chemical gradients. However, these gradients extend for only a few nano-meters and are therefore commonly ignored in mathematical models. The full complexity of the dynamics is modelled by the Poisson-Nernst-Planck (PNP) equations but these equations must be solved using temporal and spatial scales of nano-seconds and nano-meters. Here we report solutions of the PNP equations in a fraction of two abuttal cells separated by a tiny extracellular space. We show that when only the potassium channels of the two cells are open, a stationary solution is reached with the well-known Debye layer close to the membranes. When the sodium channels of one of the cells are opened, a very strong and brief electrochemical wave emanates from the channels. If the extracellular space is sufficiently small and the number of sodium channels is sufficiently high, the wave extends all the way over to the neighboring cell and may therefore explain cardiac conduction even at very low levels of gap junctional coupling.

Computational cardiac physiology for new modelers: Origins, foundations, and future

Mathematical models of the cardiovascular system have come a long way since they were first introduced in the early 19th century. Driven by a rapid development of experimental techniques, numerical methods, and computer hardware, detailed models that describe physical scales from the molecular level up to organs and organ systems have been derived and used for physiological research. Mathematical and computational models can be seen as condensed and quantitative formulations of extensive physiological knowledge and are used for formulating and testing hypotheses, interpreting and directing experimental research, and have contributed substantially to our understanding of cardiovascular physiology. However, in spite of the strengths of mathematics to precisely describe complex relationships and the obvious need for the mathematical and computational models to be informed by experimental data, there still exist considerable barriers between experimental and computational physiological research. In this review, we present a historical overview of the development of mathematical and computational models in cardiovascular physiology, including the current state of the art. We further argue why a tighter integration is needed between experimental and computational scientists in physiology, and point out important obstacles and challenges that must be overcome in order to fully realize the synergy of experimental and computational physiological research.

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.

Simulation of atrial fibrillation in a non-ohmic propagation model with dynamic gap junctions

Gap junctions exhibit nonlinear electrical properties that have been hypothesized to be relevant to arrhythmogenicity in a structurally remodeled tissue. Large-scale implementation of gap junction dynamics in 3D propagation models remains challenging. We aim to quantify the impact of nonlinear diffusion during episodes of arrhythmias simulated in a left atrial model. Homogenization of conduction properties in the presence of nonlinear gap junctions was performed by generalizing a previously developed mathematical framework. A monodomain model was solved in which conductivities were time-varying and depended on transjunctional potentials. Gap junction conductances were derived from a simplified Vogel-Weingart model with first-order gating and adjustable time constant. A bilayer interconnected cable model of the left atrium with 100  μm resolution was used. The diffusion matrix was recomputed at each time step according to the state of the gap junctions. Sinus rhythm and atrial fibrillation episodes were simulated in remodeled tissue substrates. Slow conduction was induced by reduced coupling and by diffuse or stringy fibrosis. Simulations starting from the same initial conditions were repeated with linear and nonlinear gap junctions. The discrepancy in activation times between the linear and nonlinear diffusion models was quantified. The results largely validated the linear approximation for conduction velocities >20 cm/s. In very slow conduction substrates, the discrepancy accumulated over time during atrial fibrillation, eventually leading to qualitative differences in propagation patterns, while keeping the descriptive statistics, such as cycle lengths, unchanged. The discrepancy growth rate was increased by impaired conduction, fibrosis, conduction heterogeneity, lateral uncoupling, fast gap junction time constant, and steeper action potential duration restitution.

Engineered bacterial voltage-gated sodium channel platform for cardiac gene therapy

Therapies for cardiac arrhythmias could greatly benefit from approaches to enhance electrical excitability and action potential conduction in the heart by stably overexpressing mammalian voltage-gated sodium channels. However, the large size of these channels precludes their incorporation into therapeutic viral vectors. Here, we report a platform utilizing small-size, codon-optimized engineered prokaryotic sodium channels (BacNav) driven by muscle-specific promoters that significantly enhance excitability and conduction in rat and human cardiomyocytes in vitro and adult cardiac tissues from multiple species in silico. We also show that the expression of BacNav significantly reduces occurrence of conduction block and reentrant arrhythmias in fibrotic cardiac cultures. Moreover, functional BacNav channels are stably expressed in healthy mouse hearts six weeks following intravenous injection of self-complementary adeno-associated virus (scAAV) without causing any adverse effects on cardiac electrophysiology. The large diversity of prokaryotic sodium channels and experimental-computational platform reported in this study should facilitate the development and evaluation of BacNav-based gene therapies for cardiac conduction disorders.

Deriving the Bidomain Model of Cardiac Electrophysiology From a Cell-Based Model; Properties and Comparisons

The bidomain model is considered to be the gold standard for numerical simulation of the electrophysiology of cardiac tissue. The model provides important insights into the conduction properties of the electrochemical wave traversing the cardiac muscle in every heartbeat. However, in normal resolution, the model represents the average over a large number of cardiomyocytes, and more accurate models based on representations of all individual cells have therefore been introduced in order to gain insight into the conduction properties close to the myocytes. The more accurate model considered here is referred to as the EMI model since both the extracellular space (E), the cell membrane (M) and the intracellular space (I) are explicitly represented in the model. Here, we show that the bidomain model can be derived from the cell-based EMI model and we thus reveal the close relation between the two models, and obtain an indication of the error introduced in the approximation. Also, we present numerical simulations comparing the results of the two models and thereby highlight both similarities and differences between the models. We observe that the deviations between the solutions of the models become larger for larger cell sizes. Furthermore, we observe that the bidomain model provides solutions that are very similar to the EMI model when conductive properties of the tissue are in the normal range, but large deviations are present when the resistance between cardiomyocytes is increased.

Cellular Size, Gap Junctions, and Sodium Channel Properties Govern Developmental Changes in Cardiac Conduction

Electrical conduction in cardiac ventricular tissue is regulated via sodium (Na⁺) channels and gap junctions (GJs). We and others have recently shown that Na⁺channels preferentially localize at the site of cell-cell junctions, the intercalated disc (ID), in adult cardiac tissue, facilitating coupling via the formation of intercellular Na⁺nanodomains, also termed ephaptic coupling (EpC). Several properties governing EpC vary with age, including Na⁺channel and GJ expression and distribution and cell size. Prior work has shown that neonatal cardiomyocytes have immature IDs with Na⁺channels and GJs diffusively distributed throughout the sarcolemma, while adult cells have mature IDs with preferentially localized Na⁺channels and GJs. In this study, we perform an in silico investigation of key age-dependent properties to determine developmental regulation of cardiac conduction. Simulations predict that conduction velocity (CV) biphasically depends on cell size, depending on the strength of GJ coupling. Total cell Na⁺channel conductance is predictive of CV in cardiac tissue with high GJ coupling, but not correlated with CV for low GJ coupling. We find that ephaptic effects are greatest for larger cells with low GJ coupling typically associated with intermediate developmental stages. Finally, simulations illustrate how variability in cellular properties during different developmental stages can result in a range of possible CV values, with a narrow range for both neonatal and adult myocardium but a much wider range for an intermediate developmental stage. Thus, we find that developmental changes predict associated changes in cardiac conduction.

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.

Localization of Na+ channel clusters in narrowed perinexi of gap junctions enhances cardiac impulse transmission via ephaptic coupling: a model study

Abstract

It has been proposed that when gap junctional coupling is reduced in cardiac tissue, action potential propagation can be supported via ephaptic coupling, a mechanism mediated by negative electric potentials occurring in narrow intercellular clefts of intercalated discs (IDs). Recent studies showed that sodium (Na⁺) channels form clusters near gap junction plaques in nanodomains called perinexi, where the ID cleft is even narrower. To examine the electrophysiological relevance of Na⁺ channel clusters being located in perinexi, we developed a 3D finite element model of two longitudinally abutting cardiomyocytes, with a central Na⁺ channel cluster on the ID membranes. When this cluster was located in the perinexus of a closely positioned gap junction plaque, varying perinexal width greatly modulated impulse transmission from one cell to the other, with narrow perinexi potentiating ephaptic coupling. This modulation occurred via the interplay of Na⁺ currents, extracellular potentials in the cleft and patterns of current flow within the cleft. In contrast, when the Na⁺ channel cluster was located remotely from the gap junction plaque, this modulation by perinexus width largely disappeared. Interestingly, the Na⁺ current in the ID membrane of the pre‐junctional cell switched from inward to outward during excitation, thus contributing ions to the activating channels on the post‐junctional ID membrane. In conclusion, these results indicate that the localization of Na⁺ channel clusters in the perinexi of gap junction plaques is crucial for ephaptic coupling, which is furthermore greatly modulated by perinexal width. These findings are relevant for a comprehensive understanding of cardiac excitation.

Key points

Ephaptic coupling is a cardiac conduction mechanism involving nanoscale‐level interactions between the sodium (Na⁺) current and the extracellular potential in narrow intercalated disc clefts.

When gap junctional coupling is reduced, ephaptic coupling acts in conjunction with the classical cardiac conduction mechanism based on gap junctional current flow.

In intercalated discs, Na⁺ channels form clusters that are preferentially located in the periphery of gap junction plaques, in nanodomains known as perinexi, but the electrophysiological role of these perinexi has never been examined.

In our new 3D finite element model of two cardiac cells abutting each other with their intercalated discs, a Na⁺ channel cluster located inside a narrowed perinexus facilitated impulse transmission via ephaptic coupling.

Our simulations demonstrate the role of narrowed perinexi as privileged sites for ephaptic coupling in pathological situations when gap junctional coupling is decreased.

Quantifying arrhythmic long QT effects of hydroxychloroquine and azithromycin with whole-heart optical mapping and simulations

Background

In March 2020, hydroxychloroquine (HCQ) alone or combined with azithromycin (AZM) was authorized as a treatment for COVID-19 in many countries. The therapy proved ineffective with long QT and deadly cardiac arrhythmia risks, illustrating challenges to determine the new safety profile of repurposed drugs.

Objective

To investigate proarrhythmic effects and mechanism of HCQ and AZM (combined and alone) with high doses of HCQ as in the COVID-19 clinical trials.

Methods

Proarrhythmic effects of HCQ and AZM are quantified using optical mapping with voltage-sensitive dyes in ex vivo Langendorff-perfused guinea pig (GP) hearts and with numerical simulations of a GP Luo-Rudy and a human O’Hara-Virag-Varro-Rudy models, for Epi, Endo, and M cells, in cell and tissue, incorporating the drug’s effect on cell membrane ionic currents.

Results

Experimentally, HCQ alone and combined with AZM leads to long QT intervals by prolonging the action potential duration and increased spatial dispersion of action potential (AP) repolarization across the heart, leading to proarrhythmic discordant alternans. AZM alone had a lesser arrhythmic effect with less triangulation of the AP shape. Mathematical cardiac models fail to reproduce most of the arrhythmic effects observed experimentally.

Conclusions

During public health crises, the risks and benefits of new and repurposed drugs could be better assessed with alternative experimental and computational approaches to identify proarrhythmic mechanisms. Optical mapping is an effective framework suitable to investigate the drug’s adverse effects on cardiac cell membrane ionic channels at the cellular level and arrhythmia mechanisms at the tissue and whole-organ level.

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.

It's clearly the heart! Optical transparency, cardiac tissue imaging, and computer modelling

Recent developments in clearing and microscopy enable 3D imaging with cellular resolution up to the whole organ level. These methods have been used extensively in neurobiology, but their uptake in other fields has been much more limited. Application of this approach to the human heart and effective use of the data acquired present challenges of scale and complexity. Four interlinked issues need to be addressed: 1) efficient clearing and labelling of heart tissue, 2) fast microscopic imaging of human-scale samples, 3) handling and processing of multi-terabyte 3D images, and 4) extraction of structural information in computationally tractable structure-based models of cardiac function. Preliminary studies show that each of these requirements can be achieved with the appropriate application and development of existing technologies.

Mechanisms underlying age-associated manifestation of cardiac sodium channel gain-of-function

Cardiac action potentials are initiated by sodium ion (Na+) influx through voltage-gated Na+ channels. Na+ channel gain-of-function (GOF) can arise in inherited conditions due to mutations in the gene encoding the cardiac Na+ channel, such as Long QT syndrome type 3 (LQT3). LQT3 can be a "concealed" disease, as patients with LQT3-associated mutations can remain asymptomatic until later in life; however, arrhythmias can also arise early in life in LQT3 patients, demonstrating a complex age-associated manifestation. We and others recently demonstrated that cardiac Na+ channels preferentially localize at the intercalated disc (ID) in adult cardiac tissue, which facilitates ephaptic coupling and formation of intercellular Na+ nanodomains that regulate pro-arrhythmic early afterdepolarization (EAD) formation in tissue with Na+ channel GOF. Several properties related to ephaptic coupling vary with age, such as cell size and Na+ channel and gap junction (GJ) expression and distribution: neonatal cells have immature IDs, with Na+ channels and GJs primarily diffusively distributed, while adult myocytes have mature IDs with preferentially localized Na+ channels and GJs. Here, we perform an in silico study varying critical age-dependent parameters to investigate mechanisms underlying age-associated manifestation of Na+ channel GOF in a model of guinea pig cardiac tissue. Simulations predict that total Na+ current conductance is a critical factor in action potential duration (APD) prolongation. We find a complex cell size/ Na+ channel expression relationship: increases in cell size (without concurrent increases in Na+ channel expression) suppress EAD formation, while increases in Na+ channel expression (without concurrent increases in cell size) promotes EAD formation. Finally, simulations with neonatal and early age-associated parameters predict normal APD with minimal dependence on intercellular cleft width; however, variability in cellular properties can lead to EADs presenting in early developmental stages. In contrast, for adult-associated parameters, EAD formation is highly dependent on cleft width, consistent with a mechanism underlying the age-associated manifestation of the Na+ channel GOF.

Intercellular Communication in the Heart: Therapeutic Opportunities for Cardiac Ischemia

The maintenance of tissue, organ, and organism homeostasis relies on an intricate network of players and mechanisms that assist in the different forms of cell-cell communication. Myocardial infarction, following heart ischemia and reperfusion, is associated with profound changes in key processes of intercellular communication, involving gap junctions, extracellular vesicles, and tunneling nanotubes, some of which have been implicated in communication defects associated with cardiac injury, namely arrhythmogenesis and progression into heart failure. Therefore, intercellular communication players have emerged as attractive powerful therapeutic targets aimed at preserving a fine-tuned crosstalk between the different cardiac cells in order to prevent or repair some of harmful consequences of heart ischemia and reperfusion, re-establishing myocardial function.

In vitro and In silico Models to Study SARS-CoV-2 Infection: Integrating Experimental and Computational Tools to Mimic "COVID-19 Cardiomyocyte"

The rapid dissemination of SARS-CoV-2 has made COVID-19 a tremendous social, economic, and health burden. Despite the efforts to understand the virus and treat the disease, many questions remain unanswered about COVID-19 mechanisms of infection and progression. Severe Acute Respiratory Syndrome (SARS) infection can affect several organs in the body including the heart, which can result in thromboembolism, myocardial injury, acute coronary syndromes, and arrhythmias. Numerous cardiac adverse events, from cardiomyocyte death to secondary effects caused by exaggerated immunological response against the virus, have been clinically reported. In addition to the disease itself, repurposing of treatments by using "off label" drugs can also contribute to cardiotoxicity. Over the past several decades, animal models and more recently, stem cell-derived cardiomyocytes have been proposed for studying diseases and testing treatments in vitro. In addition, mechanistic in silico models have been widely used for disease and drug studies. In these models, several characteristics such as gender, electrolyte imbalance, and comorbidities can be implemented to study pathophysiology of cardiac diseases and to predict cardiotoxicity of drug treatments. In this Mini Review, we (1) present the state of the art of in vitro and in silico cardiomyocyte modeling currently in use to study COVID-19, (2) review in vitro and in silico models that can be adopted to mimic the effects of SARS-CoV-2 infection on cardiac function, and (3) provide a perspective on how to combine some of these models to mimic "COVID-19 cardiomyocytes environment.".

Computational prediction of drug response in short QT syndrome type 1 based on measurements of compound effect in stem cell-derived cardiomyocytes

Short QT (SQT) syndrome is a genetic cardiac disorder characterized by an abbreviated QT interval of the patient's electrocardiogram. The syndrome is associated with increased risk of arrhythmia and sudden cardiac death and can arise from a number of ion channel mutations. Cardiomyocytes derived from induced pluripotent stem cells generated from SQT patients (SQT hiPSC-CMs) provide promising platforms for testing pharmacological treatments directly in human cardiac cells exhibiting mutations specific for the syndrome. However, a difficulty is posed by the relative immaturity of hiPSC-CMs, with the possibility that drug effects observed in SQT hiPSC-CMs could be very different from the corresponding drug effect in vivo. In this paper, we apply a multistep computational procedure for translating measured drug effects from these cells to human QT response. This process first detects drug effects on individual ion channels based on measurements of SQT hiPSC-CMs and then uses these results to estimate the drug effects on ventricular action potentials and QT intervals of adult SQT patients. We find that the procedure is able to identify IC50 values in line with measured values for the four drugs quinidine, ivabradine, ajmaline and mexiletine. In addition, the predicted effect of quinidine on the adult QT interval is in good agreement with measured effects of quinidine for adult patients. Consequently, the computational procedure appears to be a useful tool for helping predicting adult drug responses from pure in vitro measurements of patient derived cell lines.

Identifying Drug Response by Combining Measurements of the Membrane Potential, the Cytosolic Calcium Concentration, and the Extracellular Potential in Microphysiological Systems

Cardiomyocytes derived from human induced pluripotent stem cells (hiPSC-CMs) offer a new means to study and understand the human cardiac action potential, and can give key insight into how compounds may interact with important molecular pathways to destabilize the electrical function of the heart. Important features of the action potential can be readily measured using standard experimental techniques, such as the use of voltage sensitive dyes and fluorescent genetic reporters to estimate transmembrane potentials and cytosolic calcium concentrations. Using previously introduced computational procedures, such measurements can be used to estimate the current density of major ion channels present in hiPSC-CMs, and how compounds may alter their behavior. However, due to the limitations of optical recordings, resolving the sodium current remains difficult from these data. Here we show that if these optical measurements are complemented with observations of the extracellular potential using multi electrode arrays (MEAs), we can accurately estimate the current density of the sodium channels. This inversion of the sodium current relies on observation of the conduction velocity which turns out to be straightforwardly computed using measurements of extracellular waves across the electrodes. The combined data including the membrane potential, the cytosolic calcium concentration and the extracellular potential further opens up for the possibility of accurately estimating the effect of novel drugs applied to hiPSC-CMs.

Specific decreasing of Na+ channel expression on the lateral membrane of cardiomyocytes causes fatal arrhythmias in Brugada syndrome

Reduced cardiac sodium (Na+) channel current (INa) resulting from the loss-of-function of Na+ channel is a major cause of lethal arrhythmias in Brugada syndrome (BrS). Inspired by previous experimental studies which showed that in heart diseases INa was reduced along with expression changes in Na+ channel within myocytes, we hypothesized that the local decrease in INa caused by the alteration in Na+ channel expression in myocytes leads to the occurrence of phase-2 reentry, the major triggering mechanism of lethal arrhythmias in BrS. We constructed in silico human ventricular myocardial strand and ring models, and examined whether the Na+ channel expression changes in each myocyte cause the phase-2 reentry in BrS. Reducing Na+ channel expression in the lateral membrane of each myocyte caused not only the notch-and-dome but also loss-of-dome type action potentials and slowed conduction, both of which are typically observed in BrS patients. Furthermore, the selective reduction in Na+ channels on the lateral membrane of each myocyte together with spatial tissue heterogeneity of Na+ channel expression caused the phase-2 reentry and phase-2 reentry-mediated reentrant arrhythmias. Our data suggest that the BrS phenotype is strongly influenced by expression abnormalities as well as genetic abnormalities of Na+ channels.

Attenuating loss of cardiac conduction during no-flow ischemia through changes in perfusate sodium and calcium

Myocardial ischemia leads to conduction slowing, cell-to-cell uncoupling, and arrhythmias. We previously demonstrated that varying perfusate sodium (Na⁺) and calcium (Ca²⁺) attenuates conduction slowing and arrhythmias during simulated ischemia with continuous perfusion. Cardioprotection was selectively associated with widening of the perinexus, a gap junction adjacent nanodomain important to ephaptic coupling. It is unknown whether perfusate composition affects the perinexus or ischemic conduction during nonsimulated ischemia, when coronary flow is reduced or halted. We hypothesized that altering preischemic perfusate composition could facilitate perinexal expansion and attenuate conduction slowing during global ischemia. To test this hypothesis, ex vivo guinea pig hearts (n = 49) were Langendorff perfused with 145 or 153 mM Na⁺ and 1.25 or 2.0 mM Ca²⁺ and optically mapped during 30 min of no-flow ischemia. Altering Na⁺ and Ca²⁺ did not substantially affect baseline conduction. Increasing Na⁺ and decreasing Ca²⁺ both lowered pacing thresholds, whereas increasing Ca²⁺ narrowed perinexal width (W_p). A least squares mean estimate revealed that reduced perfusate Na⁺ and Ca²⁺ resulted in the most severe conduction slowing during ischemia. Increasing Na⁺ alone modestly attenuated conduction slowing, yet significantly delayed the median time to conduction block (10 to 16 min). Increasing both Na⁺ and Ca²⁺ selectively widened W_p during ischemia (22.7 vs. 15.7 nm) and attenuated conduction slowing to the greatest extent. Neither repolarization nor levels of total or phosphorylated connexin43 correlated with conduction slowing or block. Thus, perfusate-dependent widening of the perinexus preserved ischemic conduction and may be an adaptive response to ischemic stress.NEW & NOTEWORTHY Conduction slowing during acute ischemia creates an arrhythmogenic substrate. We have shown that extracellular ionic concentrations can alter conduction by modulating ephaptic coupling. Here, we demonstrate increased extracellular sodium and calcium significantly attenuate conduction slowing during no-flow ischemia. This effect was associated with selective widening of the perinexus, an intercalated disc nanodomain and putative cardiac ephapse. These findings suggest that acute changes in ephaptic coupling may serve as an adaptive response to ischemic stress.

Non-ohmic tissue conduction in cardiac electrophysiology: Upscaling the non-linear voltage-dependent conductance of gap junctions

Gap junctions are key mediators of intercellular communication in cardiac tissue, and their function is vital to sustaining normal cardiac electrical activity. Conduction through gap junctions strongly depends on the hemichannel arrangement and transjunctional voltage, rendering the intercellular conductance highly non-Ohmic, particularly under steady-state regimes of conduction. Despite this marked non-linear behavior, current tissue-level models of cardiac conduction are rooted in the assumption that gap-junctions conductance is constant (Ohmic), which results in inaccurate predictions of electrical propagation, particularly in the low junctional-coupling regime observed under pathological conditions. In this work, we present a novel non-Ohmic homogenization model (NOHM) of cardiac conduction that is suitable to tissue-scale simulations. Using non-linear homogenization theory, we develop a conductivity model that seamlessly upscales the voltage-dependent conductance of gap junctions, without the need of explicitly modeling gap junctions. The NOHM model allows for the simulation of electrical propagation in tissue-level cardiac domains that accurately resemble that of cell-based microscopic models for a wide range of junctional coupling scenarios, recovering key conduction features at a fraction of the computational complexity. A unique feature of the NOHM model is the possibility of upscaling the response of non-symmetric gap-junction conductance distributions, which result in conduction velocities that strongly depend on the direction of propagation, thus allowing to model the normal and retrograde conduction observed in certain regions of the heart. We envision that the NOHM model will enable organ-level simulations that are informed by sub- and inter-cellular mechanisms, delivering an accurate and predictive in-silico tool for understanding the heart function. Codes are available for download at https://github.com/dehurtado/NonOhmicConduction.

The heart relies on the propagation of electrical impulses that are mediated gap junctions, whose conduction properties vary depending on the transjunctional voltage. Despite this non-linear feature, current mathematical models assume that cardiac tissue behaves like an Ohmic (linear) material, thus delivering inaccurate results when simulated in a computer. Here we present a novel mathematical multiscale model that explicitly includes the non-Ohmic response of gap junctions in its predictions. Our results show that the proposed model recovers important conduction features modulated by gap junctions at a fraction of the computational complexity. This contribution represents an important step towards constructing computer models of a whole heart that can predict organ-level behavior in reasonable computing times.

Improved Computational Identification of Drug Response Using Optical Measurements of Human Stem Cell Derived Cardiomyocytes in Microphysiological Systems

Cardiomyocytes derived from human induced pluripotent stem cells (hiPSC-CMs) hold great potential for drug screening applications. However, their usefulness is limited by the relative immaturity of the cells' electrophysiological properties as compared to native cardiomyocytes in the adult human heart. In this work, we extend and improve on methodology to address this limitation, building on previously introduced computational procedures which predict drug effects for adult cells based on changes in optical measurements of action potentials and Ca2+ transients made in stem cell derived cardiac microtissues. This methodology quantifies ion channel changes through the inversion of data into a mathematical model, and maps this response to an adult phenotype through the assumption of functional invariance of fundamental intracellular and membrane channels during maturation. Here, we utilize an updated action potential model to represent both hiPSC-CMs and adult cardiomyocytes, apply an IC50-based model of dose-dependent drug effects, and introduce a continuation-based optimization algorithm for analysis of dose escalation measurements using five drugs with known effects. The improved methodology can identify drug induced changes more efficiently, and quantitate important metrics such as IC50 in line with published values. Consequently, the updated methodology is a step towards employing computational procedures to elucidate drug effects in adult cardiomyocytes for new drugs using stem cell-derived experimental tissues.

Modeling Depolarization Delay, Sodium Currents, and Electrical Potentials in Cardiac Transverse Tubules

T-tubules are invaginations of the lateral membrane of striated muscle cells that provide a large surface for ion channels and signaling proteins, thereby supporting excitation–contraction coupling. T-tubules are often remodeled in heart failure. To better understand the electrical behavior of T-tubules of cardiac cells in health and disease, this study addresses two largely unanswered questions regarding their electrical properties: (1) the delay of T-tubular membrane depolarization and (2) the effects of T-tubular sodium current on T-tubular potentials. Here, we present an elementary computational model to determine the delay in depolarization of deep T-tubular membrane segments as the narrow T-tubular lumen provides resistance against the extracellular current. We compare healthy tubules to tubules with constrictions and diseased tubules from mouse and human, and conclude that constrictions greatly delay T-tubular depolarization, while diseased T-tubules depolarize faster than healthy ones due to tubule widening. Increasing the tubule length non-linearly delays the depolarization. We moreover model the effect of T-tubular sodium current on intraluminal T-tubular potentials. We observe that extracellular potentials become negative during the sodium current transient (up to −40 mV in constricted T-tubules), which feedbacks on sodium channel function (self-attenuation) in a manner resembling ephaptic effects that have been described for intercalated discs where opposing membranes are very close together. The intraluminal potential and sodium current self-attenuation however greatly depend on sodium current conductance. These results show that (1) the changes in passive electrical properties of remodeled T-tubules cannot explain the excitation–contraction coupling defects in diseased cells; and (2) the sodium current may modulate intraluminal potentials. Such extracellular potentials might also affect excitation–contraction coupling and macroscopic conduction.