Cardiac alternans induced by fibroblast-myocyte coupling: mechanistic insights from computational models

Fibroblast Density is a Risk Factor for Drug-induced Arrhythmias

A recent study by Kawatou et al . has shown that the local heterogeneity of ion channel conductance is a critical substrate for focal or reentrant arrhythmias. However, the role of fibroblasts with repolarization heterogeneity in the initiation and maintenance of arrhythmias remains unknown. In this study, we investigated how diffuse fibrosis contributes to the formation of focal and reentrant arrhythmias under drug-induced heterogeneity using physiologically detailed mathematical models of the human heart. To simulate drug-induced heterogeneity, we varied the maximum conductance of transmembrane potassium and calcium currents, leading to heterogeneity in action potential duration (APD). Then, we assessed the effects of different fibrosis densities (FD) on the occurrence of premature ventricular complexes (PVCs). Fibroblasts were randomly and evenly inserted into the tissue, and various FD levels ranging from 0 to 35% were examined. We found a biphasic relationship between FD and drug-induced PVCs. Within a certain range of FD, FD positively correlated with PVC susceptibility. However, excessively high fibrosis levels were associated with reduced susceptibility to PVCs. In addition, the self-sustainability of arrhythmias exhibited a positive correlation with FD. This study demonstrates the interplay between the diffuse fibrosis and the drug-induced heterogeneity of APD in the genesis of ventricular arrhythmias.

Author summary

Sudden cardiac death remains a leading cause of death worldwide. Understanding the mechanisms underlying arrhythmia and its precursors is critical for the development of effective therapies and drugs. Repolarization heterogeneity plays a crucial role in both the initiation and maintenance of arrhythmias. Fibroblasts constitute a vital component of cardiac structure, originating from the remodeling of ventricular wall cells or the transformation of injured myocardial cells. Fibroblasts are known to couple with and alter the electrical properties of myocardial cells. However, our understanding of the role of fibroblasts in the development of arrhythmia remains limited. In this study, we employed a physiologically detailed mathematical model of cardiac tissue to investigate the roles of drug-induced heterogeneity and diffuse fibrosis in the initiation and maintenance of arrhythmias. We used 2D and 3D computational models to simulate various levels of drug-induced heterogeneity conditions with normal to pathological levels of fibroblast density (FD). We found that within a certain range of FD, fibroblasts promote PVCs under drug-induced heterogeneity. However, if FD exceeds 30%, the occurrence of PVCs decreases (biphasic relationship). On the other hand, the self-sustainability of VF (ventricular fibrillation) consistently increases with FD. This study implies that fibroblasts in cardiac tissue may play different roles in the initiation and maintenance of arrhythmia.

Modeling Fibroblast-Cardiomyocyte Interactions: Unveiling the Role of Ion Currents in Action Potential Modulation

Fibrotic cardiomyopathy represents a significant pathological condition characterized by the interaction between cardiomyocytes and fibroblasts in the heart, and it currently lacks an effective cure. In vitro platforms, such as engineered heart tissue (EHT) developed through the co-culturing of cardiomyocytes and fibroblasts, are under investigation to elucidate and manipulate these cellular interactions. We present the first integration of mathematical electrophysiological models that encapsulate fibroblast-cardiomyocyte interactions with experimental EHT studies to identify and modulate the ion channels governing these dynamics. Our findings resolve a long-standing debate regarding the effect of fibroblast coupling on cardiomyocyte action potential duration (APD). We demonstrate that these seemingly contradictory outcomes are contingent upon the specific properties of the cardiomyocyte to which the fibroblast is coupled, particularly the relative magnitudes of the fast Na+ and transient outward K+ currents within the cardiomyocyte. Our results emphasize the critical importance of detailed ionic current representation in cardiomyocytes for accurately predicting the interactions between cardiomyocytes and fibroblasts in EHT. Surprisingly, complex ion channel-based models of fibroblast electrophysiology did not outperform simplified resistance-capacitance models in this analysis. Collectively, our findings highlight the promising potential of synergizing in vitro and in silico approaches to identify therapeutic targets for cardiomyopathies.

Atrial fibroblast-derived exosomal miR-21 upregulate myocardial KCa3.1 via the PI3K-Akt pathway during rapid pacing

Background

Fibroblast-derived exosomes can regulate the electrical remodeling of cardiomyocytes, and the intermediate-conductance calcium-activated potassium channel (KCa3.1) is important in atrial electrical remodeling. However, the underlying molecular mechanisms remain unclear. This study aimed to investigate the regulation of cardiac electrophysiology by exosomes linked to KCa3.1.

Methods

Atrial myocytes (AMs) and atrial fibroblasts were isolated from Sprague-Dawley suckling rats and cultured individually. The cellular atrial fibrillation (AF) model was established via electrical stimulation (1.0 v/cm, 10 Hz), and fibroblast-derived exosomes were isolated via ultracentrifugation. Exosomes were co-cultured with AMs to investigate their influences on KCa3.1 and the underlying mechanisms. Nanoparticle tracking analysis and transmission electron microscopy were used to measure exosome particle sizes and concentrations. Whole-cell patch clamp was applied to record the current density of KCa3.1 and action potential duration (APD). The expression of miR-21-5p was detected by reverse-transcription polymerase chain reaction (RT-PCR). Western blotting or immunofluorescence was used to measure the expression of exosomal markers, Akt phosphorylation, and KCa3.1.

Results

Rapid pacing promoted the secretion of exosomes from atrial fibroblasts and miR-21-5p expression in atrial fibroblasts and exosomes. KCa3.1 protein expression and current density significantly increased, and APD50 and APD90 were sharply shortened after rapid pacing in AMs. TRAM-34 (KCa3.1 blocker) extended APD and reduced susceptibility to AF. KCa3.1 and P-AKT expressions were amplified after co-culturing AMs with exosomes secreted by atrial fibroblasts. In contrast, the increase in KCa3.1 expression was reversed after the cells were co-cultured with exosomes secreted by atrial fibroblasts that were transfected with miR-21-5p inhibitors or after the use of LY294002, a PI3K/Akt pathway inhibitor.

Conclusions

Rapid pacing promoted the secretion of exosomes from fibroblasts, and miR-21-5p was upregulated in exosomes. Moreover, the miR-21-5p-enriched exosomes upregulated KCa3.1 expression in AMs via the PI3K/Akt pathway.

Fibroblast mediated dynamics in diffusively uncoupled myocytes: a simulation study using 2-cell motifs

In healthy hearts myocytes are typically coupled to nearest neighbours through gap junctions. Under pathological conditions such as fibrosis, or in scar tissue, or across ablation lines myocytes can uncouple from their neighbours. Electrical conduction may still occur via fibroblasts that not only couple proximal myocytes but can also couple otherwise unconnected regions. We hypothesise that such coupling can alter conduction between myocytes via introduction of delays or by initiation of premature stimuli that can potentially result in reentry or conduction blocks. To test this hypothesis we have developed several 2-cell motifs and investigated the effect of fibroblast mediated electrical coupling between uncoupled myocytes. We have identified various regimes of myocyte behaviour that depend on the strength of gap-junctional conductance, connection topology, and parameters of the myocyte and fibroblast models. These motifs are useful in developing a mechanistic understanding of long-distance coupling on myocyte dynamics and enable the characterisation of interaction between different features such as myocyte and fibroblast properties, coupling strengths and pacing period. They are computationally inexpensive and allow for incorporation of spatial effects such as conduction velocity. They provide a framework for constructing scar tissue boundaries and enable linking of cellular level interactions with scar induced arrhythmia.

Fibroblasts in heart scar tissue directly regulate cardiac excitability and arrhythmogenesis

