Dataset Information

Quantitative analysis of variability in an integrated model of human ventricular electrophysiology and β-adrenergic signaling.

ABSTRACT: In ventricular myocytes, stimulation of β-adrenergic receptors activates critical cardiac signaling pathways, leading to shorter action potentials and increased contraction strength during the "fight-or-flight" response. These changes primarily result, at the cellular level, from the coordinated phosphorylation of multiple targets by protein kinase A. Although mathematical models of the intracellular signaling downstream of β-adrenergic receptor activation have previously been described, only a limited number of studies have explored quantitative interactions between intracellular signaling and electrophysiology in human ventricular myocytes. Accordingly, our objective was to develop an integrative mathematical model of β-adrenergic receptor signaling, electrophysiology, and intracellular calcium (Ca²⁺) handling in the healthy human ventricular myocyte. We combined published mathematical models of intracellular signaling and electrophysiology, then calibrated the model results against voltage clamp data and physiological changes occurring after stimulation of β-adrenergic receptors with isoproterenol. We subsequently: (1) explored how molecular variability in different categories of model parameters translated into phenotypic variability; (2) identified the most important parameters determining physiological cellular outputs in the model before and after β-adrenergic receptor stimulation; and (3) investigated which molecular level alterations can produce a phenotype indicative of heart failure with preserved ejection fraction (HFpEF). Major results included: (1) variability in parameters that controlled intracellular signaling caused qualitatively different behavior than variability in parameters controlling ion transport pathways; (2) the most important model parameters determining action potential duration and intracellular Ca²⁺ transient amplitude were generally consistent before and after β-adrenergic receptor stimulation, except for a shift in the importance of K⁺ currents in determining action potential duration; and (3) decreased Ca²⁺ uptake into the sarcoplasmic reticulum, increased Ca²⁺ extrusion through Na⁺/Ca²⁺ exchanger and decreased Ca²⁺ leak from the sarcoplasmic reticulum may contribute to HFpEF. Overall, this study provided novel insight into the phenotypic consequences of molecular variability, and our integrated model may be useful in the design and interpretation of future experimental studies of interactions between β-adrenergic signaling and cellular physiology in human ventricular myocytes.

SUBMITTER: Gong JQX

