Rate-dependent force, intracellular calcium, and action potential voltage alternans are modulated by sarcomere length and heart failure induced-remodeling of thin filament regulation in human heart failure: A myocyte modeling study

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.

A mathematical model of hiPSC cardiomyocytes electromechanics

Human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) are becoming instrumental in cardiac research, human-based cell level cardiotoxicity tests, and developing patient-specific care. As one of the principal functional readouts is contractility, we propose a novel electromechanical hiPSC-CM computational model named the hiPSC-CM-CE. This model comprises a reparametrized version of contractile element (CE) by Rice et al., 2008, with a new passive force formulation, integrated into a hiPSC-CM electrophysiology formalism by Paci et al. in 2020. Our simulated results were validated against in vitro data reported for hiPSC-CMs at matching conditions from different labs. Specifically, key action potential (AP) and calcium transient (CaT) biomarkers simulated by the hiPSC-CM-CE model were within the experimental ranges. On the mechanical side, simulated cell shortening, contraction-relaxation kinetic indices (RT50 and RT25 ), and the amplitude of tension fell within the experimental intervals. Markedly, as an inter-scale analysis, correct classification of the inotropic effects due to non-cardiomyocytes in hiPSC-CM tissues was predicted on account of the passive force expression introduced to the CE. Finally, the physiological inotropic effects caused by Verapamil and Bay-K 8644 and the aftercontractions due to the early afterdepolarizations (EADs) were simulated and validated against experimental data. In the future, the presented model can be readily expanded to take in pharmacological trials and genetic mutations, such as those involved in hypertrophic cardiomyopathy, and study arrhythmia trigger mechanisms.

Computationally efficient model of myocardial electromechanics for multiscale simulations

A model of myocardial electromechanics is suggested. It combines modified and simplified versions of previously published models of cardiac electrophysiology, excitation-contraction coupling, and mechanics. The mechano-calcium and mechano-electrical feedbacks, including the strain-dependence of the propagation velocity of the action potential, are also accounted for. The model reproduces changes in the twitch amplitude and Ca2+-transients upon changes in muscle strain including the slow response. The model also reproduces the Bowditch effect and changes in the twitch amplitude and duration upon changes in the interstimulus interval, including accelerated relaxation at high stimulation frequency. Special efforts were taken to reduce the stiffness of the differential equations of the model. As a result, the equations can be integrated numerically with a relatively high time step making the model suitable for multiscale simulation of the human heart and allowing one to study the impact of myocardial mechanics on arrhythmias.

Differentiating the effects of β-adrenergic stimulation and stretch on calcium and force dynamics using a novel electromechanical cardiomyocyte model

Cardiac electrophysiology and mechanics are strongly interconnected. Calcium is crucial in this complex interplay through its role in cellular electrophysiology and sarcomere contraction. We aim to differentiate the effects of acute β-adrenergic stimulation (β-ARS) and cardiomyocyte stretch (increased sarcomere length) on calcium-transient dynamics and force generation, using a novel computational model of cardiac electromechanics. We implemented a bidirectional coupling between the O'Hara-Rudy model of human ventricular electrophysiology and the MechChem model of sarcomere mechanics through the buffering of calcium by troponin. The coupled model was validated using experimental data from large mammals or human samples. Calcium transient and force were simulated for various degrees of β-ARS and initial sarcomere lengths. The model reproduced force-frequency, quick-release, and isotonic contraction experiments, validating the bidirectional electromechanical interactions. An increase in β-ARS increased the amplitudes of force (augmented inotropy) and calcium transient, and shortened both force and calcium-transient duration (lusitropy). An increase in sarcomere length increased force amplitude even more, but decreased calcium-transient amplitude and increased both force and calcium-transient duration. Finally, a gradient in relaxation along the thin filament may explain the nonmonotonic decay in cytosolic calcium observed with high tension. Using a novel coupled human electromechanical model, we identified differential effects of β-ARS and stretch on calcium and force. Stretch mostly contributed to increased force amplitude and β-ARS to the reduction of calcium and force duration. We showed that their combination, rather than individual contributions, is key to ensure force generation, rapid relaxation, and low diastolic calcium levels.NEW & NOTEWORTHY This work identifies the contribution of electrical and mechanical alterations to regulation of calcium and force under exercise-like conditions using a novel human electromechanical model integrating ventricular electrophysiology and sarcomere mechanics. By better understanding their individual and combined effects, this can uncover arrhythmogenic mechanisms in exercise-like situations. This publicly available model is a crucial step toward understanding the complex interplay between cardiac electrophysiology and mechanics to improve arrhythmia risk prediction and treatment.

