Excitation-contraction coupling between human atrial myocytes with fibroblasts and stretch activated channel current: a simulation study

A detailed mathematical model of the human atrial cardiomyocyte: integration of electrophysiology and cardiomechanics

Mechano-electric regulations (MER) play an important role in the maintenance of cardiac performance. Mechano-calcium and mechano-electric feedback (MCF and MEF) pathways adjust the cardiomyocyte contractile force according to mechanical perturbations and affects electro-mechanical coupling. MER integrates all these regulations in one unit resulting in a complex phenomenon. Computational modelling is a useful tool to accelerate the mechanistic understanding of complex experimental phenomena. We have developed a novel model that integrates the MER loop for human atrial cardiomyocytes with proper consideration of feedforward and feedback pathways. The model couples a modified version of the action potential (AP) Koivumäki model with the contraction model by Quarteroni group. The model simulates iso-sarcometric and isometric twitches and the feedback effects on AP and Ca2+ -handling. The model showed a biphasic response of Ca2+ transient (CaT) peak to increasing pacing rates and highlights the possible mechanisms involved. The model has shown a shift of the threshold for AP and CaT alternans from 4.6 to 4 Hz under post-operative atrial fibrillation, induced by depressed SERCA activity. The alternans incidence was dependent on a chain of mechanisms including RyRs availability time, MCF coupling, CaMKII phosphorylation, and the stretch levels. As a result, the model predicted a 10% slowdown of conduction velocity for a 20% stretch, suggesting a role of stretch in creation of substrate formation for atrial fibrillation. Overall, we conclude that the developed model provides a physiological CaT followed by a physiological twitch. This model can open pathways for the future studies of human atrial electromechanics. KEY POINTS: With the availability of human atrial cellular data, interest in atrial-specific model integration has been enhanced. We have developed a detailed mathematical model of human atrial cardiomyocytes including the mechano-electric regulatory loop. The model has gone through calibration and evaluation phases against a wide collection of available human in-vitro data. The usefulness of the model for analysing clinical problems has been preliminaryly tested by simulating the increased incidence of Ca2+ transient and action potential alternans at high rates in post-operative atrial fibrillation condition. The model determines the possible role of mechano-electric feedback in alternans incidence, which can increase vulnerability to atrial arrhythmias by varying stretch levels. We found that our physiologically accurate description of Ca2+ handling can reproduce many experimental phenomena and can help to gain insights into the underlying pathophysiological mechanisms.

Progress on role of ion channels of cardiac fibroblasts in fibrosis

Cardiac fibrosis is defined as excessive deposition of extracellular matrix (ECM) in pathological conditions. Cardiac fibroblasts (CFs) activated by injury or inflammation differentiate into myofibroblasts (MFs) with secretory and contractile functions. In the fibrotic heart, MFs produce ECM which is composed mainly of collagen and is initially involved in maintaining tissue integrity. However, persistent fibrosis disrupts the coordination of excitatory contractile coupling, leading to systolic and diastolic dysfunction, and ultimately heart failure. Numerous studies have demonstrated that both voltage- and non-voltage-gated ion channels alter intracellular ion levels and cellular activity, contributing to myofibroblast proliferation, contraction, and secretory function. However, an effective treatment strategy for myocardial fibrosis has not been established. Therefore, this review describes the progress made in research related to transient receptor potential (TRP) channels, Piezo1, Ca²⁺ release-activated Ca²⁺ (CRAC) channels, voltage-gated Ca²⁺ channels (VGCCs), sodium channels, and potassium channels in myocardial fibroblasts with the aim of providing new ideas for treating myocardial fibrosis.

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.

Regional Diversities in Fibrogenesis Weighed as a Key Determinant for Atrial Arrhythmogenesis

Atrial fibrosis plays a key role in atrial myopathy, resulting in the genesis of atrial fibrillation (AF). The abnormal distribution of fibrotic tissue, electrical coupling, paracrine interactions, and biomechanical–electrical interactions have all been suggested as causes of fibrosis-related arrhythmogenesis. Moreover, the regional difference in fibrogenesis, specifically the left atrium (LA) exhibiting a higher arrhythmogenesis and level of fibrosis than the right atrium (RA) in AF, is a key contributor to atrial arrhythmogenesis. LA fibroblasts have greater profibrotic cellular activities than RA fibroblasts, but knowledge about the regional diversity of atrial regional fibrogenesis remains limited. This article provides a comprehensive review of research findings on the association between fibrogenesis and arrhythmogenesis from laboratory to clinical evidence and updates the current understanding of the potential mechanism underlying the difference in fibrogenesis between the LA and RA.

