Coupled mathematical modeling of cisplatin electroporation.

The use of electrochemotherapy (ECT) is a well-established technique to increase the cellular uptake of cytotoxic agents within certain cancer treatment strategies. The study of the mechanisms that take part in this complex process is of high interest to gain a deeper knowledge of it, enabling the improvement of these strategies. In this work, we present a coupled multi-physics electroporation model based on a related previous one, to describe the effect of a set of electric pulses on cisplatin transport across the plasma membrane. The model applies a system of partial differential equations that includes Poisson's equation for the electric field, Nernst-Planck's equation for species transport, Maxwell's tensor and mechanical equilibrium equation for membrane deformation and Smoluchowski's equation for pore creation dynamics. Our numerical results were compared with previous numerical and experimental published data with good qualitative and quantitative agreement. These results indicate that pore aperture is favored at the cell poles by the electric field and mechanical stress forces, giving support to the dominant hypothesis of hydrophilic pore creation as the main mechanism of drug entry during an ECT treatment.

Mathematical modelling of microtumour infiltration based on in vitro experiments.

The present mathematical models of microtumours consider, in general, volumetric growth and spherical tumour invasion shapes. Nevertheless in many cases, such as in gliomas, a need for more accurate delineation of tumour infiltration areas in a patient-specific manner has arisen. The objective of this study was to build a mathematical model able to describe in a case-specific way as well as to predict in a probabilistic way the growth and the real invasion pattern of multicellular tumour spheroids (in vitro model of an avascular microtumour) immersed in a collagen matrix. The two-dimensional theoretical model was represented by a reaction-convection-diffusion equation that considers logistic proliferation, volumetric growth, a rim with proliferative cells at the tumour surface and invasion with diffusive and convective components. Population parameter values of the model were extracted from the experimental dataset and a shape function that describes the invasion area was derived from each experimental case by image processing. New possible and aleatory shape functions were generated by data mining and Monte Carlo tools by means of a satellite EGARCH model, which were fed with all the shape functions of the dataset. Then the main model is used in two different ways: to reproduce the growth and invasion of a given experimental tumour in a case-specific manner when fed with the corresponding shape function (descriptive simulations) or to generate new possible tumour cases that respond to the general population pattern when fed with an aleatory-generated shape function (predictive simulations). Both types of simulations are in good agreement with empirical data, as it was revealed by area quantification and Bland-Altman analysis. This kind of experimental-numerical interaction has wide application potential in designing new strategies able to predict as much as possible the invasive behaviour of a tumour based on its particular characteristics and microenvironment.