After heart injury, dead heart muscle is replaced by scar tissue. Fibroblasts can electrically couple with myocytes, and changes in fibroblast membrane potential can lead to myocyte excitability, which suggests that fibroblast-myocyte coupling in scar tissue may be responsible for arrhythmogenesis. However, the physiologic relevance of electrical coupling of myocytes and fibroblasts and its impact on cardiac excitability in vivo have never been demonstrated. We genetically engineered a mouse that expresses the optogenetic cationic channel ChR2 (H134R) exclusively in cardiac fibroblasts. After myocardial infarction, optical stimulation of scar tissue elicited organ-wide cardiac excitation and induced arrhythmias in these animals. Complementing computational modeling with experimental approaches, we showed that gap junctional and ephaptic coupling, in a synergistic yet functionally redundant manner, excited myocytes coupled to fibroblasts.

Human myofibroblasts increase the arrhythmogenic potential of human induced pluripotent stem cell-derived cardiomyocytes

Human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) have the potential to remuscularize infarcted hearts but their arrhythmogenicity remains an obstacle to safe transplantation. Myofibroblasts are the predominant cell-type in the infarcted myocardium but their impact on transplanted hiPSC-CMs remains poorly defined. Here, we investigate the effect of myofibroblasts on hiPSC-CMs electrophysiology and Ca2+ handling using optical mapping of advanced human cell coculture systems mimicking cell-cell interaction modalities. Human myofibroblasts altered the electrophysiology and Ca2+ handling of hiPSC-CMs and downregulated mRNAs encoding voltage channels (KV4.3, KV11.1 and Kir6.2) and SERCA2a calcium pump. Interleukin-6 was elevated in the presence of myofibroblasts and direct stimulation of hiPSC-CMs with exogenous interleukin-6 recapitulated the paracrine effects of myofibroblasts. Blocking interleukin-6 reduced the effects of myofibroblasts only in the absence of physical contact between cell-types. Myofibroblast-specific connexin43 knockdown reduced functional changes in contact cocultures only when combined with interleukin-6 blockade. This provides the first in-depth investigation into how human myofibroblasts modulate hiPSC-CMs function, identifying interleukin-6 and connexin43 as paracrine- and contact-mediators respectively, and highlighting their potential as targets for reducing arrhythmic risk in cardiac cell therapy.

Nonmyocytes as electrophysiological contributors to cardiac excitation and conduction

Although cardiac action potential (AP) generation and propagation have traditionally been attributed exclusively to cardiomyocytes (CM), other cell types in the heart are also capable of forming electrically conducting junctions. Interactions between CM and nonmyocytes (NM) enable and modulate each other's activity. This review provides an overview of the current understanding of heterocellular electrical communication in the heart. Although cardiac fibroblasts were initially thought to be electrical insulators, recent studies have demonstrated that they form functional electrical connections with CM in situ. Other NM, such as macrophages, have also been recognized as contributing to cardiac electrophysiology and arrhythmogenesis. Novel experimental tools have enabled the investigation of cell-specific activity patterns in native cardiac tissue, which is expected to yield exciting new insights into the development of novel or improved diagnostic and therapeutic strategies.

Effects of fibroblast on electromechanical dynamics of human atrial tissue—insights from a 2D discrete element model

Roughly 75% of normal myocardial tissue volume is comprised of myocytes, however, fibroblasts by number are the most predominant cells in cardiac tissue. Previous studies have shown distinctive differences in cellular electrophysiology and excitability between myocytes and fibroblasts. However, it is still unclear how the electrical coupling between the two and the increased population of fibroblasts affects the electromechanical dynamics of cardiac tissue. This paper focuses on investigating effects of fibroblast-myocyte electrical coupling (FMEC) and fibroblast population on atrial electrical conduction and mechanical contractility by using a two-dimensional Discrete Element Method (DEM) model of cardiac tissue that is different to finite element method (FEM). In the model, the electro-mechanics of atrial cells are modelled by a biophysically detailed model for atrial electrical action potentials and myofilament kinetics, and the atrial fibroblasts are modelled by an active model that considers four active membrane ionic channel currents. Our simulation results show that the FMEC impairs myocytes’ electrical action potential and mechanical contractibility, manifested by reduced upstroke velocity, amplitude and duration of action potentials, as well as cell length shortening. At the tissue level, the FMEC slows down the conduction of excitation waves, and reduces strain of the tissue produced during a contraction course. These findings provide new insights into understandings of how FMEC impairs cardiac electrical and mechanical dynamics of the heart.

Electrophysiological Consequences of Cardiac Fibrosis

For both the atria and ventricles, fibrosis is generally recognized as one of the key determinants of conduction disturbances. By definition, fibrosis refers to an increased amount of fibrous tissue. However, fibrosis is not a singular entity. Various forms can be distinguished, that differ in distribution: replacement fibrosis, endomysial and perimysial fibrosis, and perivascular, endocardial, and epicardial fibrosis. These different forms typically result from diverging pathophysiological mechanisms and can have different consequences for conduction. The impact of fibrosis on propagation depends on exactly how the patterns of electrical connections between myocytes are altered. We will therefore first consider the normal patterns of electrical connections and their regional diversity as determinants of propagation. Subsequently, we will summarize current knowledge on how different forms of fibrosis lead to a loss of electrical connectivity in order to explain their effects on propagation and mechanisms of arrhythmogenesis, including ectopy, reentry, and alternans. Finally, we will discuss a histological quantification of fibrosis. Because of the different forms of fibrosis and their diverging effects on electrical propagation, the total amount of fibrosis is a poor indicator for the effect on conduction. Ideally, an assessment of cardiac fibrosis should exclude fibrous tissue that does not affect conduction and differentiate between the various types that do; in this article, we highlight practical solutions for histological analysis that meet these requirements.

Optogenetic actuation in ChR2-transduced fibroblasts alter excitation-contraction coupling and mechano-electric feedback in coupled cardiomyocytes: a computational modeling study

With the help of the conventional electrical method and the growing optogenetic technology, cardiac fibroblasts (Fbs) have been verified to couple electrically with working myocytes and bring electrophysiological remodeling changes in them. The intrinsic properties of cardiac functional autoregulation represented by excitation-contraction coupling (ECC) and mechano-electric feedback (MEF) have also been extensively studied. However, the roles of optogenetic stimulation on the characteristics of ECC and MEF in cardiomyocytes (CMs) coupled with Fbs have been barely investigated. In this study, we proposed a combined model composed of three modules to explore these influences. Simulation results showed that (1) during ECC, an increased light duration (LD) strengthened the inflow of ChR2 current and prolonged action potential duration (APD), and extended durations of twitch and internal sarcomere deformation through the decreased dissociation of calcium with troponin C (CaTnC) complexes and the prolonged duration of Xb attachment-detachment; (2) during MEF, an increased LD was followed by a longer muscle twitch and deformation, and led to APD prolongation through the inward ChR2 current and its inward rectification kinetics, which far outweighed the effects of the delaying dissociation of CaTnC complexes and the prolonged reverse mode of Na⁺-Ca²⁺ exchange on AP shortening; (3) due to the ChR2 current's rectification feature, enhancing the light irradiance (LI) brought slight variations in peak or valley values of electrophysiological and mechanical parameters while did not change durations of AP and twitch and muscle deformation in both ECC and MEF. In conclusion, the inward ChR2 current and its inward rectification feature were found to affect significantly the durations of AP and twitch in both ECC and MEF. The roles of optogenetic actuation on both ECC and MEF should be considered in future cardiac computational optogenetics at the tissue and organ scale.

CVAR-Seg: An Automated Signal Segmentation Pipeline for Conduction Velocity and Amplitude Restitution

BACKGROUND

Rate-varying S1S2 stimulation protocols can be used for restitution studies to characterize atrial substrate, ionic remodeling, and atrial fibrillation risk. Clinical restitution studies with numerous patients create large amounts of these data. Thus, an automated pipeline to evaluate clinically acquired S1S2 stimulation protocol data necessitates consistent, robust, reproducible, and precise evaluation of local activation times, electrogram amplitude, and conduction velocity. Here, we present the CVAR-Seg pipeline, developed focusing on three challenges: (i) No previous knowledge of the stimulation parameters is available, thus, arbitrary protocols are supported. (ii) The pipeline remains robust under different noise conditions. (iii) The pipeline supports segmentation of atrial activities in close temporal proximity to the stimulation artifact, which is challenging due to larger amplitude and slope of the stimulus compared to the atrial activity.

