Transmural heterogeneity of action potentials and Ito1 in myocytes isolated from the human right ventricle

Injection of IK1 through dynamic clamp can make all the difference in patch-clamp studies on hiPSC-derived cardiomyocytes

Human-induced stem cell-derived cardiomyocytes (hiPSC-CMs) are a valuable tool for studying development, pharmacology, and (inherited) arrhythmias. Unfortunately, hiPSC-CMs are depolarized and spontaneously active, even the working cardiomyocyte subtypes such as atrial- and ventricular-like hiPSC-CMs, in contrast to the situation in the atria and ventricles of adult human hearts. Great efforts have been made, using many different strategies, to generate more mature, quiescent hiPSC-CMs with more close-to-physiological resting membrane potentials, but despite promising results, it is still difficult to obtain hiPSC-CMs with such properties. The dynamic clamp technique allows to inject a current with characteristics of the inward rectifier potassium current (IK1), computed in real time according to the actual membrane potential, into patch-clamped hiPSC-CMs during action potential measurements. This results in quiescent hiPSC-CMs with a close-to-physiological resting membrane potential. As a result, action potential measurements can be performed with normal ion channel availability, which is particularly important for the physiological functioning of the cardiac SCN5A-encoded fast sodium current (INa). We performed in vitro and in silico experiments to assess the beneficial effects of the dynamic clamp technique in dissecting the functional consequences of the SCN5A-1795insD+/- mutation. In two separate sets of patch-clamp experiments on control hiPSC-CMs and on hiPSC-CMs with mutations in ACADVL and GNB5, we assessed the value of dynamic clamp in detecting delayed afterdepolarizations and in investigating factors that modulate the resting membrane potential. We conclude that the dynamic clamp technique has highly beneficial effects in all of the aforementioned settings and should be widely used in patch-clamp studies on hiPSC-CMs while waiting for the ultimate fully mature hiPSC-CMs.

A transmurally heterogeneous model of the ventricular tissue and its application for simulation of Brugada Syndrome

We present a transmurally heterogeneous phe-nomenological model of ventricular tissue that is designed to reproduce the most important features of action potential prop-agation of endocardial, midmyocardial, and epicardial tissue. Our model consists of only 3 variables and 20 parameters. Therefore, it is highly computational efficient and easy to fit to experimental data. We exploited our myocyte model to simulate action potential propagation in a 3D slab of cardiac tissue both in healthy conditions and in presence of Brugada syndrome. The results show that our model can accurately reproduce the transmural heterogeneity of the ventricular wall and the main characteristics of electrocardiographic pattern both in healthy and pathological conditions.

A computationally efficient dynamic model of human epicardial tissue

We present a new phenomenological model of human ventricular epicardial cells and we test its reentry dynamics. The model is derived from the Rogers-McCulloch formulation of the FitzHugh-Nagumo equations and represents the total ionic current divided into three contributions corresponding to the excitatory, recovery and transient outward currents. Our model reproduces the main characteristics of human epicardial tissue, including action potential amplitude and morphology, upstroke velocity, and action potential duration and conduction velocity restitution curves. The reentry dynamics is stable, and the dominant period is about 270 ms, which is comparable to clinical values. The proposed model is the first phenomenological model able to accurately resemble human experimental data by using only 3 state variables and 17 parameters. Indeed, it is more computationally efficient than existing models (i.e., almost two times faster than the minimal ventricular model). Beyond the computational efficiency, the low number of parameters facilitates the process of fitting the model to the experimental data.

Sex Differences in Drug-Induced Arrhythmogenesis

The electrical activity in the heart varies significantly between men and women and results in a sex-specific response to drugs. Recent evidence suggests that women are more than twice as likely as men to develop drug-induced arrhythmia with potentially fatal consequences. Yet, the sex-specific differences in drug-induced arrhythmogenesis remain poorly understood. Here we integrate multiscale modeling and machine learning to gain mechanistic insight into the sex-specific origin of drug-induced cardiac arrhythmia at differing drug concentrations. To quantify critical drug concentrations in male and female hearts, we identify the most important ion channels that trigger male and female arrhythmogenesis, and create and train a sex-specific multi-fidelity arrhythmogenic risk classifier. Our study reveals that sex differences in ion channel activity, tissue conductivity, and heart dimensions trigger longer QT-intervals in women than in men. We quantify the critical drug concentration for dofetilide, a high risk drug, to be seven times lower for women than for men. Our results emphasize the importance of including sex as an independent biological variable in risk assessment during drug development. Acknowledging and understanding sex differences in drug safety evaluation is critical when developing novel therapeutic treatments on a personalized basis. The general trends of this study have significant implications on the development of safe and efficacious new drugs and the prescription of existing drugs in combination with other drugs.

Role of the rabbit whole-heart model for electrophysiologic safety pharmacology of non-cardiovascular drugs

Plenty of non-cardiovascular drugs alter cardiac electrophysiology and may ultimately lead to life-threatening arrhythmias. In clinical practice, measuring the QT interval as a marker for the repolarization period is the most common tool to assess the electrophysiologic safety of drugs. However, the sole measurement of the QT interval may be insufficient to determine the proarrhythmic risk of non-cardiovascular agents. Several other markers are considered in pre-clinical safety testing to determine potential harm on cardiac electrophysiology. Besides measuring typical electrophysiologic parameters such as repolarization duration, whole-heart models allow the determination of potential predictors for proarrhythmia. Spatial and temporal heterogeneity as well as changes of shape of the action potential can be easily assessed. In addition, provocation manoeuvers (either by electrolyte imbalances or programmed pacing protocols) may induce sustained arrhythmias and thereby determine ventricular vulnerability to arrhythmias. Compared with the human heart, the rabbit heart possesses a similar distribution of ion currents that govern cardiac repolarization, resulting in a rectangular action potential configuration in both species. In addition, similar biophysical properties of rabbit and human cardiac ion channels lead to a comparable pharmacologic response in human and rabbit hearts. Of note, arrhythmia patterns resemble in both species due to the similar effective size of human and rabbit hearts. Thus, the rabbit heart is particularly suitable for testing the electrophysiologic safety of drugs. Several experimental setups have been developed for studying cardiac electrophysiology in rabbits, ranging from single cell to tissue preparations, whole-heart setups, and in vivo models.

Canine Myocytes Represent a Good Model for Human Ventricular Cells Regarding Their Electrophysiological Properties

Due to the limited availability of healthy human ventricular tissues, the most suitable animal model has to be applied for electrophysiological and pharmacological studies. This can be best identified by studying the properties of ion currents shaping the action potential in the frequently used laboratory animals, such as dogs, rabbits, guinea pigs, or rats, and comparing them to those of human cardiomyocytes. The authors of this article with the experience of three decades of electrophysiological studies, performed in mammalian and human ventricular tissues and isolated cardiomyocytes, summarize their results obtained regarding the major canine and human cardiac ion currents. Accordingly, L-type Ca²⁺ current (I_Ca), late Na⁺ current (I_Na-late), rapid and slow components of the delayed rectifier K⁺ current (I_Kr and I_Ks, respectively), inward rectifier K⁺ current (I_K1), transient outward K⁺ current (I_to1), and Na⁺/Ca²⁺ exchange current (I_NCX) were characterized and compared. Importantly, many of these measurements were performed using the action potential voltage clamp technique allowing for visualization of the actual current profiles flowing during the ventricular action potential. Densities and shapes of these ion currents, as well as the action potential configuration, were similar in human and canine ventricular cells, except for the density of I_K1 and the recovery kinetics of I_to. I_K1 displayed a largely four-fold larger density in canine than human myocytes, and I_to recovery from inactivation displayed a somewhat different time course in the two species. On the basis of these results, it is concluded that canine ventricular cells represent a reasonably good model for human myocytes for electrophysiological studies, however, it must be borne in mind that due to their stronger I_K1, the repolarization reserve is more pronounced in canine cells, and moderate differences in the frequency-dependent repolarization patterns can also be anticipated.

A Framework for the generation of digital twins of cardiac electrophysiology from clinical 12-leads ECGs

Cardiac digital twins (Cardiac Digital Twin (CDT)s) of human electrophysiology (Electrophysiology (EP)) are digital replicas of patient hearts derived from clinical data that match like-for-like all available clinical observations. Due to their inherent predictive potential, CDTs show high promise as a complementary modality aiding in clinical decision making and also in the cost-effective, safe and ethical testing of novel EP device therapies. However, current workflows for both the anatomical and functional twinning phases within CDT generation, referring to the inference of model anatomy and parameters from clinical data, are not sufficiently efficient, robust and accurate for advanced clinical and industrial applications. Our study addresses three primary limitations impeding the routine generation of high-fidelity CDTs by introducing; a comprehensive parameter vector encapsulating all factors relating to the ventricular EP; an abstract reference frame within the model allowing the unattended manipulation of model parameter fields; a novel fast-forward electrocardiogram (Electrocardiogram (ECG)) model for efficient and bio-physically-detailed simulation required for parameter inference. A novel workflow for the generation of CDTs is then introduced as an initial proof of concept. Anatomical twinning was performed within a reasonable time compatible with clinical workflows (<4h) for 12 subjects from clinically-attained magnetic resonance images. After assessment of the underlying fast forward ECG model against a gold standard bidomain ECG model, functional twinning of optimal parameters according to a clinically-attained 12 lead ECG was then performed using a forward Saltelli sampling approach for a single subject. The achieved results in terms of efficiency and fidelity demonstrate that our workflow is well-suited and viable for generating biophysically-detailed CDTs at scale.

