Stem cell-derived cardiomyocyte heterogeneity confounds electrophysiological insights.

Human induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) offer potential as an in vitro model for studying drug cardiotoxicity and patient-specific cardiovascular disease. The inherent electrophysiological heterogeneity of these cells limits the depth of insights that can be drawn from well-designed experiments. In this review, we provide our perspective on some sources and the consequences of iPSC-CM heterogeneity. We demonstrate the extent of heterogeneity in the literature and explain how such heterogeneity is exacerbated by patch-clamp experimental artifacts in the manual and automated set-up. Finally, we discuss how this heterogeneity, caused by both intrinsic and extrinsic factors, limits our ability to build digital twins of patient-derived cardiomyocytes.

Single-cell ionic current phenotyping elucidates non-canonical features and predictive potential of cardiomyocytes during automated drug experiments.

All new drugs must go through preclinical screening tests to determine their proarrhythmic potential. While these assays effectively filter out dangerous drugs, they are too conservative, often misclassifying safe compounds as proarrhythmic. In this study, we attempt to address this shortcoming with a novel, medium-throughput drug-screening approach: we use an automated patch-clamp system to acquire optimized voltage clamp (VC) and action potential (AP) data from human induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) at several drug concentrations (baseline, 3×, 10× and 20× the effective free plasma concentrations). With our novel method, we show correlations between INa block and upstroke slowing after treatment with flecainide or quinine. Additionally, after quinine treatment, we identify significant reductions in current during voltage steps designed to isolate If and IKs. However, we do not detect any IKr block by either drug, and upon further investigation, do not see any IKr present in the iPSC-CMs when prepared for automated patch experiments (i.e. in suspension) - this is in contrast to similar experiments we have conducted with these cells using the manual patch setup. In this study, we: (1) present a proof-of-concept demonstration of a single-cell medium-throughput drug study, and (2) characterize the non-canonical electrophysiology of iPSC-CMs when prepared for experiments in a medium-throughput setting. KEY POINTS: Human induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) offer potential as an in vitro model to study the proarrhythmic potential of drugs, but insights from these cells are often limited by the low throughput of manual patch-clamp. In this study, we use a medium-throughput automated patch-clamp system to acquire action potential (AP) and complex voltage clamp (VC) data from single iPSC-CMs at multiple drug concentrations. A correlation between AP upstroke and INa transients was identified and drug-induced changes in ionic currents found. We also characterize the substantially altered physiology of iPSC-CMs when patched in an automated system, suggesting the need to investigate differences between manual and automated patch experiments.

Single-cell ionic current phenotyping explains stem cell-derived cardiomyocyte action potential morphology.

Human induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) are a promising tool to study arrhythmia-related factors, but the variability of action potential (AP) recordings from these cells limits their use as an in vitro model. In this study, we use recently published brief (10 s), dynamic voltage-clamp (VC) data to provide mechanistic insights into the ionic currents contributing to AP heterogeneity; we call this approach rapid ionic current phenotyping (RICP). Features of this VC data were correlated to AP recordings from the same cells, and we used computational models to generate mechanistic insights into cellular heterogeneity. This analysis uncovered several interesting links between AP morphology and ionic current density: both L-type calcium and sodium currents contribute to upstroke velocity, rapid delayed rectifier K+ current is the main determinant of the maximal diastolic potential, and an outward current in the activation range of slow delayed rectifier K+ is the main determinant of AP duration. Our analysis also identified an outward current in several cells at 6 mV that is not reproduced by iPSC-CM mathematical models but contributes to determining AP duration. RICP can be used to explain how cell-to-cell variability in ionic currents gives rise to AP heterogeneity. Because of its brief duration (10 s) and ease of data interpretation, we recommend the use of RICP for single-cell patch-clamp experiments that include the acquisition of APs.NEW & NOTEWORTHY We present rapid ionic current phenotyping (RICP), a current quantification approach based on an optimized voltage-clamp protocol. The method captures a rich snapshot of the ionic current dynamics, providing quantitative information about multiple currents (e.g., ICa,L, IKr) in the same cell. The protocol helped to identify key ionic determinants of cellular action potential heterogeneity in iPSC-CMs. This included unexpected results, such as the critical role of IKr in establishing the maximum diastolic potential.

Optimization of a cardiomyocyte model illuminates role of increased INa,L in repolarization reserve.