Cardiomyocyte calcium handling in health and disease: Insights from in vitro and in silico studies

Calcium (Ca²⁺) plays a central role in cardiomyocyte excitation-contraction coupling. To ensure an optimal electrical impulse propagation and cardiac contraction, Ca²⁺ levels are regulated by a variety of Ca²⁺-handling proteins. In turn, Ca²⁺ modulates numerous electrophysiological processes. Accordingly, Ca²⁺-handling abnormalities can promote cardiac arrhythmias via various mechanisms, including the promotion of afterdepolarizations, ion-channel modulation and structural remodeling. In the last 30 years, significant improvements have been made in the computational modeling of cardiomyocyte Ca²⁺ handling under physiological and pathological conditions. However, numerous questions involving the Ca²⁺-dependent regulation of different macromolecular complexes, cross-talk between Ca²⁺-dependent regulatory pathways operating over a wide range of time scales, and bidirectional interactions between electrophysiology and mechanics remain to be addressed by in vitro and in silico studies. A better understanding of disease-specific Ca²⁺-dependent proarrhythmic mechanisms may facilitate the development of improved therapeutic strategies. In this review, we describe the fundamental mechanisms of cardiomyocyte Ca²⁺ handling in health and disease, and provide an overview of currently available computational models for cardiomyocyte Ca²⁺ handling. Finally, we discuss important uncertainties and open questions about cardiomyocyte Ca²⁺ handling and highlight how synergy between in vitro and in silico studies may help to answer several of these issues.

Increased thin filament activation enhances alternans in human chronic atrial fibrillation

Action potential duration (APD) alternans (APD-ALT), defined as beat-to-beat oscillations in APD, has been proposed as an important clinical marker for chronic atrial fibrillation (cAF) risk when it occurs at pacing rates of 120-200 beats/min. Although the ionic mechanisms for occurrence of APD-ALT in human cAF at these clinically relevant rates have been investigated, little is known about the effects of myofilament protein kinetics on APD-ALT. Therefore, we used computer simulations of single cell function to explore whether remodeling in myofilament protein kinetics in human cAF alters the occurrence of APD-ALT and to uncover how these mechanisms are affected by sarcomere length and the degree of cAF-induced myofilament remodeling. Mechanistically based, bidirectionally coupled electromechanical models of human right and left atrial myocytes were constructed, incorporating both ionic and myofilament remodeling associated with cAF. By comparing results from our electromechanical model with those from the uncoupled ionic model, we found that intracellular Ca2+ concentration buffering of troponin C has a dampening effect on the magnitude of APD-ALT (APD-ANM) at slower rates (150 beats/min) due to the cooperativity between strongly bound cross-bridges and Ca2+-troponin C binding affinity. We also discovered that cAF-induced enhanced thin filament activation enhanced APD-ANM at these clinically relevant heart rates (150 beats/min). In addition, longer sarcomere lengths increased APD-ANM, suggesting that atrial stretch is an important modulator of APD-ALT. Together, these findings demonstrate that myofilament kinetics mechanisms play an important role in altering APD-ALT in human cAF. NEW & NOTEWORTHY Using a single cell simulation approach, we explored how myofilament protein kinetics alter the formation of alternans in action potential duration (APD) in human myocytes with chronic atrial fibrillation remodeling. We discovered that enhanced thin filament activation and longer sarcomere lengths increased the magnitude of APD alternans at clinically important pacing rates of 120-200 beats/min. Furthermore, we found that altered intracellular Ca2+ concentration buffering of troponin C has a dampening effect on the magnitude of APD alternans.

