- Author:
- Catherine Lloyd <c.lloyd@auckland.ac.nz>
- Date:
- 2009-12-03 12:35:32+13:00
- Desc:
- Corrected the metadata at the end - added #goldbeter_model_2006, made sure all the namespaces were 1.1 and not a mix of 1.0 and 1.1, and also corrected the keywords to bqs:subject_type, value, bag.
- Permanent Source URI:
- http://models.cellml.org/workspace/goldbeter_2006/rawfile/41027328f3bceb4b172744afb366c449fd06ab0e/goldbeter_2006.cellml
<?xml version='1.0' encoding='utf-8'?>
<model name="goldbeter_model_2006" xmlns="http://www.cellml.org/cellml/1.1#" xmlns:cellml="http://www.cellml.org/cellml/1.1#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" cmeta:id="goldbeter_model_2006">
<!-- documentation -->
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
<articleinfo>
<title>A model for the dynamics of human weight cycling</title>
<author>
<firstname>Mark</firstname>
<surname>Hanna</surname>
<affiliation>
<shortaffil>Bioengineering Institute, University of Auckland</shortaffil>
</affiliation>
</author>
</articleinfo>
<section id="dec_status">
<title>Model Status</title>
<para>Runs in COR and PCEnv to reproduce published output (figure 3). No scale is mentioned in the paper, so time is scaled arbitrarily in seconds, and variables that might be represented in arbitrary units of mass or force (weight) have been left dimensionless as their normalisation has been mentioned in the paper.</para>
</section>
<sect1 id="sec_structure">
<title>Model Structure</title>
<para>The resolution to lose weight by cognitive restraint of nutritional intake often leads to repeated bouts of weight loss and regain, a phenomenon known as weight cycling or "yo-yo dieting". A simple mathematical model for weight cycling is presented. The model is based on a feedback of psychological nature by which a subject decides to reduce dietary intake once a threshold weight is exceeded. The analysis of the model indicates that sustained oscillations in body weight occur in a parameter range bounded by critical values. Only outside this range can body weight reach a stable steay state. The model provides a theoretical framework that captures key facets of weight cycling and suggests ways to control the phenomenon. The view that weight cycling represents self-sustained oscillations has indeed specific implications. In dymamical terms, to bring weight cycling to an end, parameter values should change in such a way as to induced the transition of body weight from sustained oscillations around an unstable steady state to a stable steady state. Maintaining weight under a critical value should prevent weight cycling and allow body weight to stabilise below the oscillatory range.</para>
<para>A model for the dynamics of human weight cycling, Albert Goldbeter, 2006, <emphasis>Journal of Biosciences</emphasis>, 31, 129-136. <ulink url="http://www.ncbi.nlm.nih.gov/pubmed/16595882">PubMed ID: 16595882</ulink></para>
<informalfigure float="0" id="fig_model_flowchart">
<mediaobject>
<imageobject>
<objectinfo>
<title>Model flowchart</title>
</objectinfo>
<imagedata fileref="goldbeter_2006.png"/>
</imageobject>
</mediaobject>
<caption>PQR model for weight cycling.</caption>
</informalfigure>
</sect1>
</article>
</documentation>
<!-- units -->
<units name="s">
<!-- time: second -->
<unit units="second"/>
</units>
<units name="per_s">
<!-- per time: per second -->
<unit units="s" exponent="-1"/>
</units>
<!-- components -->
<component name="environment">
<variable units="s" public_interface="out" name="time"/>
</component>
<component name="P">
<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
<rdf:Description rdf:about="#P">Body weight, expressed as an amount in excess of a basal reference value</rdf:Description>
<rdf:Description rdf:about="#a">A proportionality constant measuring metabolic efficiency</rdf:Description>
<rdf:Description rdf:about="#b">Maximum rate of weight dissipation</rdf:Description>
<rdf:Description rdf:about="#K">The value of P yielding half-maximum rate of weight dissipation</rdf:Description>
</rdf:RDF>
<variable units="s" public_interface="in" name="time"/>
<variable units="dimensionless" public_interface="out" name="P" cmeta:id="P" initial_value="0.0"/>
<variable units="dimensionless" public_interface="in" name="Q"/>
<variable units="per_s" name="a" cmeta:id="a" initial_value="0.1"/>
<variable units="per_s" name="b" cmeta:id="b" initial_value="0.1"/>
<variable units="dimensionless" name="K" cmeta:id="K" initial_value="0.2"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="P_maths">
<apply id="dP/dt_eq"><eq/>
<apply><diff/>
<bvar><ci>time</ci></bvar>
<ci>P</ci>
</apply>
<apply><minus/>
<apply><times/>
<ci>a</ci>
<ci>Q</ci>
</apply>
<apply><times/>
<ci>b</ci>
<apply><divide/>
<ci>P</ci>
<apply><plus/>
<ci>K</ci>
<ci>P</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="P_star">
<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
<rdf:Description about="#P_star">Threshold value of P</rdf:Description>
</rdf:RDF>
<variable units="dimensionless" name="P_star" cmeta:id="P_star"/>
<variable units="per_s" public_interface="in" name="V_3"/>
<variable units="per_s" public_interface="in" name="V_4"/>
<variable units="dimensionless" public_interface="in" name="K_3"/>
<variable units="dimensionless" public_interface="in" name="K_4"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="P_star_maths">
<apply id="P_star_eq"><eq/>
<ci>P_star</ci>
<apply><times/>
<apply><divide/>
<ci>V_4</ci>
<ci>V_3</ci>
</apply>
<apply><divide/>
<apply><plus/>
<cn cellml:units="dimensionless">1</cn>
<apply><times/>
<cn cellml:units="dimensionless">2</cn>
<ci>K_3</ci>
</apply>
</apply>
<apply><plus/>
<cn cellml:units="dimensionless">1</cn>
<apply><times/>
<cn cellml:units="dimensionless">2</cn>
<ci>K_4</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="Q">
<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
<rdf:Description rdf:about="#Q">Normalised excess dietary intake, expressed as an amount in excess of a basal reference value</rdf:Description>
<rdf:Description rdf:about="#V_1">Maximum rate at which Q may increase</rdf:Description>
<rdf:Description rdf:about="#V_2">Maximum rate at which Q may decrease (scaled by cognitive restraint in equation)</rdf:Description>
</rdf:RDF>
<variable units="s" public_interface="in" name="time"/>
<variable units="dimensionless" public_interface="out" name="Q" cmeta:id="Q" initial_value="0.0"/>
<variable units="dimensionless" public_interface="in" name="R"/>
<variable units="per_s" public_interface="out" name="V_1" cmeta:id="V_1" initial_value="1"/>
<variable units="per_s" public_interface="out" name="V_2" cmeta:id="V_2" initial_value="1.5"/>
<variable units="dimensionless" public_interface="out" name="K_1" initial_value="0.01"/>
<variable units="dimensionless" public_interface="out" name="K_2" initial_value="0.01"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="Q_maths">
<apply id="dQ/dt_eq"><eq/>
<apply><diff/>
<bvar><ci>time</ci></bvar>
<ci>Q</ci>
</apply>
<apply><minus/>
<apply><times/>
<ci>V_1</ci>
<apply><divide/>
<apply><minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>Q</ci>
</apply>
<apply><plus/>
<ci>K_1</ci>
<apply><minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>Q</ci>
</apply>
</apply>
</apply>
</apply>
<apply><times/>
<ci>V_2</ci>
<ci>R</ci>
<apply><divide/>
<ci>Q</ci>
<apply><plus/>
<ci>K_2</ci>
<ci>Q</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="R">
<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
<rdf:Description rdf:about="#R">Normalised degree of cognitive restraint</rdf:Description>
<rdf:Description rdf:about="#V_3">Maximum rate of increase of cognitive restraint (scaled by excess body weight in equation)</rdf:Description>
<rdf:Description rdf:about="#V_4">Maximum rate of decrease of cognitive restraint</rdf:Description>
</rdf:RDF>
<variable units="s" public_interface="in" name="time"/>
<variable units="dimensionless" public_interface="in" name="P"/>
<variable units="dimensionless" public_interface="out" name="R" cmeta:id="R" initial_value="0.0"/>
<variable units="per_s" public_interface="out" name="V_3" cmeta:id="V_3" initial_value="6"/>
<variable units="per_s" public_interface="out" name="V_4" cmeta:id="V_4" initial_value="2.5"/>
<variable units="dimensionless" public_interface="out" name="K_3" initial_value="0.01"/>
<variable units="dimensionless" public_interface="out" name="K_4" initial_value="0.01"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="R_maths">
<apply id="dR/dt_eq"><eq/>
<apply><diff/>
<bvar><ci>time</ci></bvar>
<ci>R</ci>
</apply>
<apply><minus/>
<apply><times/>
<ci>P</ci>
<ci>V_3</ci>
<apply><divide/>
<apply><minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>R</ci>
</apply>
<apply><plus/>
<ci>K_3</ci>
<apply><minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>R</ci>
</apply>
</apply>
</apply>
</apply>
<apply><times/>
<ci>V_4</ci>
<apply><divide/>
<ci>R</ci>
<apply><plus/>
<ci>K_4</ci>
<ci>R</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="R_star">
<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
<rdf:Description rdf:about="#R_star">Threshold value of R</rdf:Description>
</rdf:RDF>
<variable units="dimensionless" name="R_star" cmeta:id="R_star"/>
<variable units="per_s" public_interface="in" name="V_1"/>
<variable units="per_s" public_interface="in" name="V_2"/>
<variable units="dimensionless" public_interface="in" name="K_1"/>
<variable units="dimensionless" public_interface="in" name="K_2"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="R_star_maths">
<apply id="R_plus_eq"><eq/>
<ci>R_star</ci>
<apply><times/>
<apply><divide/>
<ci>V_1</ci>
<ci>V_2</ci>
</apply>
<apply><divide/>
<apply><plus/>
<cn cellml:units="dimensionless">1</cn>
<apply><times/>
<cn cellml:units="dimensionless">2</cn>
<ci>K_2</ci>
</apply>
</apply>
<apply><plus/>
<cn cellml:units="dimensionless">1</cn>
<apply><times/>
<cn cellml:units="dimensionless">2</cn>
<ci>K_1</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<!-- connections -->
<connection>
<map_components component_1="environment" component_2="P"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="environment" component_2="Q"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="environment" component_2="R"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="P" component_2="Q"/>
<map_variables variable_1="Q" variable_2="Q"/>
</connection>
<connection>
<map_components component_1="P" component_2="R"/>
<map_variables variable_1="P" variable_2="P"/>
</connection>
<connection>
<map_components component_1="Q" component_2="R"/>
<map_variables variable_1="R" variable_2="R"/>
</connection>
<connection>
<map_components component_1="P_star" component_2="R"/>
<map_variables variable_1="V_3" variable_2="V_3"/>
<map_variables variable_1="V_4" variable_2="V_4"/>
<map_variables variable_1="K_3" variable_2="K_3"/>
<map_variables variable_1="K_4" variable_2="K_4"/>
</connection>
<connection>
<map_components component_1="R_star" component_2="Q"/>
<map_variables variable_1="V_1" variable_2="V_1"/>
<map_variables variable_1="V_2" variable_2="V_2"/>
<map_variables variable_1="K_1" variable_2="K_1"/>
<map_variables variable_1="K_2" variable_2="K_2"/>
</connection>
<!-- metadata -->
<rdf:RDF>
<rdf:Description rdf:about="">
<dc:creator rdf:parseType="Resource">
<vCard:N rdf:parseType="Resource">
<vCard:Family>Hanna</vCard:Family>
<vCard:Given>Mark</vCard:Given>
<vCard:Other>James</vCard:Other>
</vCard:N>
<vCard:EMAIL rdf:parseType="Resource">
<rdf:value>mark@hanna.net.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>
<dcterms:created rdf:parseType="Resource">
<dcterms:W3CDTF>2009-12-02</dcterms:W3CDTF>
</dcterms:created>
</rdf:Description>
<rdf:Description rdf:about="#goldbeter_model_2006">
<dc:title>A model for the dynamics of human weight cycling</dc:title>
<cmeta:comment rdf:parseType="Resource">
<rdf:value>
This is the CellML description of Goldbeter's mathematical model of human weight cycling
</rdf:value>
<dc:creator rdf:parseType="Resource">
<vCard:FN>Mark Hanna</vCard:FN>
</dc:creator>
</cmeta:comment>
<cmeta:species>Human</cmeta:species>
<bqs:reference rdf:parseType="Resource">
<dc:subject rdf:parseType="Resource">
<bqs:subject_type>keyword</bqs:subject_type>
<rdf:value>
<rdf:Bag>
<rdf:li>Model</rdf:li>
<rdf:li>Oscillations</rdf:li>
<rdf:li>Rythms</rdf:li>
<rdf:li>Weight cycling</rdf:li>
<rdf:li>Metabolism</rdf:li>
</rdf:Bag>
</rdf:value>
</dc:subject>
</bqs:reference>
<bqs:reference rdf:parseType="Resource">
<bqs:Pubmed_id>16595882</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>Goldbeter</vCard:Family>
<vCard:Given>Albert</vCard:Given>
</vCard:N>
</bqs:Person>
</rdf:li>
</rdf:Seq>
</dc:creator>
<dc:title>A model for the dynamics of human weight cycling</dc:title>
<dcterms:issued rdf:ParseType="Resource">
<dcterms:W3CDTF>2006-03</dcterms:W3CDTF>
</dcterms:issued>
<bqs:Journal rdf:parseType="Resource">
<dc:title>Journal of Biosciences</dc:title>
</bqs:Journal>
<bqs:volume>31</bqs:volume>
<bqs:first_page>129</bqs:first_page>
<bqs:last_page>136</bqs:last_page>
</bqs:JournalArticle>
</bqs:reference>
</rdf:Description>
</rdf:RDF>
</model>
