Location: FCU_cardiacAP_composite @ f4d1e338a31b / exposure / exposure_frontpage.rst

Author:
Shelley Fong <sfon036@UoA.auckland.ac.nz>
Date:
2022-05-19 14:40:40+12:00
Desc:
Removing Ca fluxes between cytosol and SR from this FCU
Permanent Source URI:
http://models.cellml.org/workspace/839/rawfile/f4d1e338a31bc1ed2afc4a7f02e12a36c8a18dc7/exposure/exposure_frontpage.rst

About this Functional Cell Unit (FCU)
=======================================

This is a workspace containing bond graph models of ion channels and transporters which were identified to play important roles in the generation and maintenance of an action potential within a cardiac myocyte. Simply put, an ion channel transports ion(s) from one cellular compartment to another. 

One can select components from what is available and build their own composite CellML file. 

An existing CellML file is provided, which contains the channels used by `Kernik et. al. <https://models.physiomeproject.org/e/817>`_, except for the 'funny' current channel and the 'transient outward' channel. 

The list of channels in this workspace is subject to change as literature demands it.

For a given cell model:

    **INPUTS:** 
        - Stimulus current (I\ :sub:`stim`\)
        
    **OUTPUTS:** 
        - Cytosolic and extracellular ionic concentrations (Na, K, Ca) 
        - Membrane potential (V\ :sub:`m`\)
        - Currents through each channel (I)
        
    **CHANNELS AND TRANSPORTERS:** 
        .. csv-table:: 
           :header: "Abbreviation", "Name"
           :widths: 3, 10
           
           "Na\ :sub:`B`\", "Background sodium channel"
           "fast-Na", "Fast sodium channel"
           "K1", "Inward rectifying potassium channel"
           "K\ :sub:`r`\", "Rapid delayed rectifying potassium channel"
           "K\ :sub:`s`\", "Slow delayed rectifying potassium channel"
           "Ca\ :sub:`B`\", "Background calcium channel"
           "Ca\ :sub:`leak`\", "Leakage calcium channel"
           "LCC", "L-type calcium transporter"
           "TCC", "T-type calcium transporter"
           "pCa", "Sarcolemmal calcium pump"
           "NaK", "Sodium-potassium pump"
           "NCX", "Sodium-calcium exchanger"
           "RyR", "Ryanodine receptor"
           "SERCA", "Sarcoplasmic-endoplasmic reticulum calcium pump"
           "CMDN", "Calmodulin buffer for calcium"
           "TRPN", "Troponin buffer for calcium"

Model status
=============

For individual channels, their current CellML implementation run in OpenCOR. Composite models do not run due to solver failure, including the model reproducing Kernik et. al. An exported MATLAB script is provided as an alternative.


Model overview
===================
The components of this workspace are based on existing kinetic models taken from `Kernik et. al. (2019) <https://doi.org/10.1113/JP278739>`_ or from `Luo Rudy (1994) <http://www.ncbi.nlm.nih.gov/pubmed/7514509>`_, where the mathematics are translated into the bond-graph formalism. This describes the model in energetic terms and forces adherence to the laws of thermodynamics. 

.. figure::  exposure/FCU_AP.png
   :width: 100%
   :align: center
   :alt: FCU-AP

   Fig. 1. Bond-graph formulation of the cardiac action potential FCU, following the Kernik model. *Na channels* consist of fast-Na, Na\ :sub:`B`\; *K channels* consist of K1, K\ :sub:`r`\, K\ :sub:`s`\; *Ca channels* consist of Ca\ :sub:`B`\, TCC, pCa (not voltage dependent); *CaSR channels* consist of RyR, SERCA, Ca\ :sub:`leak`\; *Ca buffers* consist of CMDN, TRPN. Subscripts *i*, *o* and *SR* refer to the cytosolic, extracellular and sarcoplasmic compartments, respectively. Adapted from `Pan et. al. <http://dx.doi.org/10.1098/rspa.2018.0106>`_.


For an individual channel:

