- Author:
- Catherine Lloyd <c.lloyd@auckland.ac.nz>
- Date:
- 2010-03-31 09:28:07+13:00
- Desc:
- Deleted duplicate file (i and j were the same). Added the orignal core model back to the workspace. Fixed broken RDF. Updated documentation. Added a curation comment to each model file.
- Permanent Source URI:
- https://models.cellml.org/workspace/lemaire_tobin_greller_cho_suva_2004/rawfile/f86ab86177e48dc60acf8bd0e7d88b2ee59809bb/lemaire_tobin_greller_cho_suva_2004.cellml
<?xml version='1.0' encoding='utf-8'?>
<!-- FILE : lemaire_model_2004.xml
CREATED : 25th July 2007
LAST MODIFIED : 25th July 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 Lemaire et al's 2004 mathematical model of the interactions between osteoblast and osteoclast activities in bone remodelling.
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#" xmlns:ns7="http://www.cellml.org/metadata/simulation/1.0#" name="lemaire_model_2004" cmeta:id="lemaire_model_2004">
<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. There is a unit inconsistency in the equation for pi_L.
</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>
<units name="day">
<unit units="second" multiplier="86400.0"/>
</units>
<units name="picomole">
<unit units="mole" prefix="pico"/>
</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 units="picomole_cells" exponent="-1"/>
</units>
<units name="picomole_day_picomole_cells">
<unit units="picomole"/>
<unit units="day" exponent="-1"/>
<unit units="picomole_cells" exponent="-1"/>
</units>
<units name="picomolar">
<unit units="mole" prefix="pico"/>
<unit units="litre" exponent="-1"/>
</units>
<units name="flux">
<unit units="picomolar"/>
<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="picomolar" exponent="-1"/>
<unit units="day" exponent="-1"/>
</units>
<component name="environment">
<variable units="day" public_interface="out" name="time" cmeta:id="environment_time"/>
</component>
<component name="R" cmeta:id="R">
<rdf:RDF 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/">
<rdf:Description rdf:about="#R">
<dc:title>R</dc:title>
<dcterms:alternative>responding osteoblasts</dcterms:alternative>
</rdf:Description>
</rdf:RDF>
<variable name="R" public_interface="out" initial_value="0.00" 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" cmeta:id="B">
<rdf:RDF 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/">
<rdf:Description rdf:about="#B">
<dc:title>B</dc:title>
<dcterms:alternative>active osteoblasts</dcterms:alternative>
</rdf:Description>
</rdf:RDF>
<variable name="B" public_interface="out" initial_value="0.00" 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" cmeta:id="C">
<rdf:RDF 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/">
<rdf:Description rdf:about="#C">
<dc:title>C</dc:title>
<dcterms:alternative>active osteoclasts</dcterms:alternative>
</rdf:Description>
</rdf:RDF>
<variable name="C" public_interface="out" initial_value="0.00" 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 name="k1" initial_value="1E-2" units="second_order_rate_constant"/>
<variable name="k2" initial_value="10.0" units="first_order_rate_constant"/>
<variable name="k3" initial_value="5.8E-4" units="second_order_rate_constant"/>
<variable name="k4" initial_value="1.7E-2" units="first_order_rate_constant"/>
<variable name="K" initial_value="10.0" units="picomolar"/>
<variable name="ko" initial_value="0.35" units="first_order_rate_constant"/>
<variable name="Io" initial_value="0.0" units="flux"/>
<variable name="IL" initial_value="0.0" units="flux"/>
<variable name="rL" initial_value="1E3" units="flux"/>
<variable name="KOP" initial_value="2E5" units="picomole_day_picomole_cells"/>
<variable name="KLP" initial_value="3E6" 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/>
<ci> k3 </ci>
<ci> k4 </ci>
</apply>
<apply><divide/>
<apply><times/>
<ci> KLP </ci>
<ci> pi_P </ci>
<ci> B </ci>
</apply>
<apply><plus/>
<cn cellml:units="dimensionless"> 1.0 </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/>
<ci> KOP </ci>
<ci> pi_P </ci>
</apply>
<ci> R </ci>
</apply>
<ci> Io </ci>
</apply>
</apply>
</apply>
</apply>
<apply><plus/>
<cn cellml:units="dimensionless"> 1.0 </cn>
<apply><divide/>
<ci> IL </ci>
<ci> rL </ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="model_parameters">
<variable name="DA" public_interface="out" initial_value="0.7" units="first_order_rate_constant"/>
<variable name="DB" public_interface="out" units="first_order_rate_constant"/>
<variable name="DC" public_interface="out" initial_value="2.1E-3" units="flux"/>
<variable name="DR" public_interface="out" initial_value="7E-4" units="flux"/>
<variable name="kB" public_interface="out" initial_value="0.189" units="first_order_rate_constant"/>
<variable name="pi_C" public_interface="out" units="dimensionless"/>
<variable name="pi_P" public_interface="out" units="dimensionless"/>
<variable name="f0" initial_value="0.05" units="dimensionless"/>
<variable name="dB" initial_value="0.7" units="first_order_rate_constant"/>
<variable name="IP" initial_value="0.0" units="flux"/>
<variable name="kP" initial_value="86.0" units="first_order_rate_constant"/>
<variable name="P" units="picomolar"/>
<variable name="P_0" units="picomolar"/>
<variable name="P_s" units="picomolar"/>
<variable name="C_s" initial_value="5E-3" units="picomolar"/>
<variable name="SP" initial_value="250.0" units="flux"/>
<variable name="k5" initial_value="0.02" units="second_order_rate_constant"/>
<variable name="k6" initial_value="3.0" 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_2="environment" component_1="R"/>
<map_variables variable_2="time" variable_1="time"/>
</connection>
<connection>
<map_components component_2="environment" component_1="B"/>
<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="model_parameters" component_1="R"/>
<map_variables variable_2="DR" variable_1="DR"/>
<map_variables variable_2="DB" variable_1="DB"/>
<map_variables variable_2="pi_C" variable_1="pi_C"/>
</connection>
<connection>
<map_components component_2="B" component_1="R"/>
<map_variables variable_2="R" variable_1="R"/>
</connection>
<connection>
<map_components component_2="model_parameters" component_1="B"/>
<map_variables variable_2="kB" variable_1="kB"/>
<map_variables variable_2="DB" variable_1="DB"/>
<map_variables variable_2="pi_C" variable_1="pi_C"/>
</connection>
<connection>
<map_components component_2="pi_L" component_1="C"/>
<map_variables variable_2="pi_L" variable_1="pi_L"/>
</connection>
<connection>
<map_components component_2="model_parameters" component_1="C"/>
<map_variables variable_2="DA" variable_1="DA"/>
<map_variables variable_2="DC" variable_1="DC"/>
<map_variables variable_2="pi_C" variable_1="pi_C"/>
<map_variables variable_2="C" variable_1="C"/>
</connection>
<connection>
<map_components component_2="pi_L" component_1="B"/>
<map_variables variable_2="B" variable_1="B"/>
</connection>
<connection>
<map_components component_2="pi_L" component_1="R"/>
<map_variables variable_2="R" variable_1="R"/>
</connection>
<connection>
<map_components component_2="model_parameters" component_1="pi_L"/>
<map_variables variable_2="pi_P" variable_1="pi_P"/>
</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#">
<!--
The following RDF block contains metadata that applies to this document
as a whole, as indicated by the empty about attribute on the
<rdf:Description> element.
-->
<rdf:Description rdf:about="">
<!--
The Model Builder Metadata. The Dublin Core "creator" element is used
to indicate the person who translated the model into CellML.
-->
<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>The Bioengineering Institute</vCard:Orgunit>
</vCard:ORG>
</dc:creator>
<!--
The Creation Date metadata. This is the date on which the model
was translated into CellML.
-->
<dcterms:created rdf:parseType="Resource">
<dcterms:W3CDTF>2007-07-26</dcterms:W3CDTF>
</dcterms:created>
<!--
The Modification History metadata. This lists the changes that have been
made to the document, who made the changes, and when they were made.
-->
<!-- The Publisher metadata. -->
<dc:publisher>
The University of Auckland, Bioengineering Institute
</dc:publisher>
</rdf:Description>
<!--
The following metadata refers to the model itself, as indicated by the
reference to the ID "lemaire_model_2004", which is declared
on the <model> element.
-->
<rdf:Description rdf:about="#lemaire_model_2004">
<!-- A human readable name for the model. -->
<dc:title>
Lemaire et al's 2004 mathematical model of the interactions between osteoblast and osteoclast activities in bone remodelling.
</dc:title>
<!-- A comment regarding the model. -->
<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>
<!-- The creator of the comment. -->
<dc:creator rdf:parseType="Resource">
<vCard:FN>Catherine Lloyd</vCard:FN>
</dc:creator>
</cmeta:comment>
<!-- Keyword(s) -->
<bqs:reference rdf:parseType="Resource">
<dc:subject 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>
<!--
The CellML Metadata Specification recommends that bibliographic metadata
is used to provide information about the original model reference. The
"identifier" attribute on the "BibliographicReference" class provides an
elegant way to identify a cited reference using a database identifier
such as Pubmed. All associated data such as author, journal title,
date, etc can be looked up on the database.
-->
<bqs:reference rdf:parseType="Resource">
<bqs:Pubmed_id>15234198</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>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 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 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 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 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>
Modeling the interactions between osteoblast and osteoclast activities in bone remodeling
</dc:title>
<dcterms:issued rdf:parseType="Resource">
<dcterms:W3CDTF>2004-08-07</dcterms:W3CDTF>
</dcterms:issued>
<bqs:Journal rdf:parseType="Resource">
<dc:title>
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>
</model>
