Genetic manipulation of calcium-handling proteins in cardiac myocytes. II. Mathematical modeling studies

Physiology of intracellular calcium buffering

Calcium signaling underlies much of physiology. Almost all the Ca2+ in the cytoplasm is bound to buffers, with typically only ∼1% being freely ionized at resting levels in most cells. Physiological Ca2+ buffers include small molecules and proteins, and experimentally Ca2+ indicators will also buffer calcium. The chemistry of interactions between Ca2+ and buffers determines the extent and speed of Ca2+ binding. The physiological effects of Ca2+ buffers are determined by the kinetics with which they bind Ca2+ and their mobility within the cell. The degree of buffering depends on factors such as the affinity for Ca2+, the Ca2+ concentration, and whether Ca2+ ions bind cooperatively. Buffering affects both the amplitude and time course of cytoplasmic Ca2+ signals as well as changes of Ca2+ concentration in organelles. It can also facilitate Ca2+ diffusion inside the cell. Ca2+ buffering affects synaptic transmission, muscle contraction, Ca2+ transport across epithelia, and the killing of bacteria. Saturation of buffers leads to synaptic facilitation and tetanic contraction in skeletal muscle and may play a role in inotropy in the heart. This review focuses on the link between buffer chemistry and function and how Ca2+ buffering affects normal physiology and the consequences of changes in disease. As well as summarizing what is known, we point out the many areas where further work is required.

Dysregulation of Calcium Handling in Duchenne Muscular Dystrophy-Associated Dilated Cardiomyopathy: Mechanisms and Experimental Therapeutic Strategies

: Duchenne muscular dystrophy (DMD) is an X-linked recessive disease resulting in the loss of dystrophin, a key cytoskeletal protein in the dystrophin-glycoprotein complex. Dystrophin connects the extracellular matrix with the cytoskeleton and stabilizes the sarcolemma. Cardiomyopathy is prominent in adolescents and young adults with DMD, manifesting as dilated cardiomyopathy (DCM) in the later stages of disease. Sarcolemmal instability, leading to calcium mishandling and overload in the cardiac myocyte, is a key mechanistic contributor to muscle cell death, fibrosis, and diminished cardiac contractile function in DMD patients. Current therapies for DMD cardiomyopathy can slow disease progression, but they do not directly target aberrant calcium handling and calcium overload. Experimental therapeutic targets that address calcium mishandling and overload include membrane stabilization, inhibition of stretch-activated channels, ryanodine receptor stabilization, and augmentation of calcium cycling via modulation of the Serca2a/phospholamban (PLN) complex or cytosolic calcium buffering. This paper addresses what is known about the mechanistic basis of calcium mishandling in DCM, with a focus on DMD cardiomyopathy. Additionally, we discuss currently utilized therapies for DMD cardiomyopathy, and review experimental therapeutic strategies targeting the calcium handling defects in DCM and DMD cardiomyopathy.

A new myofilament contraction model with ATP consumption for ventricular cell model

A new contraction model of cardiac muscle was developed by combining previously described biochemical and biophysical models. The biochemical component of the new contraction model represents events in the presence of Ca2+-crossbridge attachment and power stroke following inorganic phosphate release, detachment evoked by the replacement of ADP by ATP, ATP hydrolysis, and recovery stroke. The biophysical component focuses on Ca2+ activation and force (F b) development assuming an equivalent crossbridge. The new model faithfully incorporates the major characteristics of the biochemical and biophysical models, such as F b activation by transient Ca2+ ([Ca2+]-F b), [Ca2+]-ATP hydrolysis relations, sarcomere length-F b, and F b recovery after jumps in length under the isometric mode and upon sarcomere shortening after a rapid release of mechanical load under the isotonic mode together with the load-velocity relationship. ATP consumption was obtained for all responses. When incorporated in a ventricular cell model, the contraction model was found to share approximately 60% of the total ATP usage in the cell model.

Effects of Modified Parvalbumin EF-Hand Motifs on Cardiac Myocyte Contractile Function

