The stochastic nature of cardiac propagation at a microscopic level. Electrical description of myocardial architecture and its application to conduction

Distress-Mediated Remodeling of Cardiac Connexin-43 in a Novel Cell Model for Arrhythmogenic Heart Diseases

Gap junctions and their expression pattern are essential to robust function of intercellular communication and electrical propagation in cardiomyocytes. In healthy myocytes, the main cardiac gap junction protein connexin-43 (Cx43) is located at the intercalated disc providing a clear direction of signal spreading across the cardiac tissue. Dislocation of Cx43 to lateral membranes has been detected in numerous cardiac diseases leading to slowed conduction and high propensity for the development of arrhythmias. At the cellular level, arrhythmogenic diseases are associated with elevated levels of oxidative distress and gap junction remodeling affecting especially the amount and sarcolemmal distribution of Cx43 expression. So far, a mechanistic link between sustained oxidative distress and altered Cx43 expression has not yet been identified. Here, we propose a novel cell model based on murine induced-pluripotent stem cell-derived cardiomyocytes to investigate subcellular signaling pathways linking cardiomyocyte distress with gap junction remodeling. We tested the new hypothesis that chronic distress, induced by rapid pacing, leads to increased reactive oxygen species, which promotes expression of a micro-RNA, miR-1, specific for the control of Cx43. Our data demonstrate that Cx43 expression is highly sensitive to oxidative distress, leading to reduced expression. This effect can be efficiently prevented by the glutathione peroxidase mimetic ebselen. Moreover, Cx43 expression is tightly regulated by miR-1, which is activated by tachypacing-induced oxidative distress. In light of the high arrhythmogenic potential of altered Cx43 expression, we propose miR-1 as a novel target for pharmacological interventions to prevent the maladaptive remodeling processes during chronic distress in the heart.

Ephaptic Coupling Is a Mechanism of Conduction Reserve During Reduced Gap Junction Coupling

Many cardiac pathologies are associated with reduced gap junction (GJ) coupling, an important modulator of cardiac conduction velocity (CV). However, the relationship between phenotype and functional expression of the connexin GJ family of proteins is controversial. For example, a 50% reduction of GJ coupling has been shown to have little impact on myocardial CV due to a concept known as conduction reserve. This can be explained by the ephaptic coupling (EpC) theory whereby conduction is maintained by a combination of low GJ coupling and increased electrical fields generated in the sodium channel rich clefts between neighboring myocytes. At the same time, low GJ coupling may also increase intracellular charge accumulation within myocytes, resulting in a faster transmembrane potential rate of change during depolarization (dV/dt_max) that maintains macroscopic conduction. To provide insight into the prevalence of these two phenomena during pathological conditions, we investigated the relationship between EpC and charge accumulation within the setting of GJ remodeling using multicellular simulations and companion perfused mouse heart experiments. Conduction along a fiber of myocardial cells was simulated for a range of GJ conditions. The model incorporated intercellular variations, including GJ coupling conductance and distribution, cell-to-cell separation in the intercalated disc (perinexal width—W_P), and variations in sodium channel distribution. Perfused heart studies having conditions analogous to those of the simulations were performed using wild type mice and mice heterozygous null for the connexin gene Gja1. With insight from simulations, the relative contributions of EpC and charge accumulation on action potential parameters and conduction velocities were analyzed. Both simulation and experimental results support a common conclusion that low GJ coupling decreases and narrowing W_P increases the rate of the AP upstroke when sodium channels are densely expressed at the ends of myocytes, indicating that conduction reserve is more dependent on EpC than charge accumulation during GJ uncoupling.

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.

Multiparameter optimization of nonuniform passive diffusion properties for creating coarse-grained equivalent models of cardiac propagation

The arrhythmogenic role of discrete cardiac propagation may be assessed by comparing discrete (fine-grained) and equivalent continuous (coarse-grained) models. We aim to develop an optimization algorithm for estimating the smooth conductivity field that best reproduces the diffusion properties of a given discrete model. Our algorithm iteratively adjusts local conductivity of the coarse-grained continuous model by simulating passive diffusion from white noise initial conditions during 3-10 ms and computing the root mean square error with respect to the discrete model. The coarse-grained conductivity field was interpolated from up to 300 evenly spaced control points. We derived an approximate formula for the gradient of the cost function that required (in two dimensions) only two additional simulations per iteration regardless of the number of estimated parameters. Conjugate gradient solver facilitated simultaneous optimization of multiple conductivity parameters. The method was tested in rectangular anisotropic tissues with uniform and nonuniform conductivity (slow regions with sinusoidal profile) and random diffuse fibrosis, as well as in a monolayer interconnected cable model of the left atrium with spatially-varying fibrosis density. Comparison of activation maps served as validation. The results showed that after convergence the errors in activation time were < 1 ms for rectangular geometries and 1-3 ms in the atrial model. Our approach based on the comparison of passive properties (<10 ms simulation) avoids performing active propagation simulations (>100 ms) at each iteration while reproducing activation maps, with possible applications to investigating the impact of microstructure on arrhythmias.

Coupling between cardiac cells-An important determinant of electrical impulse propagation and arrhythmogenesis

Cardiac arrhythmias are an important cause of sudden cardiac death-a devastating manifestation of many underlying causes, such as heart failure and ischemic heart disease leading to ventricular tachyarrhythmias and ventricular fibrillation, and atrial fibrillation causing cerebral embolism. Cardiac electrical propagation is a main factor in the initiation and maintenance of cardiac arrhythmias. In the heart, gap junctions are the basic unit at the cellular level that host intercellular low-resistance channels for the diffusion of ions and small regulatory molecules. The dual voltage clamp technique enabled the direct measurement of electrical conductance between cells and recording of single gap junction channel openings. The rapid turnover of gap junction channels at the intercalated disk implicates a highly dynamic process of trafficking and internalization of gap junction connexons. Recently, non-canonical roles of gap junction proteins have been discovered in mitochondria function, cytoskeletal organization, trafficking, and cardiac rescue. At the tissue level, we explain the concepts of linear propagation and safety factor based on the model of linear cellular structure. Working myocardium is adequately represented as a discontinuous cellular network characterized by cellular anisotropy and connective tissue heterogeneity. Electrical propagation in discontinuous cellular networks reflects an interplay of three main factors: cell-to-cell electrical coupling, flow of electrical charge through the ion channels, and the microscopic tissue structure. This review provides a state-of-the-art update of the cardiac gap junction channels and their role in cardiac electrical impulse propagation and highlights a combined approach of genetics, cell biology, and physics in modern cardiac electrophysiology.

Cells, Materials, and Fabrication Processes for Cardiac Tissue Engineering