METHODS AND RESULTS

The S1 basic cycle length was estimated by time interval detection. Stimulation time windows were segmented by detecting synchronous peaks in different channels surpassing an amplitude threshold and identifying time intervals between detected stimuli. Elimination of the stimulation artifact by a matched filter allowed detection of local activation times in temporal proximity. A non-linear signal energy operator was used to segment periods of atrial activity. Geodesic and Euclidean inter electrode distances allowed approximation of conduction velocity. The automatic segmentation performance of the CVAR-Seg pipeline was evaluated on 37 synthetic datasets with decreasing signal-to-noise ratios. Noise was modeled by reconstructing the frequency spectrum of clinical noise. The pipeline retained a median local activation time error below a single sample (1 ms) for signal-to-noise ratios as low as 0 dB representing a high clinical noise level. As a proof of concept, the pipeline was tested on a CARTO case of a paroxysmal atrial fibrillation patient and yielded plausible restitution curves for conduction speed and amplitude.

CONCLUSION

The proposed openly available CVAR-Seg pipeline promises fast, fully automated, robust, and accurate evaluations of atrial signals even with low signal-to-noise ratios. This is achieved by solving the proximity problem of stimulation and atrial activity to enable standardized evaluation without introducing human bias for large data sets.

Mathematical modelling of the mechano-electric coupling in the human cardiomyocyte electrically connected with fibroblasts

Cardiac fibroblasts are interspersed within mammalian cardiac tissue. Fibroblasts are mechanically passive; however, they may communicate electrically with cardiomyocytes via gap junctions and thus affect the electrical and mechanical activity of myocytes. Several in-silico studies at both cellular (0D) and ventricular (3D) levels analysed the effects of fibroblasts on the myocardial electrical function. However, none of them addressed possible effects of fibroblast-myocyte electrical coupling to cardiomyocyte mechanical activity. In this paper, we propose a mathematical model for studying both electrical and mechanical responses of the human cardiomyocyte to its electrotonic interaction with cardiac fibroblasts. Our simulations have revealed that electrotonic interaction with fibroblasts affects not only the mechanical activity of the cardiomyocyte, comprising either moderate or significant reduction of contractility, but also the mechano-calcium and mechano-electric feedback loops, and all these effects are enhanced as the number of coupled fibroblasts is increased. Obtained results suggest that moderate values of the myocyte-fibroblast gap junction conductance (less than 1 nS) can be attributed to physiological conditions, contrasting to the higher values (2 nS and higher) proper rather for pathological situations (e.g. for infarct and/or border zones), since all mechanical indexes falls down dramatically in the case of such high conductance.

Human Cardiac Fibroblast Number and Activation State Modulate Electromechanical Function of hiPSC-Cardiomyocytes in Engineered Myocardium

Cardiac tissue engineering using hiPSC-derived cardiomyocytes is a promising avenue for cardiovascular regeneration, pharmaceutical drug development, cardiotoxicity evaluation, and disease modeling. Limitations to these applications still exist due in part to the need for more robust structural support, organization, and electromechanical function of engineered cardiac tissues. It is well accepted that heterotypic cellular interactions impact the phenotype of cardiomyocytes. The current study evaluates the functional effects of coculturing adult human cardiac fibroblasts (hCFs) in 3D engineered tissues on excitation and contraction with the goal of recapitulating healthy, nonarrhythmogenic myocardium in vitro. A small population (5% of total cell number) of hCFs in tissues improves tissue formation, material properties, and contractile function. However, two perturbations to the hCF population create disease-like phenotypes in engineered cardiac tissues. First, increasing the percentage of hCFs to 15% resulted in tissues with increased ectopic activity and spontaneous excitation rate. Second, hCFs undergo myofibroblast activation in traditional two-dimensional culture, and this altered phenotype ablated the functional benefits of hCFs when incorporated into engineered cardiac tissues. Taken together, the results of this study demonstrate that human cardiac fibroblast number and activation state modulate electromechanical function of hiPSC-cardiomyocytes and that a low percentage of quiescent hCFs are a valuable cell source to advance a healthy electromechanical response of engineered cardiac tissue for regenerative medicine applications.

Spatially Discordant Repolarization Alternans in the Absence of Conduction Velocity Restitution

Spatially discordant alternans (SDA) of action potential duration (APD) has been widely observed in cardiac tissue and is linked to cardiac arrhythmogenesis. Theoretical studies have shown that conduction velocity restitution (CVR) is required for the formation of SDA. However, this theory is not completely supported by experiments, indicating that other mechanisms may exist. In this study, we carried out computer simulations using mathematical models of action potentials to investigate the mechanisms of SDA in cardiac tissue. We show that when CVR is present and engaged, such as fast pacing from one side of the tissue, the spatial pattern of APD in the tissue undergoes either spatially concordant alternans or SDA, independent of initial conditions or tissue heterogeneities. When CVR is not engaged, such as simultaneous pacing of the whole tissue or under normal/slow heart rates, the spatial pattern of APD in the tissue can have multiple solutions, including spatially concordant alternans and different SDA patterns, depending on heterogeneous initial conditions or pre-existing repolarization heterogeneities. In homogeneous tissue, curved nodal lines are not stable, which either evolve into straight lines or disappear. However, in heterogeneous itssue, curved nodal lines can be stable, depending on their initial locations and shapes relative to the structure of the heterogeneity. Therefore, CVR-induced SDA and non-CVR-induced SDA exhibit different dynamical properties, which may be responsible for the different SDA properties observed in experimental studies and arrhythmogenesis in different clinical settings.

Mechanistic investigation of Ca2+ alternans in human heart failure and its modulation by fibroblasts

Background

Heart failure (HF) is characterized, among other factors, by a progressive loss of contractile function and by the formation of an arrhythmogenic substrate, both aspects partially related to intracellular Ca²⁺ cycling disorders. In failing hearts both electrophysiological and structural remodeling, including fibroblast proliferation, contribute to changes in Ca²⁺ handling which promote the appearance of Ca²⁺ alternans (Ca-alt). Ca-alt in turn give rise to repolarization alternans, which promote dispersion of repolarization and contribute to reentrant activity. The computational analysis of the incidence of Ca²⁺ and/or repolarization alternans under HF conditions in the presence of fibroblasts could provide a better understanding of the mechanisms leading to HF arrhythmias and contractile function disorders.

Methods and findings

The goal of the present study was to investigate in silico the mechanisms leading to the formation of Ca-alt in failing human ventricular myocytes and tissues with disperse fibroblast distributions. The contribution of ionic currents variability to alternans formation at the cellular level was analyzed and the results show that in normal ventricular tissue, altered Ca²⁺ dynamics lead to Ca-alt, which precede APD alternans and can be aggravated by the presence of fibroblasts. Electrophysiological remodeling of failing tissue alone is sufficient to develop alternans. The incidence of alternans is reduced when fibroblasts are present in failing tissue due to significantly depressed Ca²⁺ transients. The analysis of the underlying ionic mechanisms suggests that Ca-alt are driven by Ca²⁺-handling protein and Ca²⁺ cycling dysfunctions in the junctional sarcoplasmic reticulum and that their contribution to alternans occurrence depends on the cardiac remodeling conditions and on myocyte-fibroblast interactions.

Conclusion

It can thus be concluded that fibroblasts modulate the formation of Ca-alt in human ventricular tissue subjected to heart failure-related electrophysiological remodeling. Pharmacological therapies should thus consider the extent of both the electrophysiological and structural remodeling present in the failing heart.

Biomarkers Associated with Atrial Fibrosis and Remodeling

