Author Contributions
Conceptualization, L.S.; methodology, Y.H.; software, Y.H. and N.M.M.; validation, Y.H.; formal analysis, Y.H.; investigation, Y.H. and N.M.M.; resources, L.S.; data curation, Y.H.; writing—original draft preparation, Y.H., N.M.M. and L.S.; writing—review and editing, Y.H., N.M.M. and L.S.; visualization, Y.H.; supervision, L.S.; project administration, L.S.; funding acquisition, L.S. All authors have read and agreed to the published version of the manuscript.
Figure 1.
Interaction network of the tumor microenvironment in osteosarcoma. Blue arrows show activation and proliferation while inhibitions are indicated by red arrows.
Figure 1.
Interaction network of the tumor microenvironment in osteosarcoma. Blue arrows show activation and proliferation while inhibitions are indicated by red arrows.
Figure 2.
Dynamics of variables
C and
D. Part (
a) shows the value of
C at time
for the first group of osteosarcoma patients provided in [
12]. The agreement of the curve in (
b) with the green curve (results of the first group of tumors) from (
c) shows that the dynamics of
C at any point in the domain match the solution from the ODE paper. Sub-figure (
d–
f) show the dynamics of
D.
Figure 2.
Dynamics of variables
C and
D. Part (
a) shows the value of
C at time
for the first group of osteosarcoma patients provided in [
12]. The agreement of the curve in (
b) with the green curve (results of the first group of tumors) from (
c) shows that the dynamics of
C at any point in the domain match the solution from the ODE paper. Sub-figure (
d–
f) show the dynamics of
D.
Figure 3.
Evolution of the domain. The original domain is plotted in (a). At , the arrows, showing the vector field , are pointing inward, since V is decreasing in (b) whereas in (c), they are pointing outward since V is increasing, which means that the domain is growing at . Part (d–f) show the size of the domain at time , 200, and 400, respectively.
Figure 3.
Evolution of the domain. The original domain is plotted in (a). At , the arrows, showing the vector field , are pointing inward, since V is decreasing in (b) whereas in (c), they are pointing outward since V is increasing, which means that the domain is growing at . Part (d–f) show the size of the domain at time , 200, and 400, respectively.
Figure 4.
The diameter of the tumor is estimated from different sources. The red line, blue line, and black line represent the estimation from taking a square root, a linear model, and applying the displacement vector. It can be seen that the estimate obtained from applying the displacement vector is sound.
Figure 4.
The diameter of the tumor is estimated from different sources. The red line, blue line, and black line represent the estimation from taking a square root, a linear model, and applying the displacement vector. It can be seen that the estimate obtained from applying the displacement vector is sound.
Figure 5.
Illustrations for the case of more M and in the middle initially. Sub-figure (a,c) shows the initial value of and M through the domain, respectively; (d) shows at , which marks a drastic change in both the profile and values; the dynamics of M and over the whole time interval is shown in (b,e), respectively; the cancer cell concentrations at and are shown in (f,g), and the dynamics of C is shown in (h).
Figure 5.
Illustrations for the case of more M and in the middle initially. Sub-figure (a,c) shows the initial value of and M through the domain, respectively; (d) shows at , which marks a drastic change in both the profile and values; the dynamics of M and over the whole time interval is shown in (b,e), respectively; the cancer cell concentrations at and are shown in (f,g), and the dynamics of C is shown in (h).
Figure 6.
Illustrations for the case of using Robin boundary condition for 5 immune cell types. Sub-figure (a–d) show the value of , M, , and C at ; Sub-figure (e–h) show the maximum and minimum in the whole interval for , M, , and C, respectively. Using the same boundary condition for 5 immune cell types results in different profiles for these cell types.
Figure 6.
Illustrations for the case of using Robin boundary condition for 5 immune cell types. Sub-figure (a–d) show the value of , M, , and C at ; Sub-figure (e–h) show the maximum and minimum in the whole interval for , M, , and C, respectively. Using the same boundary condition for 5 immune cell types results in different profiles for these cell types.
Figure 7.
Illustrations for the case of a source of M in the middle. Sub-figure (a,b) show the value of M at and , respectively; (c) shows the dynamics of the maximum and minimum of M, from where we can see that the source has introduced a significant amount of increase in the middle; the cancer cell concentration at and are shown in (d,e); the value of the maximum and minimum of C over the whole time interval is shown in (f). These figures suggest that the cancer cells will be at the place where there are more macrophages, yet the concentration still reaches a steady state.
Figure 7.
Illustrations for the case of a source of M in the middle. Sub-figure (a,b) show the value of M at and , respectively; (c) shows the dynamics of the maximum and minimum of M, from where we can see that the source has introduced a significant amount of increase in the middle; the cancer cell concentration at and are shown in (d,e); the value of the maximum and minimum of C over the whole time interval is shown in (f). These figures suggest that the cancer cells will be at the place where there are more macrophages, yet the concentration still reaches a steady state.
Figure 8.
Illustrations for the case of a source of in the middle. Part (a,b) show the value of at and , respectively; (c) shows the dynamics of the maximum and minimum of , depicting that the value of in the middle is constantly higher than it on the boundary by a considerable amount; the cancer cell concentration at and are shown in (d,e); the value of the maximum and minimum of C over the whole time interval is shown in (f), where the growth of cancer cells in the middle are inhibited strongly by the higher concentration of .
Figure 8.
Illustrations for the case of a source of in the middle. Part (a,b) show the value of at and , respectively; (c) shows the dynamics of the maximum and minimum of , depicting that the value of in the middle is constantly higher than it on the boundary by a considerable amount; the cancer cell concentration at and are shown in (d,e); the value of the maximum and minimum of C over the whole time interval is shown in (f), where the growth of cancer cells in the middle are inhibited strongly by the higher concentration of .
Figure 9.
Illustrations for the case of a constant source of on the boundary. Su-figures (a,b) shows the value of at and , respectively; (c) shows the dynamics of the maximum and minimum of , depicting that the value of on the boundary is constantly higher than it is elsewhere; the cancer cell concentration at and are shown in (d,e); the value of the maximum and minimum of C over the whole time interval is shown in (f), where the growth of cancer cells, especially near the boundary, are inhibited strongly by the higher concentration of .
Figure 9.
Illustrations for the case of a constant source of on the boundary. Su-figures (a,b) shows the value of at and , respectively; (c) shows the dynamics of the maximum and minimum of , depicting that the value of on the boundary is constantly higher than it is elsewhere; the cancer cell concentration at and are shown in (d,e); the value of the maximum and minimum of C over the whole time interval is shown in (f), where the growth of cancer cells, especially near the boundary, are inhibited strongly by the higher concentration of .
Figure 10.
Illustrations for the case of a source of M on the boundary. Sub-figure (a,b) show the value of at and , respectively; (c) show the dynamics of the maximum and minimum of M, depicting that the value of M on the boundary is substantially higher than it is elsewhere; the cancer cell concentration at and are shown in (d,e); the value of the maximum and minimum of C over the whole time interval is shown in (f). We can see that the cancer cell concentrations approach different steady-state values, positively related to the concentrations of Macrophages.
Figure 10.
Illustrations for the case of a source of M on the boundary. Sub-figure (a,b) show the value of at and , respectively; (c) show the dynamics of the maximum and minimum of M, depicting that the value of M on the boundary is substantially higher than it is elsewhere; the cancer cell concentration at and are shown in (d,e); the value of the maximum and minimum of C over the whole time interval is shown in (f). We can see that the cancer cell concentrations approach different steady-state values, positively related to the concentrations of Macrophages.
Table 1.
Variable names corresponding to .
Table 1.
Variable names corresponding to .
| Variable Name | Biological Meaning (Concentration of) | Scaling Factor |
---|
| | Naive macrophages | |
| M | Macrophages | |
| | Naive T cells | |
| | Helper T cells | |
| | Regulatory T cells | |
| | Cytotoxic T cells and NK cells | |
| | Naive dendritic cells | |
| D | Dendritic cells | |
| C | Cancer cells | |
| N | Necrotic cells | |
| | IFN- | |
| | TGF-, IL-4, IL-10, and IL-13 | |
| | IL-6 and IL-17 | |
| | HMGB1 | |
Table 2.
Mechanical parameter values.
Table 2.
Mechanical parameter values.
Parameter | Name | Value |
---|
| Porosity | 0.2 |
K | Bulk modulus | 40,000 (Pa) [19] |
G | Shear modulus | 30,000 (Pa) [19] |
| Biot effective stress coefficient | 0.7 [20] |
| Hydraulic conductivity | () [19] |
M | Biot modulus | (Pa) [20] |
| Diffusion coefficient for cells | [21] |
| Diffusion coefficient for cytokines | [22] |
| Diffusion coefficient for HMGB1 | [23] |
Table 3.
The total number of cells V, the ratios , the crude estimate of diameters , the linear model’s prediction of diameters, and the diameter of the domain after we apply the solid displacement vector , at different times.
Table 3.
The total number of cells V, the ratios , the crude estimate of diameters , the linear model’s prediction of diameters, and the diameter of the domain after we apply the solid displacement vector , at different times.
| Cell Number V | Ratio | Crude Estimate | Linear Model Estimate | Simulated Domain Diameter |
---|
0 | 33,726 | 1 | 0.01 | 0.01 | 0.01 |
10 | 24,264 | 0.785 | 0.0089 | 0.0087 | 0.0096 |
50 | 41,447 | 1.229 | 0.0111 | 0.0112 | 0.0102 |
100 | 71,778 | 2.128 | 0.0146 | 0.0153 | 0.0116 |
150 | 118,172 | 3.504 | 0.0187 | 0.0202 | 0.014 |
200 | 184,578 | 5.473 | 0.0234 | 0.0259 | 0.0176 |
250 | 273,510 | 8.110 | 0.0285 | 0.0324 | 0.0227 |
300 | 384,156 | 11.39 | 0.0337 | 0.0392 | 0.0291 |
350 | 510,781 | 15.14 | 0.0389 | 0.0459 | 0.0367 |
400 | 642,925 | 19.06 | 0.0437 | 0.0523 | 0.0451 |
450 | 768,233 | 22.78 | 0.0477 | 0.0578 | 0.0535 |
500 | 876,630 | 25.99 | 0.0510 | 0.0622 | 0.0614 |
550 | 963,109 | 28.56 | 0.0534 | 0.0656 | 0.0684 |
600 | 1,027,711 | 30.47 | 0.0552 | 0.0680 | 0.0743 |
650 | 1,073,622 | 31.83 | 0.0564 | 0.0697 | 0.0793 |
700 | 1,105,104 | 32.77 | 0.0572 | 0.0709 | 0.0834 |
750 | 1,126,165 | 33.39 | 0.0578 | 0.0716 | 0.0868 |
800 | 1,140,022 | 33.80 | 0.0581 | 0.0721 | 0.0897 |
850 | 1,149,042 | 34.07 | 0.0584 | 0.0724 | 0.0922 |
900 | 1,154,872 | 34.24 | 0.0585 | 0.0726 | 0.0944 |
950 | 1,158,622 | 34.35 | 0.0586 | 0.0728 | 0.0963 |
1000 | 1,160,215 | 34.40 | 0.0587 | 0.0728 | 0.098 |