Location: Waugh, Sherratt, 2007 @ f93224810c0b / waugh_sherratt_2007_a.cellml

Author:
Catherine Lloyd <c.lloyd@auckland.ac.nz>
Date:
2009-11-17 15:28:47+13:00
Desc:
Removed links from the citation and added documentation ready to link to Hanne's image.
Permanent Source URI:
https://models.physiomeproject.org/workspace/waugh_sherratt_2007/rawfile/f93224810c0bf5dfa4b947c3178ab82c133bd631/waugh_sherratt_2007_a.cellml

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

CREATED :  8th November 2007

LAST MODIFIED : 8th November 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 Waugh and Sherratt's 2006 mathematical model of macrophage dynamics in diabetic wound healing.

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

<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
  <articleinfo>
  <title>Modeling the effects of treating diabetic wounds with engineered skin substitutes</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 to recreate the results in the published paper.  The model has also been checked in COR and the units are consistent.  This particular version of the model represents wound healing in a healthy individual.
          </para>
  </section>
  <sect1 id="sec_structure">
  <title>Model Structure</title>

<para>
Wounds in diabetic patients can take much longer to heal than they would in a healthy individual; some taking up to 18 months to heal properly.  This suggests the normal wound healing process is disrupted in some way in diabetic patients, but despite intensive research, the mechanisms underlying this process are still relatively poorly understood.  However, by piecing together the information which has been obtained, it is possible to develop a simple mathematical model of the wound healing process and how it is influenced by diabetes.  
</para>

<para>
The first stage of the wound healing process is inflammation, and macrophages are some of the first cells to be recruited to the site of injury.  They are derived from monocytes and they differentiate into one of three types depending on the chemical stimulus they receive, namely; cytocidal (or killer), inflammatory, and repair macrophages.  Killer macrophages remove bacteria and other debris from the wound site by phagocytosis, inflammatory macrophages secrete cytokines and chemokines to attract fibroblasts and endothelial cells to the wound and encourage their proliferation, while repair macrophages help to remodel the extracellular matrix of the wound.  The ratio of inflammatory to repair macrophages seems to be a key determinant of the rate of wound healing. 
</para>

<para>
In diabetic patients cell proliferation is often impaired, endothelial cells are more likely to undergo apoptosis, blood glucose levels are elevated, blood vessel growth is impaired, and there can be decreased collagen deposition at the wound site.  Macrophage removal to the lymph nodes may also be impaired, which might explain the presence of macrophages at the wound site long after the inflammatory stage when they would have normally been cleared in a healthy individual. 
</para>

<para>
To date, mathematical models of wound healing have tended to focus on the cell proliferation and repair stages of the wound process.  However, in diabetic patients it seems the inflammatory phase should be modelled instead.  In the paper described here, Helen Waugh and Jonathan Sherratt address this issue by developing a basic mathematical model of the inflammatory stage of the wound healing process, in particular focusing on the behaviour of the macrophage populations.  This model is an extension of their previous 2006 model (also described in CellML) and in addition to investigating the kinetics of inflammatory macrophages, repair macrophages and TGF-beta, the current model also describes changing PDGF, collagen density, hyaluronan concentration and fibroblast dynamics.  These variables were chosen, in particular, because they are components of two commercially available engineering skin substitutes; Apligraf and Dermagraft, which can be used in diabetic wound healing, and are considered in the current 2007 mathematical model described here. 
</para>

<para>
Apligraf is an artificial skin comprising a dermal layer and an epidermal layer of cells seeded onto a bioabsorbable scaffold.  Although Apligraf possesses an epidermal layer of keratinocytes the mathematical model discussed here focuses on the dermal components of the therapy - namely neonatal fibroblasts, TGF-beta, PDGF, collagen and hyaluronan.  Similarly Dermagraft comprises fibroblasts seeded onto a Vicryl scaffold and it also contains TGF-beta, PDGF, collagen and hyaluronan.  However unlike Apligraf, Dermagraft does not possess an epidermal layer of keratinocytes, and the relative density of fibroblasts and concentrations of growth factors and fibres vary slightly. 
</para>

<para>
Model simulations suggest that the key component in the successful healing of diabetic wounds appear to be hyaluronan concentration and effect wound healing therapies work by increasing the amount of hyaluronan in the wound environment.  Model simulation data agree with time-to-healing results observed in clinical trials and the model assists in our understanding of why diabetic wounds do not heal and how treatments promote wound closure.
</para>

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

<para>
Modeling the effects of treating diabetic wounds with engineered skin substitutes, Helen V. Waugh and Jonathan A. Sherratt, 2007, <emphasis>Wound Repair and Regeneration</emphasis>, 15, 556-565.  <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=pubmed&amp;cmd=Retrieve&amp;dopt=AbstractPlus&amp;list_uids=17650100&amp;query_hl=1&amp;itool=pubmed_docsum">PubMed ID: 17650100</ulink>
</para>

