Location: Sachse, Moreno, Abildskov, 2008 @ 6d1bf8b84d44 / sachse_moreno_abildskov_2008_b.cellml

Author:
Tommy Yu <tommy.yu@auckland.ac.nz>
Date:
2010-08-25 17:36:53+12:00
Desc:
e-notation fix
Permanent Source URI:
https://models.cellml.org/workspace/sachse_moreno_abildskov_2008/rawfile/6d1bf8b84d447a2915482ea244d4dee219d7420d/sachse_moreno_abildskov_2008_b.cellml

<?xml version="1.0" encoding="utf-8"?>
<!--  FILE :  sachse_model_2007.xml

CREATED :  27th August 2007

LAST MODIFIED : 27th August 2007

AUTHOR :  Catherine Lloyd
          Bioengineering Institute
          The University of Auckland
          
MODEL STATUS :  This model conforms to the CellML 1.1 Specification.

DESCRIPTION :  This file contains a CellML description of Sachse, Moreno and Abildskov's 2007 mathematical model of the electrophysiology of fibroblasts and how they interact with myocytes.

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="sachse_model_2007" name="sachse_model_2007">

<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
  <articleinfo>
  <title>Electrophysiological Modeling of Fibroblasts and Their Interaction with Myocytes</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 runs in PCEnv, COR and JSim to recreate the results in the published paper.  The units have been checked and are consistent, and we would like to thank Frank Sachse, the original model author, for his advice and assistance in getting the CellML version of this model working.  This particular CellML model represents the fibroblast.
          </para>
  </section>
  <sect1 id="sec_structure">
  <title>Model Structure</title>

<para>
Fibroblasts are the most abundant cell type present in the heart.  They are arranged in the myocardium, as a three-dimensional network, where they are responsible for the secretion of extracellular matrix proteins, and also play a role in cardiac electrophysiology.  Although fibroblasts themselves are electrically in-excitable, they can contribute to the electrophysiology of neighbouring myocytes through fibroblast-myocyte coupling.  Although the electrical coupling of fibroblasts and myocytes has been demonstrated experimentally both in vivo and in vitro, the precise role of fibroblasts as electrical bridges between myocytes still remains to be defined.   
</para>

<para>
In the paper described here, Frank Sachse, Alonso Moreno and J. A. Abildskov address this issue by presenting mathematical models of cellular and multi-cellular electrophysiology.  A novel mathematical model of cardiac fibroblasts was developed on the basis of experimental data from whole cell patch clamp and polymerase chain reaction studies.  This model was then applied, together with models of ventricular myocytes (which were based on the <ulink url="http://www.cellml.org/models/pandit_clark_giles_demir_2001_version07">Pandit <emphasis>et al.</emphasis> 2001 model</ulink> of the rat left ventricular myocyte), to determine the electrophysiological effects of heterogeneous electrical coupling of cells.
</para>

<informalfigure float="0" id="fig_reaction_diagram">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>model diagram</title>
    </objectinfo>
    <imagedata fileref="sachse_2007.png"/>
  </imageobject>
</mediaobject>
<caption>Schematic diagram of the fibroblast cell model.  A Markovian model for the time and voltage dependent outward current (I<subscript>SHkr</subscript>) is embedded within an electrophysiological cell model.  The ion channel has six different states; five closed states (C0-4) and one open state (O)</caption>
</informalfigure>

<para>
The complete original paper reference is cited below:
</para>

<para>
Electrophysiological Modeling of Fibroblasts and Their Interaction with Myocytes, F.B. Sachse, A.P. Moreno, and J.A. Abildskov, 2007, <emphasis>Annals of Biomedical Engineering</emphasis>
<ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=pubmed&amp;cmd=Retrieve&amp;dopt=AbstractPlus&amp;list_uids=17999190&amp;query_hl=1&amp;itool=pubmed_docsum">PubMed ID: 17999190</ulink>
</para>