Cardiovascular disease is the number one killer worldwide, with myocardial infarction (MI) responsible for approximately 1 in 6 deaths. The lack of endogenous regenerative capacity, added to the deleterious remodelling programme set into motion by myocardial necrosis, turns MI into a progressively debilitating disease, which current pharmacological therapy cannot halt. The advent of Regenerative Therapies over 2 decades ago kick-started a whole new scientific field whose aim was to prevent or even reverse the pathological processes of MI. As a highly dynamic organ, the heart displays a tight association between 3D structure and function, with the non-cellular components, mainly the cardiac extracellular matrix (ECM), playing both fundamental active and passive roles. Tissue engineering aims to reproduce this tissue architecture and function in order to fabricate replicas able to mimic or even substitute damaged organs. Recent advances in cell reprogramming and refinement of methods for additive manufacturing have played a critical role in the development of clinically relevant engineered cardiovascular tissues. This review focuses on the generation of human cardiac tissues for therapy, paying special attention to human pluripotent stem cells and their derivatives. We provide a perspective on progress in regenerative medicine from the early stages of cell therapy to the present day, as well as an overview of cellular processes, materials and fabrication strategies currently under investigation. Finally, we summarise current clinical applications and reflect on the most urgent needs and gaps to be filled for efficient translation to the clinical arena.

Computational modeling of cardiac growth and remodeling in pressure overloaded hearts-Linking microstructure to organ phenotype

Cardiac growth and remodeling (G&R) refers to structural changes in myocardial tissue in response to chronic alterations in loading conditions. One such condition is pressure overload where elevated wall stresses stimulate the growth in cardiomyocyte thickness, associated with a phenotype of concentric hypertrophy at the organ scale, and promote fibrosis. The initial hypertrophic response can be considered adaptive and beneficial by favoring myocyte survival, but over time if pressure overload conditions persist, maladaptive mechanisms favoring cell death and fibrosis start to dominate, ultimately mediating the transition towards an overt heart failure phenotype. The underlying mechanisms linking biological factors at the myocyte level to biomechanical factors at the systemic and organ level remain poorly understood. Computational models of G&R show high promise as a unique framework for providing a quantitative link between myocardial stresses and strains at the organ scale to biological regulatory processes at the cellular level which govern the hypertrophic response. However, microstructurally motivated, rigorously validated computational models of G&R are still in their infancy. This article provides an overview of the current state-of-the-art of computational models to study cardiac G&R. The microstructure and mechanosensing/mechanotransduction within cells of the myocardium is discussed and quantitative data from previous experimental and clinical studies is summarized. We conclude with a discussion of major challenges and possible directions of future research that can advance the current state of cardiac G&R computational modeling. STATEMENT OF SIGNIFICANCE: The mechanistic links between organ-scale biomechanics and biological factors at the cellular size scale remain poorly understood as these are largely elusive to investigations using experimental methodology alone. Computational G&R models show high promise to establish quantitative links which allow more mechanistic insight into adaptation mechanisms and may be used as a tool for stratifying the state and predict the progression of disease in the clinic. This review provides a comprehensive overview of research in this domain including a summary of experimental data. Thus, this study may serve as a basis for the further development of more advanced G&R models which are suitable for making clinical predictions on disease progression or for testing hypotheses on pathogenic mechanisms using in-silico models.

Properties of cardiac conduction in a cell-based computational model

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

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

Microheterogeneity-induced conduction slowing and wavefront collisions govern macroscopic conduction behavior: A computational and experimental study

The incidence of cardiac arrhythmias is known to be associated with tissue heterogeneities including fibrosis. However, the impact of microscopic structural heterogeneities on conduction in excitable tissues remains poorly understood. In this study, we investigated how acellular microheterogeneities affect macroscopic conduction under conditions of normal and reduced excitability by utilizing a novel platform of paired in vitro and in silico studies to examine the mechanisms of conduction. Regular patterns of nonconductive micro-obstacles were created in confluent monolayers of the previously described engineered-excitable Ex293 cell line. Increasing the relative ratio of obstacle size to intra-obstacle strand width resulted in significant conduction slowing up to 23.6% and a significant increase in wavefront curvature anisotropy, a measure of spatial variation in wavefront shape. Changes in bulk electrical conductivity and in path tortuosity were insufficient to explain these observed macroscopic changes. Rather, microscale behaviors including local conduction slowing due to microscale branching, and conduction acceleration due to wavefront merging were shown to contribute to macroscopic phenomena. Conditions of reduced excitability led to further conduction slowing and a reversal of wavefront curvature anisotropy due to spatially non-uniform effects on microscopic slowing and acceleration. This unique experimental and computation platform provided critical mechanistic insights in the impact of microscopic heterogeneities on macroscopic conduction, pertinent to settings of fibrotic heart disease.

It is well known that perturbations in the heart structure are associated with the initiation and maintenance of clinically significant cardiac arrhythmia. While previous studies have examined how single structural perturbations affect local electrical conduction, our understanding of how numerous microscopic heterogeneities act in aggregate to alter macroscopic electrical behavior is limited. In this study, we utilized simplified engineered excitable cells that contain the minimal machinery of excitability and can be directly computationally modeled. By pairing experimental and computational studies, we showed that the microscopic branching and collisions of electrical waves slow and speed conduction, respectively, resulting in macroscopic changes in the speed and pattern of electrical activation. These microscale behaviors are significantly altered under reduced excitability, resulting in exaggerated collision effects. Overall, this study helps improve our understanding of how microscopic structural heterogeneities in excitable tissue lead to abnormal action potential propagation, conducive to arrhythmias.

Virtual cardiac monolayers for electrical wave propagation

The complex structure of cardiac tissue is considered to be one of the main determinants of an arrhythmogenic substrate. This study is aimed at developing the first mathematical model to describe the formation of cardiac tissue, using a joint in silico-in vitro approach. First, we performed experiments under various conditions to carefully characterise the morphology of cardiac tissue in a culture of neonatal rat ventricular cells. We considered two cell types, namely, cardiomyocytes and fibroblasts. Next, we proposed a mathematical model, based on the Glazier-Graner-Hogeweg model, which is widely used in tissue growth studies. The resultant tissue morphology was coupled to the detailed electrophysiological Korhonen-Majumder model for neonatal rat ventricular cardiomyocytes, in order to study wave propagation. The simulated waves had the same anisotropy ratio and wavefront complexity as those in the experiment. Thus, we conclude that our approach allows us to reproduce the morphological and physiological properties of cardiac tissue.

Cell-to-cell modeling of the interface between atrial and sinoatrial anisotropic heterogeneous nets

