All-Optical Electrophysiology Refines Populations of In Silico Human iPSC-CMs for Drug Evaluation

Computational modeling of cardiac electrophysiology and arrhythmogenesis: toward clinical translation

The complexity of cardiac electrophysiology, involving dynamic changes in numerous components across multiple spatial (from ion channel to organ) and temporal (from milliseconds to days) scales, makes an intuitive or empirical analysis of cardiac arrhythmogenesis challenging. Multiscale mechanistic computational models of cardiac electrophysiology provide precise control over individual parameters, and their reproducibility enables a thorough assessment of arrhythmia mechanisms. This review provides a comprehensive analysis of models of cardiac electrophysiology and arrhythmias, from the single cell to the organ level, and how they can be leveraged to better understand rhythm disorders in cardiac disease and to improve heart patient care. Key issues related to model development based on experimental data are discussed, and major families of human cardiomyocyte models and their applications are highlighted. An overview of organ-level computational modeling of cardiac electrophysiology and its clinical applications in personalized arrhythmia risk assessment and patient-specific therapy of atrial and ventricular arrhythmias is provided. The advancements presented here highlight how patient-specific computational models of the heart reconstructed from patient data have achieved success in predicting risk of sudden cardiac death and guiding optimal treatments of heart rhythm disorders. Finally, an outlook toward potential future advances, including the combination of mechanistic modeling and machine learning/artificial intelligence, is provided. As the field of cardiology is embarking on a journey toward precision medicine, personalized modeling of the heart is expected to become a key technology to guide pharmaceutical therapy, deployment of devices, and surgical interventions.

In silico study of the mechanisms of hypoxia and contractile dysfunction during ischemia and reperfusion of hiPSC cardiomyocytes

Interconnected mechanisms of ischemia and reperfusion (IR) has increased the interest in IR in vitro experiments using human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs). We developed a whole-cell computational model of hiPSC-CMs including the electromechanics, a metabolite-sensitive sarcoplasmic reticulum Ca2+-ATPase (SERCA) and an oxygen dynamics formulation to investigate IR mechanisms. Moreover, we simulated the effect and action mechanism of levosimendan, which recently showed promising anti-arrhythmic effects in hiPSC-CMs in hypoxia. The model was validated using hiPSC-CM and in vitro animal data. The role of SERCA in causing relaxation dysfunction in IR was anticipated to be comparable to its function in sepsis-induced heart failure. Drug simulations showed that levosimendan counteracts the relaxation dysfunction by utilizing a particular Ca2+-sensitizing mechanism involving Ca2+-bound troponin C and Ca2+ flux to the myofilament, rather than inhibiting SERCA phosphorylation. The model demonstrates extensive characterization and promise for drug development, making it suitable for evaluating IR therapy strategies based on the changing levels of cardiac metabolites, oxygen and molecular pathways.

Human sodium current voltage-dependence at physiological temperature measured by coupling a patch-clamp experiment to a mathematical model

Voltage-gated Na+ channels are crucial to action potential propagation in excitable tissues. Because of the high amplitude and rapid activation of the Na+ current, voltage-clamp measurements are very challenging and are usually performed at room temperature. In this study, we measured Na+ current voltage-dependence in stem cell-derived cardiomyocytes at physiological temperature. While the apparent activation and inactivation curves, measured as the dependence of current amplitude on voltage, fall within the range reported in previous studies, we identified a systematic error in our measurements. This error is caused by the deviation of the membrane potential from the command potential of the amplifier. We demonstrate that it is possible to account for this artifact using computer simulation of the patch-clamp experiment. We obtained surprising results through patch-clamp model optimization: a half-activation of -11.5 mV and a half-inactivation of -87 mV. Although the half-activation deviates from previous research, we demonstrate that this estimate reproduces the conduction velocity dependence on extracellular potassium concentration. KEY POINTS: Voltage-gated Na+ currents play a crucial role in excitable tissues including neurons, cardiac and skeletal muscle. Measurement of Na+ current is challenging because of its high amplitude and rapid kinetics, especially at physiological temperature. We have used the patch-clamp technique to measure human Na+ current voltage-dependence in human induced pluripotent stem cell-derived cardiomyocytes. The patch-clamp data were processed by optimization of the model accounting for voltage-clamp experiment artifacts, revealing a large difference between apparent parameters of Na+ current and the results of the optimization. We conclude that actual Na+ current activation is extremely depolarized in comparison to previous studies. The new Na+ current model provides a better understanding of action potential propagation; we demonstrate that it explains propagation in hyperkalaemic conditions.

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.

