A compartmentalized mathematical model of mouse atrial myocytes

EAD Mechanisms in Hypertrophic Mouse Ventricular Myocytes: Insights from a Compartmentalized Mathematical Model

Transverse aortic constriction (TAC) is one of the experimental mouse models that are designed to investigate cardiac hypertrophy and heart failure. Most of the studies with this model are devoted to the stage of developed heart failure. However, several studies of the early stages (hypertrophy after 1 week of TAC) of this disease found significant changes in the β-adrenergic system, electrical activity, and Ca2+ dynamics in mouse ventricular myocytes. To provide a quantitative description of cardiac hypertrophy, we developed a new compartmentalized mathematical model of hypertrophic mouse ventricular myocytes for the early stage after the TAC procedure. The model described the changes in cell geometry, action potentials, [Ca2+]i transients, and β1- and β2-adrenergic signaling systems. We also showed that the hypertrophic myocytes demonstrated early afterdepolarizations (EADs) upon stimulation with isoproterenol at relatively long stimulation periods. Simulation of the hypertrophic myocyte activities revealed that the synergistic effects of the late Na+ current, the L-type Ca2+ current, and the T-type Ca2+ current were responsible for the initiation of EADs. The mechanisms of EAD and its suppression were investigated and sensitivity analysis was performed. Simulation results obtained with the hypertrophic cell model were compared to those from the normal ventricular myocytes. The developed mathematical model can be used for the explanation of the existing experimental data, for the development of the models for other hypertrophic phenotypes, and to make experimentally testable predictions of a hypertrophic myocyte's behavior.

Caveolae-associated cAMP/Ca2+-mediated mechano-chemical signal transduction in mouse atrial myocytes

Caveolae are tiny invaginations in the sarcolemma that buffer extra membrane and contribute to mechanical regulation of cellular function. While the role of caveolae in membrane mechanosensation has been studied predominantly in non-cardiomyocyte cells, caveolae contribution to cardiac mechanotransduction remains elusive. Here, we studied the role of caveolae in the regulation of Ca2+ signaling in atrial cardiomyocytes. In Langendorff-perfused mouse hearts, atrial pressure/volume overload stretched atrial myocytes and decreased caveolae density. In isolated cells, caveolae were disrupted through hypotonic challenge that induced a temporal (<10 min) augmentation of Ca2+ transients and caused a rise in Ca2+ spark activity. Similar changes in Ca2+ signaling were observed after chemical (methyl-β-cyclodextrin) and genetic ablation of caveolae in cardiac-specific conditional caveolin-3 knock-out mice. Acute disruption of caveolae, both mechanical and chemical, led to the elevation of cAMP level in the cell interior, and cAMP-mediated augmentation of protein kinase A (PKA)-phosphorylated ryanodine receptors (at Ser2030 and Ser2808). Caveolae-mediated stimulatory effects on Ca2+ signaling were abolished via inhibition of cAMP production by adenyl cyclase antagonists MDL12330 and SQ22536, or reduction of PKA activity by H-89. A compartmentalized mathematical model of mouse atrial myocytes linked the observed changes to a microdomain-specific decrease in phosphodiesterase activity, which disrupted cAMP signaling and augmented PKA activity. Our findings add a new dimension to cardiac mechanobiology and highlight caveolae-associated cAMP/PKA-mediated phosphorylation of Ca2+ handling proteins as a novel component of mechano-chemical feedback in atrial myocytes.

Regulation of cardiac function by cAMP nanodomains