Towards Mutation-Specific Precision Medicine in Atypical Clinical Phenotypes of Inherited Arrhythmia Syndromes

Most causal genes for inherited arrhythmia syndromes (IASs) encode cardiac ion channel-related proteins. Genotype-phenotype studies and functional analyses of mutant genes, using heterologous expression systems and animal models, have revealed the pathophysiology of IASs and enabled, in part, the establishment of causal gene-specific precision medicine. Additionally, the utilization of induced pluripotent stem cell (iPSC) technology have provided further insights into the pathophysiology of IASs and novel promising therapeutic strategies, especially in long QT syndrome. It is now known that there are atypical clinical phenotypes of IASs associated with specific mutations that have unique electrophysiological properties, which raises a possibility of mutation-specific precision medicine. In particular, patients with Brugada syndrome harboring an SCN5A R1632C mutation exhibit exercise-induced cardiac events, which may be caused by a marked activity-dependent loss of R1632C-Nav1.5 availability due to a marked delay of recovery from inactivation. This suggests that the use of isoproterenol should be avoided. Conversely, the efficacy of β-blocker needs to be examined. Patients harboring a KCND3 V392I mutation exhibit both cardiac (early repolarization syndrome and paroxysmal atrial fibrillation) and cerebral (epilepsy) phenotypes, which may be associated with a unique mixed electrophysiological property of V392I-Kv4.3. Since the epileptic phenotype appears to manifest prior to cardiac events in this mutation carrier, identifying KCND3 mutations in patients with epilepsy and providing optimal therapy will help prevent sudden unexpected death in epilepsy. Further studies using the iPSC technology may provide novel insights into the pathophysiology of atypical clinical phenotypes of IASs and the development of mutation-specific precision medicine.

Characterization of Spatio-Temporal Cardiac Action Potential Variability at Baseline and Under β-Adrenergic Stimulation by Combined Unscented Kalman Filter and Double Greedy Dimension Reduction

OBJECTIVE

Elevated spatio-temporal variability of human ventricular repolarization has been related to increased risk for ventricular arrhythmias and sudden cardiac death, particularly under β-adrenergic stimulation ( β-AS). This work presents a methodology for theoretical characterization of temporal and spatial repolarization variability at baseline conditions and in response to β-AS. For any measured voltage trace, the proposed methodology estimates the parameters and state variables of an underlying human ventricular action potential (AP) model by combining Double Greedy Dimension Reduction (DGDR) with automatic selection of biomarkers and the Unscented Kalman Filter (UKF). Such theoretical characterization can facilitate subsequent characterization of underlying variability mechanisms.

MATERIAL AND METHODS

Given an AP trace, initial estimates for the ionic conductances in a stochastic version of the baseline human ventricular O'Hara et al. model were obtained by DGDR. Those estimates served to initialize and update model parameter estimates by the UKF method based on formulation of an associated nonlinear state-space representation and joint estimation of model parameters and state variables. Similarly, β-AS-induced phosphorylation levels of cellular substrates were estimated by the DGDR-UKF methodology. Performance was tested by building an experimentally-calibrated population of virtual cells, from which synthetic AP traces were generated for baseline and β-AS conditions.

RESULTS

The combined DGDR-UKF methodology led to 25% reduction in the error associated with estimation of ionic current conductances at baseline conditions and phosphorylation levels under β-AS with respect to individual DGDR and UKF methods. This improvement was not at the expense of higher computational load, which was diminished by 90% with respect to the individual UKF method. Both temporal and spatial AP variability of repolarization were accurately characterized by the DGDR-UKF methodology.

CONCLUSIONS

A combined DGDR-UKF methodology is proposed for parameter and state variable estimation of human ventricular cell models from available AP traces at baseline and under β-AS. This methodology improves the estimation performance and reduces the convergence time with respect to individual DGDR and UKF methods and renders a suitable approach for computational characterization of spatio-temporal repolarization variability to be used for ascertainment of variability mechanisms and its relation to arrhythmogenesis.

Dispersion of ventricular repolarization: Temporal and spatial

Repolarization heterogeneity (RH) is an intrinsic property of ventricular myocardium and the reason for T-wave formation on electrocardiogram (ECG). Exceeding the physiologically based RH level is associated with appearance of life-threatening ventricular arrhythmias and sudden cardiac death. In this regard, an accurate and comprehensive evaluation of the degree of RH parameters is of importance for assessment of heart state and arrhythmic risk. This review is devoted to comprehensive consideration of RH phenomena in terms of electrophysiological processes underlying RH, cardiac electric field formation during ventricular repolarization, as well as clinical significance of RH and its reflection on ECG parameters. The formation of transmural, apicobasal, left-to-right and anterior-posterior gradients of action potential durations and end of repolarization times resulting from the heterogenous distribution of repolarizing ion currents and action potential morphology throughout the heart ventricles, and the different sensitivity of myocardial cells in different ventricular regions to the action of pharmacological agents, temperature, frequency of stimulation, etc., are being discussed. The review is focused on the fact that RH has different aspects – temporal and spatial, global and local; ECG reflection of various RH aspects and their clinical significance are being discussed. Strategies for comprehensive assessment of ventricular RH using different ECG indices reflecting various RH aspects are presented.

Simulating Notch-Dome Morphology of Action Potential of Ventricular Cell: How the Speeds of Positive and Negative Feedbacks on Transmembrane Voltage Can Influence the Health of a Cell?

Ventricular action potential is well-known because of its plateau phase with a spike-notch-dome morphology. As such, the morphology of action potential is necessary for ensuring a correct heart functioning. Any distraction from normal notch-dome morphology may trigger a circus movement reentry in the form of lethal ventricular fibrillation. When the epicardial action potential dome propagates from a site where it is maintained to regions where it has been lost, it gives rise to the proposed mechanism for the Brugada syndrome. Despite the impact of notch-dome dynamics on the heart function, no independent and explicit research has been performed on the simulation of notch-dome dynamics and morphology. In this paper, using a novel mathematical approach, a three-state variable model is proposed; we show that our proposed model not only can simulate morphology of action potential of ventricular cells but also can propose a biological reasonable tool for controlling of the morphology of action potential spike-notch-dome. We show that the processes of activation and inactivation of ionic gating variables (as positive or negative feedbacks on the voltage of cell membrane) and the ratio of their speeds (time constants) can be treated as a reasonable biological tool for simulating ventricular cell notch-dome. This finding may led to a new insight to the quantification of the health of a ventricular cell and may also propose a new drug therapy strategy for cardiac diseases.

In-silico human electro-mechanical ventricular modelling and simulation for drug-induced pro-arrhythmia and inotropic risk assessment

Human-based computational modelling and simulation are powerful tools to accelerate the mechanistic understanding of cardiac patho-physiology, and to develop and evaluate therapeutic interventions. The aim of this study is to calibrate and evaluate human ventricular electro-mechanical models for investigations on the effect of the electro-mechanical coupling and pharmacological action on human ventricular electrophysiology, calcium dynamics, and active contraction. The most recent models of human ventricular electrophysiology, excitation-contraction coupling, and active contraction were integrated, and the coupled models were calibrated using human experimental data. Simulations were then conducted using the coupled models to quantify the effects of electro-mechanical coupling and drug exposure on electrophysiology and force generation in virtual human ventricular cardiomyocytes and tissue. The resulting calibrated human electro-mechanical models yielded active tension, action potential, and calcium transient metrics that are in agreement with experiments for endocardial, epicardial, and mid-myocardial human samples. Simulation results correctly predicted the inotropic response of different multichannel action reference compounds and demonstrated that the electro-mechanical coupling improves the robustness of repolarisation under drug exposure compared to electrophysiology-only models. They also generated additional evidence to explain the partial mismatch between in-silico and in-vitro experiments on drug-induced electrophysiology changes. The human calibrated and evaluated modelling and simulation framework constructed in this study opens new avenues for future investigations into the complex interplay between the electrical and mechanical cardiac substrates, its modulation by pharmacological action, and its translation to tissue and organ models of cardiac patho-physiology.

A model for human action potential dynamics in vivo

