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

<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
  <articleinfo>
  <title>A Mathematical Model for Chronic Myelogenous Leukemia (CML) and T Cell Interaction</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 model runs in OpenCell to recreate the published results. This particular version of the model recreates figure 6 from the paper (where CML decreases over a 2 year period).  The units have been checked and they are consistent. The CellML model also runs in COR however due to the timescale being in days the model is not ideal for simulation in COR.
          </para>
  </section>
  <sect1 id="sec_structure">
<title>Model Structure</title>

<para>
  ABSTRACT: In this paper, we propose and analyse a mathematical model for chronic myelogenous leukemia (CML), a cancer of the blood. We model the interaction between naive T cells, effector T cells, and CML cancer cells in the body, using a system of ordinary differential equations which gives rates of change of the three cell populations. One of the difficulties in modeling CML is the scarcity of experimental data which can be used to estimate parameters values. To compensate for the resulting uncertainties, we use Latin hypercube sampling (LHS) on large ranges of possible parameter values in our analysis. A major goal of this work is the determination of parameters which play a critical role in remission or clearance of the cancer in the model. Our analysis examines 12 parameters, and identifies two of these, the growth and death rates of CML, as critical to the outcome of the system. Our results indicate that the most promising research avenues for treatments of CML should be those that affect these two significant parameters (CML growth and death rates), while altering the other parameters should have little effect on the outcome.  
</para>

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

<para>
A mathematical model for chronic myelogenous leukemia (CML) and T cell interaction, Helen Moore and Natasha K. Li, 2004, <emphasis>Journal of Theoretical Biology</emphasis>, 11, 369-391. <ulink url="http://www.ncbi.nlm.nih.gov/pubmed/15038986">PubMed ID: 15038986</ulink>
</para>

<informalfigure float="0" id="fig_reaction_diagram">
          <mediaobject>
            <imageobject>
              <objectinfo>
                <title>reaction diagram</title>
              </objectinfo>
              <imagedata fileref="moore_2004.png"/>
            </imageobject>
            <caption>
Cell population diagram showing the population dynamics of each cell type and how they interact with each other.  These dynamics and interactions are described by the mathematical model equations.
</caption>
          </mediaobject>
        </informalfigure>

