Liò2012_Modelling osteomyelitis_Control Model

  public model
Model Identifier
BIOMD0000000923
Short description
Background
This work focuses on the computational modelling of osteomyelitis, a bone pathology caused by bacteria infection (mostly Staphylococcus aureus). The infection alters the RANK/RANKL/OPG signalling dynamics that regulates osteoblasts and osteoclasts behaviour in bone remodelling, i.e. the resorption and mineralization activity. The infection rapidly leads to severe bone loss, necrosis of the affected portion, and it may even spread to other parts of the body. On the other hand, osteoporosis is not a bacterial infection but similarly is a defective bone pathology arising due to imbalances in the RANK/RANKL/OPG molecular pathway, and due to the progressive weakening of bone structure.

Results
Since both osteoporosis and osteomyelitis cause loss of bone mass, we focused on comparing the dynamics of these diseases by means of computational models. Firstly, we performed meta-analysis on a gene expression data of normal, osteoporotic and osteomyelitis bone conditions. We mainly focused on RANKL/OPG signalling, the TNF and TNF receptor superfamilies and the NF-kB pathway. Using information from the gene expression data we estimated parameters for a novel model of osteoporosis and of osteomyelitis. Our models could be seen as a hybrid ODE and probabilistic verification modelling framework which aims at investigating the dynamics of the effects of the infection in bone remodelling. Finally we discuss different diagnostic estimators defined by formal verification techniques, in order to assess different bone pathologies (osteopenia, osteoporosis and osteomyelitis) in an effective way.

Conclusions
We present a modeling framework able to reproduce aspects of the different bone remodeling defective dynamics of osteomyelitis and osteoporosis. We report that the verification-based estimators are meaningful in the light of a feed forward between computational medicine and clinical bioinformatics

Model is encoded by Ruby and submitted and curated to BioModels by Ahmad Zyoud
Format
SBML (L2V4)
Related Publication
  • Modelling osteomyelitis.
  • Liò P, Paoletti N, Moni MA, Atwell K, Merelli E, Viceconti M
  • BMC bioinformatics , 1/ 2012 , Volume 13 Suppl 14 , pages: S12 , PubMed ID: 23095605
  • Computer Laboratory, Cambridge University, William Gates Building, 15 JJ Thomson Avenue, Cambridge CB3 0FD, UK.
  • BACKGROUND: This work focuses on the computational modelling of osteomyelitis, a bone pathology caused by bacteria infection (mostly Staphylococcus aureus). The infection alters the RANK/RANKL/OPG signalling dynamics that regulates osteoblasts and osteoclasts behaviour in bone remodelling, i.e. the resorption and mineralization activity. The infection rapidly leads to severe bone loss, necrosis of the affected portion, and it may even spread to other parts of the body. On the other hand, osteoporosis is not a bacterial infection but similarly is a defective bone pathology arising due to imbalances in the RANK/RANKL/OPG molecular pathway, and due to the progressive weakening of bone structure. RESULTS: Since both osteoporosis and osteomyelitis cause loss of bone mass, we focused on comparing the dynamics of these diseases by means of computational models. Firstly, we performed meta-analysis on a gene expression data of normal, osteoporotic and osteomyelitis bone conditions. We mainly focused on RANKL/OPG signalling, the TNF and TNF receptor superfamilies and the NF-kB pathway. Using information from the gene expression data we estimated parameters for a novel model of osteoporosis and of osteomyelitis. Our models could be seen as a hybrid ODE and probabilistic verification modelling framework which aims at investigating the dynamics of the effects of the infection in bone remodelling. Finally we discuss different diagnostic estimators defined by formal verification techniques, in order to assess different bone pathologies (osteopenia, osteoporosis and osteomyelitis) in an effective way. CONCLUSIONS: We present a modeling framework able to reproduce aspects of the different bone remodeling defective dynamics of osteomyelitis and osteoporosis. We report that the verification-based estimators are meaningful in the light of a feed forward between computational medicine and clinical bioinformatics.