Computational modeling based on experimental data remains an important component in cardiac electrophysiological research, especially because clinical data such as human action potential (AP) dynamics are scarce or limited by practical or ethical concerns. Such modeling has been used to develop and test a variety of mechanistic hypotheses, with the majority of these studies involving the rate dependence of AP duration (APD) including APD restitution and conduction velocity (CV). However, there is very little information regarding the complex dynamics at the boundary of repolarization (or refractoriness) and reexcitability. Here, we developed a "minimal" ionic model of the human AP, based on in vivo human monophasic AP (MAP) recordings obtained during clinical programmed electrical stimulation (PES) to address the progressive decrease in AP take-off potential (TOP) and associated CV slowing seen during three tightly spaced extrastimuli. Recent voltage-clamp data demonstrating the effect of intracellular calcium on sodium current availability were incorporated and were required to reproduce large (>15 mV) elevations in take-off potential and progressive encroachment. Introducing clinically observed APD gradients into the model enabled us to replicate the dynamic response to PES in patients leading to conduction block and reentry formation for the positive, but not the negative, APD gradient. Finally, we modeled the dynamics of reentry and show that spiral waves follow a meandering trajectory with a period of ~180 ms. We conclude that our model reproduces a variety of electrophysiological behavior including the response to sequential premature stimuli and provides a basis for studies of the initiation of reentry in human ventricular tissue.NEW & NOTEWORTHY This work presents a new model of the action potential of the human which reproduces the complex dynamics during premature stimulation in patients.

Biowire Model of Interstitial and Focal Cardiac Fibrosis

Myocardial fibrosis is a severe global health problem due to its prevalence in all forms of cardiac diseases and direct role in causing heart failure. The discovery of efficient antifibrotic compounds has been hampered due to the lack of a physiologically relevant disease model. Herein, we present a disease model of human myocardial fibrosis and use it to establish a compound screening system. In the Biowire II platform, cardiac tissues are suspended between a pair of poly(octamethylene maleate (anhydride) citrate) (POMaC) wires. Noninvasive functional readouts are realized on the basis of the deflection of the intrinsically fluorescent polymer. The disease model is constructed to recapitulate contractile, biomechanical, and electrophysiological complexities of fibrotic myocardium. Additionally, we constructed a heteropolar integrated model with fibrotic and healthy cardiac tissues coupled together. The integrated model captures the regional heterogeneity of scar lesion, border zone, and adjacent healthy myocardium. Finally, we demonstrate the utility of the system for the evaluation of antifibrotic compounds. The high-fidelity in vitro model system combined with convenient functional readouts could potentially facilitate the development of precision medicine strategies for cardiac fibrosis modeling and establish a pipeline for preclinical compound screening.

Data-Driven Identification of Stochastic Model Parameters and State Variables: Application to the Study of Cardiac Beat-to-Beat Variability

OBJECTIVE

Enhanced spatiotemporal ventricular repolarization variability has been associated with ventricular arrhythmias and sudden cardiac death, but the involved mechanisms remain elusive. In this paper, a methodology for estimation of parameters and state variables of stochastic human ventricular cell models from input voltage data is proposed for investigation of repolarization variability.

METHODS

The proposed methodology formulates state-space representations based on developed stochastic cell models and uses the unscented Kalman filter to perform joint parameter and state estimation. Evaluation over synthetic and experimental data is presented.

RESULTS

Results on synthetically generated data show the ability of the methodology to: first, filter out measurement noise from action potential (AP) traces; second, identify model parameters and state variables from each of those individual AP traces, thus allowing robust characterization of cell-to-cell variability; and, third, replicate statistical population's distributions of input AP-based markers, including dynamic markers quantifying beat-to-beat variability. Application onto experimental data demonstrates the ability of the methodology to match input AP traces while concomitantly inferring the characteristics of underlying stochastic cell models.

CONCLUSION

A novel methodology is presented for estimation of parameters and hidden variables of stochastic cardiac computational models, with the advantage of providing a one-to-one match between each individual AP trace and a corresponding set of model characteristics.

SIGNIFICANCE

The proposed methodology can greatly help in the characterization of temporal (beat-to-beat) and spatial (cell-to-cell) variability in human ventricular repolarization and in ascertaining the corresponding underlying mechanisms, particularly in scenarios with limited available experimental data.

Comparative study of carvedilol and quinidine for inhibiting hKv4.3 channel stably expressed in HEK 293 cells

The inhibition of transient outward potassium current (I_to) is the major ionic mechanism for quinidine to treat Brugada syndrome; however, quinidine is inaccessible in many countries. The present study compared the inhibitory effect of the nonselective β-adrenergic blocker carvedilol with quinidine on human Kv4.3 (hKv4.3, encoding for I_to) channel and action potential notch using a whole-cell patch technique in HEK 293 cell line expressing KCND3 as well as in ventricular epicardial myocytes of rabbit hearts. It was found that carvedilol and quinidine inhibited hKv4.3 current in a concentration-dependent manner. The IC₅₀ of carvedilol was 1.2 μM for inhibiting hKv4.3 charge area, while the IC₅₀ of quinidine was 2.9 μM (0.2 Hz). Both carvedilol and quinidine showed typical open channel blocking properties (i.e. decreasing the time to peak of activation and increasing the inactivation of hKv4.3), negatively shifted the V_1/2 of activation and inactivation, and slowed the recovery from inactivation of the channel. Although carvedilol had weaker in use- and rate-dependent inhibition of hKv4.3 peak current than quinidine, its reduction of the charge area was more than quinidine at all frequencies (0.2-3.3 Hz). Moreover, the inhibitory effect of carvedilol on action potential notch was greater than quinidine. These results provide the novel information that carvedilol, like quinidine, significantly inhibits hKv4.3 and action potential notch, suggesting that carvedilol is likely an alternative drug for preventing malignant ventricular arrhythmias in patients with Brugada syndrome in countries where quinidine is unavailable.

A response surface optimization approach to adjust ionic current conductances of cardiac electrophysiological models. Application to the study of potassium level changes

Cardiac electrophysiological computational models are often developed from previously published models. The new models may incorporate additional features to adapt the model to a different species or may upgrade a specific ionic formulation based on newly available experimental data. A relevant challenge in the development of a new model is the estimation of certain ionic current conductances that cannot be reliably identified from experiments. A common strategy to estimate those conductances is by means of constrained non-linear least-squares optimization. In this work, a novel methodology is proposed for estimation of ionic current conductances of cardiac electrophysiological models by using a response surface approximation-based constrained optimization with trust region management. Polynomial response surfaces of a number of electrophysiological markers were built using statistical sampling methods. These markers included action potential duration (APD), triangulation, diastolic and systolic intracellular calcium concentration, and time constants of APD rate adaptation. The proposed methodology was applied to update the Carro et al. human ventricular action potential model after incorporation of intracellular potassium ([K⁺]_i) dynamics. While the Carro et al. model was well suited for investigation of arrhythmogenesis, it did not allow simulation of [K⁺]_i changes. With the methodology proposed in this study, the updated Carro et al. human ventricular model could be used to simulate [K⁺]_i changes in response to varying extracellular potassium ([K⁺]_o) levels. Additionally, it rendered values of evaluated electrophysiological markers within physiologically plausible ranges. The optimal values of ionic current conductances in the updated model were found in a notably shorter time than with previously proposed methodologies. As a conclusion, the response surface optimization-based approach proposed in this study allows estimating ionic current conductances of cardiac electrophysiological computational models while guaranteeing replication of key electrophysiological features and with an important reduction in computational cost with respect to previously published approaches. The updated Carro et al. model developed in this study is thus suitable for the investigation of arrhythmic risk-related conditions, including those involving large changes in potassium concentration.

Ca2+ Cycling Impairment in Heart Failure Is Exacerbated by Fibrosis: Insights Gained From Mechanistic Simulations

Heart failure (HF) is characterized by altered Ca²⁺ cycling, resulting in cardiac contractile dysfunction. Failing myocytes undergo electrophysiological remodeling, which is known to be the main cause of abnormal Ca²⁺ homeostasis. However, structural remodeling, specifically proliferating fibroblasts coupled to myocytes in the failing heart, could also contribute to Ca²⁺ cycling impairment. The goal of the present study was to systematically analyze the mechanisms by which myocyte–fibroblast coupling could affect Ca²⁺ dynamics in normal conditions and in HF. Simulations of healthy and failing human myocytes were performed using established mathematical models, and cells were either isolated or coupled to fibroblasts. Univariate and multivariate sensitivity analyses were performed to quantify effects of ion transport pathways on biomarkers computed from intracellular [Ca²⁺] waveforms. Variability in ion channels and pumps was imposed and populations of models were analyzed to determine effects on Ca²⁺ dynamics. Our results suggest that both univariate and multivariate sensitivity analyses are valuable methodologies to shed light into the ionic mechanisms underlying Ca²⁺ impairment in HF, although differences between the two methodologies are observed at high parameter variability. These can result from either the fact that multivariate analyses take into account ion channels or non-linear effects of ion transport pathways on Ca²⁺ dynamics. Coupling either healthy or failing myocytes to fibroblasts decreased Ca²⁺ transients due to an indirect sink effect on action potential (AP) and thus on Ca²⁺ related currents. Simulations that investigated restoration of normal physiology in failing myocytes showed that Ca²⁺ cycling can be normalized by increasing SERCA and L-type Ca²⁺ current activity while decreasing Na⁺–Ca²⁺ exchange and SR Ca²⁺ leak. Changes required to normalize APs in failing myocytes depended on whether myocytes were coupled to fibroblasts. In conclusion, univariate and multivariate sensitivity analyses are helpful tools to understand how Ca²⁺ cycling is impaired in HF and how this can be exacerbated by coupling of myocytes to fibroblasts. The design of pharmacological actions to restore normal activity should take into account the degree of fibrosis in the failing heart.

