Rendering of the source text

<?xml version="1.0" encoding="utf-8"?>
<!--
This CellML file was generated on 14/05/2010 at 11:00:54 at a.m. using:

COR (0.9.31.1371)
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 xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" name="snyder_2000" cmeta:id="snyder_2000" 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>A Mathematical Model Of Cardiocyte Ca2+ Dynamics</title>
				<author>
					<firstname>Catherine</firstname>
					<surname>Lloyd</surname>
					<affiliation>
						<shortaffil>Auckland Bioengineering Institute, The University of Auckland</shortaffil>
					</affiliation>
				</author>
			</articleinfo>
			<section id="sec_status">
				<title>Model Status</title>
				<para>
            This model runs in OpenCell and COR and the units are consistent throughout. It does not at present replicate the published results, most likely because there are a number of gating or switching processes described in the paper that were unable to be replicated precisely.  
          </para>
			</section>
			<sect1 id="sec_structure">
				<title>Model Structure</title>
				<para>
ABSTRACT: Cardiac contraction and relaxation dynamics result from a set of simultaneously interacting Ca(2+) regulatory mechanisms. In this study, cardiocyte Ca(2+) dynamics were modeled using a set of six differential equations that were based on theories, equations, and parameters described in previous studies. Among the unique features of the model was the inclusion of bidirectional modulatory interplay between the sarcoplasmic reticular Ca(2+) release channel (SRRC) and calsequestrin (CSQ) in the SR lumen, where CSQ acted as a dynamic rather than simple Ca(2+) buffer, and acted as a Ca(2+) sensor in the SR lumen as well. The inclusion of this control mechanism was central in overcoming a number of assumptions that would otherwise have to be made about SRRC kinetics, SR Ca(2+) release rates, and SR Ca(2+) release termination when the SR lumen is assumed to act as a simple, buffered Ca(2+) sink. The model was sufficient to reproduce a graded Ca(2+)-induced Ca(2+) release (CICR) response, CICR with high gain, and a system with reasonable stability. As constructed, the model successfully replicated the results of several previously published experiments that dealt with the Ca(2+) dependence of the SRRC (, J. Gen. Physiol. 85:247-289), the refractoriness of the SRRC (, Am. J. Physiol. 270:C148-C159), the SR Ca(2+) load dependence of SR Ca(2+) release (, Am. J. Physiol. 268:C1313-C1329;, J. Biol. Chem. 267:20850-20856), SR Ca(2+) leak (, J. Physiol. (Lond.). 474:463-471;, Biophys. J. 68:2015-2022), SR Ca(2+) load regulation by leak and uptake (, J. Gen. Physiol. 111:491-504), the effect of Ca(2+) trigger duration on SR Ca(2+) release (, Am. J. Physiol. 258:C944-C954), the apparent relationship that exists between sarcoplasmic and sarcoplasmic reticular calcium concentrations (, Biophys. J. 73:1524-1531), and a variety of contraction frequency-dependent alterations in sarcoplasmic [Ca(2+)] dynamics that are normally observed in the laboratory, including rest potentiation, a negative frequency-[Ca(2+)] relationship, and extrasystolic potentiation. Furthermore, under the condition of a simulated Ca(2+) overload, an alternans-like state was produced. In summary, the current model of cardiocyte Ca(2+) dynamics provides an integrated theoretical framework of fundamental cellular Ca(2+) regulatory processes that is sufficient to predict a broad array of observable experimental outcomes.
</para>
				<para>
The original paper reference is cited below:
</para>
				<para>