</sect1>
</article>
</documentation> 




  
  
  <units name="millivolt">
      <unit units="volt" prefix="milli"/>
   </units>
   <units name="microfarad">
      <unit units="farad" prefix="micro"/>
   </units>
   <units name="microsiemens">
      <unit units="siemens" prefix="micro"/>
   </units>
   <units name="nanoampere">
      <unit units="ampere" prefix="nano"/>
   </units>
   <units name="millimolar">
      <unit units="mole" prefix="milli"/>
      <unit units="litre" exponent="-1"/>
   </units>
   <units name="per_millimolar">
      <unit units="millimolar" exponent="-1"/>
   </units>
   <units name="flux">
      <unit units="millimolar"/>
      <unit units="second" exponent="-1"/>
   </units>
   <units name="first_order_rate_constant">
      <unit units="second" exponent="-1"/>
   </units>
   <units name="millijoule_per_kelvin_mole">
      <unit units="joule" prefix="milli"/>
      <unit units="kelvin" exponent="-1"/>
      <unit units="mole" exponent="-1"/>
   </units>
   <units name="coulomb_per_mole">
      <unit units="coulomb"/>
      <unit units="mole" exponent="-1"/>
   </units>
   <units name="microlitre_per_second">
      <unit units="litre" prefix="micro"/>
      <unit units="second" exponent="-1"/>
   </units>
   <component name="environment">
      <variable units="second" public_interface="out" name="time"/>
   </component>
   <component cmeta:id="eqn1" name="membrane">
      <variable units="millivolt" public_interface="out" name="Vm"/>
      <variable units="microfarad" name="Cm" initial_value="4.5e-6"/>
      <variable units="nanoampere" public_interface="in" name="I_b"/>
      <variable units="nanoampere" public_interface="in" name="I_Kir"/>
      <variable units="nanoampere" public_interface="in" name="I_Shkr"/>
      <variable units="nanoampere" public_interface="in" name="I_stim"/>
      <variable units="second" public_interface="in" name="time"/>
      <variable units="millivolt" name="VmReal" initial_value="-58.00"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>Vm</ci>
            <piecewise>
               <piece>
                  <cn type="e-notation" cellml:units="millivolt">1<sep/>-7</cn>
                  <apply>
                     <lt/>
                     <apply>
                        <abs/>
                        <ci>VmReal</ci>
                     </apply>
                     <cn type="e-notation" cellml:units="millivolt">1<sep/>-7</cn>
                  </apply>
               </piece>
               <otherwise>
                  <ci>VmReal</ci>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>VmReal</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <apply>
                     <minus/>
                     <apply>
                        <plus/>
                        <ci>I_Kir</ci>
                        <ci>I_Shkr</ci>
                        <ci>I_b</ci>
                     </apply>
                     <ci>I_stim</ci>
                  </apply>
               </apply>
               <ci>Cm</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component cmeta:id="eqn2" name="I_Kir">
      <variable units="nanoampere" public_interface="out" name="I_Kir"/>
      <variable units="microsiemens" name="GKir" initial_value="1e-3"/>
      <variable units="dimensionless" name="aKir" initial_value="0.94"/>
      <variable units="dimensionless" name="bKir" initial_value="1.26"/>
      <variable units="millivolt" name="EK"/>
      <variable units="dimensionless" name="OKir"/>
      <variable units="second" public_interface="in" name="time"/>
      <variable units="millivolt" public_interface="in" name="Vm"/>
      <variable units="millijoule_per_kelvin_mole" 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="Ko"/>
      <variable units="millimolar" public_interface="in" name="Ki"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>I_Kir</ci>
            <apply>
               <times/>
               <ci>GKir</ci>
               <ci>OKir</ci>
               <apply>
                  <root/>
                  <apply>
                     <times/>
                     <ci>Ko</ci>
                     <cn type="e-notation" cellml:units="per_millimolar">1<sep/>-3</cn>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <ci>Vm</ci>
                  <ci>EK</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>OKir</ci>
            <apply>
               <divide/>
               <cn cellml:units="dimensionless">1</cn>
               <apply>
                  <plus/>
                  <ci>aKir</ci>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <ci>bKir</ci>
                           <apply>
                              <minus/>
                              <ci>Vm</ci>
                              <ci>EK</ci>
                           </apply>
                           <ci>F</ci>
                        </apply>
                        <apply>
                           <times/>
                           <ci>R</ci>
                           <ci>T</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>EK</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>R</ci>
                     <ci>T</ci>
                  </apply>
                  <ci>F</ci>
               </apply>
               <apply>
                  <ln/>
                  <apply>
                     <divide/>
                     <ci>Ko</ci>
                     <ci>Ki</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component cmeta:id="eqn4" name="I_Shkr">
      <variable units="nanoampere" public_interface="out" name="I_Shkr"/>
      <variable units="microlitre_per_second" name="PShkr" initial_value="5.4e-9"/>
      <variable units="dimensionless" name="C0ShkrReal" initial_value="9.11e-1"/>
      <variable units="dimensionless" name="C1ShkrReal" initial_value="8.57e-2"/>
      <variable units="dimensionless" name="C2ShkrReal" initial_value="3.02e-3"/>
      <variable units="dimensionless" name="C3ShkrReal" initial_value="4.74e-5"/>
      <variable units="dimensionless" name="C4ShkrReal" initial_value="2.79e-7"/>
      <variable units="dimensionless" name="OShkrReal" initial_value="0"/>
      <variable units="dimensionless" name="C0Shkr"/>
      <variable units="dimensionless" name="C1Shkr"/>
      <variable units="dimensionless" name="C2Shkr"/>
      <variable units="dimensionless" name="C3Shkr"/>
      <variable units="dimensionless" name="C4Shkr"/>
      <variable units="dimensionless" name="OShkr"/>
      <variable units="first_order_rate_constant" name="kv"/>
      <variable units="first_order_rate_constant" name="k_v"/>
      <variable units="first_order_rate_constant" name="kv0" initial_value="30"/>
      <variable units="dimensionless" name="zv" initial_value="1.28"/>
      <variable units="first_order_rate_constant" name="k_v0" initial_value="2"/>
      <variable units="dimensionless" name="z_v" initial_value="-1.53"/>
      <variable units="first_order_rate_constant" name="ko" initial_value="77"/>
      <variable units="first_order_rate_constant" name="k_o" initial_value="18"/>
      <variable units="second" public_interface="in" name="time"/>
      <variable units="millivolt" public_interface="in" name="Vm"/>
      <variable units="millijoule_per_kelvin_mole" 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="Ko"/>
      <variable units="millimolar" public_interface="in" name="Ki"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>I_Shkr</ci>
            <apply>
               <divide/>
               <apply>
                  <times/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>PShkr</ci>
                        <ci>OShkr</ci>
                        <ci>Vm</ci>
                        <apply>
                           <power/>
                           <ci>F</ci>
                           <cn cellml:units="dimensionless">2</cn>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <ci>R</ci>
                        <ci>T</ci>
                     </apply>
                  </apply>
                  <apply>
                     <minus/>
                     <ci>Ki</ci>
                     <apply>
                        <times/>
                        <ci>Ko</ci>
                        <apply>
                           <exp/>
                           <apply>
                              <divide/>
                              <apply>
                                 <times/>
                                 <apply>
                                    <minus/>
                                    <ci>Vm</ci>
                                 </apply>
                                 <ci>F</ci>
                              </apply>
                              <apply>
                                 <times/>
                                 <ci>R</ci>
                                 <ci>T</ci>
                              </apply>
                           </apply>
                        </apply>
                     </apply>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <exp/>
                     <apply>
                        <divide/>
                        <apply>
                           <times/>
                           <apply>
                              <minus/>
                              <ci>Vm</ci>
                           </apply>
                           <ci>F</ci>
                        </apply>
                        <apply>
                           <times/>
                           <ci>R</ci>
                           <ci>T</ci>
                        </apply>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>C0Shkr</ci>
            <piecewise>
               <piece>
                  <cn cellml:units="dimensionless">0</cn>
                  <apply>
                     <lt/>
                     <ci>C0ShkrReal</ci>
                     <cn cellml:units="dimensionless">0</cn>
                  </apply>
               </piece>
               <piece>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <gt/>
                     <ci>C0ShkrReal</ci>
                     <cn cellml:units="dimensionless">1</cn>
                  </apply>
               </piece>
               <otherwise>
                  <ci>C0ShkrReal</ci>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>C1Shkr</ci>
            <piecewise>
               <piece>
                  <cn cellml:units="dimensionless">0</cn>
                  <apply>
                     <lt/>
                     <ci>C1ShkrReal</ci>
                     <cn cellml:units="dimensionless">0</cn>
                  </apply>
               </piece>
               <piece>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <gt/>
                     <ci>C1ShkrReal</ci>
                     <cn cellml:units="dimensionless">1</cn>
                  </apply>
               </piece>
               <otherwise>
                  <ci>C1ShkrReal</ci>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>C2Shkr</ci>
            <piecewise>
               <piece>
                  <cn cellml:units="dimensionless">0</cn>
                  <apply>
                     <lt/>
                     <ci>C2ShkrReal</ci>
                     <cn cellml:units="dimensionless">0</cn>
                  </apply>
               </piece>
               <piece>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <gt/>
                     <ci>C2ShkrReal</ci>
                     <cn cellml:units="dimensionless">1</cn>
                  </apply>
               </piece>
               <otherwise>
                  <ci>C2ShkrReal</ci>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>C3Shkr</ci>
            <piecewise>
               <piece>
                  <cn cellml:units="dimensionless">0</cn>
                  <apply>
                     <lt/>
                     <ci>C3ShkrReal</ci>
                     <cn cellml:units="dimensionless">0</cn>
                  </apply>
               </piece>
               <piece>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <gt/>
                     <ci>C3ShkrReal</ci>
                     <cn cellml:units="dimensionless">1</cn>
                  </apply>
               </piece>
               <otherwise>
                  <ci>C3ShkrReal</ci>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>C4Shkr</ci>
            <piecewise>
               <piece>
                  <cn cellml:units="dimensionless">0</cn>
                  <apply>
                     <lt/>
                     <ci>C4ShkrReal</ci>
                     <cn cellml:units="dimensionless">0</cn>
                  </apply>
               </piece>
               <piece>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <gt/>
                     <ci>C4ShkrReal</ci>
                     <cn cellml:units="dimensionless">1</cn>
                  </apply>
               </piece>
               <otherwise>
                  <ci>C4ShkrReal</ci>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>OShkr</ci>
            <piecewise>
               <piece>
                  <cn cellml:units="dimensionless">0</cn>
                  <apply>
                     <lt/>
                     <ci>OShkrReal</ci>
                     <cn cellml:units="dimensionless">0</cn>
                  </apply>
               </piece>
               <piece>
                  <cn cellml:units="dimensionless">1</cn>
                  <apply>
                     <gt/>
                     <ci>OShkrReal</ci>
                     <cn cellml:units="dimensionless">1</cn>
                  </apply>
               </piece>
               <otherwise>
                  <ci>OShkrReal</ci>
               </otherwise>
            </piecewise>
         </apply>
         <apply>
            <eq/>
            <ci>kv</ci>
            <apply>
               <times/>
               <ci>kv0</ci>
               <apply>
                  <exp/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>Vm</ci>
                        <ci>zv</ci>
                        <ci>F</ci>
                     </apply>
                     <apply>
                        <times/>
                        <ci>R</ci>
                        <ci>T</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>k_v</ci>
            <apply>
               <times/>
               <ci>k_v0</ci>
               <apply>
                  <exp/>
                  <apply>
                     <divide/>
                     <apply>
                        <times/>
                        <ci>Vm</ci>
                        <ci>z_v</ci>
                        <ci>F</ci>
                     </apply>
                     <apply>
                        <times/>
                        <ci>R</ci>
                        <ci>T</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>C0ShkrReal</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>k_v</ci>
                  <ci>C1Shkr</ci>
               </apply>
               <apply>
                  <times/>
                  <cn cellml:units="dimensionless">4</cn>
                  <ci>kv</ci>
                  <ci>C0Shkr</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>C1ShkrReal</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">2</cn>
                     <ci>k_v</ci>
                     <ci>C2Shkr</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">4</cn>
                     <ci>kv</ci>
                     <ci>C0Shkr</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <apply>
                     <plus/>
                     <apply>
                        <times/>
                        <cn cellml:units="dimensionless">3</cn>
                        <ci>kv</ci>
                     </apply>
                     <ci>k_v</ci>
                  </apply>
                  <ci>C1Shkr</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>C2ShkrReal</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">3</cn>
                     <ci>k_v</ci>
                     <ci>C3Shkr</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">3</cn>
                     <ci>kv</ci>
                     <ci>C1Shkr</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <apply>
                     <plus/>
                     <apply>
                        <times/>
                        <cn cellml:units="dimensionless">2</cn>
                        <ci>kv</ci>
                     </apply>
                     <apply>
                        <times/>
                        <cn cellml:units="dimensionless">2</cn>
                        <ci>k_v</ci>
                     </apply>
                  </apply>
                  <ci>C2Shkr</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>C3ShkrReal</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">4</cn>
                     <ci>k_v</ci>
                     <ci>C4Shkr</ci>
                  </apply>
                  <apply>
                     <times/>
                     <cn cellml:units="dimensionless">2</cn>
                     <ci>kv</ci>
                     <ci>C2Shkr</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <apply>
                     <plus/>
                     <ci>kv</ci>
                     <apply>
                        <times/>
                        <cn cellml:units="dimensionless">3</cn>
                        <ci>k_v</ci>
                     </apply>
                  </apply>
                  <ci>C3Shkr</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>C4ShkrReal</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <plus/>
                  <apply>
                     <times/>
                     <ci>k_o</ci>
                     <ci>OShkr</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>kv</ci>
                     <ci>C3Shkr</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <apply>
                     <plus/>
                     <ci>ko</ci>
                     <apply>
                        <times/>
                        <cn cellml:units="dimensionless">4</cn>
                        <ci>k_v</ci>
                     </apply>
                  </apply>
                  <ci>C4Shkr</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>OShkrReal</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>ko</ci>
                  <ci>C4Shkr</ci>
               </apply>
               <apply>
                  <times/>
                  <ci>k_o</ci>
                  <ci>OShkr</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component cmeta:id="eqn14" name="I_b">
      <variable units="nanoampere" public_interface="out" name="I_b"/>
      <variable units="microsiemens" name="Gb" initial_value="6.9e-6"/>
      <variable units="millivolt" name="Eb" initial_value="0"/>
      <variable units="second" public_interface="in" name="time"/>
      <variable units="millivolt" public_interface="in" name="Vm"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>I_b</ci>
            <apply>
               <times/>
               <ci>Gb</ci>
               <apply>
                  <minus/>
                  <ci>Vm</ci>
                  <ci>Eb</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="I_stim">
      <variable units="nanoampere" public_interface="out" name="I_stim"/>
      <variable units="second" name="stim_start" initial_value="0.1"/>
      <variable units="second" name="stim_end" initial_value="10"/>
      <variable units="second" name="stim_period" initial_value="1"/>
      <variable units="second" name="stim_duration" initial_value="0.001"/>
      <variable units="nanoampere" name="stim_amplitude" initial_value="0.1"/>
      <variable units="second" public_interface="in" name="time"/>
      <variable units="millivolt" public_interface="in" name="Vm"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>I_stim</ci>
            <piecewise>
               <piece>
                  <ci>stim_amplitude</ci>
                  <apply>
                     <and/>
                     <apply>
                        <geq/>
                        <ci>time</ci>
                        <ci>stim_start</ci>
                     </apply>
                  <apply>
                     <and/>
                     <apply>
                        <leq/>
                        <ci>time</ci>
                        <ci>stim_end</ci>
                     </apply>
                     <apply>
                        <leq/>
                        <apply>
                           <minus/>
                           <apply>
                              <minus/>
                              <ci>time</ci>
                              <ci>stim_start</ci>
                           </apply>
                           <apply>
                              <times/>
                              <apply>
                                 <floor/>
                                 <apply>
                                    <divide/>
                                    <apply>
                                       <minus/>
                                       <ci>time</ci>
                                       <ci>stim_start</ci>
                                    </apply>
                                    <ci>stim_period</ci>
                                 </apply>
                              </apply>
                              <ci>stim_period</ci>
                           </apply>
                        </apply>
                        <ci>stim_duration</ci>
                     </apply>
                  </apply>
                  </apply>
               </piece>
               <otherwise>
                  <cn cellml:units="nanoampere">0</cn>
               </otherwise>
            </piecewise>
         </apply>
      </math>
   </component>
   <component name="model_parameters">
      <variable units="millijoule_per_kelvin_mole" public_interface="out" name="R" initial_value="8310"/>
      <variable units="kelvin" public_interface="out" name="T" initial_value="295"/>
      <variable units="coulomb_per_mole" public_interface="out" name="F" initial_value="9.65e4"/>
      <variable units="millimolar" public_interface="out" name="Ki" initial_value="140"/>
      <variable units="millimolar" public_interface="out" name="Ko" initial_value="5"/>
   </component>
   <group>
      <relationship_ref relationship="containment"/>
      <component_ref component="membrane">
         <component_ref component="I_Kir"/>
         <component_ref component="I_Shkr"/>
         <component_ref component="I_b"/>
         <component_ref component="I_stim"/>
      </component_ref>
   </group>
   <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="I_Kir"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="I_Shkr"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="I_b"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="environment" component_1="I_stim"/>
      <map_variables variable_2="time" variable_1="time"/>
   </connection>
   <connection>
      <map_components component_2="I_Kir" component_1="membrane"/>
      <map_variables variable_2="I_Kir" variable_1="I_Kir"/>
      <map_variables variable_2="Vm" variable_1="Vm"/>
   </connection>
   <connection>
      <map_components component_2="I_Shkr" component_1="membrane"/>
      <map_variables variable_2="I_Shkr" variable_1="I_Shkr"/>
      <map_variables variable_2="Vm" variable_1="Vm"/>
   </connection>
   <connection>
      <map_components component_2="I_b" component_1="membrane"/>
      <map_variables variable_2="I_b" variable_1="I_b"/>
      <map_variables variable_2="Vm" variable_1="Vm"/>
   </connection>
   <connection>
      <map_components component_2="I_stim" component_1="membrane"/>
      <map_variables variable_2="I_stim" variable_1="I_stim"/>
      <map_variables variable_2="Vm" variable_1="Vm"/>
   </connection>
   <connection>
      <map_components component_2="model_parameters" component_1="I_Kir"/>
      <map_variables variable_2="Ko" variable_1="Ko"/>
      <map_variables variable_2="Ki" variable_1="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="model_parameters" component_1="I_Shkr"/>
      <map_variables variable_2="Ko" variable_1="Ko"/>
      <map_variables variable_2="Ki" variable_1="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>

