Location: Alexander and Wahl 2010 @ 7045a36c3bd8 / alexander_2010_system2.cellml

Author:
Hanne Nielsen <hnie010@aucklanduni.ac.nz>
Date:
2011-01-17 11:40:44+13:00
Desc:
Added xul and session file, need to fix zoom
Permanent Source URI:
https://models.cellml.org/workspace/2d8/rawfile/7045a36c3bd85666e4b738243aa195a95dabccff/alexander_2010_system2.cellml

<?xml version="1.0"?>

<model xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" name="alexander_2010" cmeta:id="alexander_2010" xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cellml="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#">

<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
  <articleinfo>
  <title>Self-tolerance and Autoimmunity in a Regulatory T Cell 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 represents system 2 (equations 2a-2d) in the original publication.  The model runs in both COR and OpenCell to replicate the published results (figure 2) with R0=0.8.  The units have been checked and they are consistent.
          </para>
  </section>
  <sect1 id="sec_structure">
<title>Model Structure</title>

<para>
ABSTRACT: The class of immunosuppressive lymphocytes known as regulatory T cells (Tregs) has been identified as a key component in preventing autoimmune diseases. Although Tregs have been incorporated previously in mathematical models of autoimmunity, we take a novel approach which emphasizes the importance of professional antigen presenting cells (pAPCs). We examine three possible mechanisms of Treg action (each in isolation) through ordinary differential equation (ODE) models. The immune response against a particular autoantigen is suppressed both by Tregs specific for that antigen and by Tregs of arbitrary specificities, through their action on either maturing or already mature pAPCs or on autoreactive effector T cells. In this deterministic approach, we find that qualitative long-term behaviour is predicted by the basic reproductive ratio R(0) for each system. When R(0) is less tHAN 1, only the trivial equilibrium exists and is stable; when R(0) is greater than 1, this equilibrium loses its stability and a stable non-trivial equilibrium appears. We interpret the absence of self-damaging populations at the trivial equilibrium to imply a state of self-tolerance, and their presence at the non-trivial equilibrium to imply a state of chronic autoimmunity. Irrespective of mechanism, our model predicts that Tregs specific for the autoantigen in question play no role in the system's qualitative long-term behaviour, but have quantitative effects that could potentially reduce an autoimmune response to sub-clinical levels. Our results also suggest an important role for Tregs of arbitrary specificities in modulating the qualitative outcome. A stochastic treatment of the same model demonstrates that the probability of developing a chronic autoimmune response increases with the initial exposure to self antigen or autoreactive effector T cells. The three different mechanisms we consider, while leading to a number of similar predictions, also exhibit key differences in both transient dynamics (ODE approach) and the probability of chronic autoimmunity (stochastic approach).
</para>

<informalfigure float="0" id="fig_reaction_diagram">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>model diagram</title>
    </objectinfo>
    <imagedata fileref="alexander_2010.png"/>
  </imageobject>
</mediaobject>
<caption>Flow chart illustrating interactions among populations and flow in/out of compartments in System 1. Populations to be modelled explicitly are in black boxes; those that are considered only implicitly or as intermediaries are in grey boxes. Movement in or out of a compartment is indicated by a black arrow; activating influences are indicated with green arrows; suppressive influences are indicated with red arrows; and other interactions or effects are indicated with blue arrows.</caption>
</informalfigure>

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