Cyclic adenosine monophosphate (cAMP) is a diffusible intracellular second messenger that plays a key role in the regulation of cardiac function. In response to the release of catecholamines from sympathetic terminals, cAMP modulates heart rate and the strength of contraction and ease of relaxation of each heartbeat. At the same time, cAMP is involved in the response to a multitude of other hormones and neurotransmitters. A sophisticated network of regulatory mechanisms controls the temporal and spatial propagation of cAMP, resulting in the generation of signaling nanodomains that enable the second messenger to match each extracellular stimulus with the appropriate cellular response. Multiple proteins contribute to this spatiotemporal regulation, including the cAMP-hydrolyzing phosphodiesterases (PDEs). By breaking down cAMP to a different extent at different locations, these enzymes generate subcellular cAMP gradients. As a result, only a subset of the downstream effectors is activated and a specific response is executed. Dysregulation of cAMP compartmentalization has been observed in cardiovascular diseases, highlighting the importance of appropriate control of local cAMP signaling. Current research is unveiling the molecular organization underpinning cAMP compartmentalization, providing original insight into the physiology of cardiac myocytes and the alteration associated with disease, with the potential to uncover novel therapeutic targets. Here, we present an overview of the mechanisms that are currently understood to be involved in generating cAMP nanodomains and we highlight the questions that remain to be answered.

A compartmentalized mathematical model of the β1- and β2-adrenergic signaling systems in ventricular myocytes from mouse in heart failure

Mouse models of heart failure are extensively used to research human cardiovascular diseases. In particular, one of the most common is the mouse model of heart failure resulting from transverse aortic constriction (TAC). Despite this, there are no comprehensive compartmentalized mathematical models that describe the complex behavior of the action potential, [Ca²⁺]_i transients, and their regulation by β₁- and β₂-adrenergic signaling systems in failing mouse myocytes. In this paper, we develop a novel compartmentalized mathematical model of failing mouse ventricular myocytes after TAC procedure. The model describes well the cell geometry, action potentials, [Ca²⁺]_i transients, and β₁- and β₂-adrenergic signaling in the failing cells. Simulation results obtained with the failing cell model are compared with those from the normal ventricular myocytes. Exploration of the model reveals the sarcoplasmic reticulum Ca²⁺ load mechanisms in failing ventricular myocytes. We also show a larger susceptibility of the failing myocytes to early and delayed afterdepolarizations and to a proarrhythmic behavior of Ca²⁺ dynamics upon stimulation with isoproterenol. The mechanisms of the proarrhythmic behavior suppression are investigated and sensitivity analysis is performed. The developed model can explain the existing experimental data on failing mouse ventricular myocytes and make experimentally testable predictions of a failing myocyte's behavior.

Evidence for a Physiological Role of T-Type Ca Channels in Ventricular Cardiomyocytes of Adult Mice

T-type Ca channels are strongly expressed and important in the developing heart. In the adult heart, these channels play a significant role in pacemaker tissues, but there is uncertainty about their presence and physiological relevance in the working myocardium. Here, we show that the T-type Ca channel isoforms Cav3.1 and Cav3.2 are expressed at a protein level in ventricular cardiomyocytes from healthy adult C57/BL6 mice. Myocytes isolated from adult wild-type and Cav3.2 KO mice showed considerable whole cell T-type Ca currents under beta-adrenergic stimulation with isoprenaline. We further show that the detectability of basal T-type Ca currents in murine wild-type cardiomyocytes depends on the applied experimental conditions. Together, these findings reveal the presence of functional T-type Ca channels in the membrane of ventricular myocytes. In addition, electrically evoked Ca release from the sarcoplasmic reticulum was significantly impaired in Cav3.2 KO compared to wild-type cardiomyocytes. Our work implies a physiological role of T-type Ca channels in the healthy adult murine ventricular working myocardium.

Bifurcations and Proarrhythmic Behaviors in Cardiac Electrical Excitations

