Rendering of the source text

<?xml version="1.0" encoding="utf-8"?>
<!--
This CellML file was generated on 11/16/2010 at 12:24:55 at PM using:

COR (0.9.31.1409)
Copyright 2002-2010 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 name="houart_1999" cmeta:id="houart_1999" 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>
			<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 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 COR and OpenCell.</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.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>
	<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#">
		<rdf:Description rdf:about="">
			<dc:creator rdf:parseType="Resource">
				<vCard:N rdf:parseType="Resource">
					<vCard:Family>Zhu</vCard:Family>
					<vCard:Given>Xin Yue</vCard:Given>
					<vCard:Other>Mona</vCard:Other>
				</vCard:N>
				<vCard:EMAIL 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 rdf:parseType="Resource">
					<vCard:Orgname>The University of Auckland</vCard:Orgname>
					<vCard:Orgunit>Auckland Bioengineering Institute</vCard:Orgunit>
				</vCard:ORG>
			</dc:creator>
			<dcterms:created rdf:parseType="Resource">
				<dcterms:W3CDTF>16-11-2010</dcterms:W3CDTF>
			</dcterms:created>
		</rdf:Description>
		<rdf:Description rdf:about="#houart_1999">
			<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 xmlns:cmeta="http://www.cellml.org/metadata/1.0#" 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 xmlns:cmeta="http://www.cellml.org/metadata/1.0#">Ca2+</cmeta:species>
			<bqs:reference rdf:parseType="Resource">
				<dc:subject 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 rdf:parseType="Resource">
				<bqs:Pubmed_id>17883229</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>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 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 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>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 rdf:parseType="Resource">
						<dcterms:W3CDTF>1999-05</dcterms:W3CDTF>
					</dcterms:issued>
					<bqs:Journal rdf:parseType="Resource">
						<dc:title>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>   <units xmlns="http://www.cellml.org/cellml/1.0#" name="minute">
      <unit units="second" multiplier="60"/>
   </units>
   <units xmlns="http://www.cellml.org/cellml/1.0#" name="per_minute">
      <unit units="minute" exponent="-1"/>
   </units>
   <units xmlns="http://www.cellml.org/cellml/1.0#" name="per_litre">
      <unit units="litre" exponent="-1"/>
   </units>
   <units xmlns="http://www.cellml.org/cellml/1.0#" name="micromolar">
      <unit units="mole" prefix="micro"/>
      <unit units="per_litre"/>
   </units>
   <units xmlns="http://www.cellml.org/cellml/1.0#" name="micromolar_per_minute">
      <unit units="micromolar"/>
      <unit units="per_minute"/>
   </units>
   <units xmlns="http://www.cellml.org/cellml/1.0#" name="per_micromolar_minute">
      <unit units="micromolar" exponent="-1"/>
      <unit units="per_minute"/>
   </units>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="environment">
      <variable name="time" units="minute" public_interface="out"/>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="Vin">
      <variable name="V_0" units="micromolar_per_minute" initial_value="2" public_interface="out"/>
      <variable name="V_1" units="micromolar_per_minute" initial_value="2" public_interface="out"/>
      <variable name="beta" units="dimensionless" initial_value="0.6" public_interface="out"/>
      <variable name="V_in" units="micromolar_per_minute" public_interface="out"/>
      <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 xmlns="http://www.cellml.org/cellml/1.0#" name="V2">
      <variable name="V_M2" units="micromolar_per_minute" initial_value="6"/>
      <variable name="Z" units="micromolar" public_interface="in" private_interface="out"/>
      <variable name="K_2" units="micromolar" initial_value="0.1"/>
      <variable name="V_2" units="micromolar_per_minute" public_interface="out"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>V_2</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>V_M2</ci>
                  <apply>
                     <power/>
                     <ci>Z</ci>
                     <cn cellml:units="dimensionless">2</cn>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <power/>
                     <ci>K_2</ci>
                     <cn cellml:units="dimensionless">2</cn>
                  </apply>
                  <apply>
                     <power/>
                     <ci>Z</ci>
                     <cn cellml:units="dimensionless">2</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="V3">
      <variable name="V_M3" units="micromolar_per_minute" initial_value="20"/>
      <variable name="K_Z" units="micromolar" initial_value="0.5"/>
      <variable name="K_A" units="micromolar" initial_value="0.2"/>
      <variable name="K_Y" units="micromolar" initial_value="0.2"/>
      <variable name="m" units="dimensionless" initial_value="2"/>
      <variable name="Y" units="micromolar" public_interface="in" private_interface="out"/>
      <variable name="A" units="micromolar" public_interface="in" private_interface="out"/>
      <variable name="Z" units="micromolar" public_interface="in" private_interface="out"/>
      <variable name="V_3" units="micromolar_per_minute" public_interface="out"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>V_3</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <ci>V_M3</ci>
                              <apply>
                                 <power/>
                                 <ci>Z</ci>
                                 <ci>m</ci>
                              </apply>
                           </apply>
                           <apply>
                              <plus/>
                              <apply>
                                 <power/>
                                 <ci>K_Z</ci>
                                 <ci>m</ci>
                              </apply>
                              <apply>
                                 <power/>
                                 <ci>Z</ci>
                                 <ci>m</ci>
                              </apply>
                           </apply>
                        </apply>
                        <apply>
                           <power/>
                           <ci>Y</ci>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                     <apply>
                        <plus/>
                        <apply>
                           <power/>
                           <ci>K_Y</ci>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                        <apply>
                           <power/>
                           <ci>Y</ci>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <power/>
                     <ci>A</ci>
                     <cn cellml:units="dimensionless">4</cn>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <power/>
                     <ci>K_A</ci>
                     <cn cellml:units="dimensionless">4</cn>
                  </apply>
                  <apply>
                     <power/>
                     <ci>A</ci>
                     <cn cellml:units="dimensionless">4</cn>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="V5">
      <variable name="V_M5" units="micromolar_per_minute" initial_value="5"/>
      <variable name="K_5" units="micromolar" initial_value="1"/>
      <variable name="K_d" units="micromolar" initial_value="0.4"/>
      <variable name="p" units="dimensionless" initial_value="2"/>
      <variable name="n" units="dimensionless" initial_value="4"/>
      <variable name="A" units="micromolar" public_interface="in" private_interface="out"/>
      <variable name="Z" units="micromolar" public_interface="in" private_interface="out"/>
      <variable name="V_5" units="micromolar_per_minute" public_interface="out"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>V_5</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>V_M5</ci>
                        <apply>
                           <power/>
                           <ci>A</ci>
                           <ci>p</ci>
                        </apply>
                     </apply>
                     <apply>
                        <plus/>
                        <apply>
                           <power/>
                           <ci>K_5</ci>
                           <ci>p</ci>
                        </apply>
                        <apply>
                           <power/>
                           <ci>A</ci>
                           <ci>p</ci>
                        </apply>
                     </apply>
                  </apply>
                  <apply>
                     <power/>
                     <ci>Z</ci>
                     <ci>n</ci>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <power/>
                     <ci>K_d</ci>
                     <ci>n</ci>
                  </apply>
                  <apply>
                     <power/>
                     <ci>Z</ci>
                     <ci>n</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="cytosol">
      <variable name="time" units="minute" public_interface="in" private_interface="out"/>
      <variable name="V_in" units="micromolar_per_minute" public_interface="in"/>
      <variable name="V_2" units="micromolar_per_minute" public_interface="in"/>
      <variable name="V_3" units="micromolar_per_minute" public_interface="in" private_interface="out"/>
      <variable name="k" units="per_minute" initial_value="10"/>
      <variable name="k_f" units="per_minute" initial_value="1" public_interface="out"/>
      <variable name="Z" units="micromolar" initial_value="0.15" public_interface="out"/>
      <variable name="Y" units="micromolar" public_interface="in" private_interface="out"/>
      <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 xmlns="http://www.cellml.org/cellml/1.0#" name="internal_pool">
      <variable name="time" units="minute" public_interface="in" private_interface="out"/>
      <variable name="V_2" units="micromolar_per_minute" public_interface="in" private_interface="out"/>
      <variable name="V_3" units="micromolar_per_minute" public_interface="in" private_interface="out"/>
      <variable name="k_f" units="per_minute" public_interface="in" private_interface="out"/>
      <variable name="Y" units="micromolar" initial_value="1" public_interface="out"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Y</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <minus/>
                  <ci>V_2</ci>
                  <ci>V_3</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>k_f</ci>
                  <ci>Y</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="InsP3_conc">
      <variable name="time" units="minute" public_interface="in"/>
      <variable name="beta" units="dimensionless" public_interface="in"/>
      <variable name="epsilon" units="per_minute" initial_value="0.1"/>
      <variable name="V_4" units="micromolar_per_minute" initial_value="2"/>
      <variable name="V_5" units="micromolar_per_minute" public_interface="in" private_interface="out"/>
      <variable name="A" units="micromolar" initial_value="0.42" public_interface="out"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>A</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>beta</ci>
                     <ci>V_4</ci>
                  </apply>
                  <ci>V_5</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>epsilon</ci>
                  <ci>A</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="Vin" component_2="cytosol"/>
      <map_variables variable_1="V_in" variable_2="V_in"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="Vin" component_2="InsP3_conc"/>
      <map_variables variable_1="beta" variable_2="beta"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <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 xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="V2" component_2="internal_pool"/>
      <map_variables variable_1="V_2" variable_2="V_2"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <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 xmlns="http://www.cellml.org/cellml/1.0#">
      <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 xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="V3" component_2="InsP3_conc"/>
      <map_variables variable_1="A" variable_2="A"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <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 xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="environment" component_2="cytosol"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="environment" component_2="internal_pool"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="environment" component_2="InsP3_conc"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="cytosol" component_2="V5"/>
      <map_variables variable_1="Z" variable_2="Z"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <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>
</model>