<para>
Self-tolerance and Autoimmunity in a Regulatory T Cell Model, Alexander HK and Wahl LM, 2010, <emphasis>Bulletin of Mathematical Biology</emphasis>.  <ulink url="http://www.ncbi.nlm.nih.gov/pubmed/20195912">PubMed ID: 20195912</ulink></para>

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



  
  <units name="day">
    <unit multiplier="86400.0" units="second"/>
  </units>
  
  <units name="first_order_rate_constant">
    <unit exponent="-1" units="day"/>
  </units>
  
  
  <component name="environment">
    <variable name="time" public_interface="out" units="day"/>
  </component>
  
  <component cmeta:id="A" name="A">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="#A">
        <dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">A</dc:title>
        <dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">mature pAPCs</dcterms:alternative> 
		<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">mature antigen presenting cells</dcterms:alternative> 
		<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">dendritic cells</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF>  
    <variable initial_value="1.0" name="A" public_interface="out" units="dimensionless"/>
    <variable name="v" public_interface="out" units="first_order_rate_constant" initial_value="1.25e5"/>
	<variable name="k" public_interface="out" units="dimensionless" initial_value="5e7"/>
    
    <variable name="f" units="dimensionless" initial_value="1e-4"/>
    <variable name="sigma1" units="first_order_rate_constant" initial_value="3e-6"/>
    <variable name="b1" units="first_order_rate_constant" initial_value="0.25"/>
	<variable name="muA" units="first_order_rate_constant" initial_value="0.25"/>
	
	<variable name="G" public_interface="in" units="dimensionless"/>
    <variable name="R" public_interface="in" units="dimensionless"/>
    <variable name="time" public_interface="in" units="day"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci> time </ci>
          </bvar> 
          <ci> A </ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <divide/>
			<apply>
              <times/>
              <ci> f </ci>
              <ci> v </ci>
              <ci> G </ci>
            </apply>
			<apply>
              <plus/>
              <ci> k </ci>
              <ci> G </ci>
            </apply>
		  </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <apply>
                <plus/>
                <apply>
                  <times/>
				  <ci> sigma1 </ci>
                  <ci> R </ci>
				</apply> 
			    <ci> b1 </ci>
		      </apply> 
			  <ci> A </ci>
		    </apply> 
            <apply>
              <times/>
              <ci> muA </ci> 
              <ci> A </ci>
            </apply> 
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component cmeta:id="R" name="R">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="#R">
        <dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">R</dc:title>
        <dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">active regulatory T cells</dcterms:alternative> 
		<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">regulatory T cells</dcterms:alternative> 
      </rdf:Description>
    </rdf:RDF>  
    <variable initial_value="0.0" name="R" public_interface="out" units="dimensionless"/>
    
    <variable name="pi1" units="first_order_rate_constant" initial_value="0.016"/>
    <variable name="beta" units="first_order_rate_constant" initial_value="200.0"/>
    <variable name="muR" units="first_order_rate_constant" initial_value="0.25"/>
    
	<variable name="A" public_interface="in" units="dimensionless"/>
    <variable name="E" public_interface="in" units="dimensionless"/>
    <variable name="time" public_interface="in" units="day"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>R</ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <apply>
              <plus/>
              <apply>
                <times/>
                <ci> pi1 </ci> 
                <ci> E </ci>
              </apply> 
			  <ci> beta </ci>
            </apply>
			<ci> A </ci>
          </apply>
          <apply>
            <times/>
            <ci> muR </ci>
            <ci> R </ci>
          </apply> 
        </apply>
      </apply>    
    </math>
  </component>
  
  <component cmeta:id="E" name="E">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="#E">
        <dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">E</dc:title>
        <dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">effector T cells</dcterms:alternative> 
      </rdf:Description>
    </rdf:RDF>  
    <variable initial_value="0.0" name="E" public_interface="out" units="dimensionless"/>
    
    <variable name="lambdaE" units="first_order_rate_constant" initial_value="1000.0"/>
    <variable name="muE" units="first_order_rate_constant" initial_value="0.25"/>
    
	<variable name="A" public_interface="in" units="dimensionless"/>
    <variable name="time" public_interface="in" units="day"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>E</ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> lambdaE </ci> 
            <ci> A </ci>
          </apply>
          <apply>
            <times/>
            <ci> muE </ci>
            <ci> E </ci>
          </apply> 
        </apply>
      </apply>    
    </math>
  </component>
  
  <component cmeta:id="G" name="G">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="#G">
        <dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">G</dc:title>
        <dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">antigen</dcterms:alternative> 
	  </rdf:Description>
    </rdf:RDF>  
    <variable initial_value="1e8" name="G" public_interface="out" units="dimensionless"/>
    
    <variable name="gamma" units="first_order_rate_constant" initial_value="2000.0"/>
    <variable name="muG" units="first_order_rate_constant" initial_value="5.0"/>
    
	<variable name="v" public_interface="in" units="first_order_rate_constant"/>
	<variable name="k" public_interface="in" units="dimensionless"/>
	<variable name="E" public_interface="in" units="dimensionless"/>
    <variable name="time" public_interface="in" units="day"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>G</ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> gamma </ci> 
            <ci> E </ci>
          </apply>
          <apply>
            <plus/>		  
			<apply>
              <divide/>
			  <apply>
                <times/>
                <ci> v </ci>
				<ci> G </ci>
			  </apply>
              <apply>
                <plus/>
                <ci> k </ci>
				<ci> G </ci>
			  </apply>
            </apply> 
			<apply>
              <times/>
              <ci> muG </ci>
              <ci> G </ci>
            </apply> 
		  </apply>
        </apply>
      </apply>    
    </math>
  </component>

  
  <connection>
    <map_components component_1="A" component_2="environment"/>
    <map_variables variable_1="time" variable_2="time"/>
  </connection>
  
  <connection>
    <map_components component_1="R" component_2="environment"/>
    <map_variables variable_1="time" variable_2="time"/>
  </connection>
  
  <connection>
    <map_components component_1="E" component_2="environment"/>
    <map_variables variable_1="time" variable_2="time"/>
  </connection>
  
  <connection>
    <map_components component_1="G" component_2="environment"/>
    <map_variables variable_1="time" variable_2="time"/>
  </connection>
 
  <connection>
    <map_components component_1="A" component_2="R"/>
    <map_variables variable_1="R" variable_2="R"/>
    <map_variables variable_1="A" variable_2="A"/>
  </connection>
  
  <connection>
    <map_components component_1="A" component_2="G"/>
    <map_variables variable_1="G" variable_2="G"/>
    <map_variables variable_1="v" variable_2="v"/>
	<map_variables variable_1="k" variable_2="k"/>
  </connection>
  
  <connection>
    <map_components component_1="R" component_2="E"/>
    <map_variables variable_1="E" variable_2="E"/>
  </connection>
  
  <connection>
    <map_components component_1="E" component_2="A"/>
    <map_variables variable_1="A" variable_2="A"/>
  </connection>
  
  <connection>
    <map_components component_1="G" component_2="E"/>
    <map_variables variable_1="E" variable_2="E"/>
  </connection>
  
  
  
  <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">
    
	<rdf:Description rdf:about="">
      <dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">
        Self-tolerance and Autoimmunity in a Regulatory T Cell Model: System 2
      </dc:title>
	  <dc:creator rdf:parseType="Resource">
        <vCard:N rdf:parseType="Resource">
          <vCard:Family>Lloyd</vCard:Family>
          <vCard:Given>Catherine</vCard:Given>
		  <vCard:Other>May</vCard:Other>
        </vCard:N>
        <vCard:EMAIL rdf:parseType="Resource">
          <rdf:value>c.lloyd@auckland.ac.nz</rdf:value>
          <rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
        </vCard:EMAIL>
        <vCard:ORG rdf:parseType="Resource">
          <vCard:Orgname>The University of Auckland</vCard:Orgname>
          <vCard:Orgunit>Auckland Bioengineering Institute</vCard:Orgunit>
        </vCard:ORG>
      </dc:creator>
    </rdf:Description>
    
    <rdf:Description rdf:about="#alexander_2010">
      <bqs:reference rdf:parseType="Resource">
        <dc:subject rdf:parseType="Resource">
          <bqs:subject_type>keyword</bqs:subject_type>
          <rdf:value>
            <rdf:Bag>
              <rdf:li>immunology</rdf:li>
			  <rdf:li>autoimmune</rdf:li>
			  <rdf:li>regulatory T cell</rdf:li>
            </rdf:Bag>
          </rdf:value>
        </dc:subject>
      </bqs:reference>

      <bqs:reference rdf:parseType="Resource">
        <bqs:Pubmed_id>20195912</bqs:Pubmed_id>
        <bqs:JournalArticle rdf:parseType="Resource">
          <dc:creator>
            <rdf:Seq>
              <rdf:li rdf:parseType="Resource">
                <bqs:Person rdf:parseType="Resource">
                  <vCard:N rdf:parseType="Resource">
                    <vCard:Family>Alexander</vCard:Family>
                    <vCard:Given>H</vCard:Given>
					<vCard:Other>K</vCard:Other>
                  </vCard:N>
                </bqs:Person>
              </rdf:li>
			  <rdf:li rdf:parseType="Resource">
                <bqs:Person rdf:parseType="Resource">
                  <vCard:N rdf:parseType="Resource">
                    <vCard:Family>Wahl</vCard:Family>
                    <vCard:Given>L</vCard:Given>
					<vCard:Other>M</vCard:Other>
                  </vCard:N>
                </bqs:Person>
              </rdf:li>
            </rdf:Seq>
          </dc:creator>
          <dc:title>
            Self-tolerance and Autoimmunity in a Regulatory T Cell Model
          </dc:title>
          <dcterms:issued rdf:parseType="Resource">
            <dcterms:W3CDTF>2010-03-02</dcterms:W3CDTF>
          </dcterms:issued>
          <bqs:Journal rdf:parseType="Resource">
            <dc:title>Bulletin of Mathematical Biology</dc:title>
          </bqs:Journal>
		</bqs:JournalArticle>
      </bqs:reference>
    </rdf:Description>
  </rdf:RDF>
</model>