<?xml version="1.0"?>
<!--
This CellML file was generated on 9/04/2010 at 12:40:37 at p.m. using:
COR (0.9.31.1371)
Copyright 2002-2010 Dr Alan Garny
http://cor.physiol.ox.ac.uk/ - cor@physiol.ox.ac.uk
CellML 1.0 was used to generate this model
http://www.cellml.org/
-->
<model 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:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" xmlns:ns7="http://www.cellml.org/metadata/simulation/1.0#" name="lemaire_model_2004" cmeta:id="lemaire_model_2004" xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cellml="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
<articleinfo>
<title>Modeling the interactions between osteoblast and osteoclast activities in bone remodeling</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 represents the core model in the published paper. The model runs in both OpenCell and COR to recreate the published results. The units are consistent throughout.
</para>
</section>
<sect1 id="sec_structure">
<title>Model Structure</title>
<para>
ABSTRACT: We propose a mathematical model explaining the interactions between osteoblasts and osteoclasts, two cell types specialized in the maintenance of the bone integrity. Bone is a dynamic, living tissue whose structure and shape continuously evolves during life. It has the ability to change architecture by removal of old bone and replacement with newly formed bone in a localized process called remodeling. The model described here is based on the idea that the relative proportions of immature and mature osteoblasts control the degree of osteoclastic activity. In addition, osteoclasts control osteoblasts differentially depending on their stage of differentiation. Despite the tremendous complexity of the bone regulatory system and its fragmentary understanding, we obtain surprisingly good correlations between the model simulations and the experimental observations extracted from the literature. The model results corroborate all behaviors of the bone remodeling system that we have simulated, including the tight coupling between osteoblasts and osteoclasts, the catabolic effect induced by continuous administration of PTH, the catabolic action of RANKL, as well as its reversal by soluble antagonist OPG. The model is also able to simulate metabolic bone diseases such as estrogen deficiency, vitamin D deficiency, senescence and glucocorticoid excess. Conversely, possible routes for therapeutic interventions are tested and evaluated. Our model confirms that anti-resorptive therapies are unable to partially restore bone loss, whereas bone formation therapies yield better results. The model enables us to determine and evaluate potential therapies based on their efficacy. In particular, the model predicts that combinations of anti-resorptive and anabolic therapies provide significant benefits compared with monotherapy, especially for certain type of skeletal disease. Finally, the model clearly indicates that increasing the size of the pool of preosteoblasts is an essential ingredient for the therapeutic manipulation of bone formation. This model was conceived as the first step in a bone turnover modeling platform. These initial modeling results are extremely encouraging and lead us to proceed with additional explorations into bone turnover and skeletal remodeling.
</para>
<informalfigure float="0" id="fig_reaction_diagram">
<mediaobject>
<imageobject>
<objectinfo>
<title>model diagram</title>
</objectinfo>
<imagedata fileref="lemaire_2004_new.png"/>
</imageobject>
</mediaobject>
<caption>Schematic diagram of the basic structure of the model.</caption>
</informalfigure>
<para>
The original paper reference is cited below:
</para>
<para>
Modeling the interactions between osteoblast and osteoclast activities in bone remodeling, Vincent Lemaire, Frank L. Tobin, Larry D. Greller, Carolyn R. Cho, and Larry J. Suva, 2004, <emphasis>Journal of Theoretical Biology</emphasis>, 229, 293-309. <ulink url="http://www.ncbi.nlm.nih.gov/pubmed/15234198">PubMed ID: 15234198</ulink>
</para>
</sect1>
</article>
</documentation>
<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
<rdf:Description rdf:about="#R">
<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">R</dc:title>
<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">responding osteoblasts</dcterms:alternative>
</rdf:Description>
</rdf:RDF>
<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
<rdf:Description rdf:about="#B">
<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">B</dc:title>
<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">active osteoblasts</dcterms:alternative>
</rdf:Description>
</rdf:RDF>
<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
<rdf:Description rdf:about="#C">
<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">C</dc:title>
<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">active osteoclasts</dcterms:alternative>
</rdf:Description>
</rdf:RDF>
<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
<rdf:Description rdf:about="">
<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">
Modeling the interactions between osteoblast and osteoclast activities in bone remodeling: Main Model
</dc:title>
<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>The Bioengineering Institute</vCard:Orgunit>
</vCard:ORG>
</dc:creator>
<dcterms:created xmlns:dcterms="http://purl.org/dc/terms/" rdf:parseType="Resource">
<dcterms:W3CDTF>2007-07-26</dcterms:W3CDTF>
</dcterms:created>
<dc:publisher xmlns:dc="http://purl.org/dc/elements/1.1/">
The University of Auckland, Bioengineering Institute
</dc:publisher>
</rdf:Description>
<rdf:Description rdf:about="#lemaire_model_2004">
<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">
Lemaire et al's 2004 mathematical model of the interactions between osteoblast and osteoclast activities in bone remodelling.
</dc:title>
<cmeta:comment rdf:parseType="Resource">
<rdf:value>
This is a CellML description of Lemaire et al's 2004 mathematical model of the interactions between osteoblast and osteoclast activities in bone remodelling.
</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#">Catherine Lloyd</vCard:FN>
</dc:creator>
</cmeta:comment>
<bqs:reference xmlns:bqs="http://www.cellml.org/bqs/1.0#" rdf:parseType="Resource">
<dc:subject xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:parseType="Resource">
<bqs:subject_type>keyword</bqs:subject_type>
<rdf:value>
<rdf:Bag>
<rdf:li>osteoblast</rdf:li>
<rdf:li>osteoclast</rdf:li>
<rdf:li>osteoporosis</rdf:li>
<rdf:li>parathyroid hormone</rdf:li>
</rdf:Bag>
</rdf:value>
</dc:subject>
</bqs:reference>
<bqs:reference xmlns:bqs="http://www.cellml.org/bqs/1.0#" rdf:parseType="Resource">
<bqs:Pubmed_id>15234198</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>Lemaire</vCard:Family>
<vCard:Given>Vincent</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>Tobin</vCard:Family>
<vCard:Given>Frank</vCard:Given>
<vCard:Other>L</vCard:Other>
</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>Greller</vCard:Family>
<vCard:Given>Larry</vCard:Given>
<vCard:Other>D</vCard:Other>
</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>Cho</vCard:Family>
<vCard:Given>Carolyn</vCard:Given>
<vCard:Other>R</vCard:Other>
</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>Suva</vCard:Family>
<vCard:Given>Larry</vCard:Given>
<vCard:Other>J</vCard:Other>
</vCard:N>
</bqs:Person>
</rdf:li>
</rdf:Seq>
</dc:creator>
<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">
Modeling the interactions between osteoblast and osteoclast activities in bone remodeling
</dc:title>
<dcterms:issued xmlns:dcterms="http://purl.org/dc/terms/" rdf:parseType="Resource">
<dcterms:W3CDTF>2004-08-07</dcterms:W3CDTF>
</dcterms:issued>
<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>229</bqs:volume>
<bqs:first_page>293</bqs:first_page>
<bqs:last_page>309</bqs:last_page>
</bqs:JournalArticle>
</bqs:reference>
</rdf:Description>
</rdf:RDF> <units name="day">
<unit multiplier="86400" units="second"/>
</units>
<units name="picomole">
<unit prefix="pico" units="mole"/>
</units>
<units base_units="yes" name="cells"/>
<units name="picomole_cells">
<unit units="picomole"/>
<unit units="cells"/>
</units>
<units name="picomole_picomole_cells">
<unit units="picomole"/>
<unit exponent="-1" units="picomole_cells"/>
</units>
<units name="picomole_day_picomole_cells">
<unit units="picomole"/>
<unit exponent="-1" units="day"/>
<unit exponent="-1" units="picomole_cells"/>
</units>
<units name="picomolar">
<unit prefix="pico" units="mole"/>
<unit exponent="-1" units="litre"/>
</units>
<units name="flux">
<unit units="picomolar"/>
<unit exponent="-1" units="day"/>
</units>
<units name="first_order_rate_constant">
<unit exponent="-1" units="day"/>
</units>
<units name="second_order_rate_constant">
<unit exponent="-1" units="picomolar"/>
<unit exponent="-1" units="day"/>
</units>
<component name="environment">
<variable cmeta:id="environment_time" name="time" public_interface="out" units="day"/>
</component>
<component name="R">
<variable cmeta:id="R_R" initial_value="0" name="R" public_interface="out" units="picomolar"/>
<variable name="DR" public_interface="in" units="flux"/>
<variable name="pi_C" public_interface="in" units="dimensionless"/>
<variable name="DB" public_interface="in" units="first_order_rate_constant"/>
<variable name="time" public_interface="in" units="day"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>R</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>DR</ci>
<ci>pi_C</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<ci>DB</ci>
<ci>pi_C</ci>
</apply>
<ci>R</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="B">
<variable cmeta:id="B_B" initial_value="0" name="B" public_interface="out" units="picomolar"/>
<variable name="kB" public_interface="in" units="first_order_rate_constant"/>
<variable name="pi_C" public_interface="in" units="dimensionless"/>
<variable name="DB" public_interface="in" units="first_order_rate_constant"/>
<variable name="R" public_interface="in" units="picomolar"/>
<variable name="time" public_interface="in" units="day"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>B</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<apply>
<divide/>
<ci>DB</ci>
<ci>pi_C</ci>
</apply>
<ci>R</ci>
</apply>
<apply>
<times/>
<ci>kB</ci>
<ci>B</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="C">
<variable cmeta:id="C_C" initial_value="0" name="C" public_interface="out" units="picomolar"/>
<variable name="DC" public_interface="in" units="flux"/>
<variable name="pi_L" public_interface="in" units="dimensionless"/>
<variable name="pi_C" public_interface="in" units="dimensionless"/>
<variable name="DA" public_interface="in" units="first_order_rate_constant"/>
<variable name="time" public_interface="in" units="day"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>C</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>DC</ci>
<ci>pi_L</ci>
</apply>
<apply>
<times/>
<ci>DA</ci>
<ci>pi_C</ci>
<ci>C</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="pi_L">
<variable name="pi_L" public_interface="out" units="dimensionless"/>
<variable initial_value="1e-2" name="k1" units="second_order_rate_constant"/>
<variable initial_value="10" name="k2" units="first_order_rate_constant"/>
<variable initial_value="5.8e-4" name="k3" units="second_order_rate_constant"/>
<variable initial_value="1.7e-2" name="k4" units="first_order_rate_constant"/>
<variable initial_value="10" name="K" units="picomolar"/>
<variable initial_value="0.35" name="ko" units="first_order_rate_constant"/>
<variable initial_value="0" name="Io" units="flux"/>
<variable initial_value="0" name="IL" units="flux"/>
<variable initial_value="1e3" name="rL" units="flux"/>
<variable initial_value="2e5" name="KOP" units="picomole_day_picomole_cells"/>
<variable initial_value="3e6" name="KLP" units="picomole_picomole_cells"/>
<variable name="pi_P" public_interface="in" units="dimensionless"/>
<variable name="B" public_interface="in" units="picomolar"/>
<variable name="R" public_interface="in" units="picomolar"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>pi_L</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<ci>k3</ci>
<ci>k4</ci>
</apply>
<ci>KLP</ci>
</apply>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="picomole_picomole_cells">1</cn>
</apply>
<ci>pi_P</ci>
<ci>B</ci>
</apply>
<apply>
<plus/>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<apply>
<times/>
<ci>k3</ci>
<ci>K</ci>
</apply>
<ci>k4</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<ci>k1</ci>
<apply>
<times/>
<ci>k2</ci>
<ci>ko</ci>
</apply>
</apply>
<apply>
<plus/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<divide/>
<ci>KOP</ci>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="picomole_picomole_cells">1</cn>
</apply>
<ci>pi_P</ci>
</apply>
<ci>R</ci>
</apply>
<ci>Io</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<plus/>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>IL</ci>
<ci>rL</ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="model_parameters">
<variable initial_value="0.7" name="DA" public_interface="out" units="first_order_rate_constant"/>
<variable name="DB" public_interface="out" units="first_order_rate_constant"/>
<variable initial_value="2.1e-3" name="DC" public_interface="out" units="flux"/>
<variable initial_value="7e-4" name="DR" public_interface="out" units="flux"/>
<variable initial_value="0.189" name="kB" public_interface="out" units="first_order_rate_constant"/>
<variable name="pi_C" public_interface="out" units="dimensionless"/>
<variable name="pi_P" public_interface="out" units="dimensionless"/>
<variable initial_value="0.05" name="f0" units="dimensionless"/>
<variable initial_value="0.7" name="dB" units="first_order_rate_constant"/>
<variable initial_value="0" name="IP" units="flux"/>
<variable initial_value="86" name="kP" units="first_order_rate_constant"/>
<variable name="P" units="picomolar"/>
<variable name="P_0" units="picomolar"/>
<variable name="P_s" units="picomolar"/>
<variable initial_value="5e-3" name="C_s" units="picomolar"/>
<variable initial_value="250" name="SP" units="flux"/>
<variable initial_value="0.02" name="k5" units="second_order_rate_constant"/>
<variable initial_value="3" name="k6" units="first_order_rate_constant"/>
<variable name="C" public_interface="in" units="picomolar"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>DB</ci>
<apply>
<times/>
<ci>f0</ci>
<ci>dB</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>pi_C</ci>
<apply>
<divide/>
<apply>
<plus/>
<ci>C</ci>
<apply>
<times/>
<ci>f0</ci>
<ci>C_s</ci>
</apply>
</apply>
<apply>
<plus/>
<ci>C</ci>
<ci>C_s</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>pi_P</ci>
<apply>
<divide/>
<apply>
<plus/>
<ci>P</ci>
<ci>P_0</ci>
</apply>
<apply>
<plus/>
<ci>P</ci>
<ci>P_s</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>P</ci>
<apply>
<divide/>
<ci>IP</ci>
<ci>kP</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>P_0</ci>
<apply>
<divide/>
<ci>SP</ci>
<ci>kP</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>P_s</ci>
<apply>
<divide/>
<ci>k6</ci>
<ci>k5</ci>
</apply>
</apply>
</math>
</component>
<connection>
<map_components component_1="R" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="B" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="C" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="R" component_2="model_parameters"/>
<map_variables variable_1="DR" variable_2="DR"/>
<map_variables variable_1="DB" variable_2="DB"/>
<map_variables variable_1="pi_C" variable_2="pi_C"/>
</connection>
<connection>
<map_components component_1="R" component_2="B"/>
<map_variables variable_1="R" variable_2="R"/>
</connection>
<connection>
<map_components component_1="B" component_2="model_parameters"/>
<map_variables variable_1="kB" variable_2="kB"/>
<map_variables variable_1="DB" variable_2="DB"/>
<map_variables variable_1="pi_C" variable_2="pi_C"/>
</connection>
<connection>
<map_components component_1="C" component_2="pi_L"/>
<map_variables variable_1="pi_L" variable_2="pi_L"/>
</connection>
<connection>
<map_components component_1="C" component_2="model_parameters"/>
<map_variables variable_1="DA" variable_2="DA"/>
<map_variables variable_1="DC" variable_2="DC"/>
<map_variables variable_1="pi_C" variable_2="pi_C"/>
<map_variables variable_1="C" variable_2="C"/>
</connection>
<connection>
<map_components component_1="B" component_2="pi_L"/>
<map_variables variable_1="B" variable_2="B"/>
</connection>
<connection>
<map_components component_1="R" component_2="pi_L"/>
<map_variables variable_1="R" variable_2="R"/>
</connection>
<connection>
<map_components component_1="pi_L" component_2="model_parameters"/>
<map_variables variable_1="pi_P" variable_2="pi_P"/>
</connection>
</model>
