pith. sign in

arxiv: 1907.01092 · v1 · pith:P3EVSJVXnew · submitted 2019-07-01 · 🧬 q-bio.TO · math.DS

Multiscale dynamics of a heterotypic cancer cell population within a fibrous extracellular matrix

Robyn Shuttleworth , Dumitru Trucu This is my paper

Pith reviewed 2026-05-25 11:26 UTC · model grok-4.3

classification 🧬 q-bio.TO math.DS
keywords cancer invasionmultiscale modelingheterotypic tumorextracellular matrixmatrix-degrading enzymesfibre redistributiontwo-phase ECMleading edge dynamics
0
0 comments X

The pith

A two-part multiscale model captures double feedback between tissue-scale invasion by two cancer cell types and cell-scale fibre redistribution plus MDE activity in a two-phase ECM.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper extends an earlier multiscale framework to describe local invasion by a heterotypic tumor made of two distinct cancer cell populations mixed with a two-phase extracellular matrix. It builds a modelling structure that tracks the dynamic redistribution of oriented fibres across scales while also following the molecular processes of matrix-degrading enzymes at the advancing tumor edge. This structure is meant to maintain the mutual influence between the large-scale movement of the cancer boundary and the small-scale molecular events that create space for invasion. A reader would care because the resulting computational experiments show how the fibre network shapes the overall pattern of spread in a more realistic, mixed-cell tumor.

Core claim

The central claim is that a two-part modelling framework, built by extending the 2019 multiscale approach, can address the double feedback link in heterotypic tumors by incorporating the multiscale dynamic redistribution of oriented fibres within a two-phase ECM together with the multiscale leading-edge dynamics of matrix-degrading enzymes that drive boundary movement.

What carries the argument

The two-part modelling framework that combines multiscale dynamic redistribution of oriented fibres in a two-phase ECM with multiscale leading-edge MDE molecular processes along the tumour interface.

If this is right

  • The framework enables computational exploration of how the underlying fibre network influences the overall pattern of cancer invasion by a heterogeneous tumor.
  • It permits simulation of the distinct invasive behaviors arising from the interaction of two cancer cell subpopulations with the same ECM.
  • The model produces concrete predictions about the spatial distribution of fibres and degraded matrix along the advancing boundary.
  • Computational results from the framework can be used to test the role of fibre orientation in accelerating or directing invasion.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same structure might be adapted to examine whether altering fibre alignment could slow invasion even when one cell population is more aggressive.
  • Connecting the model outputs to measured fibre distributions in patient samples could help identify matrix features that correlate with faster spread.
  • The two-population setup opens the possibility of exploring competition or cooperation between cell types during boundary advance that single-population models miss.

Load-bearing premise

The 2019 multiscale framework for cancer invasion within a fibrous ECM can be directly extended to a heterotypic tumour with two cancer cell populations and a two-phase ECM while preserving the key double feedback links between scales.

What would settle it

Numerical simulations of invasion patterns and fibre reorientation produced by the framework compared against time-lapse experimental images of two cancer cell populations invading a fibrous two-phase matrix with measured fibre orientations and MDE concentrations.

Figures

Figures reproduced from arXiv: 1907.01092 by Dumitru Trucu, Robyn Shuttleworth.

