Location: BG_PLC @ 98c2d94ee9d4 / exposure / exposure_frontpage.rst

Author:
Shelley Fong <s.fong@auckland.ac.nz>
Date:
2022-05-02 13:25:40+12:00
Desc:
delete sentence
Permanent Source URI:
https://models.cellml.org/workspace/848/rawfile/98c2d94ee9d4dbdfcd74f4bc2dde0c8eeb9e61a5/exposure/exposure_frontpage.rst

About this model
====================

This is a bond-graph model of the metabolism of the phospholipase-C beta variant (PLC) and the associated Gq protein in the mammalian cell.

    **INPUTS:** 
        - Gq protein (Gq) 
        
    **OUTPUTS:** 
        - Change in molar amount of diacylglycerol (DAG), inositol trisphosphate (IP3), calcium (Ca)
        
    **REACTIONS:** 
        - R0: Conversion of DAG into phosphatidylcholine (PC)
        - R1: Conversion of IP3 into inositol (Ino)
        - R2: PLC complex CaPLC catalysing the splitting of phosphatidylinositol biphosphate (PIP2) into DAG and IP3
        - R3: Similar to R2, but with enzyme complex CaGqPLC
        - R4: Binding of Ca and PLC
        - R5: Binding of Gq and PLC
        - R6: Binding of Ca and GqPLC
        - R7: Unbinding of Gq from CaGqPLC
        - R8: Unbinding of Gq (in G*GDP form) from CaGqPLC

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

The current CellML implementation runs in OpenCOR.


Model overview
===================
This model is based on existing kinetic model, 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/BG_PLC.png
   :width: 100%
   :align: center
   :alt: BG PLC reaction
   
   Fig. 1. Bond-graph formulation of the PLC network


|

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.

.. csv-table:: List of chemical species
   :header: "Abbreviation", "Name"
   :widths: 5, 15 
   
   "PLC", "Phospholipase-C (beta variant)"
   "Ca", "Calcium ion"
   "Gq", "Gq protein"
   "DAG", "Diacylglycerol"
   "IP3", "Inositol trisphosphate"
   "Ino", "Inositol"
   "PIP2", "Phosphatidylinositol biphosphate"
   "PC", "Phosphatidylcholine"
   
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/848/rawfile/6239317b4238524d2bf82a54c50849558a4d67d7/parameter_finder/kinetic_parameters_PLC.py>`_. 

Here, this solve process is performed in Python.


Original kinetic model
======================
Bhalla and Iyengar: `Emergent properties of networks of biological signaling pathways. <https://models.physiomeproject.org/e/7e3>`_