The transition between sinoatrial cells and atrial cells in the heart is not fully understood. Here we focus on cell-to-cell mathematical models involving typical sinoatrial cells and atrial cells connected with experimentally observed conductance values. We are interested mainly in the geometry of the microstructure of the conduction paths in the sinoatrial node. We show with some models that appropriate source-sink relationships between atrial and sinoatrial cells may occur according to certain geometric arrangements.

All-optical control of cardiac excitation: combined high-resolution optogenetic actuation and optical mapping

Cardiac tissue is an excitable system that can support complex spatiotemporal dynamics, including instabilities (arrhythmias) with lethal consequences. While over the last two decades optical mapping of excitation (voltage and calcium dynamics) has facilitated the detailed characterization of such arrhythmia events, until recently, no precise tools existed to actively interrogate cardiac dynamics in space and time. In this work, we discuss the combined use of new methods for space- and time-resolved optogenetic actuation and simultaneous fast, high resolution optical imaging of cardiac excitation waves. First, the mechanisms, limitations and unique features of optically induced responses in cardiomyocytes are outlined. These include the ability to bidirectionally control the membrane potential using depolarizing and hyperpolarizing opsins; the ability to induce prolonged sustained voltage changes; and the ability to control refractoriness and the shape of the cardiac action potential. At the syncytial tissue level, we discuss optogenetically enabled experimentation on cell-cell coupling, alteration of conduction properties and termination of propagating waves by light. Specific attention is given to space- and time-resolved application of optical stimulation using dynamic light patterns to perturb ongoing activation and to probe electrophysiological properties at desired tissue locations. The combined use of optical methods to perturb and to observe the system can offer new tools for precise feedback control of cardiac electrical activity, not available previously with pharmacological and electrical stimulation. These new experimental tools for all-optical electrophysiology allow for a level of precise manipulation and quantification of cardiac dynamics comparable in robustness to the computational setting, and can provide new insights into pacemaking, arrhythmogenesis and suppression or cardioversion.

Simulation of Cardiac Arrhythmias Using a 2D Heterogeneous Whole Heart Model

Simulation studies of cardiac arrhythmias at the whole heart level with electrocardiogram (ECG) gives an understanding of how the underlying cell and tissue level changes manifest as rhythm disturbances in the ECG. We present a 2D whole heart model (WHM2D) which can accommodate variations at the cellular level and can generate the ECG waveform. It is shown that, by varying cellular-level parameters like the gap junction conductance (GJC), excitability, action potential duration (APD) and frequency of oscillations of the auto-rhythmic cell in WHM2D a large variety of cardiac arrhythmias can be generated including sinus tachycardia, sinus bradycardia, sinus arrhythmia, sinus pause, junctional rhythm, Wolf Parkinson White syndrome and all types of AV conduction blocks. WHM2D includes key components of the electrical conduction system of the heart like the SA (Sino atrial) node cells, fast conducting intranodal pathways, slow conducting atriovenctricular (AV) node, bundle of His cells, Purkinje network, atrial, and ventricular myocardial cells. SA nodal cells, AV nodal cells, bundle of His cells, and Purkinje cells are represented by the Fitzhugh-Nagumo (FN) model which is a reduced model of the Hodgkin-Huxley neuron model. The atrial and ventricular myocardial cells are modeled by the Aliev-Panfilov (AP) two-variable model proposed for cardiac excitation. WHM2D can prove to be a valuable clinical tool for understanding cardiac arrhythmias.

Simulation of Ectopic Pacemakers in the Heart: Multiple Ectopic Beats Generated by Reentry inside Fibrotic Regions

The inclusion of nonconducting media, mimicking cardiac fibrosis, in two models of cardiac tissue produces the formation of ectopic beats. The fraction of nonconducting media in comparison with the fraction of healthy myocytes and the topological distribution of cells determines the probability of ectopic beat generation. First, a detailed subcellular microscopic model that accounts for the microstructure of the cardiac tissue is constructed and employed for the numerical simulation of action potential propagation. Next, an equivalent discrete model is implemented, which permits a faster integration of the equations. This discrete model is a simplified version of the microscopic model that maintains the distribution of connections between cells. Both models produce similar results when describing action potential propagation in homogeneous tissue; however, they slightly differ in the generation of ectopic beats in heterogeneous tissue. Nevertheless, both models present the generation of reentry inside fibrotic tissues. This kind of reentry restricted to microfibrosis regions can result in the formation of ectopic pacemakers, that is, regions that will generate a series of ectopic stimulus at a fast pacing rate. In turn, such activity has been related to trigger fibrillation in the atria and in the ventricles in clinical and animal studies.

Role of connexins and pannexins in cardiovascular physiology

Connexins and pannexins form connexons, pannexons and membrane channels, which are critically involved in many aspects of cardiovascular physiology. For that reason, a vast number of studies have addressed the role of connexins and pannexins in the arterial and venous systems as well as in the heart. Moreover, a role for connexins in lymphatics has recently also been suggested. This review provides an overview of the current knowledge regarding the involvement of connexins and pannexins in cardiovascular physiology.

A brief history of tissue models for cardiac electrophysiology

The last four decades have produced a number of significant advances in the developments of computer models to simulate and investigate the electrical activity of cardiac tissue. The tissue descriptions that underlie these simulations have been built from a combination of clever insight and careful comparison with measured data at multiple scales. Tissue models have not only led to greater insights into the mechanisms of life-threatening arrhythmias but have been used to engineer new therapies to treat the consequences of cardiac disease. This paper is a look back at the early years in the cardiac modeling and the challenges facing the field as models move toward the clinic.

Dynamics of propagation of premature impulses in structurally remodeled infarcted myocardium: a computational analysis

Initiation of cardiac arrhythmias typically follows one or more premature impulses either occurring spontaneously or applied externally. In this study, we characterize the dynamics of propagation of single (S2) and double premature impulses (S3), and the mechanisms of block of premature impulses at structural heterogeneities caused by remodeling of gap junctional conductance (Gj) in infarcted myocardium. Using a sub-cellular computer model of infarcted tissue, we found that |I_Na,max|, prematurity (coupling interval with the previous impulse), and conduction velocity (CV) of premature impulses change dynamically as they propagate away from the site of initiation. There are fundamental differences between the dynamics of propagation of S2 and S3 premature impulses: for S2 impulses |I_Na,max| recovers fast, prematurity decreases and CV increases as propagation proceeds; for S3 impulses low values of |I_Na,max| persist, prematurity could increase, and CV could decrease as impulses propagate away from the site of initiation. As a consequence it is more likely that S3 impulses block at sites of structural heterogeneities causing source/sink mismatch than S2 impulses block. Whether premature impulses block at Gj heterogeneities or not is also determined by the values of Gj (and the space constant λ) in the regions proximal and distal to the heterogeneity: when λ in the direction of propagation increases >40%, premature impulses could block. The maximum slope of CV restitution curves for S2 impulses is larger than for S3 impulses. In conclusion: (1) The dynamics of propagation of premature impulses make more likely that S3 impulses block at sites of structural heterogeneities than S2 impulses block; (2) Structural heterogeneities causing an increase in λ (or CV) of >40% could result in block of premature impulses; (3) A decrease in the maximum slope of CV restitution curves of propagating premature impulses is indicative of an increased potential for block at structural heterogeneities.

