Location: Thomsen, Neubig, 1989 @ 851326f149c4 / thomsen_neubig_1989.cellml

Author:
pmr2.import <nobody@models.cellml.org>
Date:
2006-07-09 07:36:31+12:00
Desc:
committing version01 of thomsen_neubig_1989
Permanent Source URI:
https://models.physiomeproject.org/workspace/thomsen_neubig_1989/rawfile/851326f149c4230e68ad4739ddc04c52cb8e409d/thomsen_neubig_1989.cellml

<?xml version='1.0' encoding='utf-8'?>
<!--  FILE :  G_protein_activation_1989.xml

CREATED :  14th November 2001

LAST MODIFIED : 9th April 2003

AUTHOR :  Catherine Lloyd
          Department of Engineering Science
          The University of Auckland
          
MODEL STATUS :  This model conforms to the CellML 1.0 Specification released on
10th August 2001, and the 16/01/2002 CellML Metadata 1.0 Specification.

DESCRIPTION :  This file contains a CellML description of a G-protein activation
pathway, based on the mathematical model of Thomsen and Neubig (1989).  

CHANGES :
  21/01/2002 - AAC - Updated metadata to conform to the 16/1/02 CellML Metadata
                     1.0 Specification.   
  22/07/2002 - CML - Added more metadata. 
  09/04/2003 - AAC - Added publication date information.  
--><model xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#" 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:cellml="http://www.cellml.org/cellml/1.0#" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:pathway_editor="http://www.physiome.com/pathway_editor/1.0#" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" pathway_editor:rendering_config_file="G protein activation pathway      1989render.xml" cmeta:id="thomsen_neubig_1989_version01" name="thomsen_neubig_1989_version01">
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
  <articleinfo>
  <title>Model Of G-Protein Activation</title>
  <author>
    <firstname>Catherine</firstname>
          <surname>Lloyd</surname>
    <affiliation>
      <shortaffil>Bioengineering Institute, University of Auckland</shortaffil>
    </affiliation>
  </author>
</articleinfo>
  <section id="sec_status">
    <title>Model Status</title>
    <para>
            This is the original unchecked version of the model imported from the previous
            CellML model repository, 24-Jan-2006.
          </para>
  </section>
  <sect1 id="sec_structure">
<title>Model Structure</title>

<para>
G-proteins, or guanine nucleotide-binding proteins, serve to couple receptors to membrane-bound enzymes such as adenylate cyclase. In turn this induces the generation of secondary messengers such as cAMP for purposes of signal transduction. The inactive form of the G-protein consists of alpha, beta, and gamma subunits with a molecule of GDP bound to the alpha subunit. When a ligand-bound receptor interacts with the G-protein, it catalyses the exchange of GDP for GTP and activates the G-protein. The G-protein is then released from the receptor and it dissociates into separate beta-gamma and alpha-GTP (active) subunits. Active G-proteins are returned to their inactive state upon the hydrolysis of GTP to GDP by the intrinsic GTPase activity of the alpha-subunit and the alpha-GDP and beta-gamma subunits can recombine. (See <xref linkend="fig_pathway_diagram"/> below).
</para>

<para>
This description of a general G-protein activation pathway is based on the mathematical model described by W.J Thomsen and R.R. Neubig in their 1989 paper. The model does not include the steps of dissociation and recombination of the alpha and beta-gamma subunits of the G-protein. In addition, it is assumed that the number of membrane bound receptors (Rt) is constant. 
</para>

<para>
The complete original paper reference is cited below:
</para>

<para>
Rapid kinetics of alpha 2-adrenergic inhibition of adenylate cyclase. Evidence for a distal rate-limiting step, W.J. Thomsen and R.R. Neubig, 1989, <ulink url="http://pubs.acs.org/journals/bichaw/">
            <emphasis>Biochemistry</emphasis>
          </ulink>, 28, 8778-8786. <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&amp;db=PubMed&amp;list_uids=2574993&amp;dopt=Abstract">PubMed ID: 2574993</ulink>
</para>

<para>
The raw CellML description of the G-protein activation pathway model can be downloaded in various formats as described in <xref linkend="sec_download_this_model"/>. For an example of a more complete documentation of another real reaction pathway, see <ulink url="${HTML_EXMPL_BI_EGF_INTRO}">The Bhalla Iyengar EGF Pathway Model, 1999</ulink>. 
</para>

<informalfigure float="0" id="fig_pathway_diagram">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>the conventional rendering of a G-protein activation pathway</title>
    </objectinfo>
    <imagedata fileref="conventional_rendering.gif"/>
  </imageobject>