CRISPRi gene modulation and all-optical electrophysiology in post-differentiated human iPSC-cardiomyocytes

Uncovering gene-phenotype relationships can be enabled by precise gene modulation in human induced pluripotent stem-cell-derived cardiomyocytes (iPSC-CMs) and follow up phenotyping using scalable all-optical electrophysiology platforms. Such efforts towards human functional genomics can be aided by recent CRISPR-derived technologies for reversible gene inhibition or activation (CRISPRi/a). We set out to characterize the performance of CRISPRi in post-differentiated iPSC-CMs, targeting key cardiac ion channel genes, KCNH2, KCNJ2, and GJA1, and providing a multiparametric quantification of the effects on cardiac repolarization, stability of the resting membrane potential and conduction properties using all-optical tools. More potent CRISPRi effectors, e.g., Zim3, and optimized viral delivery led to improved performance on par with the use of CRISPRi iPSC lines. Confirmed mild yet specific phenotype changes when CRISPRi is deployed in non-dividing differentiated heart cells is an important step towards more holistic pre-clinical cardiotoxicity testing and for future therapeutic use in vivo.

Building blocks of microphysiological system to model physiology and pathophysiology of human heart

Microphysiological systems (MPS) are drawing increasing interest from academia and from biomedical industry due to their improved capability to capture human physiology. MPS offer an advanced in vitro platform that can be used to study human organ and tissue level functions in health and in diseased states more accurately than traditional single cell cultures or even animal models. Key features in MPS include microenvironmental control and monitoring as well as high biological complexity of the target tissue. To reach these qualities, cross-disciplinary collaboration from multiple fields of science is required to build MPS. Here, we review different areas of expertise and describe essential building blocks of heart MPS including relevant cardiac cell types, supporting matrix, mechanical stimulation, functional measurements, and computational modelling. The review presents current methods in cardiac MPS and provides insights for future MPS development with improved recapitulation of human physiology.

The importance of mechanical conditions in the testing of excitation abnormalities in a population of electro-mechanical models of human ventricular cardiomyocytes

Background: Populations of in silico electrophysiological models of human cardiomyocytes represent natural variability in cell activity and are thoroughly calibrated and validated using experimental data from the human heart. The models have been shown to predict the effects of drugs and their pro-arrhythmic risks. However, excitation and contraction are known to be tightly coupled in the myocardium, with mechanical loads and stretching affecting both mechanics and excitation through mechanisms of mechano-calcium-electrical feedback. However, these couplings are not currently a focus of populations of cell models. Aim: We investigated the role of cardiomyocyte mechanical activity under different mechanical conditions in the generation, calibration, and validation of a population of electro-mechanical models of human cardiomyocytes. Methods: To generate a population, we assumed 11 input parameters of ionic currents and calcium dynamics in our recently developed TP + M model as varying within a wide range. A History matching algorithm was used to generate a non-implausible parameter space by calibrating the action potential and calcium transient biomarkers against experimental data and rejecting models with excitation abnormalities. The population was further calibrated using experimental data on human myocardial force characteristics and mechanical tests involving variations in preload and afterload. Models that passed the mechanical tests were validated with additional experimental data, including the effects of drugs with high or low pro-arrhythmic risk. Results: More than 10% of the models calibrated on electrophysiological data failed mechanical tests and were rejected from the population due to excitation abnormalities at reduced preload or afterload for cell contraction. The final population of accepted models yielded action potential, calcium transient, and force/shortening outputs consistent with experimental data. In agreement with experimental and clinical data, the models demonstrated a high frequency of excitation abnormalities in simulations of Dofetilide action on the ionic currents, in contrast to Verapamil. However, Verapamil showed a high frequency of failed contractions at high concentrations. Conclusion: Our results highlight the importance of considering mechanoelectric coupling in silico cardiomyocyte models. Mechanical tests allow a more thorough assessment of the effects of interventions on cardiac function, including drug testing.

Action potential metrics and automated data analysis pipeline for cardiotoxicity testing using optically mapped hiPSC-derived 3D cardiac microtissues