Role of the intercalated disc in cardiac propagation and arrhythmogenesis

This review article discusses mechanisms underlying impulse propagation in cardiac muscle with specific emphasis on the role of the cardiac cell-to-cell junction, called the “intercalated disc.”The first part of this review deals with the role of gap junction channels, formed by connexin proteins, as a determinant of impulse propagation. It is shown that, depending on the underlying structure of the cellular network, decreasing the conductance of gap junction channels (so-called “electrical uncoupling”) may either only slow, or additionally stabilize propagation and reverse unidirectional propagation block to bidirectional propagation. This is because the safety factor for propagation increases with decreasing intercellular electrical conductance. The role of heterogeneous connexin expression, which may be present in disease states, is also discussed. The hypothesis that so-called ephaptic impulse transmission plays a role in heart and can substitute for electrical coupling has been revived recently. Whereas ephaptic transmission can be demonstrated in theoretical simulations, direct experimental evidence has not yet been presented. The second part of this review deals with the interaction of three protein complexes at the intercalated disc: (1) desmosomal and adherens junction proteins, (2) ion channel proteins, and (3) gap junction channels consisting of connexins. Recent work has revealed multiple interactions between these three protein complexes which occur, at least in part, at the level of protein trafficking. Such interactions are likely to play an important role in the pathogenesis of arrhythmogenic cardiomyopathy, and may reveal new therapeutic concepts and targets.

Nonlinear behaviour of conduction and block in cardiac tissue with heterogeneous expression of connexin 43

Altered gap junctional coupling potentiates slow conduction and arrhythmias. To better understand how heterogeneous connexin expression affects conduction at the cellular scale, we investigated conduction in tissue consisting of two cardiomyocyte populations expressing different connexin levels. Conduction was mapped using microelectrode arrays in cultured strands of foetal murine ventricular myocytes with predefined contents of connexin 43 knockout (Cx43KO) cells. Corresponding computer simulations were run in randomly generated two-dimensional tissues mimicking the cellular architecture of the strands. In the cultures, the relationship between conduction velocity (CV) and Cx43KO cell content was nonlinear. CV first decreased significantly when Cx43KO content was increased from 0 to 50%. When the Cx43KO content was ≥60%, CV became comparable to that in 100% Cx43KO strands. Co-culturing Cx43KO and wild-type cells also resulted in significantly more heterogeneous conduction patterns and in frequent conduction blocks. The simulations replicated this behaviour of conduction. For Cx43KO contents of 10-50%, conduction was slowed due to wavefront meandering between Cx43KO cells. For Cx43KO contents ≥60%, clusters of remaining wild-type cells acted as electrical loads that impaired conduction. For Cx43KO contents of 40-60%, conduction exhibited fractal characteristics, was prone to block, and was more sensitive to changes in ion currents compared to homogeneous tissue. In conclusion, conduction velocity and stability behave in a nonlinear manner when cardiomyocytes expressing different connexin amounts are combined. This behaviour results from heterogeneous current-to-load relationships at the cellular level. Such behaviour is likely to be arrhythmogenic in various clinical contexts in which gap junctional coupling is heterogeneous.

Regulation of gap junctions by nitric oxide influences the generation of arrhythmias resulting from acute ischemia and reperfusion in vivo

Myocardial ischemia resulting from sudden occlusion of a coronary artery is one of the major causes in the appearance of severe, often life-threatening ventricular arrhythmias. Although the underlying mechanisms of these acute arrhythmias are many and varied, there is no doubt that uncoupling of gap junctions (GJs) play an important role especially in arrhythmias that are generated during phase Ib, and often terminate in sudden cardiac death. In the past decades considerable efforts have been made to explore mechanisms which regulate the function of GJs, and to find new approaches for protection against arrhythmias through the modulation of GJs. These investigations led to the development of GJ openers and inhibitors. The pharmacological modulation of GJs, however, resulted in conflicting results. It is still not clear whether opening or closing of GJs would be advantageous for the ischemic myocardium. Both maneuvers can result in protection, depending on the models, endpoints and the time of opening and closing of GJs. Furthermore, although there is substantial evidence that preconditioning decreases or delays the uncoupling of GJs, the precise mechanisms by which this attains have not yet been elucidated. In our own studies in anesthetized dogs preconditioning suppressed the ischemia and reperfusion-induced ventricular arrhythmias, and this protection was associated with the preservation of GJ function, manifested in less marked changes in electrical impedance, as well as in the maintenance of GJ permeability and phosphorylation of connexin43. Since we have substantial previous evidence that nitric oxide (NO) is an important trigger and mediator of the preconditioning-induced antiarrhythmic protection, we hypothesized that NO, among its several effects, may lead to this protection by influencing cardiac GJs. The hypotheses and theories relating to the pharmacological modulation of GJs will be discussed with particular attention to the role of NO.

Microscopic variations in interstitial and intracellular structure modulate the distribution of conduction delays and block in cardiac tissue with source-load mismatch

AIMS

Reentrant activity in the heart is often correlated with heterogeneity in both the intracellular structure and the interstitial structure surrounding cells; however, the combined effect of cardiac microstructure and interstitial resistivity in regions of source-load mismatch is largely unknown. The aim of this study was to investigate how microstructural variations in cell arrangement and increased interstitial resistivity influence the spatial distribution of conduction delays and block in poorly coupled regions of tissue.

METHODS AND RESULTS

Two-dimensional 0.6 cm × 0.6 cm computer models with idealized and realistic cellular structure were used to represent a monolayer of ventricular myocytes. Gap junction connections were distributed around the periphery of each cell at 10 μm intervals. Regions of source-load mismatch were added to the models by increasing the gap junction and interstitial resistivity in one-half of the tissue. Heterogeneity in cell shape and cell arrangement along the boundary between well-coupled and poorly coupled tissue increased variability in longitudinal conduction delays to as much as 10 ms before the onset of conduction block, resulting in wavefront breakthroughs with pronounced curvature at distinct points along the boundary. Increasing the effective interstitial resistivity reduced source-load mismatch at the transition boundary, which caused a decrease in longitudinal conduction delay and an increase in the number of wavefront breakthroughs.

CONCLUSION

Microstructural variations in cardiac tissue facilitate the formation of isolated sites of wavefront breakthrough that may enable abnormal electrical activity in small regions of diseased tissue to develop into more widespread reentrant activity.

