Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2012;8(10):e1002724.
doi: 10.1371/journal.pcbi.1002724. Epub 2012 Oct 11.

MOSAIC: a multiscale model of osteogenesis and sprouting angiogenesis with lateral inhibition of endothelial cells

Affiliations

MOSAIC: a multiscale model of osteogenesis and sprouting angiogenesis with lateral inhibition of endothelial cells

Aurélie Carlier et al. PLoS Comput Biol. 2012.

Erratum in

  • PLoS Comput Biol. 2013 Mar;9(3). doi: 10.1371/annotation/38264a13-d4b5-49cd-b54e-47330bb19fe9

Abstract

The healing of a fracture depends largely on the development of a new blood vessel network (angiogenesis) in the callus. During angiogenesis tip cells lead the developing sprout in response to extracellular signals, amongst which vascular endothelial growth factor (VEGF) is critical. In order to ensure a correct development of the vasculature, the balance between stalk and tip cell phenotypes must be tightly controlled, which is primarily achieved by the Dll4-Notch1 signaling pathway. This study presents a novel multiscale model of osteogenesis and sprouting angiogenesis, incorporating lateral inhibition of endothelial cells (further denoted MOSAIC model) through Dll4-Notch1 signaling, and applies it to fracture healing. The MOSAIC model correctly predicted the bone regeneration process and recapitulated many experimentally observed aspects of tip cell selection: the salt and pepper pattern seen for cell fates, an increased tip cell density due to the loss of Dll4 and an excessive number of tip cells in high VEGF environments. When VEGF concentration was even further increased, the MOSAIC model predicted the absence of a vascular network and fracture healing, thereby leading to a non-union, which is a direct consequence of the mutual inhibition of neighboring cells through Dll4-Notch1 signaling. This result was not retrieved for a more phenomenological model that only considers extracellular signals for tip cell migration, which illustrates the importance of implementing the actual signaling pathway rather than phenomenological rules. Finally, the MOSAIC model demonstrated the importance of a proper criterion for tip cell selection and the need for experimental data to further explore this. In conclusion, this study demonstrates that the MOSAIC model creates enhanced capabilities for investigating the influence of molecular mechanisms on angiogenesis and its relation to bone formation in a more mechanistic way and across different time and spatial scales.

PubMed Disclaimer

Conflict of interest statement

The authors have declared that no competing interests exist.

Figures

