Location: Houart, Dupont and Goldbeter 1999 @ 5a525ffac162 / houart_1999_1.1.cellml

Author:
Mona Zhu <devnull@localhost>
Date:
2011-01-19 16:25:21+13:00
Desc:
Changed import locations
Permanent Source URI:
https://models.physiomeproject.org/workspace/houart_dupont_goldbeter_1999/rawfile/5a525ffac162c4096f3a0fa1bdbcb30723b44f34/houart_1999_1.1.cellml

<?xml version="1.0"?>
<model xmlns="http://www.cellml.org/cellml/1.1#" xmlns:cmeta="http://www.cellml.org/metadata/1.1#" cmeta:id="houart_1999_1.1" name="houart_1999_1.1">
	<!-- documentation -->
	<documentation xmlns="http://cellml.org/tmp-documentation">
		<article>
			<articleinfo>
				<title>Bursting, Chaos and Birhythmicity Originating from Self-modulation of the Inositol 1,4,5-triphosphate Signal in a Model for Intracellular Ca2+ Oscillations</title>
				<author>
					<firstname>Mona</firstname>
					<surname>Zhu</surname>
					<affiliation>
						<shortaffil>Auckland Bioengineering Institute, University of Auckland</shortaffil>
					</affiliation>
				</author>
			</articleinfo>
			<section id="sec_status">
				<title>Model Status</title>
				<para>This model has been built with the differential expressions of Houart's 1999 paper for various types of intracellular Ca2+ oscillations bassed on self-modulation of the Inositol 1,4,5-triphosphate signal. The model is a 1.1 modification of the orignal 1.0 model; it uses imports that are reusable in defining equations in the model. The parameter values for this model are different for different types of intracellular Ca2+ oscillations. The parameter values used for this particular model were taken from Figure 2 on page 513 and reproduces the corresponding curves. Parameter values from pages 515 and 520 may be changed to reflect curves for Bursting, Chaos, Quasiperiodicity and Birhythmicity. The initial conditions have been set after allowing the model to run for approximately 5 minutes until steady state. This file is known to run in OpenCell only.</para>
			</section>
			<sect1 id="sec_structure">
				<title>Model Structure</title>
				<para>We investigate the various types of complex Ca2+ oscillations which can arise in a model based on the mechanism of Ca2+-induced Ca2+ release (CICR), that takes into account the Ca2+-stimulated degradation of inositol 1,4,5-trisphosphate (InsP3) by a 3-kinase. This model was previously proposed in the course of an investigation of plausible mechanisms capable of generating complex Ca2+ oscillations. Besides simple periodic behavior, this model for cytosolic Ca2+ oscillations in nonexcitable cells shows complex oscillatory phenomena like bursting or chaos. We show that the model also admits a coexistence between two stable regimes of sustained oscillations (birhythmicity). The occurrence of these various modes of oscillatory behavior is analysed by means of bifurcation diagrams. Complex oscillations are characterized by means of Poincare sections, power spectra and Lyapounov exponents. The results point to the role of self-modulation of the InsP3 signal by 3-kinase as a possible source for complex temporal patterns in Ca2+ signaling.</para>
				<para>
					Bursting, Chaos and Birhythmicity Originating from Self-modulation of the Inositol 1,4,5-triphosphate Signal in a Model for Intracellular Ca2+ Oscillations, Houart G, Dupont G, Goldbeter A, 1999, <emphasis>Bulletin of Mathematical Biology</emphasis>
					, 61, 507-530 <ulink url="http://www.ncbi.nlm.nih.gov/pubmed/17883229">PubMed ID: 17883229</ulink>
				</para>
				<informalfigure float="0" id="______">
					<mediaobject>
						<imageobject>
							<objectinfo>
								<title>Model Diagram</title>
							</objectinfo>
							<imagedata fileref="houart_1999_1.1.png"/>
						</imageobject>
					</mediaobject>
					<caption>Schematic representation of the model based on the interplay between CICR and the Ca2+-stimulated degradation of InsP3. Besides simple periodic oscillations, this model can produce complex Ca2+ oscillations including bursting, chaos, quasiperiodic behavior, as well as birhythmicity.</caption>
				</informalfigure>
			</sect1>
		</article>
	</documentation>   
	
	<!-- units -->
	<units name="minute">
      <unit multiplier="60" units="second"/>
   </units>
   <units name="per_minute">
      <unit exponent="-1" units="minute"/>
   </units>
   <units name="per_litre">
      <unit exponent="-1" units="litre"/>
   </units>
   <units name="micromolar">
      <unit prefix="micro" units="mole"/>
      <unit units="per_litre"/>
   </units>
   <units name="micromolar_per_minute">
      <unit units="micromolar"/>
      <unit units="per_minute"/>
   </units>
   
   <!-- components -->
   <component name="environment">
      <variable name="time" public_interface="out" units="minute"/>
   </component>
   <component name="Vin">
      <variable initial_value="2" name="V_0" units="micromolar_per_minute"/>
      <variable initial_value="2" name="V_1" units="micromolar_per_minute"/>
      <variable initial_value="0.6" name="beta" public_interface="out" units="dimensionless"/>
      <variable name="V_in" public_interface="out" units="micromolar_per_minute"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
		<apply>
            <eq/>
            <ci>V_in</ci>
            <apply>
               <plus/>
               <ci>V_0</ci>
               <apply>
                  <times/>
                  <ci>V_1</ci>
                  <ci>beta</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="V2">
      <variable initial_value="6" name="V_M2" private_interface="out" units="micromolar_per_minute"/>
      <variable name="Z" private_interface="out" public_interface="in" units="micromolar"/>
      <variable initial_value="0.1" name="K_2" private_interface="out" units="micromolar"/>
      <variable name="V_2" private_interface="in" public_interface="out" units="micromolar_per_minute"/>
      <variable initial_value="0" name="zero" private_interface="out" units="micromolar"/>
	  <variable initial_value="1" name="one_dm" private_interface="out" units="dimensionless"/>
	  <variable initial_value="1" name="one_mm" private_interface="out" units="micromolar"/>
	  <variable initial_value="2" name="two" private_interface="out" units="dimensionless"/>
   </component>
   <component name="V3">
      <variable initial_value="20" name="V_M3" private_interface="out" units="micromolar_per_minute"/>
      <variable initial_value="0.5" name="K_Z" private_interface="out" units="micromolar"/>
      <variable initial_value="0.2" name="K_A" private_interface="out" units="micromolar"/>
      <variable initial_value="0.2" name="K_Y" private_interface="out" units="micromolar"/>
      <variable initial_value="2" name="m" private_interface="out" units="dimensionless"/>
      <variable name="Y" private_interface="out" public_interface="in" units="micromolar"/>	
      <variable name="A" private_interface="out" public_interface="in" units="micromolar"/>
      <variable name="Z" private_interface="out" public_interface="in" units="micromolar"/>
      <variable name="V_3" private_interface="in" public_interface="out" units="micromolar_per_minute"/>
      <variable initial_value="2" name="two" private_interface="out" units="dimensionless"/>
	  <variable initial_value="4" name="four" private_interface="out" units="dimensionless"/>
   </component>
   <component name="V5">
      <variable initial_value="5" name="V_M5" private_interface="out" units="micromolar_per_minute"/>
      <variable initial_value="1" name="K_5" private_interface="out" units="micromolar"/>
      <variable initial_value="0.4" name="K_d" private_interface="out" units="micromolar"/>
      <variable initial_value="2" name="p" private_interface="out" units="dimensionless"/>
      <variable initial_value="4" name="n" private_interface="out" units="dimensionless"/>
      <variable name="A" private_interface="out" public_interface="in" units="micromolar"/>
      <variable name="Z" private_interface="out" public_interface="in" units="micromolar"/>
      <variable name="V_5" private_interface="in" public_interface="out" units="micromolar_per_minute"/>
      <variable initial_value="0" name="zero" private_interface="out" units="micromolar"/>
	  <variable initial_value="1" name="one_dm" private_interface="out" units="dimensionless"/>
	  <variable initial_value="1" name="one_mm" private_interface="out" units="micromolar"/>
   </component>
   <component name="cytosol">
      <variable name="time" private_interface="out" public_interface="in" units="minute"/>
      <variable name="V_in" public_interface="in" units="micromolar_per_minute"/>
      <variable name="V_2" public_interface="in" units="micromolar_per_minute"/>
      <variable name="V_3" public_interface="in" units="micromolar_per_minute"/>
      <variable initial_value="10" name="k" units="per_minute"/>
      <variable initial_value="1" name="k_f" public_interface="out" units="per_minute"/>
      <variable initial_value="0.15" name="Z" public_interface="out" units="micromolar"/>
      <variable name="Y" public_interface="in" units="micromolar"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Z</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <apply>
                     <minus/>
                     <ci>V_in</ci>
                     <ci>V_2</ci>
                  </apply>
                  <ci>V_3</ci>
                  <apply>
                     <times/>
                     <ci>k_f</ci>
                     <ci>Y</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>k</ci>
                  <ci>Z</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="internal_pool">
      <variable name="time" private_interface="out" public_interface="in" units="minute"/>
      <variable name="V_2" private_interface="out" public_interface="in" units="micromolar_per_minute"/>
      <variable name="V_3" private_interface="out" public_interface="in" units="micromolar_per_minute"/>
      <variable name="k_f" private_interface="out" public_interface="in" units="per_minute"/>
      <variable name="Y" private_interface="in" public_interface="out" units="micromolar"/>
      <variable initial_value="1" name="Y_init" private_interface="out" units="micromolar"/>
	  <variable initial_value="1" name="one" private_interface="out" units="dimensionless"/>
   </component>
   <component name="InsP3_conc">
      <variable name="time" private_interface="out" public_interface="in" units="minute"/>
      <variable name="beta" private_interface="out" public_interface="in" units="dimensionless"/>
      <variable initial_value="0.1" name="epsilon" private_interface="out" units="per_minute"/>
      <variable initial_value="2" name="V_4" private_interface="out" units="micromolar_per_minute"/>
      <variable name="V_5" private_interface="out" public_interface="in" units="micromolar_per_minute"/>
      <variable name="A" private_interface="in" public_interface="out" units="micromolar"/>
      <variable initial_value="0.42" name="A_init" private_interface="out" units="micromolar"/>
   </component>

   <!-- groups -->
    <group>
		<relationship_ref relationship="encapsulation"/>
			<component_ref component="V2">
				<component_ref component="V_2"/>
			</component_ref>
	</group>
	<group>
		<relationship_ref relationship="encapsulation"/>
			<component_ref component="V3">
				<component_ref component="V_3"/>
			</component_ref>
	</group>
	<group>
		<relationship_ref relationship="encapsulation"/>
			<component_ref component="V5">
				<component_ref component="V_5"/>
			</component_ref>
	</group>
	<group>
		<relationship_ref relationship="encapsulation"/>
			<component_ref component="internal_pool">
				<component_ref component="internalpool"/>
			</component_ref>
   </group>
   <group>
		<relationship_ref relationship="encapsulation"/>
			<component_ref component="InsP3_conc">
				<component_ref component="InsP3conc"/>
			</component_ref>
   </group>
      
   <!-- imports -->
   <import xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href= "constitutive_law1.cellml">
		<component component_ref="interface" name="V_2"/>
   </import>
   <import xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="constitutive_law1.cellml">
	   <component component_ref="interface" name="V_3"/>
   </import>
   <import xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="constitutive_law1.cellml">
	   <component component_ref="interface" name="V_5"/>
   </import>
   <import xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="constitutive_law2.cellml">
	   <component component_ref="interface" name="internalpool"/>
   </import>
   <import xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="constitutive_law2.cellml">
	   <component component_ref="interface" name="InsP3conc"/>
   </import>  

   <!-- connections -->
   <connection>
      <map_components component_1="Vin" component_2="cytosol"/>
      <map_variables variable_1="V_in" variable_2="V_in"/>
   </connection>
   <connection>
      <map_components component_1="Vin" component_2="InsP3_conc"/>
      <map_variables variable_1="beta" variable_2="beta"/>
   </connection>
   <connection>
      <map_components component_1="V2" component_2="cytosol"/>
      <map_variables variable_1="V_2" variable_2="V_2"/>
      <map_variables variable_1="Z" variable_2="Z"/>
   </connection>
   <connection>
      <map_components component_1="V2" component_2="internal_pool"/>
      <map_variables variable_1="V_2" variable_2="V_2"/>
   </connection>
   <connection>
      <map_components component_1="V3" component_2="cytosol"/>
      <map_variables variable_1="V_3" variable_2="V_3"/>
      <map_variables variable_1="Z" variable_2="Z"/>
   </connection>
   <connection>
      <map_components component_1="V3" component_2="internal_pool"/>
      <map_variables variable_1="Y" variable_2="Y"/>
      <map_variables variable_1="V_3" variable_2="V_3"/>
   </connection>
   <connection>
      <map_components component_1="V3" component_2="InsP3_conc"/>
      <map_variables variable_1="A" variable_2="A"/>
   </connection>
   <connection>
      <map_components component_1="V5" component_2="InsP3_conc"/>
      <map_variables variable_1="A" variable_2="A"/>
      <map_variables variable_1="V_5" variable_2="V_5"/>
   </connection>
   <connection>
      <map_components component_1="environment" component_2="cytosol"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="environment" component_2="internal_pool"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="environment" component_2="InsP3_conc"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="cytosol" component_2="V5"/>
      <map_variables variable_1="Z" variable_2="Z"/>
   </connection>
   <connection>
      <map_components component_1="internal_pool" component_2="cytosol"/>
      <map_variables variable_1="k_f" variable_2="k_f"/>
      <map_variables variable_1="Y" variable_2="Y"/>
   </connection>
   <connection>
      <map_components component_1="V_2" component_2="V2"/>
	  <map_variables variable_1="A" variable_2="V_2"/>
	  <map_variables variable_1="B" variable_2="V_M2"/>
	  <map_variables variable_1="C" variable_2="Z"/>
	  <map_variables variable_1="D" variable_2="one_mm"/>
	  <map_variables variable_1="E" variable_2="one_mm"/>
	  <map_variables variable_1="F" variable_2="K_2"/>
	  <map_variables variable_1="G" variable_2="zero"/>
	  <map_variables variable_1="H" variable_2="zero"/>
	  <map_variables variable_1="I" variable_2="two"/>
	  <map_variables variable_1="J" variable_2="one_dm"/>
	  <map_variables variable_1="K" variable_2="one_dm"/>
   </connection>
   <connection>
      <map_components component_1="V_3" component_2="V3"/>
	  <map_variables variable_1="A" variable_2="V_3"/>
	  <map_variables variable_1="B" variable_2="V_M3"/>
	  <map_variables variable_1="C" variable_2="Z"/>
	  <map_variables variable_1="D" variable_2="Y"/>
	  <map_variables variable_1="E" variable_2="A"/>
	  <map_variables variable_1="F" variable_2="K_Z"/>
	  <map_variables variable_1="G" variable_2="K_Y"/>
	  <map_variables variable_1="H" variable_2="K_A"/>
	  <map_variables variable_1="I" variable_2="m"/>
	  <map_variables variable_1="J" variable_2="two"/>
	  <map_variables variable_1="K" variable_2="four"/>
   </connection>
   <connection>
	  <map_components component_1="V_5" component_2="V5"/>
	  <map_variables variable_1="A" variable_2="V_5"/>
	  <map_variables variable_1="B" variable_2="V_M5"/>
	  <map_variables variable_1="C" variable_2="A"/>
	  <map_variables variable_1="D" variable_2="Z"/>
	  <map_variables variable_1="E" variable_2="one_mm"/>
	  <map_variables variable_1="F" variable_2="K_5"/>
	  <map_variables variable_1="G" variable_2="K_d"/>
	  <map_variables variable_1="H" variable_2="zero"/>
	  <map_variables variable_1="I" variable_2="p"/>
	  <map_variables variable_1="J" variable_2="n"/>
	  <map_variables variable_1="K" variable_2="one_dm"/>
   </connection>
   <connection>
		<map_components component_1="internalpool" component_2="internal_pool"/>
		<map_variables variable_1="A" variable_2="Y"/>
		<map_variables variable_1="A_init" variable_2="Y_init"/>
		<map_variables variable_1="time" variable_2="time"/>
		<map_variables variable_1="B" variable_2="one"/>
		<map_variables variable_1="C" variable_2="V_2"/>
		<map_variables variable_1="D" variable_2="V_3"/>
		<map_variables variable_1="E" variable_2="k_f"/>
   </connection>
   <connection>
		<map_components component_1="InsP3conc" component_2="InsP3_conc"/>
		<map_variables variable_1="A" variable_2="A"/>
		<map_variables variable_1="A_init" variable_2="A_init"/>
		<map_variables variable_1="time" variable_2="time"/>
		<map_variables variable_1="B" variable_2="beta"/>
		<map_variables variable_1="C" variable_2="V_4"/>
		<map_variables variable_1="D" variable_2="V_5"/>
		<map_variables variable_1="E" variable_2="epsilon"/>
   </connection>
   
   	<!-- metadata -->
	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		<rdf:Description rdf:about="houart_1999_1.1">
			<dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:parseType="Resource">
				<vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
					<vCard:Family>Zhu</vCard:Family>
					<vCard:Given>Xin Yue</vCard:Given>
					<vCard:Other>Mona</vCard:Other>
				</vCard:N>
				<vCard:EMAIL xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
					<rdf:value>xzhu055@aucklanduni.ac.nz</rdf:value>
					<rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
				</vCard:EMAIL>
				<vCard:ORG xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
					<vCard:Orgname>The University of Auckland</vCard:Orgname>
					<vCard:Orgunit>Auckland Bioengineering Institute</vCard:Orgunit>
				</vCard:ORG>
			</dc:creator>
			<dcterms:created xmlns:dcterms="http://purl.org/dc/terms/" rdf:parseType="Resource">
				<dcterms:W3CDTF>14-12-2010</dcterms:W3CDTF>
			</dcterms:created>
		</rdf:Description>
		<rdf:Description rdf:about="#houart_1999_1.1">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">Bursting, Chaos and birhythmicity Originating from Self-modulation of the Inositol 1,4,5-triphosphate Signal in a Model for Intracellular Ca2+ Oscillations</dc:title>
			<cmeta:comment rdf:parseType="Resource">
				<rdf:value>This is the CellML description of the 1999 mathematical model developed by Gerald Houart, Genevieve Dupont and Albert Goldbeter for the various types of complex Ca2+ oscillations based on the mechanism of CICR.
				</rdf:value>
				<dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:parseType="Resource">
					<vCard:FN xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Mona Zhu</vCard:FN>
				</dc:creator>
			</cmeta:comment>
			<cmeta:species>Ca2+</cmeta:species>
			<bqs:reference xmlns:bqs="http://www.cellml.org/bqs/1.1#" 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>Bursting</rdf:li>
							<rdf:li>Chaos</rdf:li>
							<rdf:li>birhythmicity</rdf:li>
							<rdf:li>Inositol</rdf:li>
							<rdf:li>1,4,5-triphosphate</rdf:li>
							<rdf:li>Ca2+ oscillations</rdf:li>
							<rdf:li>CICR</rdf:li>
							<rdf:li>Calcium Dynamics</rdf:li>
						</rdf:Bag>
					</rdf:value>
				</dc:subject>
			</bqs:reference>
			<bqs:reference xmlns:bqs="http://www.cellml.org/bqs/1.1#" rdf:parseType="Resource">
				<bqs:Pubmed_id>17883229</bqs:Pubmed_id>
				<bqs:JournalArticle rdf:parseType="Resource">
					<dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/">
						<rdf:Seq>
							<rdf:li rdf:parseType="Resource">
								<bqs:Person rdf:parseType="Resource">
									<vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
										<vCard:Family>Houart</vCard:Family>
										<vCard:Given>Gerald</vCard:Given>
									</vCard:N>
								</bqs:Person>
							</rdf:li>
							<rdf:li rdf:parseType="Resource">
								<bqs:Person rdf:parseType="Resource">
									<vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
										<vCard:Family>Dupont</vCard:Family>
										<vCard:Given>Genevieve</vCard:Given>
									</vCard:N>
								</bqs:Person>
							</rdf:li>
							<rdf:li rdf:parseType="Resource">
								<bqs:Person rdf:parseType="Resource">
									<vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
										<vCard:Family>Albert</vCard:Family>
										<vCard:Given>Goldbeter</vCard:Given>
									</vCard:N>
								</bqs:Person>
							</rdf:li>
						</rdf:Seq>
					</dc:creator>
					<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">Bursting, Chaos and birhythmicity Originating from Self-modulation of the Inositol 1,4,5-triphosphate Signal in a Model for Intracellular Ca2+ Oscillations</dc:title>
					<dcterms:issued xmlns:dcterms="http://purl.org/dc/terms/" rdf:parseType="Resource">
					</dcterms:issued>
					<bqs:Journal rdf:parseType="Resource">
						<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">Bulletin of Mathematical Biology</dc:title>
					</bqs:Journal>
					<bqs:volume>61</bqs:volume>
					<bqs:first_page>507</bqs:first_page>
					<bqs:last_page>530</bqs:last_page>
				</bqs:JournalArticle>
			</bqs:reference>
		</rdf:Description>
	</rdf:RDF>
</model>