Simulation of the effects of moderate stimulation/inhibition of the β1-adrenergic signaling system and its components in mouse ventricular myocytes

Depressed β-adrenergic inotropic responsiveness and intracellular calcium handling abnormalities in Duchenne Muscular Dystrophy patients' induced pluripotent stem cell-derived cardiomyocytes

Duchenne muscular dystrophy (DMD), caused by mutations in the dystrophin gene, is an X-linked disease affecting male and rarely adult heterozygous females, resulting in death by the late 20s to early 30s. Previous studies reported depressed left ventricular function in DMD patients which may result from deranged intracellular Ca2+ -handling. To decipher the mechanism(s) underlying the depressed LV function, we tested the hypothesis that iPSC-CMs generated from DMD patients feature blunted positive inotropic response to β-adrenergic stimulation. To test the hypothesis, [Ca2+ ]i transients and contractions were recorded from healthy and DMD-CMs. While in healthy CMs (HC) isoproterenol caused a prominent positive inotropic effect, DMD-CMs displayed a blunted inotropic response. Next, we tested the functionality of the sarcoplasmic reticulum (SR) by measuring caffeine-induced Ca2+ release. In contrast to HC, DMD-CMs exhibited reduced caffeine-induced Ca2+ signal amplitude and recovery time. In support of the depleted SR Ca2+ stores hypothesis, in DMD-CMs the negative inotropic effects of ryanodine and cyclopiazonic acid were smaller than in HC. RNA-seq analyses demonstrated that in DMD CMs the RNA-expression levels of specific subunits of the L-type calcium channel, the β1-adrenergic receptor (ADRβ1) and adenylate cyclase were down-regulated by 3.5-, 2.8- and 3-fold, respectively, which collectively contribute to the depressed β-adrenergic responsiveness.

Computational modeling approaches to cAMP/PKA signaling in cardiomyocytes

The cAMP/PKA pathway is a fundamental regulator of excitation-contraction coupling in cardiomyocytes. Activation of cAMP has a variety of downstream effects on cardiac function including enhanced contraction, accelerated relaxation, adaptive stress response, mitochondrial regulation, and gene transcription. Experimental advances have shed light on the compartmentation of cAMP and PKA, which allow for control over the varied targets of these second messengers and is disrupted in heart failure conditions. Computational modeling is an important tool for understanding the spatial and temporal complexities of this system. In this review article, we outline the advances in computational modeling that have allowed for deeper understanding of cAMP/PKA dynamics in the cardiomyocyte in health and disease, and explore new modeling frameworks that may bring us closer to a more complete understanding of this system. We outline various compartmental and spatial signaling models that have been used to understand how β-adrenergic signaling pathways function in a variety of simulation conditions. We also discuss newer subcellular models of cardiovascular function that may be used as templates for the next phase of computational study of cAMP and PKA in the heart, and outline open challenges which are important to consider in future models.

Mathematical model for β1-adrenergic regulation of the mouse ventricular myocyte contraction

The β1-adrenergic regulation of cardiac myocyte contraction plays an important role in regulating heart function. Activation of this system leads to an increased heart rate and stronger myocyte contraction. However, chronic stimulation of the β1-adrenergic signaling system can lead to cardiac hypertrophy and heart failure. To understand the mechanisms of action of β1-adrenoceptors, a mathematical model of cardiac myocyte contraction that includes the β1-adrenergic system was developed and studied. The model was able to simulate major experimental protocols for measurements of steady-state force-calcium relationships, cross-bridge release rate and force development rate, force-velocity relationship, and force redevelopment rate. It also reproduced quite well frequency and isoproterenol dependencies for intracellular Ca2+ concentration ([Ca2+]i) transients, total contraction force, and sarcomere shortening. The mathematical model suggested the mechanisms of increased contraction force and myocyte shortening on stimulation of β1-adrenergic receptors is due to phosphorylation of troponin I and myosin-binding protein C and increased [Ca2+]i transient resulting from activation of the β1-adrenergic signaling system. The model was used to simulate work-loop contractions and estimate the power during the cardiac cycle as well as the effects of 4-aminopyridine and tedisamil on the myocyte contraction. The developed mathematical model can be used further for simulations of contraction of ventricular myocytes from genetically modified mice and myocytes from mice with chronic cardiac diseases.

NEW & NOTEWORTHY A new mathematical model of mouse ventricular myocyte contraction that includes the β1-adrenergic system was developed. The model simulated major experimental protocols for myocyte contraction and predicted the effects of 4-aminopyridine and tedisamil on the myocyte contraction. The model also allowed for simulations of work-loop contractions and estimation of the power during the cardiac cycle.

A compartmentalized mathematical model of mouse atrial myocytes