Atrial fibrillation is the most common rhythm disturbance encountered in clinical practice. Although often considered as solely arrhythmic in nature, current evidence has established that atrial myopathy constitutes both the substrate and the outcome of atrial fibrillation, thus initiating a vicious, self-perpetuating cycle. This myopathy is triggered by stress-induced (including pressure/volume overload, inflammation, oxidative stress) responses of atrial tissue, which in the long term become maladaptive, and combine elements of both structural, especially fibrosis, and electrical remodeling, with contemporary approaches yielding potentially useful biomarkers of these processes. Biomarker value becomes greater given the fact that they can both predict atrial fibrillation occurrence and treatment outcome. This mini-review will focus on the biomarkers of atrial remodeling (both electrical and structural) and fibrosis that have been validated in human studies, including biochemical, histological and imaging approaches.

Modeling cardiac complexity: Advancements in myocardial models and analytical techniques for physiological investigation and therapeutic development in vitro

Cardiomyopathies, heart failure, and arrhythmias or conduction blockages impact millions of patients worldwide and are associated with marked increases in sudden cardiac death, decline in the quality of life, and the induction of secondary pathologies. These pathologies stem from dysfunction in the contractile or conductive properties of the cardiomyocyte, which as a result is a focus of fundamental investigation, drug discovery and therapeutic development, and tissue engineering. All of these foci require in vitro myocardial models and experimental techniques to probe the physiological functions of the cardiomyocyte. In this review, we provide a detailed exploration of different cell models, disease modeling strategies, and tissue constructs used from basic to translational research. Furthermore, we highlight recent advancements in imaging, electrophysiology, metabolic measurements, and mechanical and contractile characterization modalities that are advancing our understanding of cardiomyocyte physiology. With this review, we aim to both provide a biological framework for engineers contributing to the field and demonstrate the technical basis and limitations underlying physiological measurement modalities for biologists attempting to take advantage of these state-of-the-art techniques.

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.

Ca2+ Cycling Impairment in Heart Failure Is Exacerbated by Fibrosis: Insights Gained From Mechanistic Simulations

Heart failure (HF) is characterized by altered Ca²⁺ cycling, resulting in cardiac contractile dysfunction. Failing myocytes undergo electrophysiological remodeling, which is known to be the main cause of abnormal Ca²⁺ homeostasis. However, structural remodeling, specifically proliferating fibroblasts coupled to myocytes in the failing heart, could also contribute to Ca²⁺ cycling impairment. The goal of the present study was to systematically analyze the mechanisms by which myocyte–fibroblast coupling could affect Ca²⁺ dynamics in normal conditions and in HF. Simulations of healthy and failing human myocytes were performed using established mathematical models, and cells were either isolated or coupled to fibroblasts. Univariate and multivariate sensitivity analyses were performed to quantify effects of ion transport pathways on biomarkers computed from intracellular [Ca²⁺] waveforms. Variability in ion channels and pumps was imposed and populations of models were analyzed to determine effects on Ca²⁺ dynamics. Our results suggest that both univariate and multivariate sensitivity analyses are valuable methodologies to shed light into the ionic mechanisms underlying Ca²⁺ impairment in HF, although differences between the two methodologies are observed at high parameter variability. These can result from either the fact that multivariate analyses take into account ion channels or non-linear effects of ion transport pathways on Ca²⁺ dynamics. Coupling either healthy or failing myocytes to fibroblasts decreased Ca²⁺ transients due to an indirect sink effect on action potential (AP) and thus on Ca²⁺ related currents. Simulations that investigated restoration of normal physiology in failing myocytes showed that Ca²⁺ cycling can be normalized by increasing SERCA and L-type Ca²⁺ current activity while decreasing Na⁺–Ca²⁺ exchange and SR Ca²⁺ leak. Changes required to normalize APs in failing myocytes depended on whether myocytes were coupled to fibroblasts. In conclusion, univariate and multivariate sensitivity analyses are helpful tools to understand how Ca²⁺ cycling is impaired in HF and how this can be exacerbated by coupling of myocytes to fibroblasts. The design of pharmacological actions to restore normal activity should take into account the degree of fibrosis in the failing heart.

Quantitative systems models illuminate arrhythmia mechanisms in heart failure: Role of the Na+ -Ca2+ -Ca2+ /calmodulin-dependent protein kinase II-reactive oxygen species feedback

Quantitative systems modeling aims to integrate knowledge in different research areas with models describing biological mechanisms and dynamics to gain a better understanding of complex clinical syndromes. Heart failure (HF) is a chronic complex cardiac disease that results from structural or functional disorders impairing the ability of the ventricle to fill with or eject blood. Highly interactive and dynamic changes in mechanical, structural, neurohumoral, metabolic, and electrophysiological properties collectively predispose the failing heart to cardiac arrhythmias, which are responsible for about a half of HF deaths. Multiscale cardiac modeling and simulation integrate structural and functional data from HF experimental models and patients to improve our mechanistic understanding of this complex arrhythmia syndrome. In particular, they allow investigating how disease-induced remodeling alters the coupling of electrophysiology, Ca²⁺ and Na⁺ handling, contraction, and energetics that lead to rhythm derangements. The Ca²⁺ /calmodulin-dependent protein kinase II, which expression and activity are enhanced in HF, emerges as a critical hub that modulates the feedbacks between these various subsystems and promotes arrhythmogenesis. This article is categorized under: Physiology > Mammalian Physiology in Health and Disease Models of Systems Properties and Processes > Mechanistic Models Models of Systems Properties and Processes > Cellular Models Models of Systems Properties and Processes > Organ, Tissue, and Physiological Models.

Directed fusion of cardiac spheroids into larger heterocellular microtissues enables investigation of cardiac action potential propagation via cardiac fibroblasts

Multicellular spheroids generated through cellular self-assembly provide cytoarchitectural complexities of native tissue including three-dimensionality, extensive cell-cell contacts, and appropriate cell-extracellular matrix interactions. They are increasingly suggested as building blocks for larger engineered tissues to achieve shapes, organization, heterogeneity, and other biomimetic complexities. Application of these tissue culture platforms is of particular importance in cardiac research as the myocardium is comprised of distinct but intermingled cell types. Here, we generated scaffold-free 3D cardiac microtissue spheroids comprised of cardiac myocytes (CMs) and/or cardiac fibroblasts (CFs) and used them as building blocks to form larger microtissues with different spatial distributions of CMs and CFs. Characterization of fusing homotypic and heterotypic spheroid pairs revealed an important influence of CFs on fusion kinetics, but most strikingly showed rapid fusion kinetics between heterotypic pairs consisting of one CF and one CM spheroid, indicating that CMs and CFs self-sort in vitro into the intermixed morphology found in the healthy myocardium. We then examined electrophysiological integration of fused homotypic and heterotypic microtissues by mapping action potential propagation. Heterocellular elongated microtissues which recapitulate the disproportionate CF spatial distribution seen in the infarcted myocardium showed that action potentials propagate through CF volumes albeit with significant delay. Complementary computational modeling revealed an important role of CF sodium currents and the spatial distribution of the CM-CF boundary in action potential conduction through CF volumes. Taken together, this study provides useful insights for the development of complex, heterocellular engineered 3D tissue constructs and their engraftment via tissue fusion and has implications for arrhythmogenesis in cardiac disease and repair.

Cardiac ischemia-insights from computational models

BACKGROUND

Complementary to clinical and experimental studies, computational cardiac modeling serves to obtain a comprehensive understanding of the cardiovascular system in order to analyze dysfunction, evaluate existing, and develop novel treatment strategies.

OBJECTIVES

We describe the basics of multiscale computational modeling of cardiac electrophysiology from the molecular ion channel to the whole body scale. By modeling cardiac ischemia, we illustrate how in silico experiments can contribute to our understanding of how the pathophysiological mechanisms translate into changes observed in diagnostic tools such as the electrocardiogram (ECG).

MATERIALS AND METHODS

Quantitative in silico modeling spans a wide range of scales from ion channel biophysics to ECG signals. For each of the scales, a set of mathematical equations describes electrophysiology in relation to the other scales. Integration of ischemia-induced changes is performed on the ion channel, single-cell, and tissue level. This approach allows us to study how effects simulated at molecular scales translate to changes in the ECG.

