Caveolin-3 and Caveolin-1 Interaction Decreases Channel Dysfunction Due to Caveolin-3 Mutations.

Caveolae constitute membrane microdomains where receptors and ion channels functionally interact. Caveolin-3 (cav-3) is the key structural component of muscular caveolae. Mutations in CAV3 lead to caveolinopathies, which result in both muscular dystrophies and cardiac diseases. In cardiomyocytes, cav-1 participates with cav-3 to form caveolae; skeletal myotubes and adult skeletal fibers do not express cav-1. In the heart, the absence of cardiac alterations in the majority of cases may depend on a conserved organization of caveolae thanks to the expression of cav-1. We decided to focus on three specific cav-3 mutations (Δ62-64YTT; T78K and W101C) found in heterozygosis in patients suffering from skeletal muscle disorders. We overexpressed both the WT and mutated cav-3 together with ion channels interacting with and modulated by cav-3. Patch-clamp analysis conducted in caveolin-free cells (MEF-KO), revealed that the T78K mutant is dominant negative, causing its intracellular retention together with cav-3 WT, and inducing a significant reduction in current densities of all three ion channels tested. The other cav-3 mutations did not cause significant alterations. Mathematical modelling of the effects of cav-3 T78K would impair repolarization to levels incompatible with life. For this reason, we decided to compare the effects of this mutation in other cell lines that endogenously express cav-1 (MEF-STO and CHO cells) and to modulate cav-1 expression with an shRNA approach. In these systems, the membrane localization of cav-3 T78K was rescued in the presence of cav-1, and the current densities of hHCN4, hKv1.5 and hKir2.1 were also rescued. These results constitute the first evidence of a compensatory role of cav-1 in the heart, justifying the reduced susceptibility of this organ to caveolinopathies.

A detailed mathematical model of the human atrial cardiomyocyte: integration of electrophysiology and cardiomechanics.

Mechano-electric regulations (MER) play an important role in the maintenance of cardiac performance. Mechano-calcium and mechano-electric feedback (MCF and MEF) pathways adjust the cardiomyocyte contractile force according to mechanical perturbations and affects electro-mechanical coupling. MER integrates all these regulations in one unit resulting in a complex phenomenon. Computational modelling is a useful tool to accelerate the mechanistic understanding of complex experimental phenomena. We have developed a novel model that integrates the MER loop for human atrial cardiomyocytes with proper consideration of feedforward and feedback pathways. The model couples a modified version of the action potential (AP) Koivumäki model with the contraction model by Quarteroni group. The model simulates iso-sarcometric and isometric twitches and the feedback effects on AP and Ca2+ -handling. The model showed a biphasic response of Ca2+ transient (CaT) peak to increasing pacing rates and highlights the possible mechanisms involved. The model has shown a shift of the threshold for AP and CaT alternans from 4.6 to 4 Hz under post-operative atrial fibrillation, induced by depressed SERCA activity. The alternans incidence was dependent on a chain of mechanisms including RyRs availability time, MCF coupling, CaMKII phosphorylation, and the stretch levels. As a result, the model predicted a 10% slowdown of conduction velocity for a 20% stretch, suggesting a role of stretch in creation of substrate formation for atrial fibrillation. Overall, we conclude that the developed model provides a physiological CaT followed by a physiological twitch. This model can open pathways for the future studies of human atrial electromechanics. KEY POINTS: With the availability of human atrial cellular data, interest in atrial-specific model integration has been enhanced. We have developed a detailed mathematical model of human atrial cardiomyocytes including the mechano-electric regulatory loop. The model has gone through calibration and evaluation phases against a wide collection of available human in-vitro data. The usefulness of the model for analysing clinical problems has been preliminaryly tested by simulating the increased incidence of Ca2+ transient and action potential alternans at high rates in post-operative atrial fibrillation condition. The model determines the possible role of mechano-electric feedback in alternans incidence, which can increase vulnerability to atrial arrhythmias by varying stretch levels. We found that our physiologically accurate description of Ca2+ handling can reproduce many experimental phenomena and can help to gain insights into the underlying pathophysiological mechanisms.

A dynamic clamping approach using in silico IK1 current for discrimination of chamber-specific hiPSC-derived cardiomyocytes.

