Computational modeling of cardiac ventricular action potentials in rat and mouse: review

Developmental Changes in the Excitation-Contraction Mechanisms of the Ventricular Myocardium and Their Sympathetic Regulation in Small Experimental Animals

The developmental changes in the excitation-contraction mechanisms of the ventricular myocardium of small animals (guinea pig, rat, mouse) and their sympathetic regulation will be summarized. The action potential duration monotonically decreases during pre- and postnatal development in the rat and mouse, while in the guinea pig it decreases during the fetal stage but turns into an increase just before birth. Such changes can be attributed to changes in the repolarizing potassium currents. The T-tubule and the sarcoplasmic reticulum are scarcely present in the fetal cardiomyocyte, but increase during postnatal development. This causes a developmental shift in the Ca2+ handling from a sarcolemma-dependent mechanism to a sarcoplasmic reticulum-dependent mechanism. The sensitivity for beta-adrenoceptor-mediated positive inotropy decreases during early postnatal development, which parallels the increase in sympathetic nerve innervation. The alpha-adrenoceptor-mediated inotropy in the mouse changes from positive in the neonate to negative in the adult. This can be explained by the change in the excitation-contraction mechanism mentioned above. The shortening of the action potential duration enhances trans-sarcolemmal Ca2+ extrusion by the Na+-Ca2+ exchanger. The sarcoplasmic reticulum-dependent mechanism of contraction in the adult allows Na+-Ca2+ exchanger activity to cause negative inotropy, a mechanism not observed in neonatal myocardium. Such developmental studies would provide clues towards a more comprehensive understanding of cardiac function.

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.

Computational modeling of amylin-induced calcium dysregulation in rat ventricular cardiomyocytes

Hyperamylinemia is a condition that accompanies obesity and precedes type II diabetes, and it is characterized by above-normal blood levels of amylin, the pancreas-derived peptide. Human amylin oligomerizes easily and can deposit in the pancreas [1], brain [2], and heart [3], where they have been associated with calcium dysregulation. In the heart, accumulating evidence suggests that human amylin oligomers form moderately cation-selective [4,5] channels that embed in the cell sarcolemma (SL). The oligomers increase membrane conductance in a concentration-dependent manner [5], which is correlated with elevated cytosolic Ca²⁺. These findings motivate our core hypothesis that non-selective inward Ca²⁺ conduction afforded by human amylin oligomers increase cytosolic and sarcoplasmic reticulum (SR) Ca²⁺ load, which thereby magnifies intracellular Ca²⁺ transients. Questions remain however regarding the mechanism of amylin-induced Ca²⁺ dysregulation, including whether enhanced SL Ca²⁺ influx is sufficient to elevate cytosolic Ca²⁺ load [6], and if so, how might amplified Ca²⁺ transients perturb Ca²⁺-dependent cardiac pathways. To investigate these questions, we modified a computational model of cardiomyocytes Ca²⁺ signaling to reflect experimentally-measured changes in SL membrane permeation and decreased sarcoplasmic/endoplasmic reticulum calcium ATPase (SERCA) function stemming from acute and transgenic human amylin peptide exposure. With this model, we confirmed the hypothesis that increasing SL permeation alone was sufficient to enhance Ca²⁺ transient amplitudes. Our model indicated that amplified cytosolic transients are driven by increased Ca²⁺ loading of the SR and that greater fractional release may contribute to the Ca²⁺-dependent activation of calmodulin, which could prime the activation of myocyte remodeling pathways. Importantly, elevated Ca²⁺ in the SR and dyadic space collectively drive greater fractional SR Ca²⁺ release for human amylin expressing rats (HIP) and acute amylin-exposed rats (+Amylin) mice, which contributes to the inotropic rise in cytosolic Ca²⁺ transients. These findings suggest that increased membrane permeation induced by oligomeratization of amylin peptide in cell sarcolemma contributes to Ca²⁺ dysregulation in pre-diabetes.