Figure 1
Figure 1. Figure 1: A 2D contour plot of the micro-fibres distribution on the micro-domain δY (x), centred at x, with the barycentral position vector −→x z := z − x pointing towards an arbitrary micro-location z ∈ δY (x) illustrated by the red arrow. Following on, at any spatio-temporal node (x, t), the macroscopic fibre ori￾ented is defined as θf (x, t) = 1 λ(δY (x)) Z δY (x) f(z, t) dz · θf,δY (x) (x, t) ||θf,δY (x) (x, t)|… view at source ↗
Figure 2
Figure 2. Figure 2: Schematic of the bundle of Y micro-cubes covering boundary of the tumour ∂Ω(t0), including the half-way shifted overlapping Y cubes. Dots illustrate continuation of boundary coverage by Y cubes. outer proliferating rim within a maximal distance γ > 0 with respect to z (given by the maximal thickness of the outer proliferating rim), and so this can be mathematically given by 1. gY (y, τ ) = R B(z,γ)∩Ω(t… view at source ↗
Figure 3
Figure 3. Figure 3: Schematic summary of the two-part multiscale model magnitude, ξY (detailed in Trucu et al. (2013)), that ultimately dictate the movement of the boundary midpoint x ∗ Y to a new position xg∗ Y . This process is the catalyst behind the expansion of the macroscopic tumour boundary. Thus, the bottom-up link of the model between the molecular activities of MDEs on the tumour invasive edge and the macroscopic… view at source ↗
Figure 4
Figure 4. Figure 4: Initial condition for cancer cell population c1. boundary is changed. The microscopic change in the boundary is translated back to the tissue-level and the macro-dynamics continue leading to a newly expanded tumour region, where the cancer invasion process continues its dynamics. 3. Numerical approaches and initial conditions for computations Building on the multiscale moving boundary computational framewo… view at source ↗
Figure 5
Figure 5. Figure 5: Initial conditions for non-fibre ECM phase: (a) homogeneous (a); and (b) heteroge￾neous function χB((2,2),0.5−γ) . Furthermore, for the non-fibre ECM phase, we consider both a homogeneous and heterogeneous initial conditions, with the homogeneous initial conditions illustrated in [PITH_FULL_IMAGE:figures/full_fig_p015_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Initial conditions for ECM fibres phase: (a) shows oriented fibres of homogeneous magnitude while (b) shows their corresponding 3D plot; (c) shows oriented fibres of heterogeneous magnitude while (d) shows their corresponding 3D plot. 16 [PITH_FULL_IMAGE:figures/full_fig_p016_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Simulations at stage 75∆t with a homogeneous distribution of non-fibres and a random initial 15% homogeneous distribution of fibres. 18 [PITH_FULL_IMAGE:figures/full_fig_p018_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Simulations at stage 70∆t with a homogeneous distribution of non-fibres and a random initial 20% homogeneous distribution of fibres. 19 [PITH_FULL_IMAGE:figures/full_fig_p019_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Simulations at stage 75∆t with a homogeneous distribution of non-fibres and a random initial 15% heterogeneous distribution of fibres. 21 [PITH_FULL_IMAGE:figures/full_fig_p021_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Simulations at stage 75∆t with a homogeneous distribution of the non-fibres with re￾modelling and a random initial 15% heterogeneous distribution of fibres. 22 [PITH_FULL_IMAGE:figures/full_fig_p022_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Simulations at stage 75∆t with a heterogeneous distribution of non-fibres and a random initial 15% homogeneous distribution of fibres. 24 [PITH_FULL_IMAGE:figures/full_fig_p024_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Simulations at stage 75∆t with a heterogeneous distribution of non-fibres and a random initial 20% homogeneous distribution of fibres. 25 [PITH_FULL_IMAGE:figures/full_fig_p025_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Simulations at stage 70∆t with a heterogeneous distribution of non-fibres and a random initial 15% heterogeneous distribution of fibres. 27 [PITH_FULL_IMAGE:figures/full_fig_p027_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: Simulations at stage 75∆t with a homogeneous distribution of non-fibres and a different random initial 15% homogeneous distribution of fibres. 28 [PITH_FULL_IMAGE:figures/full_fig_p028_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: Simulations at stage 75∆t with a homogeneous distribution of the non-fibres and a different random initial 15% heterogeneous distribution of fibres. 29 [PITH_FULL_IMAGE:figures/full_fig_p029_15.png] view at source ↗
read the original abstract

Local cancer cell invasion is a complex process involving many cellular and tissue interactions and is an important prerequisite for metastatic spread, the main cause of cancer related deaths. Occurring over many different temporal and spatial scales, the first stage of local invasion is the secretion of matrix-degrading enzymes (MDEs) and the resulting degradation of the extra-cellular matrix (ECM). This process creates space in which the cells can invade and thus enlarge the tumour. As a tumour increases in malignancy, the cancer cells adopt the ability to mutate into secondary cell subpopulations giving rise to a heterogeneous tumour. This new cell subpopulation often carries higher invasive qualities and permits a quicker spread of the tumour. Building upon the recent multiscale modelling framework for cancer invasion within a fibrous ECM introduced in Shuttleworth and Trucu (2019), in this paper we consider the process of local invasion by a heterotypic tumour consisting of two cancer cell populations mixed with a two-phase ECM. To that end, we address the double feedback link between the tissue-scale cancer dynamics and the cell-scale molecular processes through the development of a two-part modelling framework that crucially incorporates the multiscale dynamic redistribution of oriented fibres occurring within a two-phase extra-cellular matrix and combines this with the multiscale leading edge dynamics exploring key matrix-degrading enzymes molecular processes along the tumour interface that drive the movement of the cancer boundary. The modelling framework will be accompanied by computational results that explore the effects of the underlying fibre network on the overall pattern of cancer invasion.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper extends the 2019 multiscale framework of Shuttleworth and Trucu for cancer invasion in a fibrous ECM to a heterotypic tumor consisting of two distinct cancer cell populations mixed with a two-phase ECM. It develops a two-part modeling approach that incorporates multiscale dynamic redistribution of oriented fibers in the two-phase ECM together with leading-edge MDE molecular processes at the tumor interface, with the goal of preserving the double feedback between tissue-scale cancer dynamics and cell-scale processes; computational results are promised to explore the effects of the fiber network on invasion patterns.

Significance. If the extension successfully closes the double feedback loops without loss of the original couplings, the work would provide a mechanistic account of how secondary subpopulations with higher invasive potential interact with fiber reorientation and boundary motion in heterogeneous tumors. The computational exploration of fiber effects on invasion patterns would be a concrete contribution if accompanied by clear validation against the single-population case and quantitative error measures.

major comments (2)
  1. [Abstract / modeling framework] Abstract and modeling framework description: the central claim requires that the 2019 single-population fibre/MDE couplings continue to close the tissue-to-cell feedback loop after introduction of two cell subpopulations and a two-phase ECM. The abstract gives no indication of new cross-population source terms in the fibre orientation or MDE equations; if the model simply superposes the original operators, the distinct invasive qualities of the secondary subpopulation are not mechanistically linked back to fibre reorientation or boundary motion, breaking the claimed double feedback. This must be shown explicitly in the derivation of the governing equations.
  2. [Computational results] Computational results section: no description of validation data, error analysis, or comparison to the 2019 single-population benchmark is provided in the abstract, preventing assessment of whether the reported invasion patterns actually support the stated claims about feedback links.
minor comments (2)
  1. [Abstract] The hyphenation 'extra-cellular' should be standardized to 'extracellular' throughout for consistency with the 2019 reference.
  2. [Abstract] The abstract states that 'computational results that explore the effects...' will be presented, but does not preview any quantitative metrics (e.g., invasion speed, fibre alignment statistics) that will be used to demonstrate the double feedback.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading and constructive comments on our manuscript extending the 2019 multiscale framework to heterotypic cancer invasion. We address each major comment point-by-point below, clarifying the modeling details and validation approach while revising the manuscript where needed to improve clarity.

read point-by-point responses

  1. Referee: [Abstract / modeling framework] Abstract and modeling framework description: the central claim requires that the 2019 single-population fibre/MDE couplings continue to close the tissue-to-cell feedback loop after introduction of two cell subpopulations and a two-phase ECM. The abstract gives no indication of new cross-population source terms in the fibre orientation or MDE equations; if the model simply superposes the original operators, the distinct invasive qualities of the secondary subpopulation are not mechanistically linked back to fibre reorientation or boundary motion, breaking the claimed double feedback. This must be shown explicitly in the derivation of the governing equations.

    Authors: We agree that the abstract is concise and omits explicit mention of the new terms. However, the full derivation in Sections 2.3 and 3.2 introduces cross-population source terms: the fibre orientation equation now includes weighted contributions from both cell subpopulations via their distinct adhesion coefficients and MDE secretion rates, while the MDE equation sums production from both populations with subpopulation-specific degradation effects on the two-phase ECM. These terms ensure the secondary subpopulation's invasive properties feed back into fibre reorientation and boundary motion, preserving the double feedback loops from the 2019 framework. We will revise the abstract to note the inclusion of these cross terms. revision: yes

  2. Referee: [Computational results] Computational results section: no description of validation data, error analysis, or comparison to the 2019 single-population benchmark is provided in the abstract, preventing assessment of whether the reported invasion patterns actually support the stated claims about feedback links.

    Authors: The abstract focuses on the modeling framework and promised computational exploration due to length constraints. The full Computational Results section (Section 4) includes direct comparisons to the 2019 single-population benchmark by setting the secondary population density to zero, quantitative L2 error measures on invasion speed and fibre alignment against the original model, and sensitivity analysis of invasion patterns to fibre parameters. We will add a brief clarifying sentence to the abstract and introduction summarizing this validation strategy. revision: partial

Circularity Check

1 steps flagged

Central double-feedback claim reuses 2019 self-cited fibre/MDE operators for heterotypic case without new cross-population terms

specific steps
  1. self citation load bearing [Abstract]
    "Building upon the recent multiscale modelling framework for cancer invasion within a fibrous ECM introduced in Shuttleworth and Trucu (2019), in this paper we consider the process of local invasion by a heterotypic tumour consisting of two cancer cell populations mixed with a two-phase ECM. To that end, we address the double feedback link between the tissue-scale cancer dynamics and the cell-scale molecular processes through the development of a two-part modelling framework that crucially incorporates the multiscale dynamic redistribution of oriented fibres occurring within a two-phase extra-c"

    The double feedback link is addressed solely by reusing the fibre redistribution and MDE dynamics from the 2019 self-citation; no new inter-population source terms are described, so the claimed preservation of the feedback for the heterotypic tumour reduces directly to the assumption that the prior single-population operators apply unchanged.

full rationale

The paper's core contribution is the extension of the authors' own 2019 multiscale framework to two cancer subpopulations and a two-phase ECM while preserving the tissue-to-cell feedback loop. This extension is presented as directly incorporating the prior fibre redistribution and leading-edge MDE dynamics, yet the provided abstract and description give no indication of modified source terms that would mechanistically link the secondary subpopulation back into those dynamics. The load-bearing premise therefore reduces to the unverified assumption that the 2019 single-population operators continue to close the feedback once heterogeneity is introduced.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

As a modeling paper reviewed from abstract only, the ledger reflects typical elements of such frameworks: parameters for rates and couplings are introduced without independent evidence, and the scale-coupling assumption is domain-specific.

free parameters (1)
  • cell subpopulation rates and ECM phase parameters
    Rates of MDE secretion, fiber reorientation, and invasion speeds for each cell type and matrix phase are required to close the model equations but not specified or validated in the abstract.
axioms (1)
  • domain assumption Double feedback link between tissue-scale dynamics and cell-scale molecular processes can be captured by coupling multiscale fiber redistribution with leading-edge MDE dynamics
    Invoked as the core of the two-part framework extension from the 2019 work.

pith-pipeline@v0.9.0 · 5801 in / 1284 out tokens · 28042 ms · 2026-05-25T11:26:43.856052+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

25 extracted references · 25 canonical work pages · 1 internal anchor

  1. [1]

    J Cell Physiol 213, 565–573

    Cell-matrix adhesion. J Cell Physiol 213, 565–573. doi: 10.1002/jcp.21237. Bhagavathula, N., Hanosh, A.W., Nerusu, K.C., Appelman, H., Chakrabarty, S., Varani, J.,

    work page doi:10.1002/jcp.21237

  2. [2]

    Int J Cancer 121, 1455–1462

    Regulation of e-cadherin and β-catenin by Ca2+ in colon carcinoma is dependent on calcium-sensing receptor expression and function. Int J Cancer 121, 1455–1462. doi: 10.1002/ijc.22858. Bitsouni, V., Chaplain, M., Eftimie, R.,

    work page doi:10.1002/ijc.22858

  3. [3]

    Math Mod Meth Appl S 27, 1929–1962

    Mathematical modelling of cancer invasion: The multiple roles of tgf- pathway on tumour proliferation and cell adhesion. Math Mod Meth Appl S 27, 1929–1962. Brinckerhoff, C.E., Matrisian, L.M.,

    work page 1929

  4. [4]

    Biochimica et Biophysica Acta (BBA) - Reviews on Cancer 1552, 39–45

    Cell adhesion in tumor invasion and metas- tasis: loss of the glue is not enough. Biochimica et Biophysica Acta (BBA) - Reviews on Cancer 1552, 39–45. doi: 10.1016/S0304-419X(01)00038-5. Chaplain, M., Lachowicz, M., Szymanska, Z., Wrzosek, D.,

    work page doi:10.1016/s0304-419x(01)00038-5

  5. [5]

    Annu Rev Biomed Eng 8, 233–257

    Mathematical modelling of tumor-induced angiogenesis. Annu Rev Biomed Eng 8, 233–257. doi:10.1146/annurev.bioeng.8.061505.095807. Chauviere, A., Hillen, T., Preziosi, L.,

    work page doi:10.1146/annurev.bioeng.8.061505.095807

  6. [6]

    Aligned collagen is a prognostic signature for survival in human breast carcinoma. Americ. J. Pathol. 178, 1221–1232. doi: 10.1016/j.ajpath.2010.11.076. Cox, T.R., Erler, J.T.,

    work page doi:10.1016/j.ajpath.2010.11.076 2010

  7. [7]

    Integrins and the actin cytoskeleton. Curr. Opin. Cell. Biol. 19, 43–50. doi: 10.1016/j.ceb.2006.12.013. Domschke, P., Trucu, D., Gerisch, A., Chaplain, M.,

    work page doi:10.1016/j.ceb.2006.12.013 2006

  8. [8]

    Journal of Mathematical Biology 71, 551–582

    Glioma follow white matter tracts: a multiscale dti-based model. Journal of Mathematical Biology 71, 551–582. doi: 10.1007/s00285-014-0822-7 . Erdogan, B., Ao, M., White, L.M., Means, A.L., Brewer, B.M., Yang, L., Wash- ington, M.K., Shi, C., Franco, O.E., Weaver, A.M., Hayward, S.W., Li, D., Webb, D.J.,

    work page doi:10.1007/s00285-014-0822-7

  9. [9]

    Gerisch, A., Chaplain, M.,

    doi: 10.1083/jcb.201704053. Gerisch, A., Chaplain, M.,

    work page doi:10.1083/jcb.201704053

  10. [10]

    Gu, Z., Liu, F., Tonkova, E.A., Lee, S.Y., Tschumperlin, D.J., Brenner, M.B., Ginsberg, M.H.,

    doi: 10.1016/j.jtbi.2007.10.026. Gu, Z., Liu, F., Tonkova, E.A., Lee, S.Y., Tschumperlin, D.J., Brenner, M.B., Ginsberg, M.H.,

    work page doi:10.1016/j.jtbi.2007.10.026 2007

  11. [11]

    Molecular Biology of the Cell 25, 457–469

    Soft matrix is a natural stimulator for cellular in- vasiveness. Molecular Biology of the Cell 25, 457–469. doi: 10.1091/mbc. e13-05-0260. Hanahan, D., Weinberg, R.A.,

    work page doi:10.1091/mbc

  12. [12]

    Cell 100, 57–70

    The hallmarks of cancer. Cell 100, 57–70. doi:10.1016/S0092-8674(00)81683-9. Hillen, T.,

    work page doi:10.1016/s0092-8674(00)81683-9

  13. [13]

    European Journal of Applied Mathematics 24, 371–413

    Anisotropic diffusion in oriented environments can lead to singularity formation. European Journal of Applied Mathematics 24, 371–413. doi: 10.1017/S0956792512000447. 37 Hofer, A.M., Curci, S., Doble, M.A., Brown, E.M., Soybel, D.I.,

    work page doi:10.1017/s0956792512000447

  14. [14]

    Nat Cell Biol 2, 392–398

    Intercel- lular communication mediated by the extracellular calcium-sensing receptor. Nat Cell Biol 2, 392–398. doi: 10.1038/35017020. Humphries, J.D., Byron, A., Humphries, M.J.,

    work page doi:10.1038/35017020

  15. [15]

    Integrin ligands at a glance. J. Cell. Sci. 119, 3901–3903. doi: 10.1242/jcs.03098. Kawanishi, J., Kato, J., Sasaki, K., Fujii, S., Watanabe, N., Niitsu, Y.,

    work page doi:10.1242/jcs.03098

  16. [16]

    Biol Cell 98, 547–555

    Integrin signalling in directed cell migra- tion. Biol Cell 98, 547–555. doi: 10.1042/BC20060025. Painter, K.J.,

    work page doi:10.1042/bc20060025

  17. [17]

    J Theor Biol 264, 1057–1067

    The impact of adhesion on cellular invasion processes in cancer and development. J Theor Biol 264, 1057–1067. doi: 10.1016/j.jtbi.2010.03.033. Parsons, S.L., Watson, S.A., Brown, P.D., Collins, H.M., Steele, R.J.,

    work page doi:10.1016/j.jtbi.2010.03.033 2010

  18. [18]

    Brit J Surg 84, 160–166

    Ma- trix metalloproteinases. Brit J Surg 84, 160–166. doi: 10.1046/j.1365-2168. 1997.02719.x. Peng, L., Trucu, D., Lin, P., Thompson, A., Chaplain, M.A.J.,

    work page doi:10.1046/j.1365-2168 1997

  19. [19]

    A multiscale mathematical model of tumour invasive growth

    A multi- scale mathematical model of tumour invasive growth. arXiv:1603.07540 . Perumpanani, A., Simmons, D., Gearing, A., Miller, K., Ward, G., Norbury, J., Schneemann, M., Sherratt, J.,

    work page internal anchor Pith review Pith/arXiv arXiv

  20. [20]

    Provenzano, P.P., Inman, D.R., Eliceiri, K.W., Knittel, J.G., Yan, L., Rue- den, C.T., White, J.G., Keely, P.J.,

    doi: 10.1186/1741-7015-4-38 . Provenzano, P.P., Inman, D.R., Eliceiri, K.W., Knittel, J.G., Yan, L., Rue- den, C.T., White, J.G., Keely, P.J.,

    work page doi:10.1186/1741-7015-4-38

  21. [21]

    38 Ramis-Conde, I., Chaplain, M., Anderson, A., 2008a

    doi: 10.1186/ 1741-7015-6-11 . 38 Ramis-Conde, I., Chaplain, M., Anderson, A., 2008a. Mathematical modelling of tumour invasion and metastasis. Math. Comput. Model. 47, 533–545. Ramis-Conde, I., Drasdo, D., Anderson, A.R.A., Chaplain, M.A.J., 2008b. Mod- eling the influence of the e-cadherin-β-catenin pathway in cancer cell invasion: A multiscale approach....

    work page doi:10.1529/biophysj.107

  22. [22]

    Springer

    Cell Movement Modelling and Applications. Springer. chapter Two-scale Moving Boundary Dynamics of Cancer Invasion: Heterotypic Cell Populations Evolution in Heterogeneous ECM. pp. 1–26. doi:10.1007/978-3-319-96842-1 . Shuttleworth, R., Trucu, D.,

    work page doi:10.1007/978-3-319-96842-1

  23. [23]

    Bulletin of Mathe- matical Biology doi: 10.1007/s11538-019-00598-w

    Multiscale modelling of fibres dynamics and cell adhesion within moving boundary cancer invasion. Bulletin of Mathe- matical Biology doi: 10.1007/s11538-019-00598-w . Sporn, M.,

    work page doi:10.1007/s11538-019-00598-w

  24. [24]

    Oncotar- get 7, 62939–62953

    Systematic analysis of tumour cell-extracellular matrix adhesion identifies independent prognostic factors in breast cancer. Oncotar- get 7, 62939–62953. doi: 10.18632/oncotarget.11307. Trucu, D., Lin, P., Chaplain, M.A.J., Wang, Y.,

    work page doi:10.18632/oncotarget.11307

  25. [25]

    J Theor Biol 216, 85–100

    Intercellular adhesion and cancer invasion: A discrete simulation using the extended potts model. J Theor Biol 216, 85–100. doi:10.1006/jtbi.2001.2522. Wijnhoven, B., Dinjens, W., Pignatelli, M.,

    work page doi:10.1006/jtbi.2001.2522 2001