clancy1

The SBML for this model was obtained from the BioModels database (BioModels ID: BIOMD0000000121) Biomodels notes: This model is according to the paper Cellular consequences of HEGR mutations in the long QT syndrome: precursors to sudden cardiac death. The author used Markovian model of cardiac Ikr in the paper. Figure4B in the paper has been reproduced using CellDesigner3.5.1. The cell is depolarized to the indicated test potential for 250ms (from 50ms to 300ms) from a holding potential of -40mV and then repolarized to -40mV. Change the value for vtest from -30,-20,-10,0,10,20,30,40 for each simulation in order to produce the different cureve in the paper. JWS Online curation: This model was curated by reproducing the figures as described in the BioModels Notes. No additional changes were made.

None

None

None

None

None

None

Cellular consequences of HERG mutations in the long QT syndrome: precursors to sudden cardiac death.

  • Colleen E Clancy
  • Yoram Rudy
Cardiovasc. Res. 2001; 50 (2): 301-313
Abstract
BACKGROUND: A variety of mutations in HERG, the major subunit of the rapidly activating component of the cardiac delayed rectifier I(Kr), have been found to underlie the congenital Long-QT syndrome, LQT2. LQT2 may give rise to severe arrhythmogenic phenotypes leading to sudden cardiac death.
OBJECTIVE: We attempt to elucidate the mechanisms by which heterogeneous LQT2 genotypes can lead to prolongation of the action potential duration (APD) and consequently the QT interval on the ECG.
METHODS: We develop Markovian models of wild-type (WT) and mutant I(Kr) channels and incorporate these models into a comprehensive model of the cardiac ventricular cell.
RESULTS: Using this virtual transgenic cell model, we describe the effects of HERG mutations on the cardiac ventricular action potential (AP) and provide insight into the mechanism by which each defect results in a net loss of repolarizing current and prolongation of APD.
CONCLUSIONS: This study demonstrates which mutations can prolong APD sufficiently to generate early afterdepolarizations (EADs), which may trigger life-threatening arrhythmias. The severity of the phenotype is shown to depend on the specific kinetic changes and how they affect I(Kr) during the time course of the action potential. Clarifying how defects in HERG can lead to impaired cellular electrophysiology can improve our understanding of the link between channel structure and cellular function.

Unit definitions have no effect on the numerical analysis of the model. It remains the responsibility of the modeler to ensure the internal numerical consistency of the model. If units are provided, however, the consistency of the model units will be checked.

Name Definition
Id Name Spatial dimensions Size
cell cell 1.0
Id Name Initial quantity Compartment Fixed
c1 0.0 cell (cell)
c2 0.0 cell (cell)
c3 1.0 cell (cell)
i 0.0 cell (cell)
o 0.06 cell (cell)

Initial assignments are expressions that are evaluated at time=0. It is not recommended to create initial assignments for all model entities. Restrict the use of initial assignments to cases where a value is expressed in terms of values or sizes of other model entities. Note that it is not permitted to have both an initial assignment and an assignment rule for a single model entity.

Definition
Id Name Objective coefficient Reaction Equation and Kinetic Law Flux bounds
reaction0000001 c3 = c2

cell*(-(b*c2) + a*c3)
reaction0000002 c2 = c1

cell*(-(bin*c1) + ain*c2)
reaction0000003 c1 = o

cell*(aa*c1 - bb*o)
reaction0000004 o = i

cell*(-(ai*i) + bi*o)
reaction0000005 c1 = i

cell*(aa*c1 - u*i)

Global parameters

Id Value
F 96.485
Gk 0.0
R 8.314
Temp 310.0
a 0.0
aa 0.0
ai 0.0
ain 2.172
b 0.0
bb 0.0
bi 0.0
bin 1.077
ik 0.0
ki 140.0
ko 5.4
u 0.0
v -40.0
vhold -40.0
vk 0.0
vtest 0.0

Local parameters

Id Value Reaction

Assignment rules

Definition
aa = 0.0655*exp(0.05547153*(-36.0 + v))
ik = Gk*(v - vk)*o
a = 0.0555*exp(0.05547153*(-12.0 + v))
b = 0.002357*exp(-0.036588*v)
Gk = 0.013500000000000002*ko^0.59
bi = 1.030072027074169*exp(0.000942*v)*(ko^(-1.0))^0.3
u = (ai*bb)/bi
bb = 0.0029357000000000003*exp(-0.02158*v)
vk = (R*Temp*log(ko/ki))/F
ai = 1.9755*exp(-0.02352*(25.0 + v))/ko
v = piecewise(vtest,50.0<time<=300.0,vhold)

Rate rules

Definition

Algebraic rules

Definition
Trigger Assignments