Recent advances in human induced pluripotent stem cell (hiPSC)-derived cardiac microtissues provide a unique opportunity for cardiotoxic assessment of pharmaceutical and environmental compounds. Here, we developed a series of automated data processing algorithms to assess changes in action potential (AP) properties for cardiotoxicity testing in 3D engineered cardiac microtissues generated from hiPSC-derived cardiomyocytes (hiPSC-CMs). Purified hiPSC-CMs were mixed with 5-25% human cardiac fibroblasts (hCFs) under scaffold-free conditions and allowed to self-assemble into 3D spherical microtissues in 35-microwell agarose gels. Optical mapping was performed to quantify electrophysiological changes. To increase throughput, AP traces from 4x4 cardiac microtissues were simultaneously acquired with a voltage sensitive dye and a CMOS camera. Individual microtissues showing APs were identified using automated thresholding after Fourier transforming traces. An asymmetric least squares method was used to correct non-uniform background and baseline drift, and the fluorescence was normalized (ΔF/F0). Bilateral filtering was applied to preserve the sharpness of the AP upstroke. AP shape changes under selective ion channel block were characterized using AP metrics including stimulation delay, rise time of AP upstroke, APD30, APD50, APD80, APDmxr (maximum rate change of repolarization), and AP triangulation (APDtri = APDmxr-APD50). We also characterized changes in AP metrics under various ion channel block conditions with multi-class logistic regression and feature extraction using principal component analysis of human AP computer simulations. Simulation results were validated experimentally with selective pharmacological ion channel blockers. In conclusion, this simple and robust automated data analysis pipeline for evaluating key AP metrics provides an excellent in vitro cardiotoxicity testing platform for a wide range of environmental and pharmaceutical compounds.

Portable low-cost macroscopic mapping system for all-optical cardiac electrophysiology

SIGNIFICANCE

All-optical cardiac electrophysiology enables the visualization and control of key parameters relevant to the detection of cardiac arrhythmias. Mapping such responses in human induced pluripotent stem-cell-derived cardiomyocytes (hiPSC-CMs) is of great interest for cardiotoxicity and personalized medicine applications.

AIM

We introduce and validate a very low-cost compact mapping system for macroscopic all-optical electrophysiology in layers of hiPSC-CMs.

APPROACH

The system uses oblique transillumination, low-cost cameras, light-emitting diodes, and off-the-shelf components (total < $ 15 , 000 ) to capture voltage, calcium, and mechanical waves under electrical or optical stimulation.

RESULTS

Our results corroborate the equivalency of electrical and optogenetic stimulation of hiPSC-CMs, and V m - [ Ca 2 + ] i similarity in conduction under pacing. Green-excitable optical sensors are combinable with blue optogenetic actuators (chanelrhodopsin2) only under very low green light ( < 0.05    mW / mm 2 ). Measurements in warmer culture medium yield larger spread of action potential duration and higher conduction velocities compared to Tyrode's solution at room temperature.

CONCLUSIONS

As multiple optical sensors and actuators are combined, our results can help handle the "spectral congestion" and avoid parameter distortion. We illustrate the utility of the system for uncovering the action of cellular uncoupling agents and show extensibility to an epi-illumination mode for future imaging of thicker native or engineered tissues.

Guidelines for assessment of cardiac electrophysiology and arrhythmias in small animals

Cardiac arrhythmias are a major cause of morbidity and mortality worldwide. Although recent advances in cell-based models, including human-induced pluripotent stem cell-derived cardiomyocytes (iPSC-CM), are contributing to our understanding of electrophysiology and arrhythmia mechanisms, preclinical animal studies of cardiovascular disease remain a mainstay. Over the past several decades, animal models of cardiovascular disease have advanced our understanding of pathological remodeling, arrhythmia mechanisms, and drug effects and have led to major improvements in pacing and defibrillation therapies. There exist a variety of methodological approaches for the assessment of cardiac electrophysiology and a plethora of parameters may be assessed with each approach. This guidelines article will provide an overview of the strengths and limitations of several common techniques used to assess electrophysiology and arrhythmia mechanisms at the whole animal, whole heart, and tissue level with a focus on small animal models. We also define key electrophysiological parameters that should be assessed, along with their physiological underpinnings, and the best methods with which to assess these parameters.