Human induced pluripotent stem cell (hiPSC)-derived cardiomyocytes (CM) constitute a mixed population of ventricular-, atrial-, nodal-like cells, limiting the reliability for studying chamber-specific disease mechanisms. Previous studies characterised CM phenotype based on action potential (AP) morphology, but the classification criteria were still undefined. Our aim was to use in silico models to develop an automated approach for discriminating the electrophysiological differences between hiPSC-CM. We propose the dynamic clamp (DC) technique with the injection of a specific IK1 current as a tool for deriving nine electrical biomarkers and blindly classifying differentiated CM. An unsupervised learning algorithm was applied to discriminate CM phenotypes and principal component analysis was used to visualise cell clustering. Pharmacological validation was performed by specific ion channel blocker and receptor agonist. The proposed approach improves the translational relevance of the hiPSC-CM model for studying mechanisms underlying inherited or acquired atrial arrhythmias in human CM, and for screening anti-arrhythmic agents.

The virtual sinoatrial node: What did computational models tell us about cardiac pacemaking?

Since its discovery, the sinoatrial node (SAN) has represented a fascinating and complex matter of research. Despite over a century of discoveries, a full comprehension of pacemaking has still to be achieved. Experiments often produced conflicting evidence that was used either in support or against alternative theories, originating intense debates. In this context, mathematical descriptions of the phenomena underlying the heartbeat have grown in importance in the last decades since they helped in gaining insights where experimental evaluation could not reach. This review presents the most updated SAN computational models and discusses their contribution to our understanding of cardiac pacemaking. Electrophysiological, structural and pathological aspects - as well as the autonomic control over the SAN - are taken into consideration to reach a holistic view of SAN activity.

Coupling and heterogeneity modulate pacemaking capability in healthy and diseased two-dimensional sinoatrial node tissue models.

Both experimental and modeling studies have attempted to determine mechanisms by which a small anatomical region, such as the sinoatrial node (SAN), can robustly drive electrical activity in the human heart. However, despite many advances from prior research, important questions remain unanswered. This study aimed to investigate, through mathematical modeling, the roles of intercellular coupling and cellular heterogeneity in synchronization and pacemaking within the healthy and diseased SAN. In a multicellular computational model of a monolayer of either human or rabbit SAN cells, simulations revealed that heterogenous cells synchronize their discharge frequency into a unique beating rhythm across a wide range of heterogeneity and intercellular coupling values. However, an unanticipated behavior appeared under pathological conditions where perturbation of ionic currents led to reduced excitability. Under these conditions, an intermediate range of intercellular coupling (900-4000 MΩ) was beneficial to SAN automaticity, enabling a very small portion of tissue (3.4%) to drive propagation, with propagation failure occurring at both lower and higher resistances. This protective effect of intercellular coupling and heterogeneity, seen in both human and rabbit tissues, highlights the remarkable resilience of the SAN. Overall, the model presented in this work allowed insight into how spontaneous beating of the SAN tissue may be preserved in the face of perturbations that can cause individual cells to lose automaticity. The simulations suggest that certain degrees of gap junctional coupling protect the SAN from ionic perturbations that can be caused by drugs or mutations.

A Novel In Silico Electromechanical Model of Human Ventricular Cardiomyocyte.

Contractility has become one of the main readouts in computational and experimental studies on cardiomyocytes. Following this trend, we propose a novel mathematical model of human ventricular cardiomyocytes electromechanics, BPSLand, by coupling a recent human contractile element to the BPS2020 model of electrophysiology. BPSLand is the result of a hybrid optimization process and it reproduces all the electrophysiology experimental indices captured by its predecessor BPS2020, simultaneously enabling the simulation of realistic human active tension and its potential abnormalities. The transmural heterogeneity in both electrophysiology and contractility departments was simulated consistent with previous computational and in vitro studies. Furthermore, our model could capture delayed afterdepolarizations (DADs), early afterdepolarizations (EADs), and contraction abnormalities in terms of aftercontractions triggered by either drug action or special pacing modes. Finally, we further validated the mechanical results of the model against previous experimental and in silico studies, e.g., the contractility dependence on pacing rate. Adding a new level of applicability to the normative models of human cardiomyocytes, BPSLand represents a robust, fully-human in silico model with promising capabilities for translational cardiology.

