Numerical simulation of electrocardiograms for full cardiac cycles in healthy and pathological conditions

ECG Localization Method Based on Volume Conductor Model and Kalman Filtering

The 12-lead electrocardiogram was invented more than 100 years ago and is still used as an essential tool in the early detection of heart disease. By estimating the time-varying source of the electrical activity from the potential changes, several types of heart disease can be noninvasively identified. However, most previous studies are based on signal processing, and thus an approach that includes physics modeling would be helpful for source localization problems. This study proposes a localization method for cardiac sources by combining an electrical analysis with a volume conductor model of the human body as a forward problem and a sparse reconstruction method as an inverse problem. Our formulation estimates not only the current source location but also the current direction. For a 12-lead electrocardiogram system, a sensitivity analysis of the localization to cardiac volume, tilted angle, and model inhomogeneity was evaluated. Finally, the estimated source location is corrected by Kalman filter, considering the estimated electrocardiogram source as time-sequence data. For a high signal-to-noise ratio (greater than 20 dB), the dominant error sources were the model inhomogeneity, which is mainly attributable to the high conductivity of the blood in the heart. The average localization error of the electric dipole sources in the heart was 12.6 mm, which is comparable to that in previous studies, where a less detailed anatomical structure was considered. A time-series source localization with Kalman filtering indicated that source mislocalization could be compensated, suggesting the effectiveness of the source estimation using the current direction and location simultaneously. For the electrocardiogram R-wave, the mean distance error was reduced to less than 7.3 mm using the proposed method. Considering the physical properties of the human body with Kalman filtering enables highly accurate estimation of the cardiac electric signal source location and direction. This proposal is also applicable to electrode configuration, such as ECG sensing systems.

Parameter variations in personalized electrophysiological models of human heart ventricles

The objectives of this study were to evaluate the accuracy of personalized numerical simulations of the electrical activity in human ventricles by comparing simulated electrocardiograms (ECGs) with real patients’ ECGs and analyzing the sensitivity of the model output to variations in the model parameters. We used standard 12-lead ECGs and up to 224 unipolar body-surface ECGs to record three patients with cardiac resynchronization therapy devices and three patients with focal ventricular tachycardia. Patient-tailored geometrical models of the ventricles, atria, large vessels, liver, and spine were created using computed tomography data. Ten cases of focal ventricular activation were simulated using the bidomain model and the TNNP 2006 cellular model. The population-based values of electrical conductivities and other model parameters were used for accuracy analysis, and their variations were used for sensitivity analysis. The mean correlation coefficient between the simulated and real ECGs varied significantly (from r = 0.29 to r = 0.86) among the simulated cases. A strong mean correlation (r > 0.7) was found in eight of the ten model cases. The accuracy of the ECG simulation varied widely in the same patient depending on the localization of the excitation origin. The sensitivity analysis revealed that variations in the anisotropy ratio, blood conductivity, and cellular apicobasal heterogeneity had the strongest influence on transmembrane potential, while variation in lung conductivity had the greatest influence on body-surface ECGs. Futhermore, the anisotropy ratio predominantly affected the latest activation time and repolarization time dispersion, while the cellular apicobasal heterogeneity mainly affected the dispersion of action potential duration, and variation in lung conductivity mainly led to changes in the amplitudes of ECGs and cardiac electrograms. We also found that the effects of certain parameter variations had specific regional patterns on the cardiac and body surfaces. These observations are useful for further developing personalized cardiac models.

A greedy classifier optimization strategy to assess ion channel blocking activity and pro-arrhythmia in hiPSC-cardiomyocytes

Novel studies conducting cardiac safety assessment using human-induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) are promising but might be limited by their specificity and predictivity. It is often challenging to correctly classify ion channel blockers or to sufficiently predict the risk for Torsade de Pointes (TdP). In this study, we developed a method combining in vitro and in silico experiments to improve machine learning approaches in delivering fast and reliable prediction of drug-induced ion-channel blockade and proarrhythmic behaviour. The algorithm is based on the construction of a dictionary and a greedy optimization, leading to the definition of optimal classifiers. Finally, we present a numerical tool that can accurately predict compound-induced pro-arrhythmic risk and involvement of sodium, calcium and potassium channels, based on hiPSC-CM field potential data.

Predicting drug-induced arrhythmias by multiscale modeling

Drugs often have undesired side effects. In the heart, they can induce lethal arrhythmias such as torsades de pointes. The risk evaluation of a new compound is costly and can take a long time, which often hinders the development of new drugs. Here, we establish a high-resolution, multiscale computational model to quickly assess the cardiac toxicity of new and existing drugs. The input of the model is the drug-specific current block from single cell electrophysiology; the output is the spatio-temporal activation profile and the associated electrocardiogram. We demonstrate the potential of our model for a low-risk drug, ranolazine, and a high-risk drug, quinidine: For ranolazine, our model predicts a prolonged QT interval of 19.4% compared with baseline and a regular sinus rhythm at 60.15 beats per minute. For quinidine, our model predicts a prolonged QT interval of 78.4% and a spontaneous development of torsades de pointes both in the activation profile and in the electrocardiogram. Our model reveals the mechanisms by which electrophysiological abnormalities propagate across the spatio-temporal scales, from specific channel blockage, via altered single cell action potentials and prolonged QT intervals, to the spontaneous emergence of ventricular tachycardia in the form of torsades de pointes. Our model could have important implications for researchers, regulatory agencies, and pharmaceutical companies on rationalizing safe drug development and reducing the time-to-market of new drugs.