fig3.py

Calcium current model optimization

Steady-state inactivation and activation iPSC-CM experimental data, with optimized dataset-specific models for different dataset (Fig. 3 A). I-V curves for ICaL. Calcium-dependent gating model formulation retained from ten Tusscher (2004) adult cardiomyocyte model with parameter optimization to fit whole cell iPSC-CM outputs (Fig. 3 B) Time constants were extracted from current recordings by Ma et al. (2011) (same as sodium current). The activation and inactivation time constant data used for model optimization, and the resulting models, are shown in Fig. 3 C and D.

Current_Ica.cellml is the main CellML files which has all the formulation for calcium current, Its associated Sedml file contains all the simulation settings.

All the CellML files and SED-ML files related to this channel need to be downloaded in a same folder (Current_Ica, gating_Ica, parameter_Ica, parameter, unit) as well as python script (fig3.py). In the python script, required Sedml file (Current_Ica.sedml) is loaded into the script and by running the code following figure is reproduced. fig3.py is used to generate the simulation and reproduces the graph shown in Figure 4 in the original study. In order to reproduce Figure 4, once all the files are downloaded to the same folder, execute the following script from the command line (command prompt):

cd [PathToThisFile]

[PathToOpenCOR]/pythonshell fig3.py

A, L-type calcium current (ICaL) steady-state inactivation and activation curves with dataset-specific model optimized to experimental data. The L-type calcium model used in the baseline cellular model is shown in black. Coloured symbols represent experimental iPSC-CM data from Ma et al. (2011) ,`Veerman et al. (2016)`_ and Es-Salah-Lamoureux et al. (2016). B, I-V curves for ICaL, C, time constants of ICaL activation gate. D, time constants of ICaL inactivation gate.