Mechanisms Underlying the Actions of Antidepressant and Antipsychotic Drugs That Cause Sudden Cardiac Arrest

A number of antipsychotic and antidepressant drugs are known to increase the risk of ventricular arrhythmias and sudden cardiac death. Based largely on a concern over the development of life-threatening arrhythmias, a number of antipsychotic drugs have been temporarily or permanently withdrawn from the market or their use restricted. While many antidepressants and antipsychotics have been linked to QT prolongation and the development of torsade de pointes arrhythmias, some have been associated with a Brugada syndrome phenotype and the development of polymorphic ventricular arrhythmias. This article examines the arrhythmic liability of antipsychotic and antidepressant drugs capable of inducing long QT and/or Brugada syndrome phenotypes. The goal of this article is to provide an update on the ionic and cellular mechanisms thought to be involved in, and the genetic and environmental factors that predispose to, the development of cardiac arrhythmias and sudden cardiac death among patients taking antidepressant and antipsychotic drugs that are in clinical use.

A Multiwell Cardiac μGMEA Platform for Action Potential Recordings from Human iPSC-Derived Cardiomyocyte Constructs

Multielectrode array (MEA) technology has been extensively used for field potential recordings from excitable cells. However, its application for action potential (AP) measurements has not been harnessed. Here, we report a novel platform for high-resolution intracellular AP recordings from induced pluripotent stem cell-cardiomyocyte constructs derived from human cardiac fibroblasts. To gain intracellular access, micro-gold MEAs were used to electroporate multiple constructs simultaneously. High-throughput AP measurements were obtained from 41 multicellular constructs. Repeated electroporations of the same cells did not affect the signal stability. Our model has the capability to distinguish subtle differences in AP morphology to characterize the network profile. Furthermore, we confirm the reliability of the system by recapitulating known drug-induced physiological and arrhythmogenic responses. Overall, the model provides a unique cardio-electronic interface for non-invasive measurements of AP dynamics for drug screening and disease modeling. This technology opens the door for identifying novel cardio-factors to enhance electrophysiological maturation.

•

Electroporation-mediated action potential (AP) recordings using MEA technology

•

Simultaneous high-throughput AP measurement from multiple cell networks

•

Multiple electroporations of the same cells over days with stable signal

•

Model validation for developmental, disease, and drug screening studies

Effect of Transmural Differences in Excitation-Contraction Delay and Contraction Velocity on Left Ventricle Isovolumic Contraction: A Simulation Study

Recent studies have shown that left ventricle (LV) exhibits considerable transmural differences in active mechanical properties induced by transmural differences in electrical activity, excitation-contraction coupling, and contractile properties of individual myocytes. It was shown that the time between electrical and mechanical activation of myocytes (electromechanical delay: EMD) decreases from subendocardium to subepicardium and, on the contrary, the myocyte shortening velocity (MSV) increases in the same direction. To investigate the physiological importance of this inhomogeneity, we developed a new finite element model of LV incorporating the observed transmural gradients in EMD and MSV. Comparative simulations with the model showed that when EMD or MSV or both were set constant across the LV wall, the LV contractility during isovolumic contraction (IVC) decreased significantly ((dp/dt)max⁡  was reduced by 2 to 38% and IVC was prolonged by 18 to 73%). This was accompanied by an increase of transmural differences in wall stress. These results suggest that the transmural differences in EMD and MSV play an important role in physiological contractility of LV by synchronising the contraction of individual layers of ventricular wall during the systole. Reduction or enhancement of these differences may therefore impair the function of LV and contribute to heart failure.

Low Resting Membrane Potential and Low Inward Rectifier Potassium Currents Are Not Inherent Features of hiPSC-Derived Cardiomyocytes

Human induced pluripotent stem cell (hiPSC) cardiomyocytes (CMs) show less negative resting membrane potential (RMP), which is attributed to small inward rectifier currents (IK1). Here, IK1 was measured in hiPSC-CMs (proprietary and commercial cell line) cultured as monolayer (ML) or 3D engineered heart tissue (EHT) and, for direct comparison, in CMs from human right atrial (RA) and left ventricular (LV) tissue. RMP was measured in isolated cells and intact tissues. IK1 density in ML- and EHT-CMs from the proprietary line was similar to LV and RA, respectively. IK1 density in EHT-CMs from the commercial line was 2-fold smaller than in the proprietary line. RMP in EHT of both lines was similar to RA and LV. Repolarization fraction and IK,ACh response discriminated best between RA and LV and indicated predominantly ventricular phenotype in hiPSC-CMs/EHT. The data indicate that IK1 is not necessarily low in hiPSC-CMs, and technical issues may underlie low RMP in hiPSC-CMs.

Patch-Clamp Recording from Human Induced Pluripotent Stem Cell-Derived Cardiomyocytes: Improving Action Potential Characteristics through Dynamic Clamp

Human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) hold great promise for studying inherited cardiac arrhythmias and developing drug therapies to treat such arrhythmias. Unfortunately, until now, action potential (AP) measurements in hiPSC-CMs have been hampered by the virtual absence of the inward rectifier potassium current (I_K1) in hiPSC-CMs, resulting in spontaneous activity and altered function of various depolarising and repolarising membrane currents. We assessed whether AP measurements in "ventricular-like" and "atrial-like" hiPSC-CMs could be improved through a simple, highly reproducible dynamic clamp approach to provide these cells with a substantial I_K1 (computed in real time according to the actual membrane potential and injected through the patch-clamp pipette). APs were measured at 1 Hz using perforated patch-clamp methodology, both in control cells and in cells treated with all-trans retinoic acid (RA) during the differentiation process to increase the number of cells with atrial-like APs. RA-treated hiPSC-CMs displayed shorter APs than control hiPSC-CMs and this phenotype became more prominent upon addition of synthetic I_K1 through dynamic clamp. Furthermore, the variability of several AP parameters decreased upon I_K1 injection. Computer simulations with models of ventricular-like and atrial-like hiPSC-CMs demonstrated the importance of selecting an appropriate synthetic I_K1. In conclusion, the dynamic clamp-based approach of I_K1 injection has broad applicability for detailed AP measurements in hiPSC-CMs.

Modeling hypothermia induced effects for the heterogeneous ventricular tissue from cellular level to the impact on the ECG

Hypothermia has a profound impact on the electrophysiological mechanisms of the heart. Experimental investigations provide a better understanding of electrophysiological alterations associated with cooling. However, there is a lack of computer models suitable for simulating the effects of hypothermia in cardio-electrophysiology. In this work, we propose a model that describes the cooling-induced electrophysiological alterations in ventricular tissue in a temperature range from 27°C to 37°C. To model the electrophysiological conditions in a 3D left ventricular tissue block it was essential to consider the following anatomical and physiological parameters in the model: the different cell types (endocardial, M, epicardial), the heterogeneous conductivities in longitudinal, transversal and transmural direction depending on the prevailing temperature, the distinct fiber orientations and the transmural repolarization sequences. Cooling-induced alterations on the morphology of the action potential (AP) of single myocardial cells thereby are described by an extension of the selected Bueno-Orovio model for human ventricular tissue using Q10 temperature coefficients. To evaluate alterations on tissue level, the corresponding pseudo electrocardiogram (pECG) was calculated. Simulations show that cooling-induced AP and pECG-related parameters, i.e. AP duration, morphology of the notch of epicardial AP, maximum AP upstroke velocity, AP rise time, QT interval, QRS duration and J wave formation are in good accordance with literature and our experimental data. The proposed model enables us to further enhance our knowledge of cooling-induced electrophysiological alterations from cellular to tissue level in the heart and may help to better understand electrophysiological mechanisms, e.g. in arrhythmias, during hypothermia.

Effect of eicosapentaenoic acid and pitavastatin on electrophysiology and anticoagulant gene expression in mice with rapid atrial pacing

