A global sensitivity tool for cardiac cell modeling: Application to ionic current balance and hypertrophic signaling

Cardiovascular diseases are the major cause of death in the developed countries. Identifying key cellular processes involved in generation of the electrical signal and in regulation of signal transduction pathways is essential for unraveling the underlying mechanisms of heart rhythm behavior. Computational cardiac models provide important insights into cardiovascular function and disease. Sensitivity analysis presents a key tool for exploring the large parameter space of such models, in order to determine the key factors determining and controlling the underlying physiological processes. We developed a new global sensitivity analysis tool which implements the Morris method, a global sensitivity screening algorithm, onto a Nimrod platform, which is a distributed resources software toolkit. The newly developed tool has been validated using the model of IP3-calcineurin signal transduction pathway model which has 30 parameters. The key driving factors of the IP3 transient behaviour have been calculated and confirmed to agree with previously published data. We next demonstrated the use of this method as an assessment tool for characterizing the structure of cardiac ionic models. In three latest human ventricular myocyte models, we examined the contribution of transmembrane currents to the shape of the electrical signal (i.e. on the action potential duration). The resulting profiles of the ionic current balance demonstrated the highly nonlinear nature of cardiac ionic models and identified key players in different models. Such profiling suggests new avenues for development of methodologies to predict drug action effects in cardiac cells.

[1]  Christopher R. Myers,et al.  Universally Sloppy Parameter Sensitivities in Systems Biology Models , 2007, PLoS Comput. Biol..

[2]  P. Hunter,et al.  Modeling hypertrophic IP3 transients in the cardiac myocyte. , 2007, Biophysical journal.

[3]  Max D. Morris,et al.  Factorial sampling plans for preliminary computational experiments , 1991 .

[4]  David Abramson,et al.  Scheduling large parametric modelling experiments on a distributed meta-computer , 1997 .

[5]  Denis Noble,et al.  Markov models for ion channels: versatility versus identifiability and speed , 2009, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[6]  Socrates Dokos,et al.  Parameter estimation in cardiac ionic models. , 2004, Progress in biophysics and molecular biology.

[7]  K. T. ten Tusscher,et al.  Alternans and spiral breakup in a human ventricular tissue model. , 2006, American journal of physiology. Heart and circulatory physiology.

[8]  Blanca Rodríguez,et al.  Impact of ionic current variability on human ventricular cellular electrophysiology. , 2009, American journal of physiology. Heart and circulatory physiology.

[9]  E. Sobie Parameter sensitivity analysis in electrophysiological models using multivariable regression. , 2009, Biophysical journal.

[10]  D. Bers,et al.  A novel computational model of the human ventricular action potential and Ca transient. , 2010, Journal of Molecular and Cellular Cardiology.

[11]  D. Noble A modification of the Hodgkin—Huxley equations applicable to Purkinje fibre action and pacemaker potentials , 1962, The Journal of physiology.

[12]  D. Noble,et al.  Systems biology and the virtual physiological human , 2009, Molecular systems biology.

[13]  Denis Noble,et al.  Contributions of HERG K+ current to repolarization of the human ventricular action potential. , 2008, Progress in biophysics and molecular biology.