Location: Friel, 1995 @ 852410459f1b / friel_1995.cellml

Author:
Catherine Lloyd <c.lloyd@auckland.ac.nz>
Date:
2010-07-10 04:21:39+12:00
Desc:
Extensively curated the model (essentially a re-write) but still can't get the CellML to replicate the published results.
Permanent Source URI:
https://staging.physiomeproject.org/workspace/friel_1995/rawfile/852410459f1bd982b0e4bfd31ccfc9c4de4d1a38/friel_1995.cellml

<?xml version='1.0' encoding='utf-8'?>
<!--  FILE :  friel_model_1995.xml

CREATED :  2nd April 2002

LAST MODIFIED : 9th April 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 D.D. Friel's 1995 model of Ca2+ oscillations in sympathetic neurons.

CHANGES:  
  22/07/2002 - CML - Added more metadata.   
  09/04/2003 - AAC - Added publication date information.  
-->
<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="friel_1995" name="friel_1995">
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
  <articleinfo>
  <title>[Ca2+]i oscillations in sympathetic neurons: an experimental test of a theoretical model</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 runs in OpenCell and COR but unfortuntately does not recreate the published figures. Specifically the model has been parameterised with the values provided in the caption of figure 4. However the CellML model does not spontaneously oscillate. This may be due to there being incorrect initial concentrations of intracellular and sub-space calcium, which are not specified in the paper.
          </para>
  </section>
  <sect1 id="sec_structure">
<title>Model Structure</title>

<para>
ABSTRACT: [Ca2+]i oscillations have been described in a variety of cells. This study focuses on caffeine-induced [Ca2+]i oscillations in sympathetic neurons. Previous work has shown that these oscillations require Ca2+ entry from the extracellular medium and Ca(2+)-induced Ca2+ release from a caffeine- and ryanodine-sensitive store. The aim of the study was to understand the mechanism responsible for the oscillations. As a starting point, [Ca2+]i relaxations were examined after membrane depolarization and exposure to caffeine. For both stimuli, post-stimulus relaxations could be described by the sum of two decaying exponential functions, consistent with a one-pool system in which Ca2+ transport between compartments is regulated by linear Ca2+ pumps and leaks. After modifying the store to include a [Ca2+]i-sensitive leak, the model also exhibits oscillations such as those observed experimentally. The model was tested by comparing measured and predicted net Ca2+ fluxes during the oscillatory cycle. Three independent fluxes were measured, describing the rates of 1) Ca2+ entry across the plasma membrane, 2) Ca2+ release by the internal store, and 3) Ca2+ extrusion across the plasma membrane and uptake by the internal store. Starting with estimates of the model parameters deduced from post-stimulus relaxations and the rapid upstroke, a set of parameter values was found that provides a good description of [Ca2+]i throughout the oscillatory cycle. With the same parameter values, there was also good agreement between the measured and simulated net fluxes. Thus, a one-pool model with a single [Ca2+]i-sensitive Ca2+ permeability is adequate to account for many of the quantitative properties of steady-state [Ca2+]i oscillations in sympathetic neurons. Inactivation of the intracellular Ca2+ permeability, cooperative nonlinear Ca2+ uptake and extrusion mechanisms, and functional links between plasma membrane Ca2+ transport and the internal store are not required.
</para>

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

<para>
[Ca2+]i oscillations in sympathetic neurons: an experimental test of a theoretical model, David D. Friel, 1995, <emphasis>Biophysical Journal</emphasis>, 68, 1752-1766.  <ulink url="http://www.ncbi.nlm.nih.gov/pubmed/7612818">PubMed ID: 7612818</ulink>
</para>


<informalfigure float="0" id="fig_cell_diagram">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>cell schematic for the model</title>
    </objectinfo>
    <imagedata fileref="friel_1995.png"/>
  </imageobject>
</mediaobject>
<caption>Schematic of the model indicating Ca<superscript>2+</superscript> compartmentalization in the extracellular matrix, cytosol and the mitochondrial matrix and pathways for Ca<superscript>2+</superscript> ion movement between the compartments.</caption>
</informalfigure>

