Mathematical modeling and numerical simulation of a multiscale cancer invasion of host tissue

Abstract

In this study, we develop an advection-reaction-difusion system of partial diferential equations (PDEs) to describe interactions between tumor cells and extracellular matrix (ECM) at the macroscopic level. At the subcellular level, we model the interaction of proteolytic enzymes and the ECM with a set of ordinary diferential equations (ODEs). A contractivity function is used to couple the macroscopic and microscopic events. The model is supplemented with nutrients supply from the underlying tissue. These PDE-ODE systems of equations model the on-set of tumor cell invasion of the host extracellular matrix. The model accounts for diferent time and spatial scales at the macroscopic and microscopic levels. Contact inhibition between the tumor cells and the tumor micro-environment are also accounted for through a nonlinear density-dependent difusion and haptotaxis coefcients. In the numerical simulations, we use a nonstandard fnite diference method to illustrate the model predictions. Qualitatively, our results confrm the three distinct layers of proliferating, quiescent and necrotic cells as observed in multicellular spheroids experiments