Contributors
Submitter of the first revision: Ahmad Zyoud
Submitter of this revision: Ahmad Zyoud
Modellers: Ahmad Zyoud

Metadata information

is (2 statements)
BioModels Database BIOMD0000000923
BioModels Database MODEL2003170001

isDescribedBy (1 statement)
PubMed 23095605

isDerivedFrom (1 statement)
PubMed 14499354

hasTaxon (1 statement)
Taxonomy Homo sapiens

hasProperty (2 statements)
Mathematical Modelling Ontology Ordinary differential equation model
NCIt Osteomyelitis

occursIn (1 statement)
Brenda Tissue Ontology bone


Curation status
Curated


Tags

Connected external resources

SBGN view in Newt Editor

Name Description Size Actions

Model files

Liò2012_Modelling osteomyelitis_Control Model.xml SBML L2V4 Liò2012 Modelling osteomyelitis_Control Model_Curated 45.18 KB Preview | Download

Additional files

Liò2012_Modelling osteomyelitis_Control Model.cps COPASI version 4.27 (Build 217) Liò2012 Modelling osteomyelitis_Control Model_Curated 85.06 KB Preview | Download
Liò2012_Modelling osteomyelitis_Control Model.sedml sed-ml L1V2 Liò2012 Modelling osteomyelitis_Control Model_Curated 10.84 KB Preview | Download

  • Model originally submitted by : Ahmad Zyoud
  • Submitted: Mar 17, 2020 4:29:01 PM
  • Last Modified: Mar 17, 2020 4:29:01 PM
Revisions
  • Version: 6 public model Download this version
    • Submitted on: Mar 17, 2020 4:29:01 PM
    • Submitted by: Ahmad Zyoud
    • With comment: Automatically added model identifier BIOMD0000000923
Legends
: Variable used inside SBML models


Species
Species Initial Concentration/Amount
Osteoclasts O c

0011016
15.0 mmol
Osteoblasts O b

0011012
316.0 mmol
Bone Density z

0016250
100.0 mmol
B

C50921
1.0 mmol
Reactions
Reactions Rate Parameters
=> Osteoclasts__O_c; B, Osteoblasts__O_b compartment*(alpha1*Osteoclasts__O_c^(g11*(1+f11*B/s))*Osteoblasts__O_b^(g21*(1+f21*B/s))-beta1*Osteoclasts__O_c) s = 100.0; g11 = 1.1; beta1 = 0.2; alpha1 = 3.0; f11 = 0.005; g21 = -0.5; f21 = 0.005
=> Osteoblasts__O_b; Osteoclasts__O_c, B compartment*(alpha2*Osteoclasts__O_c^(g12*(1+f12*B/s))*Osteoblasts__O_b^(g22-f22*B/s)-beta2*Osteoblasts__O_b) s = 100.0; beta2 = 0.02; g22 = 0.0; g12 = 1.0; f22 = 0.2; f12 = 0.0; alpha2 = 4.0
=> Bone_Density__z; Osteoclasts__O_c, Osteoblasts__O_b compartment*((-k1)*piecewise(Osteoclasts__O_c-O_cbar, (Osteoclasts__O_c-O_cbar) >= 0, 0)+k2*piecewise(Osteoblasts__O_b-O_bbar, (Osteoblasts__O_b-O_bbar) >= 0, 0)) O_bbar = 177.91; O_cbar = 1.78; k2 = 6.395E-4; k1 = 0.0748
=> B compartment*(gamma_B-V)*B*log(10, s/B) s = 100.0; V = 0.007; gamma_B = 0.005
Curator's comment:
(added: 17 Mar 2020, 16:10:22, updated: 17 Mar 2020, 16:10:22)
The paper has three models, only the control model has been curated here and the control model part of the figure 4 has been reproduced here. The rest of the models needs certain parameter change to be reproduced