</sect1>
</article>
</documentation>
  
  
  <units name="per_second">
    <unit units="second" exponent="-1"/>
  </units>
  
  <units name="nanomolar">
    <unit units="mole" prefix="nano"/>
    <unit units="litre" exponent="-1"/>
  </units>
  
  <units name="micromolar">
    <unit units="mole" prefix="micro"/>
    <unit units="litre" exponent="-1"/>
  </units>
  
  <units name="nanomolar_per_second">
    <unit units="nanomolar"/>
    <unit units="second" exponent="-1"/>
  </units>
  
  <units name="picolitre">
    <unit units="litre" prefix="pico"/>
  </units>
  
  <units name="picolitre_per_second">
    <unit units="picolitre"/>
	<unit units="second" exponent="-1"/>
  </units>
  
  <units name="picolitre_nanomolar_per_second">
    <unit units="picolitre"/>
	<unit units="nanomolar"/>
    <unit units="second" exponent="-1"/>
  </units>
  
  
  <component name="environment">
    <variable units="second" public_interface="out" name="time"/>
  </component>
  
  <component name="rate_constants">
    <variable units="per_second" public_interface="out" name="kappa_L1"/>
    <variable units="per_second" public_interface="out" name="kappa_P1" initial_value="0.132"/>
	<variable units="per_second" public_interface="out" name="kappa_L2"/>
	<variable units="per_second" public_interface="out" name="kappa_P2" initial_value="3.78"/>
    
	<variable units="per_second" name="kappa_L2_0" initial_value="0.054"/>
    <variable units="per_second" name="kappa_L2_1" initial_value="2.4"/>
    <variable units="nanomolar" name="Kd_Ca"/>
    <variable units="dimensionless" name="n" initial_value="3"/>
    
    <variable units="second" public_interface="in" name="time"/>
    <variable units="nanomolar" public_interface="in" name="Ca_i"/>
	
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="kappa_L1_calculation">
        <eq/>
        <ci> kappa_L1 </ci>
        <piecewise>
          <piece>
            <cn type="e-notation" cellml:units="per_second">5<sep/>-6</cn>
			<apply>
              <and/>
			  <apply>
                <geq/>
                <ci> time </ci>
                <cn cellml:units="second"> 0 </cn>
              </apply>
			  <apply>
                <lt/>
                <ci> time </ci>
                <cn cellml:units="second"> 30 </cn>
              </apply>
			</apply>
          </piece>
          <otherwise>
            <cn type="e-notation" cellml:units="per_second">20<sep/>-6</cn>
          </otherwise>
        </piecewise>
	  </apply>
	  
	  <apply id="kappa_L2_calculation">
        <eq/>
        <ci> kappa_L2 </ci>
        <apply>
          <divide/>
          <apply>
            <plus/>
            <ci> kappa_L2_0 </ci>
            <ci> kappa_L2_1 </ci>
          </apply>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
            <apply>
              <power/>
              <apply>
                <divide/>
                <ci> Kd_Ca </ci>
                <ci> Ca_i </ci>
              </apply>
              <ci> n </ci>
            </apply>
          </apply>
        </apply>
      </apply>
	  
	  <apply id="Kd_Ca_calculation">
        <eq/>
        <ci> Kd_Ca </ci>
        <piecewise>
          <piece>
            <cn cellml:units="micromolar"> 1.0 </cn>
            <apply>
              <and/>
			  <apply>
                <geq/>
                <ci> time </ci>
                <cn cellml:units="second"> 0 </cn>
              </apply>
			  <apply>
                <lt/>
                <ci> time </ci>
                <cn cellml:units="second"> 60 </cn>
              </apply>
			</apply>
          </piece>
          <otherwise>
            <cn cellml:units="micromolar"> 0.5 </cn>
          </otherwise>
        </piecewise>
	  </apply>
    </math>
  </component>
  
  <component name="cytosolic_calcium">
    <variable units="nanomolar" public_interface="out" name="Ca_i" initial_value="100.0"/>
    <variable units="nanomolar" public_interface="out" name="Ca_o" initial_value="1000000.0"/>
    <variable units="nanomolar" public_interface="out" name="Ca_i_ss"/>
    <variable units="dimensionless" name="gamma" public_interface="out" initial_value="0.24"/>
    
    <variable units="per_second" public_interface="in" name="kappa_L1"/>
	<variable units="per_second" public_interface="in" name="kappa_P1"/>
	<variable units="per_second" public_interface="in" name="kappa_L2"/>
	<variable units="per_second" public_interface="in" name="kappa_P2"/>
    <variable units="nanomolar" public_interface="in" name="Ca_s"/>
    <variable units="second" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="Ca_i_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Ca_i </ci>
        </apply>
        <apply>
          <plus/>
          <apply>
            <times/>
            <apply>
              <minus/>
              <apply>
                <plus/>
                <ci> kappa_L1 </ci>
                <ci> kappa_P1 </ci>
                <apply>
                  <times/>
                  <ci> gamma </ci>
                  <apply>
                    <plus/>
                    <ci> kappa_L2 </ci>
                    <ci> kappa_P2 </ci>
                  </apply>
                </apply>
              </apply>
            </apply>
            <ci> Ca_i </ci>
          </apply>
          <apply>
            <times/>
            <ci> gamma </ci>
            <ci> kappa_L2 </ci>
            <ci> Ca_s </ci> 
          </apply>
          <apply>
            <times/> 
            <ci> kappa_L1 </ci>
            <ci> Ca_o </ci> 
          </apply>
        </apply>
      </apply>
	    
      <apply id="Ca_i_ss_calculation">
        <eq/>
        <ci> Ca_i_ss </ci>
        <apply>
          <divide/>
          <ci> Ca_o </ci>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
			<apply>
              <divide/>
              <ci> kappa_P1 </ci>
              <ci> kappa_L1 </ci>
			</apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <component name="subspace_calcium">
    <variable units="nanomolar" public_interface="out" name="Ca_s" initial_value="10.0"/>
    <variable units="nanomolar" public_interface="out" name="Ca_s_ss"/>
    
	<variable units="nanomolar" public_interface="in" name="Ca_i"/>
    <variable units="nanomolar" public_interface="in" name="Ca_i_ss"/>
    <variable units="per_second" public_interface="in" name="kappa_L2"/>
    <variable units="per_second" public_interface="in" name="kappa_P2"/>
    <variable units="second" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML">
      <apply id="Ca_s_diff_eq">
        <eq/>
        <apply>
          <diff/>
          <bvar>
            <ci> time </ci>
          </bvar>
          <ci> Ca_s </ci>
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <apply>
              <plus/>
              <ci> kappa_L2 </ci>
              <ci> kappa_P2 </ci>
            </apply>
            <ci> Ca_i </ci>
          </apply>
          <apply>
            <times/>
            <ci> kappa_L2 </ci>
            <ci> Ca_s </ci>
          </apply>
        </apply>
      </apply>
      
      <apply id="Ca_s_ss_calculation">
        <eq/>
        <ci> Ca_s_ss </ci>
        <apply>
          <times/>
          <ci> Ca_i_ss </ci>
          <apply>
            <plus/>
            <cn cellml:units="dimensionless"> 1.0 </cn>
			<apply>
              <divide/>
              <ci> kappa_P2 </ci>
              <ci> kappa_L2 </ci>
			</apply>
          </apply>
        </apply>
      </apply>
    </math>
  </component>
  
  <connection>
    <map_components component_2="environment" component_1="rate_constants"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="cytosolic_calcium"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="subspace_calcium"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="rate_constants" component_1="cytosolic_calcium"/>
    <map_variables variable_2="kappa_L1" variable_1="kappa_L1"/>
	<map_variables variable_2="kappa_L2" variable_1="kappa_L2"/>
	<map_variables variable_2="kappa_P1" variable_1="kappa_P1"/>
	<map_variables variable_2="kappa_P2" variable_1="kappa_P2"/>
    <map_variables variable_2="Ca_i" variable_1="Ca_i"/>
  </connection>
  
  <connection>
    <map_components component_2="rate_constants" component_1="subspace_calcium"/>
    <map_variables variable_2="kappa_L2" variable_1="kappa_L2"/>
	<map_variables variable_2="kappa_P2" variable_1="kappa_P2"/>
  </connection>
  
  <connection>
    <map_components component_2="subspace_calcium" component_1="cytosolic_calcium"/>
    <map_variables variable_2="Ca_i" variable_1="Ca_i"/>
    <map_variables variable_2="Ca_i_ss" variable_1="Ca_i_ss"/>
    <map_variables variable_2="Ca_s" variable_1="Ca_s"/>
  </connection>