</mediaobject>
<caption>A rendering of a G-protein activation pathway. Species are represented by rounded rectangles, and reactions by arrows.</caption>
</informalfigure>

<para>
In CellML, models are thought of as connected networks of discrete components. These components may correspond to physiologically separated regions or chemically distinct objects, or may be useful modelling abstractions. This pathway model has 17 components representing chemically distinct objects (11 chemical species and six reactions) and one component for modelling convenience in which global variables such as time are declared. The CellML rendering of the G-protein activation pathway model is shown in <xref linkend="fig_cellml_rendering"/> (the different shapes in the diagram are explained in the <ulink url="${HTML_EXMPL_GRAPHICAL_NOTATION}">notation guide</ulink>). The action of a catalyst on a reaction is again represented by dashed lines. 
</para>

<informalfigure float="0" id="fig_cellml_rendering">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>the cellml rendering of a G-protein activation pathway</title>
    </objectinfo>
    <imagedata fileref="cellml_rendering.gif"/>
  </imageobject>
</mediaobject>
<caption>The CellML rendering of a G-protein activation pathway.</caption>
</informalfigure>

</sect1>
</article>
</documentation>
  
  
  <!--
    We start the model definition with a definition of some named
    sets of units for use throughout the model.
  -->
  <units name="molar">
    <unit units="mole"/>
    <unit units="litre" exponent="-1"/>
  </units>

  <units name="flux">
    <unit units="molar"/>
    <unit units="second" exponent="-1"/>
  </units>

  <units name="first_order_rate_constant_units">
    <unit units="second" exponent="-1"/>
  </units>

  <units name="second_order_rate_constant_units">
    <unit units="molar" exponent="-1"/>
    <unit units="second" exponent="-1"/>
  </units>
  
  <component name="global_variables">
    <variable units="second" public_interface="out" name="time"/>
  </component>

  <component cmeta:id="Pi" name="Pi">
    
    <variable units="molar" public_interface="out" name="Pi"/>
    <variable units="flux" public_interface="in" name="delta_Pi_rxn5"/>
    <variable units="second" public_interface="in" name="time"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>Pi</ci>
        </apply>
        <ci>delta_Pi_rxn5</ci>
      </apply>
    </math>
  </component>

  <component cmeta:id="L" name="L">
    
    <variable units="molar" public_interface="out" name="L" initial_value="1.0"/>
    <variable units="flux" public_interface="in" name="delta_L_rxn0"/>
    <variable units="second" public_interface="in" name="time"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>L</ci>
        </apply>
        <ci>delta_L_rxn0</ci>
      </apply>
    </math>
  </component>

  <component cmeta:id="R" name="R">
    
     <variable units="molar" public_interface="out" name="R" initial_value="1.0"/>
     <variable units="flux" public_interface="in" name="delta_R_rxn0"/>
     <variable units="second" public_interface="in" name="time"/>
     <math xmlns="http://www.w3.org/1998/Math/MathML">
       <apply>
        <eq/>
         <apply>
          <diff/>
           <bvar>
             <ci>time</ci>
           </bvar>
           <ci>R</ci>
         </apply>
         <ci>delta_R_rxn0</ci>
       </apply>
     </math>
  </component>

  <component cmeta:id="C" name="C">
    
    <variable units="molar" public_interface="out" name="C"/>
    <variable units="flux" public_interface="in" name="delta_C_rxn0"/>
    <variable units="flux" public_interface="in" name="delta_C_rxn4"/>
    <variable units="flux" public_interface="in" name="delta_C_rxn1"/>
    <variable units="second" public_interface="in" name="time"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>C</ci>
        </apply>
        <apply>
          <plus/>
          <ci>delta_C_rxn0</ci>
          <ci>delta_C_rxn4</ci>
          <ci>delta_C_rxn1</ci>
        </apply>
      </apply>
    </math>
  </component>

  <component cmeta:id="G_GDP" name="G_GDP">
    
    <variable units="molar" public_interface="out" name="G_GDP"/>
    <variable units="flux" public_interface="in" name="delta_G_GDP_rxn5"/>
    <variable units="flux" public_interface="in" name="delta_G_GDP_rxn1"/>
    <variable units="second" public_interface="in" name="time"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>G_GDP</ci>
        </apply>
        <apply>
          <plus/>
          <ci>delta_G_GDP_rxn5</ci>
          <ci>delta_G_GDP_rxn1</ci>
        </apply>
      </apply>
    </math>
  </component>

  <component cmeta:id="C_G_GDP" name="C_G_GDP">
    
    <variable units="molar" public_interface="out" name="C_G_GDP"/>
    <variable units="flux" public_interface="in" name="delta_C_G_GDP_rxn1"/>
    <variable units="flux" public_interface="in" name="delta_C_G_GDP_rxn2"/>
    <variable units="second" public_interface="in" name="time"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>C_G_GDP</ci>
        </apply>
        <apply>
          <plus/>
          <ci>delta_C_G_GDP_rxn1</ci>
          <ci>delta_C_G_GDP_rxn2</ci>
        </apply>
      </apply>
    </math>
  </component>

  <component cmeta:id="GDP" name="GDP">
    
    <variable units="molar" public_interface="out" name="GDP"/>
    <variable units="flux" public_interface="in" name="delta_GDP_rxn2"/>
    <variable units="second" public_interface="in" name="time"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>GDP</ci>
        </apply>
        <ci>delta_GDP_rxn2</ci>
      </apply>
    </math>
  </component>

  <component cmeta:id="C_G" name="C_G">
    
    <variable units="molar" public_interface="out" name="C_G"/>
    <variable units="flux" public_interface="in" name="delta_C_G_rxn2"/>
    <variable units="flux" public_interface="in" name="delta_C_G_rxn3"/>
    <variable units="second" public_interface="in" name="time"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>C_G</ci>
        </apply>
        <apply>
          <plus/>
          <ci>delta_C_G_rxn2</ci>
          <ci>delta_C_G_rxn3</ci>
        </apply>
      </apply>
    </math>
  </component>

  <component cmeta:id="GTP" name="GTP">
    
    <variable units="molar" public_interface="out" name="GTP" initial_value="1.0"/>
    <variable units="flux" public_interface="in" name="delta_GTP_rxn3"/>
    <variable units="second" public_interface="in" name="time"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>GTP</ci>
        </apply>
        <ci>delta_GTP_rxn3</ci>
      </apply>
    </math>
  </component>

  <component cmeta:id="C_G_GTP" name="C_G_GTP">
    
    <variable units="molar" public_interface="out" name="C_G_GTP"/>
    <variable units="flux" public_interface="in" name="delta_C_G_GTP_rxn3"/>
    <variable units="flux" public_interface="in" name="delta_C_G_GTP_rxn4"/>
    <variable units="second" public_interface="in" name="time"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>C_G_GTP</ci>
        </apply>
        <apply>
          <plus/>
          <ci>delta_C_G_GTP_rxn3</ci>
          <ci>delta_C_G_GTP_rxn4</ci>
        </apply>
      </apply>
    </math>
  </component>

  <component cmeta:id="G__GTP" name="G__GTP">
    
    <variable units="molar" public_interface="out" name="G__GTP"/>
    <variable units="flux" public_interface="in" name="delta_G__GTP_rxn4"/>
    <variable units="flux" public_interface="in" name="delta_G__GTP_rxn5"/>
    <variable units="second" public_interface="in" name="time"/>
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>G__GTP</ci>
        </apply>
        <apply>
          <plus/>
          <ci>delta_G__GTP_rxn4</ci>
          <ci>delta_G__GTP_rxn5</ci>
        </apply>
      </apply>
    </math>
  </component>

  <component name="reaction0">
    <variable units="molar" public_interface="in" name="R"/>
    <variable units="molar" public_interface="in" name="L"/>
    <variable units="molar" public_interface="in" name="C"/>
    <variable units="flux" public_interface="out" name="delta_R_rxn0"/>
    <variable units="flux" public_interface="out" name="delta_L_rxn0"/>
    <variable units="flux" public_interface="out" name="delta_C_rxn0"/>
    <variable units="second_order_rate_constant_units" name="k0" initial_value="5000000.0"/>
    <variable units="first_order_rate_constant_units" name="k0_" initial_value="0.5"/>
    <variable units="flux" name="rate"/>
    <reaction reversible="yes">
      <variable_ref variable="R">
        <role stoichiometry="1" direction="forward" delta_variable="delta_R_rxn0" role="reactant"/>
      </variable_ref>
      <variable_ref variable="L">
        <role stoichiometry="1" direction="forward" delta_variable="delta_L_rxn0" role="reactant"/>
      </variable_ref>
      <variable_ref variable="C">
        <role stoichiometry="1" direction="forward" delta_variable="delta_C_rxn0" role="product"/>
      </variable_ref>
      <variable_ref variable="rate">
        <role role="rate">
          <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
              <eq/>
              <ci>rate</ci>
              <apply>
                <plus/>
                <apply>
                  <times/>
                  <ci>k0</ci>
                  <ci>R</ci>
                  <ci>L</ci>
                </apply>
                <apply>
                  <minus/>
                  <apply>
                    <times/>
                    <ci>k0_</ci>
                    <ci>C</ci>
                  </apply>
                </apply>
              </apply>
            </apply>
          </math>
        </role>
      </variable_ref>
    </reaction>
  </component>

  <component name="reaction1">
    <variable units="molar" public_interface="in" name="G_GDP"/>
    <variable units="molar" public_interface="in" name="C"/>
    <variable units="molar" public_interface="in" name="C_G_GDP"/>
    <variable units="flux" public_interface="out" name="delta_G_GDP_rxn1"/>
    <variable units="flux" public_interface="out" name="delta_C_rxn1"/>
    <variable units="flux" public_interface="out" name="delta_C_G_GDP_rxn1"/>
    <variable units="second_order_rate_constant_units" name="k1" initial_value="0.1"/>
    <variable units="first_order_rate_constant_units" name="k1_" initial_value="0.1"/>
    <variable units="flux" name="rate"/>
    <reaction reversible="yes">
      <variable_ref variable="G_GDP">
        <role stoichiometry="1" direction="forward" delta_variable="delta_G_GDP_rxn1" role="reactant"/>
      </variable_ref>
      <variable_ref variable="C">
        <role stoichiometry="1" direction="forward" delta_variable="delta_C_rxn1" role="reactant"/>
      </variable_ref>
      <variable_ref variable="C_G_GDP">
        <role stoichiometry="1" direction="forward" delta_variable="delta_C_G_GDP_rxn1" role="product"/>
      </variable_ref>
      <variable_ref variable="rate">
        <role role="rate">
          <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
              <eq/>
              <ci>rate</ci>
              <apply>
                <plus/>
                <apply>
                  <times/>
                  <ci>k1</ci>
                  <ci>G_GDP</ci>
                  <ci>C</ci>
                </apply>
                <apply>
                  <minus/>
                  <apply>
                    <times/>
                    <ci>k1_</ci>
                    <ci>C_G_GDP</ci>
                  </apply>
                </apply>
              </apply>
            </apply>
          </math>
        </role>
      </variable_ref>
    </reaction>
  </component>

  <component name="reaction2">
    <variable units="molar" public_interface="in" name="C_G_GDP"/>
    <variable units="molar" public_interface="in" name="GDP"/>
    <variable units="molar" public_interface="in" name="C_G"/>
    <variable units="flux" public_interface="out" name="delta_C_G_GDP_rxn2"/>
    <variable units="flux" public_interface="out" name="delta_GDP_rxn2"/>
    <variable units="flux" public_interface="out" name="delta_C_G_rxn2"/>
    <variable units="first_order_rate_constant_units" name="k2" initial_value="5.0"/>
    <variable units="second_order_rate_constant_units" name="k2_" initial_value="0.0001"/>
    <variable units="flux" name="rate"/>
    <reaction reversible="yes">
      <variable_ref variable="C_G_GDP">
        <role stoichiometry="1" direction="forward" delta_variable="delta_C_G_GDP_rxn2" role="reactant"/>
      </variable_ref>
      <variable_ref variable="GDP">
        <role stoichiometry="1" direction="forward" delta_variable="delta_GDP_rxn2" role="product"/>
      </variable_ref>
      <variable_ref variable="C_G">
        <role stoichiometry="1" direction="forward" delta_variable="delta_C_G_rxn2" role="product"/>
      </variable_ref>
      <variable_ref variable="rate">
        <role role="rate">
          <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
              <eq/>
              <ci>rate</ci>
              <apply>
                <plus/>
                <apply>
                  <times/>
                  <ci>k2</ci>
                  <ci>C_G_GDP</ci>
                </apply>
                <apply>
                  <minus/>
                  <apply>
                    <times/>
                    <ci>k2_</ci>
                    <ci>GDP</ci>
                    <ci>C_G</ci>
                  </apply>
                </apply>
              </apply>
            </apply>
          </math>
        </role>
      </variable_ref>
    </reaction>
  </component>

  <component name="reaction3">
    <variable units="molar" public_interface="in" name="GTP"/>
    <variable units="molar" public_interface="in" name="C_G"/>
    <variable units="molar" public_interface="in" name="C_G_GTP"/>
    <variable units="flux" public_interface="out" name="delta_GTP_rxn3"/>
    <variable units="flux" public_interface="out" name="delta_C_G_rxn3"/>
    <variable units="flux" public_interface="out" name="delta_C_G_GTP_rxn3"/>
    <variable units="second_order_rate_constant_units" name="k3" initial_value="5000000.0"/>
    <variable units="first_order_rate_constant_units" name="k3_" initial_value="1.0"/>
    <variable units="flux" name="rate"/>
    <reaction reversible="yes">
      <variable_ref variable="GTP">
        <role stoichiometry="1" direction="forward" delta_variable="delta_GTP_rxn3" role="reactant"/>
      </variable_ref>
      <variable_ref variable="C_G">
        <role stoichiometry="1" direction="forward" delta_variable="delta_C_G_rxn3" role="reactant"/>
      </variable_ref>
      <variable_ref variable="C_G_GTP">
        <role stoichiometry="1" direction="forward" delta_variable="delta_C_G_GTP_rxn3" role="product"/>
      </variable_ref>
      <variable_ref variable="rate">
        <role role="rate">
          <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
              <eq/>
              <ci>rate</ci>
              <apply>
                <plus/>
                <apply>
                  <times/>
                  <ci>k3</ci>
                  <ci>GTP</ci>
                  <ci>C_G</ci>
                </apply>
                <apply>
                  <minus/>
                  <apply>
                    <times/>
                    <ci>k3_</ci>
                    <ci>C_G_GTP</ci>
                  </apply>
                </apply>
              </apply>
            </apply>
          </math>
        </role>
      </variable_ref>
    </reaction>
  </component>

  <component name="reaction4">
    <variable units="molar" public_interface="in" name="C_G_GTP"/>
    <variable units="molar" public_interface="in" name="G__GTP"/>
    <variable units="molar" public_interface="in" name="C"/>
    <variable units="flux" public_interface="out" name="delta_C_G_GTP_rxn4"/>
    <variable units="flux" public_interface="out" name="delta_G__GTP_rxn4"/>
    <variable units="flux" public_interface="out" name="delta_C_rxn4"/>
    <variable units="first_order_rate_constant_units" name="k4" initial_value="0.1"/>
    <variable units="flux" name="rate"/>
    <reaction reversible="no">
      <variable_ref variable="C_G_GTP">
        <role stoichiometry="1" delta_variable="delta_C_G_GTP_rxn4" role="reactant"/>
      </variable_ref>
      <variable_ref variable="G__GTP">
        <role stoichiometry="1" delta_variable="delta_G__GTP_rxn4" role="product"/>
      </variable_ref>
      <variable_ref variable="C">
        <role stoichiometry="1" delta_variable="delta_C_rxn4" role="product"/>
      </variable_ref>
      <variable_ref variable="rate">
        <role role="rate">
          <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
              <eq/>
              <ci>rate</ci>
              <apply>
                <times/>
                <ci>k4</ci>
                <ci>C_G_GTP</ci>
              </apply>
            </apply>
          </math>
        </role>
      </variable_ref>
    </reaction>
  </component>

  <component name="reaction5">
    <variable units="molar" public_interface="in" name="G__GTP"/>
    <variable units="molar" public_interface="in" name="Pi"/>
    <variable units="molar" public_interface="in" name="G_GDP"/>
    <variable units="flux" public_interface="out" name="delta_G__GTP_rxn5"/>
    <variable units="flux" public_interface="out" name="delta_Pi_rxn5"/>
    <variable units="flux" public_interface="out" name="delta_G_GDP_rxn5"/>
    <variable units="first_order_rate_constant_units" name="k5" initial_value="2.0"/>
    <variable units="flux" name="rate"/>
    <reaction reversible="no">
      <variable_ref variable="G__GTP">
        <role stoichiometry="1" delta_variable="delta_G__GTP_rxn5" role="reactant"/>
      </variable_ref>
      <variable_ref variable="Pi">
        <role stoichiometry="1" delta_variable="delta_Pi_rxn5" role="product"/>
      </variable_ref>
      <variable_ref variable="G_GDP">
        <role stoichiometry="1" delta_variable="delta_G_GDP_rxn5" role="product"/>
      </variable_ref>
      <variable_ref variable="rate">
        <role role="rate">
          <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
              <eq/>
              <ci>rate</ci>
              <apply>
                <times/>
                <ci>k5</ci>
                <ci>G__GTP</ci>
              </apply>
            </apply>
          </math>
        </role>
      </variable_ref>
    </reaction>
  </component>

  <connection>
    <map_components component_2="reaction5" component_1="Pi"/>
    <map_variables variable_2="Pi" variable_1="Pi"/>
    <map_variables variable_2="delta_Pi_rxn5" variable_1="delta_Pi_rxn5"/>
  </connection>

  <connection>
    <map_components component_2="reaction0" component_1="L"/>
    <map_variables variable_2="L" variable_1="L"/>
    <map_variables variable_2="delta_L_rxn0" variable_1="delta_L_rxn0"/>
  </connection>

  <connection>
    <map_components component_2="reaction0" component_1="R"/>
    <map_variables variable_2="R" variable_1="R"/>
    <map_variables variable_2="delta_R_rxn0" variable_1="delta_R_rxn0"/>
  </connection>

  <connection>
    <map_components component_2="reaction0" component_1="C"/>
    <map_variables variable_2="C" variable_1="C"/>
    <map_variables variable_2="delta_C_rxn0" variable_1="delta_C_rxn0"/>
  </connection>

  <connection>
    <map_components component_2="reaction4" component_1="C"/>
    <map_variables variable_2="C" variable_1="C"/>
    <map_variables variable_2="delta_C_rxn4" variable_1="delta_C_rxn4"/>
  </connection>

  <connection>
    <map_components component_2="reaction1" component_1="C"/>
    <map_variables variable_2="C" variable_1="C"/>
    <map_variables variable_2="delta_C_rxn1" variable_1="delta_C_rxn1"/>
  </connection>

  <connection>
    <map_components component_2="reaction5" component_1="G_GDP"/>
    <map_variables variable_2="G_GDP" variable_1="G_GDP"/>
    <map_variables variable_2="delta_G_GDP_rxn5" variable_1="delta_G_GDP_rxn5"/>
  </connection>

  <connection>
    <map_components component_2="reaction1" component_1="G_GDP"/>
    <map_variables variable_2="G_GDP" variable_1="G_GDP"/>
    <map_variables variable_2="delta_G_GDP_rxn1" variable_1="delta_G_GDP_rxn1"/>
  </connection>

  <connection>
    <map_components component_2="reaction1" component_1="C_G_GDP"/>
    <map_variables variable_2="C_G_GDP" variable_1="C_G_GDP"/>
    <map_variables variable_2="delta_C_G_GDP_rxn1" variable_1="delta_C_G_GDP_rxn1"/>
  </connection>

  <connection>
    <map_components component_2="reaction2" component_1="C_G_GDP"/>
    <map_variables variable_2="C_G_GDP" variable_1="C_G_GDP"/>
    <map_variables variable_2="delta_C_G_GDP_rxn2" variable_1="delta_C_G_GDP_rxn2"/>
  </connection>

  <connection>
    <map_components component_2="reaction2" component_1="GDP"/>
    <map_variables variable_2="GDP" variable_1="GDP"/>
    <map_variables variable_2="delta_GDP_rxn2" variable_1="delta_GDP_rxn2"/>
  </connection>

  <connection>
    <map_components component_2="reaction2" component_1="C_G"/>
    <map_variables variable_2="C_G" variable_1="C_G"/>
    <map_variables variable_2="delta_C_G_rxn2" variable_1="delta_C_G_rxn2"/>
  </connection>

  <connection>
    <map_components component_2="reaction3" component_1="C_G"/>
    <map_variables variable_2="C_G" variable_1="C_G"/>
    <map_variables variable_2="delta_C_G_rxn3" variable_1="delta_C_G_rxn3"/>
  </connection>

  <connection>
    <map_components component_2="reaction3" component_1="GTP"/>
    <map_variables variable_2="GTP" variable_1="GTP"/>
    <map_variables variable_2="delta_GTP_rxn3" variable_1="delta_GTP_rxn3"/>
  </connection>

  <connection>
    <map_components component_2="reaction3" component_1="C_G_GTP"/>
    <map_variables variable_2="C_G_GTP" variable_1="C_G_GTP"/>
    <map_variables variable_2="delta_C_G_GTP_rxn3" variable_1="delta_C_G_GTP_rxn3"/>
  </connection>

  <connection>
    <map_components component_2="reaction4" component_1="C_G_GTP"/>
    <map_variables variable_2="C_G_GTP" variable_1="C_G_GTP"/>
    <map_variables variable_2="delta_C_G_GTP_rxn4" variable_1="delta_C_G_GTP_rxn4"/>
  </connection>

  <connection>
    <map_components component_2="reaction4" component_1="G__GTP"/>
    <map_variables variable_2="G__GTP" variable_1="G__GTP"/>
    <map_variables variable_2="delta_G__GTP_rxn4" variable_1="delta_G__GTP_rxn4"/>
  </connection>

  <connection>
    <map_components component_2="reaction5" component_1="G__GTP"/>
    <map_variables variable_2="G__GTP" variable_1="G__GTP"/>
    <map_variables variable_2="delta_G__GTP_rxn5" variable_1="delta_G__GTP_rxn5"/>
  </connection>

  <connection>
   <map_components component_2="global_variables" component_1="Pi"/>
   <map_variables variable_2="time" variable_1="time"/>
  </connection>

  <connection>
    <map_components component_2="global_variables" component_1="L"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>

  <connection>
    <map_components component_2="global_variables" component_1="R"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>

  <connection>
    <map_components component_2="global_variables" component_1="C"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>

  <connection>
    <map_components component_2="global_variables" component_1="G_GDP"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>

  <connection>
    <map_components component_2="global_variables" component_1="C_G_GDP"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>

  <connection>
    <map_components component_2="global_variables" component_1="GDP"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>

  <connection>
    <map_components component_2="global_variables" component_1="C_G"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>

  <connection>
    <map_components component_2="global_variables" component_1="GTP"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>

  <connection>
    <map_components component_2="global_variables" component_1="C_G_GTP"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>

  <connection>
    <map_components component_2="global_variables" component_1="G__GTP"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>


