Rendering of the source text

<?xml version="1.0"?>
<model xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#" cmeta:id="morris_1981" name="morris_1981">

<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
  <articleinfo>
  <title>Voltage oscillations in the barnacle giant muscle fiber</title>
  <author>
    <firstname>Catherine</firstname>
          <surname>Lloyd</surname>
    <affiliation>
      <shortaffil>Auckland Bioengineering Institute, The University of Auckland</shortaffil>
    </affiliation>
  </author>
</articleinfo>
     <section id="sec_status">
    <title>Model Status</title>
    <para>This CellML model represents the reduced model (equations 9 and 2) with parameters values taken from the legend of figure 9.  The model runs in both OpenCell and COR to replicate figure 9a (applied current = 400).  The units have been checked and they are consistent.
</para>
  </section>
  <sect1 id="sec_structure">
<title>Model Structure</title>

<para>
ABSTRACT: Barnacle muscle fibers subjected to constant current stimulation produce a variety of types of oscillatory behavior when the internal medium contains the Ca++ chelator EGTA. Oscillations are abolished if Ca++ is removed from the external medium, or if the K+ conductance is blocked. Available voltage-clamp data indicate that the cell's active conductance systems are exceptionally simple. Given the complexity of barnacle fiber voltage behavior, this seems paradoxical. This paper presents an analysis of the possible modes of behavior available to a system of two noninactivating conductance mechanisms, and indicates a good correspondence to the types of behavior exhibited by barnacle fiber. The differential equations of a simple equivalent circuit for the fiber are dealt with by means of some of the mathematical techniques of nonlinear mechanics. General features of the system are (a) a propensity to produce damped or sustained oscillations over a rather broad parameter range, and (b) considerable latitude in the shape of the oscillatory potentials. It is concluded that for cells subject to changeable parameters (either from cell to cell or with time during cellular activity), a system dominated by two noninactivating conductances can exhibit varied oscillatory and bistable behavior.
</para>

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

<para>
Voltage oscillations in the barnacle giant muscle fiber, Morris C and Lecar H, 1981, <emphasis>Biophysical Journal</emphasis>, 35, 193-213. <ulink url="http://www.ncbi.nlm.nih.gov/pubmed/7260316">PubMed ID: 7260316</ulink> 
</para>

<informalfigure float="0" id="fig_reaction_diagram">
<mediaobject>
  <imageobject>
    <imagedata fileref="morris_1981.png"/>
  </imageobject>
</mediaobject>
<caption>Schematic diagrams of the cell model.</caption>

</informalfigure>