Stretch-activated current in human atrial myocytes and Na+ current and mechano-gated channels' current in myofibroblasts alter myocyte mechanical behavior: a computational study

Background

The activation of stretch-activated channels (SACs) in cardiac myocytes, which changes the phases of action potential repolarization, is proven to be highly efficient for the conversion of atrial fibrillation. The expression of Na⁺ current in myofibroblasts (Mfbs) regenerates myocytes’ action potentials, suggesting that Mfbs play an active role in triggering cardiac rhythm disturbances. Moreover, the excitation of mechano-gated channels (MGCs) in Mfbs depolarizes their membrane potential and contributes to the increased risk of post-infarct arrhythmia. Although these electrophysiological mechanisms have been largely known, the roles of these currents in cardiac mechanics are still debated. In this study, we aimed to investigate the mechanical influence of these currents via mathematical modeling. A novel mathematical model was developed by integrating models of human atrial myocyte (including the stretch-activated current, Ca²⁺–force relation, and mechanical behavior of a single segment) and Mfb (including our formulation of Na⁺ current and mechano-gated channels’ current). The effects of the changes in basic cycle length, number of coupled Mfbs and intercellular coupling conductance on myocyte mechanical properties were compared.

Results

Our results indicated that these three currents significantly regulated myocyte mechanical parameters. In isosarcometric contraction, these currents increased segment force by 13.8–36.6% and dropped element length by 12.1–31.5%. In isotonic contraction, there are 2.7–5.9% growth and 0.9–24% reduction. Effects of these currents on the extremum of myocyte mechanical parameters become more significant with the increase of basic cycle length, number of coupled Mfbs and intercellular coupling conductance.

Conclusions

The results demonstrated that stretch-activated current in myocytes and Na⁺ current and mechano-gated channels’ current in Mfbs significantly influenced myocyte mechanical behavior and should be considered in future cardiac mechanical mathematical modeling.

Role of Non-Myocyte Gap Junctions and Connexin Hemichannels in Cardiovascular Health and Disease: Novel Therapeutic Targets?

The heart is a complex organ composed of multiple cell types, including cardiomyocytes and different non-myocyte populations, all working closely together to determine the hearts properties and maintain normal cardiac function. Connexins are abundantly expressed proteins that form plasma membrane hemichannels and gap junctions between cells. Gap junctions are intracellular channels that allow for communication between cells, and in the heart they play a crucial role in cardiac conduction by coupling adjacent cardiomyocytes. Connexins are expressed in both cardiomyocytes and non-myocytes, including cardiac fibroblasts, endothelial cells, and macrophages. Non-myocytes are the largest population of cells in the heart, and therefore it is important to consider what roles connexins, hemichannels, and gap junctions play in these cell types. The aim of this review is to provide insight into connexin-based signalling in non-myocytes during health and disease, and highlight how targeting these proteins could lead to the development of novel therapies. We conclude that connexins in non-myocytes contribute to arrhythmias and adverse ventricular remodelling following myocardial infarction, and are associated with the initiation and development of atherosclerosis. Therefore, therapeutic interventions targeting these connexins represent an exciting new research avenue with great potential.

Simulation of Atrial Fibrosis Using Coupled Myocyte-Fibroblast Cellular and Human Atrial Models

Atrial fibrosis is characterized by expansion of extracellular matrix and increase in the number of fibroblasts which has been associated with the development and maintenance of atrial arrhythmias. However, the mechanisms how the fibrosis contributes to atrial arrhythmia remain incompletely understood. In this study, we used a proposed fibroblast model coupled with the human atrial myocyte to investigate the effects of fibrosis on atrial excitability and repolarization at both cellular and macroscopic levels. The 12-lead electrocardiogram (ECG) was also simulated to explore the index of clinical diagnosis for fibrosis. The simulation results showed that the fibrosis can modify action potential morphology of human atrial myocyte, slow down wave propagation, and have rate adaptation, thus causing the atrial electrical heterogeneity. The fibrosis alone was sufficient to cause arrhythmia, induce reentry wave, and result in low amplitude and wide P waves at normal heart rate and significant prolonged and inverse P waves at high heart rate. All these symptoms aggravated when the level of fibrosis increased. Our simulations demonstrated that fibrosis is the substrate of atrial arrhythmia and thereby may be a potential target in the treatment of atrial arrhythmias.

