- Author:
- Hanne <Hanne@hanne-nielsens-macbook.local>
- Date:
- 2010-06-08 10:07:19+12:00
- Desc:
- Added images in ai svg ang png format
- Permanent Source URI:
- https://models.cellml.org/workspace/komarova_2005/rawfile/9c11cd8d5743d130e234f7af67ed89bf71cb6744/komarova_2005.cellml
<?xml version='1.0' encoding='utf-8'?>
<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="komarova_2005" name="komarova_2005">
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
<articleinfo>
<title>Mathematical model of paracrine interactions between osteoclasts and osteoblasts predicts anabolic action of parathyroid hormone on bone</title>
<author>
<firstname>Catherine</firstname>
<surname>Lloyd</surname>
<affiliation>
<shortaffil>Auckland Bioengineering Institute, The University of Auckland</shortaffil>
</affiliation>
</author>
</articleinfo>
<section id="sec_status">
<title>Model Status</title>
<para>
This CellML model runs in both OpenCell and COR to recreate the published results (figure 2A part 3 the red line). The parameter values have been set to those listed in the main body of the paper with the exception that g12=1.4, k2=0.00075, and g21 increases to 0.15 for one day (day 1-day2). Other figures from the paper can be created by altering these parameter values in the CellML model.
</para>
</section>
<sect1 id="sec_structure">
<title>Model Structure</title>
<para>
ABSTRACT: To restore falling plasma calcium levels, PTH promotes calcium liberation from bone. PTH targets bone-forming cells, osteoblasts, to increase expression of the cytokine receptor activator of nuclear factor kappaB ligand (RANKL), which then stimulates osteoclastic bone resorption. Intriguingly, whereas continuous administration of PTH decreases bone mass, intermittent PTH has an anabolic effect on bone, which was proposed to arise from direct effects of PTH on osteoblastic bone formation. However, antiresorptive therapies impair the ability of PTH to increase bone mass, indicating a complex role for osteoclasts in the process. We developed a mathematical model that describes the actions of PTH at a single site of bone remodeling, where osteoclasts and osteoblasts are regulated by local autocrine and paracrine factors. It was assumed that PTH acts only to increase the production of RANKL by osteoblasts. As a result, PTH stimulated osteoclasts upon application, followed by compensatory osteoblast activation due to the coupling of osteoblasts to osteoclasts through local paracrine factors. Continuous PTH administration resulted in net bone loss, because bone resorption preceded bone formation at all times. In contrast, over a wide range of model parameters, short application of PTH resulted in a net increase in bone mass, because osteoclasts were rapidly removed upon PTH withdrawal, enabling osteoblasts to rebuild the bone. In excellent agreement with experimental findings, increase in the rate of osteoclast death abolished the anabolic effect of PTH on bone. This study presents an original concept for the regulation of bone remodeling by PTH, currently the only approved anabolic treatment for osteoporosis.
</para>
<informalfigure float="0" id="fig_reaction_diagram">
<mediaobject>
<imageobject>
<objectinfo>
<title>model diagram</title>
</objectinfo>
<imagedata fileref="komarova_2005.png"/>
</imageobject>
</mediaobject>
<caption>Schematic representation of the interactions between osteoclasts and osteoblasts included in the model.</caption>
</informalfigure>
<para>
The original paper reference is cited below:
</para>
<para>
Mathematical model of paracrine interactions between osteoclasts and osteoblasts predicts anabolic action of parathyroid hormone on bone, Komarova SV, 2005, <emphasis>Endocrinology</emphasis>, 146, 3589-3595. <ulink url="http://www.ncbi.nlm.nih.gov/pubmed/15860557">PubMed ID: 15860557</ulink>
</para>
</sect1>
</article>
</documentation>
<units name="day">
<unit units="second" multiplier="86400.0"/>
</units>
<units name="cell" base_units="yes" />
<units name="percent" base_units="yes" />
<units name="flux">
<unit units="cell"/>
<unit units="day" exponent="-1"/>
</units>
<units name="first_order_rate_constant">
<unit units="day" exponent="-1"/>
</units>
<units name="second_order_rate_constant">
<unit units="cell" exponent="-1"/>
<unit units="day" exponent="-1"/>
</units>
<units name="percent_per_cell_per_day">
<unit units="percent"/>
<unit units="cell" exponent="-1"/>
<unit units="day" exponent="-1"/>
</units>
<component name="environment">
<variable units="day" public_interface="out" name="time"/>
</component>
<component cmeta:id="x1" name="x1">
<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
<rdf:Description rdf:about="#x1">
<dc:title>x1</dc:title>
<dcterms:alternative>number of osteoclasts</dcterms:alternative>
</rdf:Description>
</rdf:RDF>
<variable units="cell" public_interface="out" name="x1" initial_value="10.06066"/>
<variable units="flux" public_interface="in" name="alpha1"/>
<variable units="first_order_rate_constant" public_interface="in" name="beta1"/>
<variable units="dimensionless" public_interface="in" name="g11"/>
<variable units="dimensionless" public_interface="in" name="g21"/>
<variable units="cell" public_interface="in" name="x2"/>
<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> x1 </ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci> alpha1 </ci>
<apply>
<power/>
<ci> x1 </ci>
<ci> g11 </ci>
</apply>
<apply>
<power/>
<ci> x2 </ci>
<ci> g21 </ci>
</apply>
</apply>
<apply>
<times/>
<ci> beta1 </ci>
<ci> x1 </ci>
</apply>
</apply>
</apply>
</math>
</component>
<component cmeta:id="x2" name="x2">
<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
<rdf:Description rdf:about="#x2">
<dc:title>x2</dc:title>
<dcterms:alternative>number of osteoblasts</dcterms:alternative>
</rdf:Description>
</rdf:RDF>
<variable units="cell" public_interface="out" name="x2" initial_value="212.132"/>
<variable units="first_order_rate_constant" public_interface="in" name="alpha2"/>
<variable units="first_order_rate_constant" public_interface="in" name="beta2"/>
<variable units="dimensionless" public_interface="in" name="g12"/>
<variable units="dimensionless" public_interface="in" name="g22"/>
<variable units="cell" public_interface="in" name="x1"/>
<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> x2 </ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci> alpha2 </ci>
<apply>
<power/>
<ci> x1 </ci>
<ci> g12 </ci>
</apply>
<apply>
<power/>
<ci> x2 </ci>
<ci> g22 </ci>
</apply>
</apply>
<apply>
<times/>
<ci> beta2 </ci>
<ci> x2 </ci>
</apply>
</apply>
</apply>
</math>
</component>
<component cmeta:id="z" name="z">
<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
<rdf:Description rdf:about="#z">
<dc:title>z</dc:title>
<dcterms:alternative>total bone mass</dcterms:alternative>
</rdf:Description>
</rdf:RDF>
<variable units="percent" public_interface="out" name="z" initial_value="100.0"/>
<variable units="percent_per_cell_per_day" public_interface="in" name="k1"/>
<variable units="percent_per_cell_per_day" public_interface="in" name="k2"/>
<variable units="cell" public_interface="in" name="y1"/>
<variable units="cell" public_interface="in" name="y2"/>
<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> z </ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci> k2 </ci>
<ci> y2 </ci>
</apply>
<apply>
<times/>
<ci> k1 </ci>
<ci> y1 </ci>
</apply>
</apply>
</apply>
</math>
</component>
<component cmeta:id="y1" name="y1">
<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
<rdf:Description rdf:about="#y1">
<dc:title>y1</dc:title>
<dcterms:alternative>number of actively resorbing bone cells</dcterms:alternative>
</rdf:Description>
</rdf:RDF>
<variable units="cell" public_interface="out" name="y1"/>
<variable units="cell" public_interface="in" name="x1"/>
<variable units="cell" public_interface="in" name="x1_bar"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci> y1 </ci>
<piecewise>
<piece>
<apply>
<minus/>
<ci> x1 </ci>
<ci> x1_bar </ci>
</apply>
<apply>
<gt/>
<ci> x1 </ci>
<ci> x1_bar </ci>
</apply>
</piece>
<otherwise>
<cn cellml:units="cell"> 0.0 </cn>
</otherwise>
</piecewise>
</apply>
</math>
</component>
<component cmeta:id="y2" name="y2">
<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
<rdf:Description rdf:about="#y2">
<dc:title>y2</dc:title>
<dcterms:alternative>number of actively forming bone cells</dcterms:alternative>
</rdf:Description>
</rdf:RDF>
<variable units="cell" public_interface="out" name="y2"/>
<variable units="cell" public_interface="in" name="x2"/>
<variable units="cell" public_interface="in" name="x2_bar"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci> y2 </ci>
<piecewise>
<piece>
<apply>
<minus/>
<ci> x2 </ci>
<ci> x2_bar </ci>
</apply>
<apply>
<gt/>
<ci> x2 </ci>
<ci> x2_bar </ci>
</apply>
</piece>
<otherwise>
<cn cellml:units="cell"> 0.0 </cn>
</otherwise>
</piecewise>
</apply>
</math>
</component>
<component cmeta:id="x1_bar" name="x1_bar">
<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
<rdf:Description rdf:about="#x1_bar">
<dc:title>x1_bar</dc:title>
<dcterms:alternative>number of osteoclasts at steady state</dcterms:alternative>
</rdf:Description>
</rdf:RDF>
<variable units="cell" public_interface="out" name="x1_bar"/>
<variable units="flux" public_interface="in" name="alpha1"/>
<variable units="first_order_rate_constant" public_interface="in" name="beta1"/>
<variable units="first_order_rate_constant" public_interface="in" name="alpha2"/>
<variable units="first_order_rate_constant" public_interface="in" name="beta2"/>
<variable units="dimensionless" public_interface="in" name="g22"/>
<variable units="dimensionless" public_interface="in" name="g21"/>
<variable units="dimensionless" public_interface="in" name="gamma"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci> x1_bar </ci>
<apply>
<times/>
<apply>
<power/>
<apply>
<divide/>
<ci> beta1 </ci>
<ci> alpha1 </ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<ci> g22 </ci>
</apply>
<ci> gamma </ci>
</apply>
</apply>
<apply>
<power/>
<apply>
<divide/>
<ci> beta2 </ci>
<ci> alpha2 </ci>
</apply>
<apply>
<divide/>
<ci> g21 </ci>
<ci> gamma </ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component cmeta:id="x2_bar" name="x2_bar">
<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
<rdf:Description rdf:about="#x2_bar">
<dc:title>x2_bar</dc:title>
<dcterms:alternative>number of osteoblasts at steady state</dcterms:alternative>
</rdf:Description>
</rdf:RDF>
<variable units="cell" public_interface="out" name="x2_bar"/>
<variable units="flux" public_interface="in" name="alpha1"/>
<variable units="first_order_rate_constant" public_interface="in" name="beta1"/>
<variable units="first_order_rate_constant" public_interface="in" name="alpha2"/>
<variable units="first_order_rate_constant" public_interface="in" name="beta2"/>
<variable units="dimensionless" public_interface="in" name="g11"/>
<variable units="dimensionless" public_interface="in" name="g12"/>
<variable units="dimensionless" public_interface="in" name="gamma"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci> x2_bar </ci>
<apply>
<times/>
<apply>
<power/>
<apply>
<divide/>
<ci> beta1 </ci>
<ci> alpha1 </ci>
</apply>
<apply>
<divide/>
<ci> g12 </ci>
<ci> gamma </ci>
</apply>
</apply>
<apply>
<power/>
<apply>
<divide/>
<ci> beta2 </ci>
<ci> alpha2 </ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<ci> g11 </ci>
</apply>
<ci> gamma </ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="model_parameters">
<variable units="flux" public_interface="out" name="alpha1" initial_value="3"/>
<variable units="first_order_rate_constant" public_interface="out" name="beta1" initial_value="0.2"/>
<variable units="first_order_rate_constant" public_interface="out" name="alpha2" initial_value="4"/>
<variable units="first_order_rate_constant" public_interface="out" name="beta2" initial_value="0.02"/>
<variable units="percent_per_cell_per_day" public_interface="out" name="k1" initial_value="0.24"/>
<variable units="percent_per_cell_per_day" public_interface="out" name="k2" initial_value="0.0017"/>
<variable units="dimensionless" public_interface="out" name="g11" initial_value="0.5"/>
<variable units="dimensionless" public_interface="out" name="g21"/>
<variable units="dimensionless" public_interface="out" name="g12" initial_value="1"/>
<variable units="dimensionless" public_interface="out" name="g22" initial_value="0"/>
<variable units="dimensionless" public_interface="out" name="gamma"/>
<variable units="day" public_interface="in" name="time"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci> gamma </ci>
<apply>
<minus/>
<apply>
<times/>
<ci> g12 </ci>
<ci> g21 </ci>
</apply>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<ci> g11 </ci>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<ci> g22 </ci>
</apply>
</apply>
</apply>
</apply>
</math>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci> g21 </ci>
<piecewise>
<piece>
<cn cellml:units="dimensionless"> 0.15 </cn>
<apply>
<and/>
<apply>
<geq/>
<ci> time </ci>
<cn cellml:units="day"> 1.0 </cn>
</apply>
<apply>
<lt/>
<ci> time </ci>
<cn cellml:units="day"> 2.00 </cn>
</apply>
</apply>
</piece>
<otherwise>
<cn cellml:units="dimensionless"> -0.5 </cn>
</otherwise>
</piecewise>
</apply>
</math>
</component>
<connection>
<map_components component_2="environment" component_1="x1"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="x2"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="z"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="model_parameters"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="x1" component_1="x2"/>
<map_variables variable_2="x1" variable_1="x1"/>
<map_variables variable_2="x2" variable_1="x2"/>
</connection>
<connection>
<map_components component_2="model_parameters" component_1="x1"/>
<map_variables variable_2="alpha1" variable_1="alpha1"/>
<map_variables variable_2="beta1" variable_1="beta1"/>
<map_variables variable_2="g11" variable_1="g11"/>
<map_variables variable_2="g21" variable_1="g21"/>
</connection>
<connection>
<map_components component_2="model_parameters" component_1="x2"/>
<map_variables variable_2="alpha2" variable_1="alpha2"/>
<map_variables variable_2="beta2" variable_1="beta2"/>
<map_variables variable_2="g12" variable_1="g12"/>
<map_variables variable_2="g22" variable_1="g22"/>
</connection>
<connection>
<map_components component_2="z" component_1="y1"/>
<map_variables variable_2="y1" variable_1="y1"/>
</connection>
<connection>
<map_components component_2="z" component_1="y2"/>
<map_variables variable_2="y2" variable_1="y2"/>
</connection>
<connection>
<map_components component_2="model_parameters" component_1="z"/>
<map_variables variable_2="k1" variable_1="k1"/>
<map_variables variable_2="k2" variable_1="k2"/>
</connection>
<connection>
<map_components component_2="y1" component_1="x1"/>
<map_variables variable_2="x1" variable_1="x1"/>
</connection>
<connection>
<map_components component_2="y1" component_1="x1_bar"/>
<map_variables variable_2="x1_bar" variable_1="x1_bar"/>
</connection>
<connection>
<map_components component_2="y2" component_1="x2"/>
<map_variables variable_2="x2" variable_1="x2"/>
</connection>
<connection>
<map_components component_2="y2" component_1="x2_bar"/>
<map_variables variable_2="x2_bar" variable_1="x2_bar"/>
</connection>
<connection>
<map_components component_2="model_parameters" component_1="x1_bar"/>
<map_variables variable_2="beta1" variable_1="beta1"/>
<map_variables variable_2="alpha1" variable_1="alpha1"/>
<map_variables variable_2="beta2" variable_1="beta2"/>
<map_variables variable_2="alpha2" variable_1="alpha2"/>
<map_variables variable_2="g21" variable_1="g21"/>
<map_variables variable_2="g22" variable_1="g22"/>
<map_variables variable_2="gamma" variable_1="gamma"/>
</connection>
<connection>
<map_components component_2="model_parameters" component_1="x2_bar"/>
<map_variables variable_2="beta1" variable_1="beta1"/>
<map_variables variable_2="alpha1" variable_1="alpha1"/>
<map_variables variable_2="beta2" variable_1="beta2"/>
<map_variables variable_2="alpha2" variable_1="alpha2"/>
<map_variables variable_2="g11" variable_1="g11"/>
<map_variables variable_2="g12" variable_1="g12"/>
<map_variables variable_2="gamma" variable_1="gamma"/>
</connection>
<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">
<rdf:Description rdf:about="">
<dc:creator rdf:parseType="Resource">
<vCard:N rdf:parseType="Resource">
<vCard:Family>Lloyd</vCard:Family>
<vCard:Given>Catherine</vCard:Given>
<vCard:Other>May</vCard:Other>
</vCard:N>
<vCard:EMAIL rdf:parseType="Resource">
<rdf:value>c.lloyd@auckland.ac.nz</rdf:value>
<rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
</vCard:EMAIL>
<vCard:ORG rdf:parseType="Resource">
<vCard:Orgname>The University of Auckland</vCard:Orgname>
<vCard:Orgunit>Auckland Bioengineering Institute</vCard:Orgunit>
</vCard:ORG>
</dc:creator>
</rdf:Description>
<rdf:Description rdf:about="#komarova_2005">
<bqs:reference rdf:parseType="Resource">
<dc:subject rdf:parseType="Resource">
<bqs:subject_type>keyword</bqs:subject_type>
<rdf:value>
<rdf:Bag>
<rdf:li>osteoclast</rdf:li>
<rdf:li>osteoblast</rdf:li>
<rdf:li>bone</rdf:li>
<rdf:li>parathyroid hormone</rdf:li>
<rdf:li>endocrine</rdf:li>
</rdf:Bag>
</rdf:value>
</dc:subject>
</bqs:reference>
<bqs:reference rdf:parseType="Resource">
<bqs:Pubmed_id>15860557</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>Komarova</vCard:Family>
<vCard:Given>S</vCard:Given>
<vCard:Other>V</vCard:Other>
</vCard:N>
</bqs:Person>
</rdf:li>
</rdf:Seq>
</dc:creator>
<dc:title>
Mathematical model of paracrine interactions between osteoclasts and osteoblasts predicts anabolic action of parathyroid hormone on bone
</dc:title>
<dcterms:issued rdf:parseType="Resource">
<dcterms:W3CDTF>2005-08</dcterms:W3CDTF>
</dcterms:issued>
<bqs:Journal rdf:parseType="Resource">
<dc:title>Endocrinology</dc:title>
</bqs:Journal>
<bqs:volume>146</bqs:volume>
<bqs:first_page>3589</bqs:first_page>
<bqs:last_page>3595</bqs:last_page>
</bqs:JournalArticle>
</bqs:reference>
</rdf:Description>
</rdf:RDF>
</model>