</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>
  
  <units base_units="yes" name="cells"/>
  
  <units name="cells_per_microlitre">
    <unit units="cells"/>
    <unit exponent="-1" prefix="micro" units="litre"/>
  </units>  
  
  <units name="flux">
    <unit units="cells"/>
    <unit exponent="-1" prefix="micro" units="litre"/>
    <unit exponent="-1" units="day"/>
  </units>
  
  <units name="microlitre_per_cells_day">
    <unit prefix="micro" units="litre"/>
    <unit exponent="-1" units="cells"/>
    <unit exponent="-1" units="day"/>
  </units>

  
  <component name="environment">
    <variable cmeta:id="environment_time" name="time" public_interface="out" units="day"/>
  </component>
  
  <component cmeta:id="Tn" name="Tn">
      
    <variable cmeta:id="Tn_Tn" initial_value="1510.0" name="Tn" public_interface="out" units="cells_per_microlitre"/>
    
    <variable initial_value="0.37" name="sn" units="flux"/>
    <variable initial_value="0.23" name="dn" units="first_order_rate_constant"/>
    
    <variable name="kn" public_interface="in" units="first_order_rate_constant"/>
    <variable name="eta" public_interface="in" units="cells_per_microlitre"/>
    <variable name="C" public_interface="in" units="cells_per_microlitre"/>
    <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>Tn</ci> 
        </apply>
        <apply>
          <minus/>
          <ci> sn </ci>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> dn </ci>
              <ci> Tn </ci>
            </apply>
            <apply>
              <times/>
              <ci> kn </ci>
              <ci> Tn </ci>
              <apply>
                <divide/>
                <ci> C </ci>
                <apply>
                  <plus/>
                  <ci> C </ci>
                  <ci> eta </ci>
                </apply>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component cmeta:id="Te" name="Te">
      
    <variable cmeta:id="Te_Te" initial_value="20.0" name="Te" public_interface="out" units="cells_per_microlitre"/>
    
    <variable initial_value="0.14" name="alpha_n" units="dimensionless"/>
    <variable initial_value="0.98" name="alpha_e" units="first_order_rate_constant"/>
    <variable initial_value="0.30" name="de" units="first_order_rate_constant"/>
    <variable initial_value="0.057" name="gamma_e" units="microlitre_per_cells_day"/>
    
    <variable name="kn" public_interface="in" units="first_order_rate_constant"/>
    <variable name="eta" public_interface="in" units="cells_per_microlitre"/>
    <variable name="C" public_interface="in" units="cells_per_microlitre"/>
    <variable name="Tn" public_interface="in" units="cells_per_microlitre"/>
    <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>Te</ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> alpha_n </ci>
              <ci> kn </ci>
              <ci> Tn </ci>
              <apply>
                <divide/>
                <ci> C </ci>
                <apply>
                  <plus/>
                  <ci> C </ci>
                  <ci> eta </ci>
                </apply>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> alpha_e </ci>
              <ci> Te </ci>
              <apply>
                <divide/>
                <ci> C </ci>
                <apply>
                  <plus/>
                  <ci> C </ci>
                  <ci> eta </ci>
                </apply>
              </apply>
            </apply>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> de </ci>
              <ci> Te </ci>
            </apply>
            <apply>
              <times/>
              <ci> gamma_e </ci>
              <ci> C </ci>
              <ci> Te </ci>
            </apply>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component cmeta:id="C" name="C">
      
    <variable cmeta:id="C_C" initial_value="10000.0" name="C" public_interface="out" units="cells_per_microlitre"/>
    
    <variable initial_value="230000" name="Cmax" units="cells_per_microlitre"/>
    <variable initial_value="0.0057" name="rc" units="first_order_rate_constant"/>
    <variable initial_value="0.024" name="dc" units="first_order_rate_constant"/>
    <variable initial_value="0.0034" name="gamma_c" units="microlitre_per_cells_day"/>
    
    <variable name="Te" public_interface="in" units="cells_per_microlitre"/>
    <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>C</ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> rc </ci>
            <ci> C </ci>
            <apply>
              <ln/>
              <apply>
                <divide/>
                <ci> Cmax </ci>
                <ci> C </ci>
              </apply>
            </apply>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> dc </ci>
              <ci> C </ci>
            </apply>
            <apply>
              <times/>
              <ci> gamma_c </ci>
              <ci> C </ci>
              <ci> Te </ci>
            </apply>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component name="model_parameters"> 
    <variable initial_value="0.062" name="kn" public_interface="out" units="first_order_rate_constant"/>
    <variable initial_value="720.0" name="eta" public_interface="out" units="cells_per_microlitre"/>
  </component>
  
  
  <connection>
    <map_components component_1="Tn" component_2="environment"/>
    <map_variables variable_1="time" variable_2="time"/>
  </connection>
  
  <connection>
    <map_components component_1="Te" component_2="environment"/>
    <map_variables variable_1="time" variable_2="time"/>
  </connection>
  
  <connection>
    <map_components component_1="C" component_2="environment"/>
    <map_variables variable_1="time" variable_2="time"/>
  </connection>
  
  <connection>
    <map_components component_1="model_parameters" component_2="Tn"/>
    <map_variables variable_1="kn" variable_2="kn"/>
    <map_variables variable_1="eta" variable_2="eta"/>
  </connection>
  
  <connection>
    <map_components component_1="model_parameters" component_2="Te"/>
    <map_variables variable_1="kn" variable_2="kn"/>
    <map_variables variable_1="eta" variable_2="eta"/>
  </connection>
  
  <connection>
    <map_components component_1="Tn" component_2="C"/>
    <map_variables variable_1="C" variable_2="C"/>
  </connection>
  
  <connection>
    <map_components component_1="Te" component_2="C"/>
    <map_variables variable_1="Te" variable_2="Te"/>
    <map_variables variable_1="C" variable_2="C"/>
  </connection>
  
  <connection>
    <map_components component_1="Te" component_2="Tn"/>
    <map_variables variable_1="Tn" variable_2="Tn"/>
  </connection>
  