<rdf:RDF>
  <rdf:Bag rdf:about="rdf:#55d23512-91a0-445e-9d3d-cab978e3125a">
    <rdf:li>cardiac</rdf:li>
    <rdf:li>cardiac electrophysiology</rdf:li>
    <rdf:li>fibroblast</rdf:li>
    <rdf:li>myocyte</rdf:li>
    <rdf:li>electrophysiology</rdf:li>
  </rdf:Bag>
  <rdf:Seq rdf:about="rdf:#d5c5d1cc-51af-483f-847d-3c43e0168264">
    <rdf:li rdf:resource="rdf:#a2261697-686f-46c9-a716-68dcee58cd52"/>
    <rdf:li rdf:resource="rdf:#9cff9e41-61d2-402e-9ce7-42cfe6bcb072"/>
    <rdf:li rdf:resource="rdf:#28c60c48-a93f-452c-b8b4-eb5c7de81bce"/>
  </rdf:Seq>
  <rdf:Description rdf:about="rdf:#e1e29d37-53c7-41ed-b154-097036f4f741">
    <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="">
    <dc:title>Electrophysiological Modeling of Fibroblasts and Their Interaction with Myocytes (Fibroblast Model)</dc:title>
	<dc:publisher>The University of Auckland, Bioengineering Institute</dc:publisher>
    <cmeta:comment rdf:resource="rdf:#68403c49-50e2-46bc-b477-d6d452b89e9b"/>
    <dcterms:created rdf:resource="rdf:#7af5a239-6f54-4953-8183-cbbd6672f54d"/>
    <dc:creator rdf:resource="rdf:#cd5ec6f2-3907-4d4e-ac6f-cecbbc5f1f1b"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4978d931-3f7f-4c0a-a333-20565ac8fb3c">
    <dc:creator rdf:resource="rdf:#d5c5d1cc-51af-483f-847d-3c43e0168264"/>
    <dc:title>Electrophysiological Modeling of Fibroblasts and Their Interaction with Myocytes</dc:title>
    <bqs:volume>36(1)</bqs:volume>
    <bqs:first_page>41</bqs:first_page>
    <bqs:Journal rdf:resource="rdf:#bb12bf13-c1b2-4fbf-b732-c3a309c66dc6"/>
    <dcterms:issued rdf:resource="rdf:#5f606e73-9b8e-4a61-b437-265ac93e2a7e"/>
    <bqs:last_page>56</bqs:last_page>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#7ec84fe3-d2d8-4c5d-85f5-7f20a3dfae45">
    <vCard:Given>F</vCard:Given>
    <vCard:Family>Sachse</vCard:Family>
    <vCard:Other>B</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5f606e73-9b8e-4a61-b437-265ac93e2a7e">
    <dcterms:W3CDTF>2008-01-00 00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b316a018-8d9e-40e1-899d-62283fedb013">
    <vCard:Given>A</vCard:Given>
    <vCard:Family>Moreno</vCard:Family>
    <vCard:Other>P</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="#sachse_model_2007">
    <dc:title>
        Sachse, Moreno and Abildskov's 2007 mathematical model of the electrophysiology of fibroblasts and how they interact with myocytes.
      </dc:title>
    <cmeta:comment rdf:resource="rdf:#6eaaec91-79cf-422c-8760-66f1ef593dac"/>
    <bqs:reference rdf:resource="rdf:#096cd246-65fd-498d-853a-3dac6ffd54a9"/>
    <bqs:reference rdf:resource="rdf:#6be3d549-ce0a-4950-a38c-d611e5ba99f3"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#9cff9e41-61d2-402e-9ce7-42cfe6bcb072">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#b316a018-8d9e-40e1-899d-62283fedb013"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6be3d549-ce0a-4950-a38c-d611e5ba99f3">
    <dc:subject rdf:resource="rdf:#9e8d9d25-743d-4ca9-acc0-1e61ade60a44"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#096cd246-65fd-498d-853a-3dac6ffd54a9">
    <bqs:Pubmed_id>17999190</bqs:Pubmed_id>
    <bqs:JournalArticle rdf:resource="rdf:#4978d931-3f7f-4c0a-a333-20565ac8fb3c"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4401a32d-a2a2-43ff-82d7-a608ea8194cf">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#cd5ec6f2-3907-4d4e-ac6f-cecbbc5f1f1b">
    <vCard:ORG rdf:resource="rdf:#7297cb66-a158-4c77-a767-784cb41cfdf5"/>
    <vCard:EMAIL rdf:resource="rdf:#e1e29d37-53c7-41ed-b154-097036f4f741"/>
    <vCard:N rdf:resource="rdf:#4401a32d-a2a2-43ff-82d7-a608ea8194cf"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#9e8d9d25-743d-4ca9-acc0-1e61ade60a44">
    <bqs:subject_type>keyword</bqs:subject_type>
    <rdf:value rdf:resource="rdf:#55d23512-91a0-445e-9d3d-cab978e3125a"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#7297cb66-a158-4c77-a767-784cb41cfdf5">
    <vCard:Orgname>The University of Auckland</vCard:Orgname>
    <vCard:Orgunit>The Bioengineering Institute</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#6eaaec91-79cf-422c-8760-66f1ef593dac">
    <dc:creator rdf:resource="rdf:#0e63034d-3126-4a7b-9b66-7d4dea9299fa"/>
    <rdf:value>This is a CellML description of Sachse, Moreno and Abildskov's 2007 mathematical model of the electrophysiology of fibroblasts and how they interact with myocytes.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#bb12bf13-c1b2-4fbf-b732-c3a309c66dc6">
    <dc:title>Annals of Biomedical Engineering</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#28c60c48-a93f-452c-b8b4-eb5c7de81bce">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#648f2abb-8814-4a1a-b00f-87b7c3948478"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#7af5a239-6f54-4953-8183-cbbd6672f54d">
    <dcterms:W3CDTF>2007-11-16T00:00:00+00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#0e63034d-3126-4a7b-9b66-7d4dea9299fa">
    <vCard:FN>Catherine Lloyd</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#328118c2-53e7-48f0-a6c7-49baab713d4a">
    <vCard:FN>Catherine Lloyd</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#648f2abb-8814-4a1a-b00f-87b7c3948478">
    <vCard:Given>J</vCard:Given>
    <vCard:Family>Abildskov</vCard:Family>
    <vCard:Other>A</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#68403c49-50e2-46bc-b477-d6d452b89e9b">
    <dc:creator rdf:resource="rdf:#328118c2-53e7-48f0-a6c7-49baab713d4a"/>
    <rdf:value>This CellML version of the model runs in PCEnv, COR and JSim to recreate the results in the published paper.  The units have been checked and are consistent, and we would like to thank Frank Sachse, the original model author, for his advice and assistance in getting the CellML version of this model working.  This particular CellML model represents the fibroblast.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a2261697-686f-46c9-a716-68dcee58cd52">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#7ec84fe3-d2d8-4c5d-85f5-7f20a3dfae45"/>
  </rdf:Description>
</rdf:RDF>
</model>