fig4-new.py

Calcium current model optimization

Steady-state inactivation and activation iPSC-CM experimental data, as well as optimized dataset-specific models for different dataset (Fig. 4 A). 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. 4 C and D.

act_inact.cellml is the main CellML files which show the probability of calcium channel being open or close. Its associated Sedml file contains all the simulation settings. All the CellML files and Sedml files need to be download in a same folder (act_inact, gating, parameter, unit) as well as python script (fig4-new.py). In the python script, required Sedml file is loaded into the script and by running the code following figure is reproduced. fig4-new.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 fig4-new.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). C, time constants of ICaL activation gate. D, time constants of ICaL inactivation gate.