The heart is a hierarchical dynamic system consisting of molecules, cells, and tissues, and acts as a pump for blood circulation. The pumping function depends critically on the preceding electrical activity, and disturbances in the pattern of excitation propagation lead to cardiac arrhythmia and pump failure. Excitation phenomena in cardiomyocytes have been modeled as a nonlinear dynamical system. Because of the nonlinearity of excitation phenomena, the system dynamics could be complex, and various analyses have been performed to understand the complex dynamics. Understanding the mechanisms underlying proarrhythmic responses in the heart is crucial for developing new ways to prevent and control cardiac arrhythmias and resulting contractile dysfunction. When the heart changes to a pathological state over time, the action potential (AP) in cardiomyocytes may also change to a different state in shape and duration, often undergoing a qualitative change in behavior. Such a dynamic change is called bifurcation. In this review, we first summarize the contribution of ion channels and transporters to AP formation and our knowledge of ion-transport molecules, then briefly describe bifurcation theory for nonlinear dynamical systems, and finally detail its recent progress, focusing on the research that attempts to understand the developing mechanisms of abnormal excitations in cardiomyocytes from the perspective of bifurcation phenomena.

Seeking the exclusive binding region of phenylalkylamine derivatives on human T-type calcium channels via homology modeling and molecular dynamics simulation approach

Pharmaceutical features of phenylalkylamine derivatives (PAAs) binding to calcium channels have been studied extensively in the past decades. Only a few PAAs have the binding specificity on calcium channels, for example, NNC 55‐0396. Here, we created the homology models of human Ca_v3.2, Ca_v3.3 and use them as a receptor on the rigid docking tests. The nonspecific calcium channel blocker mibefradil showed inconsistent docking preference across four domains; however, NNC 55‐0396 had a unique binding pattern on domain II specifically. The subsequent molecular dynamics (MD) simulations identified that Ca_v3.1, Ca_v3.2, and Ca_v3.3 share domain II when Ca²⁺ appearing in the neighbor region of selective filters (SFs). Moreover, free‐energy perturbation analysis suggests single mutation of lysine at P‐loop domain III, or threonine at the P‐loop domain II largely reduced the total amount of hydration‐free energy in the system. All these findings suggest that P‐loop and segment six domain II in the T‐type calcium channels (TCCs) are crucial for attracting the PAAs with specificity as the antagonist.

Evolution of mathematical models of cardiomyocyte electrophysiology

Advanced computational techniques and mathematical modeling have become more and more important to the study of cardiac electrophysiology. In this review, we provide a brief history of the evolution of cardiomyocyte electrophysiology models and highlight some of the most important ones that had a major impact on our understanding of the electrical activity of the myocardium and associated transmembrane ion fluxes in normal and pathological states. We also present the use of these models in the study of various arrhythmogenesis mechanisms, particularly the integration of experimental pharmacology data into advanced humanized models for in silico proarrhythmogenic risk prediction as an essential component of the Comprehensive in vitro Proarrhythmia Assay (CiPA) drug safety paradigm.

A Mathematical Model of the Mouse Atrial Myocyte With Inter-Atrial Electrophysiological Heterogeneity

Biophysically detailed mathematical models of cardiac electrophysiology provide an alternative to experimental approaches for investigating possible ionic mechanisms underlying the genesis of electrical action potentials and their propagation through the heart. The aim of this study was to develop a biophysically detailed mathematical model of the action potentials of mouse atrial myocytes, a popular experimental model for elucidating molecular and cellular mechanisms of arrhythmogenesis. Based on experimental data from isolated mouse atrial cardiomyocytes, a set of mathematical equations for describing the biophysical properties of membrane ion channel currents, intracellular Ca²⁺ handling, and Ca²⁺-calmodulin activated protein kinase II and β-adrenergic signaling pathways were developed. Wherever possible, membrane ion channel currents were modeled using Markov chain formalisms, allowing detailed representation of channel kinetics. The model also considered heterogeneous electrophysiological properties between the left and the right atrial cardiomyocytes. The developed model was validated by its ability to reproduce the characteristics of action potentials and Ca²⁺ transients, matching quantitatively to experimental data. Using the model, the functional roles of four K⁺ channel currents in atrial action potential were evaluated by channel block simulations, results of which were quantitatively in agreement with existent experimental data. To conclude, this newly developed model of mouse atrial cardiomyocytes provides a powerful tool for investigating possible ion channel mechanisms of atrial electrical activity at the cellular level and can be further used to investigate mechanisms underlying atrial arrhythmogenesis.