Location: Goodwin, 1965 @ 7ab656a4e509 / goodwin_1965_c.cellml

Author:
Catherine Lloyd <c.lloyd@auckland.ac.nz>
Date:
2010-08-06 20:31:06+12:00
Desc:
Fixed broken annotations and added keyword to model c.
Permanent Source URI:
http://models.cellml.org/workspace/goodwin_1965/rawfile/7ab656a4e5096d1341b3313d6e6a0e5efcb70442/goodwin_1965_c.cellml

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


--><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:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" xmlns:ns7="http://www.cellml.org/metadata/simulation/1.0#" cmeta:id="goodwin_model_1965" name="goodwin_model_1965">

<documentation xmlns="http://cellml.org/tmp-documentation">
    <article>
      <section id="sec_status">
        <title>Model Status</title>
        <para>
         This CellML model runs in both OpenCell and COR, however it does not recreate the published results - the CellML model does not oscillate.  The units have been checked and they are consistent.  This particular version of the model describes the third control system of a circuit with limit cycles.  We suspect the CellML model does not oscillate because it does not include time delays.  These are mentioned in the main body of text in the published paper, but they are not descibed mathematically in the model equations.
          </para>
      </section>
      <sect1 id="sec_structure">
        <title>Model Structure</title>
        <para>
ABSTRACT: THE demonstration of negative feedback control processes operating at the molecular level in cells is one of the most significant developments in modern biology. The phenomena of feedback inhibition and feedback repression, whereby enzymatic activities are controlled at the level of the enzyme and the gene, respectively, provide a firm experimental basis for the construction of dynamic models which represent the fundamental regulatory activity of cells. The behavior of these and other molecular control circuits thus constitutes the basis of cell physiology, and in effect provides the physiologist with his elementary units of function. The purpose of this paper is to illustrate the type of periodic behavior which can arise in model systems incorporating the essential control features of enzymatic regulatory processes, and to discuss the significance of oscillatory motion in relation to the organization of cellular processes in time.  
</para>

        <para>