<informalfigure float="0" id="fig_reaction_diagram">
					<mediaobject>
						<imageobject>
							<objectinfo>
								<title>model diagram</title>
							</objectinfo>
							<imagedata fileref="waugh_2007.png"/>
						</imageobject>
					</mediaobject>
					<caption>Schematic diagram of the mathematical model describing macrophage and fibroblast dynamics in diabetic wound healing.</caption>
				</informalfigure>
				
</sect1>
</article>
</documentation> 



  <units name="day">
    <unit units="second" multiplier="86400.0"/>
  </units>

  <units base_units="yes" name="cells"/>

  <units name="pg_per_mm3">
    <unit units="gram" prefix="pico"/>
    <unit units="metre" prefix="milli" exponent="-3"/>
  </units>

  <units name="microg_per_mm3">
    <unit units="gram" prefix="micro"/>
    <unit units="metre" prefix="milli" exponent="-3"/>
  </units>
  
  <units name="cells_per_mm3">
    <unit units="cells"/>
    <unit units="metre" prefix="milli" exponent="-3"/>
  </units>

  <units name="cells_per_mm3_per_day">
    <unit units="cells"/>
    <unit units="metre" prefix="milli" exponent="-3"/>
    <unit units="day" exponent="-1"/>
  </units>
  
  <units name="cells_per_pg_per_day">
    <unit units="cells"/>
    <unit units="gram" prefix="pico" exponent="-1"/>
    <unit units="day" exponent="-1"/>
  </units>

  <units name="mm3_per_cells">
    <unit units="metre" prefix="milli" exponent="3"/>
    <unit units="cells" exponent="-1"/>
  </units>

  <units name="mm3_per_cells_per_day">
    <unit units="metre" prefix="milli" exponent="3"/>
    <unit units="cells" exponent="-1"/>
    <unit units="day" exponent="-1"/>
  </units>

  <units name="mm3_per_microg">
    <unit units="metre" prefix="milli" exponent="3"/>
    <unit units="gram" prefix="micro" exponent="-1"/>
  </units>
  
  <units name="pg_per_cells_per_day">
    <unit units="gram" prefix="pico"/>
    <unit units="cells" exponent="-1"/>
    <unit units="day" exponent="-1"/>
  </units>

  <units name="microg_per_cells_per_day">
    <unit units="gram" prefix="micro"/>
    <unit units="cells" exponent="-1"/>
    <unit units="day" exponent="-1"/>
  </units>

  <units name="mm3_cells_per_pg2_day">
    <unit units="metre" prefix="milli" exponent="3"/>
    <unit units="cells"/>
    <unit units="gram" prefix="pico" exponent="-2"/>
    <unit units="day" exponent="-1"/>
  </units>

  <units name="mm6_cells_per_pg3_day">
    <unit units="metre" prefix="milli" exponent="6"/>
    <unit units="cells"/>
    <unit units="gram" prefix="pico" exponent="-3"/>
    <unit units="day" exponent="-1"/>
  </units>

  <units name="mm9_per_pg3">
    <unit units="metre" prefix="milli" exponent="9"/>
    <unit units="gram" prefix="pico" exponent="-3"/>
  </units>

  <units name="mm6_per_pg2">
    <unit units="metre" prefix="milli" exponent="6"/>
    <unit units="gram" prefix="pico" exponent="-2"/>
  </units>

  <units name="mm3_per_pg">
    <unit units="metre" prefix="milli" exponent="3"/>
    <unit units="gram" prefix="pico" exponent="-1"/>
  </units>

  <units name="mm9_per_microg3">
    <unit units="metre" prefix="milli" exponent="9"/>
    <unit units="gram" prefix="micro" exponent="-3"/>
  </units>

  <units name="mm6_per_microg2">
    <unit units="metre" prefix="milli" exponent="6"/>
    <unit units="gram" prefix="micro" exponent="-2"/>
  </units>

  <units name="per_day">
    <unit units="day" exponent="-1"/>
  </units>
  
  
  <component name="environment">
    <variable units="day" public_interface="out" name="time"/>
  </component>
  
  <component name="phi_I" cmeta:id="phi_I">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="phi_I">
        <dc:title>phi_I</dc:title>
        <dcterms:alternative>inflammatory macrophage cell density</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF>  
    <variable units="cells_per_mm3" public_interface="out" name="phi_I" initial_value="200.0"/>

    <variable units="dimensionless" public_interface="in" name="alpha"/>
    <variable units="dimensionless" public_interface="in" name="k1"/>
    <variable units="per_day" public_interface="in" name="k2"/>
    <variable units="mm3_per_cells" public_interface="in" name="k3"/>
    <variable units="mm3_per_cells" public_interface="in" name="k5"/>
    <variable units="mm3_per_microg" public_interface="in" name="k6"/>
    <variable units="per_day" public_interface="in" name="d1"/>
    <variable units="cells_per_mm3" public_interface="in" name="phi_R"/>
    <variable units="cells_per_mm3_per_day" public_interface="in" name="K_T"/>
    <variable units="cells_per_mm3" public_interface="in" name="F"/>
    <variable units="microg_per_mm3" public_interface="in" name="C"/>
    <variable units="day" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci> phi_I </ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> alpha </ci>
              <ci> K_T </ci>
            </apply>
            <apply>
              <times/>
              <ci> k1 </ci>
              <ci> k2 </ci>
              <ci> phi_I </ci>
              <apply>
                <minus/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <apply>
                  <plus/>
                  <apply>
                    <times/>
                    <ci> k3 </ci>
                    <apply>
                      <plus/>
                      <ci> phi_I </ci>
                      <ci> phi_R </ci>
                    </apply>
                  </apply>
                  <apply>
                    <times/>
                    <ci> k5 </ci>
                    <ci> F </ci>
                  </apply>
                  <apply>
                    <times/>
                    <ci> k6 </ci>
                    <ci> C </ci>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> d1 </ci>
            <ci> phi_I </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>

  <component name="phi_R" cmeta:id="phi_R">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="phi_R">
        <dc:title>phi_R</dc:title>
        <dcterms:alternative>repair macrophage cell density</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF>  
    <variable units="cells_per_mm3" public_interface="out" name="phi_R" initial_value="200.0"/>

    <variable units="dimensionless" public_interface="in" name="alpha"/>
    <variable units="dimensionless" public_interface="in" name="k1"/>
    <variable units="per_day" public_interface="in" name="k2"/>
    <variable units="mm3_per_cells" public_interface="in" name="k3"/>
    <variable units="mm3_per_cells" public_interface="in" name="k5"/>
    <variable units="mm3_per_microg" public_interface="in" name="k6"/>
    <variable units="per_day" public_interface="in" name="d1"/>
    <variable units="cells_per_mm3" public_interface="in" name="phi_I"/>
    <variable units="cells_per_mm3_per_day" public_interface="in" name="K_T"/>
    <variable units="cells_per_mm3" public_interface="in" name="F"/>
    <variable units="microg_per_mm3" public_interface="in" name="C"/>
    <variable units="day" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci> phi_R </ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <apply>
                <minus/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <ci> alpha </ci>
              </apply>
              <ci> K_T </ci>
            </apply>
            <apply>
              <times/>
              <ci> k1 </ci>
              <ci> k2 </ci>
              <ci> phi_R </ci>
              <apply>
                <minus/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <apply>
                  <plus/>
                  <apply>
                    <times/>
                    <ci> k3 </ci>
                    <apply>
                      <plus/>
                      <ci> phi_I </ci>
                      <ci> phi_R </ci>
                    </apply>
                  </apply>
                  <apply>
                    <times/>
                    <ci> k5 </ci>
                    <ci> F </ci>
                  </apply>
                  <apply>
                    <times/>
                    <ci> k6 </ci>
                    <ci> C </ci>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> d1 </ci>
            <ci> phi_R </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>

  <component name="T" cmeta:id="T">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="T">
        <dc:title>T</dc:title>
        <dcterms:alternative>transforming growth factor-beta concentration</dcterms:alternative>
        <dcterms:alternative>TGF-beta</dcterms:alternative>
      </rdf:Description>
    </rdf:RDF>  
    <variable units="pg_per_mm3" public_interface="out" name="T" initial_value="6.0"/>
    
    <variable units="pg_per_cells_per_day" public_interface="in" name="k4"/>
    <variable units="pg_per_cells_per_day" public_interface="in" name="k7"/>
    <variable units="per_day" public_interface="in" name="d2"/>
    <variable units="cells_per_mm3" public_interface="in" name="phi_I"/>
    <variable units="cells_per_mm3" public_interface="in" name="F"/>
    <variable units="day" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci> T </ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> k4 </ci>
              <ci> phi_I </ci>
            </apply>
            <apply>
              <times/>
              <ci> k7 </ci>
              <ci> F </ci>
            </apply> 
          </apply>
          <apply>
            <times/>
            <ci> d2 </ci>
            <ci> T </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>

  <component name="P" cmeta:id="P">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="P">
        <dc:title>P</dc:title>
        <dcterms:alternative>PDGF concentration</dcterms:alternative>
      </rdf:Description>
    </rdf:RDF>  
    <variable units="pg_per_mm3" public_interface="out" name="P" initial_value="2.0"/>
    
    <variable units="pg_per_cells_per_day" public_interface="in" name="k8"/>
    <variable units="pg_per_cells_per_day" public_interface="in" name="k9"/>
    <variable units="per_day" public_interface="in" name="d3"/>
    <variable units="cells_per_mm3" public_interface="in" name="phi_I"/>
    <variable units="cells_per_mm3" public_interface="in" name="phi_R"/>
    <variable units="cells_per_mm3" public_interface="in" name="F"/>
    <variable units="day" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci> P </ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> k8 </ci>
              <apply>
                <plus/>
                <ci> phi_I </ci>
                <ci> phi_R </ci>
              </apply>
            </apply>
            <apply>
              <times/>
              <ci> k9 </ci>
              <ci> F </ci>
            </apply> 
          </apply>
          <apply>
            <times/>
            <ci> d3 </ci>
            <ci> P </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>

  <component name="F" cmeta:id="F">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="F">
        <dc:title>F</dc:title>
        <dcterms:alternative>fibroblast density</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF>  
    <variable units="cells_per_mm3" public_interface="out" name="F" initial_value="50.0"/>

    <variable units="per_day" public_interface="in" name="k10"/>
    <variable units="mm3_per_cells" public_interface="in" name="k3"/>
    <variable units="mm3_per_cells" public_interface="in" name="k5"/>
    <variable units="mm3_per_microg" public_interface="in" name="k6"/>
    <variable units="per_day" public_interface="in" name="d4"/>
    <variable units="cells_per_mm3" public_interface="in" name="phi_I"/>
    <variable units="cells_per_mm3" public_interface="in" name="phi_R"/>
    <variable units="cells_per_mm3_per_day" public_interface="in" name="M_P"/>
    <variable units="microg_per_mm3" public_interface="in" name="C"/>
    <variable units="day" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci> F </ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/>
            <ci> M_P </ci>
            <apply>
              <times/>
              <ci> k10 </ci>
              <ci> F </ci>
              <apply>
                <minus/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <apply>
                  <plus/>
                  <apply>
                    <times/>
                    <ci> k3 </ci>
                    <apply>
                      <plus/>
                      <ci> phi_I </ci>
                      <ci> phi_R </ci>
                    </apply>
                  </apply>
                  <apply>
                    <times/>
                    <ci> k5 </ci>
                    <ci> F </ci>
                  </apply>
                  <apply>
                    <times/>
                    <ci> k6 </ci>
                    <ci> C </ci>
                  </apply>
                </apply>
              </apply>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> d4 </ci>
            <ci> F </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>

  <component name="C" cmeta:id="C">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="C">
        <dc:title>C</dc:title>
        <dcterms:alternative>collagen concentration</dcterms:alternative>
      </rdf:Description>
    </rdf:RDF>  
    <variable units="microg_per_mm3" public_interface="out" name="C" initial_value="10.0"/>
    
    <variable units="microg_per_cells_per_day" public_interface="in" name="k11"/>
    <variable units="dimensionless" public_interface="in" name="f_T"/>
    <variable units="dimensionless" public_interface="in" name="g_C"/>
    <variable units="mm3_per_cells_per_day" public_interface="in" name="d5"/>
    <variable units="cells_per_mm3" public_interface="in" name="F"/>
    <variable units="day" public_interface="in" name="time"/>
    
    <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> k11 </ci>
            <ci> F </ci>
            <ci> f_T </ci>
            <ci> g_C </ci>
          </apply>
          <apply>
            <times/>
            <ci> d5 </ci>
            <ci> F </ci>
            <ci> C </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>

  <component name="H" cmeta:id="H">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="H">
        <dc:title>H</dc:title>
        <dcterms:alternative>hyaluronan concentration</dcterms:alternative>
      </rdf:Description>
    </rdf:RDF>  
    <variable units="microg_per_mm3" public_interface="out" name="H" initial_value="0.01"/>
    
    <variable units="microg_per_cells_per_day" public_interface="in" name="k12"/>
    <variable units="per_day" public_interface="in" name="d6"/>
    <variable units="cells_per_mm3" public_interface="in" name="F"/>
    <variable units="day" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci> H </ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> k12 </ci>
            <ci> F </ci>
          </apply>
          <apply>
            <times/>
            <ci> d6 </ci>
            <ci> H </ci>
          </apply>
        </apply>
      </apply>
    </math>
  </component>

  <component name="K_T" cmeta:id="K_T">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="K_T">
        <dc:title>K_T</dc:title>
        <dcterms:alternative>monocyte migration to the wound in response to TGF-beta</dcterms:alternative>
      </rdf:Description>
    </rdf:RDF>  
    <variable units="cells_per_mm3_per_day" public_interface="out" name="K_T"/>
    
    <variable units="mm6_cells_per_pg3_day" name="tau1" initial_value="-2.47"/>
    <variable units="mm3_cells_per_pg2_day" name="tau2" initial_value="21.94"/>
    <variable units="cells_per_pg_per_day" name="tau3" initial_value="6.41"/>
    <variable units="cells_per_mm3_per_day" name="tau4" initial_value="1.75"/>
    
    <variable units="pg_per_mm3" public_interface="in" name="T"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci> K_T </ci> 
        <apply>
          <plus/>
          <apply>
            <times/>
            <ci> tau1 </ci>
            <apply>
              <power/>
              <ci> T </ci>
              <cn cellml:units="dimensionless"> 3.0 </cn>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> tau2 </ci>
            <apply>
              <power/>
              <ci> T </ci>
              <cn cellml:units="dimensionless"> 2.0 </cn>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> tau3 </ci>
            <ci> T </ci>
          </apply>
          <ci> tau4 </ci>
        </apply>
      </apply>
    </math>
  </component>

  <component name="M_P" cmeta:id="M_P">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="M_P">
        <dc:title>M_P</dc:title>
        <dcterms:alternative>fibroblast migration to the wound in response to PDGF</dcterms:alternative>
      </rdf:Description>
    </rdf:RDF>  
    <variable units="cells_per_mm3_per_day" public_interface="out" name="M_P"/>
    
    <variable units="mm6_cells_per_pg3_day" name="tau1" initial_value="15.333"/>
    <variable units="mm3_cells_per_pg2_day" name="tau2" initial_value="-167.21"/>
    <variable units="cells_per_pg_per_day" name="tau3" initial_value="452.38"/>
    <variable units="cells_per_mm3_per_day" name="tau4" initial_value="2.6593"/>
    
    <variable units="pg_per_mm3" public_interface="in" name="P"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci> M_P </ci> 
        <apply>
          <plus/>
          <apply>
            <times/>
            <ci> tau1 </ci>
            <apply>
              <power/>
              <ci> P </ci>
              <cn cellml:units="dimensionless"> 3.0 </cn>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> tau2 </ci>
            <apply>
              <power/>
              <ci> P </ci>
              <cn cellml:units="dimensionless"> 2.0 </cn>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> tau3 </ci>
            <ci> P </ci>
          </apply>
          <ci> tau4 </ci>
        </apply>
      </apply>
    </math>
  </component>

  <component name="f_T" cmeta:id="f_T">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="f_T">
        <dc:title>f_T</dc:title>
        <dcterms:alternative>fibroblast synthesis of collagen in response to TGF-beta</dcterms:alternative>
      </rdf:Description>
    </rdf:RDF>  
    <variable units="dimensionless" public_interface="out" name="f_T"/>
    
    <variable units="mm9_per_pg3" name="tau1" initial_value="0.0092"/>
    <variable units="mm6_per_pg2" name="tau2" initial_value="-0.1552"/>
    <variable units="mm3_per_pg" name="tau3" initial_value="0.6279"/>
    <variable units="dimensionless" name="tau4" initial_value="0.2527"/>
    
    <variable units="pg_per_mm3" public_interface="in" name="T"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci> f_T </ci> 
        <apply>
          <plus/>
          <apply>
            <times/>
            <ci> tau1 </ci>
            <apply>
              <power/>
              <ci> T </ci>
              <cn cellml:units="dimensionless"> 3.0 </cn>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> tau2 </ci>
            <apply>
              <power/>
              <ci> T </ci>
              <cn cellml:units="dimensionless"> 2.0 </cn>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> tau3 </ci>
            <ci> T </ci>
          </apply>
          <ci> tau4 </ci>
        </apply>
      </apply>
    </math>
  </component>

  <component name="g_C" cmeta:id="g_C">
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
      <rdf:Description rdf:about="g_C">
        <dc:title>g_C</dc:title>
        <dcterms:alternative>effect of collagen on its own synthesis by fibroblasts</dcterms:alternative>
      </rdf:Description>
    </rdf:RDF>  
    <variable units="dimensionless" public_interface="out" name="g_C"/>
    
    <variable units="mm9_per_microg3" name="tau1" initial_value="-4.33E-10"/>
    <variable units="mm6_per_microg2" name="tau2" initial_value="0.0000009"/>
    <variable units="mm3_per_microg" name="tau3" initial_value="-0.00055"/>
    <variable units="dimensionless" name="tau4" initial_value="0.13"/>
    
    <variable units="microg_per_mm3" public_interface="in" name="C"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci> g_C </ci> 
        <apply>
          <plus/>
          <apply>
            <times/>
            <ci> tau1 </ci>
            <apply>
              <power/>
              <ci> C </ci>
              <cn cellml:units="dimensionless"> 3.0 </cn>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> tau2 </ci>
            <apply>
              <power/>
              <ci> C </ci>
              <cn cellml:units="dimensionless"> 2.0 </cn>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> tau3 </ci>
            <ci> C </ci>
          </apply>
          <ci> tau4 </ci>
        </apply>
      </apply>
    </math>
  </component>

  <component name="model_parameters">
    <variable units="dimensionless" public_interface="out" name="alpha"/> 
    <variable units="dimensionless" public_interface="out" name="k1" initial_value="0.05"/>
    <variable units="per_day" public_interface="out" name="k2" initial_value="0.693"/>
    <variable units="mm3_per_cells" public_interface="out" name="k3" initial_value="0.002"/>
    <variable units="pg_per_cells_per_day" public_interface="out" name="k4" initial_value="0.07"/>
    <variable units="mm3_per_cells" public_interface="out" name="k5" initial_value="0.0025"/>
    <variable units="mm3_per_microg" public_interface="out" name="k6" initial_value="0.0004"/>
    <variable units="pg_per_cells_per_day" public_interface="out" name="k7" initial_value="0.004"/>
    <variable units="pg_per_cells_per_day" public_interface="out" name="k8" initial_value="0.015"/>
    <variable units="pg_per_cells_per_day" public_interface="out" name="k9" initial_value="0.0015"/>
    <variable units="per_day" public_interface="out" name="k10" initial_value="0.924"/>
    <variable units="microg_per_cells_per_day" public_interface="out" name="k11" initial_value="20.0"/>
    <variable units="microg_per_cells_per_day" public_interface="out" name="k12" initial_value="0.01"/>
    <variable units="per_day" public_interface="out" name="d1" initial_value="0.2"/>
    <variable units="per_day" public_interface="out" name="d2" initial_value="9.1"/>
    <variable units="per_day" public_interface="out" name="d3" initial_value="4.0"/>
    <variable units="per_day" public_interface="out" name="d4" initial_value="1.0"/>
    <variable units="mm3_per_cells_per_day" public_interface="out" name="d5" initial_value="1.5E-5"/>
    <variable units="per_day" public_interface="out" name="d6" initial_value="0.7"/>
  
    <variable units="microg_per_mm3" public_interface="in" name="H"/>

    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci> alpha </ci> 
        <apply>
          <plus/>
          <apply>
            <minus/>
            <apply>
              <times/>
              <cn cellml:units="dimensionless"> 0.197 </cn>
              <apply>
                <log/>
                <ci> H </ci>
              </apply>
            </apply>
          </apply>
          <cn cellml:units="dimensionless"> 0.4407 </cn>
        </apply>
      </apply>
    </math>
  </component>
  
  <connection>
    <map_components component_2="environment" component_1="phi_I"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>

  <connection>
    <map_components component_2="environment" component_1="phi_R"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>

  <connection>
    <map_components component_2="environment" component_1="T"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>

  <connection>
    <map_components component_2="environment" component_1="P"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>

  <connection>
    <map_components component_2="environment" component_1="F"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>

  <connection>
    <map_components component_2="environment" component_1="C"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>

  <connection>
    <map_components component_2="environment" component_1="H"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>

  <connection>
    <map_components component_2="phi_R" component_1="phi_I"/>
    <map_variables variable_2="phi_I" variable_1="phi_I"/>
    <map_variables variable_2="phi_R" variable_1="phi_R"/>
  </connection> 

  <connection>
    <map_components component_2="model_parameters" component_1="phi_I"/>
    <map_variables variable_2="alpha" variable_1="alpha"/>
    <map_variables variable_2="k1" variable_1="k1"/>
    <map_variables variable_2="k2" variable_1="k2"/>
    <map_variables variable_2="k3" variable_1="k3"/>
    <map_variables variable_2="k5" variable_1="k5"/>
    <map_variables variable_2="k6" variable_1="k6"/>
    <map_variables variable_2="d1" variable_1="d1"/>
  </connection>

  <connection>
    <map_components component_2="K_T" component_1="phi_I"/>
    <map_variables variable_2="K_T" variable_1="K_T"/>
  </connection> 

  <connection>
    <map_components component_2="F" component_1="phi_I"/>
    <map_variables variable_2="phi_I" variable_1="phi_I"/>
    <map_variables variable_2="F" variable_1="F"/>
  </connection>

  <connection>
    <map_components component_2="C" component_1="phi_I"/>
    <map_variables variable_2="C" variable_1="C"/>
  </connection>

  <connection>
    <map_components component_2="model_parameters" component_1="phi_R"/>
    <map_variables variable_2="alpha" variable_1="alpha"/>
    <map_variables variable_2="k1" variable_1="k1"/>
    <map_variables variable_2="k2" variable_1="k2"/>
    <map_variables variable_2="k3" variable_1="k3"/>
    <map_variables variable_2="k5" variable_1="k5"/>
    <map_variables variable_2="k6" variable_1="k6"/>
    <map_variables variable_2="d1" variable_1="d1"/>
  </connection>

  <connection>
    <map_components component_2="K_T" component_1="phi_R"/>
    <map_variables variable_2="K_T" variable_1="K_T"/>
  </connection>

  <connection>
    <map_components component_2="F" component_1="phi_R"/>
    <map_variables variable_2="phi_R" variable_1="phi_R"/>
    <map_variables variable_2="F" variable_1="F"/>
  </connection>

  <connection>
    <map_components component_2="C" component_1="phi_R"/>
    <map_variables variable_2="C" variable_1="C"/>
  </connection> 

  <connection>
    <map_components component_2="T" component_1="phi_I"/>
    <map_variables variable_2="phi_I" variable_1="phi_I"/>
  </connection>

  <connection>
    <map_components component_2="F" component_1="T"/>
    <map_variables variable_2="F" variable_1="F"/>
  </connection>

  <connection>
    <map_components component_2="model_parameters" component_1="T"/>
    <map_variables variable_2="k4" variable_1="k4"/>
    <map_variables variable_2="k7" variable_1="k7"/>
    <map_variables variable_2="d2" variable_1="d2"/>
  </connection>

  <connection>
    <map_components component_2="P" component_1="phi_I"/>
    <map_variables variable_2="phi_I" variable_1="phi_I"/>
  </connection>

  <connection>
    <map_components component_2="P" component_1="phi_R"/>
    <map_variables variable_2="phi_R" variable_1="phi_R"/>
  </connection>

  <connection>
    <map_components component_2="F" component_1="P"/>
    <map_variables variable_2="F" variable_1="F"/>
  </connection>

  <connection>
    <map_components component_2="model_parameters" component_1="P"/>
    <map_variables variable_2="k8" variable_1="k8"/>
    <map_variables variable_2="k9" variable_1="k9"/>
    <map_variables variable_2="d3" variable_1="d3"/>
  </connection>

  <connection>
    <map_components component_2="C" component_1="F"/>
    <map_variables variable_2="C" variable_1="C"/>
    <map_variables variable_2="F" variable_1="F"/>
  </connection>

  <connection>
    <map_components component_2="M_P" component_1="F"/>
    <map_variables variable_2="M_P" variable_1="M_P"/>
  </connection>  

  <connection>
    <map_components component_2="model_parameters" component_1="F"/>
    <map_variables variable_2="k10" variable_1="k10"/>
    <map_variables variable_2="k3" variable_1="k3"/>
    <map_variables variable_2="k5" variable_1="k5"/>
    <map_variables variable_2="k6" variable_1="k6"/>
    <map_variables variable_2="d4" variable_1="d4"/>
  </connection>

  <connection>
    <map_components component_2="f_T" component_1="C"/>
    <map_variables variable_2="f_T" variable_1="f_T"/>
  </connection>

  <connection>
    <map_components component_2="g_C" component_1="C"/>
    <map_variables variable_2="g_C" variable_1="g_C"/>
    <map_variables variable_2="C" variable_1="C"/>
  </connection>

  <connection>
    <map_components component_2="model_parameters" component_1="C"/>
    <map_variables variable_2="k11" variable_1="k11"/>
    <map_variables variable_2="d5" variable_1="d5"/>
  </connection>

  <connection>
    <map_components component_2="F" component_1="H"/>
    <map_variables variable_2="F" variable_1="F"/>
  </connection>

  <connection>
    <map_components component_2="model_parameters" component_1="H"/>
    <map_variables variable_2="k12" variable_1="k12"/>
    <map_variables variable_2="d6" variable_1="d6"/>
    <map_variables variable_2="H" variable_1="H"/>
  </connection>  

  <connection>
    <map_components component_2="T" component_1="K_T"/>
    <map_variables variable_2="T" variable_1="T"/>
  </connection>

  <connection>
    <map_components component_2="P" component_1="M_P"/>
    <map_variables variable_2="P" variable_1="P"/>
  </connection> 

  <connection>
    <map_components component_2="T" component_1="f_T"/>
    <map_variables variable_2="T" variable_1="T"/>
  </connection>