Cardiac ion currents may compensate for each other when one is compromised by a congenital or drug-induced defect. Such redundancy contributes to a robust repolarization reserve that can prevent the development of lethal arrhythmias. Most efforts made to describe this phenomenon have quantified contributions by individual ion currents. However, it is important to understand the interplay between all major ion-channel conductances, as repolarization reserve is dependent on the balance between all ion currents in a cardiomyocyte. Here, a genetic algorithm was designed to derive profiles of nine ion-channel conductances that optimize repolarization reserve in a mathematical cardiomyocyte model. Repolarization reserve was quantified using a previously defined metric, repolarization reserve current, i.e., the minimum constant current to prevent normal action potential repolarization in a cell. The optimization improved repolarization reserve current up to 84% compared to baseline in a human adult ventricular myocyte model and increased resistance to arrhythmogenic insult. The optimized conductance profiles were not only characterized by increased repolarizing current conductances but also uncovered a previously unreported behavior by the late sodium current. Simulations demonstrated that upregulated late sodium increased action potential duration, without compromising repolarization reserve current. The finding was generalized to multiple models. Ultimately, this computational approach, in which multiple currents were studied simultaneously, illuminated mechanistic insights into how the metric's magnitude could be increased and allowed for the unexpected role of late sodium to be elucidated.NEW & NOTEWORTHY Genetic algorithms are typically used to fit models or extract desired parameters from data. Here, we use the tool to produce a ventricular cardiomyocyte model with increased repolarization reserve. Since arrhythmia mitigation is dependent on multiple cardiac ion-channel conductances, study using a comprehensive, unbiased, and systems-level approach is important. The use of this optimization strategy allowed us to find robust profiles that illuminated unexpected mechanistic determinants of key ion-channel conductances in repolarization reserve.

Rapid ionic current phenotyping (RICP) identifies mechanistic underpinnings of iPSC-CM AP heterogeneity.

As a renewable, easily accessible, human-derived in vitro model, human induced pluripotent stem cell derived cardiomyocytes (iPSC-CMs) are a promising tool for studying arrhythmia-related factors, including cardiotoxicity and congenital proarrhythmia risks. An oft-mentioned limitation of iPSC-CMs is the abundant cell-to-cell variability in recordings of their electrical activity. Here, we develop a new method, rapid ionic current phenotyping (RICP), that utilizes a short (10 s) voltage clamp protocol to quantify cell-to-cell heterogeneity in key ionic currents. We correlate these ionic current dynamics to action potential recordings from the same cells and produce mechanistic insights into cellular heterogeneity. We present evidence that the L-type calcium current is the main determinant of upstroke velocity, rapid delayed rectifier K+ current is the main determinant of the maximal diastolic potential, and an outward current in the excitable range of slow delayed rectifier K+ is the main determinant of action potential duration. We measure an unidentified outward current in several cells at 6 mV that is not recapitulated by iPSC-CM mathematical models but contributes to determining action potential duration. In this way, our study both quantifies cell-to-cell variability in membrane potential and ionic currents, and demonstrates how the ionic current variability gives rise to action potential heterogeneity. Based on these results, we argue that iPSC-CM heterogeneity should not be viewed simply as a problem to be solved but as a model system to understand the mechanistic underpinnings of cellular variability.

Leak current, even with gigaohm seals, can cause misinterpretation of stem cell-derived cardiomyocyte action potential recordings.

AIMS: Human-induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) have become an essential tool to study arrhythmia mechanisms. Much of the foundational work on these cells, as well as the computational models built from the resultant data, has overlooked the contribution of seal-leak current on the immature and heterogeneous phenotype that has come to define these cells. The aim of this study is to understand the effect of seal-leak current on recordings of action potential (AP) morphology. METHODS AND RESULTS: Action potentials were recorded in human iPSC-CMs using patch clamp and simulated using previously published mathematical models. Our in silico and in vitro studies demonstrate how seal-leak current depolarizes APs, substantially affecting their morphology, even with seal resistances (Rseal) above 1 GΩ. We show that compensation of this leak current is difficult due to challenges with obtaining accurate measures of Rseal during an experiment. Using simulation, we show that Rseal measures (i) change during an experiment, invalidating the use of pre-rupture values, and (ii) are polluted by the presence of transmembrane currents at every voltage. Finally, we posit that the background sodium current in baseline iPSC-CM models imitates the effects of seal-leak current and is increased to a level that masks the effects of seal-leak current on iPSC-CMs. CONCLUSION: Based on these findings, we make recommendations to improve iPSC-CM AP data acquisition, interpretation, and model-building. Taking these recommendations into account will improve our understanding of iPSC-CM physiology and the descriptive ability of models built from such data.

