Corrias Alberto albertocorrias@nus.edu.sg National University of Singapore Division of BioEngineering 2007-09-10 National University of Singapore Corrias and Buist's 2007 quantitative model of gastric smooth muscle cellular electrophysiology. This is a CellML description of all the ion channels and intracellular processes contained in the Corrias and Buist 2007 quantitative model of gastric smooth muscle cellular activation. This encoding of the model includes a stimulus current derived from a psuedo ICC action potential. Alberto Corrias keyword electrophysiology gastric smooth muscle 17486452 Corrias Alberto Buist Martin L A quantitative model of gastric smooth muscle cellular activation 2007 Annals of Biomedical Engineering 35 1595 1607 SM with psuedo ICC stimulus bdf newton dense 1 10.0 0.0 60.0e3 Time (ms) Membrane potential (mV) #ffffff #ff0000 1 Time (ms) I_stim (pA) #ffffff #ff0000 1 Time (ms) [Ca]i (mM) #ffffff #ff0000 1 19.5e3 32.0e3 Time (ms) Current (pA) #ffffff #ff0000 1 19.5e3 32.0e3 Time (ms) Current (pA) #ffffff ICaL #ff0000 1 19.5e3 32.0e3 ILVA #0000ff 1 19.5e3 32.0e3 Time (ms) Current (pA) #ffffff #ff0000 1 19.5e3 32.0e3 Time (ms) Current (pA) #ffffff #ff0000 1 19.5e3 32.0e3 Time (ms) Current (pA) #ffffff #ff0000 1 19.5e3 32.0e3 Tidying up the electrical stimulus protocol and adding simulation and graphing metadata for my CellMLSimulator software. David Nickerson 2007-09-20 The value of Q10 for Potassium channels (Q10K) was corrected from the value in the published paper (1.5). The correct value to be used is 1.365. Alberto Corrias 2007-09-21 $\mathrm{FoRT}=\frac{F}{RT}$ $\mathrm{RToF}=\frac{RT}{F}$ Here we define a periodic stimulus based on a psuedo ICC action potential. This stimulus protocol is used to test the model under a realistic stimulus protocol, but ideally such a stimulus would originate directly from an ICC cellular electrophysiology model. This stimulus protocol simply applies the psuedo ICC stimulus a fixed number of times. $\mathrm{local\_time}=\mathrm{time}-\mathrm{stim\_start}+\mathrm{t\_ICCpeak}$ $\frac{d \mathrm{d\_Ltype\_SM}}{d \mathrm{time}}=\frac{\mathrm{d\_inf\_Ltype\_SM}-\mathrm{d\_Ltype\_SM}}{\mathrm{tau\_d\_Ltype\_SM}}$ $\frac{d \mathrm{f\_Ltype\_SM}}{d \mathrm{time}}=\frac{\mathrm{f\_inf\_Ltype\_SM}-\mathrm{f\_Ltype\_SM}}{\mathrm{tau\_f\_Ltype\_SM}}$ $\frac{d \mathrm{f\_ca\_Ltype\_SM}}{d \mathrm{time}}=\frac{\mathrm{f\_ca\_inf\_Ltype\_SM}-\mathrm{f\_ca\_Ltype\_SM}}{\mathrm{tau\_f\_ca\_Ltype\_SM}}$ $\mathrm{I\_Ltype\_SM}=\mathrm{G\_max\_Ltype}\mathrm{f\_Ltype\_SM}\mathrm{d\_Ltype\_SM}\mathrm{f\_ca\_Ltype\_SM}(\mathrm{Vm\_SM}-\mathrm{E\_Ca})$ $\frac{d \mathrm{d\_LVA\_SM}}{d \mathrm{time}}=\frac{\mathrm{d\_inf\_LVA\_SM}-\mathrm{d\_LVA\_SM}}{\mathrm{tau\_d\_LVA\_SM}}$ $\frac{d \mathrm{f\_LVA\_SM}}{d \mathrm{time}}=\frac{\mathrm{f\_inf\_LVA\_SM}-\mathrm{f\_LVA\_SM}}{\mathrm{tau\_f\_LVA\_SM}}$ $\mathrm{I\_LVA\_SM}=\mathrm{G\_max\_LVA}\mathrm{f\_LVA\_SM}\mathrm{d\_LVA\_SM}(\mathrm{Vm\_SM}-\mathrm{E\_Ca})$ $\mathrm{E\_K}=\mathrm{RToF}\ln \left(\frac{\mathrm{K\_o}}{\mathrm{K\_i}}\right)$ $\mathrm{I\_BK\_SM}=(\mathrm{G\_max\_BK}+\mathrm{T\_correction\_BK})\mathrm{d\_BK\_SM}(\mathrm{Vm\_SM}-\mathrm{E\_K})$ $\mathrm{E\_K}=\mathrm{RToF}\ln \left(\frac{\mathrm{K\_o}}{\mathrm{K\_i}}\right)$ $\mathrm{I\_bk\_SM}=\mathrm{G\_max\_bk}(\mathrm{Vm\_SM}-\mathrm{E\_K})$ $\frac{d \mathrm{xr1\_SM}}{d \mathrm{time}}=\frac{\mathrm{xr1\_inf\_SM}-\mathrm{xr1\_SM}}{\mathrm{tau\_xr1\_SM}}$ $\frac{d \mathrm{xr2\_SM}}{d \mathrm{time}}=\frac{\mathrm{xr2\_inf\_SM}-\mathrm{xr2\_SM}}{\mathrm{tau\_xr2\_SM}}$ $\mathrm{E\_K}=\mathrm{RToF}\ln \left(\frac{\mathrm{K\_o}}{\mathrm{K\_i}}\right)$ $\mathrm{I\_kr\_SM}=\mathrm{G\_max\_kr\_SM}\mathrm{xr1\_SM}\mathrm{xr2\_SM}(\mathrm{Vm\_SM}-\mathrm{E\_K})$ $\frac{d \mathrm{m\_Na\_SM}}{d \mathrm{time}}=\frac{\mathrm{m\_inf\_Na}-\mathrm{m\_Na\_SM}}{\mathrm{tau\_m\_Na}}$ $\frac{d \mathrm{h\_Na\_SM}}{d \mathrm{time}}=\frac{\mathrm{h\_inf\_Na}-\mathrm{h\_Na\_SM}}{\mathrm{tau\_h\_Na}}$ $\mathrm{E\_Na}=\mathrm{RToF}\ln \left(\frac{\mathrm{Na\_o}}{\mathrm{Na\_i}}\right)$ $\mathrm{I\_Na\_SM}=\mathrm{G\_max\_Na\_SM}\mathrm{h\_Na\_SM}\mathrm{m\_Na\_SM}(\mathrm{Vm\_SM}-\mathrm{E\_Na})$ $\frac{d \mathrm{xa1\_SM}}{d \mathrm{time}}=\frac{\mathrm{xa1\_inf\_SM}-\mathrm{xa1\_SM}}{\mathrm{tau\_xa1\_SM}}$ $\frac{d \mathrm{xa2\_SM}}{d \mathrm{time}}=\frac{\mathrm{xa2\_inf\_SM}-\mathrm{xa2\_SM}}{\mathrm{tau\_xa2\_SM}}$ $\mathrm{E\_K}=\mathrm{RToF}\ln \left(\frac{\mathrm{K\_o}}{\mathrm{K\_i}}\right)$ $\mathrm{I\_ka\_SM}=\mathrm{G\_max\_ka\_SM}\mathrm{xa1\_SM}\mathrm{xa2\_SM}(\mathrm{Vm\_SM}-\mathrm{E\_K})$ $\frac{d \mathrm{m\_NSCC\_SM}}{d \mathrm{time}}=\frac{\mathrm{m\_inf\_NSCC\_SM}-\mathrm{m\_NSCC\_SM}}{\mathrm{tau\_m\_NSCC\_SM}}$ $\mathrm{I\_NSCC\_SM}=\mathrm{G\_max\_NSCC\_SM}\mathrm{m\_NSCC\_SM}\mathrm{f\_ca\_NSCC\_SM}\mathrm{rach\_NSCC\_SM}(\mathrm{Vm\_SM}-\mathrm{E\_NSCC})$