- Author:
- WeiweiAi <wai484@aucklanduni.ac.nz>
- Date:
- 2022-07-28 12:05:55+12:00
- Desc:
- Fixed the link
- Permanent Source URI:
- http://models.cellml.org/workspace/8af/rawfile/1b1a07491a4d102536a3f7e25957ca2810d6ed76/Doc/MWC_10_experiment.rst
MWC-10 model experiment
---------------------------------
In the `MWC_10_experiment <Experiments/MWC_10_test.cellml/view>`_, the `MWC-10 model <../Components/MWC_10.cellml>`_ is configured and parameterised with an applied `voltage clamp protocol <../cellLib/Protocols/dPulse_protocol_ms.cellml>`_.
You can change the parameters of stimulation in the component ``clamp_para``. You can also bypass the parameters in the components ``free_para_10`` and ``Para``, and define your own parameters.
Steady states
----------------
We have used two approaches to get the steady states of system:
- 1. Apply the double pulse protocol and record the simulation when the system reaches plateau (after 200 ms);
- 2. Use Equation (6) and (7) in the primary paper, denoted by :math:`_{_ss}` in the following figures.
The results from these two methods are identical in our simulation.
Figure [#]_ shows the voltage dependencies and transfer function.
.. [#]
.. figure:: ../Simulation/simFig4.png
:width: 75%
:align: center
:alt: simFig4
Figure [#]_ shows the voltage dependencies and transfer function under reference conditions.
.. [#]
.. figure:: ../Simulation/simFig6_r.png
:width: 75%
:align: center
:alt: simFig6_r
Figure [#]_ shows the voltage dependencies and transfer function under perchlorate conditions, where the simulation for Fiber 911 case does not match well.
.. [#]
.. figure:: ../Simulation/simFig6_p.png
:width: 75%
:align: center
:alt: simFig6_r
Kinetic transients
-------------------
Figure [#]_ shows charge movement transients (simulation data are divided by 5) in reference simulation with different test voltages.
.. [#]
.. figure:: ../Simulation/simFig9_r.png
:width: 75%
:align: center
:alt: simFig9_r
Figure [#]_ shows charge movement transients (simulation data are divided by 5) in perchlorate simulation with different test voltages.
.. [#]
.. figure:: ../Simulation/simFig9_p.png
:width: 75%
:align: center
:alt: simFig9_p
For primary Fig. 11, we have tried different settings in Table 1 and two sets of test voltages, but we couldn't find a perfect match.
Figure [#]_ shows the open probability and charge movement current when 827 parameters are used.
.. [#]
.. figure:: ../Simulation/simFig11_827.png
:width: 75%
:align: center
:alt: simFig11_827
Figure [#]_ shows the open probability and charge movement current when 908 parameters are used.
.. [#]
.. figure:: ../Simulation/simFig11_908.png
:width: 75%
:align: center
:alt: simFig11_908
Figure [#]_ shows the open probability and charge movement current when 906 parameters are used.
.. [#]
.. figure:: ../Simulation/simFig11_906.png
:width: 75%
:align: center
:alt: simFig11_906
Figure [#]_ shows the open probability and charge movement current when 907 parameters are used.
.. [#]
.. figure:: ../Simulation/simFig11_907.png
:width: 75%
:align: center
:alt: simFig11_907
The following two figures only have simulation results, which show that the qualitative difference between reference and perchlorate conditions agrees with the original figure, that is, perchlorate conditions slow down the transitions.
Figure [#]_ shows the open probability and charge movement current under reference and perchlorate conditions when 827 parameters are used.
.. [#]
.. figure:: ../Simulation/simFig11_827-comp.png
:width: 75%
:align: center
:alt: simFig11_827-comp
Figure [#]_ shows the open probability and charge movement current under reference and perchlorate conditions when 908 parameters are used.
.. [#]
.. figure:: ../Simulation/simFig11_908-comp.png
:width: 75%
:align: center
:alt: simFig11_908-comp