<rdf:RDF>
  <rdf:Seq rdf:about="rdf:#88eb9858-6110-4e20-b6a6-c6eb5048635a">
    <rdf:li rdf:resource="rdf:#97574417-3809-43ae-968e-e8a238524b13"/>
    <rdf:li rdf:resource="rdf:#5a17a59e-3caf-4e6e-886b-1f01e2111a04"/>
  </rdf:Seq>
  <rdf:Description rdf:about="#C_G">
    <dcterms:alternative>receptor bound G-protein</dcterms:alternative>
    <dc:title>C_G</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="#R">
    <dcterms:alternative>free receptor</dcterms:alternative>
    <dc:title>R</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#c0d4e3e3-c4ec-4d0c-b8ed-fda898eabb77">
    <bqs:subject_type>keyword</bqs:subject_type>
    <rdf:value>signal transduction</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#C">
    <dcterms:alternative>receptor-ligand complex</dcterms:alternative>
    <dc:title>C</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="#L">
    <dcterms:alternative>free ligand</dcterms:alternative>
    <dc:title>L</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5a17a59e-3caf-4e6e-886b-1f01e2111a04">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#190369a6-263d-4804-9293-32e7227aebd0"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#cdc0b5b9-7de6-423c-9a6d-0e4ab161691b">
    <dcterms:W3CDTF>2001-11-14</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#1e4bcba5-2332-4a76-8bc5-bb9fea9296aa">
    <vCard:Given>W</vCard:Given>
    <vCard:Family>Thomsen</vCard:Family>
    <vCard:Other>J</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="#G__GTP">
    <dcterms:alternative>activated G-protein</dcterms:alternative>
    <dc:title>G__GTP</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="#GTP">
    <dcterms:alternative>guanosine triphosphate</dcterms:alternative>
    <dc:title>GTP</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d33bd375-b4fc-4348-8472-0e93e0997aac">
    <dc:creator rdf:resource="rdf:#2540bb20-965e-405e-b335-445805d520fe"/>
    <rdf:value>
          This is the CellML description of Thomsen and Neubig's mathematical 
          model of a G-protein activation pathway (1989).  
        </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#C_G_GDP">
    <dcterms:alternative>receptor bound G-protein-GDP</dcterms:alternative>
    <dc:title>C_G_GDP</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#190369a6-263d-4804-9293-32e7227aebd0">
    <vCard:Given>R</vCard:Given>
    <vCard:Family>Neubig</vCard:Family>
    <vCard:Other>R</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="#Pi">
    <dcterms:alternative>inorganic phosphate</dcterms:alternative>
    <dc:title>Pi</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#304b20a9-1e2f-4530-9a1c-4d730c20b649">
    <vCard:Orgname>The University of Auckland</vCard:Orgname>
    <vCard:Orgunit>The Bioengineering Research Group</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b30c1e58-0768-479f-80a5-2e81954a240e">
    <vCard:N rdf:resource="rdf:#897b445d-dabd-4ea8-8ced-b169b91f46f5"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#35357601-8fa9-42a7-969d-fa9df34a3fa6">
    <dc:creator rdf:resource="rdf:#88eb9858-6110-4e20-b6a6-c6eb5048635a"/>
    <dc:title>
           Rapid kinetics of alpha 2-adrenergic inhibition of adenylate cyclase.            Evidence for a distal rate-limiting step
          </dc:title>
    <bqs:volume>28</bqs:volume>
    <bqs:first_page>8778</bqs:first_page>
    <bqs:Journal rdf:resource="rdf:#d7ab5091-81fd-4468-9c3c-7c0975934295"/>
    <dcterms:issued rdf:resource="rdf:#ff3a88d5-66e3-4390-90aa-a6c933f0eb63"/>
    <bqs:last_page>8786</bqs:last_page>
  </rdf:Description>
  <rdf:Description rdf:about="#thomsen_neubig_1989_version01">
    <dc:title>G-protein activation pathway</dc:title>
    <cmeta:comment rdf:resource="rdf:#d33bd375-b4fc-4348-8472-0e93e0997aac"/>
    <bqs:reference rdf:resource="rdf:#2f4a9788-f7a6-40ea-8e7a-9d016d27fe99"/>
    <bqs:reference rdf:resource="rdf:#4b1239dd-4ba7-43ae-b1f0-dcf2ac614b14"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#82b9ba88-b889-477b-8856-93953dfe38a9">
    <rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
    <rdf:value>c.lloyd@auckland.ac.nz</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="#GDP">
    <dcterms:alternative>guanosine diphosphate</dcterms:alternative>
    <dc:title>GDP</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d7ab5091-81fd-4468-9c3c-7c0975934295">
    <dc:title>Biochemistry</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="#C_G_GTP">
    <dcterms:alternative>receptor bound G-protein-GTP</dcterms:alternative>
    <dc:title>C_G_GTP</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#897b445d-dabd-4ea8-8ced-b169b91f46f5">
    <vCard:Given>Autumn</vCard:Given>
    <vCard:Family>Cuellar</vCard:Family>
    <vCard:Other>A.</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#15c67f6e-2c3f-40e9-9543-8b9c792c3d65">
    <dcterms:modified rdf:resource="rdf:#8735ef6d-62ff-4237-bdf3-e3b1bffd1121"/>
    <rdf:value>
          Updated metadata to conform to the 16/1/02 CellML Metadata 1.0 
          Specification.
        </rdf:value>
    <cmeta:modifier rdf:resource="rdf:#b30c1e58-0768-479f-80a5-2e81954a240e"/>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dc:publisher>
        The University of Auckland, Bioengineering Research Group
      </dc:publisher>
    <cmeta:modification rdf:resource="rdf:#15c67f6e-2c3f-40e9-9543-8b9c792c3d65"/>
    <dcterms:created rdf:resource="rdf:#cdc0b5b9-7de6-423c-9a6d-0e4ab161691b"/>
    <dc:creator rdf:resource="rdf:#b510803b-d23a-49a2-9fd2-9e2bb3efe754"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#97574417-3809-43ae-968e-e8a238524b13">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#1e4bcba5-2332-4a76-8bc5-bb9fea9296aa"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4b1239dd-4ba7-43ae-b1f0-dcf2ac614b14">
    <dc:subject rdf:resource="rdf:#c0d4e3e3-c4ec-4d0c-b8ed-fda898eabb77"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#2540bb20-965e-405e-b335-445805d520fe">
    <vCard:FN>Catherine Lloyd</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b510803b-d23a-49a2-9fd2-9e2bb3efe754">
    <vCard:ORG rdf:resource="rdf:#304b20a9-1e2f-4530-9a1c-4d730c20b649"/>
    <vCard:EMAIL rdf:resource="rdf:#82b9ba88-b889-477b-8856-93953dfe38a9"/>
    <vCard:N rdf:resource="rdf:#7bbc3394-2fee-44a7-935e-75c75c30643f"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8735ef6d-62ff-4237-bdf3-e3b1bffd1121">
    <dcterms:W3CDTF>2002-01-21</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="#G_GDP">
    <dcterms:alternative>G-protein-GDP</dcterms:alternative>
    <dcterms:alternative>inactivated G-protein</dcterms:alternative>
    <dc:title>G_GDP</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ff3a88d5-66e3-4390-90aa-a6c933f0eb63">
    <dcterms:W3CDTF>1989-10-31</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#7bbc3394-2fee-44a7-935e-75c75c30643f">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#2f4a9788-f7a6-40ea-8e7a-9d016d27fe99">
    <bqs:Pubmed_id>2574993</bqs:Pubmed_id>
    <bqs:JournalArticle rdf:resource="rdf:#35357601-8fa9-42a7-969d-fa9df34a3fa6"/>
  </rdf:Description>
</rdf:RDF>
</model>