Association of inflammation with atrial fibrillation in hyperthyroidism

Objectives

The aim of this study was to evaluate the relationship between inflammation and development of atrial fibrillation (AF) in patients with hyperthyroidism.

Methods

A total of 65 patients with newly diagnosed hyperthyroidism, 35 of whom were in sinus rhythm and 30 of whom in AF. Thirty five age- and gender-matched patients in a control group were included in the study. Factors associated with the development of AF were evaluated by multivariate regression analysis.

Results

Factors associated with AF in multivariate analysis included high sensitivity C reactive protein (HsCRP) [odds ratio (OR): 11.19; 95% confidence interval (95% CI): 1.80-69.53; P = 0.003], free T4 (OR: 8.76; 95% CI: 2.09–36.7; P = 0.003), and left atrial diameter (OR: 1.25; 95% CI: 1.06–1.47; P = 0.008).

Conclusions

The results of the present study suggest that HsCRP, an indicator of inflammation, free T4 and left atrial diameter are associated with the development AF in patients with hyperthyroidism.

Simulations of complex and microscopic models of cardiac electrophysiology powered by multi-GPU platforms

Key aspects of cardiac electrophysiology, such as slow conduction, conduction block, and saltatory effects have been the research topic of many studies since they are strongly related to cardiac arrhythmia, reentry, fibrillation, or defibrillation. However, to reproduce these phenomena the numerical models need to use subcellular discretization for the solution of the PDEs and nonuniform, heterogeneous tissue electric conductivity. Due to the high computational costs of simulations that reproduce the fine microstructure of cardiac tissue, previous studies have considered tissue experiments of small or moderate sizes and used simple cardiac cell models. In this paper, we develop a cardiac electrophysiology model that captures the microstructure of cardiac tissue by using a very fine spatial discretization (8 μm) and uses a very modern and complex cell model based on Markov chains for the characterization of ion channel's structure and dynamics. To cope with the computational challenges, the model was parallelized using a hybrid approach: cluster computing and GPGPUs (general-purpose computing on graphics processing units). Our parallel implementation of this model using a multi-GPU platform was able to reduce the execution times of the simulations from more than 6 days (on a single processor) to 21 minutes (on a small 8-node cluster equipped with 16 GPUs, i.e., 2 GPUs per node).

Effect of heterogeneities in the cellular microstructure on propagation of the cardiac action potential

Cardiac arrhythmias are initiated in regions that undergo cellular remodeling as a result of disease. Using a sub-cellular model of myocardium, we studied the mechanism of block caused by tissue microstructure remodeling: cell geometry [quantified as length/width (L/W) cell ratio] and cell-to-cell coupling (G(j)). Heterogeneities in cell L/W ratio and G ( j ) lead to block when excitability is reduced and the corresponding space constant λ (in the direction of propagation) increases by >40 %. Tissue architectures with elongated cells (i.e. large cell L/W ratios) that are better coupled (i.e. large G(j)) are less prone to block at sites of regional heterogeneities in cell geometry and/or cell coupling than tissue architectures consisting of cells with smaller L/W ratios and/or poorer coupling. Whether an increase in tissue anisotropic ratio (ANR) is arrhythmogenic or not depends on the cellular mechanism of the increase: ANR leads to an increased risk of block when G(j) decreases, but to a decreased risk of block when cell L/W ratio increases. Our findings are useful to understand the mechanisms of block in cardiac pathologies that result in tissue architecture remodeling.

The involvement of gap junctions in the delayed phase of the protection induced by cardiac pacing in dogs

The present study has examined the role of GJ (gap junctions) in the delayed anti-arrhythmic effect of cardiac pacing, with particular reference to the time-course changes in Cx43 (connexin43) expression both after pacing (4×5 min, at a rate of 240 beats/min) and 24 h later, when the dogs were subjected to a 25 min occlusion and reperfusion of the LAD (left anterior descending coronary artery). Compared with the SP (sham-paced) controls (n=20), in dogs paced 24 h previously (n=16) there were reductions in arrhythmia severity [e.g. number of VPB (ventricular premature beats) during occlusion 294±78 compared with 63±25; survival from the combined ischaemia/reperfusion insult 20% compared with 78%], and in other ischaemic changes [epicardial ST-segment, TAT (total activation time) and tissue impedance]. Pacing also prevented the ischaemia-induced structural impairment of the intercalated discs, and preserved GJ permeability and Cx43 phosphorylation, without modifying Cx43 protein content. Following cardiac pacing the membrane and total Cx43 protein contents were unchanged up to 6 h, but were significantly reduced 12 h later (preceded by a down-regulation of Cx43 mRNA at 6 h), and returned to normal by 24 h. Interestingly, dogs that were subjected to ischaemia 12 h after cardiac pacing showed increased arrhythmia generation. We conclude that cardiac pacing results in time-dependent changes in Cx43 expression, which may alter GJ function and influence arrhythmia generation during a subsequent ischaemia/reperfusion insult. This effect is manifested in protection 24 h after pacing, but of potential clinical interest is the finding that there is a time interval after pacing during which an ischaemic event may generate severe ventricular arrhythmias.

Time frequency power profile of QRS complex obtained with wavelet transform in spontaneously hypertensive rats

We evaluated whether frequency analysis could detect the development of interstitial fibrosis in rats. SHR/Izm and age-matched WKY/Izm were used. Limb lead II electrocardiograms were recorded. Continuous wavelet transform (CWT) was applied for the time-frequency analysis. The integrated time-frequency power (ITFP) between QRS complexes was measured and compared between groups. The ITFP at low-frequency bands (≤125Hz) was significantly higher in SHR/Izm. The percent change of ITFP showed the different patterns between groups. Prominent interstitial fibrosis with an increase in TIMP-1 mRNA expression was also observed in SHR/Izm. These results were partly reproduced in a computer simulation.

ROLE OF CELLULAR COUPLING AND DISPERSION OF REFRACTORINESS IN CARDIAC ARRHYTHMIAS: A SIMULATION STUDY

Computer simulation is applied to study the role of cellular coupling, dispersion of refractoriness as well as both of them, in the mechanisms underlying cardiac arrhythmias. We first assumed that local ischemia mainly induces cell to cell dispersion in the coupling resistance (case 1), refractory period (case 2) or both (case 3). Numerical experiments, based on the van Capelle and Durrer model, showed that vortices could not be induced in these conditions. In order to be more realistic about coronary circulation we simulated a patchy dispersion of cellular properties, each patch corresponding to the zone irrigated by a small coronary artery. In these conditions, a single activation wave could give rise to abnormal activities. Probabilities of reentry, estimated for the three cases cited above, showed that a severe alteration of the coupling resistance may be an important factor in the genesis of reentry. Moreover, use of isochronal maps revealed that vortices were both stable and sustained with an alteration of coupling alone or along with reductions of action potential duration. Conversely, simulations with reduction of the refractoriness alone induced only transient patterns.