Diversity of cells and signals in the cardiovascular system.

This white paper is the outcome of the seventh UC Davis Cardiovascular Research Symposium on Systems Approach to Understanding Cardiovascular Disease and Arrhythmia. This biannual meeting aims to bring together leading experts in subfields of cardiovascular biomedicine to focus on topics of importance to the field. The theme of the 2022 Symposium was 'Cell Diversity in the Cardiovascular System, cell-autonomous and cell-cell signalling'. Experts in the field contributed their experimental and mathematical modelling perspectives and discussed emerging questions, controversies, and challenges in examining cell and signal diversity, co-ordination and interrelationships involved in cardiovascular function. This paper originates from the topics of formal presentations and informal discussions from the Symposium, which aimed to develop a holistic view of how the multiple cell types in the cardiovascular system integrate to influence cardiovascular function, disease progression and therapeutic strategies. The first section describes the major cell types (e.g. cardiomyocytes, vascular smooth muscle and endothelial cells, fibroblasts, neurons, immune cells, etc.) and the signals involved in cardiovascular function. The second section emphasizes the complexity at the subcellular, cellular and system levels in the context of cardiovascular development, ageing and disease. Finally, the third section surveys the technological innovations that allow the interrogation of this diversity and advancing our understanding of the integrated cardiovascular function and dysfunction.

An in silico-in vitro pipeline for drug cardiotoxicity screening identifies ionic pro-arrhythmia mechanisms.

BACKGROUND AND PURPOSE: Before advancing to clinical trials, new drugs are screened for their pro-arrhythmic potential using a method that is overly conservative and provides limited mechanistic insight. The shortcomings of this approach can lead to the mis-classification of beneficial drugs as pro-arrhythmic. EXPERIMENTAL APPROACH: An in silico-in vitro pipeline was developed to circumvent these shortcomings. A computational human induced pluripotent stem cell-derived cardiomyocyte (iPSC-CM) model was used as part of a genetic algorithm to design experiments, specifically electrophysiological voltage clamp (VC) protocols, to identify which of several cardiac ion channels were blocked during in vitro drug studies. Such VC data, along with dynamically clamped action potentials (AP), were acquired from iPSC-CMs before and after treatment with a control solution or a low- (verapamil), intermediate- (cisapride or quinine) or high-risk (quinidine) drug. KEY RESULTS: Significant AP prolongation (a pro-arrhythmia marker) was seen in response to quinidine and quinine. The VC protocol identified block of IKr (a source of arrhythmias) by all strong IKr blockers, including cisapride, quinidine and quinine. The protocol also detected block of ICaL by verapamil and Ito by quinidine. Further demonstrating the power of the approach, the VC data uncovered a previously unidentified If block by quinine, which was confirmed with experiments using a HEK-293 expression system and automated patch-clamp. CONCLUSION AND IMPLICATIONS: We developed an in silico-in vitro pipeline that simultaneously identifies pro-arrhythmia risk and mechanism of ion channel-blocking drugs. The approach offers a new tool for evaluating cardiotoxicity during preclinical drug screening.

Constitutively Activating GNAS Somatic Mutation in Right Ventricular Outflow Tract Tachycardia.

[Figure: see text].

Validation of quantitative measure of repolarization reserve as a novel marker of drug induced proarrhythmia.

Repolarization reserve, the robustness of a cell to repolarize even when one of the repolarization mechanisms is failing, has been described qualitatively in terms of ionic currents, but has not been quantified by a generic metric that is applicable to drug screening. Prolonged repolarization leading to repolarization failure is highly arrhythmogenic. It may lead to ventricular tachycardia caused by triggered activity from early afterdepolarizations (EADs), or it may promote the occurrence of unidirectional conduction block and reentry. Both types of arrhythmia may deteriorate into ventricular fibrillation (VF) and death. We define the Repolarization Reserve Current (RRC) as the minimum constant current necessary to prevent normal repolarization of a cell. After developing and testing RRC for nine computational ionic models of various species, we applied it experimentally to atrial and ventricular human induced pluripotent stem cell-derived cardiomyocyte (hiPSC-CM), and isolated guinea-pig ventricular cardiomyocytes. In simulations, repolarization was all-or-none with a precise, model-dependent critical RRC, resulting in a discrete shift in the Action Potential Duration (APD) - RRC relation, in the occurrence of EADs and repolarization failure. These data were faithfully reproduced in cellular experiments. RRC allows simple, fast, unambiguous quantification of the arrhythmogenic propensity in cardiac cells of various origins and species without the need of prior knowledge of underlying currents and is suitable for high throughput applications, and personalized medicine applications.

Light-Activated Dynamic Clamp Using iPSC-Derived Cardiomyocytes.

iPSC-derived cardiomyocytes (iPSC-CMs) are a potentially advantageous platform for drug screening because they provide a renewable source of human cardiomyocytes. One obstacle to their implementation is their immature electrophysiology, which reduces relevance to adult arrhythmogenesis. To address this, dynamic clamp is used to inject current representing the insufficient potassium current, IK1, thereby producing more adult-like electrophysiology. However, dynamic clamp requires patch clamp and is therefore low throughput and ill-suited for large-scale drug screening. Here, we use optogenetics to generate such a dynamic-clamp current. The optical dynamic clamp (ODC) uses outward-current-generating opsin, ArchT, to mimic IK1, resulting in more adult-like action potential morphology, similar to IK1 injection via classic dynamic clamp. Furthermore, in the presence of an IKr blocker, ODC revealed expected action potential prolongation and reduced spontaneous excitation. The ODC presented here still requires an electrode to measure Vm but provides a first step toward contactless dynamic clamp, which will not only enable high-throughput screening but may also allow control within multicellular iPSC-CM formats to better recapitulate adult in vivo physiology.

Slow [Na+]i dynamics impacts arrhythmogenesis and spiral wave reentry in cardiac myocyte ionic model.

Accumulation of intracellular Na+ is gaining recognition as an important regulator of cardiac myocyte electrophysiology. The intracellular Na+ concentration can be an important determinant of the cardiac action potential duration, can modulate the tissue-level conduction of excitation waves, and can alter vulnerability to arrhythmias. Mathematical models of cardiac electrophysiology often incorporate a dynamic intracellular Na+ concentration, which changes much more slowly than the remaining variables. We investigated the dependence of several arrhythmogenesis-related factors on [Na+]i in a mathematical model of the human atrial action potential. In cell simulations, we found that [Na+]i accumulation stabilizes the action potential duration to variations in several conductances and that the slow dynamics of [Na+]i impacts bifurcations to pro-arrhythmic afterdepolarizations, causing intermittency between different rhythms. In long-lasting tissue simulations of spiral wave reentry, [Na+]i becomes spatially heterogeneous with a decreased area around the spiral wave rotation center. This heterogeneous region forms a functional anchor, resulting in diminished meandering of the spiral wave. Our findings suggest that slow, physiological, rate-dependent variations in [Na+]i may play complex roles in cellular and tissue-level cardiac dynamics.

Regularity of beating of small clusters of embryonic chick ventricular heart-cells: experiment vs. stochastic single-channel population model.

The transmembrane potential is recorded from small isopotential clusters of 2-4 embryonic chick ventricular cells spontaneously generating action potentials. We analyze the cycle-to-cycle fluctuations in the time between successive action potentials (the interbeat interval or IBI). We also convert an existing model of electrical activity in the cluster, which is formulated as a Hodgkin-Huxley-like deterministic system of nonlinear ordinary differential equations describing five individual ionic currents, into a stochastic model consisting of a population of ∼20 000 independently and randomly gating ionic channels, with the randomness being set by a real physical stochastic process (radio static). This stochastic model, implemented using the Clay-DeFelice algorithm, reproduces the fluctuations seen experimentally: e.g., the coefficient of variation (standard deviation/mean) of IBI is 4.3% in the model vs. the 3.9% average value of the 17 clusters studied. The model also replicates all but one of several other quantitative measures of the experimental results, including the power spectrum and correlation integral of the voltage, as well as the histogram, Poincaré plot, serial correlation coefficients, power spectrum, detrended fluctuation analysis, approximate entropy, and sample entropy of IBI. The channel noise from one particular ionic current (IKs), which has channel kinetics that are relatively slow compared to that of the other currents, makes the major contribution to the fluctuations in IBI. Reproduction of the experimental coefficient of variation of IBI by adding a Gaussian white noise-current into the deterministic model necessitates using an unrealistically high noise-current amplitude. Indeed, a major implication of the modelling results is that, given the wide range of time-scales over which the various species of channels open and close, only a cell-specific stochastic model that is formulated taking into consideration the widely different ranges in the frequency content of the channel-noise produced by the opening and closing of several different types of channels will be able to reproduce precisely the various effects due to membrane noise seen in a particular electrophysiological preparation.