PROVIDER: S-EPMC7944586 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:The pig is commonly used as an experimental model of human heart disease, including for the study of mechanisms of arrhythmia. However, there exist differences between human and porcine cellular electrophysiology: The pig action potential (AP) has a deeper phase-1 notch, a longer duration at 50% repolarization, and higher plateau potentials than human. Ionic differences underlying the AP include larger rapid delayed-rectifier and smaller inward-rectifier K+-currents (IKr and IK1 respectively) in humans. AP steady-state rate-dependence and restitution is steeper in pigs. Porcine Ca2+ transients can have two components, unlike human. Although a reliable computational model for human ventricular myocytes exists, one for pigs is lacking. This hampers translation from results obtained in pigs to human myocardium. Here, we developed a computational model of the pig ventricular cardiomyocyte AP using experimental datasets of the relevant ionic currents, Ca2+-handling, AP shape, AP duration restitution, and inducibility of triggered activity and alternans. To properly capture porcine Ca2+ transients, we introduced a two-step process with a faster release in the t-tubular region, followed by a slower diffusion-induced release from a non t-tubular subcellular region. The pig model behavior was compared with that of a human ventricular cardiomyocyte (O'Hara-Rudy) model. The pig, but not the human model, developed early afterdepolarizations (EADs) under block of IK1, while IKr block led to EADs in the human but not in the pig model. At fast rates (pacing cycle length = 400 ms), the human cell model was more susceptible to spontaneous Ca2+ release-mediated delayed afterdepolarizations (DADs) and triggered activity than pig. Fast pacing led to alternans in human but not pig. Developing species-specific models incorporating electrophysiology and Ca2+-handling provides a tool to aid translating antiarrhythmic and arrhythmogenic assessment from the bench to the clinic.

Project description:A five-state model of myofilament contraction was integrated into a well-established rabbit ventricular myocyte model of ion channels, Ca(2+) transporters and kinase signaling to analyze the relative contribution of different phosphorylation targets to the overall mechanical response driven by ?-adrenergic stimulation (?-AS). ?-AS effect on sarcoplasmic reticulum Ca(2+) handling, Ca(2+), K(+) and Cl(-) currents, and Na(+)/K(+)-ATPase properties was included based on experimental data. The inotropic effect on the myofilaments was represented as reduced myofilament Ca(2+) sensitivity (XBCa) and titin stiffness, and increased cross-bridge (XB) cycling rate (XBcy). Assuming independent roles of XBCa and XBcy, the model reproduced experimental ?-AS responses on action potentials and Ca(2+) transient amplitude and kinetics. It also replicated the behavior of force-Ca(2+), release-restretch, length-step, stiffness-frequency and force-velocity relationships, and increased force and shortening in isometric and isotonic twitch contractions. The ?-AS effect was then switched off from individual targets to analyze their relative impact on contractility. Preventing ?-AS effects on L-type Ca(2+) channels or phospholamban limited Ca(2+) transients and contractile responses in parallel, while blocking phospholemman and K(+) channel (IKs) effects enhanced Ca(2+) and inotropy. Removal of ?-AS effects from XBCa enhanced contractile force while decreasing peak Ca(2+) (due to greater Ca(2+) buffering), but had less effect on shortening. Conversely, preventing ?-AS effects on XBcy preserved Ca(2+) transient effects, but blunted inotropy (both isometric force and especially shortening). Removal of titin effects had little impact on contraction. Finally, exclusion of ?-AS from XBCa and XBcy while preserving effects on other targets resulted in preserved peak isometric force response (with slower kinetics) but nearly abolished enhanced shortening. ?-AS effects on XBCa and XBcy have greater impact on isometric and isotonic contraction, respectively.

Project description:The ?1-adrenergic signaling system plays an important role in the functioning of cardiac cells. Experimental data shows that the activation of this system produces inotropy, lusitropy, and chronotropy in the heart, such as increased magnitude and relaxation rates of [Ca(2+)]i transients and contraction force, and increased heart rhythm. However, excessive stimulation of ?1-adrenergic receptors leads to heart dysfunction and heart failure. In this paper, a comprehensive, experimentally based mathematical model of the ?1-adrenergic signaling system for mouse ventricular myocytes is developed, which includes major subcellular functional compartments (caveolae, extracaveolae, and cytosol). The model describes biochemical reactions that occur during stimulation of ?1-adrenoceptors, changes in ionic currents, and modifications of Ca(2+) handling system. Simulations describe the dynamics of major signaling molecules, such as cyclic AMP and protein kinase A, in different subcellular compartments; the effects of inhibition of phosphodiesterases on cAMP production; kinetics and magnitudes of phosphorylation of ion channels, transporters, and Ca(2+) handling proteins; modifications of action potential shape and duration; magnitudes and relaxation rates of [Ca(2+)]i transients; changes in intracellular and transmembrane Ca(2+) fluxes; and [Na(+)]i fluxes and dynamics. The model elucidates complex interactions of ionic currents upon activation of ?1-adrenoceptors at different stimulation frequencies, which ultimately lead to a relatively modest increase in action potential duration and significant increase in [Ca(2+)]i transients. In particular, the model includes two subpopulations of the L-type Ca(2+) channels, in caveolae and extracaveolae compartments, and their effects on the action potential and [Ca(2+)]i transients are investigated. The presented model can be used by researchers for the interpretation of experimental data and for the developments of mathematical models for other species or for pathological conditions.

Dataset Information

Quantitative analysis of variability in an integrated model of human ventricular electrophysiology and β-adrenergic signaling.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets