Mathematical Modeling of Metastatic Cancer Migration through a Remodeling Extracellular Matrix

The spreading of cancer cells, also known as metastasis, is a lethal hallmark in cancer progression and the primary cause of cancer death. Recent cancer research has suggested that the remodeling of collagen fibers in the extracellular matrix (ECM) of the tumor microenvironment facilitates the migration of cancer cells during metastasis. ECM remodeling refers to the following two procedures: the ECM degradation caused by enzyme matrix metalloproteinases and the ECM alignment due to the cross-linking enzyme lysyl oxidase (LOX). Such modifications of ECM collagen fibers result in changes of ECM physical and biomechanical properties that affect cancer cell migration through the ECM. However, the mechanism of such cancer migration through a remodeling ECM remains not well understood. A mathematical model is proposed in this work to better describe and understand cancer migration by means of ECM remodeling. Effects of LOX are considered to enable transport of enzymes and migration of cells through a dynamic, reactive tumor microenvironment that is modulated during cell migration. For validation cases, the results obtained show comparable trends to previously established models. In novel test cases, the model predicts the impact on ECM remodeling and the overall migration of cancer cells due to the inclusion of LOX, which has not yet been included in previous cancer invasion models.


Introduction
Every year there are approximately 8.2 million cancer-related deaths worldwide.Metastasis is the primary cause of cancer death.Cancer metastasis occurs when the disease reaches its lethal stage via the uncontrolled spreading of cancer cells to invade a nearby connective tissue and other key organs in the human body.The metastasizing primary tumor cells are not the only agents that drive the progression of metastasis.Instead, metastasis is a systematic process that involves the interaction of cancer cells among a community of various biochemical and cellular factors localized in the tumor microenvironment at both the primary and secondary tumor sites [1].Many recent cancer-related studies have pointed out that the primary tumor microenvironment contains many important factors that determine whether the primary tumor progresses and proceeds to metastasize or remains dormant, staying a benign tumor [2][3][4][5][6].In the early stage of metastatic cancer invasion, migration of cancer cells first takes place by breaking away from the primary tumor site and breaching the basement membrane of the tumor.This thin barrier underlying the tumor mass is made up of mostly type IV collagen fibers that separate the tumor from the extracellular matrix (ECM) [7].After perforating the basement membrane, the escaped cancer cells then must invade the ECM and travel through a meshwork of randomly oriented collagen fibers prior to intravasation into a blood vessel [8].While maneuvering through the ECM, the direction and behavior of the migrating cells are greatly influenced by the physical and biomechanical properties of the ECM [9][10][11].
The ECM is a highly dynamic non-cellular environment present in all tissues of a human body [12].The ECM consists generally of water, polysaccharides, and proteins.Various types of fibrous proteins are present in the ECM including collagens, elastins, fibronectins, and laminins; among these, collagen is the most abundant and the main structural protein in the ECM [13].The collagen fibers are constantly being remodeled [10].The chemical remodeling of collagen fibers results in dynamic changes in the physical, chemical, and biomechanical properties of the ECM [14].In this article, ECM remodeling refers to the following two processes: (a) ECM degradation caused by enzyme matrix metalloproteinases (MMPs) and (b) ECM alignment due to the cross-linking enzyme lysyl oxidase (LOX).In a pathological condition, like cancer, such changes in the ECM properties facilitate the movement of cells [15].Via cell-ECM interaction and signaling pathways, these changes elicit cell responses to secrete chemicals like MMPs and LOX to shape the orientation of the ECM from a barrier of randomly oriented fibers into a more propitious fibrous "expressway" for migration [16][17][18].
The family of MMPs are the most prominent enzymes in ECM remodeling, responsible for proteolytic degradation of the collagen fibers in the ECM [19].Here, we focus on the collagen fibers residing in the interstitial ECM, which are mostly collagen type I and are degraded by MMP-14, one of the MMP family members secreted by cancer cells [10].Dysregulation of the activities of MMPs intervenes in almost all cancer types [20] and in the hallmarks of cancer from migration and initial invasion to angiogenesis and metastasis [21].For simplification, we lump all the family of MMPs together in this work, referring to them collectively as MMP hereafter.Due to its important role as a driving factor for cancer progression, MMP has been and continues to be an appealing and promising target of cancer therapy despite the failure of several clinical trials for MMP inhibitors conducted within the past 30 years [20,22,23].
LOX is a member of the lysyl oxidase family of copper-dependent enzymes.LOX oxidizes primary amine substrates to reactive aldehydes.Via this amine oxidase activity, collagens are cross-linked and aligned in the ECM [24].Bundles of cross-linked collagen fibers stiffen the matrix and aid the migration of cancer cells via the growth factor signaling pathway [25].Evidence of thickening and aligned collagen fibers has been observed in areas of active tissue invasion and tumor vasculature [9,16].
Mathematical modeling is a prominent practice in almost all fields of science and engineering in general and in cancer research in particular.Modeling serves as a tool for researchers to formularize and quantitatively understand physical and chemical phenomena observed in reality via the language of mathematics [26].In the last 20 years or so, many mathematical models of cancer invasion have been established and are constantly being improved upon one another .However, no model of cancer invasion has yet considered the enzyme LOX and its promising influence on metastasis via ECM cross-linking.In our mathematical model of metastatic cancer invasion, we not only added a new equation for LOX but also extended the existing models with improved features, showing how LOX affects ECM remodeling and metastatic cancer migration.
In the following sections, we derive and formulate a system of five partial differential equations (PDEs) to describe the local invasion and migration of cancer cells through a remodeling ECM.We also briefly introduce the method used to numerically solve our model in MATLAB (R2017b, MathWorks, Natick, MA, USA, 2017) and detail in Tables 1 and 2 all the parameter values applied in our simulation.We present our results obtained from the numerical simulation and interpret our in silico findings from a biological and pathological perspective.We analyze three case studies to highlight the significance of adding the enzyme LOX into the model and the impact from different modes of haptotaxis caused by LOX on the ECM and on the migration of cancer cells.Lastly, we conclude with remarks on the potential usefulness of the proposed model to further clinical understanding of the critical role of ECM remodeling in the early stages of cancer metastasis.

