Simulation of the undiseased human cardiac ventricular action potential: Model formulation and experimental validation

Thomas O'Hara, László Virág, András Varró, Yoram Rudy

Research output: Contribution to journalArticle

553 Citations (Scopus)


Cellular electrophysiology experiments, important for understanding cardiac arrhythmia mechanisms, are usually performed with channels expressed in non myocytes, or with non-human myocytes. Differences between cell types and species affect results. Thus, an accurate model for the undiseased human ventricular action potential (AP) which reproduces a broad range of physiological behaviors is needed. Such a model requires extensive experimental data, but essential elements have been unavailable. Here, we develop a human ventricular AP model using new undiseased human ventricular data: Ca2+ versus voltage dependent inactivation of L-type Ca2+ current (ICaL); kinetics for the transient outward, rapid delayed rectifier (IKr), Na+/Ca2+ exchange (INaCa), and inward rectifier currents; AP recordings at all physiological cycle lengths; and rate dependence and restitution of AP duration (APD) with and without a variety of specific channel blockers. Simulated APs reproduced the experimental AP morphology, APD rate dependence, and restitution. Using undiseased human mRNA and protein data, models for different transmural cell types were developed. Experiments for rate dependence of Ca2+ (including peak and decay) and intracellular sodium ([Na+]i) in undiseased human myocytes were quantitatively reproduced by the model. Early afterdepolarizations were induced by IKr block during slow pacing, and AP and Ca2+ alternans appeared at rates >200 bpm, as observed in the nonfailing human ventricle. Ca2+/calmodulin-dependent protein kinase II (CaMK) modulated rate dependence of Ca2+ cycling. INaCa linked Ca2+ alternation to AP alternans. CaMK suppression or SERCA upregulation eliminated alternans. Steady state APD rate dependence was caused primarily by changes in [Na+]i, via its modulation of the electrogenic Na+/K+ ATPase current. At fast pacing rates, late Na+ current and ICaL were also contributors. APD shortening during restitution was primarily dependent on reduced late Na+ and ICaL currents due to inactivation at short diastolic intervals, with additional contribution from elevated IKr due to incomplete deactivation.

Original languageEnglish
Article numbere1002061
JournalPLoS computational biology
Issue number5
Publication statusPublished - May 1 2011

ASJC Scopus subject areas

  • Ecology, Evolution, Behavior and Systematics
  • Modelling and Simulation
  • Ecology
  • Molecular Biology
  • Genetics
  • Cellular and Molecular Neuroscience
  • Computational Theory and Mathematics

Fingerprint Dive into the research topics of 'Simulation of the undiseased human cardiac ventricular action potential: Model formulation and experimental validation'. Together they form a unique fingerprint.

  • Cite this