Structural and functional plasticity in long-term cultures of adult ventricular myocytes

Cultured heart cells have long been valuable for characterizing biological mechanism and disease pathogenesis. However, these preparations have limitations, relating to immaturity in key properties like excitation-contraction coupling and β-adrenergic stimulation. Progressive attenuation of the latter is intimately related to pathogenesis and therapy in heart failure. Highly valuable would be a long-term culture system that emulates the structural and functional changes that accompany disease and development, while concurrently permitting ready access to underlying molecular events. Accordingly, we here produce functional monolayers of adult guinea-pig ventricular myocytes (aGPVMs) that can be maintained in long-term culture for several weeks. At baseline, these monolayers exhibit considerable myofibrillar organization and a significant contribution of sarcoplasmic reticular (SR) Ca(2+) release to global Ca(2+) transients. In terms of electrical signaling, these monolayers support propagated electrical activity and manifest monophasic restitution of action-potential duration and conduction velocity. Intriguingly, β-adrenergic stimulation increases chronotropy but not inotropy, indicating selective maintenance of β-adrenergic signaling. It is interesting that this overall phenotypic profile is not fixed, but can be readily enhanced by chronic electrical stimulation of cultures. This simple environmental cue significantly enhances myofibrillar organization as well as β-adrenergic sensitivity. In particular, the chronotropic response increases, and an inotropic effect now emerges, mimicking a reversal of the progression seen in heart failure. Thus, these aGPVM monolayer cultures offer a valuable platform for clarifying long elusive features of β-adrenergic signaling and its plasticity.

Pathophysiology of the cardiac late Na current and its potential as a drug target

A pathological increase in the late component of the cardiac Na(+) current, I(NaL), has been linked to disease manifestation in inherited and acquired cardiac diseases including the long QT variant 3 (LQT3) syndrome and heart failure. Disruption in I(NaL) leads to action potential prolongation, disruption of normal cellular repolarization, development of arrhythmia triggers, and propensity to ventricular arrhythmia. Attempts to treat arrhythmogenic sequelae from inherited and acquired syndromes pharmacologically with common Na(+) channel blockers (e.g. flecainide, lidocaine, and amiodarone) have been largely unsuccessful. This is due to drug toxicity and the failure of most current drugs to discriminate between the peak current component, chiefly responsible for single cell excitability and propagation in coupled tissue, and the late component (I(NaL)) of the Na(+) current. Although small in magnitude as compared to the peak Na(+) current (~1-3%), I(NaL) alters action potential properties and increases Na(+) loading in cardiac cells. With the increasing recognition that multiple cardiac pathological conditions share phenotypic manifestations of I(NaL) upregulation, there has been renewed interest in specific pharmacological inhibition of I(Na). The novel antianginal agent ranolazine, which shows a marked selectivity for late versus peak Na(+) current, may represent a novel drug archetype for targeted reduction of I(NaL). This article aims to review common pathophysiological mechanisms leading to enhanced I(NaL) in LQT3 and heart failure as prototypical disease conditions. Also reviewed are promising therapeutic strategies tailored to alter the molecular mechanisms underlying I(Na) mediated arrhythmia triggers.

Using computational modeling to predict arrhythmogenesis and antiarrhythmic therapy

The use of computational modeling to predict arrhythmia and arrhythmogensis is a relatively new field, but has nonetheless dramatically enhanced our understanding of the physiological and pathophysiological mechanisms that lead to arrhythmia. This review summarizes recent advances in the field of computational modeling approaches with a brief review of the evolution of cellular action potential models, and the incorporation of genetic mutations to understand fundamental arrhythmia mechanisms, including how simulations have revealed situation specific mechanisms leading to multiple phenotypes for the same genotype. The review then focuses on modeling drug blockade to understand how the less-than-intuitive effects some drugs have to either ameliorate or paradoxically exacerbate arrhythmia. Quantification of specific arrhythmia indicies are discussed at each spatial scale, from channel to tissue. The utility of hERG modeling to assess altered repolarization in response to drug blockade is also briefly discussed. Finally, insights gained from Ca(2+) dynamical modeling and EC coupling, neurohumoral regulation of cardiac dynamics, and cell signaling pathways are also reviewed.