Quantitative systems models illuminate arrhythmia mechanisms in heart failure: Role of the Na+ -Ca2+ -Ca2+ /calmodulin-dependent protein kinase II-reactive oxygen species feedback

Quantitative systems modeling aims to integrate knowledge in different research areas with models describing biological mechanisms and dynamics to gain a better understanding of complex clinical syndromes. Heart failure (HF) is a chronic complex cardiac disease that results from structural or functional disorders impairing the ability of the ventricle to fill with or eject blood. Highly interactive and dynamic changes in mechanical, structural, neurohumoral, metabolic, and electrophysiological properties collectively predispose the failing heart to cardiac arrhythmias, which are responsible for about a half of HF deaths. Multiscale cardiac modeling and simulation integrate structural and functional data from HF experimental models and patients to improve our mechanistic understanding of this complex arrhythmia syndrome. In particular, they allow investigating how disease-induced remodeling alters the coupling of electrophysiology, Ca²⁺ and Na⁺ handling, contraction, and energetics that lead to rhythm derangements. The Ca²⁺ /calmodulin-dependent protein kinase II, which expression and activity are enhanced in HF, emerges as a critical hub that modulates the feedbacks between these various subsystems and promotes arrhythmogenesis. This article is categorized under: Physiology > Mammalian Physiology in Health and Disease Models of Systems Properties and Processes > Mechanistic Models Models of Systems Properties and Processes > Cellular Models Models of Systems Properties and Processes > Organ, Tissue, and Physiological Models.

An integrative appraisal of mechano-electric feedback mechanisms in the heart

Mechanically-induced alterations in cardiac electrophysiology are referred to as mechano-electric feedback (MEF), and play an important role in electrical regulation of cardiac performance. The influence of mechanical stress and strain on electrophysiology has been investigated at all levels, however the role of MEF in arrhythmia remains poorly understood. During the normal contraction of the heart, mechano-sensitive processes are an implicit component of cardiac activity. Under abnormal mechanical events, stretch-activated mechanisms may contribute to local or global changes in electrophysiology (EP). While such mechanisms have been hypothesised to be involved in mechanically-initiated arrhythmias, the details of these mechanisms and their importance remain elusive. We assess the theoretical role of stretch mechanisms using coupled models of cellular electrophysiology and sarcomere contraction dynamics. Using models of single ventricular myocytes, we first investigated the potential MEF contributions of stretch-activated currents (SAC), and stretch-induced myofilament calcium release, to test how strain and fibrosis may alter cellular electrophysiology. For all models investigated, SACs were alone not sufficient to create a pro-arrhythmic perturbation of the action potential with stretch. However, when combined with stretch-induced myofilament calcium release, the action potential could be shortened depending on the timing of the strain. This effect was highly model dependent, with a canine epicardial EP model being the most sensitive. These model results suggest that known mechanisms of mechano-electric coupling in cardiac myocyte may be sufficient to be pro-arrhythmic, but only in combination and under specific strain patterns.

Dynamics of spatiotemporal line defects and chaos control in complex excitable systems

Spatiotemporal pattern formation governs dynamics and functions in various biological systems. In the heart, excitable waves can form complex oscillatory and chaotic patterns even at an abnormally higher frequency than normal heart beats, which increase the risk of fatal heart conditions by inhibiting normal blood circulation. Previous studies suggested that line defects (nodal lines) play a critical role in stabilizing those undesirable patterns. However, it remains unknown if the line defects are static or dynamically changing structures in heart tissue. Through in vitro experiments of heart tissue observation, we reveal the spatiotemporal dynamics of line defects in rotating spiral waves. We combined a novel signaling over-sampling technique with a multi-dimensional Fourier analysis, showing that line defects can translate, merge, collapse and form stable singularities with even and odd parity while maintaining a stable oscillation of the spiral wave in the tissue. These findings provide insights into a broad class of complex periodic systems, with particular impact to the control and understanding of heart diseases.