RESULTS

Ischemia induces action potential shortening and conduction slowing. Hence, ischemic myocardium has distinct and significant effects on propagation and repolarization of excitation, depending on the intramural extent of the ischemic region. For transmural and subendocardial ischemic regions, ST segment elevation and depression, respectively, were observed, whereas intermediate ischemic regions were found to be electrically silent (NSTEMI).

CONCLUSIONS

In silico modeling contributes quantitative and mechanistic insight into fundamental ischemia-related arrhythmogenic mechanisms. In addition, computational modeling can help to translate experimental findings at the (sub-)cellular level to the organ and body context (e. g., ECG), thereby providing a thorough understanding of this routinely used diagnostic tool that may translate into optimized applications.

Short-term action potential memory and electrical restitution: A cellular computational study on the stability of cardiac repolarization under dynamic pacing

Electrical restitution (ER) is a major determinant of repolarization stability and, under fast pacing rate, it reveals memory properties of the cardiac action potential (AP), whose dynamics have never been fully elucidated, nor their ionic mechanisms. Previous studies have looked at ER mainly in terms of changes in AP duration (APD) when the preceding diastolic interval (DI) changes and described dynamic conditions where this relationship shows hysteresis which, in turn, has been proposed as a marker of short-term AP memory and repolarization stability. By means of numerical simulations of a non-propagated human ventricular AP, we show here that measuring ER as APD versus the preceding cycle length (CL) provides additional information on repolarization dynamics which is not contained in the companion formulation. We focus particularly on fast pacing rate conditions with a beat-to-beat variable CL, where memory properties emerge from APD vs CL and not from APD vs DI and should thus be stored in APD and not in DI. We provide an ion-currents characterization of such conditions under periodic and random CL variability, and show that the memory stored in APD plays a stabilizing role on AP repolarization under pacing rate perturbations. The gating kinetics of L-type calcium current seems to be the main determinant of this safety mechanism. We also show that, at fast pacing rate and under otherwise identical pacing conditions, a periodically beat-to-beat changing CL is more effective than a random one in stabilizing repolarization. In summary, we propose a novel view of short-term AP memory, differentially stored between systole and diastole, which opens a number of methodological and theoretical implications for the understanding of arrhythmia development.

Reentry via high-frequency pacing in a mathematical model for human-ventricular cardiac tissue with a localized fibrotic region

Localized heterogeneities, caused by the regional proliferation of fibroblasts, occur in mammalian hearts because of diseases like myocardial infarction. Such fibroblast clumps can become sources of pathological reentrant activities, e.g., spiral or scroll waves of electrical activation in cardiac tissue. The occurrence of reentry in cardiac tissue with heterogeneities, such as fibroblast clumps, can depend on the frequency at which the medium is paced. Therefore, it is important to study the reentry-initiating potential of such fibroblast clumps at different frequencies of pacing. We investigate the arrhythmogenic effects of fibroblast clumps at high- and low-frequency pacing. We find that reentrant waves are induced in the medium more prominently at high-frequency pacing than with low-frequency pacing. We also study the other factors that affect the potential of fibroblast clumps to induce reentry in cardiac tissue. In particular, we show that the ability of a fibroblast clump to induce reentry depends on the size of the clump, the distribution and percentage of fibroblasts in the clump, and the excitability of the medium. We study the process of reentry in two-dimensional and a three-dimensional mathematical models for cardiac tissue.

Gq-activated fibroblasts induce cardiomyocyte action potential prolongation and automaticity in a three-dimensional microtissue environment

Cardiac fibroblasts (CFs) are known to regulate cardiomyocyte (CM) function in vivo and in two-dimensional in vitro cultures. This study examined the effect of CF activation on the regulation of CM electrical activity in a three-dimensional (3-D) microtissue environment. Using a scaffold-free 3-D platform with interspersed neonatal rat ventricular CMs and CFs, G_q-mediated signaling was selectively enhanced in CFs by Gα_q adenoviral infection before coseeding with CMs in nonadhesive hydrogels. After 3 days, the microtissues were analyzed by signaling assay, histological staining, quantitative PCR, Western blots, optical mapping with voltage- or Ca²⁺-sensitive dyes, and microelectrode recordings of CF resting membrane potential (RMP_CF). Enhanced G_q signaling in CFs increased microtissue size and profibrotic and prohypertrophic markers. Expression of constitutively active Gα_q in CFs prolonged CM action potential duration (by 33%) and rise time (by 31%), prolonged Ca²⁺ transient duration (by 98%) and rise time (by 65%), and caused abnormal electrical activity based on depolarization-induced automaticity. Constitutive G_q activation in CFs also depolarized RMP_CF from -33 to -20 mV and increased connexin 43 and connexin 45 expression. Computational modeling confers that elevated RMP_CF and increased cell-cell coupling between CMs and CFs in a 3-D environment could lead to automaticity. In conclusion, our data demonstrate that CF activation alone is capable of altering action potential and Ca²⁺ transient characteristics of CMs, leading to proarrhythmic electrical activity. Our results also emphasize the importance of a 3-D environment where cell-cell interactions are prevalent, underscoring that CF activation in 3-D tissue plays a significant role in modulating CM electrophysiology and arrhythmias.NEW & NOTEWORTHY In a three-dimensional microtissue model, which lowers baseline activation of cardiac fibroblasts but enables cell-cell, paracrine, and cell-extracellular matrix interactions, we demonstrate that selective cardiac fibroblast activation by enhanced G_q signaling, a pathophysiological trigger in the diseased heart, modulates cardiomyocyte electrical activity, leading to proarrhythmogenic automaticity.

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.

Stretch-Activated Current Can Promote or Suppress Cardiac Alternans Depending on Voltage-Calcium Interaction

Cardiac alternans has been linked to the onset of ventricular fibrillation and ventricular tachycardia, leading to life-threatening arrhythmias. Here, we investigated the effects of stretch-activated currents (ISAC) on alternans using a physiologically detailed model of the ventricular myocyte. We found that increasing ISAC suppresses alternans if the voltage-Ca coupling is positive or the alternans is voltage driven. However, for electromechanically discordant alternans, which occurs when the alternans is Ca driven with negative voltage-Ca coupling, increasing ISAC promotes Ca alternans. In addition, if action potential duration-Ca transients show quasiperiodicity, we observe a biphasic effect of ISAC, i.e., suppressing quasiperiodic oscillation at small stretch but promoting electromechanically discordant alternans at larger stretch. Our results demonstrate how ISAC interacts with coupled voltage-Ca dynamical systems with respect to alternans.

Computational modeling for cardiac safety pharmacology analysis: Contribution of fibroblasts

INTRODUCTION

Drug-induced proarrhythmic potential is an important regulatory criterion in safety pharmacology. The application of in silico approaches to predict proarrhythmic potential of new compounds is under consideration as part of future guidelines. Current approaches simulate the electrophysiology of a single human adult ventricular cardiomyocyte. However, drug-induced proarrhythmic potential can be different when cardiomyocytes are surrounded by non-muscle cells. Incorporating fibroblasts in models of myocardium is important particularly for predicting a drugs cardiac liability in the aging population - a growing population who take more medications and exhibit increased cardiac fibrosis. In this study, we used computational models to investigate the effects of fibroblast coupling on the electrophysiology and response to drugs of cardiomyocytes.

METHODS