Figure 1
Figure 1. Schematic overview of the multiscale model and a flowchart of its numerical implementation.
(A) Scale separation map indicating schematically the modeled processes at different spatial and temporal scales. Intracellular variables are further explained in Table 1 and govern endothelial cell (EC) behavior. Cell types that are considered at the tissue scale (MSCs: mesenchymal stem cells, CHs: chondrocytes, OBs: osteoblasts, FBs: fibroblasts) can migrate (only MSCs and FBs), proliferate (circular arrows), differentiate (vertical arrows) and produce growth factors (gb: osteogenic growth factor concentration, gc: chondrogenic growth factor concentration, gv: angiogenic growth factor concentration) and extracellular matrix (mf: fibrous tissue density, mb: bone density, mc: cartilage density, m: total tissue density). Blood vessels serve as an oxygen source (n: oxygen concentration). Variables next to an arrow indicate their mediating role for a certain tissue level process. (B) Flowchart of the numerical implementation of the MOSAIC model (formula image: vector of the continuous variables, t: time, δt: time step of the inner loop, Δt: time step of the outer loop, cv: endothelial cells).
Figure 2
Figure 2. Conceptual representation of the Dll4 distribution across the cell membranes of endothelial cells.
The dark grey ECs represent the original branch whereas Ti and Si are part of a newly formed sprout (T = tip cell, S = stalk cell). The amount of Dll4 was arbitrarily chosen, a maximal amount of 25 000 ligands per EC (for an EC with a length of 25 µm) was calculated by Bentley et al. , estimated from Liu et al. .
Figure 3
Figure 3. Geometrical domain and boundary conditions.
(left) Geometrical domain deduced from the real callus geometry at postfracture week 3 ; 1 periosteal callus; 2 intercortical callus; 3 endosteal callus; 4 cortical bone. (right) Boundary conditions. The mesenchymal stem cells (cm) and fibroblasts (cf) are released from the periosteum, the surrounding soft tissues and the bone marrow . The chondrogenic growth factors (gc) are released from the degrading bone ends whereas the cortex is modeled to be the source of osteogenic growth factors (gb) . cv indicates the initial position of the ECs.
Figure 4
Figure 4. In silico and in vivo evolution of normal fracture healing.
Temporal evolution of the bone, cartilage and fibrous tissue fractions (%) in the periosteal, intercortical and endosteal callus as predicted using the hybrid model of Peiffer et al. and the newly developed multiscale model and as measured by Harrison et al. .
Figure 5
Figure 5. Image of the amount of VEGFR-2 per EC for different Dll4 expression and Notch activity levels at post fracture day 19.
The figure is focused on the periosteal callus (area 1 in Figure 3). Remark that the tip cells have a lot of VEGFR-2 (dark brown) whereas the following stalk cells are inhibited, giving rise to a low number of VEGFR-2 (dark blue). (A) Influence of Dll4 levels. (B) Influence of Notch1 activity levels. Left: heterozygous knockout, middle: standard condition, right: overexpression.
Figure 6
Figure 6. Temporal evolution of the average amount of VEGFR-2 on ECs in the callus for different levels of VEGF addition (corresponding to Figure 7B ).
Evolution of the average VEGFR-2 concentration (day 0 corresponds to the time of fracture). The threshold line represents the first requirement that needs to be fulfilled to obtain the tip cell phenotype, i.e. V>Vmax/2. Remark that in very high VEGF concentrations (+10%) the average receptor concentration is far below the threshold.
Figure 7
Figure 7. Vasculature at 35 days post fracture for different conditions.
(A) Variation in the amount of the decoy receptor VEGFR-1 (i: Vsink = 100% (standard condition), ii: Vsink = 45%, iii: Vsink = 9%, iv: Vsink = 0.9%); (B) Variation in the amount of VEGF addition in the multiscale model with the standard tip cell selection criterion (based on V; Equation 6) (i: 0% (standard condition), ii: 0.1%, iii: 2%, iv: 10%); (C) Variation in the amount of VEGF addition in the multiscale model with an altered tip cell selection criterion (based on V′) (i: 0%, ii: 0.1%, iii: 2%, iv: 10%); (D) Variation in the amount of VEGF addition in the hybrid model (i: 0%, ii: 0.1%, iii: 2%, iv: 10%).
Figure 8
Figure 8. In silico evolution of fracture healing for various conditions.
Temporal evolution of the bone, cartilage and fibrous tissue fractions (%) in the periosteal, intercortical and endosteal callus as predicted by the MOSAIC model. (full: normal fracture healing; dashed: Dll4 overexpression (δ = 12.64, Figure 5A); dotted: VEGF-addition (+ 0.1%, Figure 7Bii).
Figure 9
Figure 9. Amount of VEGFR-2 receptors on every endothelial cell.
(A) standard condition; (B) addition of 2% VEGF.
Figure 10
Figure 10. Schematic representation of the Dll4-Notch pathway in high VEGF environments for two neighboring ECs.
Remark that although the VEGFR-2 receptor is down-regulated, the VEGF concentration is high enough to compensate for this effect resulting in high V′ levels (Equation 1).

Similar articles

Cited by

References

    1. Taguchi K, Ogawa R, Migita M, Hanawa H, Ito H, et al. (2005) The role of bone marrow-derived cells in bone fracture repair in a green fluorescent protein chimeric mouse model. Biochem Biophys Res Commun 331: 31–36. - PubMed
    1. Carmeliet P, De Smet F, Loges S, Mazzone M (2009) Branching morphogenesis and antiangiogenesis candidates: tip cells lead the way. Nat Rev Clin Oncol 6: 315–326. - PubMed
    1. Bentley K, Gerhardt H, Bates PA (2008) Agent-based simulation of notch-mediated tip cell selection in angiogenic sprout initialisation. J Theor Biol 250: 25–36 S0022-5193(07)00443-2 [pii];10.1016/j.jtbi.2007.09.015 [doi]. - PubMed
    1. De Smet F, Segura I, De Bock K, Hohensinner PJ, Carmeliet P (2009) Mechanisms of Vessel Branching Filopodia on Endothelial Tip Cells Lead the Way. Arterioscler Thromb Vasc Biol 29: 639–649. - PubMed
    1. Hellstrom M, Phng LK, Hofmann JJ, Wallgard E, Coultas L, et al. (2007) Dll4 signalling through Notch1 regulates formation of tip cells during angiogenesis. Nature 445: 776–780. - PubMed

Publication types

MeSH terms

Substances

Grants and funding

These authors were supported by the following funding: AC: PhD fellow of the Research Foundation Flanders, LG: Funded by the Special Research Fund of the University of Liège (FRS.D-10/20), KB: Funded by the Artemis network grant, Foundation Leducq, PC: This work of PC is supported by the Belgian Science Policy (IAP #P6-30); and by long-term structural funding - Methusalem funding by the Flemish Government. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.