Atrial remodeling is considered to be any persistent change in atrial structure or function, and is responsible for the development and perpetuation of atrial fibrillation (AF). Oxidative stress and intracellular pH regulation may also be linked to AF; however it remains unclear whether eicosapentaenoic acid (EPA) or statins have beneficial therapeutic effects. The aim of the present study was to investigate the effects of EPA and pitavastatin on the electrophysiology of and gene expressions in mice with rapidly-paced atria. Mice were treated with EPA (10 mg/g/day) or pitavastatin (30 ng/g/day) for 6 weeks, following which AF was simulated by 8-h atrial pacing at 1,800 bpm. The atrial electrophysiological properties and the expression of cardiac genes, potassium voltage-gated channel subfamily A member 5 (Kcna5), Kcn subfamily D member 2 (Kcnd2), Kv channel-interacting protein 2 (KChIP2), solute carrier family 9 member A1, thrombomodulin (TM) and tissue factor pathway inhibitor (TFPI) were examined using reverse transcription-quantitative polymerase chain reaction. In control mice, significant atrial electrical remodeling was observed (P<0.05); however, treatment with either EPA or pitavastatin ameliorated these electrophysiological changes (P>0.05). mRNA levels of Kcnd2, KChIP2 and Kcna5 were significantly upregulated in control mice (P<0.05), whereas treatment with EPA or pitavastatin attenuated this upregulation (P>0.05). Administration of pitavastatin significantly reduced the downregulation of both TFPI and TM (P<0.05). EPA treatment attenuated the TFPI downregulation compared with control mice (P>0.05), however no significant effect on TM expression was observed. In addition, both EPA (P>0.05) and pitavastatin (P<0.05) suppressed the overexpression of endothelial nitric oxide synthase. This was also exhibited in Ras-related C3 botulinum toxin substrate 1 genes (P<0.01 for both treatments). In conclusion, the results of the present study suggested that EPA and pitavastatin are able to prevent atrial electrical remodeling, thrombotic states and oxidative stress in rapidly-paced murine atria.

Role of abnormal repolarization in the mechanism of cardiac arrhythmia

In cardiac patients, life-threatening tachyarrhythmia is often precipitated by abnormal changes in ventricular repolarization and refractoriness. Repolarization abnormalities typically evolve as a consequence of impaired function of outward K⁺ currents in cardiac myocytes, which may be caused by genetic defects or result from various acquired pathophysiological conditions, including electrical remodelling in cardiac disease, ion channel modulation by clinically used pharmacological agents, and systemic electrolyte disorders seen in heart failure, such as hypokalaemia. Cardiac electrical instability attributed to abnormal repolarization relies on the complex interplay between a provocative arrhythmic trigger and vulnerable arrhythmic substrate, with a central role played by the excessive prolongation of ventricular action potential duration, impaired intracellular Ca²⁺ handling, and slowed impulse conduction. This review outlines the electrical activity of ventricular myocytes in normal conditions and cardiac disease, describes classical electrophysiological mechanisms of cardiac arrhythmia, and provides an update on repolarization-related surrogates currently used to assess arrhythmic propensity, including spatial dispersion of repolarization, activation-repolarization coupling, electrical restitution, TRIaD (triangulation, reverse use dependence, instability, and dispersion), and the electromechanical window. This is followed by a discussion of the mechanisms that account for the dependence of arrhythmic vulnerability on the location of the ventricular pacing site. Finally, the review clarifies the electrophysiological basis for cardiac arrhythmia produced by hypokalaemia, and gives insight into the clinical importance and pathophysiology of drug-induced arrhythmia, with particular focus on class Ia (quinidine, procainamide) and Ic (flecainide) Na⁺ channel blockers, and class III antiarrhythmic agents that block the delayed rectifier K⁺ channel (dofetilide).

Concise Review: Criteria for Chamber-Specific Categorization of Human Cardiac Myocytes Derived from Pluripotent Stem Cells

Human pluripotent stem cell‐derived cardiomyocytes (PSC‐CMs) have great potential application in almost all areas of cardiovascular research. A current major goal of the field is to build on the past success of differentiation strategies to produce CMs with the properties of those originating from the different chambers of the adult human heart. With no anatomical origin or developmental pathway to draw on, the question of how to judge the success of such approaches and assess the chamber specificity of PSC‐CMs has become increasingly important; commonly used methods have substantial limitations and are based on limited evidence to form such an assessment. In this article, we discuss the need for chamber‐specific PSC‐CMs in a number of areas as well as current approaches used to assess these cells on their likeness to those from different chambers of the heart. Furthermore, describing in detail the structural and functional features that distinguish the different chamber‐specific human adult cardiac myocytes, we propose an evidence‐based tool to aid investigators in the phenotypic characterization of differentiated PSC‐CMs. Stem Cells2017;35:1881–1897

A mathematical model for active contraction in healthy and failing myocytes and left ventricles

Cardiovascular disease is one of the leading causes of death worldwide, in particular myocardial dysfunction, which may lead to heart failure eventually. Understanding the electro-mechanics of the heart will help in developing more effective clinical treatments. In this paper, we present a multi-scale electro-mechanics model of the left ventricle (LV). The Holzapfel-Ogden constitutive law was used to describe the passive myocardial response in tissue level, a modified Grandi-Pasqualini-Bers model was adopted to model calcium dynamics in individual myocytes, and the active tension was described using the Niederer-Hunter-Smith myofilament model. We first studied the electro-mechanics coupling in a single myocyte in the healthy and diseased left ventricle, and then the single cell model was embedded in a dynamic LV model to investigate the compensation mechanism of LV pump function due to myocardial dysfunction caused by abnormality in cellular calcium dynamics. The multi-scale LV model was solved using an in-house developed hybrid immersed boundary method with finite element extension. The predictions of the healthy LV model agreed well with the clinical measurements and other studies, and likewise, the results in the failing states were also consistent with clinical observations. In particular, we found that a low level of intracellular Ca2+ transient in myocytes can result in LV pump function failure even with increased myocardial contractility, decreased systolic blood pressure, and increased diastolic filling pressure, even though they will increase LV stroke volume. Our work suggested that treatments targeted at increased contractility and lowering the systolic blood pressure alone are not sufficient in preventing LV pump dysfunction, restoring a balanced physiological Ca2+ handling mechanism is necessary.

Transmural gradients in ion channel and auxiliary subunit expression

Evolution has acted to shape the action potential in different regions of the heart in order to produce a maximally stable and efficient pump. This has been achieved by creating regional differences in ion channel expression levels within the heart as well as differences between equivalent cardiac tissues in different species. These region- and species-dependent differences in channel expression are established by regulatory evolution, evolution of the regulatory mechanisms that control channel expression levels. Ion channel auxiliary subunits are obvious targets for regulatory evolution, in order to change channel expression levels and/or modify channel function. This review focuses on the transmural gradients of ion channel expression in the heart and the role that regulation of auxiliary subunit expression plays in generating and shaping these gradients.

Ventricular stimulus site influences dynamic dispersion of repolarization in the intact human heart

Spatial variation of restitution in relation to varying stimulus site is poorly defined in the intact human heart. Repolarization gradients were shown to be dependent on site of activation with epicardial stimulation promoting significant transmural gradients. Steep restitution slopes were predominant in the normal ventricle.

The spatial variation in restitution properties in relation to varying stimulus site is poorly defined. This study aimed to investigate the effect of varying stimulus site on apicobasal and transmural activation time (AT), action potential duration (APD) and repolarization time (RT) during restitution studies in the intact human heart. Ten patients with structurally normal hearts, undergoing clinical electrophysiology studies, were enrolled. Decapolar catheters were placed apex to base in the endocardial right ventricle (RV_endo) and left ventricle (LV_endo), and an LV branch of the coronary sinus (LV_epi) for transmural recording. S1–S2 restitution protocols were performed pacing RV_endo apex, LV_endo base, and LV_epi base. Overall, 725 restitution curves were analyzed, 74% of slopes had a maximum slope of activation recovery interval (ARI) restitution (S_max) > 1 (P < 0.001); mean S_max = 1.76. APD was shorter in the LV_epi compared with LV_endo, regardless of pacing site (30-ms difference during RV_endo pacing, 25-ms during LV_endo, and 48-ms during LV_epi; 50th quantile, P < 0.01). Basal LV_epi pacing resulted in a significant transmural gradient of RT (77 ms, 50th quantile: P < 0.01), due to loss of negative transmural AT-APD coupling (mean slope 0.63 ± 0.3). No significant transmural gradient in RT was demonstrated during endocardial RV or LV pacing, with preserved negative transmural AT-APD coupling (mean slope −1.36 ± 1.9 and −0.71 ± 0.4, respectively). Steep ARI restitution slopes predominate in the normal ventricle and dynamic ARI; RT gradients exist that are modulated by the site of activation. Epicardial stimulation to initiate ventricular activation promotes significant transmural gradients of repolarization that could be proarrhythmic.

A human ventricular myocyte model with a refined representation of excitation-contraction coupling

Cardiac Ca(2+)-induced Ca(2+) release (CICR) occurs by a regenerative activation of ryanodine receptors (RyRs) within each Ca(2+)-releasing unit, triggered by the activation of L-type Ca(2+) channels (LCCs). CICR is then terminated, most probably by depletion of Ca(2+) in the junctional sarcoplasmic reticulum (SR). Hinch et al. previously developed a tightly coupled LCC-RyR mathematical model, known as the Hinch model, that enables simulations to deal with a variety of functional states of whole-cell populations of a Ca(2+)-releasing unit using a personal computer. In this study, we developed a membrane excitation-contraction model of the human ventricular myocyte, which we call the human ventricular cell (HuVEC) model. This model is a hybrid of the most recent HuVEC models and the Hinch model. We modified the Hinch model to reproduce the regenerative activation and termination of CICR. In particular, we removed the inactivated RyR state and separated the single step of RyR activation by LCCs into triggering and regenerative steps. More importantly, we included the experimental measurement of a transient rise in Ca(2+) concentrations ([Ca(2+)], 10-15 μM) during CICR in the vicinity of Ca(2+)-releasing sites, and thereby calculated the effects of the local Ca(2+) gradient on CICR as well as membrane excitation. This HuVEC model successfully reconstructed both membrane excitation and key properties of CICR. The time course of CICR evoked by an action potential was accounted for by autonomous changes in an instantaneous equilibrium open probability of couplons. This autonomous time course was driven by a core feedback loop including the pivotal local [Ca(2+)], influenced by a time-dependent decay in the SR Ca(2+) content during CICR.