The complete original paper reference is cited below:
</para>
        <para>
          Oscillatory behavior in enzymatic control processes, Brian C. Goodwin, 1965, <emphasis>Advances in Enzyme Regulation</emphasis>, 3, 425-438. (An abstract and a PDF version of the article are available to subscribers on the journal website.)  <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&amp;amp;db=PubMed&amp;amp;list_uids=5861813&amp;amp;dopt=Abstract">PubMed ID: 5861813</ulink>
        </para>
		
        <informalfigure float="0" id="fig_reaction_diagram">
          <mediaobject>
            <imageobject>
              <objectinfo>
                <title>cell diagram</title>
              </objectinfo>
              <imagedata fileref="goodwin_1965b.png"/>
            </imageobject>
          </mediaobject>
          <caption>Schematic diagram of the third control system modelled in this study: a control circuit with limit-cycle characteristics.</caption>
        </informalfigure>
      </sect1>
    </article>
  </documentation>
  
  
  <units name="first_order_rate_constant">
    <unit units="second" exponent="-1"/>
  </units>
  
  <units name="flux">
    <unit units="nanomolar"/>
    <unit units="second" exponent="-1"/>
  </units>
  
  <units name="nanomolar">
    <unit units="mole" prefix="nano"/>
    <unit units="litre" exponent="-1"/>
  </units>
  
  <units name="per_nanomolar">
    <unit units="nanomolar" exponent="-1"/>
  </units>

  
  <component name="environment">
    <variable units="second" public_interface="out" name="time" cmeta:id="environment_time"/>
  </component>
  
  
  <component cmeta:id="X1" name="X1">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="#X1">
        <dc:title>X1</dc:title>
        <dcterms:alternative>concentration of the nuclear messenger</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF>
    <variable units="nanomolar" public_interface="out" name="X1" initial_value="7"/>
	
	<variable units="flux" name="a1" initial_value="360"/>
	<variable units="first_order_rate_constant" name="b1" initial_value="1"/>
	<variable units="dimensionless" name="A1" initial_value="43"/>
	<variable units="per_nanomolar" name="k1" initial_value="1"/>
	
	<variable units="nanomolar" public_interface="in" name="Z1"/>
	<variable units="second" public_interface="in" name="time"/>
	
	<math xmlns="http://www.w3.org/1998/Math/MathML" id="1"> 
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>X1</ci>
        </apply> 
        <apply>
          <minus/>
          <apply>
            <divide/>
            <ci> a1 </ci>
		    <apply>
              <plus/>
              <ci> A1 </ci>
			  <apply>
                <times/>
			    <ci> k1 </ci>
		        <ci> Z1 </ci>
	          </apply>
			</apply>
		  </apply>
		  <apply>
            <times/>
			<ci> b1 </ci>
		    <ci> X1 </ci>
	      </apply>
		</apply>
      </apply>
    </math>
  </component>
  
  <component cmeta:id="Y1" name="Y1">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="#Y1">
        <dc:title>Y1</dc:title>
        <dcterms:alternative>concentration of the active cytoplasmic messenger</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF>
    <variable units="nanomolar" public_interface="out" name="Y1" initial_value="0"/>
	
	<variable units="first_order_rate_constant" name="beta_1" initial_value="0.6"/>
	<variable units="first_order_rate_constant" name="alpha_1" initial_value="1"/>
	
	<variable units="nanomolar" public_interface="in" name="X1"/>
	<variable units="second" public_interface="in" name="time"/>
	
	<math xmlns="http://www.w3.org/1998/Math/MathML" id="2"> 
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>Y1</ci>
        </apply> 
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> alpha_1 </ci>
		    <ci> X1 </ci>
		  </apply>
		  <apply>
            <times/>
            <ci> beta_1 </ci>
		    <ci> Y1 </ci>
		  </apply>
	    </apply>
      </apply>
    </math>
  </component>
  
  <component cmeta:id="Z1" name="Z1">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="#Z1">
        <dc:title>Z1</dc:title>
        <dcterms:alternative>concentration of the active enzyme</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF>
    <variable units="nanomolar" public_interface="out" name="Z1" initial_value="0"/>
	
	<variable units="first_order_rate_constant" name="gamma_1" initial_value="1"/>
	<variable units="first_order_rate_constant" name="delta_1" initial_value="0.8"/>
	
	<variable units="nanomolar" public_interface="in" name="Y1"/>
	<variable units="second" public_interface="in" name="time"/>
	
	<math xmlns="http://www.w3.org/1998/Math/MathML" id="3"> 
      <apply>
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci>time</ci>
          </bvar>
          <ci>Z1</ci>
        </apply> 
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> gamma_1 </ci>
		    <ci> Y1 </ci>
		  </apply>
		  <apply>
            <times/>
            <ci> delta_1 </ci>
		    <ci> Z1 </ci>
		  </apply>
	    </apply>
      </apply>
    </math>
  </component>
  
  <!-- connections -->
  
  <connection>
    <map_components component_2="environment" component_1="X1"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="Y1"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="Z1"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="Y1" component_1="X1"/>
    <map_variables variable_2="X1" variable_1="X1"/>
  </connection>
  
  <connection>
    <map_components component_2="Z1" component_1="X1"/>
    <map_variables variable_2="Z1" variable_1="Z1"/>
  </connection>
  
  <connection>
    <map_components component_2="Y1" component_1="Z1"/>
    <map_variables variable_2="Y1" variable_1="Y1"/>
  </connection>