Cardiac cells stimulated with an axial current-like waveform reproduce electrophysiological properties of tissue fibers

BACKGROUND AND OBJECTIVE

In silico electrophysiological models are generally validated by comparing simulated results with experimental data. When dealing with single-cell and tissue scales simultaneously, as occurs frequently during model development and calibration, the effects of inter-cellular coupling should be considered to ensure the trustworthiness of model predictions. The hypothesis of this paper is that the cell-tissue mismatch can be reduced by incorporating the effects of conduction into the single-cell stimulation current.

METHODS

Five different stimulation waveforms were applied to the human ventricular O'Hara-Rudy cell model. The waveforms included the commonly used monophasic and biphasic (symmetric and asymmetric) pulses, a triangular waveform and a newly proposed asymmetric waveform (stimulation A) that resembles the transmembrane current associated with AP conduction in tissue. A comparison between single-cell and fiber simulated results was established by computing the relative difference between the values of AP-derived properties at different scales, and by evaluating the differences in the contributions of ionic conductances to each evaluated property. As a proof of the benefit, we investigated multi-scale differences in the simulation of the effects induced by dofetilide, a selective I_Kr blocker with high torsadogenic risk, on ventricular repolarization at different pacing rates.

RESULTS

Out of the five tested stimulation waveforms, stimulation A produced the closest correspondence between cell and tissue simulations in terms of AP properties at steady-state and under dynamic pacing and of ionic contributors to those AP properties. Also, stimulation A reproduced the effects of dofetilide better than the other alternative waveforms, mirroring the 'beat-skipping' behavior observed at fast pacing rates in experiments with human tissue.

CONCLUSIONS

The proposed stimulation current waveform accounts for inter-cellular coupling effects by mimicking cell excitation during AP conduction. The proposed waveform improves the correspondence between simulation scales, which could improve the trustworthiness of single-cell simulations without adding computational cost.

Altered contractility in mutation-specific hypertrophic cardiomyopathy: A mechano-energetic in silico study with pharmacological insights

Introduction: Mavacamten (MAVA), Blebbistatin (BLEB), and Omecamtiv mecarbil (OM) are promising drugs directly targeting sarcomere dynamics, with demonstrated efficacy against hypertrophic cardiomyopathy (HCM) in (pre)clinical trials. However, the molecular mechanism affecting cardiac contractility regulation, and the diseased cell mechano-energetics are not fully understood yet. Methods: We present a new metabolite-sensitive computational model of human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) electromechanics to investigate the pathology of R403Q HCM mutation and the effect of MAVA, BLEB, and OM on the cell mechano-energetics. Results: We offer a mechano-energetic HCM calibration of the model, capturing the prolonged contractile relaxation due to R403Q mutation (∼33%), without assuming any further modifications such as an additional Ca²⁺ flux to the thin filaments. The HCM model variant correctly predicts the negligible alteration in ATPase activity in R403Q HCM condition compared to normal hiPSC-CMs. The simulated inotropic effects of MAVA, OM, and BLEB, along with the ATPase activities in the control and HCM model variant agree with in vitro results from different labs. The proposed model recapitulates the tension-Ca²⁺ relationship and action potential duration change due to 1 µM OM and 5 µM BLEB, consistently with in vitro data. Finally, our model replicates the experimental dose-dependent effect of OM and BLEB on the normalized isometric tension. Conclusion: This work is a step toward deep-phenotyping the mutation-specific HCM pathophysiology, manifesting as altered interfilament kinetics. Accordingly, the modeling efforts lend original insights into the MAVA, BLEB, and OM contributions to a new interfilament balance resulting in a cardioprotective effect.

Cross clinical-experimental-computational qualification of in silico drug trials on human cardiac purkinje cells for proarrhythmia risk prediction