Cardiac gene delivery of parvalbumin (Parv), an EF-hand Ca(2+) buffer, has been studied as a therapeutic strategy for diastolic heart failure, in which slow Ca(2+) reuptake is an important contributor. A limitation of wild-type (WT) Parv is the significant trade-off between faster relaxation and blunted contraction amplitude, occurring because WT-Parv sequesters Ca(2+) too early in the cardiac cycle and prematurely truncates sarcomere shortening in the facilitation of rapid relaxation. We recently demonstrated that an E → Q substitution (ParvE101Q) at amino acid 12 of the EF-hand Ca(2+)/Mg(2+) binding loop disrupts bidentate Ca(2+) binding, reducing Ca(2+) affinity by 99-fold and increasing Mg(2+) affinity twofold. ParvE101Q caused faster relaxation and not only preserved contractility, but unexpectedly increased it above untreated myocytes. To gain mechanistic insight into the increased contractility, we focused here on amino acid 12 of the EF-hand motif. We introduced an E → D substitution (ParvE101D) at this site, which converts bidentate Ca(2+) coordination to monodentate coordination. ParvE101D decreased Ca(2+) affinity by 114-fold and increased Mg(2+) affinity 28-fold compared to WT-Parv. ParvE101D increased contraction amplitude compared to both untreated myocytes and myocytes with ParvE101Q, with limited improvement in relaxation. Additionally, ParvE101D increased spontaneous contractions after pacing stress. ParvE101D also increased Ca(2+) transient peak height and was diffusely localized around the Z-line of the sarcomere, suggesting a Ca(2+)-dependent mechanism of enhanced contractility. Sarcoplasmic reticulum Ca(2+) load was not changed with ParvE101D, but postpacing Ca(2+) waves were increased. Together, these data show that inverted Ca(2+)/Mg(2+) binding affinities of ParvE101D increase myocyte contractility through a Ca(2+)-dependent mechanism without altering sarcoplasmic reticulum Ca(2+) load and by increasing unstimulated contractions and Ca(2+) waves. ParvE101D provides mechanistic insight into how changes in the Ca(2+)/Mg(2+) binding affinities of parvalbumin's EF-hand motif alter function of cardiac myocytes. These data are informative in developing new Ca(2+) buffering strategies for the failing heart.

Myokit: A simple interface to cardiac cellular electrophysiology

Myokit is a new powerful and versatile software tool for modeling and simulation of cardiac cellular electrophysiology. Myokit consists of an easy-to-read modeling language, a graphical user interface, single and multi-cell simulation engines and a library of advanced analysis tools accessible through a Python interface. Models can be loaded from Myokit's native file format or imported from CellML. Model export is provided to C, MATLAB, CellML, CUDA and OpenCL. Patch-clamp data can be imported and used to estimate model parameters. In this paper, we review existing tools to simulate the cardiac cellular action potential to find that current tools do not cater specifically to model development and that there is a gap between easy-to-use but limited software and powerful tools that require strong programming skills from their users. We then describe Myokit's capabilities, focusing on its model description language, simulation engines and import/export facilities in detail. Using three examples, we show how Myokit can be used for clinically relevant investigations, multi-model testing and parameter estimation in Markov models, all with minimal programming effort from the user. This way, Myokit bridges a gap between performance, versatility and user-friendliness.

Cardiac-specific deletion of protein phosphatase 1β promotes increased myofilament protein phosphorylation and contractile alterations

There are 3 protein phosphatase 1 (PP1) catalytic isoforms (α, β and γ) encoded within the mammalian genome. These 3 gene products share ~90% amino acid homology within their catalytic domains but each has unique N- and C-termini that likely underlie distinctive subcellular localization or functionality. In this study, we assessed the effect associated with the loss of each PP1 isoform in the heart using a conditional Cre-loxP targeting approach in mice. Ppp1ca-loxP, Ppp1cb-loxP and Ppp1cc-loxP alleles were crossed with either an Nkx2.5-Cre knock-in containing allele for early embryonic deletion or a tamoxifen inducible α-myosin heavy chain (αMHC)-MerCreMer transgene for adult and cardiac-specific deletion. We determined that while deletion of Ppp1ca (PP1α) or Ppp1cc (PP1γ) had little effect on the whole heart, deletion of Ppp1cb (PP1β) resulted in concentric remodeling of the heart, interstitial fibrosis and contractile dysregulation, using either the embryonic or adult-specific Cre-expressing alleles. However, myocytes isolated from Ppp1cb deleted hearts surprisingly showed enhanced contractility. Mechanistically we found that deletion of any of the 3 PP1 gene-encoding isoforms had no effect on phosphorylation of phospholamban, nor were Ca(2+) handling dynamics altered in adult myocytes from Ppp1cb deleted hearts. However, the loss of Ppp1cb from the heart, but not Ppp1ca or Ppp1cc, resulted in elevated phosphorylation of myofilament proteins such as myosin light chain 2 and cardiac myosin binding protein C, consistent with an enriched localization profile of this isoform to the sarcomeres. These results suggest a unique functional role for the PP1β isoform in affecting cardiac contractile function.

