- Author:
- pmr2.import <nobody@models.cellml.org>
- Date:
- 2007-05-23 05:31:15+12:00
- Desc:
- committing version02 of saucerman_brunton_michailova_mcculloch_2003
- Permanent Source URI:
- https://models.physiomeproject.org/workspace/saucerman_brunton_michailova_mcculloch_2003/rawfile/e5422b92bc7cca4e36c92f8866d76b74f22662c3/saucerman_brunton_michailova_mcculloch_2003.cellml
<?xml version='1.0' encoding='utf-8'?>
<!-- FILE : saucerman_model_2003.xml
CREATED : 18th November 2003
LAST MODIFIED : 18th November 2003
AUTHOR : Catherine Lloyd
Bioengineering Institute
The University of Auckland
MODEL STATUS : This model conforms to the CellML 1.0 Specification released on
10th August 2001, and the 16/1/02 CellML Metadata 1.0 Specification.
DESCRIPTION : This file contains a CellML description of Saucerman et al.'s 2003 mathematical model of beta-adrenergic control of cardiac myocyte contractility.
CHANGES:
--><model xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:cellml="http://www.cellml.org/cellml/1.0#" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" cmeta:id="saucerman_model_2003" name="saucerman_2003_version01">
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
<articleinfo>
<title>Modelling Beta-adrenergic Control of Cardiac Myocyte Contractility in Silico</title>
<author>
<firstname>Catherine</firstname>
<surname>Lloyd</surname>
<affiliation>
<shortaffil>Bioengineering Institute, University of Auckland</shortaffil>
</affiliation>
</author>
</articleinfo>
<section id="sec_status">
<title>Model Status</title>
<para>
This CellML version of the model can not be run in either PCEnv or COR at the moment due to the presence of circular arguments.
</para>
</section>
<sect1 id="sec_structure">
<title>Model Structure</title>
<para>
In cardiac myocytes, the beta-adrenergic signalling network is triggered in response to norepinephrine and epinephrine binding to the Gs-protein coupled beta-adrenergic receptor. Receptor-ligand binding activates the coupled Gs-protein, which in turn activates adenylate cyclase. ATP is converted into the secondary messenger cyclic AMP (cAMP). cAMP promotes dissociation of the protein kinase A (PKA) holoenzyme, whose catalytic subunits go on to phosphorylate a wide range of target proteins, two of which are the L-type calcium channel and phospholamban - which both play essential roles in the regulation of calcium dynamics and transport. In healthy cardiac myocytes, the end result of this signal transduction pathway is to provide coordinated control of contractility, metabolism and gene regulation. However, altered beta-adrenergic signalling may also play a role in the progression of heart failure.
</para>
<para>
Molecular components of this signalling pathway have been studied in detail as potential therapeutic targets in heart failure. However, the situation is complex, with intracellular compartmentation and functional integration of pathways being suggested as playing an important role in the signalling outcome. Systems level mathematical modelling will help make sense of large quantities of experimental data. Since the Hodgkin-Huxley Squid Axon Model was published in 1952, many ionic models of cell electrophysiology have been developed. Recently systems modelling has been identified as a potential method for improving understanding of signaling networks and the response to genetic and pharmaceutical perturbations.
</para>
<para>
While to date portions of the mammalian beta-adrenergic signalling network have been described by mathematical models of neuron and HEK-293 cells, no analysis has modelled and validated an entire pathway from ligand to effectors such that neuro-hormonal regulation of cell physiology can be predicted. In the publication described here, Jeffrey Saucerman <emphasis>et al.</emphasis> have developed a systems biology model to investigate the molecular mechanisms underlying the control of the beta-adrenergic signalling network over cardiac myocyte contractility (see figure below). The model was validated using a wide range of experimental data, and model simulations were used to investigate quantitatively the effects of specific molecular perturbations (with potential for pharmaceutical testing). By performing systems analysis the effects of molecular perturbations in the beta-adrenergic signalling network may be understood within the context of integrative physiology.
</para>
<para>
The model of the signal transduction pathway was embedded within an extension of the Luo-Rudy Ventricular Model II (dynamic), 1994, which was modified for the rabbit ventricular myocyte (Puglisi-Bers Rabbit Ventricular Myocyte Model, 2001). The model was further adapted by inserting the equations for the L-type calcium channel from the Jafri-Rice-Winslow Ventricular Model, 1998, and also the steady-state potassium currents from the Pandit <emphasis>et al.</emphasis> Adult Rat Left Ventricular Myocyte Model, 2001.
</para>
<para>
The complete original paper reference is cited below:
</para>
<para>
<ulink url="http://www.jbc.org/cgi/content/abstract/278/48/47997">Modeling beta-adrenergic control of cardiac myocyte contractility <emphasis>in silico</emphasis>
</ulink>, Jeffrey J. Saucerman, Laurence L. Brunton, Anushka P. Michailova, and Andrew D. McCulloch, 2003, <ulink url="http://www.jbc.org/">
<emphasis>Journal of Biological Chemistry</emphasis>
</ulink>, 48, 47997-48003. (<ulink url="http://www.jbc.org/cgi/content/full/278/48/47997">Full text (HTML)</ulink> and <ulink url="http://www.jbc.org/cgi/reprint/278/48/47997.pdf">PDF</ulink> versions of the article are available on the <emphasis>Journal of Biological Chemistry</emphasis> website.) <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&db=PubMed&list_uids=12972422&dopt=Abstract">PubMed ID: 12972422</ulink>
</para>
<informalfigure float="0" id="fig_cell_diagram">
<mediaobject>
<imageobject>
<objectinfo>
<title>cell diagram</title>
</objectinfo>
<imagedata fileref="saucerman_mcculloch_2003.png"/>
</imageobject>
</mediaobject>
<caption>Schematic diagram of the integrated model components, including the beta-adrenergic network, calcium handling, and the electrophysiology of the rat ventricular myocyte.</caption>
</informalfigure>
<informalfigure float="0" id="fig_pathway_diagram">
<mediaobject>
<imageobject>
<objectinfo>
<title>pathway diagram</title>
</objectinfo>
<imagedata fileref="saucerman_2003.png"/>
</imageobject>
</mediaobject>
<caption>Schematic diagram highlighting the signalling network component of the model.</caption>
</informalfigure>
</sect1>
</article>
</documentation>
<units name="micromolar">
<unit units="mole" prefix="micro"/>
<unit units="litre" exponent="-1"/>
</units>
<units name="millimolar">
<unit units="mole" prefix="milli"/>
<unit units="litre" exponent="-1"/>
</units>
<units name="millivolt">
<unit units="volt" prefix="milli"/>
</units>
<units name="microlitre">
<unit units="litre" prefix="micro"/>
</units>
<units name="millisecond">
<unit units="second" prefix="milli"/>
</units>
<units name="flux">
<unit units="micromolar"/>
<unit units="second" exponent="-1"/>
</units>
<units name="millimolar_per_second">
<unit units="millimolar"/>
<unit units="second" exponent="-1"/>
</units>
<units name="first_order_rate_constant">
<unit units="second" exponent="-1"/>
</units>
<units name="second_order_rate_constant">
<unit units="micromolar" exponent="-1"/>
<unit units="second" exponent="-1"/>
</units>
<units name="per_millisecond">
<unit units="second" prefix="milli" exponent="-1"/>
</units>
<units name="per_millivolt">
<unit units="volt" prefix="milli" exponent="-1"/>
</units>
<units name="per_millivolt_second">
<unit units="volt" prefix="milli" exponent="-1"/>
<unit units="second" exponent="-1"/>
</units>
<units name="cm_per_second">
<unit units="metre" prefix="centi"/>
<unit units="second" exponent="-1"/>
</units>
<units name="cm2">
<unit units="metre" prefix="centi" exponent="2"/>
</units>
<units name="milliS_per_microF">
<unit units="siemens" prefix="milli"/>
<unit units="farad" prefix="micro" exponent="-1"/>
</units>
<units name="microA_per_microF">
<unit units="ampere" prefix="micro"/>
<unit units="farad" prefix="micro"/>
</units>
<units name="millimolar3_per_second">
<unit units="millimolar" exponent="-3"/>
<unit units="second" exponent="-1"/>
</units>
<units name="millimolar4_per_second">
<unit units="millimolar" exponent="-4"/>
<unit units="second" exponent="-1"/>
</units>
<units name="millimolar_4">
<unit units="millimolar" exponent="-4"/>
</units>
<units name="joule_per_kilomole_kelvin">
<unit units="joule"/>
<unit units="mole" prefix="kilo" exponent="-1"/>
<unit units="kelvin" exponent="-1"/>
</units>
<units name="coulomb_per_mole">
<unit units="coulomb" exponent="-1"/>
<unit units="mole"/>
</units>
<units name="pico_litre">
<unit units="litre" prefix="pico"/>
</units>
<component name="environment">
<variable units="second" public_interface="out" name="time"/>
</component>
<!-- Signal transduction pathway component of the model from the Saucerman et al. 2003 paper -->
<component name="beta_adrenergic_receptor_module">
<variable units="micromolar" public_interface="out" name="LRG"/>
<variable units="micromolar" public_interface="out" name="RG"/>
<variable units="micromolar" name="LR"/>
<variable units="micromolar" name="beta1_AR" initial_value="0.01205"/>
<variable units="micromolar" name="Gs" initial_value="3.829"/>
<variable units="micromolar" name="beta1_ARact" initial_value="0.01205"/>
<variable units="micromolar" name="beta1_AR_S464" initial_value="0.0"/>
<variable units="micromolar" name="beta1_AR_S301" initial_value="1.154E-3"/>
<variable units="micromolar" name="L" initial_value="10.0"/>
<variable units="micromolar" name="Ltotmax" initial_value="1.0"/>
<variable units="micromolar" name="KL" initial_value="0.285"/>
<variable units="micromolar" name="KR" initial_value="0.062"/>
<variable units="micromolar" name="KC" initial_value="33.0"/>
<variable units="micromolar" name="Gstot" initial_value="3.83"/>
<variable units="first_order_rate_constant" name="k_betaARK_plus" initial_value="1.1E-3"/>
<variable units="first_order_rate_constant" name="k_betaARK_minus" initial_value="2.2E-3"/>
<variable units="second_order_rate_constant" name="k_PKA_plus" initial_value="3.6E-3"/>
<variable units="first_order_rate_constant" name="k_PKA_minus" initial_value="2.2E-3"/>
<variable units="second" public_interface="in" name="time"/>
<variable units="micromolar" public_interface="in" name="Gs_beta_gamma"/>
<variable units="micromolar" public_interface="in" name="PKACI"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="1">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> L </ci>
</apply>
<apply>
<minus/>
<ci> Ltotmax </ci>
<apply>
<plus/>
<ci> LR </ci>
<ci> LRG </ci>
<ci> L </ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>LR</ci>
<apply>
<divide/>
<apply>
<times/>
<ci> L </ci>
<ci> beta1_AR </ci>
</apply>
<ci> KL </ci>
</apply>
</apply>
<apply>
<eq/>
<ci>LRG</ci>
<apply>
<divide/>
<apply>
<times/>
<ci> L </ci>
<ci> beta1_AR </ci>
<ci> Gs </ci>
</apply>
<apply>
<times/>
<ci> KL </ci>
<ci> KR </ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>RG</ci>
<apply>
<divide/>
<apply>
<times/>
<ci> beta1_AR </ci>
<ci> Gs </ci>
</apply>
<ci> KC </ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> Gs </ci>
</apply>
<apply>
<minus/>
<ci> Gstot </ci>
<apply>
<plus/>
<ci> RG </ci>
<ci> LRG </ci>
<ci> Gs </ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> beta1_AR </ci>
</apply>
<apply>
<minus/>
<ci> beta1_ARact </ci>
<apply>
<plus/>
<ci> LR </ci>
<ci> LRG </ci>
<ci> RG </ci>
<ci> beta1_AR </ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>beta1_ARact</ci>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<times/>
<ci> k_betaARK_minus </ci>
<ci> beta1_AR_S464 </ci>
</apply>
<apply>
<times/>
<ci> k_betaARK_plus </ci>
<apply>
<plus/>
<ci> LR </ci>
<ci> LRG </ci>
</apply>
</apply>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci> k_PKA_minus </ci>
<ci> beta1_AR_S301 </ci>
</apply>
<apply>
<times/>
<ci> k_PKA_plus </ci>
<ci> PKACI </ci>
<ci> beta1_ARact </ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>beta1_AR_S464</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci> k_betaARK_plus </ci>
<apply>
<plus/>
<ci> LR </ci>
<ci> LRG </ci>
</apply>
</apply>
<apply>
<times/>
<ci> k_betaARK_minus </ci>
<ci> beta1_AR_S464 </ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>beta1_AR_S301</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci> k_PKA_plus </ci>
<ci> PKACI </ci>
<ci> beta1_ARact </ci>
</apply>
<apply>
<times/>
<ci> k_PKA_minus </ci>
<ci> beta1_AR_S301 </ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="Gs_activation_module">
<variable units="micromolar" public_interface="out" name="Gs_beta_gamma" initial_value="0.02569"/>
<variable units="micromolar" public_interface="out" name="Gs_alpha_GTPtot" initial_value="0.02505"/>
<variable units="micromolar" name="Gs_alpha_GDP" initial_value="6.44E-4"/>
<variable units="first_order_rate_constant" name="k_gact" initial_value="16.0"/>
<variable units="first_order_rate_constant" name="k_hyd" initial_value="0.8"/>
<variable units="second_order_rate_constant" name="k_reassoc" initial_value="1.21E3"/>
<variable units="micromolar" public_interface="in" name="RG"/>
<variable units="micromolar" public_interface="in" name="LRG"/>
<variable units="second" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="2">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Gs_alpha_GTPtot</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci> k_gact </ci>
<apply>
<plus/>
<ci> RG </ci>
<ci> LRG </ci>
</apply>
</apply>
<apply>
<times/>
<ci> k_hyd </ci>
<ci> Gs_alpha_GTPtot </ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Gs_beta_gamma</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci> k_gact </ci>
<apply>
<plus/>
<ci> RG </ci>
<ci> LRG </ci>
</apply>
</apply>
<apply>
<times/>
<ci> k_reassoc </ci>
<ci> Gs_alpha_GDP </ci>
<ci> Gs_beta_gamma </ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Gs_alpha_GDP</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci> k_hyd </ci>
<ci> Gs_alpha_GTPtot </ci>
</apply>
<apply>
<times/>
<ci> k_reassoc </ci>
<ci> Gs_alpha_GDP </ci>
<ci> Gs_beta_gamma </ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="cyclic_AMP_metabolism_module">
<variable units="micromolar" public_interface="out" name="cAMPtot" initial_value="0.8453"/>
<variable units="micromolar" name="AC" initial_value="0.04706"/>
<variable units="micromolar" name="Gs_alpha_GTP" initial_value="0.02241"/>
<variable units="micromolar" name="Gs_alpha_GTP_AC"/>
<variable units="micromolar" name="PDE" initial_value="0.0389"/>
<!-- <variable units="micromolar" name="PDE_ACT"/> -->
<variable units="micromolar" name="ACtot" initial_value="49.7E-3"/>
<variable units="micromolar" name="ATP" initial_value="5.0E3"/>
<variable units="micromolar" name="PDEtot" initial_value="38.9E-3"/>
<variable units="micromolar" name="PDEinhib"/>
<variable units="micromolar" name="IBMXtot" initial_value="0.0"/>
<variable units="micromolar" name="IBMX" initial_value="0.0"/>
<variable units="micromolar" name="fsktot" initial_value="0.0"/>
<variable units="micromolar" name="fsk" initial_value="0.0"/>
<variable units="micromolar" name="fsk_AC"/>
<variable units="first_order_rate_constant" name="kAC_basal" initial_value="0.2"/>
<variable units="first_order_rate_constant" name="kAC_fsk" initial_value="7.3"/>
<variable units="first_order_rate_constant" name="k_PDE" initial_value="5.0"/>
<variable units="first_order_rate_constant" name="kAC_Gs_alpha" initial_value="8.5"/>
<variable units="micromolar" name="Km_basal" initial_value="1.03E3"/>
<variable units="micromolar" name="Km_PDE" initial_value="1.3"/>
<variable units="micromolar" name="Km_fsk" initial_value="860.0"/>
<variable units="micromolar" name="K_fsk" initial_value="44.0"/>
<variable units="micromolar" name="Km_Gs_alpha_GTP" initial_value="0.4"/>
<variable units="micromolar" name="K_Gs_alpha" initial_value="0.4"/>
<variable units="micromolar" name="K_IBMX" initial_value="30.0"/>
<variable units="micromolar" public_interface="in" name="Gs_alpha_GTPtot"/>
<variable units="second" public_interface="in" name="time"/>
<variable units="micromolar" public_interface="in" name="cAMP"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="3">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> Gs_alpha_GTP </ci>
</apply>
<apply>
<minus/>
<ci> Gs_alpha_GTPtot </ci>
<apply>
<plus/>
<ci> Gs_alpha_GTP_AC </ci>
<ci> Gs_alpha_GTP </ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> fsk </ci>
</apply>
<apply>
<minus/>
<ci> fsktot </ci>
<apply>
<plus/>
<ci> fsk_AC </ci>
<ci> fsk </ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> AC </ci>
</apply>
<apply>
<minus/>
<ci> ACtot </ci>
<apply>
<plus/>
<ci> Gs_alpha_GTP_AC </ci>
<ci> AC </ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> PDE </ci>
</apply>
<apply>
<minus/>
<ci> PDEtot </ci>
<apply>
<plus/>
<ci> PDEinhib </ci>
<ci> PDE </ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> IBMX </ci>
</apply>
<apply>
<minus/>
<ci> IBMXtot </ci>
<apply>
<plus/>
<ci> PDEinhib </ci>
<ci> IBMX </ci>
</apply>
</apply>
</apply>
<!--
<apply><eq/>
<apply><diff/>
<bvar><ci> time </ci></bvar>
<ci> PDE_ACT </ci>
</apply>
<apply><divide/>
<apply><times/>
<ci> k_PDE </ci>
<ci> PDE </ci>
<ci> cAMP </ci>
</apply>
<apply><plus/>
<ci> Km_PDE </ci>
<ci> cAMP </ci>
</apply>
</apply>
</apply>
-->
<apply>
<eq/>
<ci>Gs_alpha_GTP_AC</ci>
<apply>
<divide/>
<apply>
<times/>
<ci> Gs_alpha_GTP </ci>
<ci> AC </ci>
</apply>
<ci> K_Gs_alpha </ci>
</apply>
</apply>
<apply>
<eq/>
<ci>fsk_AC</ci>
<apply>
<divide/>
<apply>
<times/>
<ci> fsk </ci>
<ci> AC </ci>
</apply>
<ci> K_fsk </ci>
</apply>
</apply>
<apply>
<eq/>
<ci>PDEinhib</ci>
<apply>
<divide/>
<apply>
<times/>
<ci> PDE </ci>
<ci> IBMX </ci>
</apply>
<ci> K_IBMX </ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>cAMPtot</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<times/>
<ci> kAC_basal </ci>
<ci> AC </ci>
<ci> ATP </ci>
</apply>
<apply>
<plus/>
<ci> Km_basal </ci>
<ci> ATP </ci>
</apply>
</apply>
<apply>
<divide/>
<apply>
<times/>
<ci> kAC_Gs_alpha </ci>
<ci> Gs_alpha_GTP_AC </ci>
<ci> ATP </ci>
</apply>
<apply>
<plus/>
<ci> Km_Gs_alpha_GTP </ci>
<ci> ATP </ci>
</apply>
</apply>
<apply>
<divide/>
<apply>
<times/>
<ci> kAC_fsk </ci>
<ci> fsk_AC </ci>
<ci> ATP </ci>
</apply>
<apply>
<plus/>
<ci> Km_fsk </ci>
<ci> ATP </ci>
</apply>
</apply>
</apply>
<apply>
<divide/>
<apply>
<times/>
<ci> k_PDE </ci>
<ci> PDE </ci>
<ci> cAMP </ci>
</apply>
<apply>
<plus/>
<ci> Km_PDE </ci>
<ci> cAMP </ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="PKA_activation_module">
<variable units="micromolar" public_interface="out" name="PKACI"/>
<variable units="micromolar" public_interface="out" name="PKACII"/>
<variable units="micromolar" public_interface="out" name="cAMP"/>
<variable units="micromolar" name="ARCI"/>
<variable units="micromolar" name="A2RCI"/>
<variable units="micromolar" name="A2RCII"/>
<variable units="micromolar" name="A2RI"/>
<variable units="micromolar" name="A2RII"/>
<variable units="micromolar" name="ARCII"/>
<variable units="micromolar" name="RCI"/>
<variable units="micromolar" name="PKACI_PKI"/>
<variable units="micromolar" name="PKAtemp"/>
<variable units="micromolar" name="PKI"/>
<variable units="micromolar" name="RCII"/>
<variable units="micromolar" name="PKACII_PKI"/>
<variable units="micromolar" name="Ki_pki" initial_value="0.2E-3"/>
<variable units="micromolar" name="PKAItot" initial_value="0.59"/>
<variable units="micromolar" name="PKAIItot" initial_value="0.025"/>
<variable units="micromolar" name="PKItot" initial_value="0.18"/>
<variable units="micromolar" name="KA" initial_value="9.14"/>
<variable units="micromolar" name="KB" initial_value="1.64"/>
<variable units="micromolar" name="KD" initial_value="4.375"/>
<variable units="micromolar" name="KPKI" initial_value="2.0E-4"/>
<variable units="micromolar" public_interface="in" name="cAMPtot"/>
<variable units="second" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="4">
<apply>
<eq/>
<ci>cAMP</ci>
<apply>
<minus/>
<ci> cAMPtot </ci>
<apply>
<plus/>
<apply>
<plus/>
<ci> ARCI </ci>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> A2RCI </ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> A2RI </ci>
</apply>
</apply>
<apply>
<plus/>
<ci> ARCII </ci>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> A2RCII </ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> A2RII </ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci> PKAtemp </ci>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<times/>
<ci> KA </ci>
<ci> KB </ci>
</apply>
<ci> KD </ci>
</apply>
<apply>
<divide/>
<apply>
<times/>
<ci> KA </ci>
<ci> cAMP </ci>
</apply>
<ci> KD </ci>
</apply>
<apply>
<divide/>
<apply>
<power/>
<ci> cAMP </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
<ci> KD </ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci> PKACI </ci>
<apply>
<minus/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> PKAItot </ci>
<apply>
<power/>
<ci> cAMP </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci> PKAtemp </ci>
<ci> PKACI </ci>
</apply>
<apply>
<power/>
<ci> cAMP </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci> PKACI </ci>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<divide/>
<ci> PKI </ci>
<ci> Ki_pki </ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci> PKACII </ci>
<apply>
<minus/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> PKAIItot </ci>
<apply>
<power/>
<ci> cAMP </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci> PKAtemp </ci>
<ci> PKACII </ci>
</apply>
<apply>
<power/>
<ci> cAMP </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci> PKACII </ci>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<divide/>
<ci> PKI </ci>
<ci> Ki_pki </ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>RCI</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci> KA </ci>
<ci> KB </ci>
</apply>
<apply>
<power/>
<ci> cAMP </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
</apply>
<apply>
<times/>
<apply>
<divide/>
<ci> PKACI </ci>
<ci> KD </ci>
</apply>
<ci> PKACI </ci>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<divide/>
<ci> PKI </ci>
<ci> Ki_pki </ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>ARCI</ci>
<apply>
<times/>
<apply>
<divide/>
<ci> KA </ci>
<ci> cAMP </ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<ci> PKACI </ci>
<ci> KD </ci>
</apply>
<ci> PKACI </ci>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<divide/>
<ci> PKI </ci>
<ci> Ki_pki </ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>ARCII</ci>
<apply>
<times/>
<apply>
<divide/>
<ci> KA </ci>
<ci> cAMP </ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<ci> PKACII </ci>
<ci> KD </ci>
</apply>
<ci> PKACII </ci>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<divide/>
<ci> PKI </ci>
<ci> Ki_pki </ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>A2RCI</ci>
<apply>
<times/>
<apply>
<divide/>
<ci> PKACI </ci>
<ci> KD </ci>
</apply>
<ci> PKACI </ci>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<divide/>
<ci> PKI </ci>
<ci> Ki_pki </ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>A2RI</ci>
<apply>
<times/>
<ci> PKACI </ci>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<divide/>
<ci> PKI </ci>
<ci> Ki_pki </ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>RCII</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci> KA </ci>
<ci> KB </ci>
</apply>
<apply>
<power/>
<ci> cAMP </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
</apply>
<apply>
<times/>
<apply>
<divide/>
<ci> PKACII </ci>
<ci> KD </ci>
</apply>
<ci> PKACII </ci>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<divide/>
<ci> PKI </ci>
<ci> Ki_pki </ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>A2RCII</ci>
<apply>
<times/>
<apply>
<divide/>
<ci> PKACII </ci>
<ci> KD </ci>
</apply>
<ci> PKACII </ci>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<divide/>
<ci> PKI </ci>
<ci> Ki_pki </ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>A2RII</ci>
<apply>
<times/>
<ci> PKACII </ci>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<divide/>
<ci> PKI </ci>
<ci> Ki_pki </ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>PKI</ci>
<apply>
<divide/>
<apply>
<times/>
<ci> KPKI </ci>
<ci> PKItot </ci>
</apply>
<apply>
<plus/>
<ci> KPKI </ci>
<ci> PKACI </ci>
<ci> PKACII </ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>PKACI_PKI</ci>
<apply>
<divide/>
<apply>
<times/>
<ci> PKACI </ci>
<ci> PKItot </ci>
</apply>
<apply>
<plus/>
<ci> KPKI </ci>
<ci> PKACI </ci>
<ci> PKACII </ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>PKACII_PKI</ci>
<apply>
<divide/>
<apply>
<times/>
<ci> PKACI </ci>
<ci> PKItot </ci>
</apply>
<apply>
<plus/>
<ci> KPKI </ci>
<ci> PKACI </ci>
<ci> PKACII </ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="phospholamban_regulation_module">
<variable units="dimensionless" public_interface="out" name="fracPLBp"/>
<variable units="dimensionless" public_interface="out" name="fracPLB"/>
<variable units="dimensionless" public_interface="out" name="fracPLBo" initial_value="0.9613"/>
<variable units="micromolar" name="PLBp" initial_value="4.105"/>
<variable units="micromolar" name="PLB"/>
<variable units="micromolar" name="Inhib1ptot" initial_value="0.0526"/>
<variable units="micromolar" name="Inhib1" initial_value="0.2474"/>
<variable units="micromolar" name="Inhib1p" initial_value="6.3E-5"/>
<variable units="micromolar" name="PP1_Inhib1p"/>
<variable units="micromolar" name="PP1" initial_value="0.838"/>
<variable units="micromolar" name="PP1tot" initial_value="0.89"/>
<variable units="micromolar" name="PLBtot" initial_value="106.0"/>
<variable units="micromolar" name="Inhib1tot" initial_value="0.3"/>
<variable units="dimensionless" name="epsilon" initial_value="10.0"/>
<variable units="first_order_rate_constant" name="kPKA_PLB" initial_value="54.0"/>
<variable units="micromolar" name="KmPKA_PLB" initial_value="21.0"/>
<variable units="first_order_rate_constant" name="kPKA_Inhib1" initial_value="60.0"/>
<variable units="first_order_rate_constant" name="kPP1_PLB" initial_value="8.5"/>
<variable units="micromolar" name="KmPP1_PLB" initial_value="7.0"/>
<variable units="micromolar" name="KmPKA_Inhib1" initial_value="1.0"/>
<variable units="flux" name="VmaxPP2A_Inhib1" initial_value="14.0"/>
<variable units="micromolar" name="KmPP2A_Inhib1" initial_value="1.0"/>
<variable units="micromolar" name="KInhib1" initial_value="1.0E-3"/>
<variable units="second" public_interface="in" name="time"/>
<variable units="micromolar" public_interface="in" name="PKACI"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="5">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>PP1</ci>
</apply>
<apply>
<minus/>
<ci> PP1tot </ci>
<apply>
<plus/>
<ci> PP1_Inhib1p </ci>
<ci> PP1 </ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>PP1_Inhib1p</ci>
<apply>
<divide/>
<apply>
<times/>
<ci> PP1 </ci>
<ci> Inhib1p </ci>
</apply>
<ci> KInhib1 </ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Inhib1p</ci>
</apply>
<apply>
<minus/>
<ci> Inhib1ptot </ci>
<apply>
<plus/>
<ci> PP1_Inhib1p </ci>
<ci> Inhib1p </ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>PLBp</ci>
</apply>
<apply>
<minus/>
<apply>
<divide/>
<apply>
<times/>
<ci> kPKA_PLB </ci>
<ci> PKACI </ci>
<ci> PLB </ci>
</apply>
<apply>
<plus/>
<ci> KmPKA_PLB </ci>
<ci> PLB </ci>
</apply>
</apply>
<apply>
<divide/>
<apply>
<times/>
<ci> kPP1_PLB </ci>
<ci> PP1 </ci>
<ci> PLBp </ci>
</apply>
<apply>
<plus/>
<ci> KmPP1_PLB </ci>
<ci> PLBp </ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>PLB</ci>
<apply>
<minus/>
<ci> PLBtot </ci>
<ci> PLBp </ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Inhib1ptot</ci>
</apply>
<apply>
<minus/>
<apply>
<divide/>
<apply>
<times/>
<ci> kPKA_Inhib1 </ci>
<ci> PKACI </ci>
<ci> Inhib1ptot </ci>
</apply>
<apply>
<plus/>
<ci> KmPKA_Inhib1 </ci>
<ci> Inhib1ptot </ci>
</apply>
</apply>
<apply>
<divide/>
<apply>
<times/>
<ci> VmaxPP2A_Inhib1 </ci>
<ci> Inhib1 </ci>
</apply>
<apply>
<plus/>
<ci> KmPP2A_Inhib1 </ci>
<ci> Inhib1 </ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Inhib1</ci>
</apply>
<apply>
<minus/>
<apply>
<divide/>
<apply>
<times/>
<ci> VmaxPP2A_Inhib1 </ci>
<ci> Inhib1 </ci>
</apply>
<apply>
<plus/>
<ci> KmPP2A_Inhib1 </ci>
<ci> Inhib1 </ci>
</apply>
</apply>
<apply>
<divide/>
<apply>
<times/>
<ci> kPKA_Inhib1 </ci>
<ci> PKACI </ci>
<ci> Inhib1ptot </ci>
</apply>
<apply>
<plus/>
<ci> KmPKA_Inhib1 </ci>
<ci> Inhib1ptot </ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci> fracPLBp </ci>
<apply>
<divide/>
<ci> PLBp </ci>
<ci> PLBtot </ci>
</apply>
</apply>
<apply>
<eq/>
<ci> fracPLB </ci>
<apply>
<divide/>
<ci> PLB </ci>
<ci> PLBtot </ci>
</apply>
</apply>
</math>
</component>
<component name="L_type_Ca_channel_regulation_module">
<variable units="dimensionless" public_interface="out" name="fracLCCap"/>
<variable units="dimensionless" public_interface="out" name="fracLCCbp"/>
<variable units="dimensionless" public_interface="out" name="fracLCCapo" initial_value="0.2041"/>
<variable units="dimensionless" public_interface="out" name="fracLCCbpo" initial_value="0.2336"/>
<variable units="micromolar" name="LCC_beta2_p" initial_value="1.27575E-4"/>
<variable units="micromolar" name="LCC_beta2" initial_value="1.46025E-3"/>
<!--
<variable units="micromolar" name="P1"/>
<variable units="micromolar" name="P2"/>
<variable units="micromolar" name="P3"/>
<variable units="micromolar" name="P4"/>
<variable units="micromolar" name="P5"/>
<variable units="micromolar" name="P6"/>
<variable units="micromolar" name="P7"/>
<variable units="micromolar" name="P8"/>
<variable units="micromolar" name="P9"/>
<variable units="micromolar" name="P10"/>
<variable units="micromolar" name="P11"/>
<variable units="micromolar" name="P12"/>
<variable units="micromolar" name="P13"/>
<variable units="micromolar" name="P14"/>
<variable units="micromolar" name="P15"/>
<variable units="micromolar" name="P16"/>
<variable units="micromolar" name="P17"/>
<variable units="micromolar" name="P18"/>
<variable units="micromolar" name="P19"/>
<variable units="micromolar" name="P20"/>
<variable units="micromolar" name="P21"/>
<variable units="micromolar" name="P22"/>
<variable units="micromolar" name="P23"/>
<variable units="micromolar" name="P24"/>
<variable units="micromolar" name="PHOSPH1"/>
<variable units="micromolar" name="PHOSPH2"/>
<variable units="micromolar" name="PHOSPH3"/>
<variable units="micromolar" name="PHOSPH4"/>
<variable units="micromolar" name="PHOSPH5"/>
<variable units="micromolar" name="PHOSPH6"/>
<variable units="micromolar" name="PHOSPH7"/>
<variable units="micromolar" name="PHOSPH8"/>
<variable units="micromolar" name="PHOSPH9"/>
<variable units="micromolar" name="PHOSPH10"/>
<variable units="micromolar" name="PHOSPH11"/>
<variable units="micromolar" name="PHOSPH12"/>
-->
<variable units="micromolar" name="LCCtot" initial_value="25.0E-3"/>
<variable units="micromolar" name="PP1_LCC" initial_value="25.0E-3"/>
<variable units="micromolar" name="PP2A_LCC" initial_value="25E-3"/>
<variable units="dimensionless" name="epsilon" initial_value="10.0"/>
<variable units="first_order_rate_constant" name="kPKA_LCC" initial_value="54.0"/>
<variable units="micromolar" name="KmPKA_LCC" initial_value="21.0"/>
<variable units="micromolar" name="KmPP1_LCC" initial_value="3.0"/>
<variable units="first_order_rate_constant" name="kPP1_LCC" initial_value="8.52"/>
<variable units="first_order_rate_constant" name="kPP2A_LCC" initial_value="10.1"/>
<variable units="micromolar" name="KmPP2A_LCC" initial_value="3.0"/>
<variable units="first_order_rate_constant" name="omega_LCC" initial_value="10.0"/>
<variable units="first_order_rate_constant" name="gamma0_LCC" initial_value="2.8E3"/>
<variable units="first_order_rate_constant" name="fn_LCC" initial_value="200.0"/>
<variable units="first_order_rate_constant" name="gn_LCC" initial_value="2.0E3"/>
<variable units="first_order_rate_constant" name="fc_LCC" initial_value="5.0"/>
<variable units="first_order_rate_constant" name="gc_LCC" initial_value="7.0E3"/>
<variable units="first_order_rate_constant" name="fp_LCC" initial_value="800.0"/>
<variable units="first_order_rate_constant" name="gp_LCC" initial_value="2.0E3"/>
<variable units="second" public_interface="in" name="time"/>
<variable units="micromolar" public_interface="in" name="PKACII"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="6">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>LCC_beta2_p</ci>
</apply>
<apply>
<minus/>
<apply>
<divide/>
<apply>
<times/>
<ci> epsilon </ci>
<ci> kPKA_LCC </ci>
<ci> PKACII </ci>
<apply>
<minus/>
<ci> LCCtot </ci>
<ci> LCC_beta2_p </ci>
</apply>
</apply>
<apply>
<plus/>
<ci> KmPKA_LCC </ci>
<apply>
<times/>
<ci> epsilon </ci>
<apply>
<minus/>
<ci> LCCtot </ci>
<ci> LCC_beta2_p </ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<divide/>
<apply>
<times/>
<ci> epsilon </ci>
<ci> kPP2A_LCC </ci>
<ci> PP2A_LCC </ci>
<ci> LCC_beta2_p </ci>
</apply>
<apply>
<plus/>
<ci> KmPP2A_LCC </ci>
<apply>
<times/>
<ci> epsilon </ci>
<ci> LCC_beta2_p </ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>LCC_beta2</ci>
</apply>
<apply>
<minus/>
<apply>
<divide/>
<apply>
<times/>
<ci> epsilon </ci>
<ci> kPKA_LCC </ci>
<ci> PKACII </ci>
<apply>
<minus/>
<ci> LCCtot </ci>
<ci> LCC_beta2 </ci>
</apply>
</apply>
<apply>
<plus/>
<ci> KmPKA_LCC </ci>
<apply>
<times/>
<ci> epsilon </ci>
<apply>
<minus/>
<ci> LCCtot </ci>
<ci> LCC_beta2 </ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<divide/>
<apply>
<times/>
<ci> epsilon </ci>
<ci> kPP1_LCC </ci>
<ci> PP1_LCC </ci>
<ci> LCC_beta2 </ci>
</apply>
<apply>
<plus/>
<ci> KmPP1_LCC </ci>
<apply>
<times/>
<ci> epsilon </ci>
<ci> LCC_beta2 </ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci> fracLCCbp </ci>
<apply>
<divide/>
<ci> LCC_beta2 </ci>
<ci> LCCtot </ci>
</apply>
</apply>
<apply>
<eq/>
<ci> fracLCCap </ci>
<apply>
<divide/>
<ci> LCC_beta2_p </ci>
<ci> LCCtot </ci>
</apply>
</apply>
<!-- Mode Normal -->
<!--
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>P1</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<ci> beta_LCC </ci>
<ci> P2 </ci>
</apply>
<apply>
<times/>
<ci> omega_LCC </ci>
<ci> P7 </ci>
</apply>
</apply>
<apply>
<plus/>
<apply>
<times/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 4.0 </cn>
<ci> alpha_LCC </ci>
</apply>
<ci> gamma_LCC </ci>
</apply>
<ci> P1 </ci>
</apply>
<ci> PHOSPH1 </ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>P2</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 4.0 </cn>
<ci> alpha_LCC </ci>
<ci> P1 </ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> beta_LCC </ci>
<ci> P3 </ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<ci> omega_LCC </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
<ci> P8 </ci>
</apply>
</apply>
<apply>
<plus/>
<apply>
<times/>
<apply>
<plus/>
<ci> beta_LCC </ci>
<apply>
<times/>
<cn cellml:units="dimensionless"> 3.0 </cn>
<ci> alpha_LCC </ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> gamma_LCC </ci>
</apply>
</apply>
<ci> P2 </ci>
</apply>
<ci> PHOSPH2 </ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>P3</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 3.0 </cn>
<ci> alpha_LCC </ci>
<ci> P2 </ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 3.0 </cn>
<ci> beta_LCC </ci>
<ci> P4 </ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<ci> omega_LCC </ci>
<apply>
<power/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
</apply>
<ci> P9 </ci>
</apply>
</apply>
<apply>
<plus/>
<apply>
<times/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> beta_LCC </ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> alpha_LCC </ci>
</apply>
<apply>
<times/>
<apply>
<power/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
<ci> gamma_LCC </ci>
</apply>
</apply>
<ci> P3 </ci>
</apply>
<ci> PHOSPH3 </ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>P4</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> alpha_LCC </ci>
<ci> P3 </ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 4.0 </cn>
<ci> beta_LCC </ci>
<ci> P5 </ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<ci> omega_LCC </ci>
<apply>
<power/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<cn cellml:units="dimensionless"> 3.0 </cn>
</apply>
</apply>
<ci> P10 </ci>
</apply>
</apply>
<apply>
<plus/>
<apply>
<times/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 3.0 </cn>
<ci> beta_LCC </ci>
</apply>
<ci> alpha_LCC </ci>
<apply>
<times/>
<apply>
<power/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<cn cellml:units="dimensionless"> 3.0 </cn>
</apply>
<ci> gamma_LCC </ci>
</apply>
</apply>
<ci> P4 </ci>
</apply>
<ci> PHOSPH4 </ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>P5</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<ci> alpha_LCC </ci>
<ci> P4 </ci>
</apply>
<apply>
<times/>
<ci> gn_LCC </ci>
<ci> P6 </ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<ci> omega_LCC </ci>
<apply>
<power/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<cn cellml:units="dimensionless"> 4.0 </cn>
</apply>
</apply>
<ci> P11 </ci>
</apply>
</apply>
<apply>
<plus/>
<apply>
<times/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 4.0 </cn>
<ci> beta_LCC </ci>
</apply>
<ci> fn_LCC </ci>
<apply>
<times/>
<apply>
<power/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<cn cellml:units="dimensionless"> 4.0 </cn>
</apply>
<ci> gamma_LCC </ci>
</apply>
</apply>
<ci> P5 </ci>
</apply>
<ci> PHOSPH5 </ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>P6</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci> fc_LCC </ci>
<ci> P5 </ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci> gc_LCC </ci>
<ci> P6 </ci>
</apply>
<ci> PHOSPH6 </ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>P7</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<apply>
<divide/>
<ci> beta_LCC </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
<ci> P8 </ci>
</apply>
<apply>
<times/>
<ci> gamma_LCC </ci>
<ci> P1 </ci>
</apply>
</apply>
<apply>
<plus/>
<apply>
<times/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 4.0 </cn>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> alpha_LCC </ci>
</apply>
<ci> omega_LCC </ci>
</apply>
<ci> P7 </ci>
</apply>
<ci> PHOSPH7 </ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>P8</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 4.0 </cn>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> alpha_LCC </ci>
<ci> P7 </ci>
</apply>
<apply>
<times/>
<ci> beta_LCC </ci>
<ci> P9 </ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> gamma_LCC </ci>
<ci> P2 </ci>
</apply>
</apply>
<apply>
<plus/>
<apply>
<times/>
<apply>
<plus/>
<apply>
<divide/>
<ci> beta_LCC </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 3.0 </cn>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> alpha_LCC </ci>
</apply>
<apply>
<divide/>
<ci> omega_LCC </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
</apply>
<ci> P8 </ci>
</apply>
<ci> PHOSPH8 </ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>P9</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 3.0 </cn>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> alpha_LCC </ci>
<ci> P8 </ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 3.0 </cn>
<ci> beta_LCC </ci>
</apply>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
<ci> P10 </ci>
</apply>
<apply>
<times/>
<apply>
<power/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
<ci> gamma_LCC </ci>
<ci> P3 </ci>
</apply>
</apply>
<apply>
<plus/>
<apply>
<times/>
<apply>
<plus/>
<ci> beta_LCC </ci>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> alpha_LCC </ci>
</apply>
<apply>
<divide/>
<ci> omega_LCC </ci>
<apply>
<power/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
</apply>
</apply>
<ci> P9 </ci>
</apply>
<ci> PHOSPH9 </ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>P10</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> alpha_LCC </ci>
<ci> P9 </ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> beta_LCC </ci>
<ci> P11 </ci>
</apply>
<apply>
<times/>
<apply>
<power/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<cn cellml:units="dimensionless"> 3.0 </cn>
</apply>
<ci> gamma_LCC </ci>
<ci> P4 </ci>
</apply>
</apply>
<apply>
<plus/>
<apply>
<times/>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 3.0 </cn>
<ci> beta_LCC </ci>
</apply>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> alpha_LCC </ci>
</apply>
<apply>
<divide/>
<ci> omega_LCC </ci>
<apply>
<power/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<cn cellml:units="dimensionless"> 3.0 </cn>
</apply>
</apply>
</apply>
<ci> P10 </ci>
</apply>
<ci> PHOSPH10 </ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>P11</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> alpha_LCC </ci>
<ci> P10 </ci>
</apply>
<apply>
<times/>
<ci> gc_LCC </ci>
<ci> P12 </ci>
</apply>
<apply>
<times/>
<apply>
<power/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<cn cellml:units="dimensionless"> 4.0 </cn>
</apply>
<ci> gamma_LCC </ci>
<ci> P5 </ci>
</apply>
</apply>
<apply>
<plus/>
<apply>
<times/>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 4.0 </cn>
<ci> beta_LCC </ci>
</apply>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
<ci> fc_LCC </ci>
<apply>
<divide/>
<ci> omega_LCC </ci>
<apply>
<power/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<cn cellml:units="dimensionless"> 4.0 </cn>
</apply>
</apply>
</apply>
<ci> P11 </ci>
</apply>
<ci> PHOSPH11 </ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>P12</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci> fc_LCC </ci>
<ci> P11 </ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci> gc_LCC </ci>
<ci> P12 </ci>
</apply>
<ci> PHOSPH12 </ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>P13</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<ci> beta_LCC </ci>
<ci> P14 </ci>
</apply>
<apply>
<times/>
<ci> omega_LCC </ci>
<ci> P19 </ci>
</apply>
<ci> PHOSPH1 </ci>
</apply>
<apply>
<times/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 4.0 </cn>
<ci> alpha_LCC </ci>
</apply>
<ci> gamma_LCC </ci>
</apply>
<ci> P13 </ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>P14</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 4.0 </cn>
<ci> alpha_LCC </ci>
<ci> P13 </ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> beta_LCC </ci>
<ci> P15 </ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<ci> omega_LCC </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
<ci> P20 </ci>
</apply>
<ci> PHOSPH2 </ci>
</apply>
<apply>
<times/>
<apply>
<plus/>
<ci> beta_LCC </ci>
<apply>
<times/>
<cn cellml:units="dimensionless"> 3.0 </cn>
<ci> alpha_LCC </ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> gamma_LCC </ci>
</apply>
</apply>
<ci> P14 </ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>P15</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 3.0 </cn>
<ci> alpha_LCC </ci>
<ci> P14 </ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 3.0 </cn>
<ci> beta_LCC </ci>
<ci> P16 </ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<ci> omega_LCC </ci>
<apply>
<power/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
</apply>
<ci> P21 </ci>
</apply>
<ci> PHOSPH3 </ci>
</apply>
<apply>
<times/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> beta_LCC </ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> alpha_LCC </ci>
</apply>
<apply>
<times/>
<apply>
<power/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
<ci> gamma_LCC </ci>
</apply>
</apply>
<ci> P15 </ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>P16</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> alpha_LCC </ci>
<ci> P15 </ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 4.0 </cn>
<ci> beta_LCC </ci>
<ci> P17 </ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<ci> omega_LCC </ci>
<apply>
<power/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<cn cellml:units="dimensionless"> 3.0 </cn>
</apply>
</apply>
<ci> P22 </ci>
</apply>
<ci> PHOSPH4 </ci>
</apply>
<apply>
<times/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 3.0 </cn>
<ci> beta_LCC </ci>
</apply>
<ci> alpha_LCC </ci>
<apply>
<times/>
<apply>
<power/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<cn cellml:units="dimensionless"> 3.0 </cn>
</apply>
<ci> gamma_LCC </ci>
</apply>
</apply>
<ci> P16 </ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>P17</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<ci> alpha_LCC </ci>
<ci> P16 </ci>
</apply>
<apply>
<times/>
<ci> gn_LCC </ci>
<ci> P18 </ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<ci> omega_LCC </ci>
<apply>
<power/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<cn cellml:units="dimensionless"> 4.0 </cn>
</apply>
</apply>
<ci> P23 </ci>
</apply>
<ci> PHOSPH5 </ci>
</apply>
<apply>
<times/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 4.0 </cn>
<ci> beta_LCC </ci>
</apply>
<ci> fn_LCC </ci>
<apply>
<times/>
<apply>
<power/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<cn cellml:units="dimensionless"> 4.0 </cn>
</apply>
<ci> gamma_LCC </ci>
</apply>
</apply>
<ci> P17 </ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>P18</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<ci> fc_LCC </ci>
<ci> P17 </ci>
</apply>
<ci> PHOSPH6 </ci>
</apply>
<apply>
<times/>
<ci> gc_LCC </ci>
<ci> P18 </ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>P19</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<apply>
<divide/>
<ci> beta_LCC </ci>
<cn cellml:units="dimensionless"> 4.0 </cn>
</apply>
<ci> P20 </ci>
</apply>
<apply>
<times/>
<ci> gamma_LCC </ci>
<ci> P13 </ci>
</apply>
<ci> PHOSPH7 </ci>
</apply>
<apply>
<times/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 4.0 </cn>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> alpha_LCC </ci>
</apply>
<ci> omega_LCC </ci>
</apply>
<ci> P19 </ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>P20</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 4.0 </cn>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> alpha_LCC </ci>
<ci> P19 </ci>
</apply>
<apply>
<times/>
<ci> beta_LCC </ci>
<ci> P21 </ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> gamma_LCC </ci>
<ci> P14 </ci>
</apply>
<ci> PHOSPH8 </ci>
</apply>
<apply>
<times/>
<apply>
<plus/>
<apply>
<divide/>
<ci> beta_LCC </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 3.0 </cn>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> alpha_LCC </ci>
</apply>
<apply>
<divide/>
<ci> omega_LCC </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
</apply>
<ci> P20 </ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>P21</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 3.0 </cn>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> alpha_LCC </ci>
<ci> P18 </ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 3.0 </cn>
<ci> beta_LCC </ci>
</apply>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
<ci> P20 </ci>
</apply>
<apply>
<times/>
<apply>
<power/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
<ci> gamma_LCC </ci>
<ci> P15 </ci>
</apply>
<ci> PHOSPH9 </ci>
</apply>
<apply>
<times/>
<apply>
<plus/>
<ci> beta_LCC </ci>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> alpha_LCC </ci>
</apply>
<apply>
<divide/>
<ci> omega_LCC </ci>
<apply>
<power/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
</apply>
</apply>
<ci> P21 </ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>P22</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> alpha_LCC </ci>
<ci> P17 </ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> beta_LCC </ci>
<ci> P19 </ci>
</apply>
<apply>
<times/>
<apply>
<power/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<cn cellml:units="dimensionless"> 3.0 </cn>
</apply>
<ci> gamma_LCC </ci>
<ci> P16 </ci>
</apply>
<ci> PHOSPH10 </ci>
</apply>
<apply>
<times/>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 3.0 </cn>
<ci> beta_LCC </ci>
</apply>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> alpha_LCC </ci>
</apply>
<apply>
<divide/>
<ci> omega_LCC </ci>
<apply>
<power/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<cn cellml:units="dimensionless"> 3.0 </cn>
</apply>
</apply>
</apply>
<ci> P22 </ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>P23</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> alpha_LCC </ci>
<ci> P16 </ci>
</apply>
<apply>
<times/>
<ci> gc_LCC </ci>
<ci> P24 </ci>
</apply>
<apply>
<times/>
<apply>
<power/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<cn cellml:units="dimensionless"> 4.0 </cn>
</apply>
<ci> gamma_LCC </ci>
<ci> P17 </ci>
</apply>
<ci> PHOSPH11 </ci>
</apply>
<apply>
<times/>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 4.0 </cn>
<ci> beta_LCC </ci>
</apply>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
<ci> fc_LCC </ci>
<apply>
<divide/>
<ci> omega_LCC </ci>
<apply>
<power/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<cn cellml:units="dimensionless"> 4.0 </cn>
</apply>
</apply>
</apply>
<ci> P23 </ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>P24</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<ci> fc_LCC </ci>
<ci> P23 </ci>
</apply>
<ci> PHOSPH12 </ci>
</apply>
<apply>
<times/>
<ci> gc_LCC </ci>
<ci> P24 </ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci> PHOSPH1 </ci>
<apply>
<minus/>
<apply>
<divide/>
<apply>
<times/>
<ci> epsilon </ci>
<ci> kPKA_LCC </ci>
<ci> PKACII </ci>
<ci> P1 </ci>
</apply>
<apply>
<plus/>
<ci> KmPKA_LCC </ci>
<apply>
<times/>
<ci> epsilon </ci>
<ci> LCCtot </ci>
<ci> P1 </ci>
</apply>
</apply>
</apply>
<apply>
<divide/>
<apply>
<times/>
<ci> epsilon </ci>
<ci> kPP2A_LCC </ci>
<ci> PP2A_LCC </ci>
<ci> P13 </ci>
</apply>
<apply>
<plus/>
<ci> KmPP2A_LCC </ci>
<apply>
<times/>
<ci> epsilon </ci>
<ci> LCCtot </ci>
<ci> P13 </ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci> PHOSPH2 </ci>
<apply>
<minus/>
<apply>
<divide/>
<apply>
<times/>
<ci> epsilon </ci>
<ci> kPKA_LCC </ci>
<ci> PKACII </ci>
<ci> P2 </ci>
</apply>
<apply>
<plus/>
<ci> KmPKA_LCC </ci>
<apply>
<times/>
<ci> epsilon </ci>
<ci> LCCtot </ci>
<ci> P2 </ci>
</apply>
</apply>
</apply>
<apply>
<divide/>
<apply>
<times/>
<ci> epsilon </ci>
<ci> kPP2A_LCC </ci>
<ci> PP2A_LCC </ci>
<ci> P14 </ci>
</apply>
<apply>
<plus/>
<ci> KmPP2A_LCC </ci>
<apply>
<times/>
<ci> epsilon </ci>
<ci> LCCtot </ci>
<ci> P14 </ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci> PHOSPH3 </ci>
<apply>
<minus/>
<apply>
<divide/>
<apply>
<times/>
<ci> epsilon </ci>
<ci> kPKA_LCC </ci>
<ci> PKACII </ci>
<ci> P3 </ci>
</apply>
<apply>
<plus/>
<ci> KmPKA_LCC </ci>
<apply>
<times/>
<ci> epsilon </ci>
<ci> LCCtot </ci>
<ci> P3 </ci>
</apply>
</apply>
</apply>
<apply>
<divide/>
<apply>
<times/>
<ci> epsilon </ci>
<ci> kPP2A_LCC </ci>
<ci> PP2A_LCC </ci>
<ci> P15 </ci>
</apply>
<apply>
<plus/>
<ci> KmPP2A_LCC </ci>
<apply>
<times/>
<ci> epsilon </ci>
<ci> LCCtot </ci>
<ci> P15 </ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci> PHOSPH4 </ci>
<apply>
<minus/>
<apply>
<divide/>
<apply>
<times/>
<ci> epsilon </ci>
<ci> kPKA_LCC </ci>
<ci> PKACII </ci>
<ci> P4 </ci>
</apply>
<apply>
<plus/>
<ci> KmPKA_LCC </ci>
<apply>
<times/>
<ci> epsilon </ci>
<ci> LCCtot </ci>
<ci> P4 </ci>
</apply>
</apply>
</apply>
<apply>
<divide/>
<apply>
<times/>
<ci> epsilon </ci>
<ci> kPP2A_LCC </ci>
<ci> PP2A_LCC </ci>
<ci> P16 </ci>
</apply>
<apply>
<plus/>
<ci> KmPP2A_LCC </ci>
<apply>
<times/>
<ci> epsilon </ci>
<ci> LCCtot </ci>
<ci> P16 </ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci> PHOSPH5 </ci>
<apply>
<minus/>
<apply>
<divide/>
<apply>
<times/>
<ci> epsilon </ci>
<ci> kPKA_LCC </ci>
<ci> PKACII </ci>
<ci> P5 </ci>
</apply>
<apply>
<plus/>
<ci> KmPKA_LCC </ci>
<apply>
<times/>
<ci> epsilon </ci>
<ci> LCCtot </ci>
<ci> P5 </ci>
</apply>
</apply>
</apply>
<apply>
<divide/>
<apply>
<times/>
<ci> epsilon </ci>
<ci> kPP2A_LCC </ci>
<ci> PP2A_LCC </ci>
<ci> P17 </ci>
</apply>
<apply>
<plus/>
<ci> KmPP2A_LCC </ci>
<apply>
<times/>
<ci> epsilon </ci>
<ci> LCCtot </ci>
<ci> P17 </ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci> PHOSPH6 </ci>
<apply>
<minus/>
<apply>
<divide/>
<apply>
<times/>
<ci> epsilon </ci>
<ci> kPKA_LCC </ci>
<ci> PKACII </ci>
<ci> P6 </ci>
</apply>
<apply>
<plus/>
<ci> KmPKA_LCC </ci>
<apply>
<times/>
<ci> epsilon </ci>
<ci> LCCtot </ci>
<ci> P6 </ci>
</apply>
</apply>
</apply>
<apply>
<divide/>
<apply>
<times/>
<ci> epsilon </ci>
<ci> kPP2A_LCC </ci>
<ci> PP2A_LCC </ci>
<ci> P18 </ci>
</apply>
<apply>
<plus/>
<ci> KmPP2A_LCC </ci>
<apply>
<times/>
<ci> epsilon </ci>
<ci> LCCtot </ci>
<ci> P18 </ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci> PHOSPH7 </ci>
<apply>
<minus/>
<apply>
<divide/>
<apply>
<times/>
<ci> epsilon </ci>
<ci> kPKA_LCC </ci>
<ci> PKACII </ci>
<ci> P7 </ci>
</apply>
<apply>
<plus/>
<ci> KmPKA_LCC </ci>
<apply>
<times/>
<ci> epsilon </ci>
<ci> LCCtot </ci>
<ci> P7 </ci>
</apply>
</apply>
</apply>
<apply>
<divide/>
<apply>
<times/>
<ci> epsilon </ci>
<ci> kPP2A_LCC </ci>
<ci> PP2A_LCC </ci>
<ci> P19 </ci>
</apply>
<apply>
<plus/>
<ci> KmPP2A_LCC </ci>
<apply>
<times/>
<ci> epsilon </ci>
<ci> LCCtot </ci>
<ci> P19 </ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci> PHOSPH8 </ci>
<apply>
<minus/>
<apply>
<divide/>
<apply>
<times/>
<ci> epsilon </ci>
<ci> kPKA_LCC </ci>
<ci> PKACII </ci>
<ci> P8 </ci>
</apply>
<apply>
<plus/>
<ci> KmPKA_LCC </ci>
<apply>
<times/>
<ci> epsilon </ci>
<ci> LCCtot </ci>
<ci> P8 </ci>
</apply>
</apply>
</apply>
<apply>
<divide/>
<apply>
<times/>
<ci> epsilon </ci>
<ci> kPP2A_LCC </ci>
<ci> PP2A_LCC </ci>
<ci> P20 </ci>
</apply>
<apply>
<plus/>
<ci> KmPP2A_LCC </ci>
<apply>
<times/>
<ci> epsilon </ci>
<ci> LCCtot </ci>
<ci> P20 </ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci> PHOSPH9 </ci>
<apply>
<minus/>
<apply>
<divide/>
<apply>
<times/>
<ci> epsilon </ci>
<ci> kPKA_LCC </ci>
<ci> PKACII </ci>
<ci> P9 </ci>
</apply>
<apply>
<plus/>
<ci> KmPKA_LCC </ci>
<apply>
<times/>
<ci> epsilon </ci>
<ci> LCCtot </ci>
<ci> P9 </ci>
</apply>
</apply>
</apply>
<apply>
<divide/>
<apply>
<times/>
<ci> epsilon </ci>
<ci> kPP2A_LCC </ci>
<ci> PP2A_LCC </ci>
<ci> P21 </ci>
</apply>
<apply>
<plus/>
<ci> KmPP2A_LCC </ci>
<apply>
<times/>
<ci> epsilon </ci>
<ci> LCCtot </ci>
<ci> P21 </ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci> PHOSPH10 </ci>
<apply>
<minus/>
<apply>
<divide/>
<apply>
<times/>
<ci> epsilon </ci>
<ci> kPKA_LCC </ci>
<ci> PKACII </ci>
<ci> P10 </ci>
</apply>
<apply>
<plus/>
<ci> KmPKA_LCC </ci>
<apply>
<times/>
<ci> epsilon </ci>
<ci> LCCtot </ci>
<ci> P10 </ci>
</apply>
</apply>
</apply>
<apply>
<divide/>
<apply>
<times/>
<ci> epsilon </ci>
<ci> kPP2A_LCC </ci>
<ci> PP2A_LCC </ci>
<ci> P22 </ci>
</apply>
<apply>
<plus/>
<ci> KmPP2A_LCC </ci>
<apply>
<times/>
<ci> epsilon </ci>
<ci> LCCtot </ci>
<ci> P22 </ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci> PHOSPH11 </ci>
<apply>
<minus/>
<apply>
<divide/>
<apply>
<times/>
<ci> epsilon </ci>
<ci> kPKA_LCC </ci>
<ci> PKACII </ci>
<ci> P11 </ci>
</apply>
<apply>
<plus/>
<ci> KmPKA_LCC </ci>
<apply>
<times/>
<ci> epsilon </ci>
<ci> LCCtot </ci>
<ci> P11 </ci>
</apply>
</apply>
</apply>
<apply>
<divide/>
<apply>
<times/>
<ci> epsilon </ci>
<ci> kPP2A_LCC </ci>
<ci> PP2A_LCC </ci>
<ci> P23 </ci>
</apply>
<apply>
<plus/>
<ci> KmPP2A_LCC </ci>
<apply>
<times/>
<ci> epsilon </ci>
<ci> LCCtot </ci>
<ci> P23 </ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci> PHOSPH12 </ci>
<apply>
<minus/>
<apply>
<divide/>
<apply>
<times/>
<ci> epsilon </ci>
<ci> kPKA_LCC </ci>
<ci> PKACII </ci>
<ci> P12 </ci>
</apply>
<apply>
<plus/>
<ci> KmPKA_LCC </ci>
<apply>
<times/>
<ci> epsilon </ci>
<ci> LCCtot </ci>
<ci> P12 </ci>
</apply>
</apply>
</apply>
<apply>
<divide/>
<apply>
<times/>
<ci> epsilon </ci>
<ci> kPP2A_LCC </ci>
<ci> PP2A_LCC </ci>
<ci> P24 </ci>
</apply>
<apply>
<plus/>
<ci> KmPP2A_LCC </ci>
<apply>
<times/>
<ci> epsilon </ci>
<ci> LCCtot </ci>
<ci> P24 </ci>
</apply>
</apply>
</apply>
</apply>
</apply>
-->
</math>
</component>
<!-- EP component of the model from the Saucerman et al. 2003 paper -->
<component name="membrane">
<variable units="millivolt" public_interface="out" name="V" initial_value="-85.66"/>
<variable units="joule_per_kilomole_kelvin" public_interface="out" name="R" initial_value="8314"/>
<variable units="kelvin" public_interface="out" name="T" initial_value="310"/>
<variable units="coulomb_per_mole" public_interface="out" name="F" initial_value="96485"/>
<variable units="second" public_interface="in" name="time"/>
<variable units="microA_per_microF" public_interface="in" name="I_Na"/>
<variable units="microA_per_microF" public_interface="in" name="I_Ca"/>
<variable units="microA_per_microF" public_interface="in" name="I_CaK"/>
<variable units="microA_per_microF" public_interface="in" name="I_NaCa"/>
<variable units="microA_per_microF" public_interface="in" name="I_Ki"/>
<variable units="microA_per_microF" public_interface="in" name="I_Kp"/>
<variable units="microA_per_microF" public_interface="in" name="I_p_Ca"/>
<variable units="microA_per_microF" public_interface="in" name="I_Na_b"/>
<variable units="microA_per_microF" public_interface="in" name="I_Ca_b"/>
<variable units="microA_per_microF" public_interface="in" name="I_NaK"/>
<variable units="microA_per_microF" public_interface="in" name="I_to"/>
<variable units="microA_per_microF" public_interface="in" name="I_ss"/>
<variable units="microA_per_microF" name="I_app"/>
<variable units="millivolt" name="R_clamp" initial_value="0.2"/>
<variable units="millivolt" name="V_clamp"/>
<variable units="millivolt" name="V_test" initial_value="-10.0"/>
<variable units="millivolt" name="V_hold" initial_value="-40.0"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="7">
<apply>
<eq/>
<ci> I_app </ci>
<apply>
<divide/>
<apply>
<minus/>
<ci> V_clamp </ci>
<ci> V </ci>
</apply>
<ci> R_clamp </ci>
</apply>
</apply>
<apply>
<eq/>
<ci> V_clamp </ci>
<piecewise>
<piece>
<ci> V_test </ci>
<apply>
<and/>
<apply>
<gt/>
<ci>time</ci>
<cn cellml:units="millisecond">59.1</cn>
</apply>
<apply>
<leq/>
<ci>time</ci>
<cn cellml:units="millisecond">59.5</cn>
</apply>
</apply>
</piece>
<otherwise>
<ci> V_hold </ci>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>V</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">-1E3</cn>
<apply>
<plus/>
<ci>I_Na</ci>
<ci>I_Na_b</ci>
<ci>I_NaCa</ci>
<ci>I_NaK</ci>
<ci>I_to</ci>
<ci>I_ss</ci>
<ci>I_Ki</ci>
<ci>I_Kp</ci>
<ci>I_CaK</ci>
<ci>I_Ca</ci>
<ci>I_Ca_b</ci>
<ci>I_p_Ca</ci>
<apply>
<minus/>
<ci>I_app</ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<!-- the following descriptions of the transient outward and steady state potassium currents have been taken from the Pandit et al. 2001 model. -->
<component name="Ca_independent_transient_outward_K_current">
<variable units="microA_per_microF" public_interface="out" name="I_to"/>
<variable units="milliS_per_microF" name="g_to" initial_value="0.35"/>
<variable units="dimensionless" name="a" initial_value="0.886"/>
<variable units="dimensionless" name="b" initial_value="0.114"/>
<variable units="second" public_interface="in" private_interface="out" name="time"/>
<variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
<variable units="dimensionless" private_interface="in" name="r"/>
<variable units="dimensionless" private_interface="in" name="s"/>
<variable units="dimensionless" private_interface="in" name="s_slow"/>
<variable units="millivolt" public_interface="in" name="E_K"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="8">
<apply>
<eq/>
<ci>I_to</ci>
<apply>
<times/>
<ci>g_to</ci>
<ci>r</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>a</ci>
<ci>s</ci>
</apply>
<apply>
<times/>
<ci>b</ci>
<ci>s_slow</ci>
</apply>
</apply>
<apply>
<minus/>
<ci>V</ci>
<ci>E_K</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="Ca_independent_transient_outward_K_current_r_gate">
<variable units="dimensionless" public_interface="out" name="r" initial_value="1.4E-3"/>
<variable units="second" name="tau_r"/>
<variable units="dimensionless" name="r_infinity"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="second" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="9">
<apply>
<eq/>
<ci>r_infinity</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">10.6</cn>
</apply>
<apply>
<minus/>
<cn cellml:units="millivolt">11.42</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>tau_r</ci>
<apply>
<divide/>
<cn cellml:units="second">1</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">45.16</cn>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="dimensionless">0.03577</cn>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">50</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">98.9</cn>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">0.1</cn>
</apply>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">38</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>r</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>r_infinity</ci>
<ci>r</ci>
</apply>
<ci>tau_r</ci>
</apply>
</apply>
</math>
</component>
<component name="Ca_independent_transient_outward_K_current_s_gate">
<variable units="dimensionless" public_interface="out" name="s" initial_value="1.0"/>
<variable units="second" name="tau_s"/>
<variable units="dimensionless" name="s_infinity"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="second" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="10">
<apply>
<eq/>
<ci>s_infinity</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">45.3</cn>
</apply>
<cn cellml:units="millivolt">6.8841</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>tau_s</ci>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="second">0.35</cn>
<apply>
<exp/>
<apply>
<minus/>
<apply>
<power/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">70</cn>
</apply>
<cn cellml:units="millivolt">15</cn>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
<cn cellml:units="second">0.035</cn>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>s</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>s_infinity</ci>
<ci>s</ci>
</apply>
<ci>tau_s</ci>
</apply>
</apply>
</math>
</component>
<component name="Ca_independent_transient_outward_K_current_s_slow_gate">
<variable units="dimensionless" public_interface="out" name="s_slow" initial_value="0.613"/>
<variable units="second" name="tau_s_slow"/>
<variable units="dimensionless" name="s_slow_infinity"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="second" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="11">
<apply>
<eq/>
<ci>s_slow_infinity</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">45.3</cn>
</apply>
<cn cellml:units="millivolt">6.8841</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>tau_s_slow</ci>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="second">3.7</cn>
<apply>
<exp/>
<apply>
<minus/>
<apply>
<power/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">70</cn>
</apply>
<cn cellml:units="millivolt">30</cn>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
<cn cellml:units="second">0.035</cn>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>s_slow</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>s_slow_infinity</ci>
<ci>s_slow</ci>
</apply>
<ci>tau_s_slow</ci>
</apply>
</apply>
</math>
</component>
<component name="steady_state_outward_K_current">
<variable units="microA_per_microF" public_interface="out" name="I_ss"/>
<variable units="milliS_per_microF" name="g_ss" initial_value="0.07"/>
<variable units="second" public_interface="in" private_interface="out" name="time"/>
<variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
<variable units="millivolt" public_interface="in" name="E_K"/>
<variable units="dimensionless" private_interface="in" name="r_ss"/>
<variable units="dimensionless" private_interface="in" name="s_ss"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="12">
<apply>
<eq/>
<ci>I_ss</ci>
<apply>
<times/>
<ci>g_ss</ci>
<ci>r_ss</ci>
<ci>s_ss</ci>
<apply>
<minus/>
<ci>V</ci>
<ci>E_K</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="steady_state_outward_K_current_r_ss_gate">
<variable units="dimensionless" public_interface="out" name="r_ss" initial_value="198E-3"/>
<variable units="second" name="tau_r_ss"/>
<variable units="dimensionless" name="r_ss_infinity"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="second" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>r_ss_infinity</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">11.5</cn>
</apply>
<apply>
<minus/>
<cn cellml:units="millivolt">11.82</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>tau_r_ss</ci>
<apply>
<divide/>
<cn cellml:units="second">10</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">45.16</cn>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="dimensionless">0.03577</cn>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">50</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">98.9</cn>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">0.1</cn>
</apply>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">38</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>r_ss</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>r_ss_infinity</ci>
<ci>r_ss</ci>
</apply>
<ci>tau_r_ss</ci>
</apply>
</apply>
</math>
</component>
<component name="steady_state_outward_K_current_s_ss_gate">
<variable units="dimensionless" public_interface="out" name="s_ss" initial_value="0.43"/>
<variable units="second" name="tau_s_ss"/>
<variable units="dimensionless" name="s_ss_infinity"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="second" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="13">
<apply>
<eq/>
<ci>s_ss_infinity</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">87.5</cn>
</apply>
<cn cellml:units="millivolt">10.3</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>tau_s_ss</ci>
<cn cellml:units="second">2.1</cn>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>s_ss</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci>s_ss_infinity</ci>
<ci>s_ss</ci>
</apply>
<ci>tau_s_ss</ci>
</apply>
</apply>
</math>
</component>
<!-- the following components have been taken from an adapted version of the Luo-Rudy model published by Puglisi and Bers 2001. -->
<component name="fast_sodium_current">
<variable units="microA_per_microF" public_interface="out" name="I_Na"/>
<variable units="millivolt" public_interface="out" name="E_Na"/>
<variable units="milliS_per_microF" name="g_Na" initial_value="8.0"/>
<variable units="second" public_interface="in" private_interface="out" name="time"/>
<variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
<variable units="joule_per_kilomole_kelvin" public_interface="in" name="R"/>
<variable units="coulomb_per_mole" public_interface="in" name="F"/>
<variable units="kelvin" public_interface="in" name="T"/>
<variable units="millimolar" public_interface="in" name="Nai"/>
<variable units="millimolar" public_interface="in" name="Nao"/>
<variable units="dimensionless" private_interface="in" name="m"/>
<variable units="dimensionless" private_interface="in" name="h"/>
<variable units="dimensionless" private_interface="in" name="j"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="14">
<apply>
<eq/>
<ci>E_Na</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
<ci>F</ci>
</apply>
<apply>
<log/>
<apply>
<divide/>
<ci>Nao</ci>
<ci>Nai</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>I_Na</ci>
<apply>
<times/>
<ci>g_Na</ci>
<apply>
<power/>
<ci>m</ci>
<cn cellml:units="dimensionless">3</cn>
</apply>
<ci>h</ci>
<ci>j</ci>
<apply>
<minus/>
<ci>V</ci>
<ci>E_Na</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="fast_sodium_current_m_gate">
<variable units="dimensionless" public_interface="out" name="m" initial_value="1.4E-3"/>
<variable units="first_order_rate_constant" name="alpha_m"/>
<variable units="first_order_rate_constant" name="beta_m"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="second" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="15">
<apply>
<eq/>
<ci>alpha_m</ci>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="per_millivolt_second">0.32</cn>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">47.13</cn>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_millivolt">0.1</cn>
</apply>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">47.13</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_m</ci>
<apply>
<times/>
<cn cellml:units="first_order_rate_constant">0.08</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>V</ci>
</apply>
<cn cellml:units="millivolt">11</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>m</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1E3</cn>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_m</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>m</ci>
</apply>
</apply>
<apply>
<times/>
<ci>beta_m</ci>
<ci>m</ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="fast_sodium_current_h_gate">
<variable units="dimensionless" public_interface="out" name="h" initial_value="0.99"/>
<variable units="first_order_rate_constant" name="alpha_h"/>
<variable units="first_order_rate_constant" name="beta_h"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="second" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="16">
<apply>
<eq/>
<ci>alpha_h</ci>
<piecewise>
<piece>
<apply>
<times/>
<cn cellml:units="first_order_rate_constant">0.135</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<cn cellml:units="millivolt">80</cn>
<ci>V</ci>
</apply>
<apply>
<minus/>
<cn cellml:units="millivolt">6.8</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<lt/>
<ci>V</ci>
<apply>
<minus/>
<cn cellml:units="millivolt">40</cn>
</apply>
</apply>
</piece>
<otherwise>
<cn cellml:units="first_order_rate_constant">0</cn>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>beta_h</ci>
<piecewise>
<piece>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="first_order_rate_constant">3.56</cn>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="per_millivolt">0.079</cn>
<ci>V</ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="first_order_rate_constant">3.1E5</cn>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="per_millivolt">0.35</cn>
<ci>V</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<lt/>
<ci>V</ci>
<apply>
<minus/>
<cn cellml:units="millivolt">40</cn>
</apply>
</apply>
</piece>
<otherwise>
<apply>
<divide/>
<cn cellml:units="first_order_rate_constant">1.0</cn>
<apply>
<times/>
<cn cellml:units="dimensionless">0.13</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">10.66</cn>
</apply>
<apply>
<minus/>
<cn cellml:units="millivolt">11.1</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>h</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1E3</cn>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_h</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>h</ci>
</apply>
</apply>
<apply>
<times/>
<ci>beta_h</ci>
<ci>h</ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="fast_sodium_current_j_gate">
<variable units="dimensionless" public_interface="out" name="j" initial_value="0.99"/>
<variable units="first_order_rate_constant" name="alpha_j"/>
<variable units="first_order_rate_constant" name="beta_j"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="second" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="17">
<apply>
<eq/>
<ci>alpha_j</ci>
<piecewise>
<piece>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="per_millivolt_second">1.2714E5</cn>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="per_millivolt">0.2444</cn>
<ci>V</ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="per_millivolt_second">3.474e-5</cn>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_millivolt">0.04391</cn>
</apply>
<ci>V</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">37.78</cn>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="per_millivolt">0.311</cn>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">79.23</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<lt/>
<ci>V</ci>
<apply>
<minus/>
<cn cellml:units="millivolt">40</cn>
</apply>
</apply>
</piece>
<otherwise>
<cn cellml:units="first_order_rate_constant">0</cn>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>beta_j</ci>
<piecewise>
<piece>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="first_order_rate_constant">0.1212</cn>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_millivolt">0.01052</cn>
</apply>
<ci>V</ci>
</apply>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_millivolt">0.1378</cn>
</apply>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">40.14</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<lt/>
<ci>V</ci>
<apply>
<minus/>
<cn cellml:units="millivolt">40</cn>
</apply>
</apply>
</piece>
<otherwise>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="first_order_rate_constant">0.3</cn>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_millivolt">2.535e-7</cn>
</apply>
<ci>V</ci>
</apply>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_millivolt">0.1</cn>
</apply>
<apply>
<plus/>
<ci>V</ci>
<cn cellml:units="millivolt">32</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>j</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1E3</cn>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_j</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>j</ci>
</apply>
</apply>
<apply>
<times/>
<ci>beta_j</ci>
<ci>j</ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<!-- the following description of the L-type calcium channel has been adapted from the Jafri, Rice, and Winslow 1998 model. -->
<component name="L_type_Ca_channel">
<variable units="microA_per_microF" public_interface="out" name="I_Ca"/>
<variable units="microA_per_microF" public_interface="out" name="I_CaK"/>
<variable units="first_order_rate_constant" private_interface="out" name="alpha"/>
<variable units="first_order_rate_constant" private_interface="out" name="beta"/>
<variable units="first_order_rate_constant" private_interface="out" name="gamma"/>
<variable units="cm_per_second" name="P_Ca" initial_value="1.7469E-8"/>
<variable units="cm_per_second" name="P_K" initial_value="3.234E-11"/>
<variable units="microA_per_microF" name="I_Ca05" initial_value="-0.458"/>
<variable units="microA_per_microF" name="ibarca"/>
<variable units="microA_per_microF" name="ibark"/>
<variable units="second_order_rate_constant" name="gammao" initial_value="5187.5"/>
<variable units="dimensionless" name="Nlcc" initial_value="3E5"/>
<variable units="dimensionless" name="favail"/>
<variable units="second" public_interface="in" private_interface="out" name="time"/>
<variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
<variable units="dimensionless" public_interface="in" name="fracLCCbp"/>
<variable units="dimensionless" public_interface="in" name="fracLCCbpo"/>
<variable units="dimensionless" public_interface="in" private_interface="out" name="fracLCCap"/>
<variable units="dimensionless" public_interface="in" private_interface="out" name="fracLCCapo"/>
<variable units="millimolar" public_interface="in" name="Cao"/>
<variable units="millimolar" public_interface="in" name="Cai"/>
<variable units="millimolar" public_interface="in" name="Ko"/>
<variable units="millimolar" public_interface="in" name="Ki"/>
<variable units="joule_per_kilomole_kelvin" public_interface="in" name="R"/>
<variable units="kelvin" public_interface="in" name="T"/>
<variable units="coulomb_per_mole" public_interface="in" name="F"/>
<variable units="dimensionless" private_interface="in" name="y"/>
<variable units="dimensionless" private_interface="in" name="v"/>
<variable units="dimensionless" private_interface="in" name="z"/>
<variable units="dimensionless" private_interface="in" name="x"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="18">
<apply id="I_Ca_calculation">
<eq/>
<ci> I_Ca </ci>
<apply>
<times/>
<ci> ibarca </ci>
<ci> Nlcc </ci>
<ci> favail </ci>
<apply>
<power/>
<ci> v </ci>
<cn cellml:units="dimensionless"> 4.0 </cn>
</apply>
<ci> x </ci>
<ci> y </ci>
<ci> z </ci>
</apply>
</apply>
<apply id="I_CaK_calculation">
<eq/>
<ci> I_CaK </ci>
<apply>
<times/>
<apply>
<divide/>
<ci> ibark </ci>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<divide/>
<ci> I_Ca </ci>
<ci> I_Ca05 </ci>
</apply>
</apply>
</apply>
<ci> Nlcc </ci>
<ci> favail </ci>
<apply>
<power/>
<ci> v </ci>
<cn cellml:units="dimensionless"> 4.0 </cn>
</apply>
<ci> x </ci>
<ci> y </ci>
<ci> z </ci>
</apply>
</apply>
<apply id="ibarca_calculation">
<eq/>
<ci> ibarca </ci>
<apply>
<times/>
<ci> P_Ca </ci>
<cn cellml:units="dimensionless"> 4.0 </cn>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> V </ci>
<ci> F </ci>
</apply>
<apply>
<times/>
<ci> R </ci>
<ci> T </ci>
</apply>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 0.001 </cn>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> V </ci>
<apply>
<divide/>
<ci> F </ci>
<apply>
<times/>
<ci> R </ci>
<ci> T </ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 0.341 </cn>
<ci> Cao </ci>
</apply>
</apply>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> V </ci>
<apply>
<divide/>
<ci> F </ci>
<apply>
<times/>
<ci> R </ci>
<ci> T </ci>
</apply>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless"> 1.0 </cn>
</apply>
</apply>
</apply>
</apply>
<apply id="ibark_calculation">
<eq/>
<ci> ibark </ci>
<apply>
<divide/>
<apply>
<times/>
<ci> P_K </ci>
<apply>
<divide/>
<apply>
<times/>
<ci> V </ci>
<ci> F </ci>
</apply>
<apply>
<times/>
<ci> R </ci>
<ci> T </ci>
</apply>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci> Ki </ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci> V </ci>
<ci> F </ci>
</apply>
<apply>
<times/>
<ci> R </ci>
<ci> T </ci>
</apply>
</apply>
</apply>
</apply>
<ci> Ko </ci>
</apply>
</apply>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci> V </ci>
<ci> F </ci>
</apply>
<apply>
<times/>
<ci> R </ci>
<ci> T </ci>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless"> 1.0 </cn>
</apply>
</apply>
</apply>
<apply id="alpha_calculation">
<eq/>
<ci> alpha </ci>
<apply>
<times/>
<cn cellml:units="first_order_rate_constant"> 400.0 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 2.0 </cn>
</apply>
<cn cellml:units="millivolt"> 10.0 </cn>
</apply>
</apply>
</apply>
</apply>
<apply id="beta_calculation">
<eq/>
<ci> beta </ci>
<apply>
<times/>
<cn cellml:units="first_order_rate_constant"> 50.0 </cn>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="dimensionless"> -1.0 </cn>
<apply>
<divide/>
<apply>
<plus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 2.0 </cn>
</apply>
<cn cellml:units="millivolt"> 13.0 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply id="gamma_calculation">
<eq/>
<ci> gamma </ci>
<apply>
<times/>
<ci> gammao </ci>
<ci> Cai </ci>
</apply>
</apply>
<apply id="favail_calculation">
<eq/>
<ci> favail </ci>
<apply>
<times/>
<cn cellml:units="dimensionless"> 0.5 </cn>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 0.4 </cn>
<ci> fracLCCbp </ci>
</apply>
<apply>
<plus/>
<ci> fracLCCbpo </ci>
<cn cellml:units="dimensionless"> 0.60 </cn>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="L_type_Ca_channel_v_gate">
<variable units="dimensionless" public_interface="out" name="v" initial_value="0.0"/>
<variable units="first_order_rate_constant" public_interface="in" name="alpha"/>
<variable units="first_order_rate_constant" public_interface="in" name="beta"/>
<variable units="second" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="19">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> v </ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci> alpha </ci>
<apply>
<minus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<ci> v </ci>
</apply>
</apply>
<apply>
<times/>
<ci> beta </ci>
<ci> v </ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="L_type_Ca_channel_w_gate">
<variable units="dimensionless" public_interface="out" name="w" initial_value="0.0"/>
<variable units="first_order_rate_constant" public_interface="in" name="alpha"/>
<variable units="first_order_rate_constant" public_interface="in" name="beta"/>
<variable units="second" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="20">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> w </ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> alpha </ci>
<apply>
<minus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<ci> w </ci>
</apply>
</apply>
<apply>
<divide/>
<ci> beta </ci>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> w </ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="L_type_Ca_channel_x_gate">
<variable units="dimensionless" public_interface="out" name="x" initial_value="0.13"/>
<variable units="first_order_rate_constant" public_interface="out" name="g" initial_value="2E3"/>
<variable units="first_order_rate_constant" public_interface="out" name="flcc"/>
<variable units="first_order_rate_constant" name="f" initial_value="300.0"/>
<variable units="dimensionless" public_interface="in" name="fracLCCap"/>
<variable units="dimensionless" public_interface="in" name="fracLCCapo"/>
<variable units="second" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="21">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> x </ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci> flcc </ci>
<apply>
<minus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<ci> x </ci>
</apply>
</apply>
<apply>
<times/>
<ci> g </ci>
<ci> x </ci>
</apply>
</apply>
</apply>
<apply id="flcc_calculation">
<eq/>
<ci> flcc </ci>
<apply>
<times/>
<ci> f </ci>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless"> 0.375 </cn>
<ci> fracLCCap </ci>
</apply>
<apply>
<plus/>
<ci> fracLCCapo </ci>
<cn cellml:units="dimensionless"> 0.625 </cn>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="L_type_Ca_channel_y_gate">
<variable units="dimensionless" public_interface="out" name="y" initial_value="0.96"/>
<variable units="dimensionless" name="y_infinity"/>
<variable units="dimensionless" name="tau_y"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="second" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="22">
<apply id="y_diff_eq">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> y </ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<ci> y_infinity </ci>
<ci> y </ci>
</apply>
<ci> tau_y </ci>
</apply>
</apply>
<apply id="y_infinity_calculation">
<eq/>
<ci> y_infinity </ci>
<apply>
<plus/>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 55.0 </cn>
</apply>
<cn cellml:units="millivolt"> 7.5 </cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 0.1 </cn>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<apply>
<minus/>
<ci> V </ci>
</apply>
<cn cellml:units="millivolt"> 21.0 </cn>
</apply>
<cn cellml:units="millivolt"> 6.0 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply id="tau_y_calculation">
<eq/>
<ci> tau_y </ci>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 0.02 </cn>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 0.3 </cn>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci> V </ci>
<cn cellml:units="millivolt"> 30.0 </cn>
</apply>
<cn cellml:units="millivolt"> 9.5 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="L_type_Ca_channel_z_gate">
<variable units="dimensionless" public_interface="out" name="z" initial_value="0.92"/>
<variable units="first_order_rate_constant" name="omega" initial_value="10.0"/>
<variable units="first_order_rate_constant" name="v_omega"/>
<variable units="second_order_rate_constant" name="v_gamma"/>
<variable units="first_order_rate_constant" public_interface="in" name="gamma"/>
<variable units="first_order_rate_constant" public_interface="in" name="g"/>
<variable units="dimensionless" private_interface="in" name="v"/>
<variable units="first_order_rate_constant" public_interface="in" name="flcc"/>
<variable units="dimensionless" public_interface="in" name="w"/>
<variable units="second" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="23">
<apply id="z_diff_eq">
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> z </ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci> v_omega </ci>
<apply>
<minus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<ci> z </ci>
</apply>
</apply>
<apply>
<times/>
<ci> v_gamma </ci>
<ci> z </ci>
</apply>
</apply>
</apply>
<apply id="v_omega_calculation">
<eq/>
<ci> v_omega </ci>
<apply>
<plus/>
<apply>
<times/>
<ci> omega </ci>
<apply>
<power/>
<apply>
<minus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<ci> w </ci>
</apply>
<cn cellml:units="dimensionless"> 4.0 </cn>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 0.5 </cn>
<ci> w </ci>
<apply>
<power/>
<apply>
<minus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<ci> w </ci>
</apply>
<cn cellml:units="dimensionless"> 3.0 </cn>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 0.25 </cn>
<apply>
<power/>
<ci> w </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
</apply>
<apply>
<power/>
<apply>
<minus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<ci> w </ci>
</apply>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 0.125 </cn>
<apply>
<power/>
<ci> w </ci>
<cn cellml:units="dimensionless"> 3.0 </cn>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<ci> w </ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 0.0625 </cn>
<apply>
<power/>
<ci> w </ci>
<cn cellml:units="dimensionless"> 4.0 </cn>
</apply>
</apply>
</apply>
</apply>
<apply id="v_gamma_calculation">
<eq/>
<ci> v_gamma </ci>
<apply>
<plus/>
<apply>
<times/>
<ci> gamma </ci>
<apply>
<power/>
<apply>
<minus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<ci> v </ci>
</apply>
<cn cellml:units="dimensionless"> 4.0 </cn>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci> v </ci>
<apply>
<power/>
<apply>
<minus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<ci> v </ci>
</apply>
<cn cellml:units="dimensionless"> 3.0 </cn>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 4.0 </cn>
<apply>
<power/>
<ci> v </ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
</apply>
<apply>
<power/>
<apply>
<minus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<ci> v </ci>
</apply>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 8.0 </cn>
<apply>
<power/>
<ci> v </ci>
<cn cellml:units="dimensionless"> 3.0 </cn>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<ci> v </ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless"> 16.0 </cn>
<apply>
<power/>
<ci> v </ci>
<cn cellml:units="dimensionless"> 4.0 </cn>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<divide/>
<ci> flcc </ci>
<ci> g </ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<!-- back to the Puglisi and Bers 2001 model. -->
<component name="time_independent_potassium_current">
<variable units="microA_per_microF" public_interface="out" name="I_Ki"/>
<variable units="millivolt" public_interface="out" private_interface="out" name="E_K"/>
<variable units="milliS_per_microF" name="g_Ki" initial_value="0.24"/>
<variable units="millivolt" public_interface="in" private_interface="out" name="V"/>
<variable units="millimolar" public_interface="in" name="Ko"/>
<variable units="millimolar" public_interface="in" name="Ki"/>
<variable units="joule_per_kilomole_kelvin" public_interface="in" name="R"/>
<variable units="kelvin" public_interface="in" name="T"/>
<variable units="coulomb_per_mole" public_interface="in" name="F"/>
<variable units="dimensionless" private_interface="in" name="Ki_infinity"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="24">
<apply>
<eq/>
<ci>E_K</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
<ci>F</ci>
</apply>
<apply>
<log/>
<apply>
<divide/>
<ci>Ko</ci>
<ci>Ki</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>I_Ki</ci>
<apply>
<times/>
<ci>g_Ki</ci>
<apply>
<root/>
<apply>
<divide/>
<ci>Ko</ci>
<cn cellml:units="millimolar"> 5.4 </cn>
</apply>
</apply>
<ci>Ki_infinity</ci>
<apply>
<minus/>
<ci>V</ci>
<ci>E_K</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="time_independent_potassium_current_Ki_gate">
<variable units="dimensionless" public_interface="out" name="Ki_infinity"/>
<variable units="first_order_rate_constant" name="alpha_Ki"/>
<variable units="first_order_rate_constant" name="beta_Ki"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="millivolt" public_interface="in" name="E_K"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="25">
<apply>
<eq/>
<ci>alpha_Ki</ci>
<apply>
<divide/>
<cn cellml:units="first_order_rate_constant">1.02</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="per_millivolt">0.2385</cn>
<apply>
<minus/>
<apply>
<minus/>
<ci>V</ci>
<ci>E_K</ci>
</apply>
<cn cellml:units="millivolt">59.215</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_Ki</ci>
<apply>
<divide/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="first_order_rate_constant">0.49124</cn>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="dimensionless">0.08032</cn>
<apply>
<plus/>
<apply>
<minus/>
<ci>V</ci>
<ci>E_K</ci>
</apply>
<cn cellml:units="millivolt">5.476</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<exp/>
<apply>
<times/>
<cn cellml:units="dimensionless">0.06175</cn>
<apply>
<plus/>
<apply>
<minus/>
<ci>V</ci>
<ci>E_K</ci>
</apply>
<cn cellml:units="millivolt">-594.31</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="per_millivolt">0.5143</cn>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<ci>V</ci>
<ci>E_K</ci>
</apply>
<cn cellml:units="millivolt">4.753</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Ki_infinity</ci>
<apply>
<divide/>
<ci>alpha_Ki</ci>
<apply>
<plus/>
<ci>alpha_Ki</ci>
<ci>beta_Ki</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="plateau_potassium_current">
<variable units="microA_per_microF" public_interface="out" name="I_Kp"/>
<variable units="milliS_per_microF" name="g_Kp" initial_value="0.008"/>
<variable units="dimensionless" name="Kp"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="millivolt" public_interface="in" name="E_K"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="26">
<apply>
<eq/>
<ci>Kp</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<cn cellml:units="millivolt">7.488</cn>
<ci>V</ci>
</apply>
<cn cellml:units="millivolt">5.98</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>I_Kp</ci>
<apply>
<times/>
<ci>g_Kp</ci>
<ci>Kp</ci>
<apply>
<minus/>
<ci>V</ci>
<ci>E_K</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="sarcolemmal_calcium_pump">
<variable units="microA_per_microF" public_interface="out" name="I_p_Ca"/>
<variable units="millimolar" name="K_mpCa" initial_value="0.5E-3"/>
<variable units="microA_per_microF" name="I_pCa" initial_value="1.15"/>
<variable units="millimolar" public_interface="in" name="Cai"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="27">
<apply>
<eq/>
<ci>I_p_Ca</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>I_pCa</ci>
<ci>Cai</ci>
</apply>
<apply>
<plus/>
<ci>K_mpCa</ci>
<ci>Cai</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="sodium_background_current">
<variable units="microA_per_microF" public_interface="out" name="I_Na_b"/>
<variable units="microA_per_microF" name="g_Nab" initial_value="1.18E-3"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="millivolt" public_interface="in" name="E_Na"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="28">
<apply>
<eq/>
<ci>I_Na_b</ci>
<apply>
<times/>
<ci>g_Nab</ci>
<apply>
<minus/>
<ci>V</ci>
<ci>E_Na</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="calcium_background_current">
<variable units="microA_per_microF" public_interface="out" name="I_Ca_b"/>
<variable units="microA_per_microF" name="g_Cab" initial_value="2.8E-3"/>
<variable units="millivolt" public_interface="out" name="E_Ca"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="joule_per_kilomole_kelvin" public_interface="in" name="R"/>
<variable units="kelvin" public_interface="in" name="T"/>
<variable units="coulomb_per_mole" public_interface="in" name="F"/>
<variable units="millimolar" public_interface="in" name="Cai"/>
<variable units="millimolar" public_interface="in" name="Cao"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="29">
<apply>
<eq/>
<ci>E_Ca</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>F</ci>
</apply>
</apply>
<apply>
<log/>
<apply>
<divide/>
<ci>Cao</ci>
<ci>Cai</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>I_Ca_b</ci>
<apply>
<times/>
<ci>g_Cab</ci>
<apply>
<minus/>
<ci>V</ci>
<ci>E_Ca</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="sodium_potassium_pump">
<variable units="microA_per_microF" public_interface="out" name="I_NaK"/>
<variable units="microA_per_microF" name="i_NaK" initial_value="1.1"/>
<variable units="dimensionless" name="f_NaK"/>
<variable units="millimolar" name="K_mNai" initial_value="10"/>
<variable units="millimolar" name="K_mKo" initial_value="1.5"/>
<variable units="dimensionless" name="sigma"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="joule_per_kilomole_kelvin" public_interface="in" name="R"/>
<variable units="kelvin" public_interface="in" name="T"/>
<variable units="coulomb_per_mole" public_interface="in" name="F"/>
<variable units="millimolar" public_interface="in" name="Nai"/>
<variable units="millimolar" public_interface="in" name="Nao"/>
<variable units="millimolar" public_interface="in" name="Ko"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="30">
<apply>
<eq/>
<ci>sigma</ci>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<cn cellml:units="dimensionless">7</cn>
</apply>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Nao</ci>
<cn cellml:units="dimensionless">67.3</cn>
</apply>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>f_NaK</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<times/>
<cn cellml:units="dimensionless">0.1245</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">0.1</cn>
</apply>
<ci>V</ci>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">0.0365</cn>
<ci>sigma</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<ci>V</ci>
</apply>
<ci>F</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>I_NaK</ci>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>i_NaK</ci>
<ci>f_NaK</ci>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<power/>
<apply>
<divide/>
<ci>K_mNai</ci>
<ci>Nai</ci>
</apply>
<cn cellml:units="dimensionless">1.5</cn>
</apply>
</apply>
</apply>
<ci>Ko</ci>
</apply>
<apply>
<plus/>
<ci>Ko</ci>
<ci>K_mKo</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="Na_Ca_exchanger">
<variable units="microA_per_microF" public_interface="out" name="I_NaCa"/>
<variable units="microA_per_microF" name="K_NaCa" initial_value="1483.0"/>
<variable units="millimolar" name="K_mNa" initial_value="87.5"/>
<variable units="dimensionless" name="K_sat" initial_value="0.1"/>
<variable units="dimensionless" name="eta" initial_value="0.35"/>
<variable units="millimolar" name="K_mCa" initial_value="1.38"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="joule_per_kilomole_kelvin" public_interface="in" name="R"/>
<variable units="kelvin" public_interface="in" name="T"/>
<variable units="coulomb_per_mole" public_interface="in" name="F"/>
<variable units="millimolar" public_interface="in" name="Nai"/>
<variable units="millimolar" public_interface="in" name="Nao"/>
<variable units="millimolar" public_interface="in" name="Cai"/>
<variable units="millimolar" public_interface="in" name="Cao"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="31">
<apply id="Na_Ca_exchanger">
<eq/>
<ci> I_NaCa </ci>
<apply>
<times/>
<ci> K_NaCa </ci>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<plus/>
<apply>
<power/>
<ci> K_mNa </ci>
<cn cellml:units="dimensionless"> 3.0 </cn>
</apply>
<apply>
<power/>
<ci> Nao </ci>
<cn cellml:units="dimensionless"> 3.0 </cn>
</apply>
</apply>
</apply>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<plus/>
<ci> K_mCa </ci>
<ci> Cao </ci>
</apply>
</apply>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<times/>
<ci> K_sat </ci>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<ci> eta </ci>
<cn cellml:units="dimensionless"> 1.0 </cn>
</apply>
<ci> V </ci>
<apply>
<divide/>
<ci> F </ci>
<apply>
<times/>
<ci> R </ci>
<ci> T </ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<apply>
<times/>
<apply>
<exp/>
<apply>
<times/>
<ci> eta </ci>
<ci> V </ci>
<apply>
<divide/>
<ci> F </ci>
<apply>
<times/>
<ci> R </ci>
<ci> T </ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<power/>
<ci> Nai </ci>
<cn cellml:units="dimensionless"> 3.0 </cn>
</apply>
<ci> Cao </ci>
</apply>
<apply>
<times/>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<ci> eta </ci>
<cn cellml:units="dimensionless"> 1.0 </cn>
</apply>
<ci> V </ci>
<apply>
<divide/>
<ci> F </ci>
<apply>
<times/>
<ci> R </ci>
<ci> T </ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<power/>
<ci> Nao </ci>
<cn cellml:units="dimensionless"> 3.0 </cn>
</apply>
<ci> Cai </ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="calcium_dynamics">
<variable units="flux" public_interface="out" name="I_rel"/>
<variable units="flux" public_interface="out" name="I_up"/>
<variable units="flux" public_interface="out" name="I_leak"/>
<variable units="flux" name="I_tr"/>
<variable units="millimolar" name="Ca_JSR" initial_value="1.92"/>
<variable units="millimolar" name="Ca_NSR" initial_value="1.92"/>
<variable units="dimensionless" public_interface="out" name="B_myo"/>
<variable units="first_order_rate_constant" name="G_rel"/>
<variable units="flux" name="G_rel_max" initial_value="60E3"/>
<variable units="first_order_rate_constant" name="G_rel_overload" initial_value="4000"/>
<variable units="second" name="tau_tr" initial_value="5.7E-4"/>
<variable units="millimolar" name="K_mrel" initial_value="0.8E-3"/>
<variable units="millimolar" name="delta_Ca_ith" initial_value="0.18E-3"/>
<variable units="millimolar" name="CSQN_max" initial_value="15"/>
<variable units="millimolar" name="K_mCSQN" initial_value="0.8"/>
<variable units="millimolar" name="K_mup" initial_value="3E-4"/>
<variable units="dimensionless" name="Km_up"/>
<variable units="flux" name="i_up" initial_value="4.7"/>
<variable units="millimolar" name="Ca_NSR_max" initial_value="15"/>
<variable units="millimolar" name="Km_Tn" initial_value="0.5128E-3"/>
<variable units="millimolar" name="Km_CMDN" initial_value="2.38E-3"/>
<variable units="millimolar" name="Km_INDO" initial_value="8.44E-4"/>
<variable units="millimolar" name="Tn_max" initial_value="0.07"/>
<variable units="millimolar" name="CMDN_max" initial_value="0.05"/>
<variable units="millimolar" name="INDO_max" initial_value="0.07"/>
<variable units="dimensionless" name="b_CMDN"/>
<variable units="dimensionless" name="b_Tn"/>
<variable units="dimensionless" name="b_INDO"/>
<variable units="dimensionless" name="ryr_on"/>
<variable units="dimensionless" name="ryr_off"/>
<variable units="second" name="trel"/>
<variable units="second" name="y44" initial_value="0.898"/>
<variable units="micromolar" name="SR_content"/>
<variable units="second" name="tau_on" initial_value="2E-3"/>
<variable units="second" name="tau_off" initial_value="2E-3"/>
<variable units="second" public_interface="in" name="time"/>
<variable units="millivolt" public_interface="in" name="V"/>
<variable units="coulomb_per_mole" public_interface="in" name="F"/>
<variable units="microlitre" public_interface="in" name="V_myo"/>
<variable units="microlitre" public_interface="in" name="V_JSR"/>
<variable units="microlitre" public_interface="in" name="V_NSR"/>
<variable units="cm2" public_interface="in" name="A_cap"/>
<variable units="microA_per_microF" public_interface="in" name="I_NaCa"/>
<variable units="microA_per_microF" public_interface="in" name="I_Ca"/>
<variable units="microA_per_microF" public_interface="in" name="I_p_Ca"/>
<variable units="microA_per_microF" public_interface="in" name="I_Ca_b"/>
<variable units="millimolar" public_interface="in" name="Cai"/>
<variable units="millimolar" public_interface="in" name="Cao"/>
<variable units="dimensionless" public_interface="in" name="fracPLB"/>
<variable units="dimensionless" public_interface="in" name="fracPLBo"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="32">
<apply>
<eq/>
<ci>ryr_on</ci>
<apply>
<minus/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>trel</ci>
</apply>
<ci>tau_on</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>ryr_off</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<ci>trel</ci>
</apply>
<ci>tau_off</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>I_rel</ci>
<apply>
<times/>
<ci>G_rel</ci>
<ci>ryr_on</ci>
<ci>ryr_off</ci>
<apply>
<minus/>
<ci>Ca_JSR</ci>
<ci>Cai</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>trel</ci>
<apply>
<plus/>
<ci>y44</ci>
<cn cellml:units="second"> 2E-3 </cn>
</apply>
</apply>
<apply>
<eq/>
<ci> G_rel </ci>
<apply>
<divide/>
<ci> G_rel_max </ci>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<cn cellml:units="microA_per_microF"> 5.0 </cn>
<ci>I_Ca</ci>
<ci>I_Ca_b</ci>
<ci>I_p_Ca</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>I_NaCa</ci>
</apply>
</apply>
<cn cellml:units="microA_per_microF"> 0.9 </cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>I_up</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>i_up</ci>
<apply>
<power/>
<ci>Cai</ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
</apply>
<apply>
<plus/>
<apply>
<power/>
<ci>K_mup</ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
<apply>
<power/>
<ci>Cai</ci>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Km_up</ci>
<apply>
<times/>
<ci>K_mup</ci>
<apply>
<divide/>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci>fracPLB</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<times/>
<cn cellml:units="dimensionless"> 2.0 </cn>
<ci>fracPLBo</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>I_leak</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>i_up</ci>
<ci>Ca_NSR_max</ci>
</apply>
<ci>Ca_NSR</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>I_tr</ci>
<apply>
<divide/>
<apply>
<minus/>
<ci>Ca_NSR</ci>
<ci>Ca_JSR</ci>
</apply>
<ci>tau_tr</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Ca_JSR</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<apply>
<times/>
<ci>CSQN_max</ci>
<ci>K_mCSQN</ci>
</apply>
<apply>
<power/>
<apply>
<plus/>
<ci>K_mCSQN</ci>
<ci>Ca_JSR</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<ci>I_tr</ci>
<ci>I_rel</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Ca_NSR</ci>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<ci>I_tr</ci>
</apply>
<ci>V_JSR</ci>
</apply>
<ci>V_NSR</ci>
</apply>
<ci>I_leak</ci>
</apply>
<ci>I_up</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>SR_content</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">1E3</cn>
<apply>
<plus/>
<apply>
<times/>
<apply>
<plus/>
<ci>Ca_JSR</ci>
<apply>
<divide/>
<ci>Ca_JSR</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<apply>
<times/>
<ci>CSQN_max</ci>
<ci>K_mCSQN</ci>
</apply>
<apply>
<power/>
<apply>
<plus/>
<ci>K_mCSQN</ci>
<ci>Ca_JSR</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<divide/>
<ci>V_JSR</ci>
<ci>V_myo</ci>
</apply>
</apply>
<apply>
<times/>
<ci>Ca_NSR</ci>
<apply>
<divide/>
<ci>V_NSR</ci>
<ci>V_myo</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci> B_myo </ci>
<apply>
<divide/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply>
<plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<ci> b_CMDN </ci>
<ci> b_Tn </ci>
<ci> b_Tn </ci>
<ci> b_INDO </ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci> b_CMDN </ci>
<apply>
<divide/>
<apply>
<times/>
<ci> CMDN_max </ci>
<ci> Km_CMDN </ci>
</apply>
<apply>
<power/>
<apply>
<plus/>
<ci> Km_CMDN </ci>
<ci> Cai </ci>
</apply>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci> b_Tn </ci>
<apply>
<divide/>
<apply>
<times/>
<ci> Tn_max </ci>
<ci> Km_Tn </ci>
</apply>
<apply>
<power/>
<apply>
<plus/>
<ci> Km_Tn </ci>
<ci> Cai </ci>
</apply>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci> b_INDO </ci>
<apply>
<divide/>
<apply>
<times/>
<ci> INDO_max </ci>
<ci> Km_INDO </ci>
</apply>
<apply>
<power/>
<apply>
<plus/>
<ci> Km_INDO </ci>
<ci> Cai </ci>
</apply>
<cn cellml:units="dimensionless"> 2.0 </cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci> y44 </ci>
</apply>
<piecewise>
<piece>
<apply>
<minus/>
<cn cellml:units="dimensionless"> 1 </cn>
<apply>
<times/>
<cn cellml:units="dimensionless"> 1E4 </cn>
<ci> y44 </ci>
</apply>
</apply>
<apply>
<gt/>
<ci> V </ci>
<cn cellml:units="millivolt"> 30E3 </cn>
</apply>
</piece>
<otherwise>
<cn cellml:units="second"> 1.0 </cn>
</otherwise>
</piecewise>
</apply>
</math>
</component>
<component name="ionic_concentrations">
<variable units="millimolar" public_interface="out" name="Nai" initial_value="16.0"/>
<variable units="millimolar" public_interface="out" name="Nao" initial_value="140.0"/>
<variable units="millimolar" public_interface="out" name="Cai" initial_value="1.58E-4"/>
<variable units="millimolar" public_interface="out" name="Cao" initial_value="1.8"/>
<variable units="millimolar" public_interface="out" name="Ki" initial_value="145.0"/>
<variable units="millimolar" public_interface="out" name="Ko" initial_value="5.4"/>
<variable units="dimensionless" public_interface="in" name="B_myo"/>
<variable units="second" public_interface="in" name="time"/>
<variable units="coulomb_per_mole" public_interface="in" name="F"/>
<variable units="microA_per_microF" public_interface="in" name="I_Na"/>
<variable units="microA_per_microF" public_interface="in" name="I_Na_b"/>
<variable units="microA_per_microF" public_interface="in" name="I_NaCa"/>
<variable units="microA_per_microF" public_interface="in" name="I_NaK"/>
<variable units="microA_per_microF" public_interface="in" name="I_Ca"/>
<variable units="microA_per_microF" public_interface="in" name="I_CaK"/>
<variable units="microA_per_microF" public_interface="in" name="I_p_Ca"/>
<variable units="microA_per_microF" public_interface="in" name="I_Ca_b"/>
<variable units="microA_per_microF" public_interface="in" name="I_ss"/>
<variable units="microA_per_microF" public_interface="in" name="I_Ki"/>
<variable units="microA_per_microF" public_interface="in" name="I_Kp"/>
<variable units="microA_per_microF" public_interface="in" name="I_to"/>
<variable units="flux" public_interface="in" name="I_leak"/>
<variable units="flux" public_interface="in" name="I_up"/>
<variable units="flux" public_interface="in" name="I_rel"/>
<variable units="microlitre" public_interface="in" name="V_myo"/>
<variable units="microlitre" public_interface="in" name="V_JSR"/>
<variable units="microlitre" public_interface="in" name="V_NSR"/>
<variable units="cm2" public_interface="in" name="A_cap"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="33">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Nai</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">-1.0E3</cn>
<apply>
<plus/>
<ci>I_Na</ci>
<ci>I_Na_b</ci>
<apply>
<times/>
<ci>I_NaCa</ci>
<cn cellml:units="dimensionless">3</cn>
</apply>
<apply>
<times/>
<ci>I_NaK</ci>
<cn cellml:units="dimensionless">3</cn>
</apply>
</apply>
<apply>
<divide/>
<ci>A_cap</ci>
<apply>
<times/>
<ci>V_myo</ci>
<ci>F</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Ki</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">-1.0E3</cn>
<apply>
<plus/>
<ci>I_to</ci>
<ci>I_ss</ci>
<ci>I_Ki</ci>
<ci>I_Kp</ci>
<ci>I_CaK</ci>
<ci>I_Kp</ci>
<ci>I_to</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>I_NaK</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
<apply>
<divide/>
<ci>A_cap</ci>
<apply>
<times/>
<ci>V_myo</ci>
<ci>F</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Cai</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">1E3</cn>
<apply>
<minus/>
<apply>
<plus/>
<ci>I_Ca</ci>
<ci>I_Ca_b</ci>
</apply>
<apply>
<plus/>
<ci>I_p_Ca</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>I_NaCa</ci>
</apply>
</apply>
</apply>
<apply>
<divide/>
<ci>A_cap</ci>
<apply>
<times/>
<ci>V_myo</ci>
<cn cellml:units="dimensionless">2</cn>
<ci>F</ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<apply>
<minus/>
<ci>I_up</ci>
<ci>I_leak</ci>
</apply>
<apply>
<divide/>
<ci>V_NSR</ci>
<ci>V_myo</ci>
</apply>
</apply>
</apply>
<apply>
<plus/>
<ci>B_myo</ci>
<apply>
<times/>
<ci>I_rel</ci>
<apply>
<divide/>
<ci>V_JSR</ci>
<ci>V_myo</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="model_constants">
<variable units="cm2" public_interface="out" name="A_cap" initial_value="1.534E-4"/>
<variable units="microlitre" public_interface="out" name="V_myo" initial_value="20.8E-6"/>
<variable units="microlitre" public_interface="out" name="V_JSR" initial_value="9.3E-8"/>
<variable units="microlitre" public_interface="out" name="V_NSR" initial_value="9.88E-7"/>
</component>
<group>
<relationship_ref relationship="containment"/>
<component_ref component="membrane">
<component_ref component="fast_sodium_current">
<component_ref component="fast_sodium_current_m_gate"/>
<component_ref component="fast_sodium_current_h_gate"/>
<component_ref component="fast_sodium_current_j_gate"/>
</component_ref>
<component_ref component="L_type_Ca_channel">
<component_ref component="L_type_Ca_channel_v_gate"/>
<component_ref component="L_type_Ca_channel_w_gate"/>
<component_ref component="L_type_Ca_channel_x_gate"/>
<component_ref component="L_type_Ca_channel_y_gate"/>
<component_ref component="L_type_Ca_channel_z_gate"/>
</component_ref>
<component_ref component="time_independent_potassium_current">
<component_ref component="time_independent_potassium_current_Ki_gate"/>
</component_ref>
<component_ref component="Ca_independent_transient_outward_K_current">
<component_ref component="Ca_independent_transient_outward_K_current_r_gate"/>
<component_ref component="Ca_independent_transient_outward_K_current_s_gate"/>
<component_ref component="Ca_independent_transient_outward_K_current_s_slow_gate"/>
</component_ref>
<component_ref component="steady_state_outward_K_current">
<component_ref component="steady_state_outward_K_current_r_ss_gate"/>
<component_ref component="steady_state_outward_K_current_s_ss_gate"/>
</component_ref>
<component_ref component="Na_Ca_exchanger"/>
<component_ref component="plateau_potassium_current"/>
<component_ref component="sarcolemmal_calcium_pump"/>
<component_ref component="sodium_background_current"/>
<component_ref component="calcium_background_current"/>
<component_ref component="sodium_potassium_pump"/>
<component_ref component="ionic_concentrations"/>
<component_ref component="calcium_dynamics"/>
</component_ref>
</group>
<group>
<relationship_ref relationship="encapsulation"/>
<component_ref component="fast_sodium_current">
<component_ref component="fast_sodium_current_m_gate"/>
<component_ref component="fast_sodium_current_h_gate"/>
<component_ref component="fast_sodium_current_j_gate"/>
</component_ref>
<component_ref component="L_type_Ca_channel">
<component_ref component="L_type_Ca_channel_v_gate"/>
<component_ref component="L_type_Ca_channel_w_gate"/>
<component_ref component="L_type_Ca_channel_x_gate"/>
<component_ref component="L_type_Ca_channel_y_gate"/>
<component_ref component="L_type_Ca_channel_z_gate"/>
</component_ref>
<component_ref component="time_independent_potassium_current">
<component_ref component="time_independent_potassium_current_Ki_gate"/>
</component_ref>
<component_ref component="Ca_independent_transient_outward_K_current">
<component_ref component="Ca_independent_transient_outward_K_current_r_gate"/>
<component_ref component="Ca_independent_transient_outward_K_current_s_gate"/>
<component_ref component="Ca_independent_transient_outward_K_current_s_slow_gate"/>
</component_ref>
<component_ref component="steady_state_outward_K_current">
<component_ref component="steady_state_outward_K_current_r_ss_gate"/>
<component_ref component="steady_state_outward_K_current_s_ss_gate"/>
</component_ref>
</group>
<!-- EP model connections -->
<connection>
<map_components component_2="environment" component_1="membrane"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="Ca_independent_transient_outward_K_current"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="steady_state_outward_K_current"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="fast_sodium_current"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="L_type_Ca_channel"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="ionic_concentrations"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="calcium_dynamics"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="fast_sodium_current" component_1="membrane"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="I_Na" variable_1="I_Na"/>
<map_variables variable_2="R" variable_1="R"/>
<map_variables variable_2="T" variable_1="T"/>
<map_variables variable_2="F" variable_1="F"/>
</connection>
<connection>
<map_components component_2="Ca_independent_transient_outward_K_current" component_1="membrane"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="I_to" variable_1="I_to"/>
</connection>
<connection>
<map_components component_2="steady_state_outward_K_current" component_1="membrane"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="I_ss" variable_1="I_ss"/>
</connection>
<connection>
<map_components component_2="L_type_Ca_channel" component_1="membrane"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="I_Ca" variable_1="I_Ca"/>
<map_variables variable_2="I_CaK" variable_1="I_CaK"/>
<map_variables variable_2="R" variable_1="R"/>
<map_variables variable_2="T" variable_1="T"/>
<map_variables variable_2="F" variable_1="F"/>
</connection>
<connection>
<map_components component_2="Na_Ca_exchanger" component_1="membrane"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="I_NaCa" variable_1="I_NaCa"/>
<map_variables variable_2="R" variable_1="R"/>
<map_variables variable_2="T" variable_1="T"/>
<map_variables variable_2="F" variable_1="F"/>
</connection>
<connection>
<map_components component_2="time_independent_potassium_current" component_1="membrane"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="I_Ki" variable_1="I_Ki"/>
<map_variables variable_2="R" variable_1="R"/>
<map_variables variable_2="T" variable_1="T"/>
<map_variables variable_2="F" variable_1="F"/>
</connection>
<connection>
<map_components component_2="plateau_potassium_current" component_1="membrane"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="I_Kp" variable_1="I_Kp"/>
</connection>
<connection>
<map_components component_2="sarcolemmal_calcium_pump" component_1="membrane"/>
<map_variables variable_2="I_p_Ca" variable_1="I_p_Ca"/>
</connection>
<connection>
<map_components component_2="sodium_background_current" component_1="membrane"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="I_Na_b" variable_1="I_Na_b"/>
</connection>
<connection>
<map_components component_2="calcium_background_current" component_1="membrane"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="I_Ca_b" variable_1="I_Ca_b"/>
<map_variables variable_2="R" variable_1="R"/>
<map_variables variable_2="T" variable_1="T"/>
<map_variables variable_2="F" variable_1="F"/>
</connection>
<connection>
<map_components component_2="sodium_potassium_pump" component_1="membrane"/>
<map_variables variable_2="V" variable_1="V"/>
<map_variables variable_2="I_NaK" variable_1="I_NaK"/>
<map_variables variable_2="R" variable_1="R"/>
<map_variables variable_2="T" variable_1="T"/>
<map_variables variable_2="F" variable_1="F"/>
</connection>
<connection>
<map_components component_2="ionic_concentrations" component_1="membrane"/>
<map_variables variable_2="F" variable_1="F"/>
</connection>
<connection>
<map_components component_2="calcium_dynamics" component_1="membrane"/>
<map_variables variable_2="F" variable_1="F"/>
<map_variables variable_2="V" variable_1="V"/>
</connection>
<connection>
<map_components component_2="ionic_concentrations" component_1="Ca_independent_transient_outward_K_current"/>
<map_variables variable_2="I_to" variable_1="I_to"/>
</connection>
<connection>
<map_components component_2="time_independent_potassium_current" component_1="Ca_independent_transient_outward_K_current"/>
<map_variables variable_2="E_K" variable_1="E_K"/>
</connection>
<connection>
<map_components component_2="time_independent_potassium_current" component_1="steady_state_outward_K_current"/>
<map_variables variable_2="E_K" variable_1="E_K"/>
</connection>
<connection>
<map_components component_2="ionic_concentrations" component_1="steady_state_outward_K_current"/>
<map_variables variable_2="I_ss" variable_1="I_ss"/>
</connection>
<connection>
<map_components component_2="ionic_concentrations" component_1="fast_sodium_current"/>
<map_variables variable_2="I_Na" variable_1="I_Na"/>
<map_variables variable_2="Nao" variable_1="Nao"/>
<map_variables variable_2="Nai" variable_1="Nai"/>
</connection>
<connection>
<map_components component_2="sodium_background_current" component_1="fast_sodium_current"/>
<map_variables variable_2="E_Na" variable_1="E_Na"/>
</connection>
<connection>
<map_components component_2="ionic_concentrations" component_1="L_type_Ca_channel"/>
<map_variables variable_2="Ko" variable_1="Ko"/>
<map_variables variable_2="Ki" variable_1="Ki"/>
<map_variables variable_2="Cao" variable_1="Cao"/>
<map_variables variable_2="Cai" variable_1="Cai"/>
<map_variables variable_2="I_Ca" variable_1="I_Ca"/>
<map_variables variable_2="I_CaK" variable_1="I_CaK"/>
</connection>
<connection>
<map_components component_2="calcium_dynamics" component_1="L_type_Ca_channel"/>
<map_variables variable_2="I_Ca" variable_1="I_Ca"/>
</connection>
<connection>
<map_components component_2="ionic_concentrations" component_1="Na_Ca_exchanger"/>
<map_variables variable_2="Nai" variable_1="Nai"/>
<map_variables variable_2="Nao" variable_1="Nao"/>
<map_variables variable_2="Cai" variable_1="Cai"/>
<map_variables variable_2="Cao" variable_1="Cao"/>
<map_variables variable_2="I_NaCa" variable_1="I_NaCa"/>
</connection>
<connection>
<map_components component_2="calcium_dynamics" component_1="Na_Ca_exchanger"/>
<map_variables variable_2="I_NaCa" variable_1="I_NaCa"/>
</connection>
<connection>
<map_components component_2="ionic_concentrations" component_1="time_independent_potassium_current"/>
<map_variables variable_2="Ko" variable_1="Ko"/>
<map_variables variable_2="Ki" variable_1="Ki"/>
<map_variables variable_2="I_Ki" variable_1="I_Ki"/>
</connection>
<connection>
<map_components component_2="time_independent_potassium_current" component_1="plateau_potassium_current"/>
<map_variables variable_2="E_K" variable_1="E_K"/>
</connection>
<connection>
<map_components component_2="ionic_concentrations" component_1="plateau_potassium_current"/>
<map_variables variable_2="I_Kp" variable_1="I_Kp"/>
</connection>
<connection>
<map_components component_2="calcium_dynamics" component_1="sarcolemmal_calcium_pump"/>
<map_variables variable_2="I_p_Ca" variable_1="I_p_Ca"/>
</connection>
<connection>
<map_components component_2="ionic_concentrations" component_1="sarcolemmal_calcium_pump"/>
<map_variables variable_2="I_p_Ca" variable_1="I_p_Ca"/>
<map_variables variable_2="Cai" variable_1="Cai"/>
</connection>
<connection>
<map_components component_2="ionic_concentrations" component_1="sodium_background_current"/>
<map_variables variable_2="I_Na_b" variable_1="I_Na_b"/>
</connection>
<connection>
<map_components component_2="calcium_dynamics" component_1="calcium_background_current"/>
<map_variables variable_2="I_Ca_b" variable_1="I_Ca_b"/>
</connection>
<connection>
<map_components component_2="ionic_concentrations" component_1="calcium_background_current"/>
<map_variables variable_2="Cai" variable_1="Cai"/>
<map_variables variable_2="Cao" variable_1="Cao"/>
<map_variables variable_2="I_Ca_b" variable_1="I_Ca_b"/>
</connection>
<connection>
<map_components component_2="ionic_concentrations" component_1="sodium_potassium_pump"/>
<map_variables variable_2="Nai" variable_1="Nai"/>
<map_variables variable_2="Nao" variable_1="Nao"/>
<map_variables variable_2="Ko" variable_1="Ko"/>
<map_variables variable_2="I_NaK" variable_1="I_NaK"/>
</connection>
<connection>
<map_components component_2="fast_sodium_current_m_gate" component_1="fast_sodium_current"/>
<map_variables variable_2="m" variable_1="m"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="V" variable_1="V"/>
</connection>
<connection>
<map_components component_2="fast_sodium_current_h_gate" component_1="fast_sodium_current"/>
<map_variables variable_2="h" variable_1="h"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="V" variable_1="V"/>
</connection>
<connection>
<map_components component_2="fast_sodium_current_j_gate" component_1="fast_sodium_current"/>
<map_variables variable_2="j" variable_1="j"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="V" variable_1="V"/>
</connection>
<connection>
<map_components component_2="L_type_Ca_channel_y_gate" component_1="L_type_Ca_channel"/>
<map_variables variable_2="y" variable_1="y"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="V" variable_1="V"/>
</connection>
<connection>
<map_components component_2="L_type_Ca_channel_v_gate" component_1="L_type_Ca_channel"/>
<map_variables variable_2="v" variable_1="v"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="alpha" variable_1="alpha"/>
<map_variables variable_2="beta" variable_1="beta"/>
</connection>
<connection>
<map_components component_2="L_type_Ca_channel_w_gate" component_1="L_type_Ca_channel"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="alpha" variable_1="alpha"/>
<map_variables variable_2="beta" variable_1="beta"/>
</connection>
<connection>
<map_components component_2="L_type_Ca_channel_x_gate" component_1="L_type_Ca_channel"/>
<map_variables variable_2="x" variable_1="x"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="fracLCCap" variable_1="fracLCCap"/>
<map_variables variable_2="fracLCCapo" variable_1="fracLCCapo"/>
</connection>
<connection>
<map_components component_2="L_type_Ca_channel_z_gate" component_1="L_type_Ca_channel"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="z" variable_1="z"/>
<map_variables variable_2="gamma" variable_1="gamma"/>
</connection>
<connection>
<map_components component_2="L_type_Ca_channel_x_gate" component_1="L_type_Ca_channel_z_gate"/>
<map_variables variable_2="g" variable_1="g"/>
<map_variables variable_2="flcc" variable_1="flcc"/>
</connection>
<connection>
<map_components component_2="L_type_Ca_channel_w_gate" component_1="L_type_Ca_channel_z_gate"/>
<map_variables variable_2="w" variable_1="w"/>
</connection>
<connection>
<map_components component_2="L_type_Ca_channel_v_gate" component_1="L_type_Ca_channel_z_gate"/>
<map_variables variable_2="v" variable_1="v"/>
</connection>
<connection>
<map_components component_2="time_independent_potassium_current_Ki_gate" component_1="time_independent_potassium_current"/>
<map_variables variable_2="Ki_infinity" variable_1="Ki_infinity"/>
<map_variables variable_2="E_K" variable_1="E_K"/>
<map_variables variable_2="V" variable_1="V"/>
</connection>
<connection>
<map_components component_2="Ca_independent_transient_outward_K_current_r_gate" component_1="Ca_independent_transient_outward_K_current"/>
<map_variables variable_2="r" variable_1="r"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="V" variable_1="V"/>
</connection>
<connection>
<map_components component_2="Ca_independent_transient_outward_K_current_s_gate" component_1="Ca_independent_transient_outward_K_current"/>
<map_variables variable_2="s" variable_1="s"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="V" variable_1="V"/>
</connection>
<connection>
<map_components component_2="Ca_independent_transient_outward_K_current_s_slow_gate" component_1="Ca_independent_transient_outward_K_current"/>
<map_variables variable_2="s_slow" variable_1="s_slow"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="V" variable_1="V"/>
</connection>
<connection>
<map_components component_2="steady_state_outward_K_current_r_ss_gate" component_1="steady_state_outward_K_current"/>
<map_variables variable_2="r_ss" variable_1="r_ss"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="V" variable_1="V"/>
</connection>
<connection>
<map_components component_2="steady_state_outward_K_current_s_ss_gate" component_1="steady_state_outward_K_current"/>
<map_variables variable_2="s_ss" variable_1="s_ss"/>
<map_variables variable_2="time" variable_1="time"/>
<map_variables variable_2="V" variable_1="V"/>
</connection>
<connection>
<map_components component_2="ionic_concentrations" component_1="calcium_dynamics"/>
<map_variables variable_2="I_leak" variable_1="I_leak"/>
<map_variables variable_2="I_up" variable_1="I_up"/>
<map_variables variable_2="I_rel" variable_1="I_rel"/>
<map_variables variable_2="Cai" variable_1="Cai"/>
<map_variables variable_2="Cao" variable_1="Cao"/>
<map_variables variable_2="B_myo" variable_1="B_myo"/>
</connection>
<connection>
<map_components component_2="model_constants" component_1="calcium_dynamics"/>
<map_variables variable_2="A_cap" variable_1="A_cap"/>
<map_variables variable_2="V_NSR" variable_1="V_NSR"/>
<map_variables variable_2="V_JSR" variable_1="V_JSR"/>
<map_variables variable_2="V_myo" variable_1="V_myo"/>
</connection>
<connection>
<map_components component_2="model_constants" component_1="ionic_concentrations"/>
<map_variables variable_2="A_cap" variable_1="A_cap"/>
<map_variables variable_2="V_NSR" variable_1="V_NSR"/>
<map_variables variable_2="V_JSR" variable_1="V_JSR"/>
<map_variables variable_2="V_myo" variable_1="V_myo"/>
</connection>
<!-- Signalling model connections -->
<connection>
<map_components component_2="environment" component_1="beta_adrenergic_receptor_module"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="Gs_activation_module"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="cyclic_AMP_metabolism_module"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="PKA_activation_module"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="phospholamban_regulation_module"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="L_type_Ca_channel_regulation_module"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="Gs_activation_module" component_1="beta_adrenergic_receptor_module"/>
<map_variables variable_2="Gs_beta_gamma" variable_1="Gs_beta_gamma"/>
<map_variables variable_2="LRG" variable_1="LRG"/>
<map_variables variable_2="RG" variable_1="RG"/>
</connection>
<connection>
<map_components component_2="PKA_activation_module" component_1="beta_adrenergic_receptor_module"/>
<map_variables variable_2="PKACI" variable_1="PKACI"/>
</connection>
<connection>
<map_components component_2="Gs_activation_module" component_1="cyclic_AMP_metabolism_module"/>
<map_variables variable_2="Gs_alpha_GTPtot" variable_1="Gs_alpha_GTPtot"/>
</connection>
<connection>
<map_components component_2="PKA_activation_module" component_1="cyclic_AMP_metabolism_module"/>
<map_variables variable_2="cAMP" variable_1="cAMP"/>
<map_variables variable_2="cAMPtot" variable_1="cAMPtot"/>
</connection>
<connection>
<map_components component_2="phospholamban_regulation_module" component_1="PKA_activation_module"/>
<map_variables variable_2="PKACI" variable_1="PKACI"/>
</connection>
<connection>
<map_components component_2="L_type_Ca_channel_regulation_module" component_1="PKA_activation_module"/>
<map_variables variable_2="PKACII" variable_1="PKACII"/>
</connection>
<!-- connecting the two models -->
<connection>
<map_components component_2="L_type_Ca_channel_regulation_module" component_1="L_type_Ca_channel"/>
<map_variables variable_2="fracLCCbp" variable_1="fracLCCbp"/>
<map_variables variable_2="fracLCCbpo" variable_1="fracLCCbpo"/>
<map_variables variable_2="fracLCCap" variable_1="fracLCCap"/>
<map_variables variable_2="fracLCCapo" variable_1="fracLCCapo"/>
</connection>
<connection>
<map_components component_2="calcium_dynamics" component_1="phospholamban_regulation_module"/>
<map_variables variable_2="fracPLB" variable_1="fracPLB"/>
<map_variables variable_2="fracPLBo" variable_1="fracPLBo"/>
</connection>
<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
<rdf:Bag rdf:about="rdf:#a8556578-c204-404c-b0c9-2cfca30cb5f7">
<rdf:li>Cardiac Myocyte</rdf:li>
<rdf:li>metabolism</rdf:li>
<rdf:li>electrophysiology</rdf:li>
<rdf:li>mechanics</rdf:li>
</rdf:Bag>
<rdf:Seq rdf:about="rdf:#ad72bc9a-a2c0-48d8-9ed5-a63e549fef6c">
<rdf:li rdf:resource="rdf:#be31d7e4-9620-403d-bffb-cc0554221966"/>
<rdf:li rdf:resource="rdf:#53f1d965-ecf1-4cb5-9cc6-e720ae5578d1"/>
<rdf:li rdf:resource="rdf:#e419231f-fc6c-4175-b366-6494523e15a3"/>
<rdf:li rdf:resource="rdf:#3cec4432-f09e-492b-980d-41c3f325f8cd"/>
</rdf:Seq>
<rdf:Description rdf:about="rdf:#3cec4432-f09e-492b-980d-41c3f325f8cd">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#470c1dbb-3c39-4107-ba57-b0135c5eae1a"/>
</rdf:Description>
<rdf:Description rdf:about="">
<dc:publisher>The University of Auckland, Bioengineering Institute</dc:publisher>
<cmeta:comment rdf:resource="rdf:#f78bffe1-cd0f-4ec3-86c9-3f36f876f6cd"/>
<dcterms:created rdf:resource="rdf:#49f4d8ee-41b0-4032-9264-853069efd766"/>
<dc:creator rdf:resource="rdf:#a55f9784-78de-4dfe-a483-4e05b2c0c605"/>
<cmeta:modification rdf:resource="rdf:#f7c8a1c3-6c98-41bf-85e2-1de40485d219"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#53f1d965-ecf1-4cb5-9cc6-e720ae5578d1">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#81c1dfab-72c7-40d0-b187-ccf727ed171d"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#9db2ab2d-13dd-42aa-b8ef-adb3f0a0082e">
<dc:subject rdf:resource="rdf:#219fdb75-145a-498a-bf98-ec6a2091d1db"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#0de86d89-0ae0-44d6-8038-df99a786eded">
<vCard:Given>Anushka</vCard:Given>
<vCard:Family>Michailova</vCard:Family>
<vCard:Other>P</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#f7c8a1c3-6c98-41bf-85e2-1de40485d219">
<dcterms:modified rdf:resource="rdf:#5c4e73c9-cb5e-4e08-88f9-c26ba6fe32a1"/>
<rdf:value>The model has been drastically changed. An electrophysiological excitable cell component has been added, and the author's source code was used to add initial values etc.</rdf:value>
<cmeta:modifier rdf:resource="rdf:#0b6603ba-a263-405e-9372-eddfb8b986b4"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#470c1dbb-3c39-4107-ba57-b0135c5eae1a">
<vCard:Given>Andrew</vCard:Given>
<vCard:Family>McCulloch</vCard:Family>
<vCard:Other>D</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#be31d7e4-9620-403d-bffb-cc0554221966">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#481509e2-a814-4d47-8605-2db1a7cb08bf"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#5c4e73c9-cb5e-4e08-88f9-c26ba6fe32a1">
<dcterms:W3CDTF>2007-09-21T11:43:40+12:00</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#481509e2-a814-4d47-8605-2db1a7cb08bf">
<vCard:Given>Jeffrey</vCard:Given>
<vCard:Family>Saucerman</vCard:Family>
<vCard:Other>J</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#0efe1568-9510-478e-97b5-d33d7d62e77d">
<dcterms:W3CDTF>2003-09-12</dcterms:W3CDTF>
</rdf:Description>
<rdf:Description rdf:about="rdf:#219fdb75-145a-498a-bf98-ec6a2091d1db">
<bqs:subject_type>keyword</bqs:subject_type>
<rdf:value rdf:resource="rdf:#a8556578-c204-404c-b0c9-2cfca30cb5f7"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#ec188067-d335-4afc-b707-1f3c17a2fe08">
<dc:title>Journal of Biological Chemistry</dc:title>
</rdf:Description>
<rdf:Description rdf:about="rdf:#f78bffe1-cd0f-4ec3-86c9-3f36f876f6cd">
<dc:creator rdf:resource="rdf:#ef1fa92d-7447-4331-9ed0-e1e1d093cbdb"/>
<rdf:value>Although much closer now to the published model, this CellML version of the model still can't be run in PCEnv or COR due to the presence of circular arguments. It's hoped that this problem will be addressed by the further development of the tools.</rdf:value>
</rdf:Description>
<rdf:Description rdf:about="#saucerman_model_2003">
<dc:title>
Saucerman et al.'s 2003 mathematical model of beta-adrenergic control of cardiac myocyte contractility.
</dc:title>
<cmeta:bio_entity>Cardiac Myocyte</cmeta:bio_entity>
<cmeta:comment rdf:resource="rdf:#bc9da357-61a4-490f-98fa-682c8648607a"/>
<bqs:reference rdf:resource="rdf:#9db2ab2d-13dd-42aa-b8ef-adb3f0a0082e"/>
<bqs:reference rdf:resource="rdf:#ccc52e45-21b5-4448-84e5-f5a29a280b1a"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#98fa7ace-c01f-44d2-9b9d-4874eafd1744">
<vCard:Given>Catherine</vCard:Given>
<vCard:Family>Lloyd</vCard:Family>
<vCard:Other>May</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#2d548f1a-f7ac-4be9-9a87-c694dd1435b5">
<vCard:Orgname>The University of Auckland</vCard:Orgname>
<vCard:Orgunit>The Bioengineering Institute</vCard:Orgunit>
</rdf:Description>
<rdf:Description rdf:about="rdf:#a55f9784-78de-4dfe-a483-4e05b2c0c605">
<vCard:ORG rdf:resource="rdf:#2d548f1a-f7ac-4be9-9a87-c694dd1435b5"/>
<vCard:EMAIL rdf:resource="rdf:#23ff5da7-0c81-4389-8ef9-8d71fcc139ec"/>
<vCard:N rdf:resource="rdf:#415a7e10-3261-4cb8-b735-35883e7b7f43"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#ccc52e45-21b5-4448-84e5-f5a29a280b1a">
<bqs:Pubmed_id>12972422</bqs:Pubmed_id>
<bqs:JournalArticle rdf:resource="rdf:#3b339aac-040a-47a4-be05-c0ba8b381cb9"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#4f762dc2-c157-4e7f-8ec6-d38089f38348">
<vCard:FN>Catherine Lloyd</vCard:FN>
</rdf:Description>
<rdf:Description rdf:about="rdf:#ef1fa92d-7447-4331-9ed0-e1e1d093cbdb">
<vCard:FN>Catherine Lloyd</vCard:FN>
</rdf:Description>
<rdf:Description rdf:about="rdf:#23ff5da7-0c81-4389-8ef9-8d71fcc139ec">
<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:#bc9da357-61a4-490f-98fa-682c8648607a">
<dc:creator rdf:resource="rdf:#4f762dc2-c157-4e7f-8ec6-d38089f38348"/>
<rdf:value>This is the CellML description of Saucerman et al.'s 2003 mathematical model of beta-adrenergic control of cardiac myocyte contractility.</rdf:value>
</rdf:Description>
<rdf:Description rdf:about="rdf:#81c1dfab-72c7-40d0-b187-ccf727ed171d">
<vCard:Given>Laurence</vCard:Given>
<vCard:Family>Brunton</vCard:Family>
<vCard:Other>L</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#0b6603ba-a263-405e-9372-eddfb8b986b4">
<vCard:N rdf:resource="rdf:#98fa7ace-c01f-44d2-9b9d-4874eafd1744"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#415a7e10-3261-4cb8-b735-35883e7b7f43">
<vCard:Given>Catherine</vCard:Given>
<vCard:Family>Lloyd</vCard:Family>
<vCard:Other>May</vCard:Other>
</rdf:Description>
<rdf:Description rdf:about="rdf:#3b339aac-040a-47a4-be05-c0ba8b381cb9">
<dc:creator rdf:resource="rdf:#ad72bc9a-a2c0-48d8-9ed5-a63e549fef6c"/>
<dc:title>Modeling beta-adrenergic control of cardiac myocyte contractility in silico.</dc:title>
<bqs:volume>278</bqs:volume>
<bqs:first_page>47997</bqs:first_page>
<bqs:Journal rdf:resource="rdf:#ec188067-d335-4afc-b707-1f3c17a2fe08"/>
<dcterms:issued rdf:resource="rdf:#0efe1568-9510-478e-97b5-d33d7d62e77d"/>
<bqs:last_page>48003</bqs:last_page>
</rdf:Description>
<rdf:Description rdf:about="rdf:#e419231f-fc6c-4175-b366-6494523e15a3">
<rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
<vCard:N rdf:resource="rdf:#0de86d89-0ae0-44d6-8038-df99a786eded"/>
</rdf:Description>
<rdf:Description rdf:about="rdf:#49f4d8ee-41b0-4032-9264-853069efd766">
<dcterms:W3CDTF>2003-11-18T00:00:00+00:00</dcterms:W3CDTF>
</rdf:Description>
</rdf:RDF>
</model>