<rdf:RDF>
  <rdf:Seq rdf:about="rdf:#c473fde7-55f0-4dca-9253-8f5936d50b56">
    <rdf:li rdf:resource="rdf:#4328cd8a-387e-46fc-84b9-2c3a1ac5fb15"/>
  </rdf:Seq>
  <rdf:Description rdf:about="rdf:#55bf340e-7d28-421c-a3ef-01a77e7dd6bb">
    <dc:creator rdf:resource="rdf:#a6e9c249-6083-4d8f-aeb8-544fa5f4fabf"/>
    <rdf:value>This CellML model runs in both PCEnv and COR, however it does not recreate the published results - the CellML model does not oscillate.  The units have been checked and they are consistent.  This particular version of the model describes the third control system of a circuit with limit cycles.  We suspect the CellML model does not oscillate because it does not include time delays.  These are mentioned in the main body of text in the published paper, but they are not descibed mathematically in the model equations.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#fe15a4d3-9780-4a7e-b709-08a9c359e5ae">
    <dc:creator rdf:resource="rdf:#531107a6-125f-4805-bccf-597e60966898"/>
    <rdf:value/>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dc:title>Oscillatory behavior in enzymatic control processes (Model 3)</dc:title>
	<dc:publisher/>
    <cmeta:comment rdf:resource="rdf:#55bf340e-7d28-421c-a3ef-01a77e7dd6bb"/>
    <dcterms:created rdf:resource="rdf:#02acde91-a931-430e-b152-803fd0398630"/>
    <dc:creator rdf:resource="rdf:#b9f66eee-75d4-47e0-b74f-7969550d8d5a"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#79b7bca5-b041-4eee-9a4f-4e3b0d510521">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#58167fe6-4ed4-4908-89b5-44a0ed80747f">
    <dc:creator rdf:resource="rdf:#c473fde7-55f0-4dca-9253-8f5936d50b56"/>
    <dc:title>Oscillatory behavior in enzymatic control processes</dc:title>
    <bqs:volume>3</bqs:volume>
    <bqs:first_page>425</bqs:first_page>
    <bqs:Journal rdf:resource="rdf:#21772fab-2bae-4975-90b7-736de934f274"/>
    <dcterms:issued rdf:resource="rdf:#49c32b79-2de6-4355-873a-f5ada1cdeb0c"/>
    <bqs:last_page>438</bqs:last_page>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#dd7618ea-3110-456c-94b1-f49e51984b90">
    <vCard:Given>Brian</vCard:Given>
    <vCard:Family>Goodwin</vCard:Family>
    <vCard:Other>C</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#531107a6-125f-4805-bccf-597e60966898">
    <vCard:FN/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#18cfe305-e7d3-4b1a-a59e-70773a8b0524">
    <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="rdf:#02acde91-a931-430e-b152-803fd0398630">
    <dcterms:W3CDTF>2008-09-26T00:00:00+00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="#goodwin_model_1965">
    <bqs:reference rdf:resource="rdf:#8cb1d369-65c5-414b-b706-de3c769eda0b"/>
    <cmeta:comment rdf:resource="rdf:#fe15a4d3-9780-4a7e-b709-08a9c359e5ae"/>
	<bqs:reference xmlns:bqs="http://www.cellml.org/bqs/1.0#" rdf:parseType="Resource">
  <dc:subject xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:parseType="Resource">
    <bqs:subject_type>keyword</bqs:subject_type>
    <rdf:value>
      <rdf:Bag>
        <rdf:li>circadian rhythms</rdf:li>
      </rdf:Bag>
    </rdf:value>
  </dc:subject>
</bqs:reference>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#21772fab-2bae-4975-90b7-736de934f274">
    <dc:title>Advances in Enzyme Regulation</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4328cd8a-387e-46fc-84b9-2c3a1ac5fb15">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#dd7618ea-3110-456c-94b1-f49e51984b90"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#550e6dc6-113a-48ac-abb1-d81329e99a29">
    <vCard:Orgname>The University of Auckland</vCard:Orgname>
    <vCard:Orgunit>Auckland Bioengineering Institute</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b9f66eee-75d4-47e0-b74f-7969550d8d5a">
    <vCard:ORG rdf:resource="rdf:#550e6dc6-113a-48ac-abb1-d81329e99a29"/>
    <vCard:EMAIL rdf:resource="rdf:#18cfe305-e7d3-4b1a-a59e-70773a8b0524"/>
    <vCard:N rdf:resource="rdf:#79b7bca5-b041-4eee-9a4f-4e3b0d510521"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8cb1d369-65c5-414b-b706-de3c769eda0b">
    <bqs:Pubmed_id>5</bqs:Pubmed_id>
    <bqs:JournalArticle rdf:resource="rdf:#58167fe6-4ed4-4908-89b5-44a0ed80747f"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#49c32b79-2de6-4355-873a-f5ada1cdeb0c">
    <dcterms:W3CDTF>1965-00-00 00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a6e9c249-6083-4d8f-aeb8-544fa5f4fabf">
    <vCard:FN>Catherine Lloyd</vCard:FN>
  </rdf:Description>
</rdf:RDF>
</model>