Methods
The setting of the model established in this work is an in silico metastatic tumor microenvironment (Figure 1).The model aims to unravel the interconnections between the main concepts of metastasis in the ECM: the spreading of cancer cells, the remodeling of collagen fibers, and the reactions kinetics and transport dynamics of the chemicals involved.We start the model time period just after cancer cells have penetrated the basement membrane of the primary tumor (Figure 1A).At this stage, the malignant cancer cells detach away from the primary tumor mass and squeeze through the gaps in the degraded basement membrane.These motile cancer cells proceed to invade the ECM and maneuver their way through a barrier of ECM collagen fibers undergoing remodeling (Figure 1B-D).Such modification of the ECM facilitates the migration of cancer cells via haptotactic sensing and response from cancer cells toward the degraded and cross-linked areas of ECM collagen fibers.Beyond what our model covers, the escaped cancer cells will eventually reach and intravasate nearby blood vessels or invade other connective tissues to travel to other parts of the body and initiate secondary tumors.The model established in this paper is a continuous system of five coupled PDEs describing the dynamics and interaction of cancer cells, collagen fibers, and the enzymes MMP and LOX.Two populations of collagen fibers are considered: those that are oriented randomly and those that have been cross-linked.

Cancer Cells
The population balance of cancer cells in the system is governed by contributions from three main factors: random diffusion, proliferation, and haptotaxis due to the remodeling ECM collagen fibers: where c is the number of cancer cells per volume of the spatial domain x in one dimension at time t.The first term accounts for the diffusive migration of cells with a constant diffusion coefficient D c .We assume any nutrient transport supplied to the tumor from nearby blood vessels is unaffected by the remodeling and thus is ignored in our model.The second term is included to account for cell proliferation using a first-order rate expression for c with specific reaction rate γ as an exponential decay function γ = exp(−x 2 / ), where is an adjustable parameter.The third and fourth terms in the equation implement the haptotaxis effect that directs motility of cancer cells towards the area of remodeling ECM with ρ and ρ cl denoting the haptotaxis parameters for non-cross-linked and cross-linked ECM collagen fibers, respectively.Based on studies on cancer migration using live imaging, cancer cells are prone to move more rapidly on collagen fibers in the collagen enriched area [16,25].Hence, we assume that the effect of haptotaxis on motility of cancer cells toward the cross-linked area of ECM, f cl , will be stronger than toward the area of degraded ECM that has not yet been cross-linked.We incorporate the logistic growth factor (1 in the last three terms of (1) where v 1 , v 2 , and v 3 are the occupied fractions of one unit volume of physical space by the densities of cancer cells c, regular collagen fibers f , and cross-linked collagen fibers f cl , respectively.The logistic growth factor is implemented to ensure the migration of cancer cells only into space that is not already occupied by the c, f , and f cl present in the system domain [40].The fourth term g represents the haptotactic migration of cancer cells attracted toward the cross-linked ECM collagen fibers: The term g is turned on or off (g ≡ 0) in the simulation to investigate the impact of LOX on driving the ECM stiffness via cross-linking and how that influences the overall migration behavior of cancer cells in the system.

Extracellular Matrix Collagen Fibers
Changes in the morphology and the biochemical and physical properties of the tumor-associated ECM make a substantial impact on regulating cancer progression during metastasis [12].Since ECM collagen fibers do not diffuse [52] and can be treated as a static structural support network, none of the terms incorporated in the ECM model relate to motility.Instead, the dynamics of ECM is expressed via the remodeling of its collagen fibers due to LOX and MMP.Under disease conditions like cancer, MMP and LOX are known to be dysregulated and often overexpressed, especially in carcinomas.While MMP digests the collagen fibers and generates spaces within the ECM, LOX cross-links and linearizes the fibers, further organizing the matrix of fibers.Although they have different roles in ECM remodeling, both enzymes contribute greatly to transforming ECM from a barrier of fibers into a more navigable fibrous structure that facilitates the migration of cancer cells through the ECM [10,18].
In our model, we include two species to distinguish the structures and functions of the ECM collagen fibers.The randomly oriented collagen fibers that have not been cross-linked are represented by f .The collagen fibers that are cross-linked, linearized, and aligned are denoted by f cl .Assuming a domain of constant volume, the balances for the two collagen fibers species are where the first terms in (3) and ( 4) account for the degradation of collagen fibers of either type by the concentration of MMP m present in the system domain and α f is the MMP cleavage rate of the collagen fibers.In (3), µ f is the rate constant for production due to the regular collagen fiber synthesis in the ECM.The logistic growth fraction (1 again is incorporated in the source term in (3) to check for the availability of unoccupied space for the normal synthesis of collagen fibers to take place.The term h is included in both ( 3) and ( 4) indicating the rate of converting regular collagen fibers into cross-linked fibers.The h term depends on the presence of LOX l and collagen fibers according to where β f is the reaction rate for cross-linking of collagen fibers by LOX.Like the g term in Equation ( 1), the h term in ( 3) and ( 4) is turned on and off in the simulation depending on the existence of LOX in the system.

Enzyme MMP
The evolution of the chemical signal MMP in the system is modeled using a reaction-diffusion equation where the first term accounts for MMP diffusion with a constant diffusion coefficient D m , the second term considers the natural decay of MMP with α m denoting the MMP decay rate constant, and the third term accounts for the secretion of MMP by cancer cells with a production rate constant β m .In reality, MMP is secreted by many different types of cells in the malignant tumor microenvironment including cancer-associated fibroblasts (CAFs), inflammatory cells, macrophages, and cancer cells [53].However, the collagen fibers residing in the interstitial ECM are mostly collagen type I, which is degraded mostly by MMP-1 secreted by CAFs and MMP-14 secreted by cancer cells [10,54].In the current model, we have not included the existence of CAFs.The secretion of MMP, generally, is considered from cancer cells.

Enzyme LOX
Similar to that for MMP, the reaction-diffusion equation for LOX is where D l is the chemical diffusion coefficient of LOX, α l is the natural decay rate constant for LOX, and β l is the rate constant for the secretion rate of LOX by cells.

Nondimensionalization and Parameter Values
To numerically solve the system of PDEs, nondimensionalization is first performed.Nondimensionalization helps establish fewer parameters for efficiently analyzing the system of PDEs.The dimensionless variables are where parameters L, τ, c o , f o , m o , and l o are appropriate reference values for scaling x, t, c, f and f cl , m, and l, respectively.Table 1 summarizes the reference values and parameters that are available from the literature.Considering the early stages of metastatic cancer cell invasion, the length scale L is taken to be 1 cm but could be in the range of 0.1 to 1 cm.The residence time τ is taken to be 32 h, which is a representative average in vitro doubling time for the well-established human cancer cell lines of A549 lung carcinoma cells, U87MG glioma cells, and MCF-7 and MDA-MB 231 breast cancer cells [55][56][57][58].
The diffusion coefficient of the cancer cell D c has been previously determined in an experiment of cell movement by [59].The diffusion coefficient for the MMP D m is assumed to be on the order of 10 −9 cm 2 s −1 [60] as is used in [33].The haptotaxis parameter ρ is taken from [32].The reference chemical diffusion coefficient D is from [59].The reference concentration of ECM collagen fiber f o is in the range of 10 −11 to 10 −8 M based on [61].The lower end value of f o = 10 −11 , is applied in our model.A value of 0.1 × 10 −9 M is taken from [47] as an appropriate reference chemical concentration for both MMP and LOX, m o and l o , respectively.Introducing the dimensionless quantities defined in ( 8) into ( 1)-( 7), the model can be written as The dimensionless parameters that emerge within parentheses are defined in Table 2. Values for these dimensionless parameters are obtained either through calculation (if sufficient information is provided from literature listed in Table 1, these are labeled as "Calculated") or through approximate values that are reasonable for the model (those with related values obtained from other previous models are listed as "Estimated from" followed by the sources).

Initial and Boundary Conditions
The initial condition for the cancer cells is set to the same initial condition proposed in several previous cancer invasion mathematical models [35,37,40,47].The center of the primary tumor mass resides at the left edge of the system domain x = 0. Initially, it is assumed that a fixed cluster of cancer cells already exists in the system domain from x ∈ [0, 0.25].Additionally, at x = 0.25 is the edge where the basement membrane of the primary tumor resides before being degraded.The initial distribution of the population density of cancer cells in the system is given by ĉ where here σ is a positive constant given a value of 0.01 [33].
Initially, the ECM is a mesh of randomly oriented collagen fibers that are not yet cross-linked: To ensure the physical space of the system domain is not congested, the initial condition for the concentration of ECM collagen fibers must satisfy ĉ( x, 0) We assume that MMP and LOX are absent from the system initially: Zero flux boundary conditions for symmetry are imposed on the left edge of the spatial domain ( x = 0) for all components ( ĉ, f , fcl , m, and l).On the right edge of the spatial domain, it is assumed that there are no cancer cells, cross-linked ECM fibers, MMP, or LOX ( ĉ, fcl , m, and l, respectively).
We assume that there exists a constant amount of non-cross-linked ECM fibers, f , on the right edge.Hence, the set of boundary conditions imposed on the system is

Numerical Methods and Code Repository
The system of coupled PDEs ( 9)-( 23) in this paper is numerically solved in the domain Ω = [0, 1] × (0, 20] utilizing the pdepe function, an internal PDE solver in MATLAB that discretizes the equations in space to obtain a system of ordinary differential equations in time that is then solved along the discrete grid points.This function can handle solving initial-boundary value problems for systems of parabolic and elliptic PDEs in one spatial variable x and time t [63].The PDEs that the function pdepe can solve must follow the general form Our model follows this form with n = 0, indicating rectangular coordinates for pdepe.
To enable code reuse, we wrote the model in MATLAB and shared the code including parameter values and documentation in an open-source software repository [64].

Results and Discussion
In this section, we present and analyze in silico experimental results in one dimension for three case studies of the model: (I) LOX is absent and has no effects on ECM and cancer cells (for all time steps: l ≡ 0, f cl ≡ 0, h ≡ 0, and g ≡ 0); (II) LOX is present but not coupled to the haptotactic migration of cancer cells toward the cross-linked ECM fibers (for all time steps: g ≡ 0); (III) LOX is present and is coupled to the haptotactic migration of cancer cells.Note that for the entirety of the Section 3, we have dropped the "hat" notation for convenience.As described in Section 2.7, only dimensionless quantities are used in the model equations.

Case I: LOX Is Absent
In this case study, LOX is not present in the system, so l ≡ 0. This leads to no cross-linking of ECM collagen fibers, so h ≡ 0 and f cl ≡ 0. Therefore, the term g that represents the haptotactic migration of cancer cells toward the cross-linked ECM collagen fibers is also turned off in the simulation, hence g ≡ 0. The system of PDEs is thus left with only three active species: c, f , and m.
The purpose of excluding LOX from the model in Case I is to validate our model before further exploration by comparing our results with those obtained from a highly cited published cancer invasion model [33].The model in [33] included three PDEs for c, f , and m.Our results (not shown here) in Case I for the same three components combined with the same initial and boundary conditions used in [33] show the same results as the published model.
We set the initial condition as m = 0, deviating from [33], to consider the onset of MMP secretion from cancer cells.Figure 2 shows the spatial profiles of the tumor microenvironment at different dimensionless simulation times (t = 0, 1, 10, and 20).The results obtained for the population density of cancer cells capture a slow detachment of the cells away from the primary tumor mass (the left edge of the spatial domain) as time evolves.We also observe that ECM fiber concentration is low when and where MMP concentration is high.Such trends are expected due to the fact that MMP degrades and digests ECM collagen fibers.

Case II: LOX Is Present But Not Coupled to the Haptotactic Migration of Cancer Cells
In Case II, LOX l is included in the system.Since LOX is present, the rate of cross-linking of collagen fibers h is activated in the model.However, the haptotactic migration effect on the population density of cancer cells in (9) remains inactive in the simulation (g ≡ 0).
Our primary aim for the settings in Case II is to confirm that LOX can perform its main function on the ECM via the activation of h, which transforms regular collagen fibers to cross-linked fibers.Meanwhile, in the absence of the attraction of cancer cells toward the cross-linked fibers (g ≡ 0), the migration of cancer cells is only influenced by haptotactic migration toward the pristine collagen fibers f , which are present in lower concentrations compared with Case I since they can be reduced by cross-linking.
In Figure 3, the slow detachment of cancer cells away from the primary tumor mass and the effect of ECM degradation by MMP are captured in similar trends as in the Case I results and as in [33].The main differences between the results of Cases I and II are the existence of the cross-linked ECM fibers and the f cl dynamics.With h active and g inactive, the concentration of regular ECM fibers decreases more quickly and to a larger penetration distance due to cross-linking (compare the regular ECM fiber curve in Figure 3D with that in Figure 2D).This slightly reduces the cancer cell haptotactic driving force.Additionally, the volume occupancy fractions of other species are affected by f cl in Case II.

Case III: LOX Is Present and Is Coupled to the Haptotactic Migration of Cancer Cells
In Case III, LOX is present in the system, and both the h and g terms are turned on in the model.In the simulation results (Figure 4), besides the recurrence of the expected phenomena from Cases I and II, a peak for cancer cells evolves over time away from the primary tumor (x > 0.25) (Figure 4C,D).The location of the peak of cancer cell population density corresponds to trailing the wave of cross-linked ECM fiber concentration.This suggests that cancer cells that have invaded the ECM via migration are clustered in the area where there is a high concentration of cross-linked ECM collagen fibers.Such behavior is the consequence of the haptotaxis effect caused by the directional motility of cancer cells toward the cross-linked ECM collagen fibers.Additionally, a smaller, secondary peak of regular ECM fibers appears to the left of the migrating cell front where conditions are favorable for the secretion of new regular ECM fibers.
Results of Case III have confirmed the capability of our new model to successfully implement the extended features to capture the cross-linking effect that LOX performs on the ECM.The results demonstrate how cross-linked fibers enhance the overall migration of cancer cells.

Sensitivity Analysis
Local sensitivity assesses the impact of variations in each parameter on model outputs.The Case III conditions discussed in Section 3.3 are used.All of the dimensionless parameters listed in Table 2 are considered for perturbations in the local sensitivity analysis, except for v1 , v2 , and v3 .The nominal values for the parameters are those listed in Table 2.The sensitivity analysis is conducted by varying each parameter, P j , at a time by a small perturbation ∆P j while keeping the other parameters fixed.We define the time-dependent local sensitivity index S(t) as where Y i (P, x, t) is the model prediction of the output variable i at position x and time t evaluated at parameter set P. The parameters and model equations are all dimensionless, so (25) is already normalized.Maximum deviations sorted with respect to x from the nominal case Y i (P j , P k =j , x, t) allow for the straightforward comparison of offsets over time due to parameter variations.∆P j is taken to be an increase of 10%.
We consider S = 1 to be a baseline threshold for categorizing parameter effects to be sensitive (i.e., a parameter is labeled "sensitive" if S ≥ 1).In Figure 5, we show that, besides Dc , Dm , Dl , α f , and βm , the model outputs up to simulation time t = 20 are relatively insensitive to the other parameters.c is the output most sensitive to changes in input parameters, especially Dc , Dm , Dl , α f , and βm (Figure 5A).The two local sensitivity curves of parameters βm and α f consistently overlap each other, indicating that any changes in these two parameters influence the model outputs with the same magnitudes over time (Figure 5A-C).Additionally, βm and α f impact the model outputs of c, f , and f cl more profoundly than others.The outputs of m and l are the most sensitive to the input parameters relating to chemical production by cancer cells βm and βl , respectively (Figure 5D-E).From visual inspection of Figure 5F, it is clear that the maximum deviation from the baseline with respect to x occurs for t = 20 near x = 0 and that βl is the most sensitive parameter for this scenario.The sensitivity index accounts for the same visual detection in a concise metric.The maximum with respect to x is considered because the peaks may shift in different parameter variation scenarios compared to the locations of the peaks in Figure 4, making it difficult to choose an x-position a priori at which the sensitivities can be fairly compared.

Conclusions
A model for reactions, diffusion, migration, and proliferation in the ECM undergoing dynamic remodeling has been proposed and analyzed.Case I, which only covers cancer cells, ECM density, and MMP concentration, has been verified and validated via comparison with a previous model [33].The extended features of considering the enzyme LOX and its effect on ECM and cancer migration are successfully implemented in our new model demonstrated in Cases II and III.Simulation results of Case III have confirmed the capability of the model to capture the cross-linking effect that LOX performs on ECM and how cross-linked fibers enhance the overall migration of cancer cells.Based on the current model, an additional PDE could be included to potentially aid in optimizing drug transportation into the tumor through the remodeling ECM.This could result in a better understanding of the various processes that take place within the specific microenvironment and in the determination of tissue and/or chemical factors that may inhibit an administered drug from infiltrating the tumor [65].Further, accounting for additional biological activities could allow for the model to represent a more complex and even more realistic system.Our model, in particular, has included the role of LOX in mathematical modeling of cancer migration and has raised the attention of the importance of this enzyme.The model provides a fundamental understanding that could facilitate predictions of new therapeutic development, for example, for inhibiting the effects of LOX or altering or slowing the remodeling rate of the ECM to slow down or prevent metastasis.

Figure 1 .
Figure 1.Dynamics in a metastatic tumor microenvironment: (A) The basement membrane around the edge of the primary tumor (cluster of cancer cells in red) has already been perforated by cancer cells.The surrounding collagen fibers are randomly oriented in the extracellular matrix (ECM).(B) Enzymes metalloproteinase (MMP) and lysyl oxidase (LOX) are secreted by cancer cells to degrade and cross-link collagen fibers.MMP generates spaces for cancer cells to begin detaching away from the primary tumor mass to invade the ECM.(C) Meanwhile, aligned and cross-linked collagen fibers form a fibrous pathway along which cancer cells prefer to travel.(D) Collagen fibers continue to be cross-linked, aiding the maneuvers of cancer cells further through the matrix.

Figure 2 .
Figure 2. One-dimensional numerical results for Case I when there is zero concentration of LOX (l ≡ 0) and hence, zero cross-linked ECM collagen fibers f cl in the system.Results are snapshots of the system dynamics at four simulation times: (A) t = 0, (B) t = 1, (C) t = 10, and (D) t = 20.For all four plots, the horizontal axis x indicates the dimensionless spatial position, and the vertical axis y indicates the dimensionless population density or concentration of the species listed in the legend.

Figure 3 .
Figure 3. One-dimensional numerical results for Case II when LOX (l ≡ 0) and its effect only on the ECM collagen fibers and not on the cancer haptotaxis (g ≡ 0) are considered in the modeled system.The results are snapshots of the system dynamics at four simulation times: (A) t = 0, (B) t = 1, (C) t = 10, and (D) t = 20.For all four plots, the horizontal axis x indicates the dimensionless spatial position, and the vertical axis y indicates the dimensionless population density or concentration of the species listed in the legend.

Figure 4 .
Figure 4. One-dimensional numerical results for Case III when LOX (l ≡ 0) and its effects both on the ECM collagen fibers and cancer cells motility are considered in the modeled system.Results are snapshots of the system dynamics at four simulation times: (A) t = 0, (B) t = 1, (C) t = 10, and (D) t = 20.For all four plots, the horizontal axis x indicates the dimensionless spatial position, and the vertical axis y indicates the dimensionless population density or concentration of the species listed in the legend.

Figure 5 .
Figure 5. Local sensitivity index as a function of time assessing the impacts of 10% one-at-a-time increases in dimensionless parameters listed in Table 2 on the following model output variables: (A) the population density of cancer cells, (B) the concentration of regular ECM fibers, (C) the concentration of cross-linked ECM fibers, (D) the concentration of MMP, and (E) the concentration of LOX.(F) Model output profiles for concentration of LOX as a function of position x at at simulation time t = 20 corresponding to 10% changes in each parameter input.In (A-E), the baseline marks the threshold value of S = 1.In (F) the baseline marks the nominal LOX profile.

Table 1 .
Parameter values available from the literature used in the model of metastatic invasion of cancer through remodeling ECM.

Table 2 .
Dimensionless expression and values of parameters used in MATLAB simulation of metastatic invasion of cancer through a remodeling ECM.