IK1-enhanced human-induced pluripotent stem cell-derived cardiomyocytes: an improved cardiomyocyte model to investigate inherited arrhythmia syndromes

Currently available induced pluripotent stem cell-derived cardiomyocytes (iPS-CMs) do not ideally model cellular mechanisms of human arrhythmic disease due to lack of a mature action potential (AP) phenotype. In this study, we create and characterize iPS-CMs with an electrically mature AP induced by potassium inward rectifier (IK1) enhancement. The advantages of IK1-enhanced iPS-CMs include the absence of spontaneous beating, stable resting membrane potentials at approximately -80 mV and capability for electrical pacing. Compared with unenhanced, IK1-enhanced iPS-CMs calcium transient amplitudes were larger (P < 0.05) with a typical staircase pattern. IK1-enhanced iPS-CMs demonstrated a twofold increase in cell size and membrane capacitance and increased DNA synthesis compared with control iPS-CMs (P < 0.05). Furthermore, IK1-enhanced iPS-CMs expressing the F97C-CAV3 long QT9 mutation compared with wild-type CAV3 demonstrated an increase in AP duration and late sodium current. IK1-enhanced iPS-CMs represent a more mature cardiomyocyte model to study arrhythmia mechanisms.

Human pluripotent stem cells: Prospects and challenges as a source of cardiomyocytes for in vitro modeling and cell-based cardiac repair

Human pluripotent stem cells (PSCs) represent an attractive source of cardiomyocytes with potential applications including disease modeling, drug discovery and safety screening, and novel cell-based cardiac therapies. Insights from embryology have contributed to the development of efficient, reliable methods capable of generating large quantities of human PSC-cardiomyocytes with cardiac purities ranging up to 90%. However, for human PSCs to meet their full potential, the field must identify methods to generate cardiomyocyte populations that are uniform in subtype (e.g. homogeneous ventricular cardiomyocytes) and have more mature structural and functional properties. For in vivo applications, cardiomyocyte production must be highly scalable and clinical grade, and we will need to overcome challenges including graft cell death, immune rejection, arrhythmogenesis, and tumorigenic potential. Here we discuss the types of human PSCs, commonly used methods to guide their differentiation into cardiomyocytes, the phenotype of the resultant cardiomyocytes, and the remaining obstacles to their successful translation.

Human induced pluripotent stem cell-derived versus adult cardiomyocytes: an in silico electrophysiological study on effects of ionic current block

Background and Purpose

Two new technologies are likely to revolutionize cardiac safety and drug development: in vitro experiments on human‐induced pluripotent stem cell‐derived cardiomyocytes (hiPSC‐CMs) and in silico human adult ventricular cardiomyocyte (hAdultV‐CM) models. Their combination was recently proposed as a potential replacement for the present hERG‐based QT study for pharmacological safety assessments. Here, we systematically compared in silico the effects of selective ionic current block on hiPSC‐CM and hAdultV‐CM action potentials (APs), to identify similarities/differences and to illustrate the potential of computational models as supportive tools for evaluating new in vitro technologies.

Experimental Approach

In silico AP models of ventricular‐like and atrial‐like hiPSC‐CMs and hAdultV‐CM were used to simulate the main effects of four degrees of block of the main cardiac transmembrane currents.

Key Results

Qualitatively, hiPSC‐CM and hAdultV‐CM APs showed similar responses to current block, consistent with results from experiments. However, quantitatively, hiPSC‐CMs were more sensitive to block of (i) L‐type Ca²⁺ currents due to the overexpression of the Na⁺/Ca²⁺ exchanger (leading to shorter APs) and (ii) the inward rectifier K⁺ current due to reduced repolarization reserve (inducing diastolic potential depolarization and repolarization failure).

Conclusions and Implications

In silico hiPSC‐CMs and hAdultV‐CMs exhibit a similar response to selective current blocks. However, overall hiPSC‐CMs show greater sensitivity to block, which may facilitate in vitro identification of drug‐induced effects. Extrapolation of drug effects from hiPSC‐CM to hAdultV‐CM and pro‐arrhythmic risk assessment can be facilitated by in silico predictions using biophysically‐based computational models.

Ion Channels in the Heart

Optimal cardiac function depends on proper timing of excitation and contraction in various regions of the heart, as well as on appropriate heart rate. This is accomplished via specialized electrical properties of various components of the system, including the sinoatrial node, atria, atrioventricular node, His-Purkinje system, and ventricles. Here we review the major regionally determined electrical properties of these cardiac regions and present the available data regarding the molecular and ionic bases of regional cardiac function and dysfunction. Understanding these differences is of fundamental importance for the investigation of arrhythmia mechanisms and pharmacotherapy.

In silico assessment of drug safety in human heart applied to late sodium current blockers

Drug-induced action potential (AP) prolongation leading to Torsade de Pointes is a major concern for the development of anti-arrhythmic drugs. Nevertheless the development of improved anti-arrhythmic agents, some of which may block different channels, remains an important opportunity. Partial block of the late sodium current (I(NaL)) has emerged as a novel anti-arrhythmic mechanism. It can be effective in the settings of free radical challenge or hypoxia. In addition, this approach can attenuate pro-arrhythmic effects of blocking the rapid delayed rectifying K(+) current (I(Kr)). The main goal of our computational work was to develop an in-silico tool for preclinical anti-arrhythmic drug safety assessment, by illustrating the impact of I(Kr)/I(NaL) ratio of steady-state block of drug candidates on "torsadogenic" biomarkers. The O'Hara et al. AP model for human ventricular myocytes was used. Biomarkers for arrhythmic risk, i.e., AP duration, triangulation, reverse rate-dependence, transmural dispersion of repolarization and electrocardiogram QT intervals, were calculated using single myocyte and one-dimensional strand simulations. Predetermined amounts of block of I(NaL) and I(Kr) were evaluated. "Safety plots" were developed to illustrate the value of the specific biomarker for selected combinations of IC(50)s for I(Kr) and I(NaL) of potential drugs. The reference biomarkers at baseline changed depending on the "drug" specificity for these two ion channel targets. Ranolazine and GS967 (a novel potent inhibitor of I(NaL)) yielded a biomarker data set that is considered safe by standard regulatory criteria. This novel in-silico approach is useful for evaluating pro-arrhythmic potential of drugs and drug candidates in the human ventricle.

Abnormal calcium homeostasis in heart failure with preserved ejection fraction is related to both reduced contractile function and incomplete relaxation: an electromechanically detailed biophysical modeling study

Heart failure with preserved ejection fraction (HFpEF) accounts for about 50% of heart failure cases. It has features of incomplete relaxation and increased stiffness of the left ventricle. Studies from clinical electrophysiology and animal experiments have found that HFpEF is associated with impaired calcium homeostasis, ion channel remodeling and concentric left ventricle hypertrophy (LVH). However, it is still unclear how the abnormal calcium homeostasis, ion channel and structural remodeling affect the electro-mechanical dynamics of the ventricles. In this study we have developed multiscale models of the human left ventricle from single cells to the 3D organ, which take into consideration HFpEF-induced changes in calcium handling, ion channel remodeling and concentric LVH. Our simulation results suggest that at the cellular level, HFpEF reduces the systolic calcium level resulting in a reduced systolic contractile force, but elevates the diastolic calcium level resulting in an abnormal residual diastolic force. In our simulations, these abnormal electro-mechanical features of the ventricular cells became more pronounced with the increase of the heart rate. However, at the 3D organ level, the ejection fraction of the left ventricle was maintained due to the concentric LVH. The simulation results of this study mirror clinically observed features of HFpEF and provide new insights toward the understanding of the cellular bases of impaired cardiac electromechanical functions in heart failure.

Ion channelopathies in human induced pluripotent stem cell derived cardiomyocytes: a dynamic clamp study with virtual IK1