A computational model of cardiomyocyte electrophysiology and ion handling (O'Hara, Virag, Varro, & Rudy, 2011) is coupled to a passive model of fibroblast electrophysiology to test the effects of three compounds that block cardiomyocyte ion channels. Results are compared to model results without fibroblast coupling to see how fibroblasts affect cardiomyocyte action potential duration at 90% repolarization (APD₉₀) and propensity for early after depolarization (EAD).

RESULTS

Simulation results show changes in cardiomyocyte APD₉₀ with increasing concentration of three drugs that affect cardiac function (dofetilide, vardenafil and nebivolol) when no fibroblasts are coupled to the cardiomyocyte. Coupling fibroblasts to cardiomyocytes markedly shortens APD₉₀. Moreover, increasing the number of fibroblasts can augment the shortening effect.

DISCUSSION

Coupling cardiomyocytes and fibroblasts are predicted to decrease proarrhythmic susceptibility under dofetilide, vardenafil and nebivolol block. However, this result is sensitive to parameters which define the electrophysiological function of the fibroblast. Fibroblast membrane capacitance and conductance (C_FB and G_FB) have less of an effect on APD₉₀ than the fibroblast resting membrane potential (E_FB). This study suggests that in both theoretical models and experimental tissue constructs that represent cardiac tissue, both cardiomyocytes and non-muscle cells should be considered when testing cardiac pharmacological agents.

Cardiac Fibroblast Activation Post-Myocardial Infarction: Current Knowledge Gaps

In response to myocardial infarction (MI), the wound healing response of the left ventricle (LV) comprises overlapping inflammatory, proliferative, and maturation phases, and the cardiac fibroblast is a key cell type involved in each phase. It has recently been appreciated that, early post-MI, fibroblasts transform to a proinflammatory phenotype and secrete cytokines and chemokines as well as matrix metalloproteinases (MMPs). Later post-MI, fibroblasts are activated to anti-inflammatory and proreparative phenotypes and generate anti-inflammatory and proangiogenic factors and extracellular matrix (ECM) components that form the infarct scar. Additional studies are needed to systematically examine how fibroblast activation shifts over the timeframe of the MI response and how modulation at different activation stages could alter wound healing and LV remodeling in distinct ways. This review summarizes current fibroblast knowledge as the foundation for a discussion of existing knowledge gaps.

Changes in Intracellular Na+ following Enhancement of Late Na+ Current in Virtual Human Ventricular Myocytes

The slowly inactivating or late Na+ current, INa-L, can contribute to the initiation of both atrial and ventricular rhythm disturbances in the human heart. However, the cellular and molecular mechanisms that underlie these pro-arrhythmic influences are not fully understood. At present, the major working hypothesis is that the Na+ influx corresponding to INa-L significantly increases intracellular Na+, [Na+]i; and the resulting reduction in the electrochemical driving force for Na+ reduces and (may reverse) Na+/Ca2+ exchange. These changes increase intracellular Ca2+, [Ca2+]i; which may further enhance INa-L due to calmodulin-dependent phosphorylation of the Na+ channels. This paper is based on mathematical simulations using the O'Hara et al (2011) model of baseline or healthy human ventricular action potential waveforms(s) and its [Ca2+]i homeostasis mechanisms. Somewhat surprisingly, our results reveal only very small changes (≤ 1.5 mM) in [Na+]i even when INa-L is increased 5-fold and steady-state stimulation rate is approximately 2 times the normal human heart rate (i.e. 2 Hz). Previous work done using well-established models of the rabbit and human ventricular action potential in heart failure settings also reported little or no change in [Na+]i when INa-L was increased. Based on our simulations, the major short-term effect of markedly augmenting INa-L is a significant prolongation of the action potential and an associated increase in the likelihood of reactivation of the L-type Ca2+ current, ICa-L. Furthermore, this action potential prolongation does not contribute to [Na+]i increase.

Myofibroblasts Electrotonically Coupled to Cardiomyocytes Alter Conduction: Insights at the Cellular Level from a Detailed In silico Tissue Structure Model

Fibrotic myocardial remodeling is typically accompanied by the appearance of myofibroblasts (MFBs). In vitro, MFBs were shown to slow conduction and precipitate ectopic activity following gap junctional coupling to cardiomyocytes (CMCs). To gain further mechanistic insights into this arrhythmogenic MFB-CMC crosstalk, we performed numerical simulations in cell-based high-resolution two-dimensional tissue models that replicated experimental conditions. Cell dimensions were determined using confocal microscopy of single and co-cultured neonatal rat ventricular CMCs and MFBs. Conduction was investigated as a function of MFB density in three distinct cellular tissue architectures: CMC strands with endogenous MFBs, CMC strands with coating MFBs of two different sizes, and CMC strands with MFB inserts. Simulations were performed to identify individual contributions of heterocellular gap junctional coupling and of the specific electrical phenotype of MFBs. With increasing MFB density, both endogenous and coating MFBs slowed conduction. At MFB densities of 5-30%, conduction slowing was most pronounced in strands with endogenous MFBs due to the MFB-dependent increase in axial resistance. At MFB densities >40%, very slow conduction and spontaneous activity was primarily due to MFB-induced CMC depolarization. Coating MFBs caused non-uniformities of resting membrane potential, which were more prominent with large than with small MFBs. In simulations of MFB inserts connecting two CMC strands, conduction delays increased with increasing insert lengths and block appeared for inserts >1.2 mm. Thus, electrophysiological properties of engineered CMC-MFB co-cultures depend on MFB density, MFB size and their specific positioning in respect to CMCs. These factors may influence conduction characteristics in the heterocellular myocardium.

Fibroblasts Slow Conduction Velocity in a Reconstituted Tissue Model of Fibrotic Cardiomyopathy

Myocardial function deteriorates over the course of fibrotic cardiomyopathy, due to electrophysiological and mechanical effects of myofibroblasts that are not completely understood. Although a range of experimental model systems and associated theoretical treatments exist at the levels of isolated cardiomyocytes and planar co-cultures of myofibroblasts and cardiomyocytes, interactions between these cell types at the tissue level are less clear. We studied these interactions through an engineered heart tissue (EHT) model of fibrotic myocardium and a mathematical model of the effects of cellular composition on EHT impulse conduction velocity. The EHT model allowed for modulation of cardiomyocyte and myofibroblast volume fractions, and observation of cell behavior in a three-dimensional environment that is more similar to native heart tissue than is planar cell culture. The cardiomyocyte and myofibroblast volume fractions determined the retardation of impulse conduction (spread of the action potential) in EHTs as measured by changes of the fluorescence of the Ca²⁺ probe, Fluo-2. Interpretation through our model showed retardation far in excess of predictions by homogenization theory, with conduction ceasing far below the fibroblast volume fraction associated with steric percolation. Results point to an important multiscale structural role of myofibroblasts in attenuating impulse conduction in fibrotic cardiomyopathy.

Modulation of cardiac pacemaker inter beat intervals by sinoatrial fibroblasts -a numerical study

The potential effect of sinoatrial fibroblasts on beat rate and variability of the cardiac pacemakers is not yet fully understood. Heterocellular coupling formation and fibroblast proliferation during diseased conditions may further signify the impact of those cells on sinoatrial node function. In this study we numerically modeled the impact of varying numbers of fibroblasts that are electrically coupled to a single pacemaker cell on several electrophysiological parameters. We employed cellular kinetics of the rabbit sinoatrial myocyte, and employed a range of potential gap junctional coupling between fibroblasts and myocytes. We show that increasing numbers of attached and coupled fibroblasts result in depolarization of the resting membrane potential of the pacemaker cell, as well as in attenuation in its action potential magnitude. We also demonstrate that the mean pacemaker inter-beat interval (IBI) was modulated in a non-linear, bi-phasic way by increasing numbers of attached fibroblasts, whereby an initial phase of decreasing IBIs was followed by a significant phase of exponentially increasing IBIs. These observations were more substantial for increased gap junctional coupling between the two cell types. We finally show that IBI variability exponentially increased with increasing numbers of attached and electrically coupled fibroblasts. Again, this effect was stronger with higher values of gap junctional coupling. We postulate that the last observation is related to the role of fibroblasts in amplifying membrane voltage fluctuations of attached myocytes.

The Multi-Domain Fibroblast/Myocyte Coupling in the Cardiac Tissue: A Theoretical Study

Cardiac fibroblast proliferation and concomitant collagenous matrix accumulation (fibrosis) develop during multiple cardiac pathologies. Recent studies have demonstrated direct electrical coupling between myocytes and fibroblasts in vitro, and assessed the electrophysiological implications of such coupling. However, in the living tissues, such coupling has not been demonstrated, and only indirect coupling via the extracellular space is likely to exist. In this study we employed a multi-domain model to assess the modulation of the cardiac electrophysiological properties by neighboring fibroblasts assuming only indirect coupling. Numerical simulations in 1D and 2D human atrial models showed that extracellular coupling sustains a significant impact on conduction velocity (CV) and a less significant effect on the action potential duration. Both CV and the slope of the CV restitution increased with increasing fibroblast density. This effect was more substantial for lower extracellular conductance. In 2D, spiral waves exhibited reduced frequency with increasing fibroblast density, and the propensity of wavebreaks and complex dynamics at high pacing rates significantly increased.

Computational modeling of cardiac fibroblasts and fibrosis

Altered fibroblast behavior can lead to pathologic changes in the heart such as arrhythmia, diastolic dysfunction, and systolic dysfunction. Computational models are increasingly used as a tool to identify potential mechanisms driving a phenotype or potential therapeutic targets against an unwanted phenotype. Here we review how computational models incorporating cardiac fibroblasts have clarified the role for these cells in electrical conduction and tissue remodeling in the heart. Models of fibroblast signaling networks have primarily focused on fibroblast cell lines or fibroblasts from other tissues rather than cardiac fibroblasts, specifically, but they are useful for understanding how fundamental signaling pathways control fibroblast phenotype. In the future, modeling cardiac fibroblast signaling, incorporating -omics and drug-interaction data into signaling network models, and utilizing multi-scale models will improve the ability of in silico studies to predict potential therapeutic targets against adverse cardiac fibroblast activity.

A device for co-culturing autonomic neurons and cardiomyocytes using micro-fabrication techniques

We developed a device for co-culturing sympathetic neurons, parasympathetic neurons, and cardiomyocytes using micro-fabrication techniques.

A Computational Study of the Factors Influencing the PVC-Triggering Ability of a Cluster of Early Afterdepolarization-Capable Myocytes

Premature ventricular complexes (PVCs), which are abnormal impulse propagations in cardiac tissue, can develop because of various reasons including early afterdepolarizations (EADs). We show how a cluster of EAD-generating cells (EAD clump) can lead to PVCs in a model of cardiac tissue, and also investigate the factors that assist such clumps in triggering PVCs. In particular, we study, through computer simulations, the effects of the following factors on the PVC-triggering ability of an EAD clump: (1) the repolarization reserve (RR) of the EAD cells; (2) the size of the EAD clump; (3) the coupling strength between the EAD cells in the clump; and (4) the presence of fibroblasts in the EAD clump. We find that, although a low value of RR is necessary to generate EADs and hence PVCs, a very low value of RR leads to low-amplitude EAD oscillations that decay with time and do not lead to PVCs. We demonstrate that a certain threshold size of the EAD clump, or a reduction in the coupling strength between the EAD cells, in the clump, is required to trigger PVCs. We illustrate how randomly distributed inexcitable obstacles, which we use to model collagen deposits, affect PVC-triggering by an EAD clump. We show that the gap-junctional coupling of fibroblasts with myocytes can either assist or impede the PVC-triggering ability of an EAD clump, depending on the resting membrane potential of the fibroblasts and the coupling strength between the myocyte and fibroblasts. We also find that the triggering of PVCs by an EAD clump depends sensitively on factors like the pacing cycle length and the distribution pattern of the fibroblasts.

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.

Cellular Communications in the Heart

Heart failure is one of the leading causes of morbidity and mortality worldwide. Cardiac remodelling is first an adaptive, becoming a maladaptive, compensatory mechanism that finally causes ventricular dysfunction independently of the etiology of the initial insult. In the present article the authors describe the elements of the human heart, examining their basic functions and their inter-communication under both normal and pathological circumstances. Cardiac myocytes carry out mechanical and electrical functions of the heart and cardiac fibroblasts maintain its structural integrity. Several factors can affect fibroblast activation and under pathological stress they transdifferentiate into myofibroblasts. Endothelial cells have complex biological functions, including the control of vascular permeability, vasomotion, regulation of haemostasis, immune responses and angiogenesis. The extracellular matrix is a complex architectural network consisting of a variety of proteins. Various routes using a plethora of products and mediators contribute to the cross-talk of the myocytes with endothelial cells, extracellular matrix and cardiac fibroblasts. A better understanding of the entire mechanism of cellular communication by the established or the more recently discovered agents will certainly emerge promising new perspectives when looking at the prevention of heart failure and leading to more substantial therapeutic interventions.

T-Wave Alternans Is Linked to Microvascular Obstruction and to Recurrent Coronary Ischemia After Myocardial Infarction

The purpose of this study is to investigate the relationship between T-wave alternans (TWA), infarct size and microvascular obstruction (MVO) and recurrent cardiac morbidity after ST elevation myocardial infarction (STEMI). One hundred six patients underwent TWA testing 1-12 months and 57 patients underwent cardiac magnetic resonance imaging (MRI) in the first 2-4 days after STEMI. During follow-up (3.5 ± 0.5 years), death (n = 2), ventricular tachycardia (n = 3), supraventricular tachycardia (n = 4), heart failure (n = 3) and recurrent coronary ischemia (n = 25) were observed. After multivariate analysis, positive TWA (HR2.59, CI1.10-6.11, p0.024) and larger MVO (HR1.08, CI1.01-1.16, p0.034) were associated with recurrent angina or ACS. Presence of MVO was correlated with TWA (Spearman rho 0.404, p0.002) and the impairment of LVEF (-0.524, p < 0.001). Patients after STEMI remain at a high risk of symptoms of coronary ischemia. The presence of MVO and TWA 1-12 months after STEMI is related to each other and to recurrent angina or ACS.

Turbulent states and their transitions in mathematical models for ventricular tissue: the effects of random interstitial fibroblasts

We study the dynamical behaviors of two types of spiral- and scroll-wave turbulence states, respectively, in two-dimensional (2D) and three-dimensional (3D) mathematical models, of human, ventricular, myocyte cells that are attached to randomly distributed interstitial fibroblasts; these turbulence states are promoted by (a) the steep slope of the action-potential-duration-restitution (APDR) plot or (b) early afterdepolarizations (EADs). Our single-cell study shows that (1) the myocyte-fibroblast (MF) coupling G_{j} and (2) the number N_{f} of fibroblasts in an MF unit lower the steepness of the APDR slope and eliminate the EAD behaviors of myocytes; we explore the pacing dependence of such EAD suppression. In our 2D simulations, we observe that a spiral-turbulence (ST) state evolves into a state with a single, rotating spiral (RS) if either (a) G_{j} is large or (b) the maximum possible number of fibroblasts per myocyte N_{f}^{max} is large. We also observe that the minimum value of G_{j}, for the transition from the ST to the RS state, decreases as N_{f}^{max} increases. We find that, for the steep-APDR-induced ST state, once the MF coupling suppresses ST, the rotation period of a spiral in the RS state increases as (1) G_{j} increases, with fixed N_{f}^{max}, and (2) N_{f}^{max} increases, with fixed G_{j}. We obtain the boundary between ST and RS stability regions in the N_{f}^{max}-G_{j} plane. In particular, for low values of N_{f}^{max}, the value of G_{j}, at the ST-RS boundary, depends on the realization of the randomly distributed fibroblasts; this dependence decreases as N_{f}^{max} increases. Our 3D studies show a similar transition from scroll-wave turbulence to a single, rotating, scroll-wave state because of the MF coupling. We examine the experimental implications of our study and propose that the suppression (a) of the steep slope of the APDR or (b) EADs can eliminate spiral- and scroll-wave turbulence in heterogeneous cardiac tissue, which has randomly distributed fibroblasts.

Fibroblast electrical remodeling in heart failure and potential effects on atrial fibrillation

Fibroblasts are activated in heart failure (HF) and produce fibrosis, which plays a role in maintaining atrial fibrillation (AF). The effect of HF on fibroblast ion currents and its potential role in AF are unknown. Here, we used a patch-clamp technique to investigate the effects of HF on atrial fibroblast ion currents, and mathematical computation to assess the potential impact of this remodeling on atrial electrophysiology and arrhythmogenesis. Atrial fibroblasts were isolated from control and tachypacing-induced HF dogs. Tetraethylammonium-sensitive voltage-gated fibroblast current (IKv,fb) was significantly downregulated (by ?44%), whereas the Ba(2+)-sensitive inward rectifier current (IKir,fb) was upregulated by 79%, in HF animals versus controls. The fibroblast resting membrane potential was hyperpolarized (?53 ± 2 mV vs. ?42 ± 2 mV in controls) and the capacitance was increased (29.7 ± 2.2 pF vs. 17.8 ± 1.4 pF in controls) in HF. These experimental findings were implemented in a mathematical model that included cardiomyocyte-fibroblast electrical coupling. IKir,fb upregulation had a profibrillatory effect through shortening of the action potential duration and hyperpolarization of the cardiomyocyte resting membrane potential. IKv,fb downregulation had the opposite electrophysiological effects and was antifibrillatory. Simulated pharmacological blockade of IKv,fb successfully terminated reentry under otherwise profibrillatory conditions. We conclude that HF induces fibroblast ion-current remodeling with IKv,fb downregulation and IKir,fb upregulation, and that, assuming cardiomyocyte-fibroblast electrical coupling, this remodeling has a potentially important effect on atrial electrophysiology and arrhythmogenesis, with the overall response depending on the balance of pro- and antifibrillatory contributions. These findings suggest that fibroblast K(+)-current remodeling is a novel component of AF-related remodeling that might contribute to arrhythmia dynamics.

Cardiac myocyte alternans in intact heart: Influence of cell-cell coupling and β-adrenergic stimulation

BACKGROUND

Cardiac alternans are proarrhythmic and mechanistically link cardiac mechanical dysfunction and sudden cardiac death. Beat-to-beat alternans occur when beats with large Ca(2+) transients and long action potential duration (APD) alternate with the converse. APD alternans are typically driven by Ca(2+) alternans and sarcoplasmic reticulum (SR) Ca(2+) release alternans. But the effect of intercellular communication via gap junctions (GJ) on alternans in the intact heart remains unknown.

OBJECTIVE

We assessed the effects of cell-to-cell coupling on local alternans in intact Langendorff-perfused mouse hearts, measuring single myocyte [Ca(2+)] alternans synchronization among neighboring cells, and effects of β-adrenergic receptor (β-AR) activation and reduced GJ coupling.

METHODS AND RESULTS

Mouse hearts (C57BL/6) were retrogradely perfused and loaded with Fluo8-AM to record cardiac myocyte [Ca(2+)] in situ with confocal microscopy. Single cell resolution allowed analysis of alternans within the intact organ during alternans induction. Carbenoxolone (25 μM), a GJ inhibitor, significantly increased the occurrence and amplitude of alternans in single cells within the intact heart. Alternans were concordant between neighboring cells throughout the field of view, except transiently during onset. β-AR stimulation only reduced Ca(2+) alternans in tissue that had reduced GJ coupling, matching effects seen in isolated myocytes.

CONCLUSIONS

Ca(2+) alternans among neighboring myocytes is predominantly concordant, likely because of electrical coupling between cells. Consistent with this, partial GJ uncoupling increased propensity and amplitude of Ca(2+) alternans, and made them more sensitive to reversal by β-AR activation, as in isolated myocytes. Electrical coupling between myocytes may thus limit the alternans initiation, but also allow alternans to be more stable once established.

Fibroblast proliferation alters cardiac excitation conduction and contraction: a computational study

In this study, the effects of cardiac fibroblast proliferation on cardiac electric excitation conduction and mechanical contraction were investigated using a proposed integrated myocardial-fibroblastic electromechanical model. At the cellular level, models of the human ventricular myocyte and fibroblast were modified to incorporate a model of cardiac mechanical contraction and cooperativity mechanisms. Cellular electromechanical coupling was realized with a calcium buffer. At the tissue level, electrical excitation conduction was coupled to an elastic mechanics model in which the finite difference method (FDM) was used to solve electrical excitation equations, and the finite element method (FEM) was used to solve mechanics equations. The electromechanical properties of the proposed integrated model were investigated in one or two dimensions under normal and ischemic pathological conditions. Fibroblast proliferation slowed wave propagation, induced a conduction block, decreased strains in the fibroblast proliferous tissue, and increased dispersions in depolarization, repolarization, and action potential duration (APD). It also distorted the wave-front, leading to the initiation and maintenance of re-entry, and resulted in a sustained contraction in the proliferous areas. This study demonstrated the important role that fibroblast proliferation plays in modulating cardiac electromechanical behaviour and which should be considered in planning future heart-modeling studies.

Mechanisms of ventricular arrhythmias: from molecular fluctuations to electrical turbulence

Ventricular arrhythmias have complex causes and mechanisms. Despite extensive investigation involving many clinical, experimental, and computational studies, effective biological therapeutics are still very limited. In this article, we review our current understanding of the mechanisms of ventricular arrhythmias by summarizing the state of knowledge spanning from the molecular scale to electrical wave behavior at the tissue and organ scales and how the complex nonlinear interactions integrate into the dynamics of arrhythmias in the heart. We discuss the challenges that we face in synthesizing these dynamics to develop safe and effective novel therapeutic approaches.

Nonlinear and Stochastic Dynamics in the Heart

In a normal human life span, the heart beats about 2 to 3 billion times. Under diseased conditions, a heart may lose its normal rhythm and degenerate suddenly into much faster and irregular rhythms, called arrhythmias, which may lead to sudden death. The transition from a normal rhythm to an arrhythmia is a transition from regular electrical wave conduction to irregular or turbulent wave conduction in the heart, and thus this medical problem is also a problem of physics and mathematics. In the last century, clinical, experimental, and theoretical studies have shown that dynamical theories play fundamental roles in understanding the mechanisms of the genesis of the normal heart rhythm as well as lethal arrhythmias. In this article, we summarize in detail the nonlinear and stochastic dynamics occurring in the heart and their links to normal cardiac functions and arrhythmias, providing a holistic view through integrating dynamics from the molecular (microscopic) scale, to the organelle (mesoscopic) scale, to the cellular, tissue, and organ (macroscopic) scales. We discuss what existing problems and challenges are waiting to be solved and how multi-scale mathematical modeling and nonlinear dynamics may be helpful for solving these problems.

Detection of abnormal cardiac activity using principal component analysis--a theoretical study

Electrogram-guided ablation has been recently developed for allowing better detection and localization of abnormal atrial activity that may be the source of arrhythmogeneity. Nevertheless, no clear indication for the benefit of using electrograms guided ablation over empirical ablation was established thus far, and there is a clear need of improving the localization of cardiac arrhythmogenic targets for ablation. In this paper, we propose a new approach for detection and localization of irregular cardiac activity during ablation procedures that is based on dimension reduction algorithms and principal component analysis (PCA). Using an 8×8 electrode array, our method produces manifolds that allow easy visualization and detection of possible arrhythmogenic ablation targets characterized by irregular conduction. We employ mathematical modeling and computer simulations to demonstrate the feasibility of the new approach for two well established arrhythmogenic sources for irregular conduction--spiral waves and patchy fibrosis. Our results show that the PCA method can differentiate between focal ectopic activity and spiral wave activity, as these two types of activity produce substantially different manifold shapes. Moreover, the technique allows the detection of spiral wave cores and their general meandering and drifting pattern. Fibrotic patches larger than 2 mm(2) could also be visualized using the PCA method, both for quiescent atrial tissue and for tissue exhibiting spiral wave activity. We envision that this method, contingent to further numerical and experimental validation studies in more complex, realistic geometrical configurations and with clinical data, can improve existing atrial ablation mapping capabilities, thus increasing success rates and optimizing arrhythmia management.