Interpreting genetic effects through models of cardiac electromechanics

Multiscale models of cardiac electromechanics are being increasingly focused on understanding how genetic variation and environment underpin multiple disease states. In this paper we review the current state of the art in both the development of specific models and the physiological insights they have produced. This growing research body includes the development of models for capturing the effects of changes in function in both single and multiple proteins in both specific expression systems and in vivo contexts. Finally, the potential for using this approach for ultimately predicting phenotypes from genetic sequence information is discussed.

Diastolic dysfunction and thin filament dysregulation resulting from excitation-contraction uncoupling in a mouse model of restrictive cardiomyopathy

Restrictive cardiomyopathy (RCM) has been linked to mutations in the thin filament regulatory protein cardiac troponin I (cTnI). As the pathogenesis of RCM from genotype to clinical phenotype is not fully understood, transgenic (Tg) mice were generated with cardiac specific expression of an RCM-linked missense mutation (R193H) in cTnI. R193H Tg mouse hearts with 15% stoichiometric replacement had smaller hearts and significantly elevated end diastolic pressures (EDP) in vivo. Using a unique carbon microfiber-based assay, membrane intact R193H adult cardiac myocytes generated higher passive tensions across a range of physiologic sarcomere lengths resulting in significant Ca(2+) independent cellular diastolic tone that was manifest in vivo as elevated organ-level EDP. Sarcomere relaxation and Ca(2+) decay was uncoupled in isolated R193H Tg adult myocytes due to the increase in myofilament Ca(2+) sensitivity of tension, decreased passive compliance of the sarcomere, and adaptive in vivo changes in which phospholamban (PLN) content was decreased. Further evidence of Ca(2+) and mechanical uncoupling in R193H Tg myocytes was demonstrated by the biphasic response of relaxation to increased pacing frequency versus the negative staircase seen with Ca(2+) decay. In comparison, non-transgenic myocyte relaxation closely paralleled the accelerated Ca(2+) decay. Ca(2+) transient amplitude was also significantly blunted in R193H Tg myocytes despite normal mechanical shortening resulting in myocyte hypercontractility when compared to non-transgenics. These results identify for the first time that a single point mutation in cTnI, R193H, directly causes elevated EDP due to a myocyte intrinsic loss of compliance independent of Ca(2+) cycling or altered cardiac morphology. The compound influence of impaired relaxation and elevated EDP represents a clinically severe form of diastolic dysfunction similar to the hemodynamic state documented in RCM patients.

How the Hodgkin-Huxley equations inspired the Cardiac Physiome Project

Early modelling of cardiac cells (1960-1980) was based on extensions of the Hodgkin-Huxley nerve axon equations with additional channels incorporated, but after 1980 it became clear that processes other than ion channel gating were also critical in generating electrical activity. This article reviews the development of models representing almost all cell types in the heart, many different species, and the software tools that have been created to facilitate the cardiac Physiome Project.

Parvalbumin isoforms for enhancing cardiac diastolic function

Diastolic heart failure (DHF), characterized by depressed myocardial relaxation performance and poor ventricular filling, is a distinct form of heart failure accounting for nearly half of the heart failure patients with otherwise normal systolic performance. Defective intracellular calcium (Ca2+) cycling is an important mechanism underlying impaired relaxation in DHF. Recently, genetic manipulation of Ca2+ handling proteins in cardiac myocytes has been explored for its potential therapeutic application in DHF. Specifically, ectopic expression of the skeletal muscle Ca2+ binding protein parvalbumin (Parv) has been shown to accelerate myocardial relaxation in vitro and in vivo. Parv acts as a unique "delayed" Ca2+ buffer during diastole by promoting Ca2+ transient decay and sequestration and corrects diastolic dysfunction in an energy-independent manner. This brief review summarizes the rationale and development of Parv gene transfer approaches for DHF, and in particular, discusses the divergent effects of Parv isoforms on cardiac myocyte Ca2+ handling and contractile function with the long-range goal of alleviating diastolic dysfunction in DHF.