</sect1>
</article>
</documentation>
    
  
   <units name="millisecond">
      <unit prefix="milli" units="second"/>
   </units>
   
   <units name="per_millisecond">
      <unit exponent="-1" units="millisecond"/>
   </units>
   
   <units name="millivolt">
      <unit prefix="milli" units="volt"/>
   </units>
   
   <units name="microF_per_cm2">
      <unit prefix="micro" units="farad"/>
	  <unit exponent="-2" prefix="centi" units="metre"/>
   </units>
   
   <units name="microA_per_cm2">
      <unit prefix="micro" units="ampere"/>
	  <unit exponent="-2" prefix="centi" units="metre"/>
   </units>
   
   <units name="milliS_per_cm2">
      <unit prefix="milli" units="siemens"/>
	  <unit exponent="-2" prefix="centi" units="metre"/>
   </units>
   
   
   
   <component name="environment">
      <variable cmeta:id="environment_time" name="time" public_interface="out" units="millisecond"/>
   </component>
   
   <component name="V">
      <variable cmeta:id="V_V" initial_value="-50" name="V" public_interface="out" units="millivolt"/>
      
	  <variable initial_value="20.0" name="C" units="microF_per_cm2"/>
      <variable initial_value="300.0" name="i_app" units="microA_per_cm2"/>
      
	  <variable name="i_L" public_interface="in" units="microA_per_cm2"/>
      <variable name="i_Ca" public_interface="in" units="microA_per_cm2"/>
      <variable name="i_K" public_interface="in" units="microA_per_cm2"/>
	  <variable name="time" public_interface="in" units="millisecond"/>
      
	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>V</ci>
            </apply>
            <apply>
               <divide/>
               <apply>
                  <minus/>
                  <ci>i_app</ci>
				  <apply>
                     <plus/>
                     <ci>i_L</ci>
                     <ci>i_Ca</ci>
                     <ci>i_K</ci>
                  </apply>
               </apply>
               <ci>C</ci>
            </apply>
         </apply>
      </math>
   </component>
  
   <component name="i_L">
      <variable cmeta:id="i_L_i_L" name="i_L" public_interface="out" units="microA_per_cm2"/>
      
	  <variable initial_value="2.0" name="g_L" units="milliS_per_cm2"/>
	  <variable initial_value="-50.00" name="E_L" units="millivolt"/>
      
	  <variable name="V" public_interface="in" units="millivolt"/>
      
	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_L</ci>
            <apply>
               <times/>
               <ci>g_L</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_L</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="i_Ca">
      <variable cmeta:id="i_Ca_i_Ca" name="i_Ca" public_interface="out" units="microA_per_cm2"/>
      
	  <variable initial_value="100.0" name="E_Ca" units="millivolt"/>
      <variable initial_value="4.0" name="g_Ca" units="milliS_per_cm2"/>
      <variable name="m_infinity" units="dimensionless"/>
	  <variable initial_value="10.0" name="V1" units="millivolt"/>
	  <variable initial_value="15.0" name="V2" units="millivolt"/>
      
	  <variable name="V" public_interface="in" units="millivolt"/>
     
	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_Ca</ci>
            <apply>
               <times/>
               <ci>g_Ca</ci>
               <ci>m_infinity</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_Ca</ci>
               </apply>
            </apply>
         </apply>
		 
		 <apply>
            <eq/>
            <ci>m_infinity</ci>
            <apply>
               <times/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.5</cn>
			   <apply>
                  <plus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <tanh/>
                     <apply>
                        <divide/>
						<apply>
                           <minus/>
						   <ci>V</ci>
						   <ci>V1</ci>
                        </apply>
						<ci>V2</ci>
                     </apply>
				  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="i_K">
      <variable cmeta:id="i_K_i_K" name="i_K" public_interface="out" units="microA_per_cm2"/>
      
	  <variable initial_value="-70.0" name="E_K" units="millivolt"/>
      <variable initial_value="8.0" name="g_K" units="milliS_per_cm2"/>
      
	  <variable name="time" private_interface="out" public_interface="in" units="millisecond"/>
      <variable name="V" private_interface="out" public_interface="in" units="millivolt"/>
      
	  <variable name="n" private_interface="in" units="dimensionless"/>
      
	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>i_K</ci>
            <apply>
               <times/>
               <ci>g_K</ci>
               <ci>n</ci>
               <apply>
                  <minus/>
                  <ci>V</ci>
                  <ci>E_K</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="i_K_n_gate">
      <variable initial_value="0.001" name="n" public_interface="out" units="dimensionless"/>
      
	  <variable name="n_infinity" units="dimensionless"/>
      <variable name="lambda_n" units="per_millisecond"/>
	  <variable initial_value="0.066666667" name="lambda_n_bar" units="per_millisecond"/>
      <variable initial_value="-1.0" name="V3" units="millivolt"/>
	  <variable initial_value="14.5" name="V4" units="millivolt"/>
	  
	  <variable name="V" public_interface="in" units="millivolt"/>
      <variable name="time" public_interface="in" units="millisecond"/>
      
	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>n_infinity</ci>
            <apply>
               <times/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0.5</cn>
			   <apply>
                  <plus/>
                  <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                  <apply>
                     <tanh/>
                     <apply>
                        <divide/>
						<apply>
                           <minus/>
						   <ci>V</ci>
						   <ci>V3</ci>
                        </apply>
						<ci>V4</ci>
                     </apply>
				  </apply>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <ci>lambda_n</ci>
            <apply>
               <times/>
               <ci>lambda_n_bar</ci>
               <apply>
                  <cosh/>
                  <apply>
                     <divide/>
				     <apply>
                        <minus/>
						<ci>V</ci>
						<ci>V3</ci>
                     </apply>
					 <apply>
                        <times/>
                        <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
						<ci>V4</ci>
                     </apply>               
				  </apply>
			   </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>n</ci>
            </apply>
            <apply>
               <times/>
               <ci>lambda_n</ci>
               <apply>
                  <minus/>
                  <ci>n_infinity</ci>
				  <ci>n</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>

   <group>
      <relationship_ref relationship="encapsulation"/>
      <component_ref component="i_K">
         <component_ref component="i_K_n_gate"/>
      </component_ref>
   </group>
   
   <connection>
      <map_components component_1="V" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   
   <connection>
      <map_components component_1="i_K" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>

   <connection>
      <map_components component_1="i_L" component_2="V"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_L" variable_2="i_L"/>
   </connection>
   
   <connection>
      <map_components component_1="i_Ca" component_2="V"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_Ca" variable_2="i_Ca"/>
   </connection>
   
   <connection>
      <map_components component_1="i_K" component_2="V"/>
      <map_variables variable_1="V" variable_2="V"/>
      <map_variables variable_1="i_K" variable_2="i_K"/>
   </connection>
   
   <connection>
      <map_components component_1="i_K" component_2="i_K_n_gate"/>
      <map_variables variable_1="n" variable_2="n"/>
      <map_variables variable_1="time" variable_2="time"/>
      <map_variables variable_1="V" variable_2="V"/>
   </connection>

   
   

