Li, Smith, 2009
This CellML model was written by the original model author and is known to be completely correct.
ValidateCellML confirms this model as valid CellML with consistent units
A model of Ca2+ dynamics in mouse ventricular myocytes from the C57BL/6 strain has been developed, based on experimental data at the same temperature and under consistent experimental conditions and protocols.
Parameterisation of the NCX, SERCA and LCC formulations were done sequentially. The model of NCX, based on the Luo and Rudy 1994 formulation, was first fitted from the decay of the caffeine-induced [Ca2+]i transient, which allowed us to infer the flux through SERCA from the decay of the field-stimulated [Ca2+]i transient. Frequency-dependence of the flux through SERCA was accounted for using Hund and Rudy 2002 formulation of Ca2+/calmodulin kinase II regulatory pathway, with parameters fitted to experimentally measured SERCA flux at 0.5, 1 and 3Hz pacing frequencies. The L-type Ca2+ channel model, with slight modifications to the Hinch 2004 formulation, was fitted to data from voltage-clamp experiments. All these component models were then incorporated into the whole-cell AP model of Bondarenko et al., with parameters of RyR adjusted to match the experimentally observed time-to-peak of the [Ca2+]i transient at 3Hz, and parameters for the background Ca2+ current and PMCA adjusted to match experimentally observed diastolic [Ca2+]i levels.
|Schematic diagram of the Bondarenko 2004 model.|
Caffeine - Caffeine-induced Ca2+ transients from 18 individual myocytes
Each file contains the time course of caffeine-induced Ca2+ transients from a myocyte. The traces were exported from IonOptix for the duration of the caffeine application.
Related parameter: Kncx
Fitting process: 1. The decay of individual transients were aligned and averaged. 2. The rate of decay of the average transient, after taking into account of intracellular buffering, was used to calculate Jncx. Jncx was then used to fit Kncx in the Luo and Rudy 1994 model.
3Hz - Field-stimulated Ca2+ transients at 3Hz from 19 individual myocytes
Each file contains the time course of a Ca2+ transient, field stimulated at 3Hz, from a myocyte. The time course was exported from IonOptix for a duration of 1s starting from the first stimulus.
Related parameter: Km, serca, Vserca at 3Hz
Fitting process: 1. The individual transients were averaged. The decay of the average transient, from 70ms after-peak until the end of the period was analysed. 2. Jncx was inferred from previously fitted NCX model, and the difference between the total Ca2+ flux and Jncx was taken as the net SR Ca2+ uptake Jup. 3 A small SR leak (0.0082+ through SERCA Jserca = Jup+ Jleak 4. Jserca was used to fit the Km,serca and Vserca at 3Hz.
1Hz - Field-stimulated Ca2+ transients at 1Hz from 7 individual myocytes
Each file contains the time course of a Ca2+ transient, field stimulated at 1Hz, from a myocyte.
Related parameter: Estimated Vserca at 1Hz
1. The derivation of Jserca at 1Hz is similar to above. 2. Jserca was used to constrain Vserca at 1Hz
0.5Hz - Field-stimulated Ca2+ transients at 0.5Hz from 7 individual myocytes
Each file contains the time course of a Ca2+ transient, field stimulated at 0.5Hz, from a myocyte
Related parameter: Estimated Vserca at 0.5Hz
1. The derivation of Jserca at 0.5Hz is similar to above. 2. Jserca was used to constrain Vserca at 0.5Hz
Ica - Time courses of transmembrane currents from 5 individual myocytes.
Each file contains the time course of experimentally measured L-type Ca2+ current from one myocyte, measured in nA.
Fitting process: 1. The time course of transmembrane currents individual myocytes were averaged for each test potential between -20 and 20mV (-20, -10, 0, 10, 20mV). 2. Averaged currents were converted into current densities by division by cell membrane capacitance (138.66pF). 3. Simulated sum of late K+ currents, from the Bondarenko et al 2004 formulation, for test potentials -20, -10, 0, 10, 20mV were subtracted from the decay phase of each current. 4. The decay of the current at each test potential, after correcting for late K+ currents, was fitted to a single exponential function, giving a time constant of decay tau value. 5. The IV relation and tau values at each test potential were used to fit the model of LCC.
Ica/IV.xls - Experimentally measured current-voltage relation for the L-type Ca2+ current
Average current-voltage curve.
Related parameter: PCaL,
Fitting process: as for Ica