Human induced pluripotent stem cell derived cardiomyocytes (hiPSC-CMs) are widely used in studying basic mechanisms of cardiac arrhythmias that are caused by ion channelopathies. Unfortunately, the action potential profile of hiPSC-CMs-and consequently the profile of individual membrane currents active during that action potential-differs substantially from that of native human cardiomyocytes, largely due to almost negligible expression of the inward rectifier potassium current (IK1). In the present study, we attempted to "normalize" the action potential profile of our hiPSC-CMs by inserting a voltage dependent in silico IK1 into our hiPSC-CMs, using the dynamic clamp configuration of the patch clamp technique. Recordings were made from single hiPSC-CMs, using the perforated patch clamp technique at physiological temperature. We assessed three different models of IK1, with different degrees of inward rectification, and systematically varied the magnitude of the inserted IK1. Also, we modified the inserted IK1 in order to assess the effects of loss- and gain-of-function mutations in the KCNJ2 gene, which encodes the Kir2.1 protein that is primarily responsible for the IK1 channel in human ventricle. For our experiments, we selected spontaneously beating hiPSC-CMs, with negligible IK1 as demonstrated in separate voltage clamp experiments, which were paced at 1 Hz. Upon addition of in silico IK1 with a peak outward density of 4-6 pA/pF, these hiPSC-CMs showed a ventricular-like action potential morphology with a stable resting membrane potential near -80 mV and a maximum upstroke velocity >150 V/s (n = 9). Proarrhythmic action potential changes were observed upon injection of both loss-of-function and gain-of-function IK1, as associated with Andersen-Tawil syndrome type 1 and short QT syndrome type 3, respectively (n = 6). We conclude that injection of in silico IK1 makes the hiPSC-CM a more reliable model for investigating mechanisms underlying cardiac arrhythmias.

The potassium current carried by TREK-1 channels in rat cardiac ventricular muscle

We studied the potassium current flowing through TREK-1 channels in rat cardiac ventricular myocytes. We separated the TREK-1 current from other current components by blocking most other channels with a blocker cocktail. We tried to inhibit the TREK-1 current by activating protein kinase A (PKA) with a mixture of forskolin and isobutyl-methylxanthine (IBMX). Activation of PKA blocked an outwardly rectifying current component at membrane potentials positive to -40 mV. At 37 °C, application of forskolin plus IBMX reduced the steady-state outward current measured at positive voltages by about 52 %. Application of the potassium channel blockers quinidine or tetrahexylammonium also reduced the steady-state outward current by about 50 %. Taken together, our results suggest that the increase in temperature from 22 to 37 °C increased the TREK-1 current by a factor of at least 5 and that the average density of the TREK-1 current in rat cardiomyocytes at 37 °C is about 1.5 pA/pF at +30 mV. The contribution of TREK-1 to the action potential was assessed by using a dynamic patch clamp technique. After subtraction of simulated TREK-1 currents, action potential duration at 50 or 90 % repolarisation was increased by about 12 %, indicating that TREK-1 may be functionally important in rat ventricular muscle. During sympathetic stimulation, inhibition of TREK-1 channels via PKA is expected to prolong the action potential primarily in subendocardial myocytes; this may decrease the transmural dispersion of repolarisation and thus may serve to prevent the occurrence of arrhythmias.

Dependence of phase-2 reentry and repolarization dispersion on epicardial and transmural ionic heterogeneity: a simulation study

AIMS

Phase-2 reentry (P2R) is a local arrhythmogenic phenomenon where electrotonic current propagates from a spike-and-dome action potential region to re-excite a loss-of-dome action potential region. While ionic heterogeneity has been shown to underlie P2R within the epicardium and has been hypothesized to occur transmurally, we are unaware of any study that has investigated the effects of combining these heterogeneities as they occur in the heart. Thus, we tested the hypothesis that P2R can result by either epicardial or transmural heterogeneity and that the realistic combination of the two would increase the likelihood of P2R.

METHODS AND RESULTS

We used computational ionic models of cardiac myocyte dynamics to investigate initiation and development of P2R in simulated tissues with different ionic heterogeneities. In one-dimensional transmural cable simulations, P2R occurred when the conductance of the transient outward current in the epicardial region was near the range for which epicardial action potentials switched intermittently between spike-and-dome and loss-of-dome morphologies. Phase-2 reentry was more likely in two-dimensional tissue simulations by both epicardial and transmural heterogeneity and could expand beyond its local initiation site to create a macroscopic reentry.

CONCLUSION

The characteristics and stability of action potential morphology in the epicardium are important determinants of the occurrence of both transmural and epicardial P2R and its associated arrhythmogenesis.

Cardiac potassium channel subtypes: new roles in repolarization and arrhythmia

About 10 distinct potassium channels in the heart are involved in shaping the action potential. Some of the K+ channels are primarily responsible for early repolarization, whereas others drive late repolarization and still others are open throughout the cardiac cycle. Three main K+ channels drive the late repolarization of the ventricle with some redundancy, and in atria this repolarization reserve is supplemented by the fairly atrial-specific KV1.5, Kir3, KCa, and K2P channels. The role of the latter two subtypes in atria is currently being clarified, and several findings indicate that they could constitute targets for new pharmacological treatment of atrial fibrillation. The interplay between the different K+ channel subtypes in both atria and ventricle is dynamic, and a significant up- and downregulation occurs in disease states such as atrial fibrillation or heart failure. The underlying posttranscriptional and posttranslational remodeling of the individual K+ channels changes their activity and significance relative to each other, and they must be viewed together to understand their role in keeping a stable heart rhythm, also under menacing conditions like attacks of reentry arrhythmia.

Comprehensive analyses of ventricular myocyte models identify targets exhibiting favorable rate dependence

Reverse rate dependence is a problematic property of antiarrhythmic drugs that prolong the cardiac action potential (AP). The prolongation caused by reverse rate dependent agents is greater at slow heart rates, resulting in both reduced arrhythmia suppression at fast rates and increased arrhythmia risk at slow rates. The opposite property, forward rate dependence, would theoretically overcome these parallel problems, yet forward rate dependent (FRD) antiarrhythmics remain elusive. Moreover, there is evidence that reverse rate dependence is an intrinsic property of perturbations to the AP. We have addressed the possibility of forward rate dependence by performing a comprehensive analysis of 13 ventricular myocyte models. By simulating populations of myocytes with varying properties and analyzing population results statistically, we simultaneously predicted the rate-dependent effects of changes in multiple model parameters. An average of 40 parameters were tested in each model, and effects on AP duration were assessed at slow (0.2 Hz) and fast (2 Hz) rates. The analysis identified a variety of FRD ionic current perturbations and generated specific predictions regarding their mechanisms. For instance, an increase in L-type calcium current is FRD when this is accompanied by indirect, rate-dependent changes in slow delayed rectifier potassium current. A comparison of predictions across models identified inward rectifier potassium current and the sodium-potassium pump as the two targets most likely to produce FRD AP prolongation. Finally, a statistical analysis of results from the 13 models demonstrated that models displaying minimal rate-dependent changes in AP shape have little capacity for FRD perturbations, whereas models with large shape changes have considerable FRD potential. This can explain differences between species and between ventricular cell types. Overall, this study provides new insights, both specific and general, into the determinants of AP duration rate dependence, and illustrates a strategy for the design of potentially beneficial antiarrhythmic drugs.

Several drugs intended to treat cardiac arrhythmias have failed because of unfavorable rate-dependent properties. That is, the drugs fail to alter electrical activity at fast heart rates, where this would be beneficial, but they do affect electrical activity at slow rates, where this is unwanted. In targeted studies, several agents have been shown to exhibit these unfavorable properties, suggesting that these rate-dependent responses may be intrinsic to ventricular muscle. To determine whether drugs with desirable rate-dependent properties could be rationally designed, we performed comprehensive and systematic analyses of several heart cell models. These analyses calculated the rate-dependent properties of changes in any model parameter, thereby generating simultaneously a large number of model predictions. The analyses showed that targets with favorable rate-dependent properties could indeed be identified, and further simulations uncovered the mechanisms underlying these behaviors. Moreover, a quantitative comparison of results obtained in different models provided new insight in why a given drug applied to different species, or to different tissue types, might produce different rate-dependent behaviors. Overall this study shows how a comprehensive and systematic approach to heart cell models can both identify novel targets and produce more general insight into rate-dependent alterations to cardiac electrical activity.

Overexpression of myocardin induces partial transdifferentiation of human-induced pluripotent stem cell-derived mesenchymal stem cells into cardiomyocytes

Mesenchymal stem cells (MSCs) derived from human‐induced pluripotent stem cells (iPSCs) show superior proliferative capacity and therapeutic potential than those derived from bone marrow (BM). Ectopic expression of myocardin further improved the therapeutic potential of BM‐MSCs in a mouse model of myocardial infarction. The aim was of this study was to assess whether forced myocardin expression in iPSC‐MSCs could further enhance their transdifferentiation to cardiomyocytes and improve their electrophysiological properties for cardiac regeneration. Myocardin was overexpressed in iPSC‐MSCs using viral vectors (adenovirus or lentivirus). The expression of smooth muscle cell and cardiomyocyte markers, and ion channel genes was examined by reverse transcription‐polymerase chain reaction (RT‐PCR), immunofluorescence staining and patch clamp. The conduction velocity of the neonatal rat ventricular cardiomyocytes cocultured with iPSC‐MSC monolayer was measured by multielectrode arrays recording plate. Myocardin induced the expression of α‐MHC, GATA4, α‐actinin, cardiac MHC, MYH11, calponin, and SM α‐actin, but not cTnT, β‐MHC, and MLC2v in iPSC‐MSCs. Overexpression of myocardin in iPSC‐MSC enhanced the expression of SCN9A and CACNA1C, but reduced that of KCa3.1 and Kir2.2 in iPSC‐MSCs. Moreover, BK_Ca, I_Kir, I_Cl, I_to and I_Na.TTX were detected in iPSC‐MSC with myocardin overexpression; while only BK_Ca, I_Kir, I_Cl, IK_DR, and IK_Ca were noted in iPSC‐MSC transfected with green florescence protein. Furthermore, the conduction velocity of iPSC‐MSC was significantly increased after myocardin overexpression. Overexpression of myocardin in iPSC‐MSCs resulted in partial transdifferentiation into cardiomyocytes phenotype and improved the electrical conduction during integration with mature cardiomyocytes.