The preclinical identification of drug-induced cardiotoxicity and its translation into human risk are still major challenges in pharmaceutical drug discovery. The ICH S7B Guideline and Q&A on Clinical and Nonclinical Evaluation of QT/QTc Interval Prolongation and Proarrhythmic Potential promotes human in silico drug trials as a novel tool for proarrhythmia risk assessment. To facilitate the use of in silico data in regulatory submissions, explanatory control compounds should be tested and documented to demonstrate consistency between predictions and the historic validation data. This study aims to quantify drug-induced electrophysiological effects on in silico cardiac human Purkinje cells, to compare them with existing in vitro rabbit data, and to assess their accuracy for clinical pro-arrhythmic risk predictions. The effects of 14 reference compounds were quantified in simulations with a population of in silico human cardiac Purkinje models. For each drug dose, five electrophysiological biomarkers were quantified at three pacing frequencies, and results compared with available in vitro experiments and clinical proarrhythmia reports. Three key results were obtained: 1) In silico, repolarization abnormalities in human Purkinje simulations predicted drug-induced arrhythmia for all risky compounds, showing higher predicted accuracy than rabbit experiments; 2) Drug-induced electrophysiological changes observed in human-based simulations showed a high degree of consistency with in vitro rabbit recordings at all pacing frequencies, and depolarization velocity and action potential duration were the most consistent biomarkers; 3) discrepancies observed for dofetilide, sotalol and terfenadine are mainly caused by species differences between humans and rabbit. Taken together, this study demonstrates higher accuracy of in silico methods compared to in vitro animal models for pro-arrhythmic risk prediction, as well as a high degree of consistency with in vitro experiments commonly used in safety pharmacology, supporting the potential for industrial and regulatory adoption of in silico trials for proarrhythmia prediction.

Basic Research Approaches to Evaluate Cardiac Arrhythmia in Heart Failure and Beyond

Patients with heart failure often develop cardiac arrhythmias. The mechanisms and interrelations linking heart failure and arrhythmias are not fully understood. Historically, research into arrhythmias has been performed on affected individuals or in vivo (animal) models. The latter however is constrained by interspecies variation, demands to reduce animal experiments and cost. Recent developments in in vitro induced pluripotent stem cell technology and in silico modelling have expanded the number of models available for the evaluation of heart failure and arrhythmia. An agnostic approach, combining the modalities discussed here, has the potential to improve our understanding for appraising the pathology and interactions between heart failure and arrhythmia and can provide robust and validated outcomes in a variety of research settings. This review discusses the state of the art models, methodologies and techniques used in the evaluation of heart failure and arrhythmia and will highlight the benefits of using them in combination. Special consideration is paid to assessing the pivotal role calcium handling has in the development of heart failure and arrhythmia.

Optogenetic Control of Human Induced Pluripotent Stem Cell-Derived Cardiac Tissue Models

Background

Optogenetics, using light‐sensitive proteins, emerged as a unique experimental paradigm to modulate cardiac excitability. We aimed to develop high‐resolution optogenetic approaches to modulate electrical activity in 2‐ and 3‐dimensional cardiac tissue models derived from human induced pluripotent stem cell (hiPSC)‐derived cardiomyocytes.

Methods and Results

To establish light‐controllable cardiac tissue models, opsin‐carrying HEK293 cells, expressing the light‐sensitive cationic‐channel CoChR, were mixed with hiPSC‐cardiomyocytes to generate 2‐dimensional hiPSC‐derived cardiac cell‐sheets or 3‐dimensional engineered heart tissues. Complex illumination patterns were designed with a high‐resolution digital micro‐mirror device. Optical mapping and force measurements were used to evaluate the tissues' electromechanical properties. The ability to optogenetically pace and shape the tissue's conduction properties was demonstrated by using single or multiple illumination stimulation sites, complex illumination patterns, or diffuse illumination. This allowed to establish in vitro models for optogenetic‐based cardiac resynchronization therapy, where the electrical activation could be synchronized (hiPSC‐derived cardiac cell‐sheets and engineered heart tissue models) and contractile properties improved (engineered heart tissues). Next, reentrant activity (rotors) was induced in the hiPSC‐derived cardiac cell‐sheets and engineered heart tissue models through optogenetics programmed‐ or cross‐field stimulations. Diffuse illumination protocols were then used to terminate arrhythmias, demonstrating the potential to study optogenetics cardioversion mechanisms and to identify optimal illumination parameters for arrhythmia termination.

Conclusions

By combining optogenetics and hiPSC technologies, light‐controllable human cardiac tissue models could be established, in which tissue excitability can be modulated in a functional, reversible, and localized manner. This approach may bring a unique value for physiological/pathophysiological studies, for disease modeling, and for developing optogenetic‐based cardiac pacing, resynchronization, and defibrillation approaches.