<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#"><rdf:Description rdf:nodeID="n1"><rdf:value><rdf:Description rdf:nodeID="n2"/></rdf:value><subject_type xmlns="http://www.cellml.org/bqs/1.0#" xml:lang="en">keyword</subject_type></rdf:Description><rdf:Description rdf:nodeID="n3"><rdf:type><rdf:Description rdf:about="http://www.w3.org/1999/02/22-rdf-syntax-ns#Seq"/></rdf:type><rdf:_2><rdf:Description rdf:nodeID="n4"/></rdf:_2><rdf:_1><rdf:Description rdf:nodeID="n5"/></rdf:_1></rdf:Description><rdf:Description rdf:nodeID="n6"><creator xmlns="http://purl.org/dc/elements/1.1/"><rdf:Description rdf:nodeID="n3"/></creator><title xmlns="http://purl.org/dc/elements/1.1/" xml:lang="en">
            Voltage oscillations in the barnacle giant muscle fiber
          </title><issued xmlns="http://purl.org/dc/terms/"><rdf:Description rdf:nodeID="n7"/></issued><first_page xmlns="http://www.cellml.org/bqs/1.0#" xml:lang="en">193</first_page><volume xmlns="http://www.cellml.org/bqs/1.0#" xml:lang="en">35</volume><Journal xmlns="http://www.cellml.org/bqs/1.0#"><rdf:Description rdf:nodeID="n8"/></Journal><last_page xmlns="http://www.cellml.org/bqs/1.0#" xml:lang="en">213</last_page></rdf:Description><rdf:Description rdf:nodeID="n9"><ORG xmlns="http://www.w3.org/2001/vcard-rdf/3.0#"><rdf:Description rdf:nodeID="n10"/></ORG><N xmlns="http://www.w3.org/2001/vcard-rdf/3.0#"><rdf:Description rdf:nodeID="n11"/></N><EMAIL xmlns="http://www.w3.org/2001/vcard-rdf/3.0#"><rdf:Description rdf:nodeID="n12"/></EMAIL></rdf:Description><rdf:Description rdf:about="#morris_1981"><reference xmlns="http://www.cellml.org/bqs/1.0#"><rdf:Description rdf:nodeID="n13"/></reference><reference xmlns="http://www.cellml.org/bqs/1.0#"><rdf:Description rdf:nodeID="n14"/></reference><simulation xmlns="http://www.cellml.org/metadata/simulation/1.0#"><rdf:Description rdf:nodeID="n15"/></simulation></rdf:Description><rdf:Description rdf:nodeID="n16"><rdf:first><rdf:Description rdf:nodeID="n17"/></rdf:first><rdf:rest><rdf:Description rdf:about="http://www.w3.org/1999/02/22-rdf-syntax-ns#nil"/></rdf:rest></rdf:Description><rdf:Description rdf:nodeID="n18"><Family xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">Morris</Family><Given xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">C</Given></rdf:Description><rdf:Description rdf:about=""><creator xmlns="http://purl.org/dc/elements/1.1/"><rdf:Description rdf:nodeID="n9"/></creator><title xmlns="http://purl.org/dc/elements/1.1/" xml:lang="en">
            Voltage oscillations in the barnacle giant muscle fiber: reduced model
          </title></rdf:Description><rdf:Description rdf:nodeID="n19"><Family xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">Lecar</Family><Given xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">H</Given></rdf:Description><rdf:Description rdf:nodeID="n17"><pointDensity xmlns="http://www.cellml.org/metadata/simulation/1.0#nonstandard-" xml:lang="en">10000</pointDensity><maximumStepSize xmlns="http://www.cellml.org/metadata/simulation/1.0#" xml:lang="en">10</maximumStepSize><endingValue xmlns="http://www.cellml.org/metadata/simulation/1.0#" xml:lang="en">250</endingValue></rdf:Description><rdf:Description rdf:nodeID="n20"><N xmlns="http://www.w3.org/2001/vcard-rdf/3.0#"><rdf:Description rdf:nodeID="n19"/></N></rdf:Description><rdf:Description rdf:nodeID="n8"><title xmlns="http://purl.org/dc/elements/1.1/" xml:lang="en">Biophysical Journal</title></rdf:Description><rdf:Description rdf:nodeID="n10"><Orgunit xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">Auckland Bioengineering Institute</Orgunit><Orgname xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">The University of Auckland</Orgname></rdf:Description><rdf:Description rdf:nodeID="n11"><Family xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">Lloyd</Family><Given xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">Catherine</Given><Other xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">May</Other></rdf:Description><rdf:Description rdf:nodeID="n12"><rdf:type><rdf:Description rdf:about="http://imc.org/vCard/3.0#internet"/></rdf:type><rdf:value xml:lang="en">c.lloyd@auckland.ac.nz</rdf:value></rdf:Description><rdf:Description rdf:nodeID="n2"><rdf:type><rdf:Description rdf:about="http://www.w3.org/1999/02/22-rdf-syntax-ns#Bag"/></rdf:type><rdf:_2 xml:lang="en">muscle fibre</rdf:_2><rdf:_1 xml:lang="en">electrophysiology</rdf:_1></rdf:Description><rdf:Description rdf:nodeID="n13"><subject xmlns="http://purl.org/dc/elements/1.1/"><rdf:Description rdf:nodeID="n1"/></subject></rdf:Description><rdf:Description rdf:nodeID="n5"><Person xmlns="http://www.cellml.org/bqs/1.0#"><rdf:Description rdf:nodeID="n21"/></Person></rdf:Description><rdf:Description rdf:nodeID="n7"><W3CDTF xmlns="http://purl.org/dc/terms/" xml:lang="en">1981-07</W3CDTF></rdf:Description><rdf:Description rdf:nodeID="n4"><Person xmlns="http://www.cellml.org/bqs/1.0#"><rdf:Description rdf:nodeID="n20"/></Person></rdf:Description><rdf:Description rdf:nodeID="n14"><JournalArticle xmlns="http://www.cellml.org/bqs/1.0#"><rdf:Description rdf:nodeID="n6"/></JournalArticle><Pubmed_id xmlns="http://www.cellml.org/bqs/1.0#" xml:lang="en">7260316</Pubmed_id></rdf:Description><rdf:Description rdf:nodeID="n15"><boundIntervals xmlns="http://www.cellml.org/metadata/simulation/1.0#"><rdf:Description rdf:nodeID="n16"/></boundIntervals></rdf:Description><rdf:Description rdf:nodeID="n21"><N xmlns="http://www.w3.org/2001/vcard-rdf/3.0#"><rdf:Description rdf:nodeID="n18"/></N></rdf:Description></rdf:RDF></model>