fig6-new.py

About this model

Original publication:
 Kernik et al. (2019): "A computational model of induced pluripotent stem-cell derived cardiomyocytes incorporating experimental variability from multiple data sources" J Physiol. 2019 Sep 1; 597(17): 4533-4564.
DOI:https://dx.doi.org/10.1113%2FJP277724

Figure 06

Transient outward potassium current (Ito) model optimization

The voltage dependent activation and inactivation gating variables were modeled here. experimental iPSC‐CM data collected from Ma et al. (2011), Cordeiro et al. (2013) and Veerman et al. (2016) in order to optimize data specific model. Since non of the datasets had both steady-state activation and inactivation data therefore, Inactivation time constant parameters of Veerman et al. (2016) model were optimized with time constant activation from Ma et al. (2011). Cordeiro et al. (2013) model parameters were optimized with steady-state activation parameters from Ma et al. (2011) and steady-state inactivation data from Cordeiro et al. (2013) were used to optimize the model of Ma et al. (2011). All the results are presented in Fig. 6. D. Time constant activation parameters were optimized to ten Tusscher (2004) model which is shown in Fig. 6 C.

act_inact.cellml is the main CellML file which shows the probability of transient outward potassium channel being open orclose. Its associated SED-ML file contains all the simulation settings. All the CellML files and SED-ML files need to be download in a same folder (act_inact, gating, parameter, unit) as well as python script (fig6-new.py). In the python script, required SED-ML file is loaded into the script and by running the code following figure is reproduced. fig6-new.py is used to generate the simulation and reproduces the graph shown in Figure 6 in the original study. In order to reproduce Figure 6, once all the files are downloaded to the same folder, execute the following script from the command line (command prompt):

cd [PathToThisFile]

[PathToOpenCOR]/pythonshell fig6-new.py

Transient outward potassium current model optimization

A, optimized activation with dataset-specific model fits. Different colour represent experimental iPSC-CM data from multiple laboratories. C, Time constant of Ito activation gate. D, Time constant of Ito inactivation gate.