Effects of Na+ Current and Mechanogated Channels in Myofibroblasts on Myocyte Excitability and Repolarization

Fibrotic remodeling, characterized by fibroblast phenotype switching, is often associated with atrial fibrillation and heart failure. This study aimed to investigate the effects on electrotonic myofibroblast-myocyte (Mfb-M) coupling on cardiac myocytes excitability and repolarization of the voltage-gated sodium channels (VGSCs) and single mechanogated channels (MGCs) in human atrial Mfbs. Mathematical modeling was developed from a combination of (1) models of the human atrial myocyte (including the stretch activated ion channel current, I_SAC) and Mfb and (2) our formulation of currents through VGSCs (I_{Na_Mfb}) and MGCs (I_{MGC_Mfb}) based upon experimental findings. The effects of changes in the intercellular coupling conductance, the number of coupled Mfbs, and the basic cycle length on the myocyte action potential were simulated. The results demonstrated that the integration of I_SAC, I_{Na_Mfb}, and I_{MGC_Mfb} reduced the amplitude of the myocyte membrane potential (V_max) and the action potential duration (APD), increased the depolarization of the resting myocyte membrane potential (V_rest), and made it easy to trigger spontaneous excitement in myocytes. For Mfbs, significant electrotonic depolarizations were exhibited with the addition of I_{Na_Mfb} and I_{MGC_Mfb}. Our results indicated that I_SAC, I_{Na_Mfb}, and I_{MGC_Mfb} significantly influenced myocytes and Mfbs properties and should be considered in future cardiac pathological mathematical modeling.

A review of the literature on cardiac electrical activity between fibroblasts and myocytes

Myocardial injuries often lead to fibrotic deposition. This review presents evidence supporting the concept that fibroblasts in the heart electrically couple to myocytes.

Living cardiac tissue slices: an organotypic pseudo two-dimensional model for cardiac biophysics research

Living cardiac tissue slices, a pseudo two-dimensional (2D) preparation, have received less attention than isolated single cells, cell cultures, or Langendorff-perfused hearts in cardiac biophysics research. This is, in part, due to difficulties associated with sectioning cardiac tissue to obtain live slices. With moderate complexity, native cell-types, and well-preserved cell-cell electrical and mechanical interconnections, cardiac tissue slices have several advantages for studying cardiac electrophysiology. The trans-membrane potential (Vm) has, thus far, mainly been explored using multi-electrode arrays. Here, we combine tissue slices with optical mapping to monitor Vm and intracellular Ca(2+) concentration ([Ca(2+)]i). This combination opens up the possibility of studying the effects of experimental interventions upon action potential (AP) and calcium transient (CaT) dynamics in 2D, and with relatively high spatio-temporal resolution. As an intervention, we conducted proof-of-principle application of stretch. Mechanical stimulation of cardiac preparations is well-established for membrane patches, single cells and whole heart preparations. For cardiac tissue slices, it is possible to apply stretch perpendicular or parallel to the dominant orientation of cells, while keeping the preparation in a constant focal plane for fluorescent imaging of in-slice functional dynamics. Slice-to-slice comparison furthermore allows one to assess transmural differences in ventricular tissue responses to mechanical challenges. We developed and tested application of axial stretch to cardiac tissue slices, using a manually-controlled stretching device, and recorded Vm and [Ca(2+)]i by optical mapping before, during, and after application of stretch. Living cardiac tissue slices, exposed to axial stretch, show an initial shortening in both AP and CaT duration upon stretch application, followed in most cases by a gradual prolongation of AP and CaT duration during stretch maintained for up to 50 min. After release of sustained stretch, AP duration (APD) and CaT duration reverted to shorter values. Living cardiac tissue slices are a promising experimental model for the study of cardiac mechano-electric interactions. The methodology described here can be refined to achieve more accurate control over stretch amplitude and timing (e.g. using a computer-controlled motorised stage, or by synchronising electrical and mechanical events) and through monitoring of regional tissue deformation (e.g. by adding motion tracking).