<rdf:RDF>
  <rdf:Bag rdf:about="rdf:#6e59ab32-c7a1-407e-831a-ad4efe0d40e5">
    <rdf:li>diabetes</rdf:li>
    <rdf:li>wound healing</rdf:li>
    <rdf:li>macrophage</rdf:li>
  </rdf:Bag>
  <rdf:Seq rdf:about="rdf:#39c2a0b4-8997-47cb-bfd9-5e60825ee9f4">
    <rdf:li rdf:resource="rdf:#5a3d469e-f0f6-4f17-9b5f-5f09cb03bc42"/>
    <rdf:li rdf:resource="rdf:#74380275-f841-43d5-bfc6-0c0263849170"/>
  </rdf:Seq>
  <rdf:Description rdf:about="rdf:#7a58f69d-bd49-4fff-ab20-d4084a89602b">
    <vCard:FN/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ae1ea017-e59b-4a34-99ad-5f1432e3c1bc">
    <dcterms:W3CDTF>2007-07</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ddcfc987-2cc5-4220-8c90-42e90a2dc410">
    <vCard:Given>Helen</vCard:Given>
    <vCard:Family>Waugh</vCard:Family>
    <vCard:Other>V</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#74380275-f841-43d5-bfc6-0c0263849170">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#1fd71b88-bc8e-4051-8383-c4047beffa1c"/>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dc:publisher>The University of Auckland, Bioengineering Institute</dc:publisher>
    <cmeta:comment rdf:resource="rdf:#3c8ac288-13ea-4158-a883-53e6fd7655b3"/>
    <dcterms:created rdf:resource="rdf:#fa1a96be-dba3-423a-a00b-963d13946e49"/>
    <dc:creator rdf:resource="rdf:#8c49ac41-13cd-4e2b-801d-341d7e6721a7"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#e904768f-4991-45bc-91a1-19aa328e8507">
    <dc:creator rdf:resource="rdf:#a0b2e9e1-1c61-4c33-b81a-3f230df507b8"/>
    <rdf:value>This is a CellML description of Waugh and Sherratt's 2006 mathematical model of macrophage dynamics in diabetic wound healing.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5c551607-4a40-4a3c-b4a3-93a1607dec30">
    <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:#8c49ac41-13cd-4e2b-801d-341d7e6721a7">
    <vCard:ORG rdf:resource="rdf:#fa52f2a0-576a-459f-a7f7-a05a3e942cda"/>
    <vCard:EMAIL rdf:resource="rdf:#5c551607-4a40-4a3c-b4a3-93a1607dec30"/>
    <vCard:N rdf:resource="rdf:#c8c515bb-e7e6-4897-9ee9-6274cb179b8c"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#fa1a96be-dba3-423a-a00b-963d13946e49">
    <dcterms:W3CDTF>2007-11-14T00:00:00+00:00</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3852a742-6ee2-496c-a2ea-4bc4c08601d6">
    <bqs:subject_type>keyword</bqs:subject_type>
    <rdf:value rdf:resource="rdf:#6e59ab32-c7a1-407e-831a-ad4efe0d40e5"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#9fbd5e10-333b-4779-bd0f-d53f419b6416">
    <dc:title>Wound Repair and Regeneration</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#552579d8-f713-4e12-a895-cf5886652033">
    <bqs:Pubmed_id>17650100</bqs:Pubmed_id>
    <bqs:JournalArticle rdf:resource="rdf:#d145f188-f033-4771-9b11-5c23a1219516"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#3c8ac288-13ea-4158-a883-53e6fd7655b3">
    <dc:creator rdf:resource="rdf:#7a58f69d-bd49-4fff-ab20-d4084a89602b"/>
    <rdf:value/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#5a3d469e-f0f6-4f17-9b5f-5f09cb03bc42">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#ddcfc987-2cc5-4220-8c90-42e90a2dc410"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#fa52f2a0-576a-459f-a7f7-a05a3e942cda">
    <vCard:Orgname>The University of Auckland</vCard:Orgname>
    <vCard:Orgunit>The Bioengineering Institute</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#c8c515bb-e7e6-4897-9ee9-6274cb179b8c">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#1fd71b88-bc8e-4051-8383-c4047beffa1c">
    <vCard:Given>Jonathan</vCard:Given>
    <vCard:Family>Sherratt</vCard:Family>
    <vCard:Other>A</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#dbdfad3d-2cd6-4d78-835b-05dfbd9dc72b">
    <dc:subject rdf:resource="rdf:#3852a742-6ee2-496c-a2ea-4bc4c08601d6"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#d145f188-f033-4771-9b11-5c23a1219516">
    <dc:creator rdf:resource="rdf:#39c2a0b4-8997-47cb-bfd9-5e60825ee9f4"/>
    <dc:title>Modeling the effects of treating diabetic wounds with engineered skin substitutes</dc:title>
    <bqs:volume>15</bqs:volume>
    <bqs:first_page>556</bqs:first_page>
    <bqs:Journal rdf:resource="rdf:#9fbd5e10-333b-4779-bd0f-d53f419b6416"/>
    <dcterms:issued rdf:resource="rdf:#ae1ea017-e59b-4a34-99ad-5f1432e3c1bc"/>
    <bqs:last_page>565</bqs:last_page>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#a0b2e9e1-1c61-4c33-b81a-3f230df507b8">
    <vCard:FN>Catherine Lloyd</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="#waugh_model_2007">
    <dc:title>
        Waugh and Sherratt's 2006 mathematical model of macrophage dynamics in diabetic wound healing.
      </dc:title>
    <cmeta:comment rdf:resource="rdf:#e904768f-4991-45bc-91a1-19aa328e8507"/>
    <bqs:reference rdf:resource="rdf:#552579d8-f713-4e12-a895-cf5886652033"/>
    <bqs:reference rdf:resource="rdf:#dbdfad3d-2cd6-4d78-835b-05dfbd9dc72b"/>
  </rdf:Description>
</rdf:RDF>
</model>