Myofibroblasts cause heterogeneous Cx43 reduction and are unlikely to be coupled to myocytes in the healing canine infarct

Following myocardial infarction (MI) inflammatory responses transform cardiac fibroblasts to myofibroblasts, which in vitro studies show form heterocellular gap junctions with cardiac myocytes via Connexin43 (Cx43). The ability to form heterocellular junctions in the intact heart and the impact of these junctions on propagation is unclear. We used a canine model of MI and characterized the distribution and quantity of myofibroblasts in surviving epicardial cells [epicardial border zone (EBZ)]. We found a significant increase in myofibroblasts within the EBZ and no gap junction plaques between myofibroblasts and myocytes. Because myofibroblasts produce IL-1β, which downregulates Cx43, we asked whether myofibroblast proliferation causes loss of Cx43 near myofibroblast clusters. In vitro studies showed that IL-1β caused loss of Cx43 and reduced coupling. Western blot showed a significant increase of IL-1β in the EBZ, and immunohistochemistry showed a loss of Cx43 in regions of myofibroblasts in the intact heart. Additionally, dye studies in intact heart showed no coupling between myocytes and myofibroblasts. To quantify the effect of myofibroblasts on propagation we used a two-dimensional subcellular computer model of the EBZ, which showed that heterogeneities in myofibroblast density lead to conduction abnormalities. In conclusion, an increase of myofibroblasts in the infarcted heart causes heterogeneous Cx43 levels, possibly as a result of the release of IL-1β and decreased cell-cell communication, which leads to conduction abnormalities following MI.

Conduction block in micropatterned cardiomyocyte cultures replicating the structure of ventricular cross-sections

AIMS

Structural and functional heterogeneities in cardiac tissue have been implicated in conduction block and arrhythmogenesis. However, the propensity of specific sites within the heart to initiate conduction block has not been systematically explored. We utilized cardiomyocyte cultures replicating the realistic, magnetic resonance imaging-measured tissue boundaries and fibre directions of ventricular cross-sections to investigate their roles in the development of conduction block.

METHODS AND RESULTS

The Sprague-Dawley neonatal rat cardiomyocytes were micropatterned to obtain cultures with realistic ventricular tissue boundaries and either random or realistic fibre directions. Rapid pacing was applied at multiple sites, with action potential propagation optically mapped. Excitation either failed at the stimulus site or conduction block developed remotely, often initiating reentry. The incidence of conduction block in isotropic monolayers (0% of cultures) increased with the inclusion of realistic tissue boundaries (17%) and further with realistic fibre directions (34%). Conduction block incidence was stimulus site-dependent and highest (77%) with rapid pacing from the right ventricular (RV) free wall. Furthermore, conduction block occurred exclusively at the insertion of the RV free wall into the septum, where structure-mediated current source-load mismatches acutely reduced wavefront and waveback velocity. Tissue boundaries and sharp gradients in fibre direction uniquely determined the evolution, shape, and position of conduction block lines.

CONCLUSION

Our study suggests that specific micro- and macrostructural features of the ventricle determine the incidence and spatiotemporal characteristics of conduction block, independent of spatial heterogeneities in ion channel expression.

Patterned cardiomyocytes on microelectrode arrays as a functional, high information content drug screening platform

Cardiac side effects are one of the major causes of drug candidate failures in preclinical drug development or in clinical trials and are responsible for the retraction of several already marketed therapeutics. Thus, the development of a relatively high-throughput, high information content tool to screen drugs and toxins would be important in the field of cardiac research and drug development. In this study, recordings from commercial multielectrode arrays were combined with surface patterning of cardiac myocyte monolayers to enhance the information content of the method; specifically, to enable the measurement of conduction velocity, refractory period after action potentials and to create a functional re-entry model. Two drugs, 1-Heptanol, a gap junction blocker, and Sparfloxacin, a fluoroquinone antibiotic, were tested in this system. 1-Heptanol administration resulted in a marked reduction in conduction velocity, whereas Sparfloxacin caused rapid, irregular and unsynchronized activity, indicating fibrillation. As shown in these experiments, patterning of cardiac myocyte monolayers solved several inherent problems of multielectrode recordings, increased the temporal resolution of conduction velocity measurements, and made the synchronization of external stimulation with action potential propagation possible for refractory period measurements. This method could be further developed as a cardiac side effect screening platform after combination with human cardiomyocytes.

Pathophysiological mechanisms of atrial fibrillation: a translational appraisal

Atrial fibrillation (AF) is an arrhythmia that can occur as the result of numerous different pathophysiological processes in the atria. Some aspects of the morphological and electrophysiological alterations promoting AF have been studied extensively in animal models. Atrial tachycardia or AF itself shortens atrial refractoriness and causes loss of atrial contractility. Aging, neurohumoral activation, and chronic atrial stretch due to structural heart disease activate a variety of signaling pathways leading to histological changes in the atria including myocyte hypertrophy, fibroblast proliferation, and complex alterations of the extracellular matrix including tissue fibrosis. These changes in electrical, contractile, and structural properties of the atria have been called "atrial remodeling." The resulting electrophysiological substrate is characterized by shortening of atrial refractoriness and reentrant wavelength or by local conduction heterogeneities caused by disruption of electrical interconnections between muscle bundles. Under these conditions, ectopic activity originating from the pulmonary veins or other sites is more likely to occur and to trigger longer episodes of AF. Many of these alterations also occur in patients with or at risk for AF, although the direct demonstration of these mechanisms is sometimes challenging. The diversity of etiological factors and electrophysiological mechanisms promoting AF in humans hampers the development of more effective therapy of AF. This review aims to give a translational overview on the biological basis of atrial remodeling and the proarrhythmic mechanisms involved in the fibrillation process. We pay attention to translation of pathophysiological insights gained from in vitro experiments and animal models to patients. Also, suggestions for future research objectives and therapeutical implications are discussed.

Hierarchical architecture influences calcium dynamics in engineered cardiac muscle