Human iPSC-Cardiomyocytes as an Experimental Model to Study Epigenetic Modifiers of Electrophysiology

The epigenetic landscape and the responses to pharmacological epigenetic regulators in each human are unique. Classes of epigenetic writers and erasers, such as histone acetyltransferases, HATs, and histone deacetylases, HDACs, control DNA acetylation/deacetylation and chromatin accessibility, thus exerting transcriptional control in a tissue- and person-specific manner. Rapid development of novel pharmacological agents in clinical testing—HDAC inhibitors (HDACi)—targets these master regulators as common means of therapeutic intervention in cancer and immune diseases. The action of these epigenetic modulators is much less explored for cardiac tissue, yet all new drugs need to be tested for cardiotoxicity. To advance our understanding of chromatin regulation in the heart, and specifically how modulation of DNA acetylation state may affect functional electrophysiological responses, human-induced pluripotent stem-cell-derived cardiomyocyte (hiPSC-CM) technology can be leveraged as a scalable, high-throughput platform with ability to provide patient-specific insights. This review covers relevant background on the known roles of HATs and HDACs in the heart, the current state of HDACi development, applications, and any adverse cardiac events; it also summarizes relevant differential gene expression data for the adult human heart vs. hiPSC-CMs along with initial transcriptional and functional results from using this new experimental platform to yield insights on epigenetic control of the heart. We focus on the multitude of methodologies and workflows needed to quantify responses to HDACis in hiPSC-CMs. This overview can help highlight the power and the limitations of hiPSC-CMs as a scalable experimental model in capturing epigenetic responses relevant to the human heart.

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.

Analysis of the response of human iPSC-derived cardiomyocyte tissue to ICaL block. A combined in vitro and in silico approach

The high incidence of cardiac arrythmias underlines the need for the assessment of pharmacological therapies. In this field of drug efficacy, as in the field of drug safety highlighted by the Comprehensive in Vitro Proarrhythmia Assay initiative, new pillars for research have become crucial: firstly, the integration of in-silico experiments, and secondly the evaluation of fully integrated biological systems, such as human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs). In this study, we therefore aimed to combine in-vitro experiments and in-silico simulations to evaluate the antiarrhythmic effect of L-type calcium current (I_CaL) block in hiPSC-CMs. For this, hiPSC-CM preparations were cultured and an equivalent virtual tissue was modeled. Re-entry patterns of electrical activation were induced and several biomarkers were obtained before and after I_CaL block. The virtual hiPSC-CM simulations were also reproduced using a tissue composed of adult ventricular cardiomyocytes (hAdultV-CMs). The analysis of phases, currents and safety factor for propagation showed an increased size of the re-entry core when I_CaL was blocked as a result of depressed cellular excitability. The bigger wavefront curvature yielded reductions of 12.2%, 6.9%, and 4.2% in the frequency of the re-entry for hiPSC-CM cultures, virtual hiPSC-CM, and hAdultV-CM tissues, respectively. Furthermore, I_CaL block led to a 47.8% shortening of the vulnerable window for re-entry in the virtual hiPSC-CM tissue and to re-entry vanishment in hAdultV-CM tissue. The consistent behavior between in-vitro and in-silico hiPSC-CMs and between in-silico hiPSC-CMs and hAdultV-CMs evidences that virtual hiPSC-CM tissues are suitable for assessing cardiac efficacy, as done in the present study through the analysis of I_CaL block.

CalTrack: High-Throughput Automated Calcium Transient Analysis in Cardiomyocytes

[Figure: see text].

An in silico hiPSC-Derived Cardiomyocyte Model Built With Genetic Algorithm