Computational Analysis of Mapping Catheter Geometry and Contact Quality Effects on Rotor Detection in Atrial Fibrillation.

Atrial fibrillation (AF) is the most common cardiac arrhythmia and catheter mapping has been proved to be an effective approach for detecting AF drivers to be targeted by ablation. Among drivers, the so-called rotors have gained the most attention: their identification and spatial location could help to understand which patient-specific mechanisms are acting, and thus to guide the ablation execution. Since rotor detection by multi-electrode catheters may be influenced by several structural parameters including inter-electrode spacing, catheter coverage, and endocardium-catheter distance, in this study we proposed a tool for testing the ability of different catheter shapes to detect rotors in different conditions. An approach based on the solution of the monodomain equations coupled with a modified Courtemanche ionic atrial model, that considers an electrical remodeling, was applied to simulate spiral wave dynamics on a 2D model for 7.75 s. The developed framework allowed the acquisition of unipolar signals at 2 KHz. Two high-density multipolar catheters were simulated (Advisor™ HD Grid and PentaRay®) and placed in a 2D region in which the simulated spiral wave persists longer. The configuration of the catheters was then modified by changing the number of electrodes, inter-electrodes distance, position, and atrial-wall distance for assessing how they would affect the rotor detection. In contact with the wall and at 1 mm distance from it, all the configurations detected the rotor correctly, irrespective of geometry, coverage, and inter-electrode distance. In the HDGrid-like geometry, the increase of the inter-electrode distance from 3 to 6 mm caused rotor detection failure at 2 mm distance from the LA wall. In the PentaRay-like configuration, regardless of inter-electrode distance, rotor detection failed at 3 mm endocardium-catheter distance. The asymmetry of this catheter resulted in rotation-dependent rotor detection. To conclude, the computational framework we developed is based on realistic catheter shapes designed with parameter configurations which resemble clinical settings. Results showed it is well suited to investigate how mapping catheter geometry and location affect AF driver detection, therefore it is a reliable tool to design and test new mapping catheters.

Simulation of the Effects of Extracellular Calcium Changes Leads to a Novel Computational Model of Human Ventricular Action Potential With a Revised Calcium Handling.

The importance of electrolyte concentrations for cardiac function is well established. Electrolyte variations can lead to arrhythmias onset, due to their important role in the action potential (AP) genesis and in maintaining cell homeostasis. However, most of the human AP computer models available in literature were developed with constant electrolyte concentrations, and fail to simulate physiological changes induced by electrolyte variations. This is especially true for Ca2+, even in the O'Hara-Rudy model (ORd), one of the most widely used models in cardiac electrophysiology. Therefore, the present work develops a new human ventricular model (BPS2020), based on ORd, able to simulate the inverse dependence of AP duration (APD) on extracellular Ca2+ concentration ([Ca2+]o), and APD rate dependence at 4 mM extracellular K+. The main changes needed with respect to ORd are: (i) an increased sensitivity of L-type Ca2+ current inactivation to [Ca2+]o; (ii) a single compartment description of the sarcoplasmic reticulum; iii) the replacement of Ca2+ release. BPS2020 is able to simulate the physiological APD-[Ca2+]o relationship, while also retaining the well-reproduced properties of ORd (APD rate dependence, restitution, accommodation and current block effects). We also used BPS2020 to generate an experimentally-calibrated population of models to investigate: (i) the occurrence of repolarization abnormalities in response to hERG current block; (ii) the rate adaptation variability; (iii) the occurrence of alternans and delayed after-depolarizations at fast pacing. Our results indicate that we successfully developed an improved version of ORd, which can be used to investigate electrophysiological changes and pro-arrhythmic abnormalities induced by electrolyte variations and current block at multiple rates and at the population level.

Development, Implementation and Testing of a Multicellular Dynamic Action Potential Clamp Simulator for Drug Cardiac Safety Assessment.