A Mathematical Model of Cardiocyte Ca<superscript>2+</superscript> Dynamics with a Novel Representation of Sarcoplasmic Reticular Ca<superscript>2+</superscript> Control, Steven M. Snyder, Bradley M. Palmer and Russell L. Moore, 2000, <emphasis>Biophysical Journal</emphasis>, 79, 94-115. <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&amp;db=PubMed&amp;list_uids=10866940&amp;dopt=Abstract">PubMed ID: 10866940</ulink>
				</para>
				<informalfigure float="0" id="fig_cell_diagram">
					<mediaobject>
						<imageobject>
							<objectinfo>
								<title>cell schematic for the model</title>
							</objectinfo>
							<imagedata fileref="snyder_2000.png"/>
						</imageobject>
					</mediaobject>
					<caption>Cell schematic for the model representing three intracellular compartments with buffering, the extracellular space, and the related Ca<superscript>2+</superscript> movement processes.</caption>
				</informalfigure>
			</sect1>
		</article>
	</documentation>
	<rdf:RDF>
		<rdf:Bag rdf:about="rdf:#e60dc7a6-90ff-416c-bbd7-913b0604dd8f">
			<rdf:li>calcium dynamics</rdf:li>
			<rdf:li>cardiac myocyte</rdf:li>
			<rdf:li>electrophysiology</rdf:li>
			<rdf:li>Cardiac Myocyte</rdf:li>
		</rdf:Bag>
		<rdf:Seq rdf:about="rdf:#a2464351-ed51-4f29-baf3-bf080cb2e1ba">
			<rdf:li rdf:resource="rdf:#e1cf66d3-5d67-45b6-9abd-26da004e05b6"/>
			<rdf:li rdf:resource="rdf:#3a21cb28-a446-4419-ae7a-3af777a3b642"/>
			<rdf:li rdf:resource="rdf:#34925f2b-ebf6-4a87-8e53-250a73226169"/>
		</rdf:Seq>
		<rdf:Description rdf:about="rdf:#878b537b-39ab-41bb-b57b-889d8c3be05c">
			<bqs:Pubmed>10866940</bqs:Pubmed>
			<bqs:JournalArticle rdf:resource="rdf:#28f3b3b2-8907-4f09-9ffe-63e2f5d65f2e"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#27698139-9a52-4d5f-8077-f464d06dd16d">
			<vCard:Given>Bradley</vCard:Given>
			<vCard:Family>Palmer</vCard:Family>
			<vCard:Other>M</vCard:Other>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#7c921c2f-bbb8-4156-9f27-208642f31667">
			<bqs:subject_type>keyword</bqs:subject_type>
			<rdf:value rdf:resource="rdf:#e60dc7a6-90ff-416c-bbd7-913b0604dd8f"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#d2a4ad33-5d8a-463b-b832-a7c7a10d977b">
			<dcterms:modified rdf:resource="rdf:#2bd16e58-830f-4e3e-8699-ea606e5b7455"/>
			<rdf:value>
          Added more metadata.
        </rdf:value>
			<cmeta:modifier xmlns:cmeta="http://www.cellml.org/metadata/1.0#" rdf:resource="rdf:#e16645ca-c614-4c2f-a138-6f48fbaa5b29"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#509faf33-4a7f-45fd-9f48-91e1b9db45be">
			<vCard:Given>Peter</vCard:Given>
			<vCard:Family>Villiger</vCard:Family>
			<vCard:Other>J</vCard:Other>
		</rdf:Description>
		<rdf:Description rdf:about="#snyder_2000">
			<dc:title>
        Snyder et al's 2000 mathematical model of cardiocyte Ca2+ dynamics
      </dc:title>
			<cmeta:bio_entity xmlns:cmeta="http://www.cellml.org/metadata/1.0#">Cardiac Myocyte</cmeta:bio_entity>
			<cmeta:comment xmlns:cmeta="http://www.cellml.org/metadata/1.0#" rdf:resource="rdf:#85e890e1-6713-4683-8dbe-3d4b8147d047"/>
			<bqs:reference rdf:resource="rdf:#878b537b-39ab-41bb-b57b-889d8c3be05c"/>
			<bqs:reference rdf:resource="rdf:#d5ac8eb2-c52f-480e-b220-c4f0cf97ff2d"/>
			<cmeta:species xmlns:cmeta="http://www.cellml.org/metadata/1.0#">Rat</cmeta:species>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#de896ba1-eb95-4753-a778-9b378bad3c79">
			<rdf:type rdf:resource="http://www.w3.org/2001/vcard-rdf/3.0#FN"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#2bd16e58-830f-4e3e-8699-ea606e5b7455">
			<dcterms:W3CDTF>2002-07-22</dcterms:W3CDTF>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#e16645ca-c614-4c2f-a138-6f48fbaa5b29">
			<vCard:N rdf:resource="rdf:#bd129a4b-6070-403b-9e3a-2bb1f9e86b3a"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#d5ac8eb2-c52f-480e-b220-c4f0cf97ff2d">
			<dc:subject rdf:resource="rdf:#7c921c2f-bbb8-4156-9f27-208642f31667"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#a6350034-4e4d-49eb-8231-4b4d8e643de3">
			<vCard:N rdf:resource="rdf:#509faf33-4a7f-45fd-9f48-91e1b9db45be"/>
		</rdf:Description>
		<rdf:Description rdf:about="">
			<dc:publisher>
        The University of Auckland, Auckland Bioengineering Institute
      </dc:publisher>
			<cmeta:modification xmlns:cmeta="http://www.cellml.org/metadata/1.0#" rdf:resource="rdf:#d2a4ad33-5d8a-463b-b832-a7c7a10d977b"/>
			<cmeta:modification xmlns:cmeta="http://www.cellml.org/metadata/1.0#" rdf:resource="rdf:#db2cc2b2-7166-4676-b335-f5451de91efd"/>
			<dcterms:created rdf:resource="rdf:#89e71652-c578-4a50-813d-200c81a8745e"/>
			<dc:creator rdf:resource="rdf:#80dbdf64-8014-443a-a90f-9fc9f2a871b2"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#db2cc2b2-7166-4676-b335-f5451de91efd">
			<dcterms:modified rdf:resource="rdf:#b0c42bbd-6d35-472f-8534-a29a64590121"/>
			<rdf:value>
          Removed unnecessary type. Made MathML id's unique
        </rdf:value>
			<cmeta:modifier xmlns:cmeta="http://www.cellml.org/metadata/1.0#" rdf:resource="rdf:#a6350034-4e4d-49eb-8231-4b4d8e643de3"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#4be99089-cfab-4ee2-a221-484688069f4c">
			<vCard:Orgname>The University of Auckland</vCard:Orgname>
			<vCard:Orgunit>Auckland Bioengineering Institute</vCard:Orgunit>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#8b0ea24f-97a7-420e-806a-8ca68d97b4bb">
			<vCard:Given>Russell</vCard:Given>
			<vCard:Family>Moore</vCard:Family>
			<vCard:Other>L</vCard:Other>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#bd129a4b-6070-403b-9e3a-2bb1f9e86b3a">
			<vCard:Given>Catherine</vCard:Given>
			<vCard:Family>Lloyd</vCard:Family>
			<vCard:Other>May</vCard:Other>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#b0c42bbd-6d35-472f-8534-a29a64590121">
			<dcterms:W3CDTF>2005-04-22</dcterms:W3CDTF>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#85e890e1-6713-4683-8dbe-3d4b8147d047">
			<dc:creator rdf:resource="rdf:#de896ba1-eb95-4753-a778-9b378bad3c79"/>
			<rdf:value>
          Below is a CellML description of snyder et al's 2000 mathematical model of rat cardiocyte Ca2+ dynamics with a novel representation of sarcoplasmic reticular Ca2+ control.  The mathematical representations of the movement mechanisms and the influence of Ca2+ buffering are based on experimental data taken from rat cardiac myocytes.
        </rdf:value>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#9006e114-a5df-4b72-8b5d-dc92c1c0a633">
			<dc:title>Biophysical Journal</dc:title>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#5ee1dc1f-59d1-4c38-a6f0-92dd79cfb688">
			<dcterms:W3CDTF>2000</dcterms:W3CDTF>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#34925f2b-ebf6-4a87-8e53-250a73226169">
			<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
			<vCard:N rdf:resource="rdf:#8b0ea24f-97a7-420e-806a-8ca68d97b4bb"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#e1cf66d3-5d67-45b6-9abd-26da004e05b6">
			<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
			<vCard:N rdf:resource="rdf:#2088dc91-3697-4891-acca-ac2e91dc015f"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#2088dc91-3697-4891-acca-ac2e91dc015f">
			<vCard:Given>Steven</vCard:Given>
			<vCard:Family>Snyder</vCard:Family>
			<vCard:Other>M</vCard:Other>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#3a21cb28-a446-4419-ae7a-3af777a3b642">
			<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
			<vCard:N rdf:resource="rdf:#27698139-9a52-4d5f-8077-f464d06dd16d"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#80dbdf64-8014-443a-a90f-9fc9f2a871b2">
			<vCard:ORG rdf:resource="rdf:#4be99089-cfab-4ee2-a221-484688069f4c"/>
			<vCard:EMAIL rdf:resource="rdf:#b23c8392-867a-4577-8142-94e74945a9de"/>
			<vCard:N rdf:resource="rdf:#1029f99a-ac57-4c99-a3bf-7c9e2085c073"/>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#1029f99a-ac57-4c99-a3bf-7c9e2085c073">
			<vCard:Given>Catherine</vCard:Given>
			<vCard:Family>Lloyd</vCard:Family>
			<vCard:Other>May</vCard:Other>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#28f3b3b2-8907-4f09-9ffe-63e2f5d65f2e">
			<dc:creator rdf:resource="rdf:#a2464351-ed51-4f29-baf3-bf080cb2e1ba"/>
			<dc:title>
            A Mathematical Model of Cardiocyte Ca2+ Dynamics with a Novel Representation of Sarcoplasmic Reticular Ca2+ Control
          </dc:title>
			<bqs:volume>79</bqs:volume>
			<bqs:first_page>94</bqs:first_page>
			<bqs:Journal rdf:resource="rdf:#9006e114-a5df-4b72-8b5d-dc92c1c0a633"/>
			<dcterms:issued rdf:resource="rdf:#5ee1dc1f-59d1-4c38-a6f0-92dd79cfb688"/>
			<bqs:last_page>115</bqs:last_page>
		</rdf:Description>
		<rdf:Description rdf:about="rdf:#b23c8392-867a-4577-8142-94e74945a9de">
			<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:#89e71652-c578-4a50-813d-200c81a8745e">
			<dcterms:W3CDTF>2002-03-18</dcterms:W3CDTF>
		</rdf:Description>
	</rdf:RDF>   
	
   <units xmlns="http://www.cellml.org/cellml/1.0#" name="molar">
      <unit units="mole"/>
      <unit units="litre" exponent="-1"/>
   </units>
   <units xmlns="http://www.cellml.org/cellml/1.0#" name="flux">
      <unit units="molar"/>
      <unit units="second" exponent="-1"/>
   </units>
   <units xmlns="http://www.cellml.org/cellml/1.0#" name="first_order_rate_constant">
      <unit units="second" exponent="-1"/>
   </units>
   <units xmlns="http://www.cellml.org/cellml/1.0#" name="second_order_rate_constant">
      <unit units="molar" exponent="-1"/>
      <unit units="second" exponent="-1"/>
   </units>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="environment">
      <variable name="time" units="second" public_interface="out"/>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="Ca_influx_across_the_SR">
      <variable name="Rate_Ca_influx_across_the_SR" units="flux" public_interface="out"/>
      <variable name="Ca_e" units="molar" public_interface="in"/>
      <variable name="Ca_f" units="molar" public_interface="in"/>
      <variable name="k1" units="first_order_rate_constant" initial_value="0.2"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>Rate_Ca_influx_across_the_SR</ci>
            <apply>
               <times/>
               <ci>k1</ci>
               <apply>
                  <minus/>
                  <ci>Ca_e</ci>
                  <ci>Ca_f</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="Ca_movement_through_the_SRRC">
      <variable name="Ca_2_S1" units="molar" public_interface="in"/>
      <variable name="Ca_S2" units="molar" public_interface="in"/>
      <variable name="S1" units="molar" public_interface="in"/>
      <variable name="S2" units="molar" public_interface="in"/>
      <variable name="Ca_f" units="molar" public_interface="in"/>
      <variable name="Ca_s" units="molar" public_interface="in"/>
      <variable name="dCa2_S1_dt" units="flux" public_interface="out"/>
      <variable name="dCa_S2_dt" units="flux" public_interface="out"/>
      <variable name="k_on1" units="second_order_rate_constant" initial_value="2000000000"/>
      <variable name="k_off1" units="first_order_rate_constant" initial_value="1400"/>
      <variable name="k_on2" units="second_order_rate_constant" initial_value="13000000"/>
      <variable name="k_off2" units="first_order_rate_constant" initial_value="3.9"/>
      <variable name="k_s" units="first_order_rate_constant" initial_value="9"/>
      <variable name="cas1" units="dimensionless" initial_value="0"/>
      <variable name="cas2" units="dimensionless" initial_value="0"/>
      <variable name="dcas1_dt" units="first_order_rate_constant" public_interface="out"/>
      <variable name="dcas2_dt" units="first_order_rate_constant" public_interface="out"/>
      <variable name="r_o" units="dimensionless"/>
      <variable name="Rate_Ca_movement_through_the_SRRC" units="flux" public_interface="out"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>dCa2_S1_dt</ci>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>k_on1</ci>
                  <ci>Ca_f</ci>
                  <ci>S1</ci>
               </apply>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <apply>
                        <power/>
                        <ci>k_off1</ci>
                        <cn cellml:units="dimensionless">2</cn>
                     </apply>
                     <apply>
                        <times/>
                        <ci>k_on1</ci>
                        <ci>Ca_f</ci>
                     </apply>
                  </apply>
                  <ci>Ca_2_S1</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>dCa_S2_dt</ci>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>k_on2</ci>
                  <ci>Ca_f</ci>
                  <ci>S2</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>k_off2</ci>
                  <ci>Ca_S2</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>r_o</ci>
            <apply>
               <times/>
               <ci>cas1</ci>
               <apply>
                  <minus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <ci>cas2</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>cas1</ci>
            </apply>
            <ci>dcas1_dt</ci>
         </apply>
         <apply>
            <eq/>
            <ci>dcas1_dt</ci>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>k_on1</ci>
                  <ci>Ca_f</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>cas1</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <apply>
                        <power/>
                        <ci>k_off1</ci>
                        <cn cellml:units="dimensionless">2</cn>
                     </apply>
                     <apply>
                        <times/>
                        <ci>k_on1</ci>
                        <ci>Ca_f</ci>
                     </apply>
                  </apply>
                  <ci>cas1</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>cas2</ci>
            </apply>
            <ci>dcas2_dt</ci>
         </apply>
         <apply>
            <eq/>
            <ci>dcas2_dt</ci>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>k_on2</ci>
                  <ci>Ca_f</ci>
                  <apply>
                     <minus/>
                     <cn cellml:units="dimensionless">1</cn>
                     <ci>cas2</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>k_off2</ci>
                  <ci>cas2</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>Rate_Ca_movement_through_the_SRRC</ci>
            <apply>
               <times/>
               <ci>k_s</ci>
               <ci>r_o</ci>
               <apply>
                  <minus/>
                  <ci>Ca_s</ci>
                  <ci>Ca_f</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="Ca_efflux_across_the_SR_by_NaCa_exchange">
      <variable name="Ca_f" units="molar" public_interface="in"/>
      <variable name="Km_NaCaX" units="molar" initial_value="0.000036"/>
      <variable name="Vmax_NaCaX" units="flux" initial_value="0.0012"/>
      <variable name="Rate_Ca_efflux_across_the_SR_by_NaCa_exchange" units="flux" public_interface="out"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>Rate_Ca_efflux_across_the_SR_by_NaCa_exchange</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>Vmax_NaCaX</ci>
                  <ci>Ca_f</ci>
               </apply>
               <apply>
                  <plus/>
                  <ci>Km_NaCaX</ci>
                  <ci>Ca_f</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="Ca_movement_between_the_fuzzy_space_and_cytosol">
      <variable name="Ca_c" units="molar" public_interface="in"/>
      <variable name="Ca_f" units="molar" public_interface="in"/>
      <variable name="kf" units="first_order_rate_constant" initial_value="2500"/>
      <variable name="Rate_Ca_movement_between_the_fuzzy_space_and_cytosol" units="flux" public_interface="out"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>Rate_Ca_movement_between_the_fuzzy_space_and_cytosol</ci>
            <apply>
               <times/>
               <ci>kf</ci>
               <apply>
                  <minus/>
                  <ci>Ca_f</ci>
                  <ci>Ca_c</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="Ca_uptake_by_SR_Ca_ATPase">
      <variable name="Ca_c" units="molar" public_interface="in"/>
      <variable name="Ca_s" units="molar" public_interface="in"/>
      <variable name="Km_s" units="molar" initial_value="0.00000025"/>
      <variable name="Vmax_s" units="flux" initial_value="0.000525"/>
      <variable name="Rate_Ca_uptake_by_SR_Ca_ATPase" units="flux" public_interface="out"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>Rate_Ca_uptake_by_SR_Ca_ATPase</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <ci>Vmax_s</ci>
                  <apply>
                     <minus/>
                     <apply>
                        <power/>
                        <ci>Ca_c</ci>
                        <cn cellml:units="dimensionless">2</cn>
                     </apply>
                     <apply>
                        <divide/>
                        <apply>
                           <power/>
                           <ci>Ca_s</ci>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                        <apply>
                           <power/>
                           <cn cellml:units="dimensionless">7000</cn>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <power/>
                     <ci>Km_s</ci>
                     <cn cellml:units="dimensionless">2</cn>
                  </apply>
                  <apply>
                     <power/>
                     <ci>Ca_c</ci>
                     <cn cellml:units="dimensionless">2</cn>
                  </apply>
                  <apply>
                     <divide/>
                     <apply>
                        <power/>
                        <ci>Ca_s</ci>
                        <cn cellml:units="dimensionless">2</cn>
                     </apply>
                     <apply>
                        <power/>
                        <cn cellml:units="dimensionless">7000</cn>
                        <cn cellml:units="dimensionless">2</cn>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="Ca_buffering_in_the_SR">
      <variable name="Ca_CSQ" units="molar" public_interface="in"/>
      <variable name="Ca_s" units="molar" public_interface="in"/>
      <variable name="K_ons" units="second_order_rate_constant" initial_value="8772"/>
      <variable name="K_offs" units="first_order_rate_constant" initial_value="5.596536"/>
      <variable name="Bmax_s" units="molar" initial_value="0.008"/>
      <variable name="Rate_Ca_buffering_in_the_SR" units="flux" public_interface="out"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>Rate_Ca_buffering_in_the_SR</ci>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>K_ons</ci>
                  <ci>Ca_s</ci>
                  <apply>
                     <minus/>
                     <ci>Bmax_s</ci>
                     <ci>Ca_CSQ</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>K_offs</ci>
                  <ci>Ca_CSQ</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="fuzzy_space_calcium">
      <variable name="Ca_f" units="molar" initial_value="0.12e-6" public_interface="out"/>
      <variable name="Rt" units="molar" initial_value="0.00000015"/>
      <variable name="Bmax_f1" units="molar" initial_value="0.0002"/>
      <variable name="Bmax_f2" units="molar" initial_value="0.0011"/>
      <variable name="Kb_f1" units="molar" initial_value="0.000017"/>
      <variable name="Kb_f2" units="molar" initial_value="0.000013"/>
      <variable name="V_f" units="dimensionless" initial_value="0.0013"/>
      <variable name="Rate_Ca_influx_across_the_SR" units="flux" public_interface="in"/>
      <variable name="Rate_Ca_movement_through_the_SRRC" units="flux" public_interface="in"/>
      <variable name="Rate_Ca_efflux_across_the_SR_by_NaCa_exchange" units="flux" public_interface="in"/>
      <variable name="Rate_Ca_movement_between_the_fuzzy_space_and_cytosol" units="flux" public_interface="in"/>
      <variable name="dcas1_dt" units="first_order_rate_constant" public_interface="in"/>
      <variable name="dcas2_dt" units="first_order_rate_constant" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_f</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <plus/>
                  <apply>
                     <minus/>
                     <ci>Rate_Ca_movement_through_the_SRRC</ci>
                     <apply>
                        <plus/>
                        <apply>
                           <times/>
                           <ci>Rt</ci>
                           <apply>
                              <plus/>
                              <ci>dcas1_dt</ci>
                              <ci>dcas2_dt</ci>
                           </apply>
                        </apply>
                        <ci>Rate_Ca_movement_between_the_fuzzy_space_and_cytosol</ci>
                        <ci>Rate_Ca_efflux_across_the_SR_by_NaCa_exchange</ci>
                     </apply>
                  </apply>
                  <ci>Rate_Ca_influx_across_the_SR</ci>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>Bmax_f1</ci>
                        <ci>Kb_f1</ci>
                     </apply>
                     <apply>
                        <power/>
                        <apply>
                           <plus/>
                           <ci>Ca_f</ci>
                           <ci>Kb_f1</ci>
                        </apply>
                        <cn cellml:units="dimensionless">2</cn>
                     </apply>
                  </apply>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>Bmax_f2</ci>
                        <ci>Kb_f2</ci>
                     </apply>
                     <apply>
                        <power/>
                        <apply>
                           <plus/>
                           <ci>Ca_f</ci>
                           <ci>Kb_f2</ci>
                        </apply>
                        <cn cellml:units="dimensionless">2</cn>
                     </apply>
                  </apply>
                  <ci>V_f</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="cytosolic_calcium">
      <variable name="Ca_c" units="molar" initial_value="1e-7" public_interface="out"/>
      <variable name="Bmax_c" units="molar" initial_value="0.00012"/>
      <variable name="dye_c" units="molar" initial_value="0"/>
      <variable name="Kb_c" units="molar" initial_value="0.00000096"/>
      <variable name="Kb_dye" units="molar" initial_value="2e-7"/>
      <variable name="V_c" units="dimensionless" initial_value="0.9287"/>
      <variable name="Rate_Ca_movement_between_the_fuzzy_space_and_cytosol" units="flux" public_interface="in"/>
      <variable name="Rate_Ca_uptake_by_SR_Ca_ATPase" units="flux" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_c</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>Rate_Ca_movement_between_the_fuzzy_space_and_cytosol</ci>
                  <ci>Rate_Ca_uptake_by_SR_Ca_ATPase</ci>
               </apply>
               <apply>
                  <plus/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>Bmax_c</ci>
                        <ci>Kb_c</ci>
                     </apply>
                     <apply>
                        <power/>
                        <apply>
                           <plus/>
                           <ci>Ca_c</ci>
                           <ci>Kb_c</ci>
                        </apply>
                        <cn cellml:units="dimensionless">2</cn>
                     </apply>
                  </apply>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>dye_c</ci>
                        <ci>Kb_dye</ci>
                     </apply>
                     <apply>
                        <power/>
                        <apply>
                           <plus/>
                           <ci>Ca_c</ci>
                           <ci>Kb_dye</ci>
                        </apply>
                        <cn cellml:units="dimensionless">2</cn>
                     </apply>
                  </apply>
                  <ci>V_c</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="SR_calcium">
      <variable name="Ca_s" units="molar" initial_value="201e-6" public_interface="out"/>
      <variable name="V_s" units="dimensionless" initial_value="0.07"/>
      <variable name="Rate_Ca_movement_through_the_SRRC" units="flux" public_interface="in"/>
      <variable name="Rate_Ca_uptake_by_SR_Ca_ATPase" units="flux" public_interface="in"/>
      <variable name="Rate_Ca_buffering_in_the_SR" units="flux" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_s</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <divide/>
                  <apply>
                     <minus/>
                     <ci>Rate_Ca_uptake_by_SR_Ca_ATPase</ci>
                     <ci>Rate_Ca_movement_through_the_SRRC</ci>
                  </apply>
                  <ci>V_s</ci>
               </apply>
               <ci>Rate_Ca_buffering_in_the_SR</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="extracellular_calcium">
      <variable name="Ca_e" units="molar" initial_value="0.002" public_interface="out"/>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="calsequestrin_bound_calcium">
      <variable name="Ca_CSQ" units="molar" initial_value="0" public_interface="out"/>
      <variable name="time" units="second" public_interface="in"/>
      <variable name="Rate_Ca_buffering_in_the_SR" units="flux" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_CSQ</ci>
            </apply>
            <ci>Rate_Ca_buffering_in_the_SR</ci>
         </apply>
      </math>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="SRRC_fast_activating_binding_site">
      <variable name="S1" units="molar" initial_value="0" public_interface="out"/>
      <variable name="dCa2_S1_dt" units="flux" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>S1</ci>
            </apply>
            <apply>
               <minus/>
               <ci>dCa2_S1_dt</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="SRRC_slow_inactivating_binding_site">
      <variable name="S2" units="molar" initial_value="0" public_interface="out"/>
      <variable name="dCa_S2_dt" units="flux" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>S2</ci>
            </apply>
            <apply>
               <minus/>
               <ci>dCa_S2_dt</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="Ca_bound_to_the_SRRC_fast_activating_binding_site">
      <variable name="Ca_2_S1" units="molar" initial_value="0" public_interface="out"/>
      <variable name="dCa2_S1_dt" units="flux" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_2_S1</ci>
            </apply>
            <ci>dCa2_S1_dt</ci>
         </apply>
      </math>
   </component>
   <component xmlns="http://www.cellml.org/cellml/1.0#" name="Ca_bound_to_the_SRRC_slow_inactivating_binding_site">
      <variable name="Ca_S2" units="molar" initial_value="0" public_interface="out"/>
      <variable name="dCa_S2_dt" units="flux" public_interface="in"/>
      <variable name="time" units="second" public_interface="in"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>Ca_S2</ci>
            </apply>
            <ci>dCa_S2_dt</ci>
         </apply>
      </math>
   </component>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="Ca_influx_across_the_SR" component_2="extracellular_calcium"/>
      <map_variables variable_1="Ca_e" variable_2="Ca_e"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="Ca_influx_across_the_SR" component_2="fuzzy_space_calcium"/>
      <map_variables variable_1="Ca_f" variable_2="Ca_f"/>
      <map_variables variable_1="Rate_Ca_influx_across_the_SR" variable_2="Rate_Ca_influx_across_the_SR"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="Ca_movement_through_the_SRRC" component_2="Ca_bound_to_the_SRRC_fast_activating_binding_site"/>
      <map_variables variable_1="dCa2_S1_dt" variable_2="dCa2_S1_dt"/>
      <map_variables variable_1="Ca_2_S1" variable_2="Ca_2_S1"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="Ca_movement_through_the_SRRC" component_2="Ca_bound_to_the_SRRC_slow_inactivating_binding_site"/>
      <map_variables variable_1="dCa_S2_dt" variable_2="dCa_S2_dt"/>
      <map_variables variable_1="Ca_S2" variable_2="Ca_S2"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="Ca_movement_through_the_SRRC" component_2="SRRC_fast_activating_binding_site"/>
      <map_variables variable_1="dCa2_S1_dt" variable_2="dCa2_S1_dt"/>
      <map_variables variable_1="S1" variable_2="S1"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="Ca_movement_through_the_SRRC" component_2="SRRC_slow_inactivating_binding_site"/>
      <map_variables variable_1="dCa_S2_dt" variable_2="dCa_S2_dt"/>
      <map_variables variable_1="S2" variable_2="S2"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="Ca_movement_through_the_SRRC" component_2="fuzzy_space_calcium"/>
      <map_variables variable_1="Ca_f" variable_2="Ca_f"/>
      <map_variables variable_1="dcas1_dt" variable_2="dcas1_dt"/>
      <map_variables variable_1="dcas2_dt" variable_2="dcas2_dt"/>
      <map_variables variable_1="Rate_Ca_movement_through_the_SRRC" variable_2="Rate_Ca_movement_through_the_SRRC"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="Ca_movement_through_the_SRRC" component_2="SR_calcium"/>
      <map_variables variable_1="Ca_s" variable_2="Ca_s"/>
      <map_variables variable_1="Rate_Ca_movement_through_the_SRRC" variable_2="Rate_Ca_movement_through_the_SRRC"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="Ca_efflux_across_the_SR_by_NaCa_exchange" component_2="fuzzy_space_calcium"/>
      <map_variables variable_1="Ca_f" variable_2="Ca_f"/>
      <map_variables variable_1="Rate_Ca_efflux_across_the_SR_by_NaCa_exchange" variable_2="Rate_Ca_efflux_across_the_SR_by_NaCa_exchange"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="Ca_movement_between_the_fuzzy_space_and_cytosol" component_2="fuzzy_space_calcium"/>
      <map_variables variable_1="Ca_f" variable_2="Ca_f"/>
      <map_variables variable_1="Rate_Ca_movement_between_the_fuzzy_space_and_cytosol" variable_2="Rate_Ca_movement_between_the_fuzzy_space_and_cytosol"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="Ca_movement_between_the_fuzzy_space_and_cytosol" component_2="cytosolic_calcium"/>
      <map_variables variable_1="Ca_c" variable_2="Ca_c"/>
      <map_variables variable_1="Rate_Ca_movement_between_the_fuzzy_space_and_cytosol" variable_2="Rate_Ca_movement_between_the_fuzzy_space_and_cytosol"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="Ca_uptake_by_SR_Ca_ATPase" component_2="cytosolic_calcium"/>
      <map_variables variable_1="Ca_c" variable_2="Ca_c"/>
      <map_variables variable_1="Rate_Ca_uptake_by_SR_Ca_ATPase" variable_2="Rate_Ca_uptake_by_SR_Ca_ATPase"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="Ca_uptake_by_SR_Ca_ATPase" component_2="SR_calcium"/>
      <map_variables variable_1="Ca_s" variable_2="Ca_s"/>
      <map_variables variable_1="Rate_Ca_uptake_by_SR_Ca_ATPase" variable_2="Rate_Ca_uptake_by_SR_Ca_ATPase"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="Ca_buffering_in_the_SR" component_2="SR_calcium"/>
      <map_variables variable_1="Ca_s" variable_2="Ca_s"/>
      <map_variables variable_1="Rate_Ca_buffering_in_the_SR" variable_2="Rate_Ca_buffering_in_the_SR"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="Ca_buffering_in_the_SR" component_2="calsequestrin_bound_calcium"/>
      <map_variables variable_1="Ca_CSQ" variable_2="Ca_CSQ"/>
      <map_variables variable_1="Rate_Ca_buffering_in_the_SR" variable_2="Rate_Ca_buffering_in_the_SR"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="Ca_movement_through_the_SRRC" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="fuzzy_space_calcium" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="cytosolic_calcium" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="SR_calcium" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="calsequestrin_bound_calcium" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="SRRC_fast_activating_binding_site" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="SRRC_slow_inactivating_binding_site" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="Ca_bound_to_the_SRRC_fast_activating_binding_site" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection xmlns="http://www.cellml.org/cellml/1.0#">
      <map_components component_1="Ca_bound_to_the_SRRC_slow_inactivating_binding_site" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
</model>