Ca/calmodulin kinase II differentially modulates potassium currents

BACKGROUND

Potassium currents contribute to action potential duration (APD) and arrhythmogenesis. In heart failure, Ca/calmodulin-dependent protein kinase II (CaMKII) is upregulated and can alter ion channel regulation and expression.

METHODS AND RESULTS

We examine the influence of overexpressing cytoplasmic CaMKIIdelta(C), both acutely in rabbit ventricular myocytes (24-hour adenoviral gene transfer) and chronically in CaMKIIdelta(C)-transgenic mice, on transient outward potassium current (I(to)), and inward rectifying current (I(K1)). Acute and chronic CaMKII overexpression increases I(to,slow) amplitude and expression of the underlying channel protein K(V)1.4. Chronic but not acute CaMKII overexpression causes downregulation of I(to,fast), as well as K(V)4.2 and KChIP2, suggesting that K(V)1.4 expression responds faster and oppositely to K(V)4.2 on CaMKII activation. These amplitude changes were not reversed by CaMKII inhibition, consistent with CaMKII-dependent regulation of channel expression and/or trafficking. CaMKII (acute and chronic) greatly accelerated recovery from inactivation for both I(to) components, but these effects were acutely reversed by AIP (CaMKII inhibitor), suggesting that CaMKII activity directly accelerates I(to) recovery. Expression levels of I(K1) and Kir2.1 mRNA were downregulated by CaMKII overexpression. CaMKII acutely increased I(K1), based on inhibition by AIP (in both models). CaMKII overexpression in mouse prolonged APD (consistent with reduced I(to,fast) and I(K1)), whereas CaMKII overexpression in rabbit shortened APD (consistent with enhanced I(K1) and I(to,slow) and faster I(to) recovery). Computational models allowed discrimination of contributions of different channel effects on APD.

CONCLUSIONS

CaMKII has both acute regulatory effects and chronic expression level effects on I(to) and I(K1) with complex consequences on APD.

An overview of cardiac systems biology

The cardiac system has been a major target for intensive studies in the multi-scale modeling field for many years. Reproduction of the action potential and the ionic currents of single cardiomyocytes, as well as the construction of a whole organ model is well established. Still, there are major hurdles to overcome in creating a realistic and predictive functional cardiac model due to the lack of a profound understanding of the complex molecular interactions and their outcomes controlling both normal and pathological cardiophysiology. The recent advent of systems biology offers the conceptual and practical frameworks to tackle such biological complexities. This review provides an overview of major themes in the developing field of cardiac systems biology, summarizing some of the high-throughput experiments and strategies used to integrate the datasets, and various types of computational approaches used for developing useful quantitative models capable of predicting complex biological behavior.

Hybrid duplex: a novel method to study the contractile function of heterogeneous myocardium

In an earlier study, we experimentally mimicked the effects of mechanical interaction between different regions of the ventricular wall by allowing pairs of independently maintained cardiac muscle fibers to interact mechanically in series or in parallel. This simple physiological model of heterogeneous myocardium, which has been termed “duplex,” has provided new insight into basic effects of cardiac electromechanical heterogeneity. Here, we present a novel “hybrid duplex,” where one of the elements is an isolated cardiac muscle and the other a “virtual cardiac muscle.” The virtual muscle is represented by a computational model of cardiomyocyte electromechanical activity. We present in detail the computer-based digital control system that governs the mechanical interaction between virtual and biological muscle, the software used for data analysis, and working implementations of the model. Advantages of the hybrid duplex method are discussed, and experimental recordings are presented for illustration and as proof of the principle.