The formulation of in silico biophysical models generally requires optimization strategies for reproducing experimentally observed phenomena. In electrophysiological modeling, robust nonlinear regressive methods are often crucial for guaranteeing high fidelity models. Human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs), though nascent, have proven to be useful in cardiac safety pharmacology, regenerative medicine, and in the implementation of patient-specific test benches for investigating inherited cardiac disorders. This study demonstrates the potency of heuristic techniques at formulating biophysical models, with emphasis on a hiPSC-CM model using a novel genetic algorithm (GA) recipe we proposed. The proposed GA protocol was used to develop a hiPSC-CM biophysical computer model by fitting mathematical formulations to experimental data for five ionic currents recorded in hiPSC-CMs. The maximum conductances of the remaining ionic channels were scaled based on recommendations from literature to accurately reproduce the experimentally observed hiPSC-CM action potential (AP) metrics. Near-optimal parameter fitting was achieved for the GA-fitted ionic currents. The resulting model recapitulated experimental AP parameters such as AP durations (APD₅₀, APD₇₅, and APD₉₀), maximum diastolic potential, and frequency of automaticity. The outcome of this work has implications for validating the biophysics of hiPSC-CMs in their use as viable substitutes for human cardiomyocytes, particularly in cardiac safety pharmacology and in the study of inherited cardiac disorders. This study presents a novel GA protocol useful for formulating robust numerical biophysical models. The proposed protocol is used to develop a hiPSC-CM model with implications for cardiac safety pharmacology.

Precision Medicine and cardiac channelopathies: when dreams meet reality

Precision Medicine (PM) is an innovative approach that, by relying on large populations’ datasets, patients’ genetics and characteristics, and advanced technologies, aims at improving risk stratification and at identifying patient-specific management through targeted diagnostic and therapeutic strategies. Cardiac channelopathies are being progressively involved in the evolution brought by PM and some of them are benefiting from these novel approaches, especially the long QT syndrome. Here, we have explored the main layers that should be considered when developing a PM approach for cardiac channelopathies, with a focus on modern in vitro strategies based on patient-specific human-induced pluripotent stem cells and on in silico models. PM is where scientists and clinicians must meet and integrate their expertise to improve medical care in an innovative way but without losing common sense. We have indeed tried to provide the cardiologist’s point of view by comparing state-of-the-art techniques and approaches, including revolutionary discoveries, to current practice. This point matters because the new approaches may, or may not, exceed the efficacy and safety of established therapies. Thus, our own eagerness to implement the most recent translational strategies for cardiac channelopathies must be tempered by an objective assessment to verify whether the PM approaches are indeed making a difference for the patients. We believe that PM may shape the diagnosis and treatment of cardiac channelopathies for years to come. Nonetheless, its potential superiority over standard therapies should be constantly monitored and assessed before translating intellectually rewarding new discoveries into clinical practice.

Comparison of the Simulated Response of Three in Silico Human Stem Cell-Derived Cardiomyocytes Models and in Vitro Data Under 15 Drug Actions

Objectives: Improvements in human stem cell-derived cardiomyocyte (hSC-CM) technology have promoted their use for drug testing and disease investigations. Several in silico hSC-CM models have been proposed to augment interpretation of experimental findings through simulations. This work aims to assess the response of three hSC-CM in silico models (Koivumäki2018, Kernik2019, and Paci2020) to simulated drug action, and compare simulation results against in vitro data for 15 drugs.

Methods: First, simulations were conducted considering 15 drugs, using a simple pore-block model and experimental data for seven ion channels. Similarities and differences were analyzed in the in silico responses of the three models to drugs, in terms of Ca²⁺ transient duration (CTD₉₀) and occurrence of arrhythmic events. Then, the sensitivity of each model to different degrees of blockage of Na⁺ (I_Na), L-type Ca²⁺ (I_CaL), and rapid delayed rectifying K⁺ (I_Kr) currents was quantified. Finally, we compared the drug-induced effects on CTD₉₀ against the corresponding in vitro experiments.

Results: The observed CTD₉₀ changes were overall consistent among the in silico models, all three showing changes of smaller magnitudes compared to the ones measured in vitro. For example, sparfloxacin 10 µM induced +42% CTD₉₀ prolongation in vitro, and +17% (Koivumäki2018), +6% (Kernik2019), and +9% (Paci2020) in silico. Different arrhythmic events were observed following drug application, mainly for drugs affecting I_Kr. Paci2020 and Kernik2019 showed only repolarization failure, while Koivumäki2018 also displayed early and delayed afterdepolarizations. The spontaneous activity was suppressed by Na⁺ blockers and by drugs with similar effects on I_CaL and I_Kr in Koivumäki2018 and Paci2020, while only by strong I_CaL blockers, e.g. nisoldipine, in Kernik2019. These results were confirmed by the sensitivity analysis.