As drugs can be multichannel blockers it is important to assess their cardiac safety taking into account multiple currents. In silico action potential (AP) models have been proposed for being able to integrate drugs effect on ionic currents and generate the resulting AP. However, a mathematical description of drug effects is required, which could be inaccurate. Dynamic Clamp has been proposed for drug cardiac safety assessment. In the dynamic action potential clamp (dAPC) configuration it creates an hybrid model connecting a real cell with a computer simulation. This way, drugs could be administrated directly to real cells, and effects on currents can be taken into account when generating the AP. Here we design and simulate a parallel multichannel dAPC system. The system includes the real cells overexpressing the currents of interest, the voltage clamp acquisition system, and the AP in silico model.

Development, calibration, and validation of a novel human ventricular myocyte model in health, disease, and drug block.

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.

Decades of intensive experimental and clinical research have revealed much about how the human heart works. Though incomplete, this knowledge has been used to construct computer models that represent the activity of this organ as a whole, and of its individual chambers (the atria and ventricles), tissues and cells. Such models have been used to better understand life-threatening irregular heartbeats; they are also beginning to be used to guide decisions about the treatment of patients and the development of new drugs by the pharmaceutical industry. Yet existing computer models of the electrical activity of the human heart are sometimes inconsistent with experimental data. This problem led Tomek et al. to try to create a new model that was consistent with established biophysical knowledge and experimental data for a wide range of conditions including disease and drug action. Tomek et al. designed a strategy that explicitly separated the construction and validation of a model that could recreate the electrical activity of the ventricles in a human heart. This model was able to integrate and explain a wide range of properties of both healthy and diseased hearts, including their response to different drugs. The development of the model also uncovered and resolved theoretical inconsistencies that have been present in almost all models of the heart from the last 25 years. Tomek et al. hope that their new human heart model will enable more basic, translational and clinical research into a range of heart diseases and accelerate the development of new therapies.

Acute effect of a peritoneal dialysis exchange on electrolyte concentration and QT interval in uraemic patients.

BACKGROUND: Hemodialysis (HD) sessions induce changes in plasma electrolytes that lead to modifications of QT interval, virtually associated with dangerous arrhythmias. It is not known whether such a phenomenon occurs even during peritoneal dialysis (PD). The aim of the study is to analyze the relationship between dialysate and plasma electrolyte modifications and QT interval during a PD exchange. METHODS: In 15 patients, two manual PD 4-h exchanges were performed, using two isotonic solutions with different calcium concentration (Ca++1.25 and Ca1.75++ mmol/L). Dialysate and plasma electrolyte concentration and QT interval (ECG Holter recording) were monitored hourly. A computational model simulating the ventricular action potential during the exchange was also performed. RESULTS: Dialysis exchange induced a significant plasma alkalizing effect (p < 0.001). Plasma K+ significantly decreased at the third hour (p < 0.05). Plasma Na+ significantly decreased (p < 0.001), while plasma Ca++ slightly increased only when using the Ca 1.75++ mmol/L solution (p < 0.01). The PD exchange did not induce modifications of clinical relevance in the QT interval, while a significant decrease in heart rate (p < 0.001) was observed. The changes in plasma K+ values were significantly inversely correlated to QT interval modifications (p < 0.001), indicating that even small decreases of K+ were consistently paralleled by small QT prolongations. These results were perfectly confirmed by the computational model. CONCLUSIONS: The PD exchange guarantees a greater cardiac electrical stability compared to the HD session and should be preferred in patients with a higher arrhythmic risk. Moreover, our study shows that ventricular repolarization is extremely sensitive to plasma K+ changes, also in normal range.

Action potential contour contributes to species differences in repolarization response to ß-adrenergic stimulation.