Applications of Dynamic Clamp to Cardiac Arrhythmia Research: Role in Drug Target Discovery and Safety Pharmacology Testing.

Dynamic clamp, a hybrid-computational-experimental technique that has been used to elucidate ionic mechanisms underlying cardiac electrophysiology, is emerging as a promising tool in the discovery of potential anti-arrhythmic targets and in pharmacological safety testing. Through the injection of computationally simulated conductances into isolated cardiomyocytes in a real-time continuous loop, dynamic clamp has greatly expanded the capabilities of patch clamp outside traditional static voltage and current protocols. Recent applications include fine manipulation of injected artificial conductances to identify promising drug targets in the prevention of arrhythmia and the direct testing of model-based hypotheses. Furthermore, dynamic clamp has been used to enhance existing experimental models by addressing their intrinsic limitations, which increased predictive power in identifying pro-arrhythmic pharmacological compounds. Here, we review the recent advances of the dynamic clamp technique in cardiac electrophysiology with a focus on its future role in the development of safety testing and discovery of anti-arrhythmic drugs.

Global Optimization of Ventricular Myocyte Model to Multi-Variable Objective Improves Predictions of Drug-Induced Torsades de Pointes.

In silico cardiac myocyte models present powerful tools for drug safety testing and for predicting phenotypical consequences of ion channel mutations, but their accuracy is sometimes limited. For example, several models describing human ventricular electrophysiology perform poorly when simulating effects of long QT mutations. Model optimization represents one way of obtaining models with stronger predictive power. Using a recent human ventricular myocyte model, we demonstrate that model optimization to clinical long QT data, in conjunction with physiologically-based bounds on intracellular calcium and sodium concentrations, better constrains model parameters. To determine if the model optimized to congenital long QT data better predicts risk of drug-induced long QT arrhythmogenesis, in particular Torsades de Pointes risk, we tested the optimized model against a database of known arrhythmogenic and non-arrhythmogenic ion channel blockers. When doing so, the optimized model provided an improved risk assessment. In particular, we demonstrate an elimination of false-positive outcomes generated by the baseline model, in which simulations of non-torsadogenic drugs, in particular verapamil, predict action potential prolongation. Our results underscore the importance of currents beyond those directly impacted by a drug block in determining torsadogenic risk. Our study also highlights the need for rich data in cardiac myocyte model optimization and substantiates such optimization as a method to generate models with higher accuracy of predictions of drug-induced cardiotoxicity.

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

KEY POINTS: Arrhythmias result from disruptions to cardiac electrical activity, although the factors that control cellular action potentials are incompletely understood. We combined mathematical modelling with experiments in heart cells from guinea pigs to determine how cellular electrical activity is regulated. A mismatch between modelling predictions and the experimental results allowed us to construct an improved, more predictive mathematical model. The balance between two particular potassium currents dictates how heart cells respond to perturbations and their susceptibility to arrhythmias. ABSTRACT: Imbalances of ionic currents can destabilize the cardiac action potential and potentially trigger lethal cardiac arrhythmias. In the present study, we combined mathematical modelling with information-rich dynamic clamp experiments to determine the regulation of action potential morphology in guinea pig ventricular myocytes. Parameter sensitivity analysis was used to predict how changes in ionic currents alter action potential duration, and these were tested experimentally using dynamic clamp, a technique that allows for multiple perturbations to be tested in each cell. Surprisingly, we found that a leading mathematical model, developed with traditional approaches, systematically underestimated experimental responses to dynamic clamp perturbations. We then re-parameterized the model using a genetic algorithm, which allowed us to estimate ionic current levels in each of the cells studied. This unbiased model adjustment consistently predicted an increase in the rapid delayed rectifier K+ current and a drastic decrease in the slow delayed rectifier K+ current, and this prediction was validated experimentally. Subsequent simulations with the adjusted model generated the clinically relevant prediction that the slow delayed rectifier is better able to stabilize the action potential and suppress pro-arrhythmic events than the rapid delayed rectifier. In summary, iterative coupling of simulations and experiments enabled novel insight into how the balance between cardiac K+ currents influences ventricular arrhythmia susceptibility.

Illuminating Myocyte-Fibroblast Homotypic and Heterotypic Gap Junction Dynamics Using Dynamic Clamp.

