Cellular and ionic causes of variability in the electrophysiological activity of hearts from individuals of the same species are unknown. However, improved understanding of this variability is key to enable prediction of the response of specific hearts to disease and therapies. Limitations of current mathematical modeling and experimental techniques hamper our ability to provide insight into variability. Here, we describe a methodology to unravel the ionic determinants of intersubject variability exhibited in experimental recordings, based on the construction and calibration of populations of models. We illustrate the methodology through its application to rabbit Purkinje preparations, because of their importance in arrhythmias and safety pharmacology assessment. We consider a set of equations describing the biophysical processes underlying rabbit Purkinje electrophysiology, and we construct a population of over 10,000 models by randomly assigning specific parameter values corresponding to ionic current conductances and kinetics. We calibrate the model population by closely comparing simulation output and experimental recordings at three pacing frequencies. We show that 213 of the 10,000 candidate models are fully consistent with the experimental dataset. Ionic properties in the 213 models cover a wide range of values, including differences up to ±100% in several conductances. Partial correlation analysis shows that particular combinations of ionic properties determine the precise shape, amplitude, and rate dependence of specific action potentials. Finally, we demonstrate that the population of models calibrated using data obtained under physiological conditions quantitatively predicts the action potential duration prolongation caused by exposure to four concentrations of the potassium channel blocker dofetilide.cardiac electrophysiology | computational biology | mathematical modeling | systems biology | drug B iological variability is exhibited at all levels in all organs of living organisms. It manifests itself as differences in physiological function between individuals of the same species and often more drastically by significant differences in the outcome of their exposure to pathological conditions. Thus, healthy cardiac cells of the same species and location exhibit a qualitatively similar response to a stimulus, i.e., the action potential (AP). However, significant quantitative intersubject differences exist in AP morphology and duration, which may explain the different individual response of each of the cells (and patients) to disease or drug action.The variability underlying the physiological and pathological responses of different individuals has often been ignored in experimental and theoretical research, ultimately hampering the extrapolation of results to a population level. Experimentalists often average the results obtained in different preparations to reduce experimental error, therefore also averaging out the effects of intersubject variation and resulting in an important loss of information. T...

Early prediction of cardiotoxicity is critical for drug development. Current animal models raise ethical and translational questions, and have limited accuracy in clinical risk prediction. Human-based computer models constitute a fast, cheap and potentially effective alternative to experimental assays, also facilitating translation to human. Key challenges include consideration of inter-cellular variability in drug responses and integration of computational and experimental methods in safety pharmacology. Our aim is to evaluate the ability of in silico drug trials in populations of human action potential (AP) models to predict clinical risk of drug-induced arrhythmias based on ion channel information, and to compare simulation results against experimental assays commonly used for drug testing. A control population of 1,213 human ventricular AP models in agreement with experimental recordings was constructed. In silico drug trials were performed for 62 reference compounds at multiple concentrations, using pore-block drug models (IC50/Hill coefficient). Drug-induced changes in AP biomarkers were quantified, together with occurrence of repolarization/depolarization abnormalities. Simulation results were used to predict clinical risk based on reports of Torsade de Pointes arrhythmias, and further evaluated in a subset of compounds through comparison with electrocardiograms from rabbit wedge preparations and Ca2+-transient recordings in human induced pluripotent stem cell-derived cardiomyocytes (hiPS-CMs). Drug-induced changes in silico vary in magnitude depending on the specific ionic profile of each model in the population, thus allowing to identify cell sub-populations at higher risk of developing abnormal AP phenotypes. Models with low repolarization reserve (increased Ca2+/late Na+ currents and Na+/Ca2+-exchanger, reduced Na+/K+-pump) are highly vulnerable to drug-induced repolarization abnormalities, while those with reduced inward current density (fast/late Na+ and Ca2+ currents) exhibit high susceptibility to depolarization abnormalities. Repolarization abnormalities in silico predict clinical risk for all compounds with 89% accuracy. Drug-induced changes in biomarkers are in overall agreement across different assays: in silico AP duration changes reflect the ones observed in rabbit QT interval and hiPS-CMs Ca2+-transient, and simulated upstroke velocity captures variations in rabbit QRS complex. Our results demonstrate that human in silico drug trials constitute a powerful methodology for prediction of clinical pro-arrhythmic cardiotoxicity, ready for integration in the existing drug safety assessment pipelines.

beta-Adrenergic challenges trigger TdP in a reproducible manner in this model of drug-induced LQT1. Paradoxical prolongation and increased temporal and spatial dispersion of repolarization precipitate TdP. Incremental LV systolic aftercontractions precede TdP, suggesting abnormal cellular Ca(2+) handling contributes to the arrhythmogenic mechanism.

No abstract