Myofilament protein dynamics modulate EAD formation in human hypertrophic cardiomyopathy

Patients with hypertrophic cardiomyopathy (HCM), a disease associated with sarcomeric protein mutations, often suffer from sudden cardiac death (SCD) resulting from arrhythmia. In order to advance SCD prevention strategies, our understanding of how sarcomeric mutations in HCM patients contribute to enhanced arrhythmogenesis needs to be improved. Early afterdepolarizations (EADs) are an important mechanism underlying arrhythmias associated with HCM-SCD. Although the ionic mechanisms underlying EADs have been studied in general, whether myofilament protein dynamics mechanisms also underlie EADs remains unknown. Thus, our goals were to investigate if myofilament protein dynamics mechanisms underlie EADs and to uncover how those mechanisms are affected by pacing rate, sarcomere length (SL), and different levels of HCM-induced myofilament remodeling. To achieve this, a mechanistically-based bidirectionally coupled human electrophysiology-force myocyte model under the conditions of HCM was constructed. HCM ionic remodeling included a reduced repolarization reserve, while HCM myofilament modeling involved altered thin filament activation. We found that the mechanoelectric feedback (MEF) on calcium dynamics in the bidirectionally coupled model, via Troponin C buffering of cytoplasmic Ca²⁺, was the myofilament mechanism underlying EADs. Incorporating MEF diminished the degree of repolarization reserve reduction necessary for EADs to emerge and increased the frequency of EAD occurrence, especially at faster pacing rates. Longer SLs and enhanced thin filament activation diminished the effects of MEF on EADs. Together these findings demonstrate that myofilament protein dynamics mechanisms play an important role in EAD formation.

The link between abnormal calcium handling and electrical instability in acquired long QT syndrome--Does calcium precipitate arrhythmic storms?

Release of Ca(2+) ions from sarcoplasmic reticulum (SR) into myocyte cytoplasm and their binding to troponin C is the final signal form myocardial contraction. Synchronous contraction of ventricular myocytes is necessary for efficient cardiac pumping function. This requires both shuttling of Ca(2+) between SR and cytoplasm in individual myocytes, and organ-level synchronization of this process by means of electrical coupling among ventricular myocytes. Abnormal Ca(2+) release from SR causes arrhythmias in the setting of CPVT (catecholaminergic polymorphic ventricular tachycardia) and digoxin toxicity. Recent optical mapping data indicate that abnormal Ca(2+) handling causes arrhythmias in models of both repolarization impairment and profound bradycardia. The mechanisms involve dynamic spatial heterogeneity of myocardial Ca(2+) handling preceding arrhythmia onset, cell-synchronous systolic secondary Ca(2+) elevation (SSCE), as well as more complex abnormalities of intracellular Ca(2+) handling detected by subcellular optical mapping in Langendorff-perfused hearts. The regional heterogeneities in Ca(2+) handling cause action potential (AP) heterogeneities through sodium-calcium exchange (NCX) activation and eventually overwhelm electrical coupling of the tissue. Divergent Ca(2+) dynamics among different myocardial regions leads to temporal instability of AP duration and - on the patient level - in T wave lability. Although T-wave alternans has been linked to cardiac arrhythmias, non-alternans lability is observed in pre-clinical models of the long QT syndrome (LQTS) and CPVT, and in LQTS patients. Analysis of T wave lability may provide a real-time window on the abnormal Ca(2+) dynamics causing specific arrhythmias such as Torsade de Pointes (TdP).