.. figure::  exposure/BG_individual_channel.png
   :width: 50%
   :align: center
   :alt: BG-individual-channel

   Fig. 2. Bond-graph formulation of an individual gated channel, for the transport of ion **q** from compartment *i* to *o*. **G** represents the gate, for gated channels. Membrane voltage **mem** modulates the GHK conductance of the channel and the gate. Adapted from Pan et. al.
   
For the above bond-graphs, a '0' node refers to a junction where all chemical potentials are the same. A '1' node refers to all fluxes being the same going in and out of the junction.

The following is a list of all available channels of this workspace. Unless otherwise stated, a source index of '1' indicates the model was based on `Kernik <https://models.physiomeproject.org/e/817>`_; a source of '2' on the `Luo-Rudy 1994 <https://models.physiomeproject.org/e/81>`_ model; '3' indicates that the bond-graph model was made by `Pan et. al. <http://dx.doi.org/10.1098/rspa.2018.0106>`_ (and mostly based on Luo Rudy). For cases '2' and '3', models have been reparameterised with Kernik cell dimensions. 

.. csv-table:: List of all available channels
   :header: "Abbreviation", "Name", "Kinetic source"
   :widths: 3, 10, 3
   
   "Ca\ :sub:`B`\", "Background calcium channel", "2"
   "Ca\ :sub:`leak`\", "Leakage calcium channel", "1"
   "CMDN", "Calmodulin buffer for calcium", "3"
   "fast-Na", "Fast sodium channel", "3"
   "K", "Time-dependent potassium channel", "3"
   "K1", "Inward rectifying potassium channel", "3"
   "K\ :sub:`ATP`\", "ATP-sensitive potassium channel", "`Clancy and Rudy <https://models.physiomeproject.org/e/514>`_" 
   "K\ :sub:`r`\", "Rapid delayed rectifying potassium channel", "1"
   "K\ :sub:`s`\", "Slow delayed rectifying potassium channel", "1"
   "LCC", "L-type calcium transporter", "3"
   "Na\ :sub:`B`\", "Background sodium channel", "2"
   "NaK", "Sodium-potassium pump", "3"
   "NCX", "Sodium-calcium exchanger", "3"
   "NS", "Non-specific calcium activated channel", "2"
   "pCa", "Sarcolemmal calcium pump", "1"
   "RyR", "Ryanodine receptor", "`Stern et. al. <https://models.physiomeproject.org/exposure/b060fdbcfae8c7d85e595c24d36ab11b>`_"
   "SERCA", "Sarcoplasmic-endoplasmic reticulum calcium pump", "3"
   "TCC", "T-type calcium transporter", "1"
   "TRPN", "Troponin buffer for calcium", "3"
   
Parameter finding
~~~~~~~~~~~~~~~~~
A description of the process to find bond-graph parameters is shown in the folder `parameter_finder <parameter_finder>`_, which relies on the:

1. stoichiometry of system
    
2. kinetic constants for forward/reverse reactions

  - If not already, all reactions are made reversible by assigning a small value to the reverse  direction.
  
3. `linear algebra script <https://models.physiomeproject.org/workspace/839/rawfile/174ee5da185a9a8254e6c9affcbeb67b30c81190/parameter_finder/find_BG_parameters_composite.py>`_, which outputs the bond-graph parameters K and :math:`{\kappa}`.

Here, this solve process is performed in Python. For a given composite model, the user must specify which components to include through the variable :code:`subsystem_names`. The stoichiometry for the whole system is automatically found.


Notes on creating a composite model
-------------------------------------
  - Step 3 (above) must be run each time a change is made to the composite model (e.g. addition or removal of a component).
  - Kinetic parameters fitted to a kinetic model must be found anew for a cell with different volumes and membrane capacitance values for each individual channel. This would apply to the GHK fits to a linear (ohmic) current, and to finding transition parameters for gating subcomponents of a gated channel. Following this, the bond-graph parameters are then recalculated.
  - New channels *must* have species names consistent with names already present in the model.

|

For further reading into the derivation of this process, please see `Pan et. al. <http://dx.doi.org/10.1098/rspa.2018.0106>`_.