<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#"><rdf:Description rdf:about="rdf:#273c0ad7-3920-40ec-b95c-d1bc272dfc2b"><W3CDTF xmlns="http://purl.org/dc/terms/" xml:lang="en">2004-04-21</W3CDTF></rdf:Description><rdf:Description rdf:nodeID="n1"><rdf:first><rdf:Description rdf:nodeID="n2"/></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:about="rdf:#19a39640-f856-4b14-8202-2059539fc516"><rdf:value><rdf:Description rdf:about="rdf:#2d21808a-1b91-4741-85c7-1f95235a667d"/></rdf:value><subject_type xmlns="http://www.cellml.org/bqs/1.0#" xml:lang="en">keyword</subject_type></rdf:Description><rdf:Description rdf:about="rdf:#4a94e2ab-7c70-4fb9-8e1c-f058ba3fa8f8"><title xmlns="http://purl.org/dc/elements/1.1/" xml:lang="en">Journal of Theoretical Biology</title></rdf:Description><rdf:Description rdf:nodeID="n2"><endingValue xmlns="http://www.cellml.org/metadata/simulation/1.0#" xml:lang="en">800</endingValue><pointDensity xmlns="http://www.cellml.org/metadata/simulation/1.0#nonstandard-" xml:lang="en">100000</pointDensity></rdf:Description><rdf:Description rdf:about="rdf:#ee7b2cab-b0ca-407f-8f40-7fbab6824c86"><rdf:type><rdf:Description rdf:about="http://www.cellml.org/bqs/1.0#Person"/></rdf:type><N xmlns="http://www.w3.org/2001/vcard-rdf/3.0#"><rdf:Description rdf:about="rdf:#2a0b1381-b6b5-4b08-9644-2c5ff91679de"/></N></rdf:Description><rdf:Description rdf:about=""><title xmlns="http://purl.org/dc/elements/1.1/" xml:lang="en">A Mathematical Model for Chronic Myelogenous Leukemia (CML) and T Cell Interaction (Decreasing CML Variant)</title><creator xmlns="http://purl.org/dc/elements/1.1/"><rdf:Description rdf:about="rdf:#9c4ff941-7a32-4cdd-889a-89e9b2b5f78f"/></creator><publisher xmlns="http://purl.org/dc/elements/1.1/" xml:lang="en">
        The University of Auckland, Bioengineering Institute
      </publisher><created xmlns="http://purl.org/dc/terms/"><rdf:Description rdf:about="rdf:#36b349de-8251-41b4-9b46-f651b1f32228"/></created></rdf:Description><rdf:Description rdf:about="#moore_2004"><title xmlns="http://purl.org/dc/elements/1.1/" xml:lang="en">
        Moore and Li's 2004 mathematical model for chronic myelogenous leukemia and T cell interaction.
      </title><reference xmlns="http://www.cellml.org/bqs/1.0#"><rdf:Description rdf:about="rdf:#4d88ecab-9800-4d91-a070-e022b200c3ae"/></reference><reference xmlns="http://www.cellml.org/bqs/1.0#"><rdf:Description rdf:about="rdf:#c96a07ad-0ffb-4f96-9589-883fc3d997df"/></reference><simulation xmlns="http://www.cellml.org/metadata/simulation/1.0#"><rdf:Description rdf:nodeID="n3"/></simulation><cmeta:bio_entity xml:lang="en">T lymphocyte</cmeta:bio_entity><cmeta:comment><rdf:Description rdf:about="rdf:#8649b4d6-f8a9-4736-94c6-e98872398295"/></cmeta:comment></rdf:Description><rdf:Description rdf:about="rdf:#9c4ff941-7a32-4cdd-889a-89e9b2b5f78f"><N xmlns="http://www.w3.org/2001/vcard-rdf/3.0#"><rdf:Description rdf:about="rdf:#b784e183-a2d7-4b9a-b535-56193bf324a1"/></N><EMAIL xmlns="http://www.w3.org/2001/vcard-rdf/3.0#"><rdf:Description rdf:about="rdf:#4bd57f3b-d972-415f-b1c5-dcabcc9078a8"/></EMAIL><ORG xmlns="http://www.w3.org/2001/vcard-rdf/3.0#"><rdf:Description rdf:about="rdf:#772c5b02-2381-443d-93f1-dc052b5d5003"/></ORG></rdf:Description><rdf:Description rdf:about="rdf:#772c5b02-2381-443d-93f1-dc052b5d5003"><Orgname xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">The University of Auckland</Orgname><Orgunit xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">The Bioengineering Institute</Orgunit></rdf:Description><rdf:Description rdf:about="rdf:#b784e183-a2d7-4b9a-b535-56193bf324a1"><Other xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">May</Other><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></rdf:Description><rdf:Description rdf:about="rdf:#67e7c4c0-4be2-458d-97f0-980f4e27eb19"><rdf:type><rdf:Description rdf:about="http://www.w3.org/1999/02/22-rdf-syntax-ns#Seq"/></rdf:type><rdf:_2><rdf:Description rdf:about="rdf:#8d45c880-8322-4699-a01d-b600eaf81253"/></rdf:_2><rdf:_1><rdf:Description rdf:about="rdf:#ee7b2cab-b0ca-407f-8f40-7fbab6824c86"/></rdf:_1></rdf:Description><rdf:Description rdf:about="rdf:#4bd57f3b-d972-415f-b1c5-dcabcc9078a8"><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:about="rdf:#c96a07ad-0ffb-4f96-9589-883fc3d997df"><Pubmed_id xmlns="http://www.cellml.org/bqs/1.0#" xml:lang="en">15038986</Pubmed_id><JournalArticle xmlns="http://www.cellml.org/bqs/1.0#"><rdf:Description rdf:about="rdf:#7d3ae7e2-66a9-4468-8433-2fd9368938d1"/></JournalArticle></rdf:Description><rdf:Description rdf:about="/Tn"><title xmlns="http://purl.org/dc/elements/1.1/" xml:lang="en">Tn</title><alternative xmlns="http://purl.org/dc/terms/" xml:lang="en">naive T cells</alternative></rdf:Description><rdf:Description rdf:nodeID="n3"><boundIntervals xmlns="http://www.cellml.org/metadata/simulation/1.0#"><rdf:Description rdf:nodeID="n1"/></boundIntervals></rdf:Description><rdf:Description rdf:about="rdf:#8649b4d6-f8a9-4736-94c6-e98872398295"><rdf:value xml:lang="en">
          This is the CellML description of Moore and Li's 2004 mathematical model for chronic myelogenous leukemia and T cell interaction.
        </rdf:value><creator xmlns="http://purl.org/dc/elements/1.1/"><rdf:Description rdf:about="rdf:#07facc3f-7e93-4f47-9b8f-9dc281783370"/></creator></rdf:Description><rdf:Description rdf:about="rdf:#4d88ecab-9800-4d91-a070-e022b200c3ae"><subject xmlns="http://purl.org/dc/elements/1.1/"><rdf:Description rdf:about="rdf:#19a39640-f856-4b14-8202-2059539fc516"/></subject></rdf:Description><rdf:Description rdf:about="rdf:#8d45c880-8322-4699-a01d-b600eaf81253"><rdf:type><rdf:Description rdf:about="http://www.cellml.org/bqs/1.0#Person"/></rdf:type><N xmlns="http://www.w3.org/2001/vcard-rdf/3.0#"><rdf:Description rdf:about="rdf:#a7e9d350-4382-45b6-bd3b-01256a0b0e18"/></N></rdf:Description><rdf:Description rdf:about="rdf:#36b349de-8251-41b4-9b46-f651b1f32228"><W3CDTF xmlns="http://purl.org/dc/terms/" xml:lang="en">2004-03-27</W3CDTF></rdf:Description><rdf:Description rdf:about="rdf:#2d21808a-1b91-4741-85c7-1f95235a667d"><rdf:type><rdf:Description rdf:about="http://www.w3.org/1999/02/22-rdf-syntax-ns#Bag"/></rdf:type><rdf:_3 xml:lang="en">t cell</rdf:_3><rdf:_4 xml:lang="en">cancer</rdf:_4><rdf:_2 xml:lang="en">immunology</rdf:_2><rdf:_1 xml:lang="en">T cell</rdf:_1></rdf:Description><rdf:Description rdf:about="rdf:#2a0b1381-b6b5-4b08-9644-2c5ff91679de"><Family xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">Moore</Family><Given xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">Helen</Given></rdf:Description><rdf:Description rdf:about="rdf:#07facc3f-7e93-4f47-9b8f-9dc281783370"><FN xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">Catherine Lloyd</FN></rdf:Description><rdf:Description rdf:about="/Te"><title xmlns="http://purl.org/dc/elements/1.1/" xml:lang="en">Te</title><alternative xmlns="http://purl.org/dc/terms/" xml:lang="en">effector T cells specific to CML</alternative></rdf:Description><rdf:Description rdf:about="rdf:#7d3ae7e2-66a9-4468-8433-2fd9368938d1"><title xmlns="http://purl.org/dc/elements/1.1/" xml:lang="en">
            A mathematical model for chronic myelogenous leukemia (CML) and T cell interaction 
          </title><issued xmlns="http://purl.org/dc/terms/"><rdf:Description rdf:about="rdf:#273c0ad7-3920-40ec-b95c-d1bc272dfc2b"/></issued><first_page xmlns="http://www.cellml.org/bqs/1.0#" xml:lang="en">513</first_page><Journal xmlns="http://www.cellml.org/bqs/1.0#"><rdf:Description rdf:about="rdf:#4a94e2ab-7c70-4fb9-8e1c-f058ba3fa8f8"/></Journal><last_page xmlns="http://www.cellml.org/bqs/1.0#" xml:lang="en">523</last_page><creator xmlns="http://purl.org/dc/elements/1.1/"><rdf:Description rdf:about="rdf:#67e7c4c0-4be2-458d-97f0-980f4e27eb19"/></creator><volume xmlns="http://www.cellml.org/bqs/1.0#" xml:lang="en">227</volume></rdf:Description><rdf:Description rdf:about="rdf:#a7e9d350-4382-45b6-bd3b-01256a0b0e18"><Other xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">K</Other><Family xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">Li</Family><Given xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">Natasha</Given></rdf:Description><rdf:Description rdf:about="/C"><title xmlns="http://purl.org/dc/elements/1.1/" xml:lang="en">C</title><alternative xmlns="http://purl.org/dc/terms/" xml:lang="en">chronic myelogenous leukemia (CML) cancer cells</alternative></rdf:Description></rdf:RDF></model>