Tumor growth model of ductal carcinoma: from in situ phase to stroma invasion.

This paper aims at modeling breast cancer transition from the in situ stage -when the tumor is confined to the duct- to the invasive phase. Such a transition occurs thanks to the degradation of the duct membrane under the action of specific enzymes so-called matrix metalloproteinases (MMPs). The model consists of advection-reaction equations that hold in the duct and in the surrounding tissue, in order to describe the proliferation and the necrosis of the cancer cells in each subdomain. The divergence of the velocity is given by the increase of the cell densities. Darcy law is imposed in order to close the system. The key-point of the modeling lies in the description of the transmission conditions across the duct. Nonlinear Kedem-Katchalsky transmission conditions across the membrane describe the discontinuity of the pressure as a linear function of the flux. These transmission conditions make it possible to describe the transition from the in situ stage to the invasive phase at the macroscopic level. More precisely, the membrane permeability increases with respect to the local concentration of MMPs. The cancer cells are no more confined to the duct and the tumor invades the surrounding tissue. The model is enriched by the description of nutrients concentration, tumor necrosis factors, and MMPs production. The mathematical model is implemented in a 3D C++ -code, which is based on well-adapted finite difference schemes on Cartesian grid. The membrane interface is described by a level-set, and the transmission conditions are precisely approached at the second order thanks to well-suited sharp stencils. Our continuous approach provides new significant insights in the macroscopic modeling of the breast cancer phase transition, due to the membrane degradation by MMP enzymes.