- Author:
- WeiweiAi <wai484@aucklanduni.ac.nz>
- Date:
- 2022-06-13 12:35:09+12:00
- Desc:
- Update the documentation
- Permanent Source URI:
- https://models.physiomeproject.org/workspace/83f/rawfile/309a2cd1708eb79e1e9f8b928b41c9d36bcb94fc/Doc/Single_KCl_stim_experiment.rst
Single-cell response to KCl stimulation
-----------------------------------------
In the `single-cell response to KCl stimulation <Experiments/PASMC_test.cellml/view>`_, the `isolated cell model <../Components/PASMC.cellml>`_ is configured and parameterised with an applied `concentration clamp protocol <../cellLib/Protocols/X_clamp_protocol.cellml>`_.
You can change the parameters of stimulation in the component ``clamp_para``. You can also bypass the parameters in the component ``Para`` and define your own parameters.
The simulation experiment can be obtained by loading the `corresponding SED-ML document <../Experiments/PASMC_test.sedml>`_ into OpenCOR and executing the simulation. We use the Python script `simFig4 <../Simulation/src/simFig4.py>`_ to perform the simulation for the following figures.
We set the intracellular :math:`K^{+}` concentration :math:`K_{in}=198.7874 mM` and the initial value of the open state probability of :math:`K^{+}` channels :math:`w_i=0.5` during the simulation.
Dynamics when :math:`K_{out}=5.9 mM`
------------------------------------------------------
Figure [#]_ shows the oscillating dynamics of :math:`Ca^{2+}` concentrations and the cell membrane potential, when the extracellular :math:`K^{+}` concentration :math:`K_{out}=5.9 mM`.
.. [#]
.. figure:: ../Simulation/simFig4_1.png
:width: 75%
:align: center
:alt: simFig4_1
Single-cell response to KCl stimulation
------------------------------------------------------
We set the extracellular :math:`K^{+}` concentration :math:`K_{out}=5.9 mM` before :math:`t=10s` and :math:`K_{out}` steps up to 81.2 mM at :math:`t=10s`. Since the reversal potential for :math:`K^{+}` current is associated with :math:`K_{out}`, the reversal potential :math:`v_K` in original equation (13) is also changed at :math:`t=10s`.
Figure [#]_ shows the simulation result where there is an initial surge of :math:`c_i` when we use the initial values from the cited model (Parthimos et al. 1999).
.. [#]
.. figure:: ../Simulation/simFig4_0.png
:width: 75%
:align: center
:alt: simFig4_0
We tried different initial values and Figure [#]_ shows a few traces compared to the original figure. Although the closest orange trace has some overlap with the original data, there is still discrepancy between them. We think this might have something to do with the intracellular :math:`K^{+}` concentration value :math:`K_{in}=198.7874 mM` or the holding extracellular :math:`K^{+}` concentration :math:`K_{out}=5.9 mM` before :math:`t=10s`.
.. [#]
.. figure:: ../Simulation/simFig4_ci.png
:width: 75%
:align: center
:alt: simFig4_ci
Simulated intracellular :math:`Ca^{2+}` concentration :math:`c_i` response to KCl stimulation.
