We propose a mathematical model for the transport of DNA plasmids from the extracellular matrix up to the cell nucleus. The model couples two phenomena: the electroporation process, describing the cell membrane permeabilization to plasmids and the intracellular transport enhanced by the presence of microtubules. Numerical simulations of cells with arbitrary geometry, in 2D and 3D, and a network of microtubules show numerically the importance of the microtubules and the electroporation on the effectiveness of the DNA transfection, as observed by previous biological data. The paper proposes efficient numerical tools for forthcoming optimized procedures of cell transfection.

We consider two models which were both designed to describe the movement of eukaryotic cells responding to chemical signals. Besides a common standard parabolic equation for the diffusion of a chemoattractant, like chemokines or growth factors, the two models differ for the equations describing the movement of cells. The first model is based on a quasilinear hyperbolic system with damping, the other one on a degenerate parabolic equation. The two models have the same stationary solutions, which may contain some regions with vacuum. We first explain in details how to discretize the quasilinear hyperbolic system through an upwinding technique, which uses an adapted reconstruction, which is able to deal with the transitions to vacuum. Then we concentrate on the analysis of asymptotic preserving properties of the scheme towards a discretization of the parabolic equation, obtained in the large time and large damping limit, in order to present a numerical comparison between the asymptotic behavior of these two models. Finally we perform an accurate numerical comparison of the two models in the time asymptotic regime, which shows that the respective solutions have a quite different behavior for large times.

We introduce a numerical scheme to approximate a quasilinear hyperbolic system which models the movement of cells under the influence of chemotaxis. Since we expect to find solutions which contain vacuum parts, we propose an upwinding scheme which properly handles the presence of vacuum and which gives a good approximation of the time asymptotic states of the system. For this scheme we prove some basic analytical properties and study its stability near some of the steady states of the system. Finally, we present some numerical simulations which show the dependence of the asymptotic behavior of the solutions upon the parameters of the system.