Aim: Repolarization response to ß-adrenergic (ß-AR) stimulation differs between guinea-pig and canine myocytes and, within the latter, between myocardial layers. Correlative analysis suggests that this may be due to differences in action potential (AP) contour. Here we tested whether AP contour may set the response of current and of repolarization to ß-AR stimulation (10 nM isoproterenol, ISO). Methods and results: The responses of AP and current to ISO were measured under I-clamp and "AP-clamp" in guinea-pig (GP), dog epicardial (DEPI) and dog subendocardial (DENDO) myocytes. Dynamic-clamp (DC) was used to evaluate the impact of AP features on AP response to ISO. ISO prolonged AP duration (APD) in GP myocytes, did not affect it in DENDO and shortened it in DEPI ones. The current induced by ISO (IISO) sharply differed between GP and canine myocytes and, to a lesser extent, between DENDO and DEPI ones. Differences in IISO profile likely important in setting APD response (time-to-peak, time-to-reversal), were minimized when canine myocytes where clamped with GP AP-waveforms and vice versa. Introduction of a "notch" in GP AP (by DC) was alone insufficient to affect the APD response to ISO; nevertheless, when incorporated in a GP AP-waveform, the main "canine" AP features ("notch" and low plateau potential) caused IISO of GP myocytes to acquire canine features. Conclusion: Early repolarization contour and level of plateau potential contribute to species-specificity of IISO profile. Changes in AP contour, also when generated by modulation of ISO-insensitive currents, may be crucial in setting APD response to ß-AR stimulation.

The Pioneering Work of Enrico Morselli (1852-1929) in Light of Modern Scientific Research on Hypnosis and Suggestion.

This article examines research on hypnosis and suggestion, starting with the nineteenth-century model proposed by Enrico Morselli (1852-1929), an illustrious Italian psychiatrist and psychologist. The authors conducted an original psychophysiological analysis of hypnosis, distancing the work from the neuropathological concept of the time and proposing a model based on a naturalistic approach to investigating mental processes. The issues investigated by Morselli, including the definition of hypnosis and analysis of specific mental processes such as attention and memory, are reviewed in light of modern research. From the view of modern neuroscientific concepts, some problems that originated in the nineteenth century still appear to be present and pose still-open questions.

D242N, a KV7.1 LQTS mutation uncovers a key residue for IKs voltage dependence.

KV7.1 and KCNE1 co-assemble to give rise to the IKs current, one of the most important repolarizing currents of the cardiac action potential. Its relevance is underscored by the identification of >500 mutations in KV7.1 and, at least, 36 in KCNE1, that cause Long QT Syndrome (LQTS). The aim of this study was to characterize the biophysical and cellular consequences of the D242N KV7.1 mutation associated with the LQTS. The mutation is located in the S4 transmembrane segment, within the voltage sensor of the KV7.1 channel, disrupting the conserved charge balance of this region. Perforated patch-clamp experiments show that, unexpectedly, the mutation did not disrupt the voltage-dependent activation but it removed the inactivation and slowed the activation kinetics of D242N KV7.1 channels. Biotinylation of cell-surface protein and co-immunoprecipitation experiments revealed that neither plasma membrane targeting nor co-assembly between KV7.1 and KCNE1 was altered by the mutation. However, the association of D242N KV7.1 with KCNE1 strongly shifted the voltage dependence of activation to more depolarized potentials (+50mV), hindering IKs current at physiologically relevant membrane potentials. Both functional and computational analysis suggest that the clinical phenotype of the LQTS patients carrying the D242N mutation is due to impaired action potential adaptation to exercise and, in particular, to increase in heart rate. Moreover, our data identify D242 aminoacidic position as a potential residue involved in the KCNE1-mediated regulation of the voltage dependence of activation of the KV7.1 channel.

The archive of the History of Psychology at the University of Rome, Sapienza.