Conclusion: To conclude, The CTD₉₀ changes observed in silico are qualitatively consistent with our in vitro data, although our simulations show differences in drug responses across the hSC-CM models, which could stem from variability in the experimental data used in their construction.

In vitro and In silico Models to Study SARS-CoV-2 Infection: Integrating Experimental and Computational Tools to Mimic "COVID-19 Cardiomyocyte"

The rapid dissemination of SARS-CoV-2 has made COVID-19 a tremendous social, economic, and health burden. Despite the efforts to understand the virus and treat the disease, many questions remain unanswered about COVID-19 mechanisms of infection and progression. Severe Acute Respiratory Syndrome (SARS) infection can affect several organs in the body including the heart, which can result in thromboembolism, myocardial injury, acute coronary syndromes, and arrhythmias. Numerous cardiac adverse events, from cardiomyocyte death to secondary effects caused by exaggerated immunological response against the virus, have been clinically reported. In addition to the disease itself, repurposing of treatments by using "off label" drugs can also contribute to cardiotoxicity. Over the past several decades, animal models and more recently, stem cell-derived cardiomyocytes have been proposed for studying diseases and testing treatments in vitro. In addition, mechanistic in silico models have been widely used for disease and drug studies. In these models, several characteristics such as gender, electrolyte imbalance, and comorbidities can be implemented to study pathophysiology of cardiac diseases and to predict cardiotoxicity of drug treatments. In this Mini Review, we (1) present the state of the art of in vitro and in silico cardiomyocyte modeling currently in use to study COVID-19, (2) review in vitro and in silico models that can be adopted to mimic the effects of SARS-CoV-2 infection on cardiac function, and (3) provide a perspective on how to combine some of these models to mimic "COVID-19 cardiomyocytes environment.".

Cardiac optogenetics: a decade of enlightenment

The electromechanical function of the heart involves complex, coordinated activity over time and space. Life-threatening cardiac arrhythmias arise from asynchrony in these space-time events; therefore, therapies for prevention and treatment require fundamental understanding and the ability to visualize, perturb and control cardiac activity. Optogenetics combines optical and molecular biology (genetic) approaches for light-enabled sensing and actuation of electrical activity with unprecedented spatiotemporal resolution and parallelism. The year 2020 marks a decade of developments in cardiac optogenetics since this technology was adopted from neuroscience and applied to the heart. In this Review, we appraise a decade of advances that define near-term (immediate) translation based on all-optical electrophysiology, including high-throughput screening, cardiotoxicity testing and personalized medicine assays, and long-term (aspirational) prospects for clinical translation of cardiac optogenetics, including new optical therapies for rhythm control. The main translational opportunities and challenges for optogenetics to be fully embraced in cardiology are also discussed.

Syncytium cell growth increases Kir2.1 contribution in human iPSC-cardiomyocytes

Human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) enable cardiotoxicity testing and personalized medicine. However, their maturity is of concern, including relatively depolarized resting membrane potential and more spontaneous activity compared with adult cardiomyocytes, implicating low or lacking inward rectifier potassium current (I_k1). Here, protein quantification confirms Kir2.1 expression in hiPSC-CM syncytia, albeit several times lower than in adult heart tissue. We find that hiPSC-CM culture density influences Kir2.1 expression at the mRNA level (potassium inwardly rectifying channel subfamily J member 2) and at the protein level and its associated electrophysiology phenotype. Namely, all-optical cardiac electrophysiology and pharmacological treatments reveal reduction of spontaneous and irregular activity and increase in action potential upstroke in denser cultures. Blocking I_k1-like currents with BaCl₂ increased spontaneous frequency and blunted action potential upstrokes during pacing in a dose-dependent manner only in the highest-density cultures, in line with I_k1’s role in regulating the resting membrane potential. Our results emphasize the importance of syncytial growth of hiPSC-CMs for more physiologically relevant phenotype and the power of all-optical electrophysiology to study cardiomyocytes in their multicellular setting.

NEW & NOTEWORTHY We identify cell culture density and cell-cell contact as an important factor in determining the expression of a key ion channel at the transcriptional and the protein levels, KCNJ2/Kir2.1, and its contribution to the electrophysiology of human induced pluripotent stem cell-derived cardiomyocytes. Our results indicate that studies on isolated cells, out of tissue context, may underestimate the cellular ion channel properties being characterized.