Forced myocardin expression in human‐induced pluripotent stem cell (hiPSC)‐derived mesenchymal stem cells lead to partial transdifferentiation into cardiomyocytes and smooth muscle cells phenotypes through modification in ion channel expression profile and electrical conduction velocity.

A quantitative comparison of the behavior of human ventricular cardiac electrophysiology models in tissue

Numerical integration of mathematical models of heart cell electrophysiology provides an important computational tool for studying cardiac arrhythmias, but the abundance of available models complicates selecting an appropriate model. We study the behavior of two recently published models of human ventricular action potentials, the Grandi-Pasqualini-Bers (GPB) and the O'Hara-Virág-Varró-Rudy (OVVR) models, and compare the results with four previously published models and with available experimental and clinical data. We find the shapes and durations of action potentials and calcium transients differ between the GPB and OVVR models, as do the magnitudes and rate-dependent properties of transmembrane currents and the calcium transient. Differences also occur in the steady-state and S1–S2 action potential duration and conduction velocity restitution curves, including a maximum conduction velocity for the OVVR model roughly half that of the GPB model and well below clinical values. Between single cells and tissue, both models exhibit differences in properties, including maximum upstroke velocity, action potential amplitude, and minimum diastolic interval. Compared to experimental data, action potential durations for the GPB and OVVR models agree fairly well (although OVVR epicardial action potentials are shorter), but maximum slopes of steady-state restitution curves are smaller. Although studies show alternans in normal hearts, it occurs only in the OVVR model, and only for a narrow range of cycle lengths. We find initiated spiral waves do not progress to sustained breakup for either model. The dominant spiral wave period of the GPB model falls within clinically relevant values for ventricular tachycardia (VT), but for the OVVR model, the dominant period is longer than periods associated with VT. Our results should facilitate choosing a model to match properties of interest in human cardiac tissue and to replicate arrhythmia behavior more closely. Furthermore, by indicating areas where existing models disagree, our findings suggest avenues for further experimental work.

Effect of ion concentration changes in the limited extracellular spaces on sarcolemmal ion transport and Ca2+ turnover in a model of human ventricular cardiomyocyte

We have developed a computer model of human cardiac ventricular myocyte (CVM), including t-tubular and cleft spaces with the aim of evaluating the impact of accumulation-depletion of ions in restricted extracellular spaces on transmembrane ion transport and ionic homeostasis in human CVM. The model was based on available data from human CVMs. Under steady state, the effect of ion concentration changes in extracellular spaces on [Ca²⁺]_i-transient was explored as a function of critical fractions of ion transporters in t-tubular membrane (not documented for human CVM). Depletion of Ca²⁺ and accumulation of K⁺ occurring in extracellular spaces slightly affected the transmembrane Ca²⁺ flux, but not the action potential duration (APD₉₀). The [Ca²⁺]_i-transient was reduced (by 2%–9%), depending on the stimulation frequency, the rate of ion exchange between t-tubules and clefts and fractions of ion-transfer proteins in the t-tubular membrane. Under non-steady state, the responses of the model to changes of stimulation frequency were analyzed. A sudden increase of frequency (1–2.5 Hz) caused a temporal decrease of [Ca²⁺] in both extracellular spaces, a reduction of [Ca²⁺]_i-transient (by 15%) and APD₉₀ (by 13 ms). The results reveal different effects of activity-related ion concentration changes in human cardiac t-tubules (steady-state effects) and intercellular clefts (transient effects) in the modulation of membrane ion transport and Ca²⁺ turnover.

Effects of deformation on transmural dispersion of repolarization using in silico models of human left ventricular wedge

Mechanical deformation affects the electrical activity of the heart through multiple feedback loops. The purpose of this work is to study the effect of deformation on transmural dispersion of repolarization and on surface electrograms using an in silico human ventricular wedge. To achieve this purpose, we developed a strongly coupled electromechanical cell model by coupling a human left ventricle electrophysiology model and an active contraction model reparameterized for human cells. This model was then embedded in tissue simulations on the basis of bidomain equations and nonlinear solid mechanics. The coupled model was used to evaluate effects of mechanical deformation on important features of repolarization and electrograms. Our results indicate an increase in the T-wave amplitude of the surface electrograms in simulations that account for the effects of cardiac deformation. This increased T-wave amplitude can be explained by changes to the coupling between neighboring myocytes, also known as electrotonic effect. The thickening of the ventricular wall during repolarization contributes to the decoupling of cells in the transmural direction, enhancing action potential heterogeneity and increasing both transmural repolarization dispersion and T-wave amplitude of surface electrograms. The simulations suggest that a considerable percentage of the T-wave amplitude (15%) may be related to cardiac deformation.

Sodium overload due to a persistent current that attenuates the arrhythmogenic potential of a novel LQT3 mutation

Long QT syndrome (LQTS) is a congenital abnormality of cardiac repolarization that manifests as a prolonged QT interval on 12-lead electrocardiograms (ECGs). The syndrome may lead to syncope and sudden death from ventricular tachyarrhythmias known as torsades de pointes. An increased persistent Na⁺ current is known to cause a Ca²⁺ overload in case of ischemia for example. Such increased Na⁺ persistent current is also usually associated to the LQT3 syndrome. The purpose of this study was to investigate the pathological consequences of a novel mutation in a family affected by LQTS. The impact of biophysical defects on cellular homeostasis are also investigated. Genomic DNA was extracted from blood samples, and a combination of PCR and DNA sequencing of several LQTS-linked genes was used to identify mutations. The mutation was reproduced in vitro and was characterized using the patch clamp technique and in silico quantitative analysis. A novel mutation (Q1476R) was identified on the SCN5A gene encoding the cardiac Na⁺ channel. Cells expressing the Q1476R mutation exhibited biophysical alterations, including a shift of SS inactivation and a significant increase in the persistent Na⁺ current. The in silico analysis confirmed the arrhythmogenic character of the Q1476R mutation. It further revealed that the increase in persistent Na⁺ current causes a frequency-dependent Na⁺ overload in cardiomyocytes co-expressing WT and mutant Na_v1.5 channels that, in turn, exerts a moderating effect on the lengthening of the action potential (AP) duration caused by the mutation. The Q1476R mutation in SCN5A results in a three-fold increase in the window current and a persistent inward Na⁺ current. These biophysical defects may expose the carrier of the mutation to arrhythmias that occur preferentially in the patient at rest or during tachycardia. However, the Na⁺ overload counterbalances the gain-of-function of the mutation and is beneficial in that it prevents severe arrhythmias at intermediate heart rates.

Electronic "expression" of the inward rectifier in cardiocytes derived from human-induced pluripotent stem cells

BACKGROUND

Human-induced pluripotent stem cell (h-iPSC)-derived cardiac myocytes are a unique model in which human myocyte function and dysfunction are studied, especially those from patients with genetic disorders. They are also considered a major advance for drug safety testing. However, these cells have considerable unexplored potential limitations when applied to quantitative action potential (AP) analysis. One major factor is spontaneous activity and resulting variability and potentially anomalous behavior of AP parameters.

OBJECTIVE

To demonstrate the effect of using an in silico interface on electronically expressed I(K1), a major component lacking in h-iPSC-derived cardiac myocytes.

METHODS

An in silico interface was developed to express synthetic I(K1) in cells under whole-cell voltage clamp.

RESULTS

Electronic I(K1) expression established a physiological resting potential, eliminated spontaneous activity, reduced spontaneous early and delayed afterdepolarizations, and decreased AP variability. The initiated APs had the classic rapid upstroke and spike and dome morphology consistent with data obtained with freshly isolated human myocytes as well as the readily recognizable repolarization attributes of ventricular and atrial cells. The application of 1 µM of BayK-8644 resulted in anomalous AP shortening in h-iPSC-derived cardiac myocytes. When I(K1) was electronically expressed, BayK-8644 prolonged the AP, which is consistent with the existing results on native cardiac myocytes.

CONCLUSIONS

The electronic expression of I(K1) is a simple and robust method to significantly improve the physiological behavior of the AP and electrical profile of h-iPSC-derived cardiac myocytes. Increased stability enables the use of this preparation for a controlled quantitative analysis of AP parameters, for example, drug responsiveness, genetic disorders, and dynamic behavior restitution profiles.