Location: Model of excitation-contraction in uterine myocytes from a pregnant rat @ 3a34f2fe0e1b / Simulation / src / Figs_sim.py

Author:
WeiweiAi <wai484@aucklanduni.ac.nz>
Date:
2021-10-11 15:26:25+13:00
Desc:
Revert ICaT_gating to Tong 2011; I_K1, -V+7.7 in qinf; Add gating variables in individualCurrents
Permanent Source URI:
https://staging.physiomeproject.org/workspace/701/rawfile/3a34f2fe0e1b6f676409d4cfa1c6726b71b1961e/Simulation/src/Figs_sim.py

# To reproduce the data needed for clamping experiments in associated original paper,

import clampExp 
import imp
imp.reload(clampExp)
'''
simFign='FigA1F'
V_clamp=[-60, 120, 120, -50, 50, 10]
paraKeys=['Clamp_parameters/Cai']
paras=[9.43e-5]
outKeys=['outputs/I_CaL']
clampExp.clampExp(simFign, V_clamp,paraKeys,paras,outKeys)

simFign='FigA2F'
V_clamp=[-80, 100, 100, -70, 50, 10]
paraKeys=['Clamp_parameters/Cai']
paras=[9.43e-5]
outKeys=['outputs/I_CaT']
clampExp.clampExp(simFign, V_clamp,paraKeys,paras,outKeys)

simFign='FigA3F'
V_clamp=[-90, 50, 50, -50, 50, 10]
paraKeys=['Clamp_parameters/Nai']
paras=[4]
outKeys=['outputs/I_Na']
clampExp.clampExp(simFign, V_clamp,paraKeys,paras,outKeys)

simFign='FigA4F'
V_clamp=[-80, 8000, 8000, -60, 60, 10]
paraKeys=['Clamp_parameters/Ki']
paras=[140]
outKeys=['outputs/I_K1']
clampExp.clampExp(simFign, V_clamp,paraKeys,paras,outKeys)
'''
simFign='FigA5F'
V_clamp=[-80, 8000, 8000, -60, 60, 10]
paraKeys=['Clamp_parameters/Ki']
paras=[140]
outKeys=['outputs/I_K2']
clampExp.clampExp(simFign, V_clamp,paraKeys,paras,outKeys)

simFign='FigA6F'
V_clamp=[-80, 150, 150, -80, 60, 10]
paraKeys=['Clamp_parameters/Ki']
paras=[140]
outKeys=['outputs/I_KA']
clampExp.clampExp(simFign, V_clamp,paraKeys,paras,outKeys)

simFign='FigA7F'
V_clamp=[-10, 3000, 3000, -20, 130, 10]
paraKeys=['Clamp_parameters/Nai','Clamp_parameters/Ki']
paras=[4,140]
outKeys=['outputs/I_h']
clampExp.clampExp(simFign, V_clamp,paraKeys,paras,outKeys)

simFign='FigA8F'
V_clamp=[0, 100, 100, -60, 40, 10]
paraKeys=['Clamp_parameters/Cai']
paras=[9.43e-5]
outKeys=['outputs/I_Cl']
clampExp.clampExp(simFign, V_clamp,paraKeys,paras,outKeys)