Rendering of the source text

<?xml version="1.0" encoding="utf-8"?>
<!--
This CellML file was generated on 21/08/2009 at 9:19:39 at a.m. using:

COR (0.9.31.1309)
Copyright 2002-2009 Dr Alan Garny
http://cor.physiol.ox.ac.uk/ - cor@physiol.ox.ac.uk

CellML 1.0 was used to generate this model
http://www.cellml.org/
-->

<model xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" name="lemon_2003" cmeta:id="lemon_2003" xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cellml="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#">

<documentation xmlns="http://cellml.org/tmp-documentation">
    <article>
      <section id="sec_status">
        <title>Model Status</title>
        <para>
	This CellML model is the model which was used to produce the original results in the paper. The units have been checked and are consistent and this model is known to be valid CellML.
          </para>
      </section>
      <sect1 id="sec_structure">
        <title>Model Structure</title>
        
		<para>ABSTRACT:  A mathematical account is given of the processes governing the time courses of calcium ions (Ca2+), inositol 1,4,5-trisphosphate (IP(3)) and phosphatidylinositol 4,5-bisphosphate (PIP(2)) in single cells following the application of external agonist to metabotropic receptors. A model is constructed that incorporates the regulation of metabotropic receptor activity, the G-protein cascade and the Ca2+ dynamics in the cytosol. It is subsequently used to reproduce observations on the extent of desensitization and sequestration of the P(2)Y(2) receptor following its activation by uridine triphosphate (UTP). The theory predicts the dependence on agonist concentration of the change in the number of receptors in the membrane as well as the time course of disappearance of receptors from the plasmalemma, upon exposure to agonist. In addition, the extent of activation and desensitization of the receptor, using the calcium transients in cells initiated by exposure to agonist, is also predicted. Model predictions show the significance of membrane PIP(2) depletion and resupply on the time course of IP(3) and Ca2+ levels. Results of the modelling also reveal the importance of receptor recycling and PIP(2) resupply for maintaining Ca2+ and IP(3) levels during sustained application of agonist.
        </para>
       
	<para>
	The original paper is cited below:
	</para>
	<para>
	Metabotropic receptor activation, desensitization and sequestration-I: modelling calcium and inositol 1,4,5-trisphosphate dynamics following receptor activation. G.Lemon, W.G. Gibson and M.R. Bennett, 2003, Journal of Theoretical Biology, 223, 93-111. <ulink url="http://www.ncbi.nlm.nih.gov/pubmed/12782119">PubMed ID: 12782119</ulink>
	</para>

<informalfigure float="0" id="fig_reaction_diagram">
	<mediaobject>
		<imageobject>
			<objectinfo>
			<title>Diagram of model</title>
			</objectinfo>
			<imagedata fileref="lemon_2003.png"/>
		</imageobject>
	</mediaobject>