The History of Psychology Archive at the University of Rome, Sapienza was founded in 2008 in the Department of Dynamic and Clinical Psychology. The archive aspires to become an indispensable tool to (a) understand the currents, schools, and research traditions that have marked the path of Italian psychology, (b) focus on issues of general and applied psychology developed in each university, (c) identify experimental and clinical-differential methodologies specific to each lab, (d) reconstruct the genesis and consolidation of psychology institutions and, ultimately, (e) write a "story," set according to the most recent historiographical criteria. The archive is designed according to scholarship on the history of Italian psychology from the past two decades. The online archive is divided into five sections for ease of access. The Sapienza archive is a work in progress and it has plans for expansion. (PsycINFO Database Record

A new KCNQ1 mutation at the S5 segment that impairs its association with KCNE1 is responsible for short QT syndrome.

AIMS: KCNQ1 and KCNE1 encode Kv7.1 and KCNE1, respectively, the pore-forming and the accessory subunits of the slow delayed rectifier potassium current, IKs. KCNQ1 mutations are associated with long and short QT syndrome. The aim of this study was to characterize the biophysical and cellular phenotype of a KCNQ1 missense mutation, F279I, found in a 23-year-old man with a corrected QT interval (QTc) of 356 ms and a family history of sudden cardiac death. METHODS AND RESULTS: Experiments were performed using perforated patch-clamp, western blot, co-immunoprecipitation, biotinylation, and immunocytochemistry techniques in HEK293, COS7 cells and in cardiomyocytes transfected with WT Kv7.1/KCNE1 or F279I Kv7.1/KCNE1 channels. In the absence of KCNE1, F279I Kv7.1 current exhibited a lesser degree of inactivation than WT Kv7.1. Also, functional analysis of F279I Kv7.1 in the presence of KCNE1 revealed a negative shift in the activation curve and an acceleration of the activation kinetics leading to a gain of function in IKs. The co-assembly between F279I Kv7.1 channels and KCNE1 was markedly decreased compared with WT Kv7.1 channels, as revealed by co-immunoprecipitation and Föster Resonance Energy Transfer experiments. All these effects contribute to the increase of IKs when channels incorporate F279I Kv7.1 subunits, as shown by a computer model simulation of these data that predicts a shortening of the action potential (AP) consistent with the patient phenotype. CONCLUSION: The F279I mutation induces a gain of function of IKs due to an impaired gating modulation of Kv7.1 induced by KCNE1, leading to a shortening of the cardiac AP.

IKr Impact on Repolarization and Its Variability Assessed by Dynamic Clamp.

BACKGROUND: Repolarization and its stability are exquisitely sensitive to IKr features. Information on the relative importance of specific IKr abnormalities is missing and would assist in the evaluation of arrhythmogenic risk. METHODS AND RESULTS: In single guinea-pig myocytes, endogenous IKr was replaced by modeled IKr (mIKr) by dynamic clamp (DC) at a cycle length of 1 s. mIKr parameters were systematically modified, and the resulting changes in action potential duration (APD) and its short term variability (SD1) were measured. We observed that (1) IKr blockade increased SD1 more than expected by its dependency on APD; (2) mIKr completely reversed APD and SD1 changes caused by IKr blockade; (3) repolarization was most sensitive to inactivation shifts, which affected APD and SD1 concordantly; (4) activation shifts of the same magnitude had marginal impact on APD, but only when reducing mIKr, they significantly increased SD1; (5) changes in maximal conductance resulted in a pattern similar to that of activation shifts. CONCLUSIONS: The largest effect on repolarization and its stability are expected from changes in IKr inactivation. APD is less sensitive to changes in other IKr gating parameters, which are better revealed by SD1 changes. SD1 may be more sensitive than APD in detecting IKr-dependent repolarization abnormalities.

Combined action potential- and dynamic-clamp for accurate computational modelling of the cardiac IKr current.

In the present work Action-Potential clamp (APC) and Dynamic clamp (DC) were used in combination in order to optimize the Luo-Rudy (LRd) mathematical formulation of the guinea-pig rapid delayed rectifier K(+) current (IKr), and to validate the optimized model. To this end, IKr model parameters were adjusted to fit the experimental E4031-sensitive current (IE4031) recorded under APC in guinea-pig myocytes. Currents generated by LRd model (ILRd) and the optimized one (IOpt) were then compared by testing their suitability to replace IE4031 under DC. Under APC, ILRd was significantly larger than IE4031 (mean current densities 0.51±0.01 vs 0.21±0.05pA/pF; p<0.001), mainly because of different rectification. IOpt mean density (0.17±0.01pA/pF) was similar to the IE4031 one (NS); moreover, IOpt accurately reproduced IE4031 distribution along the different AP phases. Models were then compared under DC by blocking native IKr (5µM E4031) and replacing it with ILRd or IOpt. Whereas injection of ILRd overshortened AP duration (APD90) (by 25% of its pre-block value), IOpt injection restored AP morphology and duration to overlap pre-block values. This study highlights the power of APC and DC for the identification of reliable formulations of ionic current models. An optimized model of IKr has been obtained which fully reversed E4031 effects on the AP. The model strongly diverged from the widely used Luo-Rudy formulation; this can be particularly relevant to the in silico analysis of AP prolongation caused by IKr blocking or alterations.