Fibroblasts play a significant role in the development of electrical and mechanical dysfunction of the heart; however, the underlying mechanisms are only partially understood. One widely studied mechanism suggests that fibroblasts produce excess extracellular matrix, resulting in collagenous septa that slow propagation, cause zig-zag conduction paths, and decouple cardiomyocytes, resulting in a substrate for cardiac arrhythmia. An emerging mechanism suggests that fibroblasts promote arrhythmogenesis through direct electrical interactions with cardiomyocytes via gap junction (GJ) channels. In the heart, three major connexin (Cx) isoforms, Cx40, Cx43, and Cx45, form GJ channels in cell-type-specific combinations. Because each Cx is characterized by a unique time- and transjunctional voltage-dependent profile, we investigated whether the electrophysiological contributions of fibroblasts would vary with the specific composition of the myocyte-fibroblast (M-F) GJ channel. Due to the challenges of systematically modifying Cxs in vitro, we coupled native cardiomyocytes with in silico fibroblast and GJ channel electrophysiology models using the dynamic-clamp technique. We found that there is a reduction in the early peak of the junctional current during the upstroke of the action potential (AP) due to GJ channel gating. However, effects on the cardiomyocyte AP morphology were similar regardless of the specific type of GJ channel (homotypic Cx43 and Cx45, and heterotypic Cx43/Cx45 and Cx45/Cx43). To illuminate effects at the tissue level, we performed multiscale simulations of M-F coupling. First, we developed a cell-specific model of our dynamic-clamp experiments and investigated changes in the underlying membrane currents during M-F coupling. Second, we performed two-dimensional tissue sheet simulations of cardiac fibrosis and incorporated GJ channels in a cell type-specific manner. We determined that although GJ channel gating reduces junctional current, it does not significantly alter conduction velocity during cardiac fibrosis relative to static GJ coupling. These findings shed more light on the complex electrophysiological interplay between cardiac fibroblasts and myocytes.

Improving cardiomyocyte model fidelity and utility via dynamic electrophysiology protocols and optimization algorithms.

Mathematical models of cardiac electrophysiology are instrumental in determining mechanisms of cardiac arrhythmias. However, the foundation of a realistic multiscale heart model is only as strong as the underlying cell model. While there have been myriad advances in the improvement of cellular-level models, the identification of model parameters, such as ion channel conductances and rate constants, remains a challenging problem. The primary limitations to this process include: (1) such parameters are usually estimated from data recorded using standard electrophysiology voltage-clamp protocols that have not been developed with model building in mind, and (2) model parameters are typically tuned manually to subjectively match a desired output. Over the last decade, methods aimed at overcoming these disadvantages have emerged. These approaches include the use of optimization or fitting tools for parameter estimation and incorporating more extensive data for output matching. Here, we review recent advances in parameter estimation for cardiomyocyte models, focusing on the use of more complex electrophysiology protocols and global search heuristics. We also discuss future applications of such parameter identification, including development of cell-specific and patient-specific mathematical models to investigate arrhythmia mechanisms and predict therapy strategies.

Computational Approaches to Understanding the Role of Fibroblast-Myocyte Interactions in Cardiac Arrhythmogenesis.

The adult heart is composed of a dense network of cardiomyocytes surrounded by nonmyocytes, the most abundant of which are cardiac fibroblasts. Several cardiac diseases, such as myocardial infarction or dilated cardiomyopathy, are associated with an increased density of fibroblasts, that is, fibrosis. Fibroblasts play a significant role in the development of electrical and mechanical dysfunction of the heart; however the underlying mechanisms are only partially understood. One widely studied mechanism suggests that fibroblasts produce excess extracellular matrix, resulting in collagenous septa. These collagenous septa slow propagation, cause zig-zag conduction paths, and decouple cardiomyocytes resulting in a substrate for arrhythmia. Another emerging mechanism suggests that fibroblasts promote arrhythmogenesis through direct electrical interactions with cardiomyocytes via gap junctions. Due to the challenges of investigating fibroblast-myocyte coupling in native cardiac tissue, computational modeling and in vitro experiments have facilitated the investigation into the mechanisms underlying fibroblast-mediated changes in cardiomyocyte action potential morphology, conduction velocity, spontaneous excitability, and vulnerability to reentry. In this paper, we summarize the major findings of the existing computational studies investigating the implications of fibroblast-myocyte interactions in the normal and diseased heart. We then present investigations from our group into the potential role of voltage-dependent gap junctions in fibroblast-myocyte interactions.