Human-based modelling and simulations are becoming ubiquitous in biomedical science due to their ability to augment experimental and clinical investigations. Cardiac electrophysiology is one of the most advanced areas, with cardiac modelling and simulation being considered for virtual testing of pharmacological therapies and medical devices. Current models present inconsistencies with experimental data, which limit further progress. In this study, we present the design, development, calibration and independent validation of a human-based ventricular model (ToR-ORd) for simulations of electrophysiology and excitation-contraction coupling, from ionic to whole-organ dynamics, including the electrocardiogram. Validation based on substantial multiscale simulations supports the credibility of the ToR-ORd model under healthy and key disease conditions, as well as drug blockade. In addition, the process uncovers new theoretical insights into the biophysical properties of the L-type calcium current, which are critical for sodium and calcium dynamics. These insights enable the reformulation of L-type calcium current, as well as replacement of the hERG current model.

Hypertrophic cardiomyopathy (HCM) is a cause of sudden arrhythmic death, but the understanding of its pro-arrhythmic mechanisms and an effective pharmacological treatment are lacking. HCM electrophysiological remodelling includes both increased inward and reduced outward currents, but their role in promoting repolarisation abnormalities remains unknown. The goal of this study is to identify key ionic mechanisms driving repolarisation abnormalities in human HCM, and to evaluate anti-arrhythmic effects of single and multichannel inward current blocks. Experimental ionic current, action potential (AP) and Ca(2+)-transient (CaT) recordings were used to construct populations of human non-diseased and HCM AP models (n=9118), accounting for inter-subject variability. Simulations were conducted for several degrees of selective and combined inward current block. Simulated HCM cardiomyocytes exhibited prolonged AP and CaT, diastolic Ca(2+) overload and decreased CaT amplitude, in agreement with experiments. Repolarisation abnormalities in HCM models were consistently driven by L-type Ca(2+) current (ICaL) re-activation, and ICaL block was the most effective intervention to normalise repolarisation and diastolic Ca(2+), but compromised CaT amplitude. Late Na(+) current (INaL) block partially abolished repolarisation abnormalities, with small impact on CaT. Na(+)/Ca(2+) exchanger (INCX) block effectively restored repolarisation and CaT amplitude, but increased Ca(2+) overload. Multichannel block increased efficacy in normalising repolarisation, AP biomarkers and CaT amplitude compared to selective block. Experimentally-calibrated populations of human AP models identify ICaL re-activation as the key mechanism for repolarisation abnormalities in HCM, and combined INCX, INaL and ICaL block as effective anti-arrhythmic therapies also able to partially reverse the HCM electrophysiological phenotype.

Calmodulin (CaM) is a small protein, encoded by three genes (CALM1-3), exerting multiple Ca²⁺-dependent modulatory roles. A mutation (F142L) affecting only one of the six CALM alleles is associated with long QT syndrome (LQTS) characterized by recurrent cardiac arrests. This phenotypic severity is unexpected from the predicted allelic balance. In this work, the effects of heterozygous CALM1-F142L have been investigated in human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) obtained from a LQTS patient carrying the F142L mutation, i.e. in the context of native allelic ratio and potential gene modifiers. Skin fibroblasts of the mutation carrier and two unrelated healthy subjects (controls) were reprogrammed to hiPSC and differentiated into hiPSC-CMs. Scanty I_K1 expression, an hiPSC-CMs feature potentially biasing repolarization, was corrected by addition of simulated I_K1 (Dynamic-Clamp). Abnormalities in repolarization rate-dependency (in single cells and cell aggregates), membrane currents and intracellular Ca²⁺dynamics were evaluated as putative arrhythmogenic factors. CALM1-F142L prolonged repolarization, altered its rate-dependency and its response to isoproterenol. This was associated with severe impairment of Ca²⁺-dependent inactivation (CDI) of I_CaL, resulting in augmented inward current during the plateau phase. As a result, the repolarization of mutant cells failed to adapt to high pacing rates, a finding well reproduced by using a recent hiPSC-CM action potential model. The mutation failed to affect I_Ks and I_NaL and changed I_f only marginally. Intracellular Ca²⁺dynamics and Ca²⁺store stability were not significantly modified. Mutation-induced repolarization abnormalities were reversed by verapamil. The main functional derangement in CALM1-F142L was prolonged repolarization with altered rate-dependency and sensitivity to β-adrenergic stimulation. Impaired CDI of I_CaL underlined the electrical abnormality, which was sensitive to I_CaL blockade. High mutation penetrance was confirmed in the presence of the native genotype, implying strong dominance of effects.

This white paper presents principles for validating proarrhythmia risk prediction models for regulatory use as discussed at the In Silico Breakout Session of a Cardiac Safety Research Consortium/Health and Environmental Sciences Institute/US Food and Drug Administration–sponsored Think Tank Meeting on May 22, 2018. The meeting was convened to evaluate the progress in the development of a new cardiac safety paradigm, the Comprehensive in Vitro Proarrhythmia Assay (CiPA). The opinions regarding these principles reflect the collective views of those who participated in the discussion of this topic both at and after the breakout session. Although primarily discussed in the context of in silico models, these principles describe the interface between experimental input and model‐based interpretation and are intended to be general enough to be applied to other types of nonclinical models for proarrhythmia assessment. This document was developed with the intention of providing a foundation for more consistency and harmonization in developing and validating different models for proarrhythmia risk prediction using the example of the CiPA paradigm.