Changes in myocyte cell shape and tissue structure are concurrent with changes in electromechanical function in both the developing and diseased heart. While the anisotropic architecture of cardiac tissue is known to influence the propagation of the action potential, the influence of tissue architecture and its potential role in regulating excitation-contraction coupling (ECC) are less well defined. We hypothesized that changes in the shape and the orientation of cardiac myocytes induced by spatial arrangement of the extracellular matrix (ECM) affects ECC. To test this hypothesis, we isolated and cultured neonatal rat ventricular cardiac myocytes on various micropatterns of fibronectin where they self-organized into tissues with varying degrees of anisotropy. We then measured the morphological features of these engineered myocardial tissues across several hierarchical dimensions by measuring cellular aspect ratio, myocyte area, nuclear density and the degree of cytoskeletal F-actin alignment. We found that when compared with isotropic tissues, anisotropic tissues have increased cellular aspect ratios, increased nuclear densities, decreased myocyte cell areas and smaller variances in actin alignment. To understand how tissue architecture influences cardiac function, we studied the role of anisotropy on intracellular calcium ([Ca(2+)](i)) dynamics by characterizing the [Ca(2+)](i)-frequency relationship of electrically paced tissues. When compared with isotropic tissues, anisotropic tissues displayed significant differences in [Ca(2+)](i) transients, decreased diastolic baseline [Ca(2+)](i) levels and greater [Ca(2+)](i) influx per cardiac cycle. These results suggest that ECM cues influence tissue structure at cellular and subcellular levels and regulate ECC.

A computer model of engineered cardiac monolayers

Engineered monolayers created using microabrasion and micropatterning methods have provided a simplified in vitro system to study the effects of anisotropy and fiber direction on electrical propagation. Interpreting the behavior in these culture systems has often been performed using classical computer models with continuous properties. However, such models do not account for the effects of random cell shapes, cell orientations, and cleft spaces inherent in these monolayers on the resulting wavefront conduction. This work presents a novel methodology for modeling a monolayer of cardiac tissue in which the factors governing cell shape, cell-to-cell coupling, and degree of cleft space are not constant but rather are treated as spatially random with assigned distributions. This modeling approach makes it possible to simulate wavefront propagation in a manner analogous to performing experiments on engineered monolayer tissues. Simulated results are compared to previously published measured data from monolayers used to investigate the role of cellular architecture on conduction velocities and anisotropy ratios. We also present an estimate for obtaining the electrical properties from these networks and demonstrate how variations in the discrete cellular architecture affect the macroscopic conductivities. The simulations support the common assumption that under normal ranges of coupling strength, tissues with relatively uniform distributions of cell shapes and connectivity can be represented using continuous models with conductivities derived from random discrete cellular architecture using either global or local estimates. The results also reveal that in the presence of abrupt changes in cell orientation, local estimates of tissue properties predict smoother changes in conductivity that may not adequately predict the discrete nature of propagation at the transition sites.

Effect of cell geometry on conduction velocity in a subcellular model of myocardium

We have studied the effect of cell geometry on propagation velocity of the cardiac impulse using a subcellular computer model of myocardium. Variation of cell size has only small effects on longitudinal and transverse conduction velocities, when the ratio of cell length/width is constant, for cell sizes (length x width) between (60 microm x 20 microm) and (120 microm x 40 microm). The results were not dependent on gap-junction conductance (range 0.25-1 microS), gap-junction distribution, or the specific tissue architecture. Longitudinal conduction velocity increased with the cell length/width ratio and transverse velocity decreased. The cell length/width ratio was a good estimator of the anisotropic ratio. In conclusion, cell length/width ratio is more important than cell size in determining conduction velocity.

Extracellular space attenuates the effect of gap junctional remodeling on wave propagation: a computational study

Ionic channels and gap junctions are remodeled in cells from the 5-day epicardial border zone (EBZ) of the healing canine infarct. The main objective of the study was to determine the effect of gap junctional conductance (Gj) remodeling and Cx43 redistribution to the lateral membrane on conduction velocity (theta) and anisotropic ratio, and how gap junctional remodeling is modulated by the extracellular space. We first implemented subcellular monodomain and two-domain computer models of normal epicardium (NZ) to understand how extracellular space modulates the relationship between Gj and theta in NZ. We found that the extracellular space flattens the Gj-theta relationship, thus theta becomes less sensitive to changes in Gj. We then investigated the functional consequences of Gj remodeling and Cx43 distribution in subcellular computer models of cells of the outer pathway (IZo) and central pathway (IZc) of reentrant circuits. In IZo cells, side-to-side (transverse) Gj is 10% the value in NZ cells. Such Gj remodeling causes a 45% decrease in transverse theta (theta(T)). Inclusion of an extracellular space reduces the decrease in theta(T) to 31%. In IZc cells, Cx43 redistribution along the lateral membrane results in a 29% increase in theta(T). That increase in theta(T) is a consequence of the decrease in access resistance to the Cx43 plaques that occur with the Cx43 redistribution. Extracellular space reduces the increase in theta(T) to 10%.

IN CONCLUSION

1), The extracellular space included in normal epicardial simulations flattens the Gj-theta relationship with theta becoming less sensitive to changes in Gj. 2), The extracellular space attenuates the effects of gap junction epicardial border zone remodeling (i.e., Gj reduction and Cx43 lateralization) on theta(T).

Uniform action potential repolarization within the sarcolemma of in situ ventricular cardiomyocytes

Previous studies have speculated, based on indirect evidence, that the action potential at the transverse (t)-tubules is longer than at the surface membrane in mammalian ventricular cardiomyocytes. To date, no technique has enabled recording of electrical activity selectively at the t-tubules to directly examine this hypothesis. We used confocal line-scan imaging in conjunction with the fast response voltage-sensitive dyes ANNINE-6 and ANNINE-6plus to resolve action potential-related changes in fractional dye fluorescence (DeltaF/F) at the t-tubule and surface membranes of in situ mouse ventricular cardiomyocytes. Peak DeltaF/F during action potential phase 0 depolarization averaged -21% for both dyes. The shape and time course of optical action potentials measured with the water-soluble ANNINE-6plus were indistinguishable from those of action potentials recorded with intracellular microelectrodes in the absence of the dye. In contrast, optical action potentials measured with the water-insoluble ANNINE-6 were significantly prolonged compared to the electrical recordings obtained from dye-free hearts, suggesting electrophysiological effects of ANNINE-6 and/or its solvents. With either dye, the kinetics of action potential-dependent changes in DeltaF/F during repolarization were found to be similar at the t-tubular and surface membranes. This study provides what to our knowledge are the first direct measurements of t-tubule electrical activity in ventricular cardiomyocytes, which support the concept that action potential duration is uniform throughout the sarcolemma of individual cells.

Loading effect of fibroblast-myocyte coupling on resting potential, impulse propagation, and repolarization: insights from a microstructure model

