- Author:
- Tommy Yu <tommy.yu@auckland.ac.nz>
- Date:
- 2015-02-12 17:33:39+13:00
- Desc:
- Pruning off that extra head
- Permanent Source URI:
- https://models.physiomeproject.org/workspace/halloy_bernard_loussouarn_goldbeter_2002/rawfile/9b24d1c9771a149c9dc7af9134d4f35b0ad36629/halloy_2002_1.1.cellml
<?xml version="1.0"?>
<model xmlns="http://www.cellml.org/cellml/1.1#" xmlns:cmeta="http://www.cellml.org/metadata/1.1#" cmeta:id="halloy_2002_1.1" name="halloy_2002_1.1">
<!-- documentation -->
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
<articleinfo>
<title>The Follicular Automaton Model for Hair Cycles</title>
<author>
<firstname>Catherine</firstname>
<surname>Lloyd</surname>
<affiliation>
<shortaffil>Auckland Bioengineering Institute, The University of Auckland</shortaffil>
</affiliation>
</author>
<author>
<firstname>Mona</firstname>
<surname>Zhu</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 to recreate the same steady state output of the original published model (as seen in in figure 5(d)). The units have been checked and they are consistent. This model has been re-coded to remove the reaction element and replace it with a simple MathML description of the model reaction kinetics. This is thought to be truer to the original publication, and information regarding the enzyme kinetics etc will later be added to the metadata through use of an ontology. The CellMl model represents the deterministic version of the follicular automaton model. Currently CellML cannot be used to describe stochastic variables.
</para>
</section>
<sect1 id="sec_structure">
<title>Model Structure</title>
<para>
ABSTRACT: Human scalp hair consists of a set of about 10 to the 5 follicles which progress independently through developmental cycles. Each hair follicle successively goes through the anagen (A), catagen (C), telogen (T) and latency (L) phases that correspond, respectively, to growth, arrest and hair shedding before a new anagen phase is initiated. Long-term experimental observations in a group of ten male, alopecic and non-alopecic volunteers allowed determination of the characteristics of hair follicle cycles. On the basis of these observations, we previously proposed a follicular automaton model to simulate the dynamics of human hair cycles and the development of different patterns of alopecia [Halloy et al. (2000) Proc. Natl Acad. Sci. U.S.A.97, 8328-8333]. The automaton model is defined by a set of rules that govern the stochastic transitions of each follicle between the successive states A, T, L and the subsequent return to A. These transitions occur independently for each follicle, after time intervals given stochastically by a distribution characterized by a mean and a standard deviation. The follicular automaton model was shown to account both for the dynamical transitions observed in a single follicle, and for the behaviour of an ensemble of independently cycling follicles. Here, we extend these results and investigate additional properties of the model. We present a deterministic version of the follicular automaton. We show that numerical simulations of the stochastic version of the automaton yield steady-state level of follicles in the different phases which approach the levels predicted by the deterministic equations as the number of follicles progressively increases. Only the stochastic version can successfully reproduce the fluctuations of the fractions of follicles in each of the three phases, observed in small follicle populations. When the standard deviation is reduced or when the follicles become otherwise synchronized, e.g. by a periodic external signal inducing the transition of anagen follicles into telogen phase, large-amplitude oscillations occur in the fractions of follicles in the three phases. These oscillations are not observed in humans but are reminiscent of the phenomenon of moulting observed in a number of mammalian species. Copyright 2002 Elsevier Science Ltd.
</para>
<para>
The original paper reference is cited below:
</para>
<para>The Follicular Automaton Model: Effect of Stochasticity and of Synchronization of Hair Cycles, J. Halloy, B.A. Bernard, G. Loussouarn and A. Goldbeter, 2002, <emphasis>Journal of Theoretical Biology</emphasis>
, 214, 469-479. <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&db=PubMed&list_uids=11846603&dopt=Abstract">PubMed ID: 11846603</ulink>
</para>
<informalfigure float="0" id="fig_reaction_diagram">
<mediaobject>
<imageobject>
<objectinfo>
<title>reaction_diagram</title>
</objectinfo>
<imagedata fileref="halloy_2002.png"/>
</imageobject>
</mediaobject>
<caption>The above diagram represents the transition of a model hair follicle from anagen (A) to telogen (T) to latency (L) phase, successively. After phase T, the follicle may either die or miniaturise (transition to M; this usually occurs after a critical number of cycles), or complete a cycle by entering a new A phase.</caption>
</informalfigure>
</sect1>
</article>
</documentation>
<!-- units -->
<units name="month">
<unit multiplier="2592000" units="second"/>
</units>
<units name="flux">
<unit units="fraction"/>
<unit exponent="-1" units="month"/>
</units>
<units name="first_order_rate_constant">
<unit exponent="-1" units="month"/>
</units>
<!-- components -->
<component name="environment">
<variable cmeta:id="environment_time" name="time" public_interface="out" units="month"/>
</component>
<component cmeta:id="L" name="L">
<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
<rdf:Description rdf:about="#L">
<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">L</dc:title>
<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">latency phase</dcterms:alternative>
</rdf:Description>
</rdf:RDF>
<variable cmeta:id="L_L" name="L" private_interface="in" public_interface="out" units="fraction"/>
<variable initial_value="0.33" name="L_init" private_interface="out" units="fraction"/>
<variable name="mu_T" private_interface="out" public_interface="in" units="month"/>
<variable name="mu_L" private_interface="out" public_interface="in" units="month"/>
<variable name="mu_A" private_interface="out" public_interface="in" units="month"/>
<variable name="epsilon" private_interface="out" public_interface="in" units="first_order_rate_constant"/>
<variable name="M" private_interface="out" public_interface="in" units="fraction"/>
<variable name="T" private_interface="out" public_interface="in" units="fraction"/>
<variable name="time" private_interface="out" public_interface="in" units="month"/>
<variable name="L_0" private_interface="in" units="dimensionless"/>
</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/">anagen phase</dcterms:alternative>
</rdf:Description>
</rdf:RDF>
<variable cmeta:id="A_A" name="A" private_interface="in" public_interface="out" units="fraction"/>
<variable initial_value="0.34" name="A_init" private_interface="out" units="fraction"/>
<variable name="mu_A" private_interface="out" public_interface="in" units="month"/>
<variable name="mu_L" private_interface="out" public_interface="in" units="month"/>
<variable name="mu_T" private_interface="out" public_interface="in" units="month"/>
<variable name="L" private_interface="out" public_interface="in" units="fraction"/>
<variable name="time" private_interface="out" public_interface="in" units="month"/>
<variable name="A_0" private_interface="in" units="dimensionless"/>
<variable initial_value="0" name="eps_zero" private_interface="out" units="first_order_rate_constant"/>
<variable initial_value="0" name="M_zero" private_interface="out" units="fraction"/>
</component>
<component cmeta:id="T" name="T">
<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
<rdf:Description rdf:about="#T">
<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">T</dc:title>
<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">telogen phase</dcterms:alternative>
</rdf:Description>
</rdf:RDF>
<variable cmeta:id="T_T" name="T" private_interface="in" public_interface="out" units="fraction"/>
<variable initial_value="0.33" name="T_init" private_interface="out" units="fraction"/>
<variable name="mu_A" private_interface="out" public_interface="in" units="month"/>
<variable name="mu_T" private_interface="out" public_interface="in" units="month"/>
<variable name="mu_L" private_interface="out" public_interface="in" units="month"/>
<variable name="A" private_interface="out" public_interface="in" units="fraction"/>
<variable name="time" private_interface="out" public_interface="in" units="month"/>
<variable name="T_0" private_interface="in" units="dimensionless"/>
<variable initial_value="0" name="eps_zero" private_interface="out" units="first_order_rate_constant"/>
<variable initial_value="0" name="M_zero" private_interface="out" units="fraction"/>
</component>
<component cmeta:id="M" name="M">
<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
<rdf:Description rdf:about="#M">
<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">M</dc:title>
<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">miniaturisation phase</dcterms:alternative>
</rdf:Description>
</rdf:RDF>
<variable cmeta:id="M_M" initial_value="0" name="M" public_interface="out" units="fraction"/>
<variable name="epsilon" public_interface="in" units="first_order_rate_constant"/>
<variable name="T" public_interface="in" units="fraction"/>
<variable name="time" public_interface="in" units="month"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>M</ci>
</apply>
<apply>
<times/>
<ci>epsilon</ci>
<ci>T</ci>
</apply>
</apply>
</math>
</component>
<component name="reaction_constants">
<variable initial_value="2" name="mu_L" public_interface="out" units="month"/>
<variable initial_value="3.6" name="mu_A" public_interface="out" units="month"/>
<variable initial_value="2.4" name="mu_T" public_interface="out" units="month"/>
<variable initial_value="0" name="epsilon" public_interface="out" units="first_order_rate_constant"/>
</component>
<!-- groups -->
<group>
<relationship_ref relationship="encapsulation"/>
<component_ref component="A">
<component_ref component="Ai"/>
<component_ref component="A_0i"/>
</component_ref>
</group>
<group>
<relationship_ref relationship="encapsulation"/>
<component_ref component="T">
<component_ref component="Ti"/>
<component_ref component="T_0i"/>
</component_ref>
</group>
<group>
<relationship_ref relationship="encapsulation"/>
<component_ref component="L">
<component_ref component="Li"/>
<component_ref component="L_0i"/>
</component_ref>
</group>
<!-- imports -->
<import xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="Y.cellml">
<component component_ref="interface" name="Ai"/>
</import>
<import xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="Y.cellml">
<component component_ref="interface" name="Ti"/>
</import>
<import xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="Y.cellml">
<component component_ref="interface" name="Li"/>
</import>
<import xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="Y_0.cellml">
<component component_ref="interface" name="A_0i"/>
</import>
<import xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="Y_0.cellml">
<component component_ref="interface" name="T_0i"/>
</import>
<import xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="Y_0.cellml">
<component component_ref="interface" name="L_0i"/>
</import>
<import xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="fraction.cellml">
<units name="fraction" units_ref="fraction"/>
</import>
<!-- connections -->
<connection>
<map_components component_1="L" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="A" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="T" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="M" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="L" component_2="T"/>
<map_variables variable_1="T" variable_2="T"/>
</connection>
<connection>
<map_components component_1="L" component_2="A"/>
<map_variables variable_1="L" variable_2="L"/>
</connection>
<connection>
<map_components component_1="A" component_2="T"/>
<map_variables variable_1="A" variable_2="A"/>
</connection>
<connection>
<map_components component_1="L" component_2="M"/>
<map_variables variable_1="M" variable_2="M"/>
</connection>
<connection>
<map_components component_1="T" component_2="M"/>
<map_variables variable_1="T" variable_2="T"/>
</connection>
<connection>
<map_components component_1="A" component_2="reaction_constants"/>
<map_variables variable_1="mu_L" variable_2="mu_L"/>
<map_variables variable_1="mu_A" variable_2="mu_A"/>
<map_variables variable_1="mu_T" variable_2="mu_T"/>
</connection>
<connection>
<map_components component_1="T" component_2="reaction_constants"/>
<map_variables variable_1="mu_T" variable_2="mu_T"/>
<map_variables variable_1="mu_A" variable_2="mu_A"/>
<map_variables variable_1="mu_L" variable_2="mu_L"/>
</connection>
<connection>
<map_components component_1="L" component_2="reaction_constants"/>
<map_variables variable_1="mu_L" variable_2="mu_L"/>
<map_variables variable_1="mu_T" variable_2="mu_T"/>
<map_variables variable_1="mu_A" variable_2="mu_A"/>
<map_variables variable_1="epsilon" variable_2="epsilon"/>
</connection>
<connection>
<map_components component_1="M" component_2="reaction_constants"/>
<map_variables variable_1="epsilon" variable_2="epsilon"/>
</connection>
<connection>
<map_components component_1="Ai" component_2="A"/>
<map_variables variable_1="Y" variable_2="A"/>
<map_variables variable_1="Y_init" variable_2="A_init"/>
<map_variables variable_1="mu_Z" variable_2="mu_L"/>
<map_variables variable_1="mu_Y" variable_2="mu_A"/>
<map_variables variable_1="epsilon" variable_2="eps_zero"/>
<map_variables variable_1="M" variable_2="M_zero"/>
<map_variables variable_1="Z" variable_2="L"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="Ti" component_2="T"/>
<map_variables variable_1="Y" variable_2="T"/>
<map_variables variable_1="Y_init" variable_2="T_init"/>
<map_variables variable_1="mu_Z" variable_2="mu_A"/>
<map_variables variable_1="mu_Y" variable_2="mu_T"/>
<map_variables variable_1="epsilon" variable_2="eps_zero"/>
<map_variables variable_1="M" variable_2="M_zero"/>
<map_variables variable_1="Z" variable_2="A"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="Li" component_2="L"/>
<map_variables variable_1="Y" variable_2="L"/>
<map_variables variable_1="Y_init" variable_2="L_init"/>
<map_variables variable_1="mu_Z" variable_2="mu_T"/>
<map_variables variable_1="mu_Y" variable_2="mu_L"/>
<map_variables variable_1="epsilon" variable_2="epsilon"/>
<map_variables variable_1="M" variable_2="M"/>
<map_variables variable_1="Z" variable_2="T"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="A_0i" component_2="A"/>
<map_variables variable_1="Y_0" variable_2="A_0"/>
<map_variables variable_1="mu_T" variable_2="mu_T"/>
<map_variables variable_1="mu_L" variable_2="mu_L"/>
<map_variables variable_1="mu_A" variable_2="mu_A"/>
</connection>
<connection>
<map_components component_1="T_0i" component_2="T"/>
<map_variables variable_1="Y_0" variable_2="T_0"/>
<map_variables variable_1="mu_T" variable_2="mu_T"/>
<map_variables variable_1="mu_L" variable_2="mu_L"/>
<map_variables variable_1="mu_A" variable_2="mu_A"/>
</connection>
<connection>
<map_components component_1="L_0i" component_2="L"/>
<map_variables variable_1="Y_0" variable_2="L_0"/>
<map_variables variable_1="mu_T" variable_2="mu_T"/>
<map_variables variable_1="mu_L" variable_2="mu_L"/>
<map_variables variable_1="mu_A" variable_2="mu_A"/>
</connection>
<!-- metadata -->
<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
<rdf:Description rdf:about="halloy_2002_1.1.cellml">
<dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:parseType="Resource">
<vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
<vCard:Family>Zhu</vCard:Family>
<vCard:Given>Xin Yue</vCard:Given>
<vCard:Other>Mona</vCard:Other>
</vCard:N>
<vCard:EMAIL xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
<rdf:value>xzhu055@aucklanduni.ac.nz</rdf:value>
<rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
</vCard:EMAIL>
<vCard:ORG xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
<vCard:Orgname>The University of Auckland</vCard:Orgname>
<vCard:Orgunit>Auckland Bioengineering Institute</vCard:Orgunit>
</vCard:ORG>
</dc:creator>
<dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:parseType="Resource">
<vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
<vCard:Family>Lloyd</vCard:Family>
<vCard:Given>Catherine</vCard:Given>
<vCard:Other>May</vCard:Other>
</vCard:N>
<vCard:EMAIL xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" 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 xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
<vCard:Orgname>The University of Auckland</vCard:Orgname>
<vCard:Orgunit>Auckland Bioengineering Institute</vCard:Orgunit>
</vCard:ORG>
</dc:creator>
<dcterms:created xmlns:dcterms="http://purl.org/dc/terms/" rdf:parseType="Resource">
<dcterms:W3CDTF>19-01-2011</dcterms:W3CDTF>
</dcterms:created>
</rdf:Description>
<rdf:Description rdf:about="#halloy_2002_1.1.cellml">
<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">The Follicular Automaton Model: Effect of Stochasticity and of
Synchronization of Hair Cycles</dc:title>
<cmeta:comment rdf:parseType="Resource">
<rdf:value>This is the CellML description of Halloy et al's 2002 follicular automaton model.
</rdf:value>
<dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:parseType="Resource">
<vCard:FN xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Mona Zhu</vCard:FN>
</dc:creator>
</cmeta:comment>
<bqs:reference xmlns:bqs="http://www.cellml.org/bqs/1.1#" rdf:parseType="Resource">
<bqs:Pubmed_id>11846603</bqs:Pubmed_id>
<bqs:JournalArticle rdf:parseType="Resource">
<dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/">
<rdf:Seq>
<rdf:li rdf:parseType="Resource">
<bqs:Person rdf:parseType="Resource">
<vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
<vCard:Family>Halloy</vCard:Family>
<vCard:Given>J</vCard:Given>
</vCard:N>
</bqs:Person>
</rdf:li>
<rdf:li rdf:parseType="Resource">
<bqs:Person rdf:parseType="Resource">
<vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
<vCard:Family>Goldbeter</vCard:Family>
<vCard:Given>A</vCard:Given>
</vCard:N>
</bqs:Person>
</rdf:li>
<rdf:li rdf:parseType="Resource">
<bqs:Person rdf:parseType="Resource">
<vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
<vCard:Family>Bernard</vCard:Family>
<vCard:Given>B</vCard:Given>
</vCard:N>
</bqs:Person>
</rdf:li>
<rdf:li rdf:parseType="Resource">
<bqs:Person rdf:parseType="Resource">
<vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource">
<vCard:Family>Loussouarn</vCard:Family>
<vCard:Given>G</vCard:Given>
</vCard:N>
</bqs:Person>
</rdf:li>
</rdf:Seq>
</dc:creator>
<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">The Follicular Automaton Model: Effect of Stochasticity and of
Synchronization of Hair Cycles</dc:title>
<bqs:Journal rdf:parseType="Resource">
<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">Journal of Theoretical Biology</dc:title>
</bqs:Journal>
<bqs:volume>214</bqs:volume>
<bqs:first_page>469</bqs:first_page>
<bqs:last_page>479</bqs:last_page>
</bqs:JournalArticle>
</bqs:reference>
</rdf:Description>
</rdf:RDF>
</model>