<rdf:RDF>
  <rdf:Bag rdf:about="rdf:#38593c84-2b8a-4e76-9f16-2450acbdce46">
    <rdf:li>calcium dynamics</rdf:li>
    <rdf:li>sympathetic</rdf:li>
    <rdf:li>electrophysiology</rdf:li>
    <rdf:li>neuron</rdf:li>
	<rdf:li>neurobiology</rdf:li>
  </rdf:Bag>
  <rdf:Description rdf:about="rdf:#a523d839-d741-4ac6-8d0d-a940a448694b">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#67948363-0da2-4b6d-9fcc-0d233741ddd2">
    <bqs:Person rdf:resource="rdf:#0d6247e7-8a49-4371-a91d-0f1b62e4519b"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#783f94d7-a602-4efe-a009-70789daf65c9">
    <dcterms:modified rdf:resource="rdf:#097b3f66-46c3-4f5d-8533-55b2f977180f"/>
    <rdf:value>
          Added publication date information.
        </rdf:value>
    <cmeta:modifier rdf:resource="rdf:#91828417-0569-4499-b066-fc8c6af1379b"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#370c9a52-4b33-4c61-a3d0-795baa44bf25">
    <dcterms:modified rdf:resource="rdf:#330f494b-128e-4382-b84f-7063551353bb"/>
    <rdf:value>
          Added more metadata.
        </rdf:value>
    <cmeta:modifier rdf:resource="rdf:#f360b3e3-ad5a-4519-b3cf-c28fd8addb44"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#364d5d63-5b8f-4d5b-b27d-874df4234b04">
    <dc:creator rdf:resource="rdf:#67948363-0da2-4b6d-9fcc-0d233741ddd2"/>
    <dc:title>
            [Ca2+]i oscillations in sympathetic neurons: an experimental test of a theoretical model
          </dc:title>
    <bqs:volume>68</bqs:volume>
    <bqs:first_page>1752</bqs:first_page>
    <bqs:Journal rdf:resource="rdf:#d0e39266-dde4-4a8a-81c5-66d11b4e003e"/>
    <dcterms:issued rdf:resource="rdf:#1b72924c-01df-4725-b754-69fa4f3496d8"/>
    <bqs:last_page>1766</bqs:last_page>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#209e8eea-1e47-4185-b87d-d3704e8d95f7">
    <vCard:Given>David</vCard:Given>
    <vCard:Family>Friel</vCard:Family>
    <vCard:Other>D</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#59e4443d-4c3b-40eb-a8ab-b9029bd9ea90">
    <dcterms:W3CDTF>2002-04-02</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#1b72924c-01df-4725-b754-69fa4f3496d8">
    <dcterms:W3CDTF>1995-05</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#404066de-08f8-46f1-937e-9d554b7d5480">
    <bqs:Pubmed_id>7612818</bqs:Pubmed_id>
    <bqs:JournalArticle rdf:resource="rdf:#364d5d63-5b8f-4d5b-b27d-874df4234b04"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5f1a4948-db95-4dbd-bdd6-85879e1d13f4">
    <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:#097b3f66-46c3-4f5d-8533-55b2f977180f">
    <dcterms:W3CDTF>2003-04-09</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="#friel_1995">
    <dc:title>
        A Model Of Calcium Oscillations in Sympathetic Neurons 
      </dc:title>
    <cmeta:bio_entity>Sympathetic Neuron</cmeta:bio_entity>
    <cmeta:comment rdf:resource="rdf:#c0ab5e71-5262-443b-aba7-7db714a2affe"/>
    <bqs:reference rdf:resource="rdf:#404066de-08f8-46f1-937e-9d554b7d5480"/>
    <bqs:reference rdf:resource="rdf:#9d759697-4299-4616-b8fb-d7c18b622781"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#0d6247e7-8a49-4371-a91d-0f1b62e4519b">
    <vCard:N rdf:resource="rdf:#209e8eea-1e47-4185-b87d-d3704e8d95f7"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#113e461c-7f67-4e70-b08f-998352b084f3">
    <vCard:ORG rdf:resource="rdf:#cf2f6d31-4136-41fb-9590-2d8ada1b579c"/>
    <vCard:EMAIL rdf:resource="rdf:#5f1a4948-db95-4dbd-bdd6-85879e1d13f4"/>
    <vCard:N rdf:resource="rdf:#a523d839-d741-4ac6-8d0d-a940a448694b"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#aac2b6ce-2e80-4cb5-b672-8a1c263e01e4">
    <vCard:Given>Autumn</vCard:Given>
    <vCard:Family>Cuellar</vCard:Family>
    <vCard:Other>A</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#c0ab5e71-5262-443b-aba7-7db714a2affe">
    <dc:creator rdf:resource="rdf:#1442df3c-4a09-4214-9c05-42c0a1765b54"/>
    <rdf:value>
          This is the CellML description of David Friel's 1995 model of calcium 
          oscillations in sympathetic neurons.
        </rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#f360b3e3-ad5a-4519-b3cf-c28fd8addb44">
    <vCard:N rdf:resource="rdf:#acac6a63-a306-4fe0-beb7-c549e7998a2d"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#9d759697-4299-4616-b8fb-d7c18b622781">
    <dc:subject rdf:resource="rdf:#fdb08dac-74db-4f52-9895-16f933ba43e5"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#330f494b-128e-4382-b84f-7063551353bb">
    <dcterms:W3CDTF>2002-07-22</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d0e39266-dde4-4a8a-81c5-66d11b4e003e">
    <dc:title>Biophysical Journal</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#fdb08dac-74db-4f52-9895-16f933ba43e5">
    <bqs:subject_type>keyword</bqs:subject_type>
    <rdf:value rdf:resource="rdf:#38593c84-2b8a-4e76-9f16-2450acbdce46"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#91828417-0569-4499-b066-fc8c6af1379b">
    <vCard:N rdf:resource="rdf:#aac2b6ce-2e80-4cb5-b672-8a1c263e01e4"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#cf2f6d31-4136-41fb-9590-2d8ada1b579c">
    <vCard:Orgname>The University of Auckland</vCard:Orgname>
    <vCard:Orgunit>Auckland Bioengineering Institute</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#acac6a63-a306-4fe0-beb7-c549e7998a2d">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dc:publisher>
        The University of Auckland, Auckland Bioengineering Institute
      </dc:publisher>
    <cmeta:modification rdf:resource="rdf:#370c9a52-4b33-4c61-a3d0-795baa44bf25"/>
    <cmeta:modification rdf:resource="rdf:#783f94d7-a602-4efe-a009-70789daf65c9"/>
    <dcterms:created rdf:resource="rdf:#59e4443d-4c3b-40eb-a8ab-b9029bd9ea90"/>
    <dc:creator rdf:resource="rdf:#113e461c-7f67-4e70-b08f-998352b084f3"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#1442df3c-4a09-4214-9c05-42c0a1765b54">
    <vCard:FN>Catherine Lloyd</vCard:FN>
  </rdf:Description>
</rdf:RDF>
</model>