The numerous nonmyocytes present within the myocardium may establish electrical connections with myocytes through gap junctions, formed naturally or as a result of a cell therapy. The strength of the coupling and its potential impact on action potential characteristics and conduction are not well understood. This study used computer simulation to investigate the load-induced electrophysiological consequences of the coupling of myocytes with fibroblasts, where the fibroblast resting potential, density, distribution, and coupling strength were varied. Conduction velocity (CV), upstroke velocity, and action potential duration (APD) were analyzed for longitudinal and transverse impulse propagation in a two-dimensional microstructure tissue model, developed to represent a monolayer culture of cardiac cells covered by a layer of fibroblasts. The results show that 1) at weak coupling (<0.25 nS), the myocyte resting potential was elevated, leading to CV up to 5% faster than control; 2) at intermediate coupling, the myocyte resting potential elevation saturated, whereas the current flowing from the myocyte to the fibroblast progressively slowed down both CV and upstroke velocity; 3) at strong couplings (>8 nS), all of the effects saturated; and 4) APD at 90% repolarization was usually prolonged by 0-20 ms (up to 60-80 ms for high fibroblast density and coupling) by the coupling to fibroblasts. The changes in APD depended on the fibroblast resting potential. This complex, coupling-dependent interaction of fibroblast and myocytes also has relevance to the integration of other nonmyocytes in the heart, such as those used in cellular therapies.

Study of myocardial fiber pathway using magnetic resonance diffusion tensor imaging

The purpose of this study was to investigate myocardial fiber pathway distribution in order to provide supplemental information on myocardial fiber architecture and cardiac mechanics. Diffusion tensor imaging (DTI) with medium diffusion resolution (15 directions) was performed on normal canine heart samples (N=6) fixed in formalin. With the use of diffusion tensor fiber tracking, left ventricle (LV) myocardial fiber pathways and helix angles were computed pixel by pixel at short-axis slices from base to apex. Distribution of DTI-tracked fiber pathway length and number was analyzed quantitatively as a function of fiber helix angle in step of 9 degrees . The long fiber pathways were found to have small helix angles. They are mostly distributed in the middle myocardium and run circumferentially. Fiber pathways tracked at the middle and upper LV are generally longer than those near the apex. Majority of fiber pathways have small helix angles between -20 degrees and 20 degrees , dominating the fiber architecture in myocardium. Likely, such myocardial fiber pathway measurement by DTI may reflect the spatial connectiveness or connectivity of elastic myofiber bundles along their preferential pathway of electromechanical activation. The dominance of the long and circumferentially running fiber pathways found in the study may explain the circumferential predominance in left ventricular contraction.

Dynamics of sustained reentry in a loop model with discrete gap junction resistances

The dynamics of reentry is studied in a one-dimensional loop of model cardiac cells with discrete intercellular gap junction resistance (R). Each cell is represented by a continuous cable with ionic current given by a modified Beeler-Reuter formulation. For R below a limiting value, propagation is found to change from period-1 to quasiperiodic (QP) at a critical loop length (L(crit)) that decreases with R. Quasiperiodic reentry exists from L(crit) to a minimum length (L(min)), which also shortens with R. The decrease of L(crit) (R) is not a simple scaling, but the bifurcation can still be predicted from the slope of the restitution curve giving the duration of the action potential as a function of the diastolic interval. However, the shape of the restitution curve changes with R. An increase of R does not seem to increase the number of possible QP solutions since, as in the continuous cable, only two QP modes of propagation were found despite an extensive search through alternative initial conditions.

Oblique propagation of activation allows the detection of uncoupling microstructures from cardiac near field behavior

Wave fronts of cardiac activation, when propagating oblique to the fiber axis, reveal small fractionations and distortions in extracellular potential waveforms Phi(e) as well as in parameters derived from Phi(e) like dPhi(e)/dt and the gradient of Phi(e), the cardiac near field E. dPhi(e)/dt shows multiple deflections and E changes its morphology forming abnormal or even two or multiple subsequent loops. We analyze segments of such irregular loops of E obtained from in-vitro experiments and from computer simulation of a 2D-tissue sheet with a longitudinal oriented obstacle. In computer simulations we found that the individual sections of E reflect fairly well individual pathways of activation separated and delayed by the presence of a structural obstacle similar to connective tissue embedded longitudinally in heart tissue. Electrophysiological experiments with papillary muscles confirm this near field behavior.

The helical ventricular myocardial band of Torrent-Guasp

We live in an era of substantial progress in understanding myocardial structure and function at genetic, molecular, and microscopic levels. Yet, ventricular myocardium has proven remarkably resistant to macroscopic analyses of functional anatomy. Pronounced and practically indefinite global and local structural anisotropy of its fibers and other ventricular wall constituents produces electrical and mechanical properties that are nonlinear, anisotropic, time varying, and spatially inhomogeneous. The helical ventricular myocardial band of Torrent-Guasp is a revolutionary new concept in understanding global, 3-dimensional, functional architecture of the ventricular myocardium. This concept defines the principal, cumulative vectors, integrating the tissue architecture (ie, form) and net forces developed (ie, function) within the ventricular mass. The primary purpose of this review is to emphasize the importance of this concept, in the light of collaborative efforts to establish an integrative approach, defining ventricular form and function by linking across multiple scales of biological organization, as explained in the ongoing Physiome project. Because one of the most important scientific missions in this century is integration of basic research with clinical medicine, we believe that this knowledge is not of merely academic importance, but is also the essential prerequisite in clinical evaluation and treatment of different heart diseases.

Mechanism of origin of conduction disturbances in aging human atrial bundles: experimental and model study

BACKGROUND

Aging is associated with a significant increase in atrial tachyarrhythmias, especially atrial fibrillation. A macroscopic repolarization gradient created artificially by a stimulus at one site before a premature stimulus from a second site is widely considered to be part of the experimental protocol necessary for the initiation of such arrhythmias in the laboratory. How such gradients occur naturally in aging atrial tissue is unknown.

OBJECTIVE

The objective of this study was to determine if the pattern of cellular connectivity in aging human atrial bundles produces a mechanism for variable early premature responses.

METHODS

Extracellular and intracellular potentials were recorded after control and premature stimuli at a single site in aging human atrial bundles. We also measured cellular geometry, the distribution of connexins, and the distribution of collagenous septa. A model of the atrial bundles was constructed based on the morphological results. Action potential propagation and the sodium current were analyzed after premature stimuli in the model.

RESULTS

Similar extracellular potential waveform responses occurred after early premature stimuli in the aging bundles and in the model. Variable premature conduction patterns were accounted for by the single model of aging atrial structure. A major feature of the model results was that the conduction events and the magnitude of the sodium current at multiple sites were very sensitive to small changes in the location and the timing of premature stimuli.

CONCLUSION

In aging human atrial bundles stimulated from only a single site, premature stimuli induce variable arrhythmogenic conduction responses. The generation of these responses is greatly enhanced by remodeling of cellular connectivity during aging. The results provide insight into sodium current structural interactions as a general mechanism of arrhythmogenic atrial responses to premature stimuli.