Various experimental mouse models are extensively used to research human diseases, including atrial fibrillation, the most common cardiac rhythm disorder. Despite this, there are no comprehensive mathematical models that describe the complex behavior of the action potential and [Ca²⁺]_i transients in mouse atrial myocytes. Here, we develop a novel compartmentalized mathematical model of mouse atrial myocytes that combines the action potential, [Ca²⁺]_i dynamics, and β-adrenergic signaling cascade for a subpopulation of right atrial myocytes with developed transverse-axial tubule system. The model consists of three compartments related to β-adrenergic signaling (caveolae, extracaveolae, and cytosol) and employs local control of Ca²⁺ release. It also simulates ionic mechanisms of action potential generation and describes atrial-specific Ca²⁺ handling as well as frequency dependences of the action potential and [Ca²⁺]_i transients. The model showed that the T-type Ca²⁺ current significantly affects the later stage of the action potential, with little effect on [Ca²⁺]_i transients. The block of the small-conductance Ca²⁺-activated K⁺ current leads to a prolongation of the action potential at high intracellular Ca²⁺. Simulation results obtained from the atrial model cells were compared with those from ventricular myocytes. The developed model represents a useful tool to study complex electrical properties in the mouse atria and could be applied to enhance the understanding of atrial physiology and arrhythmogenesis.NEW & NOTEWORTHY A new compartmentalized mathematical model of mouse right atrial myocytes was developed. The model simulated action potential and Ca²⁺ dynamics at baseline and after stimulation of the β-adrenergic signaling system. Simulations showed that the T-type Ca²⁺ current markedly prolonged the later stage of atrial action potential repolarization, with a minor effect on [Ca²⁺]_i transients. The small-conductance Ca²⁺-activated K⁺ current block resulted in prolongation of the action potential only at the relatively high intracellular Ca²⁺.

Oxidative and glycolytic skeletal muscles show marked differences in gene expression profile in Chinese Qingyuan partridge chickens

Oxidative and glycolytic myofibers have different structures and metabolic characteristics and their ratios are important in determining poultry meat quality. However, the molecular mechanisms underlying their differences are unclear. In this study, global gene expression profiling was conducted in oxidative skeletal muscle (obtained from the soleus, or SOL) and glycolytic skeletal muscle (obtained from the extensor digitorum longus, or EDL) of Chinese Qingyuan partridge chickens, using the Agilent Chicken Gene Expression Chip. A total of 1224 genes with at least 2-fold differences were identified (P < 0.05), of which 654 were upregulated and 570 were downregulated in SOL. GO, KEGG pathway, and co-expressed gene network analyses suggested that PRKAG3, ATP2A2, and PPARGC1A might play important roles in myofiber composition. The function of PPARGC1A gene was further validated. PPARGC1A mRNA expression levels were higher in SOL than in EDL muscles throughout the early postnatal development stages. In myoblast cells, shRNA knockdown of PPARGC1A significantly inhibited some muscle development and transition-related genes, including PPP3CA, MEF2C, and SM (P < 0.01 or P < 0.05), and significantly upregulated the expression of FWM (P < 0.05). Our study demonstrates strong transcriptome differences between oxidative and glycolytic myofibers, and the results suggest that PPARGC1A is a key gene involved in chicken myofiber composition and transition.

Distinct physiological effects of β1- and β2-adrenoceptors in mouse ventricular myocytes: insights from a compartmentalized mathematical model

The β1- and β2-adrenergic signaling systems play different roles in the functioning of cardiac cells. Experimental data show that the activation of the β1-adrenergic signaling system produces significant inotropic, lusitropic, and chronotropic effects in the heart, whereas the effects of the β2-adrenergic signaling system is less apparent. In this paper, a comprehensive compartmentalized experimentally based mathematical model of the combined β1- and β2-adrenergic signaling systems in mouse ventricular myocytes is developed to simulate the experimental findings and make testable predictions of the behavior of the cardiac cells under different physiological conditions. Simulations describe the dynamics of major signaling molecules in different subcellular compartments; kinetics and magnitudes of phosphorylation of ion channels, transporters, and Ca2+ handling proteins; modifications of action potential shape and duration; and [Ca2+]i and [Na+]i dynamics upon stimulation of β1- and β2-adrenergic receptors (β1- and β2-ARs). The model reveals physiological conditions when β2-ARs do not produce significant physiological effects and when their effects can be measured experimentally. Simulations demonstrated that stimulation of β2-ARs with isoproterenol caused a marked increase in the magnitude of the L-type Ca2+ current, [Ca2+]i transient, and phosphorylation of phospholamban only upon additional application of pertussis toxin or inhibition of phosphodiesterases of type 3 and 4. The model also made testable predictions of the changes in magnitudes of [Ca2+]i and [Na+]i fluxes, the rate of decay of [Na+]i concentration upon both combined and separate stimulation of β1- and β2-ARs, and the contribution of phosphorylation of PKA targets to the changes in the action potential and [Ca2+]i transient.