<caption>A schematic diagram of the processes described by the mathematical model.</caption>
</informalfigure>

      </sect1>
    </article>
  </documentation>

	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#"><rdf:Description rdf:nodeID="n1"><boundIntervals xmlns="http://www.cellml.org/metadata/simulation/1.0#"><rdf:Description rdf:nodeID="n2"/></boundIntervals></rdf:Description><rdf:Description rdf:about="#All_Model_3_BACK_ORIGINAL"><simulation xmlns="http://www.cellml.org/metadata/simulation/1.0#"><rdf:Description rdf:nodeID="n1"/></simulation></rdf:Description><rdf:Description rdf:nodeID="n2"><rdf:rest><rdf:Description rdf:about="http://www.w3.org/1999/02/22-rdf-syntax-ns#nil"/></rdf:rest><rdf:first><rdf:Description rdf:nodeID="n3"/></rdf:first></rdf:Description><rdf:Description rdf:nodeID="n3"><endingValue xmlns="http://www.cellml.org/metadata/simulation/1.0#" xml:lang="en">1000</endingValue></rdf:Description></rdf:RDF>   <units name="micromolar">
      <unit units="mole" prefix="micro"/>
      <unit units="litre" exponent="-1"/>
   </units>
   <units name="per_micromolar">
      <unit units="micromolar" exponent="-1"/>
   </units>
   <units name="per_micromole">
      <unit units="mole" prefix="micro" exponent="-1"/>
   </units>
   <units name="flux">
      <unit units="micromolar"/>
      <unit units="second" exponent="-1"/>
   </units>
   <units name="per_second">
      <unit units="second" exponent="-1"/>
   </units>
   <units name="per_micromolar_per_second">
      <unit units="micromolar" exponent="-1"/>
      <unit units="per_second"/>
   </units>
   <units name="micromolar_square_per_second">
      <unit units="micromolar" exponent="2"/>
      <unit units="per_second"/>
   </units>
   <component name="environment">
      <variable name="time" units="second" public_interface="out"/>
   </component>
   <component name="Parameters">
      <variable name="R_T" units="dimensionless" initial_value="2e4" public_interface="out"/>
      <variable name="K_1" units="micromolar" initial_value="5" public_interface="out"/>
      <variable name="K_2" units="micromolar" initial_value="100" public_interface="out"/>
      <variable name="k_r" units="per_second" initial_value="1.75e-4" public_interface="out"/>
      <variable name="k_p" units="per_second" initial_value="0.03" public_interface="out"/>
      <variable name="k_e" units="per_second" initial_value="6e-3" public_interface="out"/>
      <variable name="xi" units="dimensionless" initial_value="0.85" public_interface="out"/>
      <variable name="G_T" units="dimensionless" initial_value="1e5" public_interface="out"/>
      <variable name="k_deg" units="per_second" initial_value="1.25" public_interface="out"/>
      <variable name="k_a" units="per_second" initial_value="0.017" public_interface="out"/>
      <variable name="k_d" units="per_second" initial_value="0.15" public_interface="out"/>
      <variable name="PIP_2_T" units="dimensionless" initial_value="5e7" public_interface="out"/>
      <variable name="r_r" units="per_second" initial_value="10" public_interface="out"/>
      <variable name="delta" units="dimensionless" initial_value="1.238e-3" public_interface="out"/>
      <variable name="K_c" units="micromolar" initial_value="0.4" public_interface="out"/>
      <variable name="alpha" units="per_second" initial_value="2.781e-5" public_interface="out"/>
      <variable name="N_a" units="per_micromole" initial_value="6.02252e17" public_interface="out"/>
      <variable name="v" units="litre" initial_value="5e-13" public_interface="out"/>
      <variable name="epsilon_r" units="dimensionless" initial_value="0.185" public_interface="out"/>
      <variable name="d_1" units="micromolar" initial_value="0.13" public_interface="out"/>
      <variable name="d_2" units="micromolar" initial_value="1.05" public_interface="out"/>
      <variable name="d_3" units="micromolar" initial_value="0.943" public_interface="out"/>
      <variable name="d_5" units="micromolar" initial_value="0.0823" public_interface="out"/>
      <variable name="a_2" units="per_micromolar_per_second" initial_value="0.2" public_interface="out"/>
      <variable name="B_e" units="micromolar" initial_value="150" public_interface="out"/>
      <variable name="K_e" units="micromolar" initial_value="10" public_interface="out"/>
      <variable name="B_ER" units="micromolar" initial_value="120000" public_interface="out"/>
      <variable name="K_ER" units="micromolar" initial_value="1200" public_interface="out"/>
      <variable name="B_x" units="micromolar" initial_value="50" public_interface="out"/>
      <variable name="K_x" units="micromolar" initial_value="0.2" public_interface="out"/>
      <variable name="k_3" units="micromolar" initial_value="0.4" public_interface="out"/>
      <variable name="eta_1" units="per_second" initial_value="575" public_interface="out"/>
      <variable name="eta_2" units="per_second" initial_value="5.2" public_interface="out"/>
      <variable name="eta_3" units="flux" initial_value="45" public_interface="out"/>
      <variable name="C_T" units="micromolar" initial_value="67" public_interface="out"/>
   </component>
   <component name="ligand">
      <variable name="L" units="micromolar" initial_value="1000" public_interface="out"/>
   </component>
   <component name="RS">
      <variable name="RS" units="dimensionless" initial_value="1.7e4" public_interface="out"/>
      <variable name="R_T" units="dimensionless" public_interface="in"/>
      <variable name="RS_p" units="dimensionless" public_interface="in"/>
      <variable name="L" units="micromolar" public_interface="in"/>
      <variable name="K_1" units="micromolar" public_interface="in"/>
      <variable name="k_r" units="per_second" public_interface="in"/>
      <variable name="k_p" units="per_second" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>RS</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>k_r</ci>
                     <ci>R_T</ci>
                  </apply>
                  <apply>
                     <times/>
                     <apply>
                        <plus/>
                        <ci>k_r</ci>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <ci>k_p</ci>
                              <ci>L</ci>
                           </apply>
                           <apply>
                              <plus/>
                              <ci>K_1</ci>
                              <ci>L</ci>
                           </apply>
                        </apply>
                     </apply>
                     <ci>RS</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>k_r</ci>
                  <ci>RS_p</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="RS_p">
      <variable name="RS_p" units="dimensionless" initial_value="0" public_interface="out"/>
      <variable name="R_T" units="dimensionless" public_interface="in"/>
      <variable name="RS" units="dimensionless" public_interface="in"/>
      <variable name="L" units="micromolar" public_interface="in"/>
      <variable name="K_1" units="micromolar" public_interface="in"/>
      <variable name="K_2" units="micromolar" public_interface="in"/>
      <variable name="k_r" units="per_second" public_interface="in"/>
      <variable name="k_p" units="per_second" public_interface="in"/>
      <variable name="k_e" units="per_second" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>RS_p</ci>
            </apply>
            <apply>
               <times/>
               <ci>L</ci>
               <apply>
                  <minus/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>k_p</ci>
                        <ci>RS</ci>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>K_1</ci>
                        <ci>L</ci>
                     </apply>
                  </apply>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>k_e</ci>
                        <ci>RS_p</ci>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>K_2</ci>
                        <ci>L</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="G_GTP">
      <variable name="G" units="dimensionless" initial_value="0" public_interface="out"/>
      <variable name="G_T" units="dimensionless" public_interface="in"/>
      <variable name="delta" units="dimensionless" public_interface="in"/>
      <variable name="k_a" units="per_second" public_interface="in"/>
      <variable name="k_d" units="per_second" public_interface="in"/>
      <variable name="rho_r" units="dimensionless" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>G</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>k_a</ci>
                  <apply>
                     <plus/>
                     <ci>delta</ci>
                     <ci>rho_r</ci>
                  </apply>
                  <apply>
                     <minus/>
                     <ci>G_T</ci>
                     <ci>G</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>k_d</ci>
                  <ci>G</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="IP_3">
      <variable name="IP_3" units="micromolar" initial_value="0.01" public_interface="out"/>
      <variable name="PIP_2" units="dimensionless" public_interface="in"/>
      <variable name="r_h" units="per_second" public_interface="in"/>
      <variable name="k_deg" units="per_second" public_interface="in"/>
      <variable name="N_a" units="per_micromole" public_interface="in"/>
      <variable name="v" units="litre" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>IP_3</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>r_h</ci>
                     <ci>PIP_2</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>N_a</ci>
                     <ci>v</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>k_deg</ci>
                  <ci>IP_3</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="PIP_2">
      <variable name="PIP_2" units="dimensionless" initial_value="49997000" public_interface="out"/>
      <variable name="IP_3" units="micromolar" public_interface="in"/>
      <variable name="r_h" units="per_second" public_interface="in"/>
      <variable name="r_r" units="per_second" public_interface="in"/>
      <variable name="PIP_2_T" units="dimensionless" public_interface="in"/>
      <variable name="k_deg" units="per_second" public_interface="in"/>
      <variable name="N_a" units="per_micromole" public_interface="in"/>
      <variable name="v" units="litre" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>PIP_2</ci>
            </apply>
            <apply>
               <plus/>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
                        <apply>
                           <plus/>
                           <ci>r_h</ci>
                           <ci>r_r</ci>
                        </apply>
                     </apply>
                     <ci>PIP_2</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>r_r</ci>
                     <ci>N_a</ci>
                     <ci>v</ci>
                     <ci>IP_3</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>r_r</ci>
                  <ci>PIP_2_T</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="r_h">
      <variable name="r_h" units="per_second" public_interface="out"/>
      <variable name="alpha" units="per_second" public_interface="in"/>
      <variable name="C" units="micromolar" public_interface="in"/>
      <variable name="K_c" units="micromolar" public_interface="in"/>
      <variable name="G" units="dimensionless" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>r_h</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>alpha</ci>
                     <ci>C</ci>
                  </apply>
                  <apply>
                     <plus/>
                     <ci>K_c</ci>
                     <ci>C</ci>
                  </apply>
               </apply>
               <ci>G</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="rho_r">
      <variable name="rho_r" units="dimensionless" public_interface="out"/>
      <variable name="L" units="micromolar" public_interface="in"/>
      <variable name="RS" units="dimensionless" public_interface="in"/>
      <variable name="R_T" units="dimensionless" public_interface="in"/>
      <variable name="K_1" units="micromolar" public_interface="in"/>
      <variable name="xi" units="dimensionless" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>rho_r</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>L</ci>
                  <ci>RS</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>xi</ci>
                  <ci>R_T</ci>
                  <apply>
                     <plus/>
                     <ci>K_1</ci>
                     <ci>L</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="C">
      <variable name="C" units="micromolar" initial_value="0.0961" public_interface="out"/>
      <variable name="C_ER" units="micromolar" public_interface="in"/>
      <variable name="beta" units="dimensionless" public_interface="in"/>
      <variable name="epsilon_r" units="dimensionless" public_interface="in"/>
      <variable name="eta_1" units="per_second" public_interface="in"/>
      <variable name="eta_2" units="per_second" public_interface="in"/>
      <variable name="eta_3" units="flux" public_interface="in"/>
      <variable name="k_3" units="micromolar" public_interface="in"/>
      <variable name="h" units="dimensionless" public_interface="in"/>
      <variable name="m_infinit" units="dimensionless" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>C</ci>
            </apply>
            <apply>
               <times/>
               <ci>beta</ci>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>epsilon_r</ci>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <ci>eta_1</ci>
                           <apply>
                              <power/>
                              <ci>m_infinit</ci>
                              <cn cellml:units="dimensionless">3</cn>
                           </apply>
                           <apply>
                              <power/>
                              <ci>h</ci>
                              <cn cellml:units="dimensionless">3</cn>
                           </apply>
                        </apply>
                        <ci>eta_2</ci>
                     </apply>
                     <apply>
                        <minus/>
                        <ci>C_ER</ci>
                        <ci>C</ci>
                     </apply>
                  </apply>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>eta_3</ci>
                        <apply>
                           <power/>
                           <ci>C</ci>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                     <apply>
                        <plus/>
                        <apply>
                           <power/>
                           <ci>k_3</ci>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                        <apply>
                           <power/>
                           <ci>C</ci>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="h">
      <variable name="h" units="dimensionless" initial_value="0.6155" public_interface="out"/>
      <variable name="h_infinit" units="dimensionless" public_interface="in"/>
      <variable name="tau_h" units="second" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>h</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>h_infinit</ci>
                  <ci>h</ci>
               </apply>
               <ci>tau_h</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="tau_h">
      <variable name="tau_h" units="second" public_interface="out"/>
      <variable name="a_2" units="per_micromolar_per_second" public_interface="in"/>
      <variable name="zeta" units="micromolar" public_interface="in"/>
      <variable name="C" units="micromolar" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>tau_h</ci>
            <apply>
               <power/>
               <apply>
                  <times/>
                  <ci>a_2</ci>
                  <apply>
                     <plus/>
                     <ci>zeta</ci>
                     <ci>C</ci>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <cn cellml:units="dimensionless">1</cn>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="h_infinit">
      <variable name="h_infinit" units="dimensionless" public_interface="out"/>
      <variable name="zeta" units="micromolar" public_interface="in"/>
      <variable name="C" units="micromolar" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>h_infinit</ci>
            <apply>
               <divide/>
               <ci>zeta</ci>
               <apply>
                  <plus/>
                  <ci>zeta</ci>
                  <ci>C</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="zeta">
      <variable name="zeta" units="micromolar" public_interface="out"/>
      <variable name="d_1" units="micromolar" public_interface="in"/>
      <variable name="d_2" units="micromolar" public_interface="in"/>
      <variable name="d_3" units="micromolar" public_interface="in"/>
      <variable name="IP_3" units="micromolar" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>zeta</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>d_2</ci>
                  <apply>
                     <plus/>
                     <ci>IP_3</ci>
                     <ci>d_1</ci>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <ci>IP_3</ci>
                  <ci>d_3</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="m_infinit">
      <variable name="m_infinit" units="dimensionless" public_interface="out"/>
      <variable name="d_1" units="micromolar" public_interface="in"/>
      <variable name="d_5" units="micromolar" public_interface="in"/>
      <variable name="IP_3" units="micromolar" public_interface="in"/>
      <variable name="C" units="micromolar" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>m_infinit</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <ci>IP_3</ci>
                     <apply>
                        <plus/>
                        <ci>d_1</ci>
                        <ci>IP_3</ci>
                     </apply>
                  </apply>
                  <ci>C</ci>
               </apply>
               <apply>
                  <plus/>
                  <ci>d_5</ci>
                  <ci>C</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="beta">
      <variable name="beta" units="dimensionless" public_interface="out"/>
      <variable name="K_e" units="micromolar" public_interface="in"/>
      <variable name="B_e" units="micromolar" public_interface="in"/>
      <variable name="K_x" units="micromolar" public_interface="in"/>
      <variable name="B_x" units="micromolar" public_interface="in"/>
      <variable name="C" units="micromolar" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>beta</ci>
            <apply>
               <power/>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>K_e</ci>
                        <ci>B_e</ci>
                     </apply>
                     <apply>
                        <power/>
                        <apply>
                           <plus/>
                           <ci>K_e</ci>
                           <ci>C</ci>
                        </apply>
                        <cn cellml:units="dimensionless">2</cn>
                     </apply>
                  </apply>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>K_x</ci>
                        <ci>B_x</ci>
                     </apply>
                     <apply>
                        <power/>
                        <apply>
                           <plus/>
                           <ci>K_x</ci>
                           <ci>C</ci>
                        </apply>
                        <cn cellml:units="dimensionless">2</cn>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <cn cellml:units="dimensionless">1</cn>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="gamma">
      <variable name="gamma" units="dimensionless" public_interface="out"/>
      <variable name="K_e" units="micromolar" public_interface="in"/>
      <variable name="B_e" units="micromolar" public_interface="in"/>
      <variable name="K_x" units="micromolar" public_interface="in"/>
      <variable name="B_x" units="micromolar" public_interface="in"/>
      <variable name="C" units="micromolar" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>gamma</ci>
            <apply>
               <power/>
               <apply>
                  <plus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <divide/>
                     <ci>B_e</ci>
                     <apply>
                        <plus/>
                        <ci>K_e</ci>
                        <ci>C</ci>
                     </apply>
                  </apply>
                  <apply>
                     <divide/>
                     <ci>B_x</ci>
                     <apply>
                        <plus/>
                        <ci>K_x</ci>
                        <ci>C</ci>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <cn cellml:units="dimensionless">1</cn>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="C_ER">
      <variable name="C_ER" units="micromolar" public_interface="out"/>
      <variable name="K_ER" units="micromolar" public_interface="in"/>
      <variable name="B_ER" units="micromolar" public_interface="in"/>
      <variable name="epsilon_r" units="dimensionless" public_interface="in"/>
      <variable name="C_T" units="micromolar" public_interface="in"/>
      <variable name="C" units="micromolar" public_interface="in"/>
      <variable name="gamma" units="dimensionless" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>C_ER</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <ci>K_ER</ci>
                  <apply>
                     <times/>
                     <ci>B_ER</ci>
                     <ci>epsilon_r</ci>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <ci>C_T</ci>
                  <apply>
                     <divide/>
                     <ci>C</ci>
                     <ci>gamma</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="RS_E">
      <variable name="RS_E" units="dimensionless" public_interface="out"/>
      <variable name="R_T" units="dimensionless" public_interface="in"/>
      <variable name="L" units="micromolar" public_interface="in"/>
      <variable name="K_1" units="micromolar" public_interface="in"/>
      <variable name="K_2" units="micromolar" public_interface="in"/>
      <variable name="k_r" units="per_second" public_interface="in"/>
      <variable name="k_p" units="per_second" public_interface="in"/>
      <variable name="k_e" units="per_second" public_interface="in"/>
      <variable name="xi" units="dimensionless" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>RS_E</ci>
            <apply>
               <plus/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>k_r</ci>
                        <apply>
                           <plus/>
                           <cn cellml:units="dimensionless">1</cn>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <apply>
                                    <divide/>
                                    <ci>k_p</ci>
                                    <ci>k_e</ci>
                                 </apply>
                                 <apply>
                                    <plus/>
                                    <ci>K_2</ci>
                                    <ci>L</ci>
                                 </apply>
                              </apply>
                              <apply>
                                 <plus/>
                                 <ci>K_1</ci>
                                 <ci>L</ci>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                     <apply>
                        <plus/>
                        <ci>k_r</ci>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <ci>k_p</ci>
                              <ci>L</ci>
                           </apply>
                           <apply>
                              <plus/>
                              <ci>K_1</ci>
                              <ci>L</ci>
                           </apply>
                        </apply>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <apply>
                                 <divide/>
                                 <apply>
                                    <times/>
                                    <ci>k_r</ci>
                                    <ci>k_p</ci>
                                 </apply>
                                 <ci>k_e</ci>
                              </apply>
                              <apply>
                                 <plus/>
                                 <ci>K_2</ci>
                                 <ci>L</ci>
                              </apply>
                           </apply>
                           <apply>
                              <plus/>
                              <ci>K_1</ci>
                              <ci>L</ci>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
                  <ci>xi</ci>
                  <ci>R_T</ci>
               </apply>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>xi</ci>
                  </apply>
                  <ci>R_T</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <connection>
      <map_components component_1="RS" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="RS_p" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="environment" component_2="G_GTP"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="environment" component_2="IP_3"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="environment" component_2="PIP_2"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="environment" component_2="C"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="h" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="RS" component_2="ligand"/>
      <map_variables variable_1="L" variable_2="L"/>
   </connection>
   <connection>
      <map_components component_1="RS_p" component_2="ligand"/>
      <map_variables variable_1="L" variable_2="L"/>
   </connection>
   <connection>
      <map_components component_1="rho_r" component_2="ligand"/>
      <map_variables variable_1="L" variable_2="L"/>
   </connection>
   <connection>
      <map_components component_1="RS" component_2="Parameters"/>
      <map_variables variable_1="R_T" variable_2="R_T"/>
      <map_variables variable_1="K_1" variable_2="K_1"/>
      <map_variables variable_1="k_r" variable_2="k_r"/>
      <map_variables variable_1="k_p" variable_2="k_p"/>
   </connection>
   <connection>
      <map_components component_1="RS" component_2="rho_r"/>
      <map_variables variable_1="RS" variable_2="RS"/>
   </connection>
   <connection>
      <map_components component_1="RS_p" component_2="Parameters"/>
      <map_variables variable_1="R_T" variable_2="R_T"/>
      <map_variables variable_1="K_1" variable_2="K_1"/>
      <map_variables variable_1="K_2" variable_2="K_2"/>
      <map_variables variable_1="k_r" variable_2="k_r"/>
      <map_variables variable_1="k_p" variable_2="k_p"/>
      <map_variables variable_1="k_e" variable_2="k_e"/>
   </connection>
   <connection>
      <map_components component_1="RS" component_2="RS_p"/>
      <map_variables variable_1="RS" variable_2="RS"/>
      <map_variables variable_1="RS_p" variable_2="RS_p"/>
   </connection>
   <connection>
      <map_components component_1="G_GTP" component_2="rho_r"/>
      <map_variables variable_1="rho_r" variable_2="rho_r"/>
   </connection>
   <connection>
      <map_components component_1="G_GTP" component_2="Parameters"/>
      <map_variables variable_1="G_T" variable_2="G_T"/>
      <map_variables variable_1="delta" variable_2="delta"/>
      <map_variables variable_1="k_a" variable_2="k_a"/>
      <map_variables variable_1="k_d" variable_2="k_d"/>
   </connection>
   <connection>
      <map_components component_1="IP_3" component_2="PIP_2"/>
      <map_variables variable_1="IP_3" variable_2="IP_3"/>
      <map_variables variable_1="PIP_2" variable_2="PIP_2"/>
   </connection>
   <connection>
      <map_components component_1="IP_3" component_2="r_h"/>
      <map_variables variable_1="r_h" variable_2="r_h"/>
   </connection>
   <connection>
      <map_components component_1="IP_3" component_2="Parameters"/>
      <map_variables variable_1="k_deg" variable_2="k_deg"/>
      <map_variables variable_1="N_a" variable_2="N_a"/>
      <map_variables variable_1="v" variable_2="v"/>
   </connection>
   <connection>
      <map_components component_1="PIP_2" component_2="r_h"/>
      <map_variables variable_1="r_h" variable_2="r_h"/>
   </connection>
   <connection>
      <map_components component_1="PIP_2" component_2="Parameters"/>
      <map_variables variable_1="r_r" variable_2="r_r"/>
      <map_variables variable_1="PIP_2_T" variable_2="PIP_2_T"/>
      <map_variables variable_1="k_deg" variable_2="k_deg"/>
      <map_variables variable_1="N_a" variable_2="N_a"/>
      <map_variables variable_1="v" variable_2="v"/>
   </connection>
   <connection>
      <map_components component_1="G_GTP" component_2="r_h"/>
      <map_variables variable_1="G" variable_2="G"/>
   </connection>
   <connection>
      <map_components component_1="r_h" component_2="Parameters"/>
      <map_variables variable_1="alpha" variable_2="alpha"/>
      <map_variables variable_1="K_c" variable_2="K_c"/>
   </connection>
   <connection>
      <map_components component_1="rho_r" component_2="Parameters"/>
      <map_variables variable_1="R_T" variable_2="R_T"/>
      <map_variables variable_1="K_1" variable_2="K_1"/>
      <map_variables variable_1="xi" variable_2="xi"/>
   </connection>
   <connection>
      <map_components component_1="C" component_2="r_h"/>
      <map_variables variable_1="C" variable_2="C"/>
   </connection>
   <connection>
      <map_components component_1="C" component_2="Parameters"/>
      <map_variables variable_1="epsilon_r" variable_2="epsilon_r"/>
      <map_variables variable_1="eta_1" variable_2="eta_1"/>
      <map_variables variable_1="eta_2" variable_2="eta_2"/>
      <map_variables variable_1="eta_3" variable_2="eta_3"/>
      <map_variables variable_1="k_3" variable_2="k_3"/>
   </connection>
   <connection>
      <map_components component_1="C" component_2="h"/>
      <map_variables variable_1="h" variable_2="h"/>
   </connection>
   <connection>
      <map_components component_1="C" component_2="tau_h"/>
      <map_variables variable_1="C" variable_2="C"/>
   </connection>
   <connection>
      <map_components component_1="C" component_2="h_infinit"/>
      <map_variables variable_1="C" variable_2="C"/>
   </connection>
   <connection>
      <map_components component_1="C" component_2="m_infinit"/>
      <map_variables variable_1="m_infinit" variable_2="m_infinit"/>
      <map_variables variable_1="C" variable_2="C"/>
   </connection>
   <connection>
      <map_components component_1="C" component_2="beta"/>
      <map_variables variable_1="beta" variable_2="beta"/>
      <map_variables variable_1="C" variable_2="C"/>
   </connection>
   <connection>
      <map_components component_1="C" component_2="C_ER"/>
      <map_variables variable_1="C_ER" variable_2="C_ER"/>
      <map_variables variable_1="C" variable_2="C"/>
   </connection>
   <connection>
      <map_components component_1="h" component_2="h_infinit"/>
      <map_variables variable_1="h_infinit" variable_2="h_infinit"/>
   </connection>
   <connection>
      <map_components component_1="h" component_2="tau_h"/>
      <map_variables variable_1="tau_h" variable_2="tau_h"/>
   </connection>
   <connection>
      <map_components component_1="tau_h" component_2="Parameters"/>
      <map_variables variable_1="a_2" variable_2="a_2"/>
   </connection>
   <connection>
      <map_components component_1="tau_h" component_2="zeta"/>
      <map_variables variable_1="zeta" variable_2="zeta"/>
   </connection>
   <connection>
      <map_components component_1="h_infinit" component_2="zeta"/>
      <map_variables variable_1="zeta" variable_2="zeta"/>
   </connection>
   <connection>
      <map_components component_1="zeta" component_2="Parameters"/>
      <map_variables variable_1="d_1" variable_2="d_1"/>
      <map_variables variable_1="d_2" variable_2="d_2"/>
      <map_variables variable_1="d_3" variable_2="d_3"/>
   </connection>
   <connection>
      <map_components component_1="zeta" component_2="IP_3"/>
      <map_variables variable_1="IP_3" variable_2="IP_3"/>
   </connection>
   <connection>
      <map_components component_1="m_infinit" component_2="Parameters"/>
      <map_variables variable_1="d_1" variable_2="d_1"/>
      <map_variables variable_1="d_5" variable_2="d_5"/>
   </connection>
   <connection>
      <map_components component_1="m_infinit" component_2="IP_3"/>
      <map_variables variable_1="IP_3" variable_2="IP_3"/>
   </connection>
   <connection>
      <map_components component_1="beta" component_2="Parameters"/>
      <map_variables variable_1="K_e" variable_2="K_e"/>
      <map_variables variable_1="K_x" variable_2="K_x"/>
      <map_variables variable_1="B_e" variable_2="B_e"/>
      <map_variables variable_1="B_x" variable_2="B_x"/>
   </connection>
   <connection>
      <map_components component_1="gamma" component_2="Parameters"/>
      <map_variables variable_1="K_e" variable_2="K_e"/>
      <map_variables variable_1="K_x" variable_2="K_x"/>
      <map_variables variable_1="B_e" variable_2="B_e"/>
      <map_variables variable_1="B_x" variable_2="B_x"/>
   </connection>
   <connection>
      <map_components component_1="gamma" component_2="C"/>
      <map_variables variable_1="C" variable_2="C"/>
   </connection>
   <connection>
      <map_components component_1="C_ER" component_2="Parameters"/>
      <map_variables variable_1="K_ER" variable_2="K_ER"/>
      <map_variables variable_1="B_ER" variable_2="B_ER"/>
      <map_variables variable_1="epsilon_r" variable_2="epsilon_r"/>
      <map_variables variable_1="C_T" variable_2="C_T"/>
   </connection>
   <connection>
      <map_components component_1="C_ER" component_2="gamma"/>
      <map_variables variable_1="gamma" variable_2="gamma"/>
   </connection>
   <connection>
      <map_components component_1="RS_E" component_2="Parameters"/>
      <map_variables variable_1="K_1" variable_2="K_1"/>
      <map_variables variable_1="K_2" variable_2="K_2"/>
      <map_variables variable_1="k_p" variable_2="k_p"/>
      <map_variables variable_1="k_e" variable_2="k_e"/>
      <map_variables variable_1="k_r" variable_2="k_r"/>
      <map_variables variable_1="xi" variable_2="xi"/>
      <map_variables variable_1="R_T" variable_2="R_T"/>
   </connection>
   <connection>
      <map_components component_1="RS_E" component_2="ligand"/>
      <map_variables variable_1="L" variable_2="L"/>
   </connection>
   
   <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">
    <foo xmlns="http://www.example.org/foo/bar/baz" />
    <!--
      The following RDF block contains metadata that applies to this document
      as a whole, as indicated by the empty about attribute on the 
      <rdf:Description> element.
    -->
    <rdf:Description rdf:about="">
      <!--
        The Model Builder Metadata.  The Dublin Core "creator" element is used  
        to indicate the person who translated the model into CellML.
      -->
      <dc:creator rdf:parseType="Resource">
        <vCard:N rdf:parseType="Resource">
          <vCard:Family>Liu</vCard:Family>
          <vCard:Given>Wei</vCard:Given>
		</vCard:N>
        <vCard:EMAIL rdf:parseType="Resource">
          <rdf:value>wliu052@aucklanduni.ac.nz</rdf:value>
          <rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
        </vCard:EMAIL>
        <vCard:ORG rdf:parseType="Resource">
          <vCard:Orgname>The University of Auckland</vCard:Orgname>
        </vCard:ORG>
      </dc:creator>
      
      <!--
        The Creation Date metadata. This is the date on which the model
        was translated into CellML.
      -->
      <dcterms:created rdf:parseType="Resource">
        <dcterms:W3CDTF>2009-09-29</dcterms:W3CDTF>
      </dcterms:created>
      
      <!--
        The Modification History metadata. This lists the changes that have been
        made to the document, who made the changes, and when they were made.
      -->
     
    </rdf:Description>
    
    <!--
      The following metadata refers to the model itself, as indicated by the
      reference to the ID "lemon_2003", which is    
      declared on the <model> element.
    -->
    <rdf:Description rdf:about="#lemon_2003">
      <!-- A human readable name for the model. -->
      <dc:title>
        The Lemon et al. 2003 model of calcium and inositol 1,4,5-trisphosphate dynamics following receptor activation
      </dc:title>
      
      <!-- A comment regarding the model. -->
      <cmeta:comment rdf:parseType="Resource">
        <rdf:value>
          This is the CellML description of Lemon et al.'s mathematical model of calcium and inositol 1,4,5-trisphosphate dynamics following receptor activation
        </rdf:value>
        <!-- The creator of the comment. -->
        <dc:creator rdf:parseType="Resource">
          <vCard:FN>Catherine Lloyd</vCard:FN>
        </dc:creator>
      </cmeta:comment>
      
      <!--  Keyword(s) -->
      <bqs:reference rdf:parseType="Resource">
        <dc:subject rdf:parseType="Resource">
          <bqs:subject_type>keyword</bqs:subject_type>
          <rdf:value>
            <rdf:Bag>
              <rdf:li>calcium dynamics</rdf:li>
            </rdf:Bag>
          </rdf:value>
        </dc:subject>
      </bqs:reference>

      <bqs:reference rdf:parseType="Resource">
        <bqs:Pubmed_id>12782119</bqs:Pubmed_id>
        <bqs:JournalArticle rdf:parseType="Resource">
          <dc:creator>
            <rdf:Seq>
              <rdf:li rdf:parseType="Resource">
                <bqs:Person rdf:parseType="Resource">
                  <vCard:N rdf:parseType="Resource">
                    <vCard:Family>Lemon</vCard:Family>
                    <vCard:Given>G</vCard:Given>
                  </vCard:N>
                </bqs:Person>
              </rdf:li>
              <rdf:li rdf:parseType="Resource">
                <bqs:Person rdf:parseType="Resource">
                  <vCard:N rdf:parseType="Resource">
                    <vCard:Family>Gibson</vCard:Family>
                    <vCard:Given>W</vCard:Given>
					<vCard:Other>G</vCard:Other>
				  </vCard:N>
                </bqs:Person>
              </rdf:li>
			  <rdf:li rdf:parseType="Resource">
                <bqs:Person rdf:parseType="Resource">
                  <vCard:N rdf:parseType="Resource">
                    <vCard:Family>Bennett</vCard:Family>
                    <vCard:Given>M</vCard:Given>
					<vCard:Other>R</vCard:Other>
                  </vCard:N>
                </bqs:Person>
              </rdf:li>
            </rdf:Seq>
          </dc:creator>
          <dc:title>
           Metabotropic receptor activation, desensitization and sequestration-I: modelling calcium and inositol 1,4,5-trisphosphate dynamics following receptor activation
          </dc:title>
          <dcterms:issued rdf:parseType="Resource">
            <dcterms:W3CDTF>2003-07-07</dcterms:W3CDTF>
          </dcterms:issued>
          <bqs:Journal rdf:parseType="Resource">
            <dc:title>Journal of Theoretical Biology</dc:title>
          </bqs:Journal>
          <bqs:volume>223</bqs:volume>
          <bqs:first_page>93</bqs:first_page>
          <bqs:last_page>111</bqs:last_page>
        </bqs:JournalArticle>
      </bqs